Abstract
Background:
Motor and cognitive decline as part of the normal aging process is linked to alterations in synaptic plasticity and reduction of adult neurogenesis in the dorsal striatum. Neuroinflammation, particularly in the form of microglial activation, is suggested to contribute to these age-associated changes.
Objective and Methods:
To explore the molecular basis of alterations in striatal function during aging we analyzed RNA-Seq data for 117 postmortem human dorsal caudate samples and 97 putamen samples acquired through GTEx.
Results:
Increased expression of neuroinflammatory transcripts including TREM2, MHC II molecules HLA-DMB, HLA-DQA2, HLA-DPA1, HLA-DPB1, HLA-DMA and HLA-DRA, complement genes C1QA, C1QB, CIQC and C3AR1, and MHCI molecules HLA-B and HLA-F was identified. We also identified down-regulation of transcripts involved in neurogenesis, synaptogenesis, and synaptic pruning, including DCX, CX3CL1, and CD200, and the canonical WNTs WNT7A, WNT7B, and WNT8A. The canonical WNT signaling pathway has previously been shown to mediate adult neurogenesis and synapse formation and growth. Recent findings also highlight the link between WNT/β-catenin signaling and inflammation pathways.
Conclusions:
These findings suggest that age-dependent attenuation of canonical WNT signaling plays a pivotal role in regulating striatal plasticity during aging. Dysregulation of WNT/β-catenin signaling via astrocyte–microglial interactions is suggested to be a novel mechanism that drives the decline of striatal neurogenesis and altered synaptic connectivity and plasticity, leading to a subsequent decrease in motor and cognitive performance with age. These findings may aid in the development of therapies targeting WNT/β-catenin signaling to combat cognitive and motor impairments associated with aging.
Keywords
Introduction
Normal aging is associated with a progressive decline in motor and cognitive functions. The dorsal striatum, which includes the caudate nucleus and putamen, receives dopaminergic input from the substantia nigra and glutamatergic input from the cortex and thalamus (Do et al., 2012). The dorsal striatum is involved in controlling motor movement and action-contingent learning (Delgado, Miller, Inati, & Phelps, 2005; Delgado, Nystrom, Fissell, Noll, & Fiez, 2000; Haber, 2016; Haruno et al., 2004; Kreitzer & Malenka, 2008; Tricomi, Delgado, & Fiez, 2004; C. Wu, Garamszegi, Xie, & Mash, 2017). Altered striatal plasticity has been implicated in the age-related deterioration in motor performance indicators such as balance deficits and slowed movement (Cham, Perera, Studenski, & Bohnen, 2007; Seidler et al., 2010; van Dyck et al., 2008), and cognitive impairment such as deficits in working memory, episodic memory, and executive function (Erixon-Lindroth et al., 2005; Klostermann, Braskie, Landau, O’Neil, & Jagust, 2012). The molecular mechanisms involved in the process of aging, however, remain poorly understood.
Many studies have demonstrated structural brain changes in normal aging (Fjell & Walhovd, 2010). MRI studies showed longitudinal shrinkage of the caudate nucleus and putamen in healthy adults (Gunning-Dixon, Head, McQuain, Acker, & Raz, 1998; Raz et al., 2003). A study in aged cats (Levine et al., 1986) found decreased spine density and dendritic length in medium-sized spiny neurons of the caudate nucleus. Recently, using 14C birth-dating, adult neurogenesis has been verified in the adult human subventricular zone (SVZ) of the lateral ventricle wall and the subgranular zone (SGZ) of the hippocampal dentate gyrus (Ernst et al., 2014; Spalding et al., 2013). In contrast to rodents, SVZ-derived neuroblasts in humans migrated only to the caudate nucleus and putamen, but not to the olfactory bulb (Ernst et al., 2014). Since adult-generated striatal neurons are preferentially destroyed in the neurodegenerative Huntington’s disease (HD) and adult neurogenesis in the human striatum has been shown to decline during aging (Ernst et al., 2014), age-associated reductions in adult striatal neurogenesis, to the extent that those adult-generated neurons are inadequate for normal striatal function, may contribute to neurodegenerative disorders.
Rodent experiments have suggested that inflammation plays a role in age-associated changes in adult neurogenesis and neuroplasticity, and these often relate to aberrant microglial activation (Gebara, Sultan, Kocher-Braissant, & Toni, 2013; L’Episcopo et al., 2013; Zhang et al., 2013). Notably, up-regulation of genes reflecting microglial activation has been reported in the aging human brain (Cribbs et al., 2012). Microglia, the resident immune cells of the brain, survey and shape neural circuits through phagocytosis of unnecessary synaptic elements and newborn neurons in the developing and adult brain (Brown & Neher, 2014; Paolicelli & Gross, 2011; Sierra et al., 2010; Wake, Moorhouse, Miyamoto, & Nabekura, 2013). On the other hand, inflammatory activation of microglia has been implicated in the pathogenesis of various neurological disorders (Frick, Williams, & Pittenger, 2013; Mondelli, Vernon, Turkheimer, Dazzan, & Pariante, 2017; Mosher & Wyss-Coray, 2014; Wyss-Coray, 2016).
In aged mice, activated microglia may contribute to the decline of SVZ neurogenesis through the up-regulation of microglial pro-inflammatory mediators (L’Episcopo et al., 2013). Recent studies indicate that astrocytes respond to microglia-derived inflammatory molecules by increasing expression of secreted WNT ligands which act to restrain microglial activation by attenuating pro-inflammatory cytokine overexpression (L’Episcopo et al., 2014; L’Episcopo et al., 2011; L’Episcopo et al., 2012; Marchetti & Pluchino, 2013). Microglia overactivation in aged mice has, however, been shown to impair the neurogenic capacity of neural progenitor cells (NPCs) in the SVZ through aberrant astrocyte–microglia crosstalk and subsequent down-regulation of WNT/β-catenin signaling in both astrocytes and microglia (L’Episcopo et al., 2012; L’Episcopo et al., 2013). Nevertheless, the role of WNT-glial interactions in the aged human brain is not known.
To identify age-dependent changes in gene expression networks we performed differential expression and pathway enrichment analyses on 117 postmortem human caudate and 97 putamen samples. These samples were from the Genotype-Tissue Expression (GTEx) database, a NIH founded project that provides extensive gene expression and clinical data from a large cohort of postmortem tissues unbiased for any diseases (Consortium, 2013, 2015; Mele et al., 2015).
Materials and methods
GTEx database and pre-processing of RNA-Sequencing dataset
The GTEx database (v6) collected 1632 brain samples from 13 brain regions. Detailed information on sample collection, RNA sequencing, and data processing pipelines can be found in the GTEx Consortium paper (Mele et al., 2015). Specifically, only optimal (SMTORMVE not equal to “FLAGGED”) TrueSeq.v1 (SMGEBTCHT = “TrueSeq.v1”) RNA-Seq samples in the caudate (N = 117) and putamen (N = 97) were selected for this study (Table S1). Gene expression datasets encompassing raw read counts and RPKM (Reads Per Kilobase of transcript per Million mapped reads) were generated using RNA-SeQC. In addition, GTEx dataset with GENECODE v19 as a reference transcriptome resulted in 56,318 annotated genes.
We evaluated similarities in gene expression between the caudate and putamen samples for each individual by performing hierarchical clustering and inter-array correlation (IAC) (Oldham et al., 2008). Read counts were transformed using the voom function in Limma before IAC analysis (Law, Chen, Shi, & Smyth, 2014). Relationships among RNA-Seq samples were measured using average linkage hierarchical clustering with a distance metric of 1-IAC.
Identification of significant age-associated genes
Genes differentially expressed based on age were identified using a procedure similar to that described previously (C. Wu, Bendriem, Garamszegi, Song, & Lee, 2017; C. Wu, Chen, Shu, & Lee, 2017; C. Wu, Xu, Tsai, Freed, & Lee, 2017). We performed differential gene expression analysis with DESeq2 (v1.12.3), and considered genes to be expressed if they had more than five reads in at least 25%of total samples (Love, Huber, & Anders, 2014). Next, we used the “svaseq” function from the sva R package to identify hidden batch effects and cofounding factors (Leek, 2014). In addition to gender, 15 surrogate variables were added to the formula in DESeq2. In either caudate or putamen samples, the residuals for normalized read counts, after gender and surrogate variables correction, were tested against continuous age using the following negative binomial (NB) generalized linear regression model (GLM):
where Kij is the read counts for gene j in sample i, fitted with a negative binomial distribution. αj is a gene-specific dispersion parameter. μij represents fitted mean, containing a sample-specific size factor si and a covariate-dependent part qij (Love et al., 2014). In Equation 3, β0 is the regression intercept for gene j, ɛij is the error term. β1j, β2j, and βkj (k = 1, . . . , 15) denote the regression coefficients of age, gender and kth surrogate variables for gene j respectively.
The Wald test P-values were adjusted for multiple testing using the Benjamini-Hochberg (BH) algorithm. Linear association with the continuous “AGE” variable was used as the primary strategy for detecting age-associated genes. To further validate our results, we arbitrarily divided RNA-Seq samples into three age groups: young (AGE≤55), middle-aged (55 < AGE < 65), and old (AGE≥65). The same formula shown in Equation 3 above was applied with the exception of the primary variable being changed to categorical age groups (“AGEGROUP”) with three-factor levels:
The list of differentially expressed genes that differed between old- and young-age groups was compared to that of the age-associated genes derived from linear association with continuous age. The regression coefficients β1j indicate the log2 fold changes for gene j from the contrast between the age groups, whereas log2 fold changes for continuous age means change per unit (year) (Love et al., 2014). The threshold for significance of differential expression of genes between old- and young-age groups was at a false discovery rate (FDR) level of 0.2. We defined significant age-associated genes at FDR < 0.2 and uncorrected P-value < 0.01 from linear association with continuous age as the primary strategy.
Gene set enrichment analysis
Gene set enrichment analysis (GSEA) was performed by using clusterProfiler (v3.0.2) in R (Yu, Wang, Han, & He, 2012), which implements the algorithm developed by the Broad Institute (Subramanian et al., 2005). Specifically, we constructed a pre-ranked gene list of all expressed genes ordered by log2 fold changes for continuous age (beta regression coefficients) from DESeq2 package. Enrichment score (ES) and the significance level of ES (nominal P-value) were calculated by 1000 phenotype-based permutation test. Pre-defined pathways from Kyoto Encyclopedia of Genes and Genomes (KEGG PATHWAY database) with minimal gene size of 10 and maximal of 500 were used in GSEA. q-values were calculated for FDR control. Significant pathways with q-value less than 0.05 were reported.
Overrepresentation test
We performed overrepresentation tests on 1832 age-associated genes by using meta-analysis. The function “enrichKEGG” or “enrichGO” of the cluster Profiler package was used to test whether up- or down-regulated age-associated genes are significantly overrepresented in given pathways from KEGG or gene ontology (GO) categories (Biological Process, BP). Hypergeometric distribution was used to evaluate the P-value. q-values were reported for FDR control. For pathway-based analyses, we defined significant age-associated biological functions at the level of q-value less than 0.05. Selected significant age-associated KEGG pathways were visualized through “pathview” R package (W. Luo & Brouwer, 2013).
Expression correlation analysis of synaptic and immune genes
To evaluate the gene-gene correlations within synaptic and immune gene sets, Pearson correlation matrix based on gene expression in RPKM was constructed and displayed by corrplot package in R. The function “rcorr” from Hmisc package was used to compute the significance levels for Pearson correlations. Any correlation that did not pass the significance level (P-value < 0.01) was left blank.
Replication by using human putamen dataset from UKBEC and meta-analysis
The UK Human Brain Expression Consortium (UKBEC) dataset is the largest control brain gene expression dataset generated from microarray analysis. This dataset contains 129 putamen samples from Caucasians aged 16–102 years (Table S1). A detailed description of tissue collection and dissection, RNA isolation, and processing of brain samples was described in Trabzuni’s study (Trabzuni et al., 2011). Briefly, the expression profiles of postmodern brain tissues were performed on the Affymetrix GeneChip® Human Exon 1.0 ST Arrays (Affymetrix, High Wycombe, UK). The gene-level expression dataset (GSE60862) was already normalized using Robust Multichip Analysis (RMA) with GC background correction. We extracted 129 putamen samples with 26493 gene signatures. For each gene signature, we performed linear regression of age against gender and surrogate variable corrected expression values by using limma R package (Leek & Storey, 2007; Ritchie et al., 2015):
Where yij is RMA-normalized gene expression value for gene j in sample i. The “logFC” for continuous age variable produced by limma package is the beta regression coefficient (β1j). Multiple testing P-values were adjusted using the BH method.
To identify common age-associated genes in the dorsal striatum, we performed a meta-analysis by using P-value combination approaches (Fisher and inverse normal) to increase detection power (Rau, Marot, & Jaffrezic, 2014). All significant genes (FDR < 0.05) without differential expression conflicts were reported.
All statistical computing was performed in the R (v3.3, https://www.r-project.org/). Distribution of clinical characteristics and Pearson/Spearman correlation were performed with JMP v12.2.0 (SAS).
RNA-Seq dataset can be downloaded from dbGaP, study accession no. phs000424.v6.p1. Details about the GTEx resource can be found in the GTEx portal (http://www.gtexportal.org). Microarray dataset from UKBEC can be downloaded from GEO, accession no. GSE60862.
Results
Identifying significant age-associated genes in the human caudate and putamen
The inter-array correlation (IAC) analysis showed that RNA-Seq samples from either the caudate (mean IAC = 0.93) or putamen (mean IAC = 0.927) were highly correlated (Fig. S1A-D). This confirms an overall consistency amongst samples in each of the two brain regions.
Genes that were differentially expressed between young- (AGE≤55) and old-age (AGE≥65) groups with an FDR level of 0.2 were compared with age-associated genes derived from linear association with continuous age (Fig. 1B). A large overlap of the age-associated genes between these two strategies was observed (Fig. 1B), indicating consistency between the results of the two analyses. Since group comparison ignores variances within each age group, we opted to use the linear association with continuous age strategy for the subsequent analyses. Volcano plots (Fig. 1C, D) display significant age-associated gene expression changes in the caudate and putamen. At the level of FDR < 0.2 and P-value < 0.01, we identified 2494 age-associated genes in the caudate and 1004 age-associated genes in the putamen. Lists of these genes are presented in Tables S2 and S3. A potential explanation for the larger list of age-associated genes in the caudate samples compared to those of the putamen is the small sample size and lower IAC of the putamen dataset.

Transcriptome analysis reveals a large number of significant age-associated genes in caudate and putamen. (A) Histograms showing age distribution of subjects in caudate and putamen. (B) Shared significant age-associated genes between two different data analysis strategies. Red circles represent age-associated genes derived from linear association with continuous age (FDR < 0.2 and uncorrected P-value < 0.01), and green circles show the significant differentially expressed (DE) genes in old age group contrasted with young age group (FDR < 0.2). (C, D) Volcano plots show the -log10 (P-value) and beta regression coefficient of age for all expressed genes in the caudate (C) and putamen (D). Age-associated genes at the FDR level of 0.05, 0.1, and 0.2 are indicated by red, green, and yellow, respectively.
The dorsal striatum houses the caudate and putamen which have similar organization of their neuronal circuitry (Do, Kim, Bakes, Lee, & Kaang, 2012). The question then arises whether age-related gene expression changes in these two regions are coordinated. We have shown that caudate and putamen share 472 common significant age-associated genes (Fig. 2A), suggesting a substantial overlap between the caudate and putamen. Moreover, beta regression coefficients for age were highly correlated (Pearson r = 0.56) between the caudate and putamen (Fig. 2B). A negative regression coefficient of age for gene j signifies that for every unit (year) increase, we expect a beta unit decrease in gene expression (log2 normalized read counts; Equation 3), holding all other variables constant. Amongst a total of 472 shared genes, all directions of age-related changes (beta coefficient) were identical between the caudate and putamen, wherein 255 were down-regulated and 217 were up-regulated (Table S4). These results suggest consistent age-related gene expression changes between these two brain regions.

Correlation of age-related gene expression changes between caudate and putamen. (A) Venn diagram showing 472 shared significant age-associated genes between caudate and putamen. (B) Each point shows the beta regression coefficient of age in both brain regions. Both the caudate and putamen showing a high correlation in age-related gene expression changes (P-value < 0.0001 and Pearson’s r = 0.56).
Gene Set Enrichment Analysis (GSEA) revealed a total of 65 and 21 significant pathways with a q-value less than 0.05 in the caudate and putamen, respectively (Tables S5 and S6). Notably, 86%(18 out of 21) of enriched pathways in the putamen overlapped with those in the caudate (Fig. 3A). All directions of these age-associated pathway changes (Normalized Enrichment Score) were identical. This substantial overlap between the caudate and putamen indicates that although expression changes are weaker in the putamen, they follow a trend similar to those identified in the caudate.

Gene set enrichment analysis (GSEA) reveals a high similarity between the caudate and putamen in age-associated pathways. (A) Venn diagram showing the overlap of the significant age-associated KEGG pathways (FDR q-value < 0.05) between the caudate and putamen. (B) Top 10 significant age-associated KEGG pathways in the caudate. (C) Common significant age-associated pathways between the caudate and putamen.
Amongst eighteen pathways found in both the caudate and putamen, the pathways which were up-regulated throughout aging pertained to functions of the immune system and inflammation, including asthma, allograft rejection, autoimmune thyroid disease, staphylococcus aureus infection, viral myocarditis, complement and coagulation cascades, graft-versus-host disease, and antigen processing and presentation (Fig. 3C). The pathways which were down-regulated throughout the aging process were involved in proper synapse maintenance and function, including glutamatergic synapse, morphine addiction, retrograde endocannabinoid signaling, GABAergic synapse, and cholinergic synapse (Fig. 3C). These results are consistent with previous gene expression studies in the aging human brain (Mohan, Mather, Thalamuthu, Baune, & Sachdev, 2016; Sibille, 2013).
Notably, the RAP1 signaling pathway was down-regulated throughout aging in both caudate and putamen (Fig. 3C). RAP1 is essential in maintaining genome stability by protecting telomeric DNA ends from being shortened or altered via non-homologous end joining or homologous recombination (Bae & Baumann, 2007; Martinez et al., 2010; Sarthy, Bae, Scrafford, & Baumann, 2009; Sfeir, Kabir, van Overbeek, Celli, & de Lange, 2010). Telomeres are naturally shortened with age and telomere length has been previously linked to life expectancy in humans (Heidinger et al., 2012; Shammas, 2011). RAP1 down-regulation throughout aging; therefore, suggests that RAP1 signaling may partially drive observed age-related changes in the brain.
Recent evidence using histological and 14C dating approaches indicated that the turnover rates of neuronal populations in the striatum of the adult human brain decline with age (Ernst et al., 2014). In agreement with these findings, we found that expression of DCX, a marker of newly born neurons in the adult brain, was decreased in both caudate and putamen (Table S4), suggesting that striatal neurogenesis decreases in the aging human brain. WNT signaling has been identified as an essential regulatory pathway in adult neurogenesis and synaptic plasticity (Inestrosa & Arenas, 2010; Lie et al., 2005). Importantly, abnormal crosstalk between the WNT signaling and inflammatory response has been shown to impact SVZ plasticity and cause nigrostriatal injury (L’Episcopo et al., 2012; L’Episcopo et al., 2013). The WNT signaling pathway was found to be down-regulated throughout aging and is one of the top enriched pathways in the caudate (Fig. 3B, Table S5). GSEA in the putamen also showed a down-regulated WNT signaling pathway (P-value = 0.02, q-value = 0.11), but it was not included in our final list due to stringent cutoff criteria (Table S6). Moreover, the PI3K-Akt signaling pathway, which has been shown to positively modulate WNT signaling (Fukumoto et al., 2001; G. Lee et al., 2010), was down-regulated in the putamen (Table S6), further supporting the claim that altered WNT signaling participates in the aging process in the striatum.
We next used the “pathview” R package to visualize molecular networks that cause progressive loss of functional integrity of the striatum during the aging process (W. Luo & Brouwer, 2013). As shown in Fig. 4, expression of WNT signaling genes revealed a concordant decrease (negative regression coefficient) in both caudate and putamen. Moreover, concordant repression of age-associated GABAergic and glutamatergic synapse pathway genes was observed in both caudate and putamen (Fig. 5, Fig. S2). The putamen exhibited a slightly lower concordance decrease in expression of age-associated dopaminergic synapse pathway genes than those of the caudate (Fig. S3). Nevertheless, a concordant age-associated decrease in expression of genes required for long-term potentiation in both caudate and putamen (Fig. S4) is consistent with loss of synaptic function during aging.

Gene expression changes (beta regression coefficients) in the WNT signaling pathway visualized using Pathview. Expression of WNT signaling pathway genes showing concordant negative associations with age in the caudate (A) and putamen (B). Up-regulated genes (positive regression coefficient for age) appear in red and down-regulated genes (negative regression coefficient for age) appear in green.

Gene expression changes (beta regression coefficients) in the GABAergic synapse pathway visualized using Pathview. Expression of GABAergic synapse pathway genes showing concordant negative associations with age in the caudate (A) and putamen (B). Up-regulated genes (positive regression coefficient for age) appear in red, and down-regulated genes (negative regression coefficient for age) appear in green.
Expression of antigen processing and presentation genes, including both class I and class II MHC molecules (MHCI and MHCII, respectively), showed a concordant increase in both caudate and putamen (Fig. 6), suggesting activation of immune responses. Correlation analysis of expression of age-related synaptic and immune genes showed a highly coordinated network within functionally similar gene groups in both caudate (Fig. S5) and putamen (Fig. S6), implicating an association between impaired synaptic function and neuroinflammation in aging.

Gene expression changes (beta regression coefficients) in the antigen processing and presentation pathway visualized using Pathview. Expression of antigen processing and presentation pathway genes showing concordant positive associations with age in the caudate (A) and putamen (B). Up-regulated genes (positive regression coefficient for age) appear in red, and down-regulated genes (negative regression coefficient for age) appear in green.
To validate the results, we compared the RNA-seq data to the UKBEC microarray-based gene expression dataset for the putamen. The UKBEC dataset contains 129 putamen samples from Caucasian individuals aged 16–102 (Fig. 7A) (Trabzuni et al., 2011). We performed a similar linear regression approach based on the concepts previously described to analyze gene expression changes associated with age (Equation 5). Our results consistently showed that a similarly large number of age-associated genes (2462) were identified at the level of FDR < 0.2 and uncorrected P-value < 0.01 (Fig. 7B). All shared expressed genes in three datasets of UKBEC-putamen, GTEx-putamen, and GTEx-caudate (n = 15382) were merged by gene symbols. Beta regression coefficients of age for all 15382 genes were significantly correlated between these three datasets (p < 0.0001, Fig. S7). Moreover, the Venn diagrams (Fig. 7C, 7D) show an overlap of the significant age-associated genes and pathways amongst these three datasets. Notably, both WNT and PI3K-Akt signaling pathways were down-regulated in the putamen using the UKBEC dataset (Table S7), suggesting consistent age-related changes in WNT signaling in the human dorsal striatum across the three different datasets.

Validation using microarray-based UKBEC putamen gene expression dataset. (A) Age distribution of 129 donors from UKBEC. (B) Volcano plot shows the -log10 (P-value) and beta regression coefficient of age for all expressed genes in the putamen from UKBEC. Age-associated genes at the FDR level of 0.05, 0.1, and 0.2 are indicated by red, green, and yellow, respectively. (C) Three-way Venn diagram showing shared significant age-associated genes among three datasets (FDR < 0.2 and uncorrected P-value < 0.01). (D) Three-way Venn diagram showing shared significant age-related KEGG pathways among three datasets using GSEA (q-value < 0.05). (E) Five-way Venn diagram comparing the list of age-associated genes at a 5%BH threshold using three datasets and two P-value combination methods. (F) Each point shows the beta regression coefficient of age for all significant age-associated genes from the meta-analysis (n = 1832). A correlation probability matrix showing that three datasets are highly correlated. Histograms showing distributions of beta coefficients. Red lines represent the linear fit and 95%density contours.
We next performed a meta-analysis of age-related gene expression profiles on three datasets of UKBEC-putamen, GTEx-putamen, and GTEx-caudate using P-value combination approaches (Fisher and inverse normal) to identify shared molecular changes in the aging human dorsal striatum (Fig. 7E, S8). Significant age-associated genes at a 5%FDR threshold by meta-analysis were identified (Table S8). Among these 1832 genes, 1390 were negatively associated with age and 442 were positively associated with age. Figure 7F showed that beta regression coefficients of age for all significant age-associated genes from the meta-analysis (n = 1832) were significantly correlated between these three datasets (p < 0.0001). Notably, there was an age-related decrease in expression of the marker for newborn neurons DCX amongst the three datasets, further suggesting a reduction in newly-generated neurons (Table S8).
Overrepresentation tests for up- and down-regulated genes using GSEA revealed similar transcriptional signatures to those from the GTEx dataset. As shown in Table S9, down-regulated aging genes were significantly enriched in the PI3K-Akt signaling pathway and pathways acting on synaptic communication in the brain, whereas up-regulated genes were mainly enriched in immune and inflammatory pathways. Overrepresentation tests based on GO biological process database provided additional insights into biological functions of these significant age-associated genes from the meta-analysis. Among the down-regulated genes, the most enriched categories were cognition, learning or memory, WNT signaling, neurogenesis, neuronal migration, axonal and synaptic plasticity, and synaptic transmission. Among the up-regulated genes, the most enriched categories were antigen processing and presentation, and T-cell activation and aggregation (Table S10). These results showed consistent aging-related gene expression profiles amongst three datasets.
We identified a comprehensive list of aging-related genes and pathways in two human dorsal striatal regions, caudate and putamen. Genes pertaining to neurogenesis, synaptic plasticity, and striatal functions were significantly down-regulated, whereas genes within immune response and inflammation pathways were significantly up-regulated. Notably WNT signaling, which functions at the neuroimmune interface and plays pivotal roles in regulating neural progenitors, post-mitotic neurons, and microglia in the adult human brain, was down-regulated in association with aging. Other prior studies have examined gene expression patterns in aged human brain, but none have focused specifically on the dorsal striatum (Dillman et al., 2017; Donertas et al., 2017). The present study identifies a potential link between WNT signaling and age-mediated changes in neurogenesis, synaptic function, and immune cell activation in the dorsal striatum.
Age-mediated down-regulation of neurogenesis
Adult neurogenesis is the continuous generation and integration of neurons into existing neural circuitry. These adult-born neurons are structurally plastic and functional. Adult hippocampal neurogenesis has been implicated in cognitive processes such as hippocampal-dependent learning, memory, and mood regulation (Anacker & Hen, 2017; Goncalves, Schafer, & Gage, 2016). Age-dependent impairment of adult neurogenesis in the dentate gyrus (DG) has been linked to hippocampal dysfunction, and counteracting this aging-induced pathophysiology can reduce age-related cognitive decline (Casadesus et al., 2004; Drapeau et al., 2003; Driscoll et al., 2006; Mayo et al., 2003; Sun, Evans, Hsieh, Panici, & Bartke, 2005). The full spectrum of functions of striatal neurogenesis remains unknown; yet, the lower turnover rate of adult-generated striatal neurons in HD patients suggests that the interruption of striatal neurons turnover might be a cause of the motor and cognitive deficits that these patients exhibit (Ernst et al., 2014). It is therefore reasonable to assume that newly generated striatal neurons have a functional role in the striatum. Our results indicate an age-dependent reduction in the neuronal marker DCX transcript in both caudate and putamen (Table S4) and significant enrichment of down-regulated genes in neurogenesis-related pathways (Table S10), both of which are consistent with reports of a decline in striatal neuronal cell turnover during human aging (Ernst et al., 2014). The above observations suggest that reduced striatal neurogenesis in the aged brain relates to a gradual decay in the functionality of the dorsal striatum.
Age-induced changes in inflammatory and synaptic pathways
Correlation analysis of expression of age-mediated synaptic and immune response genes implicates the association of synaptic deficits and immune activation with aging in both caudate and putamen (Fig. S5, S6). The neuroinflammatory response following microglial activation is thought to contribute to age-related brain plasticity impairments (Dickson et al., 1992; X. G. Luo, Ding, & Chen, 2010). Microglia can be neuroprotective in the developing brain but have also been shown to become neurotoxic and destructive in the aged brain when activated by adverse stimuli, including neurodegeneration (X. G. Luo et al., 2010; Perry & Holmes, 2014; Streit, Miller, Lopes, & Njie, 2008). Microglia-derived proinflammatory mediators, including TNF-α and IL-1β, have been shown to impact adult neurogenesis (Iosif et al., 2006; L’Episcopo et al., 2013; Zunszain et al., 2012); however, transcript expression for these proinflammatory cytokines did not change with age in the dorsal striatum in the current study.
On the other hand, CX3CL1, a chemokine ligand expressed in neurons that plays a vital role in neuroprotection when bound to microglia-specific CX3CR1 (Cardona et al., 2006; Paolicelli & Gross, 2011), was down-regulated with age (Table S8). CX3CR1-knockout mice exhibit decreased adult hippocampal neurogenesis, impaired LTP in the hippocampus, and impaired learning and memory (Bachstetter et al., 2011; Rogers et al., 2011), and exogenous CX3CR1 administration was able to reverse these phenotypes (Bachstetter et al., 2011)(Bachstetter et al., 2011). Decreased levels of CX3CL1 in the aging hippocampus have also been reported (Bachstetter et al., 2011; Vukovic, Colditz, Blackmore, Ruitenberg, & Bartlett, 2012). Recent studies also showed that the neuron-microglia interface via CX3CL1/R1 signaling plays a role in synaptic pruning, a process crucial in eliminating malfunctioning synapses in order to establish efficient or functional neuronal circuits (Paolicelli & Gross, 2011). CX3CR1-knockout mice exhibited reduced microglia numbers and subsequent synaptic pruning deficits, including problems with synaptic transmission and brain connectivity (Paolicelli & Gross, 2011; Zhan et al., 2014). Decreased expression of CX3CL1 with age in the dorsal striatum (Table S8) therefore suggests that neuron-microglia interaction mediated by CX3CL1/CX3CR1 signaling may be impaired during the aging process, possibly contributing to deficits in neurogenesis and synaptic plasticity.
Similar to CX3CL1, the membrane glycoprotein CD200 expressed by neurons inhibits microglia activation when bound to its receptor CD200R (Hoek et al., 2000). CD200 was down-regulated in the aging caudate (Table S2), further pointing to a potential link between increased microglial activation and aging. In addition, the microglia-specific cell surface receptor TREM2 (Colonna & Wang, 2016) was up-regulated with age in the caudate (Table 2), which can either be explained by excessive microglial proliferation or increased TREM2 expression by individual microglia (Frank et al., 2008; Raha et al., 2017). Interestingly, a rare mutant variant of TREM2 was identified as a risk factor for Alzheimer’s disease (Guerreiro et al., 2013; Jonsson et al., 2013), and TREM2 transcription was increased in the frontal cortex and hippocampus of various Alzheimer’s disease mouse models (Matarin et al., 2015; Raha et al., 2017). TREM2 is now well-known as an inductor of microglial activation and enhancer of phagocytic abilities, while also suppressing induction of proinflammatory cytokines (Neumann & Takahashi, 2007), suggesting an age-related increase of TREM2 in the caudate may influence the capacity of microglia to react to aging-related processes.
Microglia are the chief MHCII-expressing antigen-presenting cells (APCs) in the brain (X. G. Luo et al., 2010; Schetters, Gomez-Nicola, Garcia-Vallejo, & Van Kooyk, 2017). Expression of MHCII in the brain can be rapidly up-regulated on the surface of microglia and is often used as a marker of their activation (Wyss-Coray & Mucke, 2002). Previous studies have shown that, in non-human primates, microglial expression of MHCII increases with age in parallel with increases of microglial phagocytic activity, suggesting that microglial will respond more efficiently to stimulation (Peters, Josephson, & Vincent, 1991; Sheffield & Berman, 1998; Sloane, Hollander, Moss, Rosene, & Abraham, 1999). Moreover, post-mortem brain tissue transcriptomics revealed significant up-regulation of MHCII genes in the healthy aging human brain and across neurodegenerative diseases (Bossers et al., 2010; Durrenberger et al., 2015; Katsel, Tan, & Haroutunian, 2009; Parachikova et al., 2007). Our results showed significant up-regulation of MHCII antigen processing and presentation machinery in GO enrichment results in dorsal striatum from meta-analysis, including HLA-DMB, HLA-DQA2, HLA-DPA1, HLA-DPB1, HLA-DMA, and HLA-DRA (Table S10). It is known that MHCII enables antigen presentation to CD4 + T cells; however, little is known about interactions between microglia and CD4 + T cells in the aged brain. Nevertheless, age-dependent induction of MHCII transcripts in the dorsal striatum indicates functional changes in microglia, which may have detrimental effects in brain aging.
Age-dependent shortening of telomeres has been previously reported in the brain (Flanary & Streit, 2003). Given that microglia proliferate substantially throughout development and adulthood, as opposed to post-mitotic neurons, which no longer replicate once mature, telomere shortening in microglia likely significantly contributes to the overall reduction of brain telomere length. Studies have previously identified significant telomere shortening and reduction of telomerase activity in microglia in aged rats (Flanary, Sammons, Nguyen, Walker, & Streit, 2007; Flanary & Streit, 2004). Telomere shorting in humans has been correlated with cognitive decline and dementia (Eitan, Hutchison, & Mattson, 2014; Grodstein et al., 2008; Honig, Kang, Schupf, Lee, & Mayeux, 2012; Honig, Schupf, Lee, Tang, & Mayeux, 2006; Valdes et al., 2010). Age-associated down-regulation of the telomere binding protein RAP1 signaling pathway in the dorsal striatum identified in the present study (Fig. 3C; Table S7, S9) suggests that telomere dysfunction in microglia could influence age-related changes in brain function.
Recent studies suggest that immune-system related molecules, including complement and MHC class I (MHCI) proteins, which play a role in synapse formation and pruning in the young brain, may erroneously eliminate functional synapses in late adulthood, contributing to neurodegenerative diseases (Lazarczyk et al., 2016; Y. Wu, Dissing-Olesen, MacVicar, & Stevens, 2015). C3, an essential player in both classical and alternative complement cascades, localizes to immature synapses and recruits complement receptor 3 (CR3)-expressing microglia for synapse elimination (Cribbs et al., 2012; Stephan, Barres, & Stevens, 2012; Stevens et al., 2007). Microglia-mediated synapse elimination is decreased in C3 and CR3-deficient mice (Schafer et al., 2012), suggesting that complement proteins located at immature synapses result in the targeted removal of these synapses by microglia. Interestingly, CR3 expression levels are increased with age in humans (Cribbs et al., 2012) and play a role in the erroneous elimination of functional synapses, thereby contributing to age-associated cognitive decline (Shi et al., 2015). C1q, another player in the complement cascade (Stephan et al., 2012; Stevens et al., 2007), is also overexpressed in aged human and mouse brains but does not appear to mediate synapse elimination (Stephan et al., 2013). Our results revealed significant up-regulation of C1QA, C1QB, CIQC, and C3AR1 in the caudate with age (Table S2), and meta-analysis showed an age-dependent increase of C1QB and C3AR1 expression in the dorsal striatum (Table S8)
Like complement cascade proteins, MHCI proteins have been implicated in synapse maintenance and synaptic plasticity in the developing CNS throughout non-pathological aging (Goddard, Butts, & Shatz, 2007; Huh et al., 2000; Lazarczyk et al., 2016; H.Lee et al., 2014). MHCI protein expression was significantly increased in the cognitively normal aging human brain (Katsel et al., 2009), indicating that increase in MHCI protein levels might be important in regulating and preserving cognitive function. The age-associated up-regulation of antigen processing and presentation pathway in both caudate and putamen (Fig. 3C), together with up-regulation of antigen processing and presentation machinery in GO enrichment results in dorsal striatum from meta-analysis, including MHCI molecules HLA-B and HLA-F (Table S10), suggest that enhanced MHCI gene expression in the dorsal striatum might contribute to synaptic formation and plasticity during normal aging. However, whether up-regulation of MHCI genes in healthy aging subjects is functionally associated with antigen presentation or is involved in interactions between neural and glial cells needs further investigation.
In contrast to the complement and coagulation cascades pathways that were up-regulated, glutamatergic, GABAergic, and cholinergic synapse proteins were down-regulated in both caudate and putamen in aging striatum (Fig. 3C). Age-associated down-regulation of dopaminergic synaptic pathways and long-term potentiation and long-term depression (LTP/LTD) pathways were identified as two of the most significantly enriched pathways in the caudate (Fig. 3B).
Canonical WNT signaling in age-mediated neuroinflammation and plasticity
The WNT pathway recently has been linked with the regulation of multiple processes in the adult brain, including neurogenesis, axonal growth, synapse formation and plasticity, and cognitive health (Inestrosa & Arenas, 2010; Maiese, Li, Chong, & Shang, 2008; Oliva, Vargas, & Inestrosa, 2013). In addition, WNT molecules act as anti-inflammatory factors, commonly secreted by astrocytes to suppress pro-inflammatory cytokines released by neighboring microglial and prevent their activation (L’Episcopo et al., 2014; Marchetti & Pluchino, 2013). The extracellular WNT molecules are compartmentalized in the cell into three separate major pathways: the canonical WNT/β-catenin pathway and the non-canonical WNT/planar cell polarity (PCP) or WNT/Ca2+ pathways (Inestrosa & Varela-Nallar, 2014). Our results showed significant down-regulation of WNT7A, WNT7B, and WNT8A (Table S10), which are activators of the canonical WNT/β-catenin pathway. Moreover, the PI3K-Akt signaling pathway, a positive regulator of the canonical WNT pathway, was down-regulated with age (Table S6, S7, S9) (Fukumoto et al., 2001), suggesting age-mediated down-regulation of canonical WNT signaling. Many studies have shown that canonical WNT/β-catenin pathway participates in adult neurogenesis, axon remodeling, and synaptic plasticity (Lie et al., 2005; Oliva et al., 2013). In addition, β-catenin down-regulation was shown to be involved in inflammation-induced neurogenic impairment (L’Episcopo et al., 2013). Therefore, age-associated inhibition of canonical WNT/β-catenin signaling may both exacerbate a deterioration in the neuroprotective milieu and promote inflammatory cascades.
The blood-brain barrier (BBB) is a semi-permeable membrane made up of a specialized endothelial lining whose integrity is sustained by microglia, astrocytes, pericytes, and neurons. In a healthy state, the BBB restricts access of immune cells into the CNS and maintains homeostasis. During inflammation, BBB integrity is compromised and immune cells can enter the CNS parenchyma. An impaired BBB is a hallmark of the aging brain (Erdo, Denes, & de Lange, 2017). Recent studies using MRI in healthy subjects with no or mild cognitive impairment found that BBB permeability and overall loss of integrity in the hippocampus increased linearly with age, and was highest in subjects with mild cognitive impairment. This suggests that loss of BBB integrity is an early event in the aging human hippocampus that may contribute to cognitive impairment (Montagne et al., 2015). Notably, the death receptor 6 (TNFRSF21) and TROY (TNFRSF19) were shown to act downstream of WNT/β-catenin signaling in BBB endothelial cells to regulate CNS angiogenesis and barrier genesis (Tam et al., 2012). Dysregulation of TNFRSF21/TNFRSF19 signaling was shown to disrupt the BBB (Sonar & Lal, 2015). Our results identified down-regulation of both these TNFRSF members with age (Table S8). Thus down-regulation of the WNT/β-catenin signaling pathway may influence TNFRSF21 and TNFRSF19 expression and contribute to structural and functional impairment in the BBB during aging.
Conclusion
Down-regulation of transcripts involved in neurogenesis, synaptogenesis, synapse pruning, and canonical WNTs, with contrasting up-regulation of neuroinflammatory transcripts was identified in the aging caudate and putamen. These data suggest that age-dependent attenuation of canonical WNT signaling in the dorsal striatum may play a pivotal role at the neuron-glial interface to regulate striatal plasticity during aging. WNT/β-catenin signaling dysregulation via astrocyte–microglial interactions during aging sheds light on novel mechanisms driving the decline of striatal neurogenesis and altered synaptic connectivity and plasticity, ultimately leading to impaired motor and cognitive abilities. These findings suggest promising therapeutic opportunities for targeting WNT/β-catenin signaling in order to enhance endogenous restoration or to promote preservation of brain function during normal aging.
Author contributions
C.W. conducted data analysis, C.W., R.M.B., W.J.F., and C-T.L. wrote the manuscript and contributed to the interpretation of data. C.W. and C-T.L. are the guarantors of this work and take responsibility for the accuracy of the data analysis.
