Original Research ARTICLE
Long non-coding RNA expression during aging in the human subependymal zone
- 1Garvan Institute of Medical Research, Sydney, NSW, Australia
- 2St Vincent’s Clinical School and School of Biotechnology and Biomolecular Sciences, University of New South Wales, Sydney, NSW, Australia
- 3Schizophrenia Research Institute, Sydney, NSW, Australia
- 4Schizophrenia Research Laboratory, Neuroscience Research Australia, Sydney, NSW, Australia
- 5School of Psychiatry, University of New South Wales, Sydney, NSW, Australia
The human subependymal zone (SEZ) is debatably a source of newly born neurons throughout life and neurogenesis is a multi-step process requiring distinct transcripts during cell proliferation and early neuronal maturation, along with orchestrated changes in gene expression during cell state/fate transitions. Furthermore, it is becoming increasingly clear that the majority of our genome that results in production of non-protein-coding RNAs plays vital roles in the evolution, development, adaptation, and region-specific function of the human brain. We predicted that some transcripts expressed in the SEZ may be unique to this specialized brain region, and that a comprehensive transcriptomic analysis of this region would aid in defining expression changes during neuronal birth and growth in adult humans. Here, we used deep RNA sequencing of human SEZ tissue during adulthood and aging to characterize the transcriptional landscape with a particular emphasis on long non-coding RNAs (lncRNAs). The data show predicted age-related changes in mRNAs encoding proliferation, progenitor, and inflammatory proteins as well as a unique subset of lncRNAs that are highly expressed in the human SEZ, many of which have unknown functions. Our results suggest the existence of robust proliferative and neuronal differentiation potential in the adult human SEZ and lay the foundation for understanding the involvement of lncRNAs in postnatal neurogenesis and potentially associated neurodevelopmental diseases that emerge after birth.
The human subependymal zone (SEZ; also known as the subventricular zone; SVZ), found adjacent to the lateral ventricle, is the largest reservoir of newly born neuronal and glial cells in the adult brain and, in primates, is estimated to produce ~10 times the number of new neurons compared to hippocampal neurogenesis (1). In late embryonic development, the area lying juxtaposed to the lateral wall of the lateral ventricle and just medial to the ganglionic eminence is responsible for the generation of inhibitory interneurons destined for the cortex and caudate. In adult rodents, proliferative capacity is retained in the SEZ and is critically important to the genesis and replacement of olfactory bulb interneurons that travel from the SEZ via the rostral migratory stream to the olfactory bulb throughout life (2). The extent of SEZ proliferation, neurogenesis, and the existence of a rostral migratory stream is debated in adult humans with some evidence supporting active neurogenesis and migration postnatally, especially in the first few years of human life (3–7). There are also examples of reactive neurogenesis in the human SEZ in response to injury or neural degeneration indicating that this region harbors precursor cells throughout life and implicating activation of this region following brain damage (1). This demonstrates that the adult brain may not only be capable of generating new cells from the SEZ but also that these cells could travel throughout the brain if the necessary instructive cues were present. If the appropriate temporal and spatial requirements of developing cells within the SEZ are not met, altered production of interneuron precursors in the SEZ may also contribute to interneuron deficiencies and the development of brain disorders, such as a host of psychiatric disorders including schizophrenia and autism spectrum disorder (8). There is also evidence that the neurogenic potential of the mammalian brain is attenuated with the normal aging process (9) suggesting that transcriptional events involved in SEZ neurogenesis may also change as humans age.
To more fully understand functional genomic expression in specialized regions of the human brain, such as the SEZ, requires the exploration of the entire transcriptome as the human genome contains only ~1.5% protein-coding sequences (10) but at least 80% of the genome is dynamically transcribed resulting in a transcriptional landscape dominated by long and short non-protein-coding RNAs (11, 12). Long non-coding RNAs (lncRNAs) are expressed from conventional promoters and share many features with protein-coding transcripts, including intron–exon boundaries and alternative splicing; they also commonly possess some distinct features in that many are nuclear-localized and are far more cell- and tissue-specific than protein-coding transcripts (13). There are over 56,000 lncRNAs currently annotated in the human genome (14), many of which have appeared relatively recently in evolution, with about one-third being primate-specific (15).
Long non-coding RNAs display highly specific spatial and temporal expression patterns (13, 16–18), implying that their functions lie in the refinement of regulatory circuits specific to particular cell types and activities, especially the epigenetic control of differentiation and development (19, 20). The majority of lncRNAs are expressed in the brain, many exclusively (13, 15). It is, therefore, unsurprising that they have critical roles in brain development (21) and are emerging as significant contributors to brain disease, such as schizophrenia (22), when dysregulated.
Long non-coding RNAs have been implicated in neurogenesis (23) and functional studies are accumulating demonstrating the mechanisms underpinning their involvement. For example, knockdown studies have identified lncRNAs that are specifically required for neurogenesis, with three lncRNAs (designated N1–N3) shown to be critical for proper neuronal differentiation (24). These lncRNAs also associate with chromatin remodeling complexes, a common theme among lncRNAs such as REST and SUZ12, to control cell fate. Moreover, Ng et al. have demonstrated that RMST, an lncRNA regulated by REST and induced during neurogenesis, interacts with the transcription factor SOX2 and is required for its recruitment to key neurogenesis promoting genes (25). Furthermore, distal-less homeobox 1 (DLX1) while a protein-coding gene, plays an important role in the specification of interneuron subtypes in the SEZ (26, 27), the activity of the sense strand mRNA may be modulated by an antisense lncRNA in mice, termed Dlx1as (28). Although, Dlx1as is thought not to exist in humans (28).
A deeper transcriptomic analysis of human neurogenic regions may be required to determine if there are known and unique lncRNAs expressed in the human SEZ that may be involved in interneuron genesis. Although the SEZ is strongly linked to neurodevelopmental programs and may be a source of adult neurogenesis, there has been no study to date investigating the non-coding RNA molecular landscape of this region and transcriptional changes during aging. Therefore, we sought to characterize age-dependent expression of protein-coding RNA and lncRNA transcripts to gain insight into SEZ neurogenesis.
Materials and Methods
Human Subependymal Zone Tissue
Tissue from the anterior caudate of 21 normal individuals between the ages of 21 and 81 years was obtained from the New South Wales Tissue Resource Centre (Sydney, NSW, Australia; HREC 07261). Cases were screened for neuropathology and toxicology and were free of psychiatric or neurological disease. The cause of death for most cases was cardiac (cardiac n = 17, cardiac/respiratory n = 1, unknown n = 2, cancer n = 1) and the cohort consisted of 3 females and 18 males, with an average age of 53.95 years, average PMI of 30.1, average pH of 6.68, and average RIN of 7.3 (Table S1 in Supplementary Material). Fresh-frozen caudate tissue was sectioned on a Leica CM3050 S cryostat, taking sets of 20μm × 60 μm sections interspersed with 10μm × 14 μm sections. Tissue for RNA extraction was dissected from 60 μm thick sections while frozen. Cuts were made ~2 mm deep to the lateral surface of the lateral ventricle to include the SEZ using Wescott spring scissors (T106, ProSci Tech). For each case, tissue was dissected from 3 sets of 3–4 adjacent 60 μm sections spaced ~1340 μm to give 10 sections/case (~40 mg tissue total/case).
Total RNA was extracted from SEZ-containing tissue for each case using Trizol (Invitrogen), and RNA quality was assessed on an Agilent Technologies 2100 Bioanalyzer with an RNA 6000 Nano kit (Agilent Technologies, USA) according to the manufacturer’s instructions. cDNA was synthesized from 3 μg total RNA per case using SuperScript® III First-Strand Synthesis kit and random hexamers (Invitrogen). For LINC00657, NEAT1, Dlx1, and Dlx1 enhancer, 2 μg of total RNA was DNase treated using TURBO DNA-free kit (Ambion) prior to cDNA synthesis to minimize genomic DNA amplification.
Quantitative Real-Time PCR
TaqMan gene expression assays
Transcript levels were measured by quantitative real-time PCR (qPCR) on an ABI Prism 7900HT Fast Real-time PCR system with 384-well format and TaqMan Gene Expression Assay [Applied Biosystems; Ki67, Hs010324433; see Ref. (29)]. Changes in expression of mRNAs for one coding (DCX) and two expressed lncRNAs, one novel (LINC00657) and one known (NEAT1), with age were validated with qPCR in the entire cohort with a total of 21 cases (Table S1 in Supplementary Material). TaqMan gene expression assays Hs01035496_m1 and Hs01008264_s1 were used to detect DCX and NEAT1 mRNAs, respectively, and custom probe AI6RO2N was used to detect LINC00657. All measurements from each subject were performed in duplicate and relative quantities determined from a seven-point standard curve of pooled cDNA. Transcript quantities for Ki67 were normalized by the geometric mean of four housekeeping genes: ubiquitin C (Hs00824723_m1), actin β (Hs99999903_m1), glyceraldehyde-3-phosphate dehydrogenase (Hs99999905_m1), TATA box binding protein (Hs00427620_m1) that did not correlate with age (r = −0.02, p > 0.05). Quantity means for DCX, LINC00657, and NEAT1 were normalized to glyceraldehyde-3-phosphate dehydrogenase (Hs99999905_m1) transcript expression.
Quantitative real-time PCR was also conducted using KAPA SYBR Fast qPCR universal kit (KAPA Biosystems, USA) according to manufacturer’s instructions in 384-well plates. Primers used were DLX1 (forward: CTCAGGTCAAGATCTGGTTC; reverse: GGATGAAGAGTTAGGGTTCC), DLX1 enhancer (forward: CGAGGATTAACACTTCCTGAA; reverse: GGGAGTGATTATGTATGCACC), and GAPDH (forward: CAGCCTCAAGATCATCAGCA; reverse: ATGGACTGTGGTCATGAGTC). Ten microliter qPCR reactions were performed in triplicate using 2 μL of diluted cDNA per reaction. Quantitative PCR reagent master mixes included 2× KAPA SYBR FAST qPCR Master Mix (2×) Universal, ROX-high reference dye and 200 nM final concentration of each forward and reverse primer (Integrated DNA Technologies, Inc., USA). An Applied Biosystems 7900HT Fast Real-Time PCR machine fitted with a 384-well thermal block was used for qPCR (Life Technologies, USA). Thermal cycling conditions consisted of initial 3 min enzyme activation at 95°C followed by 40 cycles of 95°C for 5 s, 60°C for 30 s, and finished with the default instrument dissociation protocol. Melt curve (dissociation) analysis was performed to identify the presence of primer–dimers and to analyze the specificity of the reaction. Amplicon sizes were verified for all primer pairs by agarose gel electrophoresis. Initial analysis of raw qPCR data utilized SDS v2.4.1 software (Life Technologies, USA).
Eleven samples were chosen for deep sequencing based on cost limitations while maintaining statistical relevance and the best possible spread of age and sex.
RNA isolation for library preparation
Prior to library preparation, 1.5 μg of each total RNA sample was DNase-treated using TURBO DNase (Ambion, USA) according to the manufacturer’s instructions followed by purification with Agencourt RNAClean XP beads (Beckman Coulter, USA) also according to manufacturer’s instructions. RNA concentration was measured using a Nanodrop 2000 spectrophotometer (Thermo Fisher Scientific, USA).
Five hundred nanograms of total RNA were used as input material for library preparation using the TruSeq Stranded Total RNA Sample Prep Kit (Illumina, USA) according to manufacturer’s instructions including the recommended modification to the RNA fragmentation duration to account for partially degraded RNA. Individual libraries were indexed as recommended by Illumina.
Quantification and quality control of DNA library templates
Indexed DNA libraries were analyzed individually using an Agilent Technologies 2100 Bioanalyzer with the DNA 1000 kit according to the manufacturer’s instructions (Agilent Technologies, USA). Libraries were diluted and pooled to a final concentration of 10 nM each in nuclease-free H2O (Ambion, USA). Pooled libraries were quantitated using a Life Technologies Qubit 2.0 Fluorometer with the Qubit dsDNA HS Assay Kit (Life Technologies, USA) and further diluted to 2 nM. Final DNA library concentration was confirmed using a Qubit dsDNA HS Assay Kit. PCR-competent library DNA concentration was verified using the universal KAPA Library Quantification Kit for Illumina Sequencing Platforms according to manufacturer’s instructions (KAPA Biosystems, USA).
Total RNA sequencing was performed using the Illumina HiSeq2500 platform with 100 bp paired-end sequencing with a fragment size of ~295 bp. Illumina TruSeq version 3 chemistry was used for cluster generation and sequencing.
Initial and post trimming quality control was performed with FastQC (version 0.10.1). Adapter and quality trimming was performed with TrimGalore [version 0.3.3; (30)] including adapter cutting, a minimal length of 20 bp and a Quality Phred score cutoff of 20. Trimmed paired-end reads were aligned against assembly GRCh37 of the human genome with TopHat [version 2.0.10; (31)] and Bowtie [version 2.1.0; (32)] using a pre-built transcriptome based on Gencode (version 19). Post processing and quality control of the alignment was performed with Bedtools [version 2.17.0; (33)], Samtools [version 0.1.19; (34)], Samstat [version 1.08; (35)], and RNA-SeQC [version 1.1.7; (36)]. The pipeline incorporating the aforementioned algorithms was built with NGSane [version 0.4.0; (37)]. Transcript normalization was performed using upper quartile normalization applied through the edgeR package (38). The average million reads per sample are ~16 million paired-end reads (16,106,463) and the average mapping of these reads is 94.38%. For transcriptome analysis, alignment counts were resolved on gene and transcript identifiers based on Gencode (version 19) annotation. Quantification was performed post alignment with HTSeq (version 0.5.4p5) in union mode considering only uniquely mapped reads. Transcriptome assembly was carried out with Cufflinks (39). Age regression analysis was performed using GraphPad Prism Version 6.0B with a linear regression fit and a 95% confidence interval with a p value of <0.05 deemed significant.
Genome-Wide Transcriptomic Analysis of the Adult Human SEZ
First, we wanted to confirm that proliferating cells within the cell cycle are likely to exist in the human SEZ. By qPCR we found that mRNA for the proliferation marker Ki67, required in all active phases of cell division while absent in the non-proliferating state, was present and that the levels of this transcript decreased sharply with increasing age (Figure 1). Next, to determine the molecular identity of lncRNAs in the adult human SEZ and to understand the mechanisms underpinning neurogenesis that may be applicable during adulthood, we used next generation sequencing (NGS) for genome-wide analysis. Transcript expression levels were ranked according to average HTSeq counts (expression counts; Table S2 and Figure S1 in Supplementary Material) over the 11 samples. Interrogating the top 100 expressed transcripts using DAVID 6.7 analysis for functional annotation (40) revealed highly significant cluster annotation (high stringency) for neuron development and differentiation (Enrichment score 6.18; Table S3 in Supplementary Material) providing confidence that this dataset may be valuable for investigating neurogenesis-related mechanisms in the human SEZ and supporting that this region is potentially enriched in immature neurons. We confirmed a decrease in neuronal differentiation in the human SEZ with advancing age as our NGS data revealed significant declines in the early neuronal differentiation and migration marker doublecortin [DCX; (41); Figure 2A; confirmed by qPCR (Figure S2A in Supplementary Material) in a larger cohort of samples (Table S1 in Supplementary Material)] and predicting additional mRNAs, which can be involved in proliferation, such as β-catenin [CTNNB1; (42); Figure 2B]. Interestingly, however, the neuronal progenitor marker paired box 6 [PAX6; (43); Figure 2C] and the early glial specification factor nuclear factor I A [NFIA; (44); Figure 2D], were predicted to be not significantly altered with age and trended upwards with increasing age implying that the progenitor pool for neurogenesis may persist during aging. These data suggest that significant neurogenesis may still be possible in the adult brain but proliferative levels of certain precursor cells may decline with regards to age.
Figure 1. Proliferation significantly decreases in the human SEZ during aging. Quantitative PCR reveals that mRNA for the proliferation marker Ki67 sharply decreases in an age-dependent manner. (Significance is defined as p value <0.05).
Figure 2. Next generation sequencing shows age-dependent decreases in early neuron differentiation and putative proliferation markers but stable expression of neural progenitor markers in the human SEZ. Significant decreases during aging are observed in the expression levels of (A) the immature neuronal marker doublecortin (DCX; validated by quantitative PCR; Figure S2A in Supplementary Material) and (B) a key promoter of neuronal proliferation β-catenin (CTNNB1). Non-significant age-dependent alterations are observed in the expression of radial glial markers (C) paired box 6 (PAX6) and (D) nuclear factor I A (NFIA). (Significance is defined as p value <0.05).
Increase in Inflammatory Markers in the Human SEZ during Aging
Increased brain inflammation with aging is thought to underlie, at least in part, the gradual decline of human brain function (45, 46). Our data predict that mRNA-encoding receptors of inflammatory mediators interleukin (IL)-1 and IL-6, IL1R, and IL6R, are significantly upregulated during aging in the human SEZ (Figures 3A,B). These changes were in contrast to the predicted decrease in proliferative (Figures 1 and 2B), progenitor (Figures 2C,D), and immature neuron (Figure 2A) markers, and could reflect indicators of brain injury, senescence, or neurogenic impairment. However, as our samples were from donors with no overt phenotype these increases are likely due to the involvement of IL-1 and IL-6 pathways in “normative” senescence (47).
Figure 3. Inflammatory receptor mRNAs related to senescence increase in the human SEZ with age. Significant increases during aging are observed in the expression levels of (A) interleukin-1 receptor 1 (IL1R1) and (B) interleukin-6 receptor (IL6R). (Significance is defined as p value <0.05).
A Subset of lncRNAs are Highly Expressed and Upregulated during Aging in the Human SEZ
To begin to characterize more fully the molecular signature of lncRNAs in the human SEZ, we determined the highest expressed lncRNAs from the human SEZ NGS data (Table 1). Expectedly, due to the generally low expression of lncRNAs in any particular brain region, we found that only 30 lncRNAs fell into the top 2262 of the most highly expressed transcripts, i.e., the 30th most highly expressed lncRNA was 2262nd (Table S2 in Supplementary Material). We chose the top 30 as they represent a reasonably high level of expression. The fairly well-known, widely and robustly expressed lncRNAs, such as MALAT1, involved in neuronal function (48) and proliferation (49), and NEAT1, implicated in paraspeckle formation and alternative splicing (50), were among the most highly expressed lncRNAs in the human SEZ. Additional lncRNAs that have been linked to differentiation were in the highly expressed group (Table 1) and were predicted to increase in expression during aging [Figures 4A–D; NEAT1 examined by qPCR increases with age; however, the relationship is not significant (Figure S2B in Supplementary Material) in a larger cohort of samples (Table S1 in Supplementary Material)]. For example, the lncRNA Gomafu is implicated in neuronal differentiation (18) and alternative splicing in human cortical neurons (22) and we find that Gomafu expression is suggested to increase in the human SEZ during aging. This is also true for the lncRNA TUG1 that is a direct target of the tumor suppressor p53 and reported to regulate proliferation (51).
Figure 4. Highly expressed lncRNAs that increase in expression during aging in the human SEZ. A significant increase during aging in the human SEZ is observed in the expression levels of (A) MALAT1. Insignificant, although upward trending, increases in expression during aging are also seen for the highly expressed lncRNAs (B) Gomafu, (C) NEAT1, and (D) TUG1. (Significance is defined as p value <0.05).
A Unique Subset of lncRNAs of Unknown Function Decreased during Aging in the Human SEZ
Our analysis of the SEZ reveals that there are highly expressed and uncharacterized lncRNAs that display significant regulation that parallel decreased age-related proliferation such as LINC00657 (Figure 5A) and SNORD3A (Figure 5B). The function of LINC00657 is unknown but this lncRNA appears to decline rapidly during aging (Figure 5A), a result validated through qPCR in a larger cohort of samples (see Figure S2C and Table S1 in Supplementary Material). This lncRNA is conserved in vertebrates (Figure S3 in Supplementary Material) and, although widely expressed, is enriched in the brain (Figure S3 in Supplementary Material). SNORD3A, an lncRNA encoding a small nucleolar RNA (C/D Box 3A), is predicted to be similarly downregulated during aging in the human SEZ (Figure 5B). This lncRNA is mammalian-specific (Figure S4 in Supplementary Material) and its expression levels are lower in brain tissue than elsewhere (Figure S4 in Supplementary Material) and may reflect expression specificity in particular cell types or regions, or dependence on activity in the human brain.
Figure 5. Examples of highly expressed and uncharacterized lncRNAs that are downregulated during aging in the human SEZ. Significant age-related decreases in the human SEZ are observed for the expression of the lncRNAs (A) LINC00657 (Validated by quantitative PCR; Figure S2C in Supplementary Material) and (B) SNORD3A. (Significance is defined as p value <0.05).
Possible Enhancer Transcript for DLX1 Detected in the Human SEZ
A lncRNA, antisense to the mouse Dlx1 gene (Dlx1as), may regulate Dlx1 function and hence interneuron specification but is proposed not to exist in humans (28). Interestingly, when we examined our human SEZ RNAseq data specifically examining the DLX locus, our analysis revealed reads that aligned to the syntenic region of Dlx1as. The position and structure of the human antisense transcript appears altered compared with mouse (Figures 6A,B). The total size of the antisense transcripts (excluding intronic regions) between mouse and human are comparable; however, the human transcript seems to have lost the DLX1-overlapping portion and may not possess an intron as seen in the mouse (at least in the adult human SEZ transcript pool). In a comprehensive study of enhancer regions, Andersson and colleagues designated this transcript as a possible actively transcribed enhancer (52). The putative human enhancer DLX1 (eDLX1) transcript is expressed at much lower levels in the SEZ than DLX1 (Figure 6C) but is of interest as neither ENCODE/GENCODE, RefSeq, UCSC (Figure S5 in Supplementary Material) nor the Illumina Body Atlas (data not shown) could find evidence of the human transcript emanating from the opposite DNA from the DLX1 sense strand. This may be due to the fact that it is uniquely expressed in the SEZ, or other discrete regions, and no other study has sequenced this region in depth. It remains to be seen whether this transcript indeed acts as an enhancer or traditional antisense regulator. Both DLX1 (Figure 6C) and eDLX1 (Figure 6D) are predicted to be expressed at relatively stable levels throughout the age range studied.
Figure 6. Detection of a human DLX1 antisense transcript that may function as an enhancer in the human SEZ. Graphical depiction of the corresponding (A) human (DLX1, antisense transcripts, and DLX2) and (B) mouse (Dlx1, Dlx1as, and Dlx2) regions of the genome. Both (C) DLX1 and (D) the putative DLX1 enhancer transcript expression are predicted to remain relatively constant during aging in the human SEZ. (Significance is defined as p value <0.05).
In the case of neurodevelopmental diseases, interneuron maldevelopment may underlie conditions such as schizophrenia and autism due to the fact that interneurons are specified and still developing well into postnatal life (7, 53–55) and accumulating evidence supports an inhibitory neuron dysfunction that may be common in both conditions and is, in fact, well replicated in schizophrenia (8, 56, 57). However, the role of interneurons in brain disorders first diagnosed in older adults or in normal brain aging is not as clear. Although different types of interneurons are present in the brain, a general deficit in interneuron progenitor pools could affect all interneuron subtypes or it may be that interneuronal specification potential changes with age such that abnormal molecular environments arising at a particular life stage could lead to distinct outcomes. As the human SEZ is a major source of interneuron production during development (1) and may contribute to adult neurogenesis, an understanding of the molecular control of this region is fundamental to our understanding of human brain diseases. In this study, we provide transcriptional evidence to suggest that interneuronal genesis may continue throughout human life in the SEZ, including the detection of mRNA encoding a molecule (DCX) known to be expressed by migrating interneurons (6, 58–60), and an mRNA encoding a transcription factor capable of directing cell fate toward the interneuron phenotype, DLX, and finally a putative DLX1 enhancer transcript, perhaps providing a window into unique transcripts that operate to induce proliferation, cell differentiation, and early neuronal migration even into advanced years of life in humans.
Our transcriptomic analysis of the human SEZ during adulthood revealed a bias for neuronal differentiation and an lncRNA “signature” that includes a small subset of highly expressed functionally known and unknown lncRNA transcripts. Interestingly, our data predict that the expression of some putative proliferative mRNA markers (Ki67, CNTTB1) and immature neuron markers (DCX) decrease during aging but intriguingly, progenitor markers (PAX6, NFIA) remain stable suggesting that the transition from multipotent progenitor cells to more committed precursors and neuroblasts (those in a transitional state) is reduced with age. An age-dependent decrease in neurogenesis seemed to coincide with the predicted increase in some of the highly expressed and proliferation-associated lncRNAs, such as Gomafu, MALAT1, and NEAT1 that form nuclear speckles and require activity to perform their functions. We did not anticipate this result, and it may be that during aging these lncRNAs become unresponsive to mitogenic or environmental stimulation or the stimuli themselves are diminished or not relayed correctly.
Several studies suggest that genesis of new neurons may be induced in the adult brain under trauma such as ischemia and injury (61). Injury or disease resulting in inflammation increases levels of pro-inflammatory cytokines such as IL-6, and we have previously reported that increased cortical cytokines are related to increased density of neurons in the white matter in adult humans, which can be interpreted as a possible increase in immature migratory neurons (56, 59). Furthermore, in this study, we predicted significant increases in IL1R and IL6R in the human SEZ with advancing age (Figure 3). These observations coupled with other reports of increased pro-inflammatory cytokines in brain aging (62) would suggest that we may expect to find increases in cell proliferation and early neuronal differentiation in the aging SEZ. However, this was not the case. While it may be that increases in inflammatory cytokines can be involved with inducing neurogenesis in a younger brain (63), inflammatory pathways are also implicated in senescence (47) and IL-6 can inhibit neuronal differentiation of neural stem cells (64) and astrocytic overexpression of IL-6 inhibits hippocampal neurogenesis (65). We find that the mRNA-encoding cytokine receptors for both IL1 and IL6 are suggested to increase in the aging human SEZ while the proliferation marker Ki67 is decreased. Our results would be consistent with a speculative model where receptors for cytokines may be involved in sensing and translating the demise of the human brain and stimulating senescence of the SEZ. It is also possible that although inflammation-induced neurogenesis is potentially viable in the young brain this system may become more sensitized to inflammation in the aging brain making it unable to respond to injury- or degeneration-induced inflammation in the same manner (66).
Our data suggest that neurogenesis in the adult human brain continues across the typical life span. However, normal levels of neurogenesis may be hampered in aging perhaps by an altered ability of resident SEZ progenitor cells to successfully initiate proliferation in response to extrinsic cues, i.e., up regulation of cytokines. Additionally, we have identified that many lncRNAs are expressed in this region, some of which may be unique. Future experiments aimed at a greater understanding of the role of lncRNAs in human neurogenesis and how cell intrinsic molecular factors may interact with cell extrinsic molecular cues in cell specific patterns over aging are needed to advance our understanding of how to harness and control new neuronal birth. Also, separating these possibilities of change in proliferative response of neuronal precursors to cytokines over aging would be of great relevance to the field. Furthermore, mechanistic insight into the function of highly expressed and uncharacterized lncRNAs in the human SEZ could uncover neurogenic regulatory processes central to this region. We believe that this study will inform future analyses of neurodevelopmental and neurodegenerative diseases and highlight the need to examine lncRNAs in larger sample sizes in future transcriptomic studies to reveal region-specific function.
Conflict of Interest Statement
Cynthia Shannon Weickert is a member of the advisory board for Lundbeck and consultant for Roche Pharmaceuticals. The other co-authors declare no conflict of interest.
This work was supported by the Schizophrenia Research Institute (utilizing infrastructure funding from the New South Wales Ministry of Health and the Macquarie Group Foundation), the University of New South Wales, and Neuroscience Research Australia. CW is a recipient of a National Health and Medical Research Council (Australia) Senior Research Fellowship (#1021970). The New South Wales Tissue Resource Centre at the University of Sydney is supported by the National Health and Medical Research Council of Australia, Schizophrenia Research Institute, and National Institute of Alcohol Abuse and Alcoholism (grant number R24AA012725). We gratefully acknowledge the Swiss National Science Foundation for the funding of BG (P2EZP3_152143).
The Supplementary Material for this article can be found online at http://www.frontiersin.org/Journal/10.3389/fneur.2015.00045/abstract
1. Curtis MA, Low VF, Faull RL. Neurogenesis and progenitor cells in the adult human brain: a comparison between hippocampal and subventricular progenitor proliferation. Dev Neurobiol (2012) 72:990–1005. doi: 10.1002/dneu.22028
3. Curtis MA, Kam M, Nannmark U, Anderson MF, Axell MZ, Wikkelso C. Human neuroblasts migrate to the olfactory bulb via a lateral ventricular extension. Science (2007) 315:1243–9. doi:10.1126/science.1136281
6. Wang C, Liu F, Liu YY, Zhao CH, You Y, Wang L. Identification and characterization of neuroblasts in the subventricular zone and rostral migratory stream of the adult human brain. Cell Res (2011) 21:1534–50. doi:10.1038/cr.2011.83
7. Weickert CS, Webster MJ, Colvin SM, Herman MM, Hyde TM, Weinberger DR. Localization of epidermal growth factor receptors and putative neuroblasts in human subependymal zone. J Comp Neurol (2000) 423:359–72. doi:10.1002/1096-9861(20000731)423:3<359::AID-CNE1>3.0.CO;2-0
9. Encinas JM, Michurina TV, Peunova N, Park JH, Tordo J, Peterson DA. Division-coupled astrocytic differentiation and age-related depletion of neural stem cells in the adult hippocampus. Cell Stem Cell (2011) 8:566–79. doi:10.1016/j.stem.2011.03.010
15. Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H. The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res (2012) 22:1775–89. doi:10.1101/gr.132159.111
16. Dinger ME, Amaral PP, Mercer TR, Pang KC, Bruce SJ, Gardiner BB. Long noncoding RNAs in mouse embryonic stem cell pluripotency and differentiation. Genome Res (2008) 18:1433–45. doi:10.1101/gr.078378.108
18. Mercer TR, Qureshi IA, Gokhan S, Dinger ME, Li G, Mattick JS. Long noncoding RNAs in neuronal-glial fate specification and oligodendrocyte lineage maturation. BMC Neurosci (2010) 11:14. doi:10.1186/1471-2202-11-14
22. Barry G, Briggs JA, Vanichkina DP, Poth EM, Beveridge NJ, Ratnu VS. The long non-coding RNA Gomafu is acutely regulated in response to neuronal activation and involved in schizophrenia-associated alternative splicing. Mol Psychiatry (2014) 19:486–94. doi:10.1038/mp.2013.45
24. Ng SY, Johnson R, Stanton LW. Human long non-coding RNAs promote pluripotency and neuronal differentiation by association with chromatin modifiers and transcription factors. EMBO J (2012) 31:522–33. doi:10.1038/emboj.2011.459
26. Cobos I, Calcagnotto ME, Vilaythong AJ, Thwin MT, Noebels JL, Baraban SC. Mice lacking Dlx1 show subtype-specific loss of interneurons, reduced inhibition and epilepsy. Nat Neurosci (2005) 8:1059–68. doi:10.1038/nn1499
29. Weickert CS, Sheedy D, Rothmond DA, Dedova I, Fung S, Garrick T. Selection of reference gene expression in a schizophrenia brain cohort. Aust N Z J Psychiatry (2010) 44:59–70. doi:10.3109/00048670903393662
36. DeLuca DS, Levin JZ, Sivachenko A, Fennell T, Nazaire MD, Williams C. RNA-SeQC: RNA-seq metrics for quality control and process optimization. Bioinformatics (2012) 28:1530–2. doi:10.1093/bioinformatics/bts196
37. Buske FA, French HJ, Smith MA, Clark SJ, Bauer DC. NGSANE: a lightweight production informatics framework for high-throughput data analysis. Bioinformatics (2014) 30:1471–2. doi:10.1093/bioinformatics/btu036
38. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics (2010) 26:139–40. doi:10.1093/bioinformatics/btp616
39. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol (2010) 28:511–5. doi:10.1038/nbt.1621
41. Couillard-Despres S, Winner B, Schaubeck S, Aigner R, Vroemen M, Weidner N. Doublecortin expression levels in adult brain reflect neurogenesis. Eur J Neurosci (2005) 21:1–14. doi:10.1111/j.1460-9568.2004.03813.x
44. Piper M, Barry G, Hawkins J, Mason S, Lindwall C, Little E. NFIA controls telencephalic progenitor cell differentiation through repression of the Notch effector Hes1. J Neurosci (2010) 30:9127–39. doi:10.1523/JNEUROSCI.6167-09.2010
45. Cribbs DH, Berchtold NC, Perreau V, Coleman PD, Rogers J, Tenner AJ. Extensive innate immune gene activation accompanies brain aging, increasing vulnerability to cognitive decline and neurodegeneration: a microarray study. J Neuroinflammation (2012) 9:179. doi:10.1186/1742-2094-9-179
47. Orjalo AV, Bhaumik D, Gengler BK, Scott GK, Campisi J. Cell surface-bound IL-1alpha is an upstream regulator of the senescence-associated IL-6/IL-8 cytokine network. Proc Natl Acad Sci U S A (2009) 106:17031–6. doi:10.1073/pnas.0905299106
48. Bernard D, Prasanth KV, Tripathi V, Colasse S, Nakamura T, Xuan Z. A long nuclear-retained non-coding RNA regulates synaptogenesis by modulating gene expression. EMBO J (2010) 29:3082–93. doi:10.1038/emboj.2010.199
49. Tripathi V, Shen Z, Chakraborty A, Giri S, Freier SM, Wu X. Long noncoding RNA MALAT1 controls cell cycle progression by regulating the expression of oncogenic transcription factor B-MYB. PLoS Genet (2013) 9:e1003368. doi:10.1371/journal.pgen.1003368
51. Zhang EB, Yin DD, Sun M, Kong R, Liu XH, You LH. P53-regulated long non-coding RNA TUG1 affects cell proliferation in human non-small cell lung cancer, partly through epigenetically regulating HOXB7 expression. Cell Death Dis (2014) 5:e1243. doi:10.1038/cddis.2014.201
53. Fung SJ, Webster MJ, Sivagnanasundaram S, Duncan C, Elashoff M, Weickert CS. Expression of interneuron markers in the dorsolateral prefrontal cortex of the developing human and in schizophrenia. Am J Psychiatry (2010) 167:1479–88. doi:10.1176/appi.ajp.2010.09060784
56. Fung SJ, Fillman SG, Webster MJ, Shannon Weickert C. Schizophrenia and bipolar disorder show both common and distinct changes in cortical interneuron markers. Schizophr Res (2014) 155:26–30. doi:10.1016/j.schres.2014.02.021
57. Joshi D, Fullerton JM, Weickert CS. Elevated ErbB4 mRNA is related to interneuron deficit in prefrontal cortex in schizophrenia. J Psychiatr Res (2014) 53:125–32. doi:10.1016/j.jpsychires.2014.02.014
58. Francis F, Koulakoff A, Boucher D, Chafey P, Schaar B, Vinet MC. Doublecortin is a developmentally regulated, microtubule-associated protein expressed in migrating and differentiating neurons. Neuron (1999) 23:247–56. doi:10.1016/S0896-6273(00)80777-1
59. Fung SJ, Joshi D, Allen KM, Sivagnanasundaram S, Rothmond DA, Saunders R. Developmental patterns of doublecortin expression and white matter neuron density in the postnatal primate prefrontal cortex and schizophrenia. PLoS One (2011) 6:e25194. doi:10.1371/journal.pone.0025194
60. Gleeson JG, Lin PT, Flanagan LA, Walsh CA. Doublecortin is a microtubule-associated protein and is expressed widely by migrating neurons. Neuron (1999) 23:257–71. doi:10.1016/S0896-6273(00)80778-3
Keywords: long non-coding RNA, subependymal zone, aging, neurodevelopmental disease, neurogenesis, interneuron
Citation: Barry G, Guennewig B, Fung S, Kaczorowski D and Weickert CS (2015) Long non-coding RNA expression during aging in the human subependymal zone. Front. Neurol. 6:45. doi: 10.3389/fneur.2015.00045
Received: 29 August 2014; Accepted: 21 February 2015;
Published online: 09 March 2015.
Edited by:Owen Murray Rennert, National Institutes of Health, USA
Reviewed by:Khyobeni Mozhui, University of Tennessee Health Science Center, USA
Tin-Lap Lee, The Chinese University of Hong Kong, Hong Kong
Copyright: © 2015 Barry, Guennewig, Fung, Kaczorowski and Weickert. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guy Barry, Garvan Institute of Medical Research, 384 Victoria Street, Darlinghurst, NSW 2010, Australia e-mail: firstname.lastname@example.org