Int J Med Sci 2020; 17(4):437-448. doi:10.7150/ijms.39667
Plasma Exosome-derived MicroRNAs as Novel Biomarkers of Traumatic Brain Injury in Rats
1. Emergency Center, Zhongnan Hospital of Wuhan University, Wuhan 430071, China.
2. Hubei Clinical Research Center for Emergency and Resuscitation, Zhongnan Hospital of Wuhan University, Wuhan 430071, China.
Wang P, Ma H, Zhang Y, Zeng R, Yu J, Liu R, Jin X, Zhao Y. Plasma Exosome-derived MicroRNAs as Novel Biomarkers of Traumatic Brain Injury in Rats. Int J Med Sci 2020; 17(4):437-448. doi:10.7150/ijms.39667. Available from http://www.medsci.org/v17p0437.htm
Traumatic brain injury (TBI) is a widespread central nervous system (CNS) condition and a leading cause of death, disability, and long-term disability including seizures and emotional and behavioral issues. To date, applicable diagnostic biomarkers have not been elucidated. MicroRNAs (miRNAs) are enriched and stable in exosomes in plasma. Therefore, we speculated that miRNAs in plasma exosomes might serve as novel biomarkers for TBI diagnosis and are also involved in the pathogenesis of TBI. In this study, we first isolated exosomes from peripheral blood plasma in rats with TBI and then investigated the alterations in miRNA expression in exosomes by high-throughput RNA sequencing. As a result, we identified 50 significantly differentially expressed miRNAs, including 31 upregulated and 19 downregulated miRNAs. Then, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that the most highly correlated pathways that were identified were the MAPK signaling pathway, regulation of actin cytoskeleton, Rap1 signaling pathway and Ras signaling pathway. This study provides novel perspectives on miRNAs in peripheral blood plasma exosomes, which not only could be used as biomarkers of TBI diagnosis but could also be manipulated as therapeutic targets of TBI.
Keywords: Traumatic brain injury, exosomes, microRNA, plasma, biomarker, rat.
Due to the increasing incidence of traffic accidents and falls, an increasing number of people every year seek medical care for traumatic brain injury (TBI). Among these individuals, 75% or more are considered to have a mild or moderate TBI . Although severe TBI is one of the main causes of death or disability in patients, mild or moderate TBI is a leading cause of long-term disability including seizures and emotional and behavioural issues . The damage to the cranium and the intracranial contents caused by a TBI can be divided into the primary injury and delayed secondary injury, which is more profound and includes axonal shearing , neuroinflammation [4,5], neurochemical changes [6, 7], brain edema , vascular injury , cell apoptosis  and mitochondrial dysfunction .
Currently, strategies for diagnosing and treating TBI in patients are limited . Imaging tests (computed tomography, CT or magnetic resonance imaging, MRI) are generally applied to diagnose and classify TBI in clinical practice but are limited because of low sensitivity in patients with mild injuries [13, 14]. Cerebrospinal fluid (CSF) or blood tests are useful in the assessment of the severity of the injury and subsequent recovery of TBI patients by detecting protein or miRNA components in the peripheral blood . For instance, protein tests that can detect cytokines, chemokines, tau, calcium-binding protein S100B, glial fibrillary acidic protein (GFAP), neuron-specific enolase (NSE), ubiquitin carboxy-terminal hydrolase L1 (UCH-L1), and myelin basic protein (MBP) are candidates for classification and prognostication of TBI patients . However, these tests have detection sensitivity problems due to the low concentration of the substances (fM-pM). MiRNAs, which are small (19-23 nt) endogenous RNA molecules, are post-transcriptional regulators of gene expression and protein synthesis, bind to target sites in the noncoding region of mRNA, are enriched and stable in exosomes and are involved in the onset of many diseases. Microarray and sequencing are the most often used methods for miRNA quantitation. We chose the latter method because sequencing yields more precise accuracy and the ability to investigate the novel structure of miRNA. Furthermore, microRNAs (miRNA) could be indicators of the pathology of TBI. Serum or plasma miRNAs have been reported as specific and sensitive biomarkers of many CNS conditions [16, 17], including TBI [18-21]. For example, a subset of ten unique miRNAs was shown to have the ability to diagnose and distinguish between mild/moderate TBI and severe TBI [22, 23]. However, one constraint is that miRNAs can easily dissolve or degrade in plasma samples.
Recently, several studies have suggested that extracellular vehicles (EVs) might play a vital role in TBI pathogenesis . Exosomes are EVs with diameters that typically range from 40-150 nm. These particles have a phospholipid bilayer that carries proteins, metabolites, lipids, mRNAs, and miRNAs that are released from cells and could be an important medium for cell-cell or cell-internal environment communication . After trauma, brain cell-derived exosomes containing intracranial substances cross the blood-brain barrier (BBB) and are transferred to the peripheral blood. Therefore, they might serve as diagnostic and prognostic biomarkers. MiRNAs have been shown to be stable in exosomes and EVs . Alteration of miRNA expression has also shown to occur with disease onset and progression. Therefore, we conjectured that miRNAs in exosomes and EVs in the peripheral blood plasma after TBI might play a vital role in regulating gene expression, protein synthesis and specific signaling pathways, which could constitute a good biomarker for diagnosis.
The present study aimed to identify vital miRNAs, provide novel targets for further study, and evaluate the alterations in miRNA expression in exosomes in the plasma after TBI by high-throughput RNA sequencing. Bioinformatics analysis was used to predict the function of differentially expressed miRNAs by gene ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway analysis. Pathway relationship networks was constructed based on differentially expressed miRNAs to explore the interactions among them. In addition, 11 miRNAs were randomly selected to confirm the reliability of the microRNA sequencing results by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). This study revealed several novel potential roles of these differentially expressed miRNAs in plasma exosomes from the peripheral blood after TBI.
Materials and Methods
Experiments were performed using male adult SPF-grade SD rats (12--16 weeks old and 300-400 g, purchased from the Animal Center of Wuhan University) that were acclimated at least one week before surgery. All animal experiment protocols were approved by the Animal Experiment Center and Ethics Committee of Zhongnan Hospital of Wuhan University and followed the National Institutes of Health Guide for the Care and Use of Laboratory Animals.
TBI model and experimental groups
In this study, we used the weight-drop (WD) method as described in previous studies to induce moderate TBI in rats [26, 27]. Briefly, the rats were intraperitoneally anaesthetized with 1% pentobarbital (30 mg/kg). After anaesthesia, the rats were fixed, the hair was shaved, and the skin was sterilized and cut. The bregma was exposed, and a 5-mm-diameter hole (2.0 mm posterior from the bregma and 2.0 mm right of the sagittal suture) was drilled into the skull; then, the dura mater was exposed. A 40 g -weight was used for a vertical strike at 25 cm, resulting in moderate TBI. The sham group received surgery but did not undergo the WD.
After 24 hours, all rats were anaesthetized and peripheral blood was collected in a tube containing EDTA as an anticoagulant. The samples were stored at 4°C for a short time and immediately centrifuged at 4°C and 3000 rpm for 10 min. The supernatant was carefully collected and transferred to a new tube without disturbing the intermediate buffy coat layer, and then the sample was passed through a 0.8 μm filter (SLAA033SB, Millex, Millipore) to remove additional cellular fragments and cell debris. Then the samples were collected for further steps. In total, 14 samples (seven from the TBI group and seven from the sham group) were used to isolate exosomes and extract miRNAs. Among them, six samples (three from the sham group and three from the TBI group) were used for miRNA sequencing, and eight samples (four from the sham group and four from the TBI group) were chosen for qRT-PCR validation.
Isolation and detection of exosomes
Exosomes derived from plasma were isolated from the plasma using an exoEasy Maxi Kit (Qiagen, 76064) following the manufacturer's instructions. Exosomes were observed under a transmission electron microscope (TEM; FEI Tecnai G2 spirit, The Czech Republic), stained with phosphotungstic acid and detected at a voltage of 80kV. A nanoparticle tracking analyser (NTA) (ZetaView) was used to determine the size distribution of the plasma-derived exosomes.
The proteins were extracted from plasma exosomes of both groups. Western blotting was performed to detect typical exosomal markers using WES, an automated capillary-based size sorting system (ProteinSimple, San Jose CA). The data were analysed using Compass software (ProteinSimple, San Jose CA). The primary antibodies used included CD63 (SBI, EXOAB-CD63A-1, 1:50), HSP70 (SBI, EXOAB-Hsp70A-1, 1:50), and Transferrin (Abcam, ab82411, 1:50), which was used as a reference control.
microRNA sequencing and subsequent bioinformatics analyses
The miRNA was isolated from exosome samples using an exoRNeasy Serum/Plasma Maxi Kit (Qiagen, 77064) following the manufacturer's instructions. The sequence library was determined by an Agilent 2100 bioanalyser. Sequencing and analysis of small RNAs were conducted by OE biotech Co. Ltd. (Shanghai, China). The basic reads were converted into raw data by base calling and were obtained from an Illumine Hiseq4000 sequencer. Low-quality reads were filtered, and reads with 5' primer contaminants or poly (A) sequences were removed. Reads without a 3'adapter and insert tag and reads shorter than 15 nt or longer than 41 nt were filtered from the raw data, and the clean reads were obtained.
For the primary analysis, the length distribution of the clean sequences in the reference genome was determined. Non-coding RNAs were annotated as rRNAs, tRNAs, small nuclear RNAs (snRNAs), and small nucleolar RNAs (snoRNAs). These RNAs were aligned and then subjected to BLAST  searches against the Rfam v.10.1 (http://www.sanger.ac.uk/software/Rfam) [29, 30] and GenBank databases (http://www.ncbi.nlm.nih.gov/genbank/). Known miRNAs were identified by alignment against the miRBase v.21 database (http://www.mirbase.org/) , and the expression patterns of known miRNAs in different samples were analysed. Unannotated small RNAs were analysed by miRDeep2  to predict novel miRNAs. The corresponding miRNA star sequence was also identified based on the hairpin structure of the pre-miRNA and the miRbase database.
Differentially expressed miRNAs were identified with a threshold p-value < 0.05. The p-value was calculated with the DEG algorithm in the R package. The targets of differentially expressed miRNAs were predicted using Targetscan, miRDB, and Diana software, and the interaction was selected as the final decision.
GO enrichment (http://www.geneontology.org) and KEGG (http://www.genome.jp/kegg) pathway enrichment analyses of differentially expressed miRNA-target -genes were performed using R based on the hypergeometric distribution.
To verify the accuracy of the miRNA sequencing data, qRT-PCR was performed as a validation method. The target RNA and housekeeping genes of each sample were subjected to RT-PCR reactions on a CFX Connect Real-Time PCR System (Bio-Rad; CFX Maestro 1.0 software) with TB Green Premix Ex Taq II (Tli RNaseH Plus) (RR820A, TaKaRa, Japan).
RNA was reverse transcribed to synthesize cDNA using a Mir-X miRNA First Strand Synthesis Kit (638313, TaKaRa, Osaka, Japan); The PCR and premier design were performed according to the manufacturer's instructions. Briefly, RNA samples were mixed with mRQ buffer and mRQ Enzyme, preferably in a thermal cycler and incubated for 1 h at 37°C then the reaction was terminated at 85°C for 5min to inactivate the enzymes. When finished, a 9-fold volume of RNase-free water was added.
All procedures were performed as follows: 95°C, 30 sec; 45 PCR cycles (95°C, 5 sec; 60°C, 30 sec [plate read]); 95°C, 10 sec. The temperature was slowly increased from 65°C to 95°C to establish the melting curve of the PCR product (automatic instrument, ramp rate 0.05°C /s). The primers of the randomly selected miRNAs and the internal parameters of U6 are shown in Table 3.
All statistical analyses were performed in SPSS 23.0. Data are expressed as the mean ± SD. Student's t-test was used to evaluate differences between the two groups. A p-value<0.05 was considered to indicate a significant difference.
Detection of exosomes
Transmission electron microscope (TEM) analysis, nanoparticle tracking analysis (NTA) and western blot analysis were used to identify the exosomes that were obtained. In TEM images, the exosomes were irregular spheres ranging from 30-150 nm in diameter with a clearly defined and relatively intact membrane (Figure 1A). The NTA results showed a continuous range of exosomes, with a peak of 140.3nm (Figure 1B), suggesting a significant range of exosomes. The western blot analysis showed that the exosomes of both groups expressed the typical exosomal markers such as CD63 and HSP70 (Figure 1C).
Differential expression of miRNAs in exosomes from the plasma after TBI
The miRNAs were sequenced in three TBI and three sham exosome samples from plasma. The miRNA-seq reads of each sample are shown in Table 1. The heatmap shows sample-to-sample distances, and samples with high similarity are preferentially clustered together (Figure 2A). We identified 50 significantly differentially expressed miRNAs that met the criteria of an expression fold change ≥ 2.0 and p value ≤ 0.05. Of these, 31 were upregulated and 19 were downregulated (Table 2). A heat map of 50 differentially expressed miRNAs was generated to illustrate the distinguishable miRNA expression profile of the samples (Figure 2B). A volcano plot was created according to the means of the expression values of the two groups of samples and the threshold of the expression fold change was 2.0 (Figure 2C).
|Sample||Raw reads||Clean reads||miRNA nunber|
miRNA, micro ribonucleic acid; TBI, traumatic brain injury.
Detection of exosomes. (A) Representative transmission electron microscope (TEM) image of exosomes. TEM image of exosomes with negative staining to enhance the view of membrane structures (scale bar=100nm). Black arrow showed a typical exosome structure. (B) Representative nanoparticle tracking (NTA) analysis result of exosomes. The peak scale was measures as 140.3nm. (C) Representative western blot (WES) images of exosomes. Exosomes were probed for CD63, HSP70, and transferrin. CD, cluster of differentiation; HSP, heat shock protein.(Click on the image to enlarge.)
Correlation of samples, Differential expression of micro ribonucleic acids (miRNAs) in exosomes from plasma after traumatic brain injury (TBI). (A) The heatmap showed sample to sample distances, samples with high similarity were clustered together preferentially. (B) Heat plot of 50 differentially expressed miRNAs. Each column represents one sample; each row represents one probe set. The dendrogram on the top reveals the sample clustering; the dendrogram on the left reveals the gene clustering. (C) Volcano plot of expressions of miRNAs. Red spots were up-regulated miRNAs with significantly differential expressions, green spots were down-regulated, and gray spots were miRNAs with non-differential expressions. X-axis represented for log2 Fold Change; Y-axis represented for -log10 P value.(Click on the image to enlarge.)
Top 50 of differentially expressed miRNAs.
FC, fold change.
Primers sequences used in this study.
|miRNA||Primer type||Primer Sequence (5'-3')|
The reverse sequence and U6 sequence were provided in Mir-X miRNA First Strand Synthesis Kit (638313, TaKaRa, Osaka, Japan) as identified by previous studies.
GO analysis of differentially expressed miRNA targeted genes
The genes were subjected to GO function analysis to annotate and speculate the function of these miRNAs (Figure 3). GO analysis is divided into three parts: biological process (BP), cell component (CC), and molecular function (MF). GO analysis of the BP showed that differentially expressed miRNAs were significantly associated with the following terms: regulation of histamine secretion by mast cell; positive regulation of cell proliferation by VEGF-activated platelet derived growth factor receptor signaling pathway; signal transduction by trans-phosphorylation; metanephric glomerular capillary formation; polyamine catabolic process; positive regulation of transcription by RNA polymerase II; positive regulation of transcription, DNA-templated; negative regulation of transcription by RNA polymerase II; transcription by RNA polymerase II; and regulation of transcription by RNA polymerase II. The significant GO CC terms of differentially expressed miRNAs included the followings: cytoplasm, nucleus, cytosol, membrane, neuronal cell body, nucleoplasm, Golgi apparatus, cell-cell junction, postsynaptic density, and dendritic spine. For MF, differentially expressed miRNA were associated with the following terms: protein binding; protein kinase binding; protein serine/threonine kinase activity; transcriptional activator activity, RNA polymerase II proximal promoter sequence-specific DNA binding; RNA polymerase II proximal promoter sequence-specific DNA binding; protein domain specific binding; ATP binding; DNA binding; kinase activity; and platelet-derived growth factor binding.
Gene ontology (GO) analysis of differentially expressed micro ribonucleic acids (miRNAs) with top 30 gene of -log10 P value in each GO terms.(Click on the image to enlarge.)
KEGG pathway analysis of differentially expressed miRNA targeted genes
A KEGG pathway dot plot was constructed to show the significant enrichment pathway with the top 20 enrichment score values. KEGG pathway analysis predicted the pathways affected by the variation of miRNAs in exosomes from the plasma after TBI (Figure 4A). The top 20 pathways included MAPK signaling pathway, Proteoglycans in cancer, Tight junction, Focal adhesion, Mitophagy - animal, Regulation of actin cytoskeleton, Insulin resistance, Rap1 signaling pathway, FoxO signaling pathway, EGFR tyrosine kinase inhibitor resistance, Gap junction, Endocrine resistance, Axon guidance, Glioma, Neurotrophin signaling pathway, Longevity regulating pathway - mammal, Ras signaling pathway, Non-small cell lung cancer, Endocytosis, Pathways in cancer. Of all of these pathways, Mitophagy - animal seemed to have the highest enrichment score.
Similar to GO classification, KEGG classification counted the number or components of proteins or genes at the following function levels: cellular processes (transport and catabolism, cellular community-eukaryotes, cell motility, cell growth and death), environmental information processing (signaling molecules and interaction, signal transduction, membrane transport), genetic information processing (translation, transcription, replication and repair, fold, sorting and degradation), human diseases (substance dependence, neurodegenerative diseases, infectious diseases, immune diseases, endocrine and metabolic diseases, drug resistance, cardiovascular diseases, cancers), metabolism (Xenobiotics biodegradation and metabolism, nucleotide metabolism, metabolism of terpenoids and polyketides, metabolism of other amino acids, metabolism of cofactors and vitamins, lipid metabolism, glycan biosynthesis and metabolism, energy metabolism, carbohydrate metabolism, biosynthesis of other secondary metabolites, amino acid metabolism), and organismal systems (sensory system, nervous system, immune system, excretory system environmental adaptation, endocrine system, digestive system, development, circulatory system, aging) (Figure 4B).
Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of differentially expressed micro ribonucleic acids (miRNAs) (A) KEGG pathways of top20 enrichment score. (B) Distribution of KEGG Level2 of different miRNAs.(Click on the image to enlarge.)
A pathway relationship network of the top 20 significant pathways of differentially expressed miRNAs was constructed (Figure 5), showing that the MAPK signaling pathway, Regulation of actin cytoskeleton, Rap1 signaling pathway, and Ras signaling pathway may have a regulatory effect on upstream levels.
Validation of the accuracy of miRNA-seq data by qRT-PCR
Eleven differentially expressed miRNAs were randomly selected for validation. The target gene expression levels of these RNAs were normalized to their internal control, U6. Comparisons of the expression levels of these miRNAs by sequencing and qRT-PCR are shown in Figure 6. The expression levels of these miRNAs were calculated using the 2-δδCT method. The data are expressed as the means ± SD. There were significant differences in nine of all the selected miRNAs (p<0.05). Among them, the variation tendencies of miR-124-3p, miR-142-3p, miR-145-3p, miR-374-5p, miR-532-5p, miR-29b-3p and miR-106b-5p were in accordance with the micro sequencing; however, the variation tendencies of miR-92a-3p and miR-451-5p were not in accordance with the sequencing results. Moreover, although no significant differences were found (p > 0.05), the other miRNAs (miR-181c-3p and miR-9a-3p) exhibited concordant results with the miRNA sequencing results (data not shown). In total, the objective miRNA validation rate was 7/9, showing that the miRNA expression profiles were reliable.
Pathway relation network analysis of the top 20 enrichment of miRNAs. The network was established based on the results of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and KEGG database search.(Click on the image to enlarge.)
The altered expression profile of miRNAs in exosomes from the plasma of peripheral blood after TBI in rats was studied for the first time by high-throughput whole transcriptome sequencing and subsequent bioinformatics analysis. We made several novel observations that improve the understanding of post-TBI molecular and intermolecular interactions as well as the key signaling pathways in which they participate. First, we demonstrated that the components of exosomes from plasma change after TBI. Exosomes were isolated and then confirmed by TEM, NTA and WES. The change in biological material in exosomes from plasma can reflect the physiological and pathological processes after TBI and may be a potential target of TBI therapy.
In this study, we found a total of 50 miRNAs in exosomes that were differentially expressed after TBI, including 31 upregulated and 19 downregulated miRNAs. These altered miRNAs might participate in the progression of TBI. Several previous studies have demonstrated the involvement of various spectrums of miRNAs from different sources in the TBI process, including miRNAs from cerebrospinal fluid (CSF) , serum , brain tissue , and saliva . Compared to these samples, CSF or brain tissue is difficult to obtain in an emergency room (ER), and preparation of serum usually takes hours, while plasma can be acquired relatively easily along with blood tests in ER. In the present study, miRNA samples were acquired from exosomes in plasma. We found a different spectrum of miRNAs compared to those found in previous studies, but this difference coincided with a previous study , which found that measurement of exosome miRNAs cannot be replaced by measurement of plasma miRNAs and vice versa. We suggest that this discrepancy may be due to sample disparity, miRNAs in plasma, plasma-exosomes, or the fact that other biofluids might be regulated by different mechanisms.
Quantitative reverse transcription polymerase chain reaction validation results of 9 selected micro ribonucleic acids. *, p<0.05; **, p<0.01, ***, p<0.001.(Click on the image to enlarge.)
We found some miRNAs that were reported in previous studies. For example, miR124-3p has long been studied and reported to have a role in regulating hepatocellular carcinoma , breast cancer [35, 36], acute lymphoblastic leukaemia  and non-small cell lung cancer . It has also been reported in CNS diseases such as Alzheimer's disease  and glioma [40, 41]. Previous studies  have shown that miR-124-3p is upregulated after TBI in brain tissue in microglial exosomes from the acute to the chronic phase of TBI, and increases in miR-124-3p in microglia promote an anti-inflammatory response in both in vitro and in vivo experiments. MiR-124-3p may also be a chronic regulator of gene expression after brain injury , therefore, miR-124-3p is a promising therapeutic target for interventions in neuroinflammation after TBI. Additionally, miR-142-3p, which has also long been studied, has been reported in several tumours [44-46] and other diseases [47, 48]. Its expression in rat hippocampus after TBI has been reported , indicating that it may be a sensitive and informative biomarker for forensic assessment of TBI. As to miR-106-5p and miR-532-5p, they have been reported in ischemic stroke as biomarkers, and they may serve a novel role in the pathogenesis . Furthermore, miR-181c-3p has been shown in previous studies to be regulated by several circRNAs after TBI in mice, including chr17:39846653-39847022+, chr17:39846695-39847059+, and chr17:39846692-39847021+ . Although involvement of this mechanism has not been reported in TBI, it has been reported in other fields [52-54], indicating that it may have a therapeutic role in TBI treatment or assessment.
Moreover, to gain insight into the potential function of the differentially expressed miRNAs, KEGG pathway analysis was performed, which revealed that the differentially expressed miRNAs participated in several biological pathways. Of these pathways, the MAPK signaling pathway, Ras signaling pathway, and Rap1 signaling pathway are related to environmental information processing and signal transduction. Based on pathway analysis, we constructed a pathway relationship network, which indicated that the most important upstream pathway was MAPK signaling pathway and Ras signaling pathway, and Rap1 signaling pathway were downstream pathways also involved in this network. Many studies have indicated the role and the potential therapy of activation of RAS or RAP1 in traumatic brain injury [55-57]. There were also several studies about the actin cytoskeleton, which showed the mechanism in TBI procedure [58, 59]. As the upstream pathway, MAPK signaling pathway may serve as another important role.
The MAPK signaling pathway, which is a group of serine/threonine protein kinases involved in the cellular signal transduction , which regulated many cellular functions, affects a large number of human diseases. Previous studies have demonstrated the effect of MAPK signals on tumor metastasis [61, 62], osteogenesis [63, 64], osteoarthritis [65, 66], or certain inflammatory diseases . For example, it can regulate the proliferation and differentiation of brain cells . Furthermore, this pathway was also found to be involved in neuroinflammation in traumatic brain injury; according to a study that showed that p38α/MAPK regulates microglial to TBI mice . In other diseases, a research on osteoarthritis showed a suppressed miR-92-3p expression up activation of MAPK in chondrocytes ; miR-221-5p regulates MAPK/ERK signaling pathway in prostate cancer cells, which suggested miR-221-5p plays an important role in prostate cancer progression ; a research in glioma showed the upregulation of miR-124-3p induced cellular apoptosis via MAPK activation, and promote cell autophagy . These combined evidences showed several miRNAs in our studies are according with MAPK signaling pathway.
Last, we predicted the miRNAs that target the mRNA network according to the database and bioinformatics criteria. Briefly, the network showed that rno-miR-181c-3p, rno-miR-28-3p and rno-miR-451-5p regulated few mRNAs. In contrast, rno-miR-124-3p, rno-miR-106b-5p, rno-miR-29b-3p, and rno-miR-9a-5p, regulated a large number of mRNAs and seemed to be more important. According to our network and previous studies, these mRNAs and the related mRNAs may merit further research in the TBI process.
In the present study, for the first time, we identified a series of differentially expressed miRNAs in exosomes obtained from the plasma after WD-induced TBI in rats, which may be related to the physiological and pathological processes that occur after TBI. We also predicted the potential roles of pathways and the interactions of miRNAs with mRNAs. The present study demonstrated that miRNA in plasma exosomes are altered after TBI, which suggests that they may serve as biomarkers for the diagnosis or merit consideration as therapeutic targets in treating TBI.
TBI: Traumatic Brain Injury; CNS: Central Nervous System; miRNA: Micro Ribonucleic Acid; CD: Cluster of Differentiation; HSP: Hat Sock Protein; GO: Gene Oncology; KEGG: Kyoto Encyclopedia of Gene and Genomes; BP: Biological Process; CC: Cellular Component; MF: Molecular Function.
This study was supported by Emergency Diagnostic and Therapeutic Center of Central China (Wuhan, Hubei, China), Science, Technology and Innovation Seed Fund of Zhongnan Hospital of Wuhan University (cxpy2017039, WJ2019H010/znpy2018043), and the Fundamental Research Funds for the Central Universities (2042019kf0149).
Xiaoqing Jin and Yan Zhao conceived of and designed the research plans; Pengcheng Wang performed most of the experiments and analyzed the data; Haoli Ma, Yuxian Zhang, Rong Zeng, Jiangtao Yu, and Ruining Liu provided technical assistance to Pengcheng Wang; Pengcheng Wang wrote the article with contributions from all the authors; Xiaoqing Jin and Yan Zhao supervised and supported the writing.
The authors have declared that no competing interest exists.
1. Voss JD, Connolly J, Schwab KA, Scher AI. Update on the Epidemiology of Concussion/Mild Traumatic Brain Injury. Curr Pain Headache Rep. 2015;19:32
2. Albrecht JS, Peters ME, Smith GS, Rao V. Anxiety and Posttraumatic Stress Disorder Among Medicare Beneficiaries After Traumatic Brain Injury. J Head Trauma Rehabil. 2017;32:178-84
3. van Eijck MM, Schoonman GG, van der Naalt J, de Vries J, Roks G. Diffuse axonal injury after traumatic brain injury is a prognostic factor for functional outcome: a systematic review and meta-analysis. Brain Inj. 2018;32:395-402
4. Yoo D, Magsam AW, Kelly AM, Stayton PS, Kievit FM, Convertine AJ. Core-Cross-Linked Nanoparticles Reduce Neuroinflammation and Improve Outcome in a Mouse Model of Traumatic Brain Injury. ACS Nano. 2017;11:8600-11
5. White TE, Ford GD, Surles-Zeigler MC, Gates AS, Laplaca MC, Ford BD. Gene expression patterns following unilateral traumatic brain injury reveals a local pro-inflammatory and remote anti-inflammatory response. BMC Genomics. 2013;14:282
6. Shahim P, Tegner Y, Gustafsson B, Gren M, Arlig J, Olsson M. et al. Neurochemical Aftermath of Repetitive Mild Traumatic Brain Injury. JAMA Neurol. 2016;73:1308-15
7. Birnie M, Morrison R, Camara R, Strauss KI. Temporal changes of cytochrome P450 (Cyp) and eicosanoid-related gene expression in the rat brain after traumatic brain injury. BMC Genomics. 2013;14:303
8. Tucker B, Aston J, Dines M, Caraman E, Yacyshyn M, McCarthy M. et al. Early Brain Edema is a Predictor of In-Hospital Mortality in Traumatic Brain Injury. J Emerg Med. 2017;53:18-29
9. Tao XG, Shi JH, Hao SY, Chen XT, Liu BY. Protective Effects of Calpain Inhibition on Neurovascular Unit Injury through Downregulating Nuclear Factor-kappaB-related Inflammation during Traumatic Brain Injury in Mice. Chin Med J (Engl). 2017;130:187-98
10. White TE, Surles-Zeigler MC, Ford GD, Gates AS, Davids B, Distel T. et al. Bilateral gene interaction hierarchy analysis of the cell death gene response emphasizes the significance of cell cycle genes following unilateral traumatic brain injury. BMC Genomics. 2016;17:130
11. Gonzalez-Dominguez R. Medium-chain Fatty Acids as Biomarkers of Mitochondrial Dysfunction in Traumatic Brain Injury. EBioMedicine. 2016;12:8-9
12. Galgano M, Toshkezi G, Qiu X, Russell T, Chin L, Zhao LR. Traumatic Brain Injury: Current Treatment Strategies and Future Endeavors. Cell Transplant. 2017;26:1118-30
13. Mutch CA, Talbott JF, Gean A. Imaging Evaluation of Acute Traumatic Brain Injury. Neurosurg Clin N Am. 2016;27:409-39
14. Posti JP, Takala RSK, Lagerstedt L, Dickens AM, Hossain I, Mohammadian M. et al. Correlation of Blood Biomarkers and Biomarker Panels with Traumatic Findings on Computed Tomography after Traumatic Brain Injury. J Neurotrauma. 2019
15. Karnati HK, Garcia JH, Tweedie D, Becker RE, Kapogiannis D, Greig NH. Neuronal Enriched Extracellular Vesicle Proteins as Biomarkers for Traumatic Brain Injury. J Neurotrauma. 2019;36:975-87
16. Jin XF, Wu N, Wang L, Li J. Circulating microRNAs: a novel class of potential biomarkers for diagnosing and prognosing central nervous system diseases. Cellular and molecular neurobiology. 2013;33:601-13
17. Kichukova TM, Popov NT, Ivanov HY, Vachev TI. Circulating microRNAs as a Novel Class of Potential Diagnostic Biomarkers in Neuropsychiatric Disorders. Folia medica. 2015;57:159-72
18. Mitra B, Rau TF, Surendran N, Brennan JH, Thaveenthiran P, Sorich E. et al. Plasma micro-RNA biomarkers for diagnosis and prognosis after traumatic brain injury: A pilot study. J Clin Neurosci. 2017;38:37-42
19. Yang T, Song J, Bu X, Wang C, Wu J, Cai J. et al. Elevated serum miR-93, miR-191, and miR-499 are noninvasive biomarkers for the presence and progression of traumatic brain injury. J Neurochem. 2016;137:122-9
20. Di Pietro V, Ragusa M, Davies D, Su Z, Hazeldine J, Lazzarino G. et al. MicroRNAs as Novel Biomarkers for the Diagnosis and Prognosis of Mild and Severe Traumatic Brain Injury. J Neurotrauma. 2017;34:1948-56
21. Taheri S, Tanriverdi F, Zararsiz G, Elbuken G, Ulutabanca H, Karaca Z. et al. Circulating MicroRNAs as Potential Biomarkers for Traumatic Brain Injury-Induced Hypopituitarism. J Neurotrauma. 2016;33:1818-25
22. Bhomia M, Balakathiresan NS, Wang KK, Papa L, Maheshwari RK. A Panel of Serum MiRNA Biomarkers for the Diagnosis of Severe to Mild Traumatic Brain Injury in Humans. Sci Rep. 2016;6:28148
23. Chandran R, Sharma A, Bhomia M, Balakathiresan NS, Knollmann-Ritschel BE, Maheshwari RK. Differential expression of microRNAs in the brains of mice subjected to increasing grade of mild traumatic brain injury. Brain Inj. 2017;31:106-19
24. Kumar A, Stoica BA, Loane DJ, Yang M, Abulwerdi G, Khan N. et al. Microglial-derived microparticles mediate neuroinflammation after traumatic brain injury. J Neuroinflammation. 2017;14:47
25. Jeppesen DK, Fenix AM, Franklin JL, Higginbotham JN, Zhang Q, Zimmerman LJ. et al. Reassessment of Exosome Composition. Cell. 2019;177:428-45 e18
26. Zhao L, Liu Q, Ma S, Zhang Y, Liang P. TPEN Attenuates Neural Autophagy Induced by Synaptically-released Zinc Translocation and Improves Histological Outcomes after Traumatic Brain Injury in Rats. Ann Clin Lab Sci. 2018;48:446-52
27. Buchele F, Morawska MM, Schreglmann SR, Penner M, Muser M, Baumann CR. et al. Novel Rat Model of Weight Drop-Induced Closed Diffuse Traumatic Brain Injury Compatible with Electrophysiological Recordings of Vigilance States. J Neurotrauma. 2016;33:1171-80
28. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403-10
29. Griffiths-Jones S, Bateman A, Marshall M, Khanna A, Eddy SR. Rfam: an RNA family database. Nucleic Acids Res. 2003;31:439-41
30. Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ. miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008;36:D154-8
31. Friedlander MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40:37-52
32. Di Pietro V, Porto E, Ragusa M, Barbagallo C, Davies D, Forcione M. et al. Salivary MicroRNAs: Diagnostic Markers of Mild Traumatic Brain Injury in Contact-Sport. Front Mol Neurosci. 2018;11:290
33. Xie JX, Fan X, Drummond CA, Majumder R, Xie Y, Chen T. et al. MicroRNA profiling in kidney disease: Plasma versus plasma-derived exosomes. Gene. 2017;627:1-8
34. Hu XX, Feng J, Huang XW, Lu PZ, Wang ZX, Dai HQ. et al. Histone deacetylases up-regulate C/EBPalpha expression through reduction of miR-124-3p and miR-25 in hepatocellular carcinoma. Biochem Biophys Res Commun. 2019
35. Yan G, Li Y, Zhan L, Sun S, Yuan J, Wang T. et al. Decreased miR-124-3p promoted breast cancer proliferation and metastasis by targeting MGAT5. Am J Cancer Res. 2019;9:585-96
36. Duan H, Li X, Wang C, Hao P, Song W, Li M. et al. Functional hyaluronate collagen scaffolds induce NSCs differentiation into functional neurons in repairing the traumatic brain injury. Acta Biomater. 2016;45:182-95
37. Yang T, Jin X, Lan J, Wang W. Long non-coding RNA SNHG16 has Tumor suppressing effect in acute lymphoblastic leukemia by inverse interaction on hsa-miR-124-3p. IUBMB Life. 2019;71:134-42
38. Tang LX, Chen GH, Li H, He P, Zhang Y, Xu XW. Long non-coding RNA OGFRP1 regulates LYPD3 expression by sponging miR-124-3p and promotes non-small cell lung cancer progression. Biochem Biophys Res Commun. 2018;505:578-85
39. Zhou Y, Deng J, Chu X, Zhao Y, Guo Y. Role of Post-Transcriptional Control of Calpain by miR-124-3p in the Development of Alzheimer's Disease. J Alzheimers Dis. 2019;67:571-81
40. Yang JX, Liu B, Yang BY, Meng Q. Long non-coding RNA homeobox (HOX) A11-AS promotes malignant progression of glioma by targeting miR-124-3p. Neoplasma. 2018;65:505-14
41. Hu D, Zhang Y. Circular RNA HIPK3 promotes glioma progression by binding to miR-124-3p. Gene. 2019;690:81-9
42. Huang S, Ge X, Yu J, Han Z, Yin Z, Li Y. et al. Increased miR-124-3p in microglial exosomes following traumatic brain injury inhibits neuronal inflammation and contributes to neurite outgrowth via their transfer into neurons. FASEB J. 2018;32:512-28
43. Vuokila N, Lukasiuk K, Bot AM, van Vliet EA, Aronica E, Pitkanen A. et al. miR-124-3p is a chronic regulator of gene expression after brain injury. Cell Mol Life Sci. 2018;75:4557-81
44. Jia XP, Meng LL, Fang JC, Wang HW, Chen J, Zhou J. et al. Aberrant Expression of miR-142-3p and its Target Gene HMGA1 and FZD7 in Breast Cancer and its Clinical Significance. Clin Lab. 2018;64:915-21
45. Jia AY, Castillo-Martin M, Domingo-Domenech J, Bonal DM, Sanchez-Carbayo M, Silva JM. et al. A common MicroRNA signature consisting of miR-133a, miR-139-3p, and miR-142-3p clusters bladder carcinoma in situ with normal umbrella cells. Am J Pathol. 2013;182:1171-9
46. Ghanbari R, Mosakhani N, Asadi J, Nouraee N, Mowla SJ, Yazdani Y. et al. Downregulation of Plasma MiR-142-3p and MiR-26a-5p in Patients With Colorectal Carcinoma. Iran J Cancer Prev. 2015;8:e2329
47. Makino K, Jinnin M, Kajihara I, Honda N, Sakai K, Masuguchi S. et al. Circulating miR-142-3p levels in patients with systemic sclerosis. Clin Exp Dermatol. 2012;37:34-9
48. Kleinsteuber K, Heesch K, Schattling S, Kohns M, Sander-Julch C, Walzl G. et al. Decreased expression of miR-21, miR-26a, miR-29a, and miR-142-3p in CD4(+) T cells and peripheral blood from tuberculosis patients. PLoS One. 2013;8:e61609
49. Sun TY, Chen XR, Liu ZL, Zhao LL, Jiang YX, Qu GQ. et al. Expression profiling of microRNAs in hippocampus of rats following traumatic brain injury. J Huazhong Univ Sci Technolog Med Sci. 2014;34:548-53
50. Li P, Teng F, Gao F, Zhang M, Wu J, Zhang C. Identification of circulating microRNAs as potential biomarkers for detecting acute ischemic stroke. Cell Mol Neurobiol. 2015;35:433-47
51. Zhao RT, Zhou J, Dong XL, Bi CW, Jiang RC, Dong JF. et al. Circular Ribonucleic Acid Expression Alteration in Exosomes from the Brain Extracellular Space after Traumatic Brain Injury in Mice. J Neurotrauma. 2018;35:2056-66
52. Rajarajan D, Selvarajan S, Charan Raja MR, Kar Mahapatra S, Kasiappan R. Genome-wide analysis reveals miR-3184-5p and miR-181c-3p as a critical regulator for adipocytes-associated breast cancer. J Cell Physiol. 2019
53. Shen X, Li Y, Sun G, Guo D, Bai X. miR-181c-3p and -5p promotes high-glucose-induced dysfunction in human umbilical vein endothelial cells by regulating leukemia inhibitory factor. Int J Biol Macromol. 2018;115:509-17
54. Chen S, Ainiwaer B, Qing S, Liu T, Ma Z, Shi Y. et al. Expression levels of miR-181c-3p and miR-5692b in esophageal cancer and their clinical significance. Zhonghua Bing Li Xue Za Zhi. 2015;44:905-9
55. Xiong Y, Mahmood A, Chopp M. Remodeling dendritic spines for treatment of traumatic brain injury. Neural Regen Res. 2019;14:1477-80
56. Dai X, Yi M, Wang D, Chen Y, Xu X. Changqin NO. 1 inhibits neuronal apoptosis via suppressing GAS5 expression in a traumatic brain injury mice model. Biol Chem. 2019;400:753-63
57. Shi Y, Yi C, Li X, Wang J, Zhou F, Chen X. Overexpression of Mitofusin2 decreased the reactive astrocytes proliferation in vitro induced by oxygen-glucose deprivation/reoxygenation. Neurosci Lett. 2017;639:68-73
58. Berrout J, Jin M, O'Neil RG. Critical role of TRPP2 and TRPC1 channels in stretch-induced injury of blood-brain barrier endothelial cells. Brain Res. 2012;1436:1-12
59. Zhang X, Li Z, Zhang Q, Chen L, Huang X, Zhang Y. et al. Mechanisms Underlying H2O2-Evoked Carbonyl Modification of Cytoskeletal Protein and Axon Injury in PC-12 Cells. Cell Physiol Biochem. 2018;48:1088-98
60. Seger R, Krebs EG. The MAPK signaling cascade. FASEB J. 1995;9:726-35
61. Stephanie M. J. Fliedner TE, Nikoletta K. Lendvai,Uma Shankavaram,Svenja Nolting,Robert Wesley,Abdel G. Elkahloun,Hendrik Ungefroren,Angela Oldoerp,Gary Lampert,Hendrik Lehnert,Henri Timmers, Karel Pacak. Anti-Cancer Potential of MAPK Pathway Inhibition in Paragangliomas-Effect of Different Statins on Mouse Pheochromocytoma Cells. plos one. 2014:9
62. Ge C, Yang Q, Zhao G, Yu H, Kirkwood KL, Franceschi RT. Interactions between extracellular signal-regulated kinase 1/2 and P38 MAP Kinase pathways in the control of Runx2 phosphorylation and transcriptional activity. J Bone Miner Res. 2012;27:538-51
63. Li W, Li G, Zhang Y, Wei S, Song M, Wang W. et al. Role of P2 x 7 receptor in the differentiation of bone marrow stromal cells into osteoblasts and adipocytes. Experimental cell research. 2015;339:367-79
64. Harrison EB, Hochfelder CG, Lamberty BG, Meays BM, Morsey BM, Kelso ML. et al. Traumatic brain injury increases levels of miR-21 in extracellular vesicles: implications for neuroinflammation. FEBS Open Bio. 2016;6:835-46
65. Wang P, Mao Z, Pan Q, Lu R, Huang X, Shang X. et al. Histone deacetylase-4 and histone deacetylase-8 regulate interleukin-1beta-induced cartilage catabolic degradation through MAPK/JNK and ERK pathways. International journal of molecular medicine. 2018;41:2117-27
66. Mao Z, Wang P, Pan Q, Huang X, Zhang R, Shang X. et al. Pravastatin alleviates interleukin 1beta-induced cartilage degradation by restoring impaired autophagy associated with MAPK pathway inhibition. International immunopharmacology. 2018;64:308-18
67. Ismail HM, Yamamoto K, Vincent TL, Nagase H, Troeberg L, Saklatvala J. Interleukin-1 Acts via the JNK-2 Signaling Pathway to Induce Aggrecan Degradation by Human Chondrocytes. Arthritis & rheumatology. 2015;67:1826-36
68. Uhlen P, Fritz N, Smedler E, Malmersjo S, Kanatani S. Calcium signaling in neocortical development. Dev Neurobiol. 2015;75:360-8
69. Bachstetter AD, Rowe RK, Kaneko M, Goulding D, Lifshitz J, Van Eldik LJ. The p38alpha MAPK regulates microglial responsiveness to diffuse traumatic brain injury. J Neurosci. 2013;33:6143-53
70. Mao G, Wu P, Zhang Z, Zhang Z, Liao W, Li Y. et al. MicroRNA-92a-3p Regulates Aggrecanase-1 and Aggrecanase-2 Expression in Chondrogenesis and IL-1beta-Induced Catabolism in Human Articular Chondrocytes. Cell Physiol Biochem. 2017;44:38-52
71. Shao N, Ma G, Zhang J, Zhu W. miR-221-5p enhances cell proliferation and metastasis through post-transcriptional regulation of SOCS1 in human prostate cancer. BMC Urol. 2018;18:14
72. You Y, Wang R, Shao N, Zhi F, Yang Y. Luteolin suppresses tumor proliferation through inducing apoptosis and autophagy via MAPK activation in glioma. Onco Targets Ther. 2019;12:2383-96
Corresponding authors: Xiaoqing Jin, MD; E-mail: redjinedu.cn. Yan Zhao, MD, PhD; E-mail: doctoryanzhaoedu.cn. Tel.: +86-27-67811471. Fax: +86-27-67811472.