Skip to main content
Download PDF

Construction of a lncRNA-associated competing endogenous RNA regulatory network after traumatic brain injury in mouse

Molecular Brain volume 15, Article number: 40 (2022) Cite this article

A Correction to this article was published on 13 July 2022

This article has been updated

Abstract

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.

Introduction

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 [3]. TBI can be divided into initial injury and secondary brain injury [4]. 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 [5]. 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 [6]. 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 [7]. 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 [11], cerebral infarction [12], type 2 diabetes [13] 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

Animals

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 [22], 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) [23]. The generated SAM (Sequence Alignment/Map) files were converted to the BAM (Binary Alignment/Map) format by using samtools (version 1.10) [24]. The BAM files were converted to counts files by using featureCounts in subread software (version 2.0.0) [25]. 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 [26], 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/) [27]. 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 [28]. Then, miRNA-targeted mRNAs were identified using miRTarBase, a highly reliable miRNA reference database [29]. 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/) [30].

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 [31]. 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.

Statistical analysis

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).

Results

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.

Fig. 1
figure 1

LncRNA profile based on RNA-seq data. A Distribution of all identified lncRNAs in TBI samples. B Volcano plot of DElncRNAs in mice cortex 24 h post-TBI compared with control cortex. Red points represent upregulated DElncRNAs and blue points represent downregulated DElncRNAs in TBI group. The genesymbols of ten upregulated lncRNAs and ten downregulated lncRNAs with the most significant expression differences are shown in the plot. C Hierarchical cluster heatmaps of DElncRNAs. Each row represents an RNA, and each column represents a sample. Red indicates relatively high expression, and blue indicates relatively low expression. D Distribution of DElncRNAs in TBI, showing upregulated (gray) and downregulated (black) lncRNAs in each chromosome (chr)

Full size image

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.

Fig. 2
figure 2

mRNA profile based on RNA-seq data. A Distribution of all identified mRNA in TBI samples. B Volcano plot of DEmRNAs in mice cortex 24 h post-TBI compared with control cortex. Red points represent upregulated DEmRNAs and blue points represent downregulated DEmRNAs in TBI group. The genesymbols of ten upregulated mRNAs and ten downregulated mRNAs with the most significant expression differences are shown in the plot. C Hierarchical cluster heatmaps of DEmRNAs. Each row represents an RNA, and each column represents a sample. Red indicates relatively high expression, and blue indicates relatively low expression. D Distribution of DEmRNAs in TBI, showing upregulated (gray) and downregulated (black) mRNAs in each chromosome (chr)

Full size image

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.

Fig. 3
figure 3

A The top 20 enriched terms in the GO biological process analysis. B The top 20 enriched pathways in the KEGG pathway analysis. C, D ceRNA network and subnetwork related to TBI. C lncRNA-associated ceRNA regulatory network related to TBI. D Subnetwork of lncRNA Neat1. Blue squares represent lncRNAs, tan triangles represent miRNAs, and red circles represent mRNAs. E–I Validation of expression alternations of miRNAs after TBI via qRT-PCR. All the experiments were performed and analyzed in triplicate. (*P < 0.05)

Full size image

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 [36]. 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 [42]. 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 [43]. 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 [44]. In addition, upregulated miR-31-5p was reported to inhibit cell apoptosis [45]. Yang et al. revealed that lncRNA Neat1 was able to modulate inflammatory responses via miR-31-5p/POU2F1 axis [46]. 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.

Fig. 4
figure 4

A The top 20 enriched terms in the GO biological process analysis. B PPI network of mRNAs in ceRNA network. Circles represent protein-coding genes, and edges represent the interactions between two proteins. The purple circles represent top 10 hub genes with highest degrees. C–J Validation of expression differences of DElncRNAs and DEmRNAs via qRT-PCR. C–D P2ry12 and Hes5 expression significantly decrease in injury cortex 24 h post-TBI compared to control cortex. E–F Cxcr2 and Mmp12 expression significantly increase in injury cortex 24 h post-TBI compared to control cortex. G–H lncRNA C030018K13Rik and Gm36823 expression significantly decrease in injury cortex 24 h post-TBI compared to control cortex. I–J lncRNA H19 and Mir155hg expression significantly increase in injury cortex 24 h post-TBI compared to control cortex. All the experiments were performed and analyzed in triplicate. (*P < 0.05, **P < 0.01, ***P < 0.001)

Full size image

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).

Discussion

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 [12], type 2 diabetes [13] and cardiac hypertrophy [49]. 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. [42] 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 [53]. 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 [44]. Furthermore, Myd88 was significantly upregulated after TBI and was able to activate NF-κB and then regulate inflammation during TBI process [43]. 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 [56]. 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 [65]. 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 [66]. 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.

Change history

References

  1. Corrigan JD, Selassie AW, Orman JA. The epidemiology of traumatic brain injury. J Head Trauma Rehabil. 2010;25:72–80.

    Article  PubMed  Google Scholar 

  2. 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.

    Article  PubMed  Google Scholar 

  3. (2010). Traumatic brain injury: time to end the silence. Lancet Neurol 9, 331.

  4. Stoica BA, Faden AI. Cell death mechanisms and modulation in traumatic brain injury. Neurotherapeutics. 2010;7:3–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. 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.

    Article  PubMed  Google Scholar 

  6. 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.

    Article  Google Scholar 

  7. Wilusz JE, Sunwoo H, Spector DL. Long noncoding RNAs: functional surprises from the RNA world. Genes Dev. 2009;23:1494–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Caley DP, Pink RC, Trujillano D, Carter DR. Long noncoding RNAs, chromatin, and development. Sci World J. 2010;10:90–102.

    Article  CAS  Google Scholar 

  9. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. 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.

    Article  CAS  Google Scholar 

  16. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  17. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. 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.

    Article  CAS  Google Scholar 

  25. 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.

    Article  CAS  Google Scholar 

  26. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. 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.

    Article  CAS  Google Scholar 

  28. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. 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.

    Article  PubMed Central  CAS  Google Scholar 

  30. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. 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.

    Article  PubMed  Google Scholar 

  34. 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.

    Article  CAS  PubMed  Google Scholar 

  35. 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.

    Article  CAS  PubMed  Google Scholar 

  36. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  37. 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.

    Article  CAS  PubMed  Google Scholar 

  38. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  39. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. 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.

    Article  PubMed  Google Scholar 

  41. 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.

    Article  CAS  PubMed  Google Scholar 

  42. 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.

    Article  CAS  PubMed  Google Scholar 

  43. 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.

    Article  CAS  Google Scholar 

  44. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. 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.

    Article  CAS  PubMed  Google Scholar 

  46. 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.

    Article  CAS  PubMed  Google Scholar 

  47. Pan Y-B, Sun Z-L, Feng D-F. The Role of MicroRNA in Traumatic Brain Injury. Neuroscience. 2017;367:189–99.

    Article  CAS  PubMed  Google Scholar 

  48. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. 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.

    Google Scholar 

  50. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. 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.

    Article  CAS  PubMed  Google Scholar 

  52. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. 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.

    Article  CAS  PubMed  Google Scholar 

  54. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. 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.

    PubMed  PubMed Central  Google Scholar 

  56. Li N, Liu Y, Cai J. LncRNA MIR155HG regulates M1/M2 macrophage polarization in chronic obstructive pulmonary disease. Biomed Pharmacother. 2019;117: 109015.

    Article  CAS  PubMed  Google Scholar 

  57. Filipovic R, Jakovcevski I, Zecevic N. GRO-α and CXCR2 in the human fetal brain and multiple sclerosis lesions. Dev Neurosci. 2003;25:279–90.

    Article  CAS  PubMed  Google Scholar 

  58. 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.

    Article  CAS  PubMed  Google Scholar 

  59. 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.

    Article  CAS  PubMed  Google Scholar 

  60. 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.

    Article  CAS  PubMed  Google Scholar 

  61. 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.

    Article  CAS  PubMed  Google Scholar 

  62. 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.

    Article  CAS  PubMed  Google Scholar 

  63. 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.

    Article  CAS  PubMed  Google Scholar 

  64. 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.

    Article  CAS  PubMed  Google Scholar 

  65. 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.

    Article  CAS  PubMed  Google Scholar 

  66. 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.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

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.

Funding

This work was funded by the Medical Science and Technology Project of Zhejiang province (2018ZD022).

Author information

Author notes

  1. Siqi Wang, Yiyu Sun and Shaobo Hu contributed equally to this manuscript

Authors and Affiliations

  1. Department of Nuclear Medicine, Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Hangzhou, 310000, Zhejiang, China

    Siqi Wang & Cen Lou

  2. Department of Plastic and Reconstructive Surgery, Shanghai 9th People’s Hospital, Shanghai Jiao Tong University, 639 Zhi Zao Ju Road, Shanghai, 200011, China

    Yiyu Sun

  3. Department of Neurosurgery, Li Hui Li Hospital of Medical Center of Ningbo, Ningbo, Zhejiang, China

    Shaobo Hu

  4. Department of Neurosurgery, Second Affiliated Hospital, School of Medicine, Zhejiang University, NO. 88 Jiefang Rd, Hangzhou, 310009, Zhejiang, China

    Yuan-Bo Pan

Authors
  1. Siqi Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yiyu Sun
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Shaobo Hu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Cen Lou
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yuan-Bo Pan
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Y-BP designed the study. SW, Y-YS, SH, Y-BP and CL analyzed and interpreted the data. SW wrote the article. All authors reviewed and approved the final manuscript.

Corresponding authors

Correspondence to Cen Lou or Yuan-Bo Pan.

Ethics declarations

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

Not applicable.

Competing interests

The authors declare no conflicts of interest.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1.

PCR primes used in this study.

Additional file 2.

The list of DElncRNAs and DEmRNAs.

Additional file 3.

Interactions between lncRNA and miRNA in the ceRNA network.

Additional file 4.

Interactions between miRNA and mRNA in the ceRNA network.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13041-022-00925-8

Keywords

  • Bioinformatics analysis
  • Traumatic brain injury
  • ceRNA
  • lncRNA
  • Neat1
  • Myd88

Molecular Brain

ISSN: 1756-6606