Research Paper Volume 11, Issue 19 pp 8313—8328
LncRNAs down-regulate Myh1, Casr, and Mis18a expression in the Substantia Nigra of aged male rats
- 1 Department of Neurobiology, Hebei Medical University, Hebei Province, Shijiazhuang 050017, China
- 2 Department of Human Anatomy, Hebei Medical University, Hebei Province, Shijiazhuang 050017, China
- 3 Hebei Laboratory Animal Center, Hebei Medical University, Hebei Province, Shijiazhuang 050017, China
- 4 Grade 2015 Eight-year Clinical Medicine Program, School of Basic Medical Sciences, Hebei Medical University, Hebei Province, Shijiazhuang 050017, China
- 5 Neuroscience Research Center, Hebei Medical University, Hebei Province, Shijiazhuang 050017, China
- 6 Hebei Key Laboratory of Forensic Medicine, Department of Forensic Medicine, Hebei Medical University, Hebei Province, Shijiazhuang 050017, China
Received: May 24, 2019 Accepted: September 21, 2019 Published: October 2, 2019https://doi.org/10.18632/aging.102321
How to Cite
Copyright © 2019 Zhang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
In this study, we used high-throughput RNA sequencing to identify mRNAs, long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) that are differentially expressed in the Substantia Nigra (SN) of aged and young rats. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were used to perform functional annotation of mRNAs that were either differentially expressed themselves (DEMs), targeted by differentially expressed lncRNAs (DELs), or the parents of differentially expressed circRNAs (DECs). A total of 112 DEMs, 163 DELs, and 98 DECs were found in the SN of aged rats. The down-regulated lncRNA NONRATT010417.2 targeted the down-regulated mRNA Myh1, while the down-regulated lncRNA NONRATT015586.2 and the up-regulated lncRNAs NONRATT000490.2 and NONRATT007029.2 all targeted the down-regulated mRNAs Casr and Mis18a. Western blots and RT-qPCR revealed that Myh1, Casr, and Mis18a protein and mRNA expression were significantly reduced in aged rats compared to young rats. This study improves our understanding of the transcriptional alterations underlying aging-related changes in the SN and provides a foundation for future studies of associated molecular mechanisms.
Aging is associated with the accumulation of deficits in cognition and motor control and with increased risk of neurodegenerative diseases, such as Parkinson’s disease [1–3]. The nigrostriatal dopaminergic pathway consists of projections from dopaminergic cell bodies in the Substantia Nigra (SN) to the corpus striatum, including the caudate, putamen, and globus pallidus . The dopaminergic neurons of the SN are vulnerable to age-related degenerative processes . In normal aging, the number and dopamine content of dopaminergic neurons in the SN progressively decreases, which leads to aging-related behavioral deficits [6, 7]. The dopaminergic neurons of the SN are vulnerable to oxidative stress  and inflammatory attack . Reduced dopamine activity during normal aging is also linked to age-related mitochondrial DNA damage . The mechanisms responsible for dysfunction and degeneration of the SN dopamine neurons in normal aging are complex and are not yet fully understood.
Age-related transcriptome changes may contribute to physiological changes in the dopaminergic neurons of the SN during aging. High-throughput RNA sequencing technologies have helped demonstrate that less than 3% of human transcripts are protein-coding RNAs; the rest are non-protein-coding RNAs, including long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) [11, 12]. Changes in protein expression contribute to changes in cellular function that characterize the aging process. It is therefore necessary to examine the expression of age-associated genes. Increasing evidence indicates that lncRNAs and circRNAs play crucial roles in various biological processes and diseases [13–19]. LncRNAs modulate gene expression via transcriptional, post-transcriptional, translational, and post-translational regulation [15, 20]. CircRNAs can also play important roles during and after transcription by functioning as microRNA (miRNA) sponges, regulating splicing, and interacting with RNA-binding proteins [21–23]. It is therefore important to examine changes in lncRNAs and circRNAs in addition to messenger RNAs (mRNAs) in the SN during normal aging.
In this study, we used high-throughput RNA sequencing to investigate mRNA, lncRNA, and circRNAs expression profiles. We identified RNAs that are differentially expressed in the SNs of 6-month-old rats and 24-month-old rats. We then identified Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with the differentially expressed mRNAs (DEMs), the mRNAs targeted by differentially expressed lncRNAs (DELs), and the parental mRNAs of the differentially expressed circRNAs (DECs). Finally, we evaluated the functional annotations of the DEMs targeted by DELs or that were parents of DECs to identify putative regulatory roles using bioinformatics approaches. Our findings help elucidate the molecular mechanisms by which mRNAs, lncRNAs, and circRNAs contribute to the cognitive and motor impairments as well as the functional deficits observed in SN dopaminergic neurons that occur during normal aging.
Age-dependent DEMs, DELs, and DECs in the SN
A total of 112 DEMs (Supplementary Table 1), 163 DELs (Supplementary Table 2), and 98 DECs (Supplementary Table 3) were found in the SN of 24-month old (24Mon) rats comparted to six-month old (6Mon) rats. Of these, 56 DEMs (50.0% of all DEMs), 80 DELs (49.1% of all DELs), and 52 DECs (53.1% of all DECs) were up-regulated, while 56 DEMs (50.0%), 83 DELs (50.9%), and 46 DECs (46.9%) were down-regulated, in the SN of 24Mon rats. Scatter diagrams (A-C), volcano diagrams (D-F), and hierarchical clustering analysis (G-I) showing differences in SN mRNA, lncRNA, and circRNA expression between 24Mon and 6Mon rats are shown in Figure 1. The ten DEMs for which expression changed the most in the SN of 24Mon rats compared to 6Mon rats are listed in Table 1.
Figure 1. Scatter diagrams (A–C), volcano diagrams (D–F), and hierarchical clustering analysis (G–I) of relative mRNA (A, D, G), lncRNA (B, E, H), and circRNA (C, F, I) expression signals in the SN of 6Mon and 24Mon rats.
Table 1. Top 10 up- and down-regulated DEMs in the SN of 24Mon rats.
|Gene id||Gene name||Description||Fold change||P-value||Regulation|
|ENSRNOG00000047706||RT1-CE16||RT1 class I, locus CE16||Inf||2.95E-06||DOWN|
|ENSRNOG00000001865||Serpind1||serpin peptidase inhibitor, clade D (heparin cofactor), member 1||Inf||0.00018||DOWN|
|ENSRNOG00000021555||Mis18a||MIS18 kinetochore protein A||-4.75492||0.000454||DOWN|
|ENSRNOG00000060523||Rps2-ps2||ribosomal protein S2, pseudogene 2||-3.75752||0.000683||DOWN|
|ENSRNOG00000049695||Myh1||myosin, heavy chain 1, skeletal muscle, adult||-3.37537||3.84E-06||DOWN|
|ENSRNOG00000020953||Ms4a7||membrane-spanning 4-domains, subfamily A, member 7||-2.88917||0.000289||DOWN|
|ENSRNOG00000025946||Igf2bp2||insulin-like growth factor 2 mRNA binding protein 2||-2.40051||0.000306||DOWN|
|ENSRNOG00000013426||Mrgprf||MAS-related GPR, member F||-2.13363||2.38E-07||DOWN|
|ENSRNOG00000037167||Rtp3||receptor (chemosensory) transporter protein 3||Inf||0.000422||UP|
|ENSRNOG00000020400||Mef2b||myocyte enhancer factor 2B||4.379482||0.000155||UP|
|ENSRNOG00000050325||Emr4||EGF-like module containing, mucin-like, hormone receptor-like sequence 4||4.295064||0.000351||UP|
|ENSRNOG00000002031||Naa11||N(alpha)-acetyltransferase 11, NatA catalytic subunit||3.278133||0.000379||UP|
|ENSRNOG00000042220||Vsx1||visual system homeobox 1||3.162037||0.000119||UP|
|ENSRNOG00000007664||Tnfrsf13c||tumor necrosis factor receptor superfamily, member 13c||2.702566||0.000239||UP|
|ENSRNOG00000032246||Acsm3||acyl-CoA synthetase medium-chain family member 3||2.653876||0.000266||UP|
|ENSRNOG00000000451||RT1-Ba||RT1 class II, locus Ba||2.639609||1.30E-07||UP|
The 163 DELs were divided among six biotype categories: three (1.8%) were bidirectional, eight (4.9%) were exonic-antisense, 85 (52.1%) were exonic-sense, 49 (30.1%) were intergenic, two (1.2%) were intronic-antisense, and 16 (9.8%) were intronic-sense (Supplementary Table 4 and Figure 2A). The 98 DECs were divided among three categories: 77 (78.6%) were exon, 18 (18.4%) were intergenic-region, and three (3.1%) were intron (Supplementary Table 5 and Figure 2B).
Functional annotation of DEMs
A total of 112 DEMs were obviously responsive to 23 biological processes (BP), 14 cellular components (CC), and eight molecular functions (MF) in the SN of 24Mon rats (Figure 3A and Supplementary Table 6). These DEMs were associated with 290 functional enrichments (Supplementary Table 7). The top 30 GO enrichments are shown in Figure 3B. KEGG classifications indicated that these DEMs were responsive to four cellular processes and to two environmental information processing, one genetic information processing, four metabolism, and five organismal systems (Figure 3C and Supplementary Table 8). The DEMs were significantly enriched in 77 KEGG pathways (Supplementary Table 9). The top 30 KEGG pathway enrichments are shown in Figure 3D.
Functional annotation of mRNAs targeted by DELs
A total of 737 targeting relationships were identified between DELs and mRNAs (Supplementary Table 10). The mRNAs targeted by DELs were obviously responsive to 25 BPs, 17 CCs, and 12 MFs in the SN of 24Mon rats (Figure 4A and Supplementary Table 11) and were associated with 858 functional enrichments (Supplementary Table 12). The top 30 GO enrichments are shown in Figure 4B. KEGG classifications indicated that the mRNAs targeted by DELs were responsive to four cellular processes and to three environmental information processing, four genetic information processing, ten metabolism, and eight organismal systems (Figure 4C and Supplementary Table 13). The mRNAs targeted by DELs were significantly enriched in 170 KEGG pathways (Supplementary Table 14). The top 30 KEGG pathway enrichments are shown in Figure 4D.
Functional annotation of mRNAs that were parents of DECs
A total of 67 parental relationships were identified between DECs and mRNAs (Supplementary Table 15). The parental mRNAs of DECs were obviously responsive to 21 BPs, 15 CCs, and ten MFs in the SN of 24Mon rats (Figure 5A and Supplementary Table 16) and were associated with 258 functional enrichments (Supplementary Table 17). The top 30 GO enrichments are shown in Figure 5B. KEGG classifications indicate that these parental mRNAs are responsive to three cellular processes and to two environmental information processing, three genetic information processing, six metabolism, and five organismal systems (Figure 5C and Supplementary Table 18). The parental mRNAs were significantly enriched in 60 KEGG pathways (Supplementary Table 19). The top 30 KEGG pathway enrichments are shown in Figure 5D.
Analysis of DEMs targeted by DELs and DEMs that are parents of DECs
Among the 163 DELs, four (NONRATT010417.2, NONRATT015586.2, NONRATT000490.2, and NONRATT007029.2) targeted three DEMs (Myh1, Casr, and Mis18a) (Table 2). None of the parental mRNAs of the 98 DECs were differentially expressed. In the SN of 24Mon rats, down-regulated lncRNA NONRATT010417.2 targeted the down-regulated mRNA Myh1 (ENSRNOG00000049695), while the down-regulated lncRNA NONRATT015586.2 and the up-regulated lncRNAs NONRATT000490.2 and NONRATT007029.2 targeted the down-regulated mRNAs Casr (ENSRNOG00000002265) and Mis18a (ENSRNOG00000021555) (Table 2). GO analysis indicated that these three DEMs (Myh1, Casr, and Mis18a) were obviously responsive to 14 BPs, ten CCs, and four MFs in the SN of 24Mon rats (Figure 6A and Supplementary Table 20) and were not significantly associated with any GO functional enrichments. KEGG pathway analysis revealed that Myh1 is responsive to cellular processes and significantly enriched in the tight junction (rno04530). Myh1 was related to cell polarity (Figure 6B). Casr and Mis18a were neither responsive to any categories nor enriched in any pathways in the KEGG analysis.
Table 2. DELs target DEMs in the SN of 24Mon rats.
|LncRNA_id||Fold change||Pvalue||Qvalue||Regulation||Target gene_id||Target gene name||Fold change||Pvalue||Qvalue||Target gene Regulation|
Myh1, Casr, and Mis18a protein and mRNA expression in the SN of 6Mon and 24Mon rats
Western blot analysis was used to measure levels of Myh1, Casr, and Mis18a proteins extracted from the SN (Figure 7A). Expression of all three proteins was significantly reduced in 24Mon rats compared to 6Mon rats (Figure 7B). RT-qPCR revealed that Myh1, Casr, and Mis18a mRNA expression were significantly reduced in 24Mon rats compared to 6Mon rats (Figure 7C).
In this study, we identified a total of 112 mRNAs that are differentially expressed in the SNs of 24- and six-month old rats. Of these DEMs, 56 were down-regulated and 56 were up-regulated in the SN of aged rats. GO analysis indicated that the DEMs were obviously responsive to 23 BPs, 14 CCs, and eight MFs and were associated with 290 functional enrichments. KEGG pathway analysis revealed that the DEMs were responsive to four cellular processes and to two environmental information processing, one genetic information processing, four metabolism, and five organismal systems, and were significantly enriched in 77 KEGG pathways. Overall, our results suggest that a large number of genes and pathways could play pivotal roles in normal aging processes in the SN.
While the DEMs identified here may be particularly important for understanding aging-related dysfunction and degeneration of SN dopamine neurons, their roles in the SN are poorly understood. Compared to young rats, the ten most strongly down-regulated mRNAs in aged rats were RT1-CE16, AABR07042599.1, Serpind1, Mis18a, Rps2-ps2, Myh1, Ms4a7, Igf2bp2, Casr, and Mrgprf, and the most strongly up-regulated mRNAs were Rtp3, Mef2b, Emr4, AABR07001734.1, Sell, Naa11, Vsx1, Tnfrsf13c, Acsm3, and RT1-Ba (Table 1). Major histocompatibility complex (MHC) class I (known as RT1-CE in the rat) protein expression affects neuroinflammatory processes and immune-mediated neurodegeneration both in Parkinson’s disease and in general [24, 25]. Serpin peptidase inhibitor plays an important role in synapse development and regulates synaptic plasticity [26, 27]. Serpind1 is very weakly expressed in brain tissues . MS4A family members (including Ms4a7) are likely involved in signal transduction in many different cell lineages . Igf2bp2, which regulates oxidative phosphorylation, binds several mRNAs that encode mitochondrial respiratory chain complex subunits and interacts with complex I (NADH:ubiquinone oxidoreductase) proteins . Igf2bp2 is strikingly enriched in developing axon tracts, including spinal commissural axons . MEF2 family proteins are involved in cell differentiation, proliferation, migration, and apoptosis . DNA methylation, which is an epigenetic mechanism commonly associated with gene silencing, regulates Naa11 . Mrgprf is classified as a MAS-related G-protein coupled receptor, and its function is currently unknown . The abnormal expression of these genes observed in this study might help illuminate new regulatory mechanisms in the aged SN.
Regulation of gene and protein expression by non-coding RNAs (ncRNAs) has become a popular research topic in recent years [35, 36]. NcRNAs, including lncRNAs and circRNAs, play an important role in brain biology and molecular pathologies associated with neurodegeneration [22, 37]. For example, lncRNAs and circRNAs affect cell proliferation, differentiation, quiescence, senescence, stress and immune response, and many other cellular functions associated with aging [12, 38–40]. LncRNAs can regulate gene expression in the following ways: inter-chromosomal interactions, formation of nuclear structures or R-loops, regulating post-transcriptional mRNA decay, regulating cellular localization of RNA- and DNA-binding proteins, acting as guides or decoys for transcription factors, providing scaffolds for chromatin modifying complexes, and acting as miRNA sponges . Large numbers of circRNAs that are biosynthesized from transcripts expressed in neuronal cells have been identified in neuronal tissues . CircRNAs can modulate miRNA function, regulate protein production and gene transcription, and act as miRNA sponges [22, 42]. The lncRNA HOTAIR affects differentiation of mesenchymal stem cells and is implicated in senescence-associated DNA-methylation . The muscleblind-like-3 (MBNL3) splicing factor promotes hepatocellular carcinoma (HCC) by increasing paxillin (PXN) expression through the alternative splicing of lncRNA-PXN-AS1. In HCC tissues, overexpression of MBNL3 promotes the inclusion of exon 4 of lncRNA-PXN-AS1, shifts the PXN-AS1-S isoform to PXN-AS1-L, and upregulates PXN mRNA and protein expression . The lncRNA uc.173 regulates growth of the intestinal mucosa and stimulates intestinal epithelial renewal by reducing levels of miRNA195 . Down-regulation of lncRNA MALAT1 attenuates neuronal cell death by regulating miR-30a expression, which in turn suppresses Beclin1-dependent autophagy, in cerebral ischemic stroke . The circRNA circAGFG1 exerts oncogenic effects by sponging miR-195-5p, which modulates cyclin E1 expression and promotes tumorigenesis and the development of triple-negative breast cancer . Here, we identified 163 DELs and 98 DECs in the SN of 24Mon rats. A single DEL might target multiple mRNAs, and multiple DELs can also target the same mRNA (Supplementary Table 10). In addition, while each DEC is derived from a single parental mRNA, a single mRNA might produce multiple DECs (Supplementary Table 15). GO and KEGG pathway analysis was performed to investigate the biological functions of the genes targeted by DELs and the parental genes of DECs. Our results indicate that those genes are involved in many processes and pathways. Thus, the complicated changes in regulatory mechanisms that occur in the SN during aging require further investigation.
We also used bioinformatics approaches to predict the potential regulatory roles of DEMs targeted by DELs and DEMs that were parents of DECs. Seven target relationships between four DELs and three DEMs were identified in the aged SN. In aged rats, the down-regulated lncRNA NONRATT010417.2 targeted the down-regulated mRNA Myh1, while the down-regulated lncRNA NONRATT015586.2 and the up-regulated lncRNAs NONRATT000490.2 and NONRATT007029.2 targeted the down-regulated mRNAs Casr and Mis18a. Immunoblotting and RT-qPCR revealed that expression of the Myh1, Casr, and Mis18a proteins and mRNAs was also significantly reduced in aged male rats. In addition, GO analysis revealed that the DELs targeting Myh1, Casr, and Mis18a are involved in behavior (GO:0007610), binding (GO:0005488), developmental processes (GO:0032502), metabolic processes (GO:0008152), signal transducer activity (GO:0004871), molecular transducer activity (GO:0060089), and other classifications. KEGG pathway analysis indicated that Myh1 is responsive to cellular processes, significantly enriched in the tight junction (rno04530), and related to cell polarity. Myosins constitute a superfamily of ubiquitous, multifunctional, and structurally related motor proteins that catalyze the hydrolysis of ATP to power directed movement on filamentous actin. The myosin superfamily plays fundamental roles in neuronal function, plasticity, morphogenesis, and survival . The non-muscle myosins are likely involved in neuronal motility and the extension and maintenance of axons and dendrites . Casr is a widely expressed G-protein coupled receptor. In neurons, Casr responds to changes in calcium levels that are involved in neuronal growth and migration, synaptic plasticity, and a variety of other processes . Down-regulation of Casr expression is associated with neuronal injury, calcium disturbances, increased production of reactive oxygen species, and decreased nitric oxide release . Mis18a, which is down-regulated in the amygdala during exposure to low-novelty environments, was enriched for GO terms relating to chromatin and nucleosome assembly and remodeling . Mis18a knockout causes severe chromosomal missegregation, inhibition of CENPA, and ultimately cell death . Reduced Mis18a expression is associated with problems in the formation of the nuclear spindle and chromosome segregation . Therefore, we hypothesized that low expression of Myh1, Casr, and Mis18a in the SN of aged male rats is likely regulated by lncRNA. However, the specific mechanism, such as the direct connection and regulation of the indicated lncRNAs with Myh1, Casr and Mis18a, remains unsolved in this study and should be investigated in vitro using cell culture model in future experiments.
Mitochondria play key roles in bioenergetics, generation of reactive oxygen species, apoptosis, signal transduction, calcium homeostasis, and cell survival and death [55, 56]. Mitochondrial dysfunction has emerged as a hallmark of the aging process and is linked to the development of numerous age-related pathologies, including neurodegenerative disorders [57–59]. Both nuclear genes and mitochondrial genes code for proteins that are involved in mitochondrial functions [60, 61]. Whether the expression of proteins related to mitochondrial function is also regulated by ncRNAs remains largely unknown. In our study, lncRNA NONRATT015586.2, which was down-regulated, and lncRNAs NONRATT000490.2 and NONRATT007029.2, which were up-regulated, all targeted the ATP synthase mitochondrial F1 complex assembly factor 1 (Atpaf1), mitochondrial ribosomal protein 63 (Mrp63), and mitochondrial ribosomal protein L51 (Mrpl51) mRNAs; there were no age-dependent differences in the expression of these mRNAs (Supplementary Table 10). Unfortunately, we were not able to examine expression of the Atpaf1, Mrp63, and Mrpl51 proteins. Studies have shown that abnormal expression of Atpaf1 and Mrpl51 is closely related to mitochondrial dysfunction [62, 63]. It is possible that differential expression of these mRNAs was not observed in this study due to the small number of samples examined (n=3) or because the identified lncRNAs affected translation rather than transcription.
In conclusion, we characterized mRNAs, lncRNAs, and circRNAs that were differentially expressed in the SN of aged and young rats. We found that the indicated lncRNAs down-regulated the expression of Myh1, Casr, and Mis18a in the aged rat SN. These findings help to elucidate transcriptional underpinnings of aging in the SN and provide a foundation for future studies on associated molecular mechanisms.
Materials and Methods
Animals housing and ethics
A total of six male Wistar rats were assigned to either the 6Mon (n = 3) or the 24Mon (n = 3) group. All rats were supplied by the Experimental Animal Center of Hebei Medical University and were housed in an air-conditioned room (22 ± 2°C) on a 12-h light–dark cycle (lights on 06:00 h). Food and water were available ad libitum. All experimental procedures were conducted in accordance with the rules in the “Guidelines for the Care and Use of Mammals in Neuroscience and Behavioral Research” and were approved by the Committee of Ethics in Animal Experiments at Hebei Medical University.
All rats in each group were sacrificed by decapitation. The brains were removed quickly and the blocks containing the SN in the ventral midbrain (between 2.96 and 3.70 mm, Paxinos and Watson, 1998) was dissected on an ice-cold plate using a scalpel for ophthalmic surgery and a stereomicroscope. Tissue blocks containing the SN were processed for high-throughput RNA sequencing and western blot analysis.
RNA extraction and purification
Total RNA was extracted from the SN using the RNeasy Mini Kit (Cat 74106, Qiagen) according the manufacturer’s instructions; RNA integrity number (RIN, from 0 to 10) was determined using an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, US). Higher RIN scores indicate greater RNA integrity; only qualified total RNA samples with RIN ≥ 7.0 were further purified using the RNAClean XP Kit (Cat A63987, Beckman Coulter, Inc., Kraemer Boulevard, Brea, CA, USA) and the RNase-Free DNase Set (Cat 79254, QIAGEN, GmBH, Germany).
Sequencing library construction and sequencing
Purified SN RNA samples were used for library preparation. The TruSeq RNA Sample Preparation Kit (RS-122-2001, Illumina, San Francisco, CA, USA) was then used to synthesize paired-end libraries according to the manufacturer’s guidelines. Library construction, Illumina sequencing, sequence analysis of high-quality reads, and bioinformatic data analysis were performed by the Shanghai Biotechnology Corporation.
RNA sequence mapping and bioinformatic analysis of differentially expressed mRNAs, lncRNAs, and circRNAs
Before read mapping, clean reads were isolated from the raw reads using Seqtk (https://github.com/lh3/seqtk). Hisat2 (version 2.0.4) was used for alignment of clean reads and genome mapping . Relative expression levels of all matched unigenes were normalized by transforming the clean reads into fragment kilobase of exon per million fragments mapped (FPKM) for mRNAs and lncRNAs and spliced reads per billion mapping (SRPBM) for circRNAs. RNAs were considered significantly differentially expressed under the following criteria: (1) fold-change (FC) value > 2 (|log2FC|>1), and (2) adjusted false discovery rate (FDR), which was a q-value < 0.05 for mRNAs and lncRNAs and a p-value < 0.05 for circRNAs.
GO and KEGG analysis
To better understand the biological functions and potential regulatory roles of dysregulated mRNAs, lncRNAs, and circRNAs, we performed GO (http://geneontology.org) and KEGG (https://www. kegg.jp) pathway and enrichment analysis. GO analysis involved classifying associations between DEMs, mRNAs targeted by DELs, and parental mRNAs of DECs and the following three categories: biological process (BP), cellular component (CC), and molecular function (MF). GO terms were widely used as functional enrichment analysis for many genes. KEGG pathway analysis included classifications and pathway enrichment analysis. A q-value < 0.05 denoted significant enrichments in the GO and KEGG pathways.
Western blot analysis
The tissue blocks used for detection of Myh1, Casr, and Mis18a proteins were homogenized in radioimmunoprecipitation assay (RIPA) buffer containing 1% Triton X-100, 0.1% SDS, 0.5% sodium deoxycholate, and protease inhibitors (100 μg/mL phenylmethanesulfonyl fluoride, 30 μg/mL aprotinin, and 1 mM sodium orthovanadate), and then sonicated four times for 10 s each. After centrifugation at 12,000×g for 20 min at 4°C, the supernatant was collected and centrifuged again as described above. The resulting final supernatant was stored at -80°C until use. Samples from the SN (50 μg) were diluted with 2× sample buffer (50 mM Tris, pH 6.8, 2% SDS, 10% glycerol, 0.1% bromophenol blue, and 5% β-mercaptoethanol), heated for 5 min at 95°C before sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) on a 10% gel, and subsequently transferred to a polyvinylidene fluoride (PVDF) membrane (Millipore). The membrane was incubated for 2 h with 5% nonfat dry milk in Tris-buffered saline (TBS) containing 0.05% Tween 20 (TBST) (20 mM Tris-Cl, 137 mM NaCl, 0.1% Tween 20, pH 7.6) at room temperature. The membrane was rinsed in TBST three times and then incubated overnight with rabbit anti-Myh1 polyclonal antibody (Sigma, SAB2104768, 1:100), rabbit anti-Casr polyclonal antibody (Sigma, SAB4503369, 1:500), or rabbit anti-Mis18a polyclonal antibody (Biorbyt, orb325235, 1:200) at 4°C. After three washes, the membrane was incubated for 2 h in horseradish peroxidase-conjugated anti-rabbit IgG (1:6000). Immunoreactive bands were visualized using enhanced chemiluminescence (ECL, Amersham Biosciences) and film exposure (Kodak XBT-1). Following stripping, each PVDF membrane was subsequently immunoblotted with mouse anti-β-actin monoclonal antibody (Santa Cruz Biotechnology, Santa Cruz, CA, USA; sc-47778). The labeling densities for Myh1, Casr, and Mis18a were compared with those of β-actin, which served as the endogenous control. The specific bands were recorded by a scanner and analyzed by transmittance densitometry using Gel-Pro Analyzer Analysis software (Media Cybernetics).
Real-time quantitative polymerase chain reaction (RT-PCR)
Total RNA (2 μg) from the SN was subjected to reverse transcription using random primer to obtain the first-strand cDNA template. Real-time fluorescence quantitative PCR was performed with 0.8-μl cDNA (diluted 1:10), 2-μl specific primers, and 2X GoTaq® Green Master Mix (Promega, USA) with a final volume of 20 μl. PCR was performed as follows: an initial cycle at 95 °C for 10min, followed by 40 cycles at 95 °C for 15 s, 58 °C for 20 s, and 72 °C for 27 s. The products were then analyzed by melting curve to confirm the specificity of amplification. Expression of Myh1, Casr and Mis18a genes was analyzed using glyceraldehyde-3-phosphate dehydrogenase (GAPDH) as an internal control. The sets of primers were as follows: Myh1 (5′-CCTGGATGA TCTACACCTACTC-3′ and 5′-GTCAGAGATAGAGA AGATGTGGG-3′), Casr (5′-TGACGAGCCTCAGAAG AATGC-3′ and 5′-CCTCCACCACTAATGACGAAGC-3), Mis18a (5′-TGGAACCCTGCTCTACCTCT-3′ and 5′-CCAATGGCGGCTCGAATCTT-3′), GAPDH (5′-TG AACGGGAAGCTCACTG-3′ and 5′-GCTTCACCAC CTTCTTGATG-3′).
Data are reported as means ± SD. Independent-samples t-tests were used for statistical analysis. A value of p < 0.05 was considered significant.
Conflicts of Interest
The authors have no conflicts of interest to declare.
This project was financially supported by the Natural Science Foundation of China (No. 81200252, 81871119), the Natural Science Foundation of Hebei Province of China (No. C2017206072), the Growth Fund of the Chunyu Plan of Hebei Medical University (CYCZ2019004), and the University Student’s Innovation Project of Hebei Medical University (USIP2018199).
- 1. Błaszczyk JW. Nigrostriatal interaction in the aging brain: new therapeutic target for Parkinson’s disease. Acta Neurobiol Exp (Wars). 2017; 77:106–112. https://doi.org/10.21307/ane-2017-041 [PubMed]
- 2. Beitz JM. Parkinson’s disease: a review. Front Biosci (Schol Ed). 2014; 6:65–74. https://doi.org/10.2741/s415 [PubMed]
- 3. Paladino S, Conte A, Caggiano R, Pierantoni GM, Faraonio R. Nrf2 Pathway in Age-Related Neurological Disorders: Insights into MicroRNAs. Cell Physiol Biochem. 2018; 47:1951–76. https://doi.org/10.1159/000491465 [PubMed]
- 4. Freundlieb N, François C, Tandé D, Oertel WH, Hirsch EC, Höglinger GU. Dopaminergic substantia nigra neurons project topographically organized to the subventricular zone and stimulate precursor cell proliferation in aged primates. J Neurosci. 2006; 26:2321–5. https://doi.org/10.1523/jneurosci.4859-05.2006 [PubMed]
- 5. Costa KM. The effects of aging on substantia nigra dopamine neurons. J Neurosci. 2014; 34:15133–4. https://doi.org/10.1523/jneurosci.3739-14.2014 [PubMed]
- 6. Zhang GL, Wang W, Kang YX, Xue Y, Yang H, Zhou CM, Shi GM. Chronic testosterone propionate supplement could activated the Nrf2-ARE pathway in the brain and ameliorated the behaviors of aged rats. Behav Brain Res. 2013; 252:388–95. https://doi.org/10.1016/j.bbr.2013.05.063 [PubMed]
- 7. Zhang G, Li S, Kang Y, Che J, Cui R, Song S, Cui H, Shi G. Enhancement of dopaminergic activity and region-specific activation of Nrf2-ARE pathway by intranasal supplements of testosterone propionate in aged male rats. Horm Behav. 2016; 80:103–116. https://doi.org/10.1016/j.yhbeh.2016.02.001 [PubMed]
- 8. Rekaik H, Blaudin de Thé FX, Fuchs J, Massiani-Beaudoin O, Prochiantz A, Joshi RL. Engrailed homeoprotein protects mesencephalic dopaminergic neurons from oxidative stress. Cell Rep. 2015; 13:242–250. https://doi.org/10.1016/j.celrep.2015.08.076 [PubMed]
- 9. Pal R, Tiwari PC, Nath R, Pant KK. Role of neuroinflammation and latent transcription factors in pathogenesis of Parkinson’s disease. Neurol Res. 2016; 38:1111–1122. https://doi.org/10.1080/01616412.2016.1249997 [PubMed]
- 10. Parkinson GM, Dayas CV, Smith DW. Increased mitochondrial DNA deletions in substantia nigra dopamine neurons of the aged rat. Curr Aging Sci. 2014; 7:155–60. https://doi.org/10.2174/1874609808666150122150850 [PubMed]
- 11. Guennewig B, Cooper AA. The central role of noncoding RNA in the brain. Int Rev Neurobiol. 2014; 116:153–94. https://doi.org/10.1016/b978-0-12-801105-8.00007-2 [PubMed]
- 12. Arendt T, Ueberham U, Janitz M. Non-coding transcriptome in brain aging. Aging (Albany NY). 2017; 9:1943–44. https://doi.org/10.18632/aging.101290 [PubMed]
- 13. Kour S, Rath PC. Long noncoding RNAs in aging and age-related diseases. Ageing Res Rev. 2016; 26:1–21. https://doi.org/10.1016/j.arr.2015.12.001 [PubMed]
- 14. Grammatikakis I, Panda AC, Abdelmohsen K, Gorospe M. Long noncoding RNAs (lncRNAs) and the molecular hallmarks of aging. Aging (Albany NY). 2014; 6:992–1009. https://doi.org/10.18632/aging.100710 [PubMed]
- 15. Kim C, Kang D, Lee EK, Lee JS. Long Noncoding RNAs and RNA-Binding Proteins in Oxidative Stress, Cellular Senescence, and Age-Related Diseases. Oxid Med Cell Longev. 2017; 2017:2062384. https://doi.org/10.1155/2017/2062384 [PubMed]
- 16. Kumar L, Shamsuzzama, Haque R, Baghel T, Nazir A. Circular RNAs: the Emerging Class of Non-coding RNAs and Their Potential Role in Human Neurodegenerative Diseases. Mol Neurobiol. 2017; 54:7224–7234. https://doi.org/10.1007/s12035-016-0213-8 [PubMed]
- 17. Haque S, Harries LW. Circular RNAs (circRNAs) in health and disease. Genes (Basel). 2017; 8:1513. https://doi.org/10.3390/genes8120353 [PubMed]
- 18. Pan L, Ma J, Pan F, Zhao D, Gao J. Long non-coding RNA expression profiling in aging rats with erectile dysfunction. Cell Physiol Biochem. 2015; 37:1513–26. https://doi.org/10.1159/000438519 [PubMed]
- 19. Shen S, Jiang H, Bei Y, Xiao J, Li X. Long non-coding RNAs in cardiac remodeling. Cell Physiol Biochem. 2017; 41:1830–37. https://doi.org/10.1161/circresaha.113.300226 [PubMed]
- 20. Marchese FP, Raimondi I, Huarte M. The multidimensional mechanisms of long noncoding RNA function. Genome Biol. 2017; 18:206. https://doi.org/10.1186/s13059-017-1348-2 [PubMed]
- 21. Qu S, Zhong Y, Shang R, Zhang X, Song W, Kjems J, Li H. The emerging landscape of circular RNA in life processes. RNA Biol. 2017; 14:992–9. https://doi.org/10.1080/15476286.2016.1220473 [PubMed]
- 22. Floris G, Zhang L, Follesa P, Sun T. Regulatory role of circular RNAs and neurological disorders. Mol Neurobiol. 2017; 54:5156–65. https://doi.org/10.1007/s12035-016-0055-4 [PubMed]
- 23. Li X, Yang L, Chen LL. The Biogenesis, functions, and challenges of circular RNAs. Mol Cell. 2018; 71:428–42. https://doi.org/10.1016/j.molcel.2018.06.034 [PubMed]
- 24. Cebrián C, Zucca FA, Mauri P, Steinbeck JA, Studer L, Scherzer CR, Kanter E, Budhu S, Mandelbaum J, Vonsattel JP, Zecca L, Loike JD, Sulzer D. MHC-I expression renders catecholaminergic neurons susceptible to T-cell-mediated degeneration. Nat Commun. 2014; 5:3633. https://doi.org/10.1038/ncomms4633 [PubMed]
- 25. Cebrián C, Loike JD, Sulzer D. Neuronal MHC-I expression and its implications in synaptic function, axonal regeneration and Parkinson’s and other brain diseases. Front Neuroanat. 2014; 8:114. https://doi.org/10.3389/fnana.2014.00114 [PubMed]
- 26. Teesalu T, Kulla A, Simisker A, Sirén V, Lawrence DA, Asser T, Vaheri A. Tissue plasminogen activator and neuroserpin are widely expressed in the human central nervous system. Thromb Haemost. 2004; 92:358–68. https://doi.org/10.1160/th02-12-0310 [PubMed]
- 27. Kinghorn KJ, Crowther DC, Sharp LK, Nerelius C, Davis RL, Chang HT, Green C, Gubb DC, Johansson J, Lomas DA. Neuroserpin binds Abeta and is a neuroprotective component of amyloid plaques in Alzheimer disease. J Biol Chem. 2006; 281:29268–77. https://doi.org/10.1074/jbc.m600690200 [PubMed]
- 28. Woodward KJ, Stampalia J, Vanyai H, Rijhumal H, Potts K, Taylor F, Peverall J, Grumball T, Sivamoorthy S, Alinejad-Rokny H, Wray J, Whitehouse A, Nagarajan L, Scurlock J, Afchani S, Edwards M, Murch A, Beilby J, Baynam G, Kiraly-Borri C, McKenzie F, Heng JIT. Atypical nested 22q11.2 duplications between LCR22B and LCR22D are associated with neurodevelopmental phenotypes including autism spectrum disorder with incomplete penetrance. Mol Genet Genomic Med. 2019; 7:e00507. https://doi.org/10.1002/mgg3.507 [PubMed]
- 29. Liang Y, Buckley TR, Tu L, Langdon SD, Tedder TF. Structural organization of the human MS4A gene cluster on Chromosome 11q12. Immunogenetics. 2001; 53:357–68. https://doi.org/10.1007/s002510100339 [PubMed]
- 30. Janiszewska M, Suvà ML, Riggi N, Houtkooper RH, Auwerx J, Clément-Schatlo V, Radovanovic I, Rheinbay E, Provero P, Stamenkovic I. Imp2 controls oxidative phosphorylation and is crucial for preserving glioblastoma cancer stem cells. Genes Dev. 2012; 26:1926–44. https://doi.org/10.1101/gad.188292.112 [PubMed]
- 31. Preitner N, Quan J, Li X, Nielsen FC, Flanagan JG. IMP2 axonal localization, RNA interactome, and function in the development of axon trajectories. Development. 2016; 143:2753–9. https://doi.org/10.1242/dev.128348 [PubMed]
- 32. Pon JR, Marra MA. MEF2 transcription factors: developmental regulators and emerging cancer genes. Oncotarget. 2016; 7:2297–312. https://doi.org/10.18632/oncotarget.6223 [PubMed]
- 33. Pang AL, Clark J, Chan WY, Rennert OM. Expression of human NAA11 (ARD1B) gene is tissue-specific and is regulated by DNA methylation. Epigenetics. 2011; 6:1391–9. https://doi.org/10.4161/epi.6.11.18125 [PubMed]
- 34. Bader M, Alenina N, Andrade-Navarro MA, Santos RA. MAS and its related G protein-coupled receptors, Mrgprs. Pharmacol Rev. 2014; 66:1080–105. https://doi.org/10.1124/pr.113.008136 [PubMed]
- 35. Subhramanyam CS, Hu Q. Non-Coding RNA in brain development and disorder. Curr Med Chem. 2017; 24:1983–97. https://doi.org/10.2174/0929867324666170124151436 [PubMed]
- 36. Zhou S, Ding F, Gu X. Non-coding RNAs as Emerging regulators of neural injury responses and regeneration. Neurosci Bull. 2016; 32:253–64. https://doi.org/10.1007/s12264-016-0028-7 [PubMed]
- 37. Idda ML, Munk R, Abdelmohsen K, Gorospe M. Noncoding RNAs in Alzheimer’s disease. Wiley Interdiscip Rev RNA. 2018; 9:10.1002/wrna.1463. https://doi.org/10.1002/wrna.1463 [PubMed]
- 38. Montes M, Lund AH. Emerging roles of lncRNAs in senescence. FEBS J. 2016; 283:2414–26. https://doi.org/10.1111/febs.13679 [PubMed]
- 39. Abdelmohsen K, Panda AC, De S, Grammatikakis I, Kim J, Ding J, Noh JH, Kim KM, Mattison JA, de Cabo R, Gorospe M. Circular RNAs in monkey muscle: age-dependent changes. Aging (Albany NY). 2015; 7:903–10. https://doi.org/10.18632/aging.100834 [PubMed]
- 40. Panda AC, Grammatikakis I, Munk R, Gorospe M, Abdelmohsen K. Emerging roles and context of circular RNAs. Wiley Interdiscip Rev RNA. 2017; 8:10.1002/wrna.1386. https://doi.org/10.1002/wrna.1386 [PubMed]
- 41. Ashwal-Fluss R, Meyer M, Pamudurti NR, Ivanov A, Bartok O, Hanan M, Evantal N, Memczak S, Rajewsky N, Kadener S. circRNA biogenesis competes with pre-mRNA splicing. Mol Cell. 2014; 56:55–66. https://doi.org/10.1016/j.molcel.2014.08.019 [PubMed]
- 42. Qu S, Yang X, Li X, Wang J, Gao Y, Shang R, Sun W, Dou K, Li H. Circular RNA: A new star of noncoding RNAs. Cancer Lett. 2015; 365:141–8. https://doi.org/10.1016/j.canlet.2015.06.003 [PubMed]
- 43. Kalwa M, Hänzelmann S, Otto S, Kuo CC, Franzen J, Joussen S, Fernandez-Rebollo E, Rath B, Koch C, Hofmann A, Lee SH, Teschendorff AE, Denecke B, Lin Q, Widschwendter M, Weinhold E, Costa IG, Wagner W. The lncRNA HOTAIR impacts on mesenchymal stem cells via triple helix formation. Nucleic Acids Res. 2016; 44:10631–43. https://doi.org/10.1093/nar/gkw802 [PubMed]
- 44. Yuan JH, Liu XN, Wang TT, Pan W, Tao QF, Zhou WP, Wang F, Sun SH. The MBNL3 splicing factor promotes hepatocellular carcinoma by increasing PXN expression through the alternative splicing of lncRNA-PXN-AS1. Nat Cell Biol. 2017; 19:820–32. https://doi.org/10.1038/ncb3538 [PubMed]
- 45. Xiao L, Wu J, Wang JY, Chung HK, Kalakonda S, Rao JN, Gorospe M, Wang JY. Long Noncoding RNA uc.173 Promotes Renewal of the Intestinal Mucosa by Inducing Degradation of MicroRNA 195. Gastroenterology. 2018; 154:599–611. https://doi.org/10.1053/j.gastro.2017.10.009 [PubMed]
- 46. Guo D, Ma J, Yan L, Li T, Li Z, Han X, Shui S. Down-Regulation of Lncrna MALAT1 Attenuates Neuronal Cell Death Through Suppressing Beclin1-Dependent Autophagy by Regulating Mir-30a in Cerebral Ischemic Stroke. Cell Physiol Biochem. 2017; 43:182–194. https://doi.org/10.1159/000480337 [PubMed]
- 47. Yang R, Xing L, Zheng X, Sun Y, Wang X, Chen J. The circRNA circAGFG1 acts as a sponge of miR-195-5p to promote triple-negative breast cancer progression through regulating CCNE1 expression. Mol Cancer. 2019; 18:4. https://doi.org/10.1186/s12943-018-0933-7 [PubMed]
- 48. Hirokawa N, Niwa S, Tanaka Y. Molecular motors in neurons: transport mechanisms and roles in brain function, development, and disease. Neuron. 2010; 68:610–38. https://doi.org/10.1016/j.neuron.2010.09.039 [PubMed]
- 49. Miller M, Bower E, Levitt P, Li D, Chantler PD. Myosin II distribution in neurons is consistent with a role in growth cone motility but not synaptic vesicle mobilization. Neuron. 1992; 8:25–44. https://doi.org/10.1016/0896-6273(92)90106-n [PubMed]
- 50. Ruat M, Traiffort E. Roles of the calcium sensing receptor in the central nervous system. Best Pract Res Clin Endocrinol Metab. 2013; 27:429–42. https://doi.org/10.1016/j.beem.2013.03.001 [PubMed]
- 51. Dong S, Li G, Zheng D, Wu J, Sun D, Yang F, Yu X, Li T, Sun A, Liu J, Zhong X, Xu C, Lu F, Zhang W. A novel role for the calcium sensing receptor in rat diabetic encephalopathy. Cell Physiol Biochem. 2015; 35:38–50. https://doi.org/10.1159/000369673 [PubMed]
- 52. Cohen JL, Ata AE, Jackson NL, Rahn EJ, Ramaker RC, Cooper S, Kerman IA, Clinton SM. Differential stress induced c-Fos expression and identification of region-specific miRNA-mRNA networks in the dorsal raphe and amygdala of high-responder/low-responder rats. Behav Brain Res. 2017; 319:110–23. https://doi.org/10.1016/j.bbr.2016.11.015 [PubMed]
- 53. Baumann K. Keeping centromeric identity. Nat Rev Mol Cell Biol. 2012; 13:340. https://doi.org/10.1038/nrm3356 [PubMed]
- 54. Pussila M, Törönen P, Einarsdottir E, Katayama S, Krjutškov K, Holm L, Kere J, Peltomäki P, Mäkinen MJ, Linden J, Nyström M. Mlh1 deficiency in normal mouse colon mucosa associates with chromosomally unstable colon cancer. Carcinogenesis. 2018; 39:788–97. https://doi.org/10.1093/carcin/bgy056 [PubMed]
- 55. Li PA, Hou X, Hao S. Mitochondrial biogenesis in neurodegeneration. J Neurosci Res. 2017; 95:2025–9. https://doi.org/10.1002/jnr.24042 [PubMed]
- 56. Aufschnaiter A, Kohler V, Diessl J, Peselj C, Carmona-Gutierrez D, Keller W, Büttner S. Mitochondrial lipids in neurodegeneration. Cell Tissue Res. 2017; 367:125–40. https://doi.org/10.1007/s00441-016-2463-1 [PubMed]
- 57. Kong Y, Trabucco SE, Zhang H. Oxidative stress, mitochondrial dysfunction and the mitochondria theory of aging. Interdiscip Top Gerontol. 2014; 39:86–107. https://doi.org/10.1159/000358901 [PubMed]
- 58. Srivastava S. The mitochondrial basis of aging and age-related disorders. Genes (Basel). 2017; 8:398. https://doi.org/10.3390/genes8120398 [PubMed]
- 59. Jha SK, Jha NK, Kumar D, Ambasta RK, Kumar P. Linking mitochondrial dysfunction, metabolic syndrome and stress signaling in Neurodegeneration. Biochim Biophys Acta Mol Basis Dis. 2017; 1863:1132–46. https://doi.org/10.1016/j.bbadis.2016.06.015 [PubMed]
- 60. Koopman WJ, Distelmaier F, Smeitink JA, Willems PH. OXPHOS mutations and neurodegeneration. EMBO J. 2013; 32:9–29. https://doi.org/10.1038/emboj.2012.300 [PubMed]
- 61. Ladd AC, Brohawn DG, Thomas RR, Keeney PM, Berr SS, Khan SM, Portell FR, Shakenov MZ, Antkowiak PF, Kundu B, Tustison N, Bennett JP. RNA-seq analyses reveal that cervical spinal cords and anterior motor neurons from amyotrophic lateral sclerosis subjects show reduced expression of mitochondrial DNA-encoded respiratory genes, and rhTFAM may correct this respiratory deficiency. Brain Res. 2017; 1667:74–83. https://doi.org/10.1016/j.brainres.2017.05.010 [PubMed]
- 62. Smiles WJ, Parr EB, Coffey VG, Lacham-Kaplan O, Hawley JA, Camera DM. Protein coingestion with alcohol following strenuous exercise attenuates alcohol-induced intramyocellular apoptosis and inhibition of autophagy. Am J Physiol Endocrinol Metab. 2016; 311:E836–E849. https://doi.org/10.1152/ajpendo.00303.2016 [PubMed]
- 63. Haque ME, Elmore KB, Tripathy A, Koc H, Koc EC, Spremulli LL. Properties of the C-terminal tail of human mitochondrial inner membrane protein Oxa1L and its interactions with mammalian mitochondrial ribosomes. J Biol Chem. 2010; 285:28353–62. https://doi.org/10.1074/jbc.m110.148262 [PubMed]
- 64. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015; 12:357–60. https://doi.org/10.1038/nmeth.3317 [PubMed]