- Open Access
Construction of a lncRNA-associated competing endogenous RNA regulatory network after traumatic brain injury in mouse
Molecular Brain volume 15, Article number: 40 (2022)
Traumatic brain injury (TBI) is a major public health problem worldwide which causes high mortality and disability. Functioning as microRNA (miRNA) sponges, long non-coding RNA (lncRNA) regulates the expression of protein-coding genes in a competing endogenous RNA (ceRNA) network. However, the lncRNA-associated ceRNA in TBI remains unclear. In this study, we processed the raw SRR files of mice cortex samples of sham injury (n = 3) and TBI groups (n = 3) to count files. Then, the expression profiles of lncRNAs and mRNAs were identified, and 86 differentially expressed (DE) lncRNAs and 1201 DEmRNAs between sham and TBI groups were identified. The DEmRNAs were used to perform enrichment analyses. Next, a lncRNA-miRNA-mRNA regulatory ceRNA network was constructed. The network consisted of 23 mRNAs, 5 miRNAs and 2 lncRNAs. The expression alternations of the 5 miRNAs were validated via qRT-PCR. The subnetwork of hub lncRNA Neat1 was extracted. We identified a potential inflammatory associated regulatory axis: Neat1/miR-31-5p/Myd88 axis. The PPI network based on DEmRNA involved in ceRNA network was constructed PPI networks to identify the hub genes. Finally, DElncRNAs and DEmRNAs were selected randomly and validated by qRT-PCR. In conclusion, with the lncRNA-miRNA-mRNA ceRNA network provided above, we can improve our understanding of the regulatory mechanisms and interaction among lncRNAs, miRNAs and mRNAs in TBI process.
Traumatic brain injury (TBI) is a major public health problem all over the world. It is estimated that 10 million people worldwide suffer from TBI annually [1, 2]. WHO estimates that TBI will be the third most common cause of death and disability in global terms by 2020 . TBI can be divided into initial injury and secondary brain injury . Initial injury causes cells to die in the moment of mechanical shock. Secondary injury refers to the biochemical and physiological changes which would cause the apoptosis and death of neural cells, the formation of brain edema, the disruption of the blood–brain barrier, and the subsequent neurological disorders . Various efforts have been made over the past decades, however, TBI remains a disease with high mortality and disability, which brings a huge burden on the families and society . The complex molecular mechanisms behind brain injury need to be further explained.
Long non-coding RNAs (lncRNAs) belong to the non-coding RNA family with more than 200 nucleotides in length . Although lncRNAs do not encode any protein products, many studies have proved that they could be involved in the regulation of gene expression at epigenetic, transcriptional or post-transcriptional levels [8, 9]. Several studies have shown that lncRNAs can act as miRNA sponges in the ceRNA regulatory network and further regulate protein expression 10. In recent years, more and more studies have focused on the role of lncRNA-related ceRNA regulatory networks in intracranial aneurysm , cerebral infarction , type 2 diabetes  and various types of cancer [14, 15, 16]. Several lncRNAs, such as lncRNA GM12371, lncRNA Evf2 and lncRNA Pinky, were reported to be involved in regulation of synaptic functions and neurodevelopment [17, 18, 19]. Moreover, some lncRNAs were found to be specifically expressed in nervous system, such as cerebellar cortex [20, 21]. These findings indicated that lncRNAs could serve as a major role in the nervous system. However, the potential role of the lncRNA-associated ceRNA regulatory network in TBI remains unclear.
In this study, we investigated the expression of lncRNAs and mRNAs in the injury cortex 24 h after TBI in mice via RNA-seq data analyses and further constructed a lncRNA-miRNA-mRNA ceRNA regulatory network. Firstly, we identified differentially expressed lncRNAs (DElncRNA) and mRNAs (DEmRNAs) between normal brain cortex and injury cortex. The DEmRNAs were used to perform enrichment analyses. Then, we developed a lncRNA-miRNA-mRNA ceRNA regulatory network using integrated bioinformatics analysis. The network contained 23 mRNAs, 5 miRNAs and 2 lncRNAs. The expression alternations of the 5 miRNAs were validated via qRT-PCR. Next, Subnetwork of Neat1 was further analyzed, and inflammatory related Neat1/miR-31-5p/Myd88 axis was identified. Gene Oncology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted to evaluate the mRNAs in the network to identify additional biological functions. Furthermore, protein–protein interaction (PPI) network was constructed, and hub genes with key roles in the PPI network were identified. Finally, the difference in expression of four DElncRNA and four DEmRNA was validated via qRT-PCR. The construction of the lncRNA-miRNA-mRNA ceRNA network may provide insight into the regulatory mechanism of TBI.
Materials and methods
All animal procedures were performed in accordance with the Guidelines for Care and Use of Laboratory Animals of the Zhejiang University and approved by the Animal Ethics Committee of the Zhejiang University. Adult male C57BL/6 mice were purchased from Shanghai Slac Experimental Animal Center.
Controlled cortical impact (CCI) model of TBI
Both sham-injury and TBI mice (12–16 weeks old and 22–25 g) were anesthetized by intraperitoneal injection of 1% pentobarbital sodium solution (0.1 ml/20 g). Six mice were divided into sham group (n = 3) and TBI group (n = 3) according to the random number table. The mice heads were shaved and disinfected by wiping with iodophors. The animals were mounted in a prone position on a stereotaxic instrument (RWD Life Science, China) and fixed with auxiliary ear and incisor bars. Using the sterile surgical procedures, the mice received a midline cranial skin incision, and the scalp was retracted to expose the skull. Then, the mice received a right lateral craniotomy (3.5 mm in diameter) with 1.5 mm lateral to sagittal suture and 2.0 mm posterior to the bregma remaining the dura mater intact, which was performed with a motorized drill. The cortical impact was performed at a velocity of 5.0 m/s, a depth of 2.0 mm below the cortical surface, and an impact duration of 180 ms, which would cause a moderately severe contusion in the sensorimotor cortex. Medical grade cyanoacrylate gel was applied to the exposed dura mater and skull surface after CCI. The hole on the skull was filled with bone wax, and the skin incision was sutured with an absorbable suture. The antibiotic ointment was applied to the suture area. The mice were wrapped in an electric blanket to maintain the body temperature and transferred to a clean cage. Sham injury mice underwent the same craniotomy and postoperative care procedures. All the mice were anesthetized with 1% pentobarbital sodium solution (0.1 ml/20 g) 24 h post-TBI. The mice were transcardially perfused with 10 mL 4℃ 0.9% saline. The ipsilateral cortex around the contusion site was dissected rapidly and stored at − 80℃. All the 6 cortex samples were used for qRT-PCR.
RNA-seq data processing
The SRR files (SRR3271216-SRR3271221) of mice cortex samples of sham injury (n = 3) and TBI groups (n = 3) were retrieved from Gene Expression Omnibus (GEO) NCBI , which were converted to ‘fastq’ format data using sratoolkit (version 2.9.6). Using FastQC (version 0.11.9), quality control of the raw sequence data was performed. Low-quality reads and adaptors were removed by the trimmomatic software (version 0.36). The reference genome and genome annotation files of mouse (Release M24 GRCm38.p6) were downloaded from GENCODE (www.gencodegenes.org). Then, the high-quality reads were aligned to the reference genome with hisat2 (version 2.0.0) . The generated SAM (Sequence Alignment/Map) files were converted to the BAM (Binary Alignment/Map) format by using samtools (version 1.10) . The BAM files were converted to counts files by using featureCounts in subread software (version 2.0.0) . Then, the “.csv” files were generated, each consisting of all gene counts for a particular sample. All of the “.csv” files were combined, and a single file with sample names depicted as columns and gene names depicted as rows was obtained. LncRNAs and mRNAs were annotated according to the genome annotation files.
Identification of differentially expressed lncRNAs and mRNAs
Using the DESeq2 package , the counts files were normalized and differentially expressed analyses were performed. lncRNAs and mRNAs with |log2(fold change)|> 1 and adjusted P-values < 0.05 were considered DElncRNAs and DEmRNAs.
Functional enrichment analysis
Both the Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using the Database for Annotation, Visualization, and Integrated Discovery platform (DAVID 6.8; https://david.ncifcrf.gov/) . DEmRNAs or the mRNAs involved in ceRNA network were used for GO and KEGG pathway enrichment analyses. The GO terms and pathways, with corrected P-values < 0.05 using the Benjamini method, were considered significant functional categories.
Construction of ceRNA network
Firstly, the interactions of DElncRNAs and targeted miRNAs with high stringency (> = 3) were identified using starBase v2.0 . Then, miRNA-targeted mRNAs were identified using miRTarBase, a highly reliable miRNA reference database . The miRTarBase database has accumulated more than three hundred and sixty thousand miRNA-targeted interactions which were experimentally validated by western blot, reporter assay, microarray and next-generation sequencing experiments. Moreover, the targeted mRNAs that were not differentially expressed between control cortex and TBI cortex samples were filtered out. The ceRNA network was constructed and viewed using Cytoscape (version 3.7.2; http://www.cytoscape.org/) .
Construction of protein–protein interaction (PPI) network
The Search Tool for the Retrieval of Interacting Genes (STRING version 11.0; www.string-db.org) was used to construct the PPI network of DEmRNAs involved in the ceRNA network . The interactions with a score more than 0.4 were included. The PPI network was constructed and viewed using Cytoscape (version 3.7.2).
RNA preparation and qRT-PCR
The brain cortex tissues of sham and TBI groups were cut into small chunks and rinsed with PBS, which were further dissolved with TRIzol reagent (Invitrogen, USA) to acquire total RNA. Then, using Reverse Transcription Kit (TaKaRa, Osaka, Japan), the isolated RNA was reverse-transcribed to cDNA. The qRT-PCR analyses were performed on an Applied Biosystems 7500 Fast Real-Time PCR System (Roche, Basel, Switzerland) with SYBR Green supermix (172-5150, Bio-Rad, Shanghai, China). All the experiments were performed and analyzed in triplicate. Five miRNAs (mmu-miR-377-3p, mmu-miR-185-5p, mmu-miR-107-3p, mmu-miR-31-5p, and mmu-miR-130a-3p), four DElncRNAs (lncRNA C030018K13Rik, lncRNA Gm36823, lncRNA H19 and lncRNA Mir155hg) and four DEmRNA (P2ry12, Hes5, Cxcr2 and Mmp12) were randomly selected to perform qRT-PCR for validation. The primer sequences used in this study were shown in Additional file 1. The lncRNA, miRNAs and mRNA expression levels were calculated according to the 2−ΔΔCt method.
GraphPad Prism (version 6.0, GraphPad Software, San Diego, CA, USA) and R language (3.4.3) were used for statistical analysis. Statistical differences were determined by Student’s t test for two-group comparisons. A P-value < 0.05 was considered statistically significant. Barplots were generated by GraphPad Prism. Other plots, including bubble charts, volcano plots and heat maps, were produced by R language (3.4.3).
LncRNA expression profile
The total detected reads for samples from sham (n = 3) and TBI (n = 3) groups were 100,290,212 and 96,819,024, respectively. The RNA-seq analysis of 3 normal brain cortex tissues and 3 TBI cortex tissues identified 9945 lncRNAs. After deleting lncRNAs with average read count below 1 across all samples, there were 5700 lncRNAs remaining. The distribution of these lncRNAs in the chromosomes was shown in Fig. 1A. Of these lncRNAs, 86 lncRNAs were differentially expressed (|log2FC|> 1 and adjusted P-value < 0.05) in mouse cortex 24 h post-TBI relative to normal brain cortex, which contained 47 upregulated lncRNAs and 39 downregulated lncRNAs (Fig. 1B, Additional file 2). Hierarchical clustering analysis showed that lncRNA expression profiles in the injured cortex were significantly different from those in the control group (Fig. 1C). The differentially expressed lncRNAs (DElncRNAs) were distributed on all the chromosomes, although the distribution on the chromosomes was not equal (Fig. 1D). Chromosome 7 had the largest number of DElncRNAs, of which 6 were upregulated and 2 were downregulated, accounting for 9.3% (8/86) of all DElncRNAs.
mRNA expression profile
The RNA-seq analysis of 3 normal brain cortex tissues and 3 TBI cortex tissues also identified 21,807 mRNAs. After deleting mRNAs with average read count below 1 across all samples, there were 17,711 mRNAs remaining. The distribution of these mRNAs on all the chromosomes was shown in Fig. 2A. 1201 differentially expressed mRNAs (DEmRNAs) were identified between mouse cortex 24 h post-TBI and normal brain cortex, which contained 1076 upregulated mRNAs and 125 downregulated mRNAs (Fig. 2B, Additional file 2). Hierarchical clustering analysis showed that mRNA expression profiles in the injured cortex were significantly different from those in the normal cortex (Fig. 2C). The DEmRNAs were distributed on all the chromosomes, although the distribution on the chromosomes was not equal (Fig. 2D). Chromosome 7 had the largest number of DEmRNAs, of which 103 were upregulated and 8 were downregulated, accounting for 9.24% (111/1201) of all DEmRNAs.
Gene ontology and pathway analysis of DEmRNAs
To explore the potential functional implication of the 1,201 DEmRNAs, GO enrichment and KEGG pathway analyses were performed. In the GO enrichment analysis, a total of 466 enriched GO terms in the Biological Process (BP) were identified. The top 20 significantly enriched terms were shown (Fig. 3A). The DEmRNAs were primarily enriched in immune inflammatory response-related BPs, such as “immune system process”, “inflammatory response”, “neutrophil chemotaxis”, “immune response” and “response to lipopolysaccharide”. Furthermore, we also found that these DEmRNAs were also enriched in “cell response to interferon-beta”, “cellular response to tumor necrosis factor”, “cellular response to interleukin-1” and “cellular response to interferon-gamma”, which indicated that the neural cells may be actively adapting to the drastically changing microenvironment. In addition, a total 62 enriched pathways were identified via the KEGG pathway analysis. The top 20 significantly enriched pathways were shown (Fig. 3B). Among these pathways, “TNF signaling pathway” and “NF-kappa B signaling pathway” were corresponding to the enriched GO terms, which indicated that the two pathways may play significant roles in acute phase of TBI.
Construction of a ceRNA regulatory network in TBI
To elucidate the regulatory mechanism of TBI, a lncRNA‐miRNA‐mRNA ceRNA network of TBI was developed. As for integrated ceRNA network, lncRNA‐miRNA‐mRNA axis consists of two forms: (1) downregulated lncRNAs, upregulated miRNAs, and downregulated mRNAs; (2) upregulated lncRNAs, downregulated miRNAs, and upregulated mRNAs. According to the changes in the expression of previously reported miRNAs after TBI [32, 33], we included the five miRNAs (mmu-miR-377-3p, mmu-miR-185-5p, mmu-miR-107-3p, mmu-miR-31-5p, and mmu-miR-130a-3p) with the correct expression trends to develop the lncRNA‐miRNA‐mRNA ceRNA network (Fig. 3C). The ceRNA network in this study only contained upregulated lncRNAs, downregulated miRNAs and upregulated mRNAs. The expression levels of mmu-miR-377-3p, mmu-miR-185-5p, mmu-miR-107-3p, mmu-miR-31-5p, and mmu-miR-130a-3p in brain tissues were reported to downregulate after TBI in previous studies [32, 33]. Moreover, we used qRT-PCR to further confirm the expression alternations of the five miRNAs in network (Fig. 3E–I). The lncRNA‐miRNA and miRNA‐mRNA relationship pairs were listed in Additional file 3 and Additional file 4. The network was constructed with 30 nodes (23 mRNAs, 5 miRNAs and 2 lncRNAs) and 28 edges. Blue squares represented lncRNAs, tan triangles represented miRNAs, and red circles represented mRNAs. Among the network, lncRNA Neat1 directly connected with 4 miRNAs and indirectly connected with 21 mRNAs, which indicated that lncRNA Neat1 could play an important role in acute phase of TBI.
Previous studies reported that downregulation of miR-377 could promote angiogenesis and inhibit inflammation to alleviate ischemic brain injury and renal ischemia/reperfusion injury [34, 35]. Feng et al. found that lncRNA ADAMTS9-AS2 could upregulate IGFBP-2 expression via decreasing miR-185-5p expression, further promoting angiogenesis . Moreover, several studies have confirmed that lncRNA Neat1 could sponge miR-185-5p to regulate IGF2 expression or DNMT1/mTOR signaling, promoting cancer progression [37, 38]. In addition, the lncRNA Neat1 was reported to regulate CDK6 and CDK14 via sponging miR107 and further promote tumor growth [39, 40, 41]. Previous study found that over-expressed lncRNA Neat1 could promote axon growth of neurons in vitro, inhibit cell apoptosis and restrict inflammation after TBI . To investigate the downstream mechanism of Neat1, we further analyzed subnetwork of Neat1 (Fig. 3D). For example, in this ceRNA network, an upregulated DEmRNA, myeloid differentiation factor 88 (Myd88), were reported to be overexpressed dramatically in brain tissues after TBI and activate NF-κB as well as promote the production of proinflammatory cytokines . Moreover, Cai et al. found that miR-31-5p/Myd88/ NF-κB pathway plays an important role in blood–brain barrier damage after subarachnoid hemorrhage . In addition, upregulated miR-31-5p was reported to inhibit cell apoptosis . Yang et al. revealed that lncRNA Neat1 was able to modulate inflammatory responses via miR-31-5p/POU2F1 axis . Thus, we suggested that Neat1/miR-31-5p/Myd88 axis may be a crucial pathway in TBI process.
Functional enrichment analyses of TBI ceRNA network
To explore the biological functions of ceRNA network in TBI, we performed the GO and KEGG pathway analyses for mRNAs involved in the network. The top 20 enriched GO terms were shown in Fig. 4A. The terms of signaling pathway regulation was significantly enriched, including “positive regulation of I-kappaB kinase/NF-kappaB signaling”, “positive regulation of JNK cascade” and “positive regulation of ERK1 and ERK2 cascade”. The significant enriched terms of “angiogenesis” and “positive regulation of smooth muscle cell proliferation” indicated that the ceRNA network may be involved in blood vessel reconstruction of cortex post-TBI. Moreover, the ceRNA network was also involved in immune and inflammatory response.
Construction of PPI network and identification of hub genes
To further explore the relationship among the DEmRNA involved in ceRNA network, we constructed PPI networks via the STRING database. The PPI network was shown in Fig. 4B. Furthermore, using Cytoscape App cytoHubba, we identified the hub genes with highest degrees (bigger purple circles), including Myd88, Tnfrsf1b, Cd1d1, Itgam and Irf1. These hub genes could play important roles in TBI ceRNA network.
Validation of DElncRNAs and DEmRNAs using qRT-PCR
To validate the reliability of the RNA-seq data, we randomly selected two up-regulated lncRNAs (lncRNA H19 and lncRNA Mir155hg), two down-regulated lncRNAs (lncRNA Gm36823 and lncRNA C030018K13Rik), two up-regulated mRNAs (Cxcr2 and Mmp12) and two down-regulated mRNAs (Hes5 and P2ry12) that were abundantly expressed and exhibited significant changes. We used qRT-PCR to analyze expression differences in control and injury cortex. The qRT-PCR analysis results were mostly consistent with the RNA-seq data (Fig. 4C–J).
In recent years, an increasing number of studies have revealed that ncRNAs, including lncRNAs and miRNAs, play a key role in TBI [47, 48]. More and more studies have focused on the role of lncRNA-related ceRNA regulatory networks in various diseases, such as cerebral infarction , type 2 diabetes  and cardiac hypertrophy . However, the potential role of the lncRNA-associated ceRNA regulatory network in TBI still remains unclear. In the present study, 9945 lncRNAs and 21,807 mRNAs were identified via RNA-seq analysis. 86 lncRNAs (47 upregulated and 39 downregulated) and 1201 mRNAs (1076 upregulated and 125 downregulated) were detected to dysregulated in the cortex of mice 24 h after TBI. Functional enrichment analyses indicated that dysregulated genes were involved in immune inflammatory processes and cell responses to tumor necrosis factor/interferon-gamma/interferon-beta/interleukin-1. Moreover, ceRNA regulatory network was constructed based on 23 mRNAs, 5 miRNAs and 2 lncRNAs. Function prediction indicated that the network was involved in angiogenesis, immune and inflammatory response, and activation of several signaling pathways. The subnetwork of Neat1 was further analyzed. Furthermore, PPI network was constructed and 5 hub genes were identified. Finally, randomly selected DEmRNAs and DElncRNAs were validated via qRT-PCR. This study could provide a comprehensive perspective on the underlying lncRNA regulatory mechanism in TBI.
The upregulation of LncRNA Neat1 was found to function in early apoptosis. Zhong et al.  explored the relationship between Neat1 and bexarotene in TBI treatment in mice. They found that bexarotene can upregulated Neat1, which further inhibited apoptosis and inflammation, contributing to better motor and cognitive function after TBI. In this study, we found that 5 genes in Neat1 subnetwork were involved in regulation of apoptosis, such as Igfbp3, Myd88, Itgam, Tnfrsf1b and Spp1 [50, 51, 52]. The potential miRNAs between these apoptosis-related genes and Neat1 were also shown in the subnetwork (Fig. 3D). Besides, we revealed that 5 genes in Neat1 subnetwork were involved in regulation of inflammatory response, such as Myd88, Krt16, Spp1, Tnfrsf1b and Sele. Likewise, the potential miRNAs between these inflammatory response-related genes and Neat1 were also shown in the subnetwork (Fig. 3D). Interestingly, miR-31-5p was reported to be able to regulate the angiogenesis of vascular endothelial cells . In addition, previous study found that miR-31-5p was involved in blood–brain barrier repair and Myd88/NF-κB pathway-mediated inflammation after subarachnoid hemorrhage . Furthermore, Myd88 was significantly upregulated after TBI and was able to activate NF-κB and then regulate inflammation during TBI process . Thus, we suggested that Neat1/miR-31-5p/Myd88 axis might play an important role in regulation inflammation cytokines and blood–brain barrier repair in brain tissues after TBI, which needs to be further confirmed in future study.
The expression of lncRNA H19 in brain tissues was upregulated after intracerebral hemorrhage (ICH), which could activate NF-κB and enhance inflammatory responses to aggravate brain edema and neurological injury [54, 55]. The role of lncRNA Mir155hg in brain injury has rarely been reported. Li et al. reported that upregulation of lncRNA Mir155hg was able to promote M1 phenotype macrophage polarization and the release of proinflammatory cytokines in chronic obstructive pulmonary disease . However, the function of lncRNA Gm36823 and lncRNA C030018K13Rik has not been reported so far. Cxcr2, a chemokine receptor on cellular surface, could upregulate and activate microglia in cerebral stoke, Alzheimer’s disease and multiple sclerosis [57,58,59, 60]. Moreover, Cxcr2 might be involved in neutrophil infiltration and subsequent neurodegeneration following TBI [61, 62]. The Hes5 was reported to be downregulated after brain injury [63, 64], while the Mmp12 expression was significantly increased in intracerebral hemorrhage . However, the function of Mmp12 and Hes5 in brain injury needs to further studied. P2ry12 as a microglia marker could be utilized to evaluate the microglial cells. Previous study found that P2ry12 positive microglia were significantly increased after TBI . However, in this study, we found the P2ry12 expression was decreased after TBI. The different expression trends might be attributed to injury severity, brain regions or time points for sampling.
To summary, a lncRNA-associated ceRNA regulatory network of TBI was successfully constructed, which may provide a comprehensive view of the underlying mechanisms of gene regulation and interaction in TBI. Moreover, we proposed that the regulatory network centred on lncRNA Neat1 may play a critical part in TBI. Additionally, we further suggested that the Neat1/miR-31-5p/Myd88 axis might be potential downstream molecular bases for Neat1 to regulate apoptosis, inflammation and blood–brain barrier damage. However, these regulatory axes require further studies to confirm the molecular mechanisms.
Availability of data and materials
The datasets used or analyzed during the current study are available from the corresponding author on reasonable request.
Corrigan JD, Selassie AW, Orman JA. The epidemiology of traumatic brain injury. J Head Trauma Rehabil. 2010;25:72–80.
Feigin VL, Theadom A, Barker-Collo S, Starkey NJ, McPherson K, Kahan M, Dowell A, Brown P, Parag V, Kydd R, Jones K, Jones A, Ameratunga S. Incidence of traumatic brain injury in New Zealand: a population-based study. Lancet Neurol. 2013;12:53–64.
(2010). Traumatic brain injury: time to end the silence. Lancet Neurol 9, 331.
Stoica BA, Faden AI. Cell death mechanisms and modulation in traumatic brain injury. Neurotherapeutics. 2010;7:3–12.
Zweckberger K, Eros C, Zimmermann R, Kim SW, Engel D, Plesnila N. Effect of early and delayed decompressive craniectomy on secondary brain damage after controlled cortical impact in mice. J Neurotrauma. 2006;23:1083–93.
Abou-Abbass H, Bahmad H, Ghandour H, Fares J, Wazzi-Mkahal R, Yacoub B, Darwish H, Mondello S, Harati H, El Sayed MJ, Tamim H, Kobeissy F. Epidemiology and clinical characteristics of traumatic brain injury in Lebanon: a systematic review. Medicine (Baltimore). 2016;95: e5342.
Wilusz JE, Sunwoo H, Spector DL. Long noncoding RNAs: functional surprises from the RNA world. Genes Dev. 2009;23:1494–504.
Caley DP, Pink RC, Trujillano D, Carter DR. Long noncoding RNAs, chromatin, and development. Sci World J. 2010;10:90–102.
Guttman M, Donaghey J, Carey BW, Garber M, Grenier JK, Munson G, Young G, Lucas AB, Ach R, Bruhn L, Yang X, Amit I, Meissner A, Regev A, Rinn JL, Root DE, Lander ES. lincRNAs act in the circuitry controlling pluripotency and differentiation. Nature. 2011;477:295–300.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011;146:353–8.
Pan Y-B, Lu J, Yang B, Lenahan C, Zhang J, Shao A. Construction of competitive endogenous RNA network reveals regulatory role of long non-coding RNAs in intracranial aneurysm. BMC Neurosci. 2021;22:15.
Zou J-B, Chai H-B, Zhang X-F, Guo D-Y, Tai J, Wang Y, Liang Y-L, Wang F, Cheng J-X, Wang J, Shi Y-J. Reconstruction of the lncRNA-miRNA-mRNA network based on competitive endogenous RNA reveal functional lncRNAs in Cerebral Infarction. Sci Rep. 2019;9:12176.
Lin Z, Li X, Zhan X. Construction of competitive endogenous RNA network reveals regulatory role of long non-coding RNAs in type 2 diabetes mellitus. J Cell Mol Med. 2017;21:3204–13.
Pan H, Pan J, Song S, Ji L, Lv H, Yang Z. Identification and development of long non-coding RNA-associated regulatory network in colorectal cancer. J Cell Mol Med. 2019;23:5200–10.
Wang JD, Zhou HS, Tu XX, He Y, Liu QF, Liu Q, Long ZJ. Prediction of competing endogenous RNA coexpression network as prognostic markers in AML. Aging (Albany NY). 2019;11:3333–47.
Long J, Xiong J, Bai Y, Mao J, Lin J, Xu W, Zhang H, Chen S, Zhao H. Construction and investigation of a lncRNA-associated ceRNA regulatory network in cholangiocarcinoma. Front Oncol. 2019;9:649.
Raveendra BL, Swarnkar S, Avchalumov Y, Liu X-A, Grinman E, Badal K, Reich A, Pascal BD, Puthanveettil SV. Long noncoding RNA GM12371 acts as a transcriptional regulator of synapse function. Proc Natl Acad Sci. 2018;115:E10197–205.
Bond AM, Vangompel MJ, Sametsky EA, Clark MF, Savage JC, Disterhoft JF, Kohtz JD. Balanced gene regulation by an embryonic brain ncRNA is critical for adult hippocampal GABA circuitry. Nat Neurosci. 2009;12:1020–7.
Ramos AD, Andersen RE, Liu SJ, Nowakowski TJ, Hong SJ, Gertz C, Salinas RD, Zarabi H, Kriegstein AR, Lim DA. The long noncoding RNA Pnky regulates neuronal differentiation of embryonic and postnatal neural stem cells. Cell Stem Cell. 2015;16:439–47.
Mercer TR, Dinger ME, Sunkin SM, Mehler MF, Mattick JS. Specific expression of long noncoding RNAs in the mouse brain. Proc Natl Acad Sci U S A. 2008;105:716–21.
Wang F, Wang Q, Liu B, Mei L, Ma S, Wang S, Wang R, Zhang Y, Niu C, Xiong Z, Zheng Y, Zhang Z, Shi J, Song X. The long noncoding RNA Synage regulates synapse stability and neuronal function in the cerebellum. Cell Death Differ. 2021;28:2634–50.
Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.
Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019;37:907–15.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The sequence alignment/map format and SAMtools. Bioinformatics (Oxford, England). 2009;25:2078–9.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics (Oxford, England). 2014;30:923–30.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Li J-H, Liu S, Zhou H, Qu L-H, Yang J-H. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein–RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2013;42:D92–7.
Chou C-H, Shrestha S, Yang C-D, Chang N-W, Lin Y-L, Liao K-W, Huang W-C, Sun T-H, Tu S-J, Lee W-H, Chiew M-Y, Tai C-S, Wei T-Y, Tsai T-R, Huang H-T, Wang C-Y, Wu H-Y, Ho S-Y, Chen P-R, Chuang C-H, Hsieh P-J, Wu Y-S, Chen W-L, Li M-J, Wu Y-C, Huang X-Y, Ng FL, Buddhakosai W, Huang P-C, Lan K-C, Huang C-Y, Weng S-L, Cheng Y-N, Liang C, Hsu W-L, Huang H-D. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2017;46:D296–302.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, Santos A, Doncheva NT, Roth A, Bork P, Jensen LJ, von Mering C. The STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessible. Nucleic Acids Res. 2016;45:D362–8.
Redell JB, Liu Y, Dash PK. Traumatic brain injury alters expression of hippocampal microRNAs: potential regulators of multiple pathophysiological processes. J Neurosci Res. 2009;87:1435–48.
Meissner L, Gallozzi M, Balbi M, Schwarzmaier S, Tiedt S, Terpolilli NA, Plesnila N. Temporal profile of microRNA expression in contused cortex after traumatic brain injury in mice. J Neurotrauma. 2016;33:713–20.
Fan Y, Ding S, Sun Y, Zhao B, Pan Y, Wan J. MiR-377 regulates inflammation and angiogenesis in rats after cerebral ischemic injury. J Cell Biochem. 2018;119:327–37.
Liu Z, Yang Q, Wei Q, Chang Y, Qu M, Yu L. The protective effect of miR-377 inhibitor against renal ischemia-reperfusion injury through inhibition of inflammation and oxidative stress via a VEGF-dependent mechanism in mice. Mol Immunol. 2019;106:153–8.
Feng N, Wang Z, Wu Y, Zheng H, Jiang X, Wang Z, Qu F, Zhang Z. ADAMTS9-AS2 promotes angiogenesis of brain microvascular endothelial cells through regulating miR-185–5p/IGFBP-2 axis in ischemic stroke. Mol Neurobiol. 2022. https://doi.org/10.1007/s12035-021-02641-1.
Yu H, Xu A, Wu B, Wang M, Chen Z. Long noncoding RNA NEAT1 promotes progression of glioma as a ceRNA by sponging miR-185-5p to stimulate DNMT1/mTOR signaling. J Cell Physiol. 2021;236:121–30.
Zhuang S-T, Cai Y-J, Liu H-P, Qin Y, Wen J-F. LncRNA NEAT1/miR-185-5p/IGF2 axis regulates the invasion and migration of colon cancer. Mol Genet Genomic Med. 2020;8: e1125.
Wang P, Wu T, Zhou H, Jin Q, He G, Yu H, Xuan L, Wang X, Tian L, Sun Y, Liu M, Qu L. Long noncoding RNA NEAT1 promotes laryngeal squamous cell cancer through regulating miR-107/CDK6 pathway. J Exp Clin Cancer Res. 2016;35:22.
Yang X, Xiao Z, Du X, Huang L, Du G. Silencing of the long non-coding RNA NEAT1 suppresses glioma stem-like properties through modulation of the miR-107/CDK6 pathway. Oncol Rep. 2017;37:555–62.
Zhen Y, Nan Y, Guo S, Zhang L, Li G, Yue S, Liu X. Knockdown of NEAT1 repressed the malignant progression of glioma through sponging miR-107 and inhibiting CDK14. J Cell Physiol. 2019;234:10671–9.
Zhong J, Jiang L, Huang Z, Zhang H, Cheng C, Liu H, He J, Wu J, Darwazeh R, Wu Y, Sun X. The long non-coding RNA Neat1 is an important mediator of the therapeutic effect of bexarotene on traumatic brain injury in mice. Brain Behav Immun. 2017;65:183–94.
Zhu H-T, Bian C, Yuan J-C, Chu W-H, Xiang X, Chen F, Wang C-S, Feng H, Lin J-K. Curcumin attenuates acute inflammatory injury by inhibiting the TLR4/MyD88/NF-κB signaling pathway in experimental traumatic brain injury. J Neuroinflamm. 2014;11:59.
Cai L, Ge B, Xu S, Chen X, Yang H. Up-regulation of circARF3 reduces blood-brain barrier damage in rat subarachnoid hemorrhage model via miR-31-5p/MyD88/NF-κB axis. Aging. 2021;13:21345–63.
Peng H, Wang L, Su Q, Yi K, Du J, Wang Z. MiR-31-5p promotes the cell growth, migration and invasion of colorectal cancer cells by targeting NUMB. Biomed Pharmacother. 2019;109:208–16.
Yang Y, Xue J, Qin L, Zhang J, Liu J, Yu J. LncRNA NEAT1 Promotes Inflammatory Response in Sepsis via the miR-31-5p/POU2F1 Axis. Inflammation. 2021;44:1518–28.
Pan Y-B, Sun Z-L, Feng D-F. The Role of MicroRNA in Traumatic Brain Injury. Neuroscience. 2017;367:189–99.
Li Z, Han K, Zhang D, Chen J, Xu Z, Hou L. The role of long noncoding RNA in traumatic brain injury. Neuropsychiatr Dis Treat. 2019;15:1671–7.
Gao S, Chen Y, Xu B, Yu C, Yue M, Tan X, Zhang J, Feng C, Song C, Ai B, Bai X, Liu Y, Li X, Zhao J, Shi S, Chen X, Su X, Li C. Identification and analysis of a key long non-coding RNAs (lncRNAs)-associated module reveal functional lncRNAs in cardiac hypertrophy. J Cell Mol Med. 2018;22:892–903.
Bolisetty S, Traylor AM, Kim J, Joseph R, Ricart K, Landar A, Agarwal A. Heme oxygenase-1 inhibits renal tubular macroautophagy in acute kidney injury. J Am Soc Nephrol. 2010;21:1702–12.
Li S, Zhao Y, He X, Kim TH, Kuharsky DK, Rabinowich H, Chen J, Du C, Yin XM. Relief of extrinsic pathway inhibition by the Bid-dependent mitochondrial release of Smac in Fas-mediated hepatocyte apoptosis. J Biol Chem. 2002;277:26912–20.
Shibue T, Takeda K, Oda E, Tanaka H, Murasawa H, Takaoka A, Morishita Y, Akira S, Taniguchi T, Tanaka N. Integral role of Noxa in p53-mediated apoptotic response. Genes Dev. 2003;17:2233–8.
Ye J, Zhu J, Chen H, Qian J, Zhang L, Wan Z, Chen F, Sun S, Li W, Luo C. A novel lncRNA-LINC01116 regulates tumorigenesis of glioma by targeting VEGFA. Int J Cancer. 2020;146:248–61.
Chen B, Wang H, Lv C, Mao C, Cui Y. Long non-coding RNA H19 protects against intracerebral hemorrhage injuries via regulating microRNA-106b-5p/acyl-CoA synthetase long chain family member 4 axis. Bioengineered. 2021;12:4004–15.
Mao S, Huang H, Chen X. lncRNA H19 aggravates brain injury in rats following experimental intracerebral hemorrhage via NF-κB pathway. Comput Math Methods Med. 2022;2022:3017312.
Li N, Liu Y, Cai J. LncRNA MIR155HG regulates M1/M2 macrophage polarization in chronic obstructive pulmonary disease. Biomed Pharmacother. 2019;117: 109015.
Filipovic R, Jakovcevski I, Zecevic N. GRO-α and CXCR2 in the human fetal brain and multiple sclerosis lesions. Dev Neurosci. 2003;25:279–90.
Popivanova BK, Koike K, Tonchev AB, Ishida Y, Kondo T, Ogawa S, Mukaida N, Inoue M, Yamashima T. Accumulation of microglial cells expressing ELR motif-positive CXC chemokines and their receptor CXCR2 in monkey hippocampus after ischemia-reperfusion. Brain Res. 2003;970:195–204.
Watson K, Fan G-H. Macrophage inflammatory protein 2 inhibits β-amyloid peptide (1–42)-mediated hippocampal neuronal apoptosis through activation of mitogen-activated protein kinase and phosphatidylinositol 3-kinase signaling pathways. Mol Pharmacol. 2005;67:757–65.
Bye N, Habgood MD, Callaway JK, Malakooti N, Potter A, Kossmann T, Morganti-Kossmann MC. Transient neuroprotection by minocycline following traumatic brain injury is associated with attenuated microglial activation but no changes in cell apoptosis or neutrophil infiltration. Exp Neurol. 2007;204:220–33.
Semple BD, Bye N, Ziebell JM, Morganti-Kossmann MC. Deficiency of the chemokine receptor CXCR2 attenuates neutrophil infiltration and cortical damage following closed head injury. Neurobiol Dis. 2010;40:394–403.
Semple BD, Kossmann T, Morganti-Kossmann MC. Role of chemokines in CNS health and pathology: a focus on the CCL2/CCR2 and CXCL8/CXCR2 networks. J Cereb Blood Flow Metab. 2010;30:459–73.
Xiao X, Bai P, Cao S, Jiang Y, Liang W, Wang T, Luo X, Guan Q, Gao L, Zhang L. Integrated bioinformatics analysis for the identification of key molecules and pathways in the hippocampus of rats after traumatic brain injury. Neurochem Res. 2020;45:928–39.
Zhang Y-M, Chen S-X, Dai Q-F, Jiang S-T, Chen A-L, Tang C-Z, Zhang Y-Q. Effect of acupuncture on the notch signaling pathway in rats with brain injury. Chin J Integr Med. 2018;24:537–44.
Power C, Henry S, Del Bigio MR, Larsen PH, Corbett D, Imai Y, Yong VW, Peeling J. Intracerebral hemorrhage induces macrophage activation and matrix metalloproteinases. Ann Neurol. 2003;53:731–42.
Saglam E, Zırh S, Aktas CC, Muftuoglu SF, Bilginer B. Papaverine provides neuroprotection by suppressing neuroinflammation and apoptosis in the traumatic brain injury via RAGE-NF-κB pathway. J Neuroimmunol. 2021;352: 577476.
We would like to express our thanks to the National Library of Medicine for giving us the privilege to access the raw data of various GEO series.
This work was funded by the Medical Science and Technology Project of Zhejiang province (2018ZD022).
Ethics approval and consent to participate
The animal study was reviewed and approved by the Animal Ethics Committee of the Zhejiang University.
Consent for publication
The authors declare no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original version of this article was revised: an error in the affiliations assignment for Yiyu Sun was corrected
About this article
Cite this article
Wang, S., Sun, Y., Hu, S. et al. Construction of a lncRNA-associated competing endogenous RNA regulatory network after traumatic brain injury in mouse. Mol Brain 15, 40 (2022). https://doi.org/10.1186/s13041-022-00925-8
- Bioinformatics analysis
- Traumatic brain injury