Original Research ARTICLE
Transcription Profile of Aging and Cognition-Related Genes in the Medial Prefrontal Cortex
- 1Department of Neuroscience, McKnight Brain Institute, University of Florida, Gainesville, FL, USA
- 2Genetics and Genomics Program, Genetics Institute, University of Florida, Gainesville, FL, USA
Cognitive function depends on transcription; however, there is little information linking altered gene expression to impaired prefrontal cortex function during aging. Young and aged F344 rats were characterized on attentional set shift and spatial memory tasks. Transcriptional differences associated with age and cognition were examined using RNA sequencing to construct transcriptomic profiles for the medial prefrontal cortex (mPFC), white matter, and region CA1 of the hippocampus. The results indicate regional differences in vulnerability to aging. Age-related gene expression in the mPFC was similar to, though less robust than, changes in the dorsolateral PFC of aging humans suggesting that aging processes may be similar. Importantly, the pattern of transcription associated with aging did not predict cognitive decline. Rather, increased mPFC expression of genes involved in regulation of transcription, including transcription factors that regulate the strength of excitatory and inhibitory inputs, and neural activity-related immediate-early genes was observed in aged animals that exhibit delayed set shift behavior. The specificity of impairment on a mPFC-dependent task, associated with a particular mPFC transcriptional profile indicates that impaired executive function involves altered transcriptional regulation and neural activity/plasticity processes that are distinct from that described for impaired hippocampal function.
The extent to which cognition declines over the course of aging varies across individuals. Neuroimaging studies indicate that the pattern of cognitive decline is related to changes in the structure and activity of the prefrontal cortex (PFC) and hippocampus (Grady et al., 2005; Persson et al., 2006; Dennis et al., 2008; Park and Reuter-Lorenz, 2009; Migo et al., 2016), suggesting that individual differences in cognitive aging may result from vulnerability and reorganization of these neural systems. In addition, altered white matter integrity could influence connectivity of the PFC with other brain regions (O'Sullivan et al., 2001; Pfefferbaum et al., 2005; Salat et al., 2005; Andrews-Hanna et al., 2007; Bennett et al., 2011; Borghesani et al., 2013). The molecular mechanisms for vulnerability and adaptive reorganization during aging have been the subject of speculation (Jackson et al., 2009; Kumar et al., 2009; McEwen and Morrison, 2013; Gray and Barnes, 2015). Previous work using microarray technology indicates that over the course of aging, the transcription of genes linked to inflammation and synaptic function increases and decreases, respectively, within a number of brain regions (Prolla, 2002; Blalock et al., 2003; Verbitsky et al., 2004; Erraji-Benchekroun et al., 2005; Loerch et al., 2008; Burger, 2010; Bordner et al., 2011; Haberman et al., 2011; VanGuilder et al., 2011; Zeier et al., 2011; Cribbs et al., 2012; Yuan et al., 2012; Berchtold et al., 2013; Primiani et al., 2014), suggesting possible mechanisms for variability in cognitive decline.
While it is widely thought that transcription is linked to cognitive function, there is relatively little information on the PFC transcriptional profile, which attempts to link altered gene expression to an age-related decline in behaviors that depend on the PFC. Indeed, the PFC provides several unique challenges for examining the relationship of transcription to age-related cognitive impairment. The PFC can be divided into several sub-regions and there is a long-standing debate over the equivalence of anatomical regions within the PFC across species (Uylings and van Eden, 1990; Preuss, 1995; Brown and Bowman, 2002; Vertes, 2004; Hoover and Vertes, 2007). In addition, the PFC is involved in executive function, which encompasses a number of cognitive processes including attention, response inhibition, working memory, and mental flexibility (Robbins, 1996; Bizon et al., 2012).
In the current study, we exploit individual differences in behavior to examine the relationship between age-related changes in cognition and transcription. Young and aged rats were characterized on two tasks that are age-sensitive, including an attentional set shift task that depends on the mPFC (Brown and Bowman, 2002; Kesner and Churchwell, 2011) and on a hippocampal-dependent spatial episodic memory task (Foster, 2012; Foster et al., 2012). RNA sequencing (RNA-seq) was used to construct transcriptomic profiles for the mPFC, white matter, and CA1 region of the hippocampus. Expression differences associated with aging and cognition, defined by variability in set shift or spatial memory behavior, were examined. Finally, the aging and cognition mPFC gene sets were compared to microarray data from other studies to test specific hypotheses. The results indicate that expression of immediate-early genes (IEGs) related to neural activity and synaptic plasticity decline with age in the mPFC; however, within the group of aged animals, expression of IEGs is up regulated in animals that exhibit delayed set shift behavior.
Procedures involving animal subjects have been reviewed and approved by the Institutional Animal Care and Use Committee and were in accordance with guidelines established by the U.S. Public Health Service Policy on Humane Care and Use of Laboratory. Male Fischer 344 rats of two ages, young (5–6 months, n = 11) and aged (17–22 months, n = 20) were obtained from National Institute on Aging colony (Taconic) through the University of Florida Animal Care and Service facility. Animals were maintained on a 12:12 h light schedule, and provided ad lib access to food and water prior to the set shifting task.
Set Shifting Operant Task
Testing in the set shifting task was conducted in standard rat behavioral test chambers (30.5 × 25.4 × 30.5 cm, Coulbourn Instruments, Whitehall, PA) with metal front and back walls, transparent Plexiglas side walls, and a floor composed of steel rods (0.4 cm in diameter) spaced 1.1 cm apart. Each test chamber was housed in a sound-attenuating cubicle, and was equipped with a recessed food pellet delivery trough located 2 cm above the floor in the center of the front wall. The trough was fitted with a photobeam to detect head entries and a 1.12 W lamp for illumination. Food rewards consisted of one 45 mg grain-based food pellet for each correct response (PJAI, Test Diet, Richmond, IN). Two retractable levers were located to the left and right of the food trough (11 cm above the floor), and a 1.12 W cue lamp was located 3.8 cm above each lever. An additional 1.12 W house light was mounted near the top of the rear wall of the sound-attenuating cubicle. An activity monitor was positioned above each test chamber to monitor locomotor activity throughout each session. This monitor consisted of an array of infrared (body heat) detectors focused over the entire test chamber. Movement in the test chamber (in x, y, or z planes) was defined as a relative change in the infrared energy falling on the different detectors. A computer interfaced with the behavioral test chambers and equipped with Graphic State 3.01 software (Coulbourn Instruments) was used to control experiments and collect data.
The design of the set shifting task was based previously published methods (Floresco et al., 2008; Beas et al., 2013). Prior to the start of behavioral testing, rats were reduced to 85% of their free feeding weights over the course of 5 days and maintained at this weight for the duration of the experiments. Rats were trained and tested in the same behavioral testing chamber during the course of the experiment. Rats progressed through four stages of shaping prior to the start of the set shifting task, with new stages beginning on the day immediately following completion of the previous stage. On the day prior to Shaping Stage 1, each rat was given five 45 mg food pellets in its home cage to reduce neophobia to the food reward used in the task. Shaping Stage 1 consisted of a 64-min session of magazine training, involving 38 deliveries of a single food pellet with an inter-trial interval (ITI) of 100 ± 40 s. Shaping Stage 2 consisted of lever press training, in which a single lever (left or right, counterbalanced across groups) was extended and a press resulted in delivery of a single food pellet. After reaching a criterion of 50 lever presses in 30 min, rats were then trained on the opposite lever using the same procedures.
Shaping Stage 3 consisted of 90 trials that were designed to train rats to press the levers after their insertion into the test chamber. Each 20 s trial began with illumination of the house light and insertion of a single lever (either left or right, randomly selected within each pair of trials) into the test chamber where it remained for a maximum of 10 s. A response on the lever within this time window resulted in retraction of the lever, delivery of a single food pellet, and continued illumination of the house light for an additional 4 s. If a rat failed to respond on the lever within 10 s, the lever was retracted and the house light turned off, and the trial was scored as an omission. Rats received at least 4 daily sessions in this stage, and were trained until reaching criterion performance of fewer than 10 omissions out of the 90 trials.
Shaping Stage 4 was designed to determine each rat's side bias (i.e., preference for one lever over the other). Each trial consisted of multiple phases. In the first phase of a trial, the house light was illuminated and both levers were inserted into the test chamber. A response on either lever resulted in retraction of both levers and delivery of a single food pellet. In the second phase of a trial, both levers were again inserted, but only a response on the lever opposite to that chosen in the first phase resulted in food delivery. A response on the same lever chosen in the first phase (i.e., “incorrect”) resulted in the levers being retracted and the house light being extinguished. After a “correct” response in this second phase of a trial, a new trial was initiated, whereas after an “incorrect” response, the second phase of the trial was repeated. The second phase was repeated until rats made a “correct” response. The session ended after a total of 45 completed trials. The side associated with the greatest number of total responses across this phase of testing was considered a rat's biased side.
Visual cue discrimination
Following shaping stage 4, rats were trained to press the lever signaled by the illumination of a cue light over the lever. Each 20 s trial began with illumination of one of the cue lights (left or right, randomly selected in each pair of trials). After 3 s, the house light was illuminated and both levers were inserted into the chamber (the cue light remained illuminated while the levers were extended). A response on the lever corresponding to the cue light (a correct response) resulted in the house light remaining on for 4 s, during which time the levers were retracted, the cue light was extinguished, and a single food pellet was delivered. A response on the opposite lever (an incorrect response) or failure to respond within 10 s (omission) resulted in retraction of both levers and all lights being extinguished. Rats were considered to have acquired the task upon reaching criterion performance of eight consecutive correct trials (and at least 30 total trials, excluding omissions), with the maximum number of trials per session set at 120. Rats that failed to acquire the task within a single session (young: n = 5; aged: n = 13) received additional sessions on subsequent days.
Left/Right discrimination (Set shift)
After reaching criterion performance on the visual cue discrimination, rats were tested the next day in the set shift condition, in which the task contingencies were altered. In this condition, rats were required to ignore the visual cue and instead to consistently choose the left or right lever (whichever was not their biased side as determined in Shaping Stage 4). Hence, accurate performance required rats to “shift” their attention away from the visual cue and toward the left/right position of the lever. Beyond the shift in reward contingencies, trials were identical in presentation to those in the visual cue discrimination (i.e., on each trial, both levers were presented; with the cue light illuminated over one lever). As in the visual cue discrimination, the location of the illuminated cue light was randomized (left or right) in each pair of trials. Rats were considered to have acquired the task upon reaching criterion performance of eight consecutive correct trials, excluding omissions. The maximum number of trials per session was set at 120 and all rats acquired the task within a single session.
Morris Water Maze
Following completion of set shifting, animals were again provided ad lib access to food and water for ~7 weeks prior to testing on the water maze. Animals were trained in a black tank, 1.7 m in diameter, positioned in a well-lit room. The pool was surrounded by black walls and black curtain. For spatial training an assortment of two- and three-dimensional cues were hung on the walls and curtain. Water (27 ± 2°C) was maintained at a level ~8 cm below the surface of the tank. Methods employed to assess sensory-motor deficits and impaired episodic spatial memory on the water maze have been published previously (Foster et al., 1991; Kumar and Foster, 2013; Guidi et al., 2014). For cue and spatial tasks, training consisted of five blocks with three trials per block and training on each task was massed into a single day. Inter-trial intervals were 20 s and inter-block intervals were ~15 min. Rats remained on the platform between trials and in home cages under the heat lamp after each block. Behavioral data was acquired with Noldus EthoVision computer tracking software (Noldus Information Technology, Leesburg, VA, USA) and included path-length and latency to escape to the platform, platform crossing and time in the goal and opposite quadrants.
Rats were first trained on the cue discrimination version of the water escape task. The escape platform was extended ~1 cm above the water level and a white Styrofoam flag was attached. For each trial, the platform position and start location were randomized. If an animal did not escape the water maze within 60 s, the rat was gently guided to the platform. Three days following cue training, animals were trained on the spatial discrimination task. For spatial discrimination, the escape platform was hidden ~1.5 cm beneath the water level and remained in the same location relative to the distal cues in the room for the duration of the initial spatial training. Fifteen minutes following the end of training on block 5, a free-swim probe trial was administered as a measure of learning. For the probe trial, the platform was removed and the animal placed in the tank for 60 s. A spatial discrimination index was computed according to the formula (G − O)/(G + O) where G and O represent the percent of time spent in the goal quadrant and quadrant opposite the goal, respectively.
Statistical analysis of behavior
The total numbers of trials required to achieve criterion (TTC) on the visual cue discrimination and on the left/right discrimination (set shift) were used as the indices of performance. Mean distance to find the platform during each training block for the water maze cue and spatial tasks and probe trial data, platform crossing, and discrimination index, were employed to examine learning on the water maze. For the distance measures, repeated measures analyses of variance (ANOVAs) were used to examine age and training effects. One way ANOVAs were used to examine aged effects for the water maze probe trial data and TTC measures from the operant tasks. Fisher's protected least significant difference comparisons, with the p-value set at 0.05, were used to localize differences.
Two weeks following water maze testing, rats were anesthetized with isoflurane (Piramal Healthcare), decapitated and the brain was rapidly removed. The PFC was blocked into 1 mm coronal slices. The mPFC including the prelimbic and infralimbic regions were collected from two sections (between +5.0 and +2.5 anterior to bregma; Paxinos and Watson, 1986). For a subset of animals (young = 8, aged = 9), white matter was collected adjacent to the mPFC (Figure 1). For region CA1, the hippocampus was isolated, a 1–2 mm slice was removed from the dorsal hippocampus, and the CA1 region was dissected (Blalock et al., 2003; Zeier et al., 2011). The collected tissue was immediately frozen in liquid nitrogen and stored in −80°C until processed.
Figure 1. Region of the mPFC and white matter (WM) collected for RNA-seq. The right panel provides a schematic of a coronal slice +2.7 anterior to bregma diagram as adapted from Paxinos and Watson (1986) and illustrates the region of the mPFC and white matter collected for RNA-seq. The left panel shows a coronal slice from this same region.
RNA and Library Preparation
RNA was isolated from the mPFC, CA1, and white matter using the RNeasy Lipid Tissue Mini kit (Qiagen, catalog number 74804) and DNase digestion was performed with the RNase-Free DNase Set (Qiagen, catalog number 79254). The concentration was measured with the NanoDrop 2000 spectrophotometer and the RNA integrity number (RIN) was quantified by the University of Florida Interdisciplinary Center for Biotechnology Research using the High Sensitivity RNA Screen Tape in an Agilent 2200 Tapestation system. The average RIN across all regions was 8.02 (±SEM 0.05). External RNA Controls Consortium (ERCC) spike-in controls (Thermo Fisher, catalog number 4456740) were added to a subset of samples and the mRNA was selected with the Dynabeads mRNA DIRECT Micro kit (Thermo Fisher, catalog number 61021). Whole transcriptome libraries were prepared with the Ion Total RNA-seq Kit v2 (Thermo Fisher, catalog number 4475936) with the addition of the Ion Xpress barcodes for multiplex sequencing (Thermo Fisher, catalog number 4475485). The concentration of the libraries was quantified by the Qubit dsDNA HS Assay (Thermo Fisher, catalog number Q32851) and size distribution was evaluated with the High Sensitivity D1000 Screen Tape in the Tapestation system.
Reverse Transcription Quantitative Polymerase Chain Reaction
Reverse transcription quantitative polymerase chain reaction (RT-qPCR) was performed in a subset of samples to validate RNA-seq results. cDNA was prepared using the QuantiTect Reverse Transcription kit (Qiagen, catalog number 205311) and quantitative PCR was completed with the TaqMan Gene Expression Assays (Arc: Rn00571208_g1, Egr1: Rn00561138_m1, Egr2: Rn00586224_m1, Egr4: Rn00569509_g1, Fos: Rn02396759_m1, Lin7b: Rn00572781_m1, Gapdh: Rn01775763_g1) in a 7300 Real-Time PCR system with SDS software version 1.3.1 (Applied Biosystems). The ΔΔCT method (Livak and Schmittgen, 2001) was used to determine the relative cDNA levels. Differences in the subset of RNA-seq and RT-qPCR were confirmed using t-tests between young and aged rats and between age impaired and unimpaired animals.
Sequencing, Bioinformatics, and Statistical Analysis
Template preparation was performed in the Ion Chef system and sequencing was completed in the Ion Proton (Thermo Fisher). ERCC analysis was executed in the Torrent Server with the ERCC analysis plugin. Spiked samples contained R2 above 0.9 with at least 60 transcripts. On average, each sample contained 18.8 million reads of 131 base pair length. Low quality reads were removed from the FASTQ files and the data was aligned to the rn5 genome using the two step alignment method for Ion Proton transcriptome data with TopHat2 and Bowtie2 in the Partek Flow servers (Partek Inc.). Gene-level counts were generated from BAM files using the featureCounts function in the R package Rsubread (Liao et al., 2014). The Rnor_5.0.78.gtf file was used for annotation and count normalization was performed with the DESeq package in R. The data for this study has been uploaded to NCBI's Gene Expression Omnibus under the accession number: GSE75772.
Gene filtering and initial statistical analysis was performed according to our previously published work (Blalock et al., 2003; Aenlle et al., 2009; Aenlle and Foster, 2010; Zeier et al., 2011). Gene lists were initially filtered to remove those genes with counts of 5 or less, which resulted in the detection of over 18,000 Ensembl database genes from each area. Gene lists were further filtered such that only genes with annotation of at least one gene ontology (GO) term were considered for gene enrichment analysis. Filtering for GO terms resulted in 15075 mPFC genes, 15235 genes, and 15084 white matter genes. For differential expression analysis associated with age, a statistical filter was performed in each tissue type independently, using a one-way ANOVA generated in Partek Genomics Suite 6.6, with p < 0.025 according to our previous work (Blalock et al., 2003; Aenlle et al., 2009; Aenlle and Foster, 2010; Zeier et al., 2011). For examination of mPFC gene expression related to cognition, Pearson's correlations were calculated between the behavioral TTC measures for the set shifting task or discrimination index score for spatial learning and the expression of each gene in the mPFC transcriptome. Correlations were limited to aged animals in order to remove age as a confound. Due to multiple comparisons, the confidence in any single gene is low; therefore, gene enrichment analysis was performed under the assumption that changes in biological process with age or cognition would result in a shift in the expression of clusters of genes related to the biological process. For gene enrichment and functional annotation clustering analysis, data sets of genes that exhibited an increase or decrease in expression were separately submitted to the NIH database for annotation, visualization, and integrated discovery (DAVID; Huang et al., 2007a,b). Enrichment analysis was limited to gene ontology for biological processes and cellular components with the Benjamini False Discovery Rate (FDR) p < 0.05 as a cut-off for cluster selection. The heat map were generated in Partek Genomics Suite 6.6 using the genes that were identified in DAVID with counts which were standardized to z-scores. In other cases, specific hypotheses were tested by comparing gene expression with previous published work using microarrays. In this case, we determined whether the previously published genes were detected by our procedures. This set of genes that were common across studies represented the total data set. Next, we used a fold change and chi squared test to determine if the genes were altered in the same direction, relative to what would be expected by chance. The number of genes that were significantly altered in the same direction was determined using one-tailed t-tests with the direction specified by the microarray studies. A FDR was calculated by calculating the number of genes expected to change in a specified direction using the formula T*p, where T = total number of genes tested and p is the significance level (0.05). The expected number of genes was divided by the number of genes that were significantly different in the predicted direction to obtain the FDR.
Young (n = 11) and aged (n = 20) animals were trained on the visual discrimination operant task followed by set shift testing. Aged animals exhibited considerable variability in performance on each task (Figure 2), with some aged animals performing in a range similar to young. Examination of the TTC for visual discrimination indicated no effect of age [F(1, 29) = 2.15, p = 0.15; Figure 2A]. In contrast, an age difference in TTC was observed for the set shifting (left/right discrimination) behavior [F(1, 29) = 7.95, p < 0.01], with a subset of aged rats exhibiting an increase in TTC relative to young (Figure 2B).
Figure 2. Performance on the visual discrimination and set shift operant tasks. Trials to criteria (TTC) are illustrated for individual aged (filled circles, n = 20) and young (open circles, n = 11) animals during performance of the (A) initial visual discrimination and (B) set shift tasks. Asterisk indicates that aged animals exhibited more trails to criteria for the set shift task (p < 0.01). The open bars indicate the mean TTC for each group.
For the cue discrimination version of the water maze task, all animals were able to find the visible platform during the 60 s time limit during the last three trials (block 5). A repeated measures ANOVA for the cue discrimination task indicated an effect of training [F(4, 116) = 7.84, p < 0.0001] and age [F(1, 29) = 11.43, p < 0.005] and an interaction of age and training [F(4, 116) = 2.92, p < 0.05] due to superior performance by young animals during the final training blocks (Figure 3A). A repeated measures ANOVA for the spatial discrimination task indicated an effect of training [F(4, 116) = 17.82, p < 0.0001] and an interaction of age and training [F(4, 116) = 3.13, p < 0.05; Figure 3B]. The results of the probe trial indicated a decrease in platform crossings for aged animals [F(1, 29) = 8.52, p < 0.01; Figure 3C]. No age effect was observed for the discrimination index; however, consistent with previous reports (Blalock et al., 2003; Foster, 2012; Kumar and Foster, 2013; Guidi et al., 2014), there was considerable variability (Figure 3D). For aged animals, the relationship between behavioral measures on the operant tasks (visual discrimination TTC, set shift TTC) and on the water maze (distance to escape on block 5 of the cue task, distance to escape on block 5 of the spatial task, platform crossings, and discrimination index of the probe trial) was examined (Table 1). The results indicated a correlation between the distance on block 5 of the spatial task and the discrimination index (r = −0.484, p < 0.05), such that longer escape distances on the spatial task were associated with poorer discrimination index scores.
Figure 3. Performance on the water maze task. Symbols indicate the mean (±SEM) escape path length to the escape platform during five training blocks on the (A) cue and (B) spatial discrimination tasks for young (open symbols) and aged (filled symbols) animals. Individual (C) platform crossing and (D) discrimination index scores for young (open symbols) and aged (filled symbols) animals. The open bars indicate the means for each group.
Gene Expression Related to Aging
RNA-seq libraries from the mPFC and region CA1 of the hippocampus were prepared from all animals (young = 11, aged = 20). In addition, libraries were prepared from white matter, collected from a subset of animals (young = 8, aged = 9). Gene expression data was statistically filtered for age differences using ANOVAs for each tissue type with a cut-off set at p < 0.025. The largest and smallest number of age-related genes was observed for white matter (1529 genes) and region CA1 (286 genes), with 731 genes altered in the mPFC. Although, expression changes were relatively distinct, there was some overlap. Furthermore, for genes that overlapped across any two regions, the number of up regulated genes was 2–8 times larger than genes that decreased expression (Figure 4).
Figure 4. Number of genes altered during aging across regions. Graphic summary of the total number of genes whose expression either significantly (p < 0.025) increased (up arrow) or decreased (down arrow) in the mPFC, region CA1, white matter, and across regions.
Based on functional studies, it has been suggested that the mPFC of the rat may be analogous to the dorsolateral PFC of humans and non-human primates (Dias et al., 1997; Birrell and Brown, 2000; Kesner and Churchwell, 2011). To determine possible expression changes in the rat mPFC that may correspond to changes observed in humans, we compared our results with a microarray study examining age-related changes in the dorsolateral (Brodmann's area 9) PFC of humans (Erraji-Benchekroun et al., 2005). For the 414 age-related probes from the dorsolateral PFC of the original human data set, we were able to identify 318 rat genes in the mPFC and the direction of change, increasing or decreasing, was determined for the rat. From the set of 318 genes, we observed 203 genes (64%) that changed in the same direction as predicted in humans and a chi square analysis indicated that the directional changes were different than that expected by chance (x2 = 11.86, p < 0.001). One-tailed t-tests (p < 0.05 with direction specified by age-related changes in humans) were conducted to examine age differences in expression of these 318 genes. A total of 60 genes (39 decreased, 21 increased) reached significance (FDR: 0.27; Table 2).
The data sets for differentially expressed up regulated and down regulated genes were separately submitted to NIH DAVID for enrichment analysis based on gene ontology for biological processes and cellular components. Cluster selection cut-off was restricted to clusters with a Benjamini FDR p < 0.05. In each case, several genes were observed in multiple related clusters, which are shown below.
In general, age-related changes in transcription involve up regulation of genes linked to immune response, oxidative stress, and the lysosome (Lee et al., 2000; Blalock et al., 2003; Verbitsky et al., 2004; Fraser et al., 2005; Rowe et al., 2007; de Magalhães et al., 2009; Kadish et al., 2009; Lipinski et al., 2010; VanGuilder et al., 2011; Zeier et al., 2011; Yuan et al., 2012). Expression of immune response related genes was particularly evident in the white matter and mPFC. For the 940 white matter genes, enrichment was observed for immune response (GO:0006955, 85 genes, FDR p = 4.1−25), defense response (GO:0006952, 70 genes, FDR p = 1.3−16), and the lysosome (GO:0005764, 33 genes, FDR p = 2.2−9). Similarly, for the 344 genes that increased expression in the mPFC, gene enrichment was observed for biological processes linked to oxidation reduction (GO:0055114, 30 genes, FDR p = 0.02), response to wounding (GO:0009611, 24 genes, FDR p = 0.02), and adaptive immune response (GO:0002250, 10 genes, FDR p = 0.01; Table 3). Increased expression in the mPFC was also observed for carboxylic and catabolic process (GO:0046395, 13 genes, FDR p = 0.002; Figure 5). Finally, for 182 CA1 genes, enrichment was observed for the lysosome (GO:0005764, 10 genes, FDR p = 0.007). Up regulation was observed for response to wounding (GO:0009611, 11 genes, FDR p = 0.38) and defense response (GO:0006952, 10 genes, FDR p = 0.44); however, the clusters did not reach the FDR cut-off.
Figure 5. Heat map of age-related changes in gene expression for the mPFC. Each row represents a differentially expressed gene (p < 0.025) associated with aging. Expression for each gene was converted to a standardized score and the color represents the standard deviation increasing (red) or decreasing (blue) relative to the mean (gray). The age-related gene enrichment clusters are indicated (FDR p < 0.05). Top clusters: Genes that exhibit a decrease from young (left) to aged (right) animals. Bottom: Genes that exhibit an increase from young to aged animals.
Previous work indicates that aging is associated with decreased expression of genes linked to neuronal/synaptic genes (Blalock et al., 2003; Lu et al., 2004; Verbitsky et al., 2004; Burger et al., 2007; Aenlle and Foster, 2010; Zeier et al., 2011; Berchtold et al., 2013; Primiani et al., 2014). In the case of the 589 genes that were decreased in white matter, enrichment was observed for cell division (GO:0051301, 16 genes, FDR p = 0.04). For the 104 genes that decreased in area CA1, clustering did not pass our cut-off. For the 387 mPFC genes that decreased with age, common genes were observed for clusters related to the synapse (GO:0045202, 25 genes, FDR p = 1.2−4) and postsynaptic membrane (GO:0045211, 14 genes, FDR p = 4.1−4; Table 4). In addition, decreased expression was observed for neuron projection (GO:0043005, 22 genes, FDR p = 0.015; Figure 5). Thus, the age-related decrease in neuronal genes was particularly evident in the mPFC.
mPFC Gene Expression Related to Behavior
To examine behavioral specificity of transcriptional changes, Pearson's correlations were run comparing expression of mPFC genes with the TTC for set shifting and visual discrimination, and the discrimination index score for the water maze. Correlations were limited to aged animals in order to remove age as a confound and correlations were performed across all genes in the mPFC transcriptome. Using a cut-off set at p < 0.025 (r = 0.499), a total of 416 genes were correlated with behavioral flexibility. Most of the genes (73%) were positively correlated with the set shift TTC score (303 genes increasing expression with impairment), with 113 genes negatively correlated (decreased in impaired animals). A similar analysis for the visual discrimination TTC indicated many fewer mPFC genes (135 total) correlated with acquisition of the visual discrimination, with 103 mPFC genes positively correlated (increased in animals with more trials to criteria) and 32 genes that were negatively correlated (decreased in animals with more trials to criteria). Unlike the set shift behavior, more mPFC genes exhibited decreased expression (62%) in animals with poorer spatial learning. A total of 401 mPFC genes correlated with the water maze discrimination index, with 249 genes positively correlated with the discrimination index (decreased in animals with poor spatial learning) and 152 mPFC genes negatively correlated (increased in spatial learning impaired animals; Figure 6).
Figure 6. Number of genes correlated with behavioral measures. Graphic summary for the number of mPFC genes whose expression either significantly increased (filled) or decreased (open) in relation to performance on the set shift and visual discrimination operant tasks, and a spatial memory task (Pearson's correlation p < 0.025). A large number of genes exhibited increased expression associated with impaired set shifting compared to visual discrimination and spatial discrimination index (DI).
Gene enrichment analysis was conducted to determine which mPFC genes and biological processes might be good markers for the age-related impairment in set shift behavior. For mPFC genes that were negatively correlated with the set shift TTC (decreased expression in impaired animals that delay shifting), no gene enrichment clusters were observed to pass the cut-off for cluster selection. For gene expression that positively correlated with the TTC score (increasing in animals that delayed shifting), clusters were observed for response to organic substance (GO:0010033, 34 genes, FDR p = 0.01), regulation of apoptosis (GO:0042981, 27 genes, FDR p = 0.02) and the regulation of transcription (GO:0045449, 46 genes, FDR p = 0.027), which contained the largest number of genes (Table 5). An index score describing expression of transcription regulation genes was generated for each animal by first standardizing the expression of the 46 transcription genes and the standard scores were averaged within each animal. The standard scores were plotted against the standardized TTC set shift score (r = 0.89) to illustrate the correspondence of transcription regulators with set shift behavior (Figure 7). Interestingly, several IEGs linked to neuronal activity (Arc, Egr1, Egr2, Egr3, Egr4, Fos, Fosb, Fosl2, Junb) were observed to increase in association with delayed set shift behavior.
Figure 7. Impaired set shifting is associated with increased expression of genes involved in transcription regulation. The z-scores for the cluster of 46 transcription regulation genes were averaged for each animal. Individual mean z-scores (y-axes) are plotted relative to z-scores for set shift TTC (x-axes). The correlation is illustrated as a regression line (dashed line).
Next we compared our results to a previously published microarray study that examined gene expression across rat species that differ on tests of attention (Qiu et al., 2010). The spontaneously hypertensive-rat (SHR), a well characterized model of impaired attention, exhibiting impairment on the five-choice serial reaction time task and attentional set shift relative to the Wistar-Kyoto rat (De Bruin et al., 2003; Sagvolden et al., 2005). Using RNA-seq, we were able to detect 48 genes that were correlated with set shift behavior and were previously reported to differ in the mPFC of SHR and Wistar-Kyoto rats (Qiu et al., 2010). In order to compare our data with the results of Qiu et al. (2010), a mean split for the set shift TTC scores of aged animals (mean TTC = 51.7) was used to separate aged animals into those that delayed shifting and were considered aged-impaired (AI) and aged-unimpaired (AU). One-tailed t-tests were run on the 48 genes using a p < 0.05 and the direction specified by the results of Qiu et al. (2010). For the 26 genes that were predicted to increase, 16 genes exhibited a significant increase (FDR: 0.08; Table 6). Interestingly, for the 16 genes that exhibited a significant increase in expression in AI rats, these genes also exhibited decreased expression during aging with seven genes (Arc, Egr1, Egr3, Egr4, Junb, Klf10, Nr4a3) exhibiting a significant (p < 0.05) decrease with age. Finally, for the 22 genes that were previously reported to decrease expression in SHR animals, no genes were significantly decreased in AI animals compared to AU animals.
In order to provide some validation of the findings, RT-qPCR was performed on 5 IEGs from a subset of young (n = 9), AI (n = 6), and AU (n = 6) animals. The animals were selected based on set shifting performance to insure group differences in behavior (Figure 8) and an ANOVA confirmed a group difference [F(2, 18) = 17.15, p < 0.0001]. Post-hoc tests confirmed that AI animals exhibited an increase in the TTC relative to young and AU animals. Figure 9 shows a comparison of RT-qPCR results relative to the gene counts for the same genes and animals. The genes selected were IEGs that were increased in impaired animals (Arc, Egr1, Egr2, Egr4, Fos) and the comparisons (i.e., t-tests) were limited to confirmation of results observed for the whole data set. For both the RNA-seq and RT-qPCR measures of expression, t-tests indicated a difference in expression (p < 0.05) between the subset of AI and AU animals (Figure 9). In addition, Lin7b and Egr4 were expected to decrease with age, which was confirmed for the RNA-seq and RT-qPCR using t-tests to compare the subset of young and aged animals (Figure 9).
Figure 8. Set shifting performance for animals used in RT-qPCR validations. The animals were selected based on set shifting performance to insure group differences in behavior. The bars illustrate this difference as the mean + SEM TTC for young (n = 9) and aged animals classified as unimpaired (AU, n = 6) and impaired (AI, n = 6) on the set shift task. Asterisks indicate a significant (p < 0.05) difference relative to AI animals.
Figure 9. Comparison between RT-qPCR and RNA-seq. Six genes were selected for validation experiments using a subset of animals. Each panel provides the mPFC expression determined by RT-qPCR (left, ΔΔCT values) and RNA-seq (right, counts). Two-tailed t-tests confirmed increased expression of Arc, Fos, Egr1, Egr2, and Egr4 in AI, relative to AU rats. Gene expression for young animals is provided for comparison to aged animals. For two genes, Lin7b and Egr4, age differences were confirmed (***p < 0.005, **p < 0.025, *p < 0.05).
Neuroimaging indicates that the pattern of cognitive decline is related to the structure or activity of specific brain regions (Grady et al., 2005; Persson et al., 2006; Dennis et al., 2008; Park and Reuter-Lorenz, 2009; Migo et al., 2016). Furthermore, altered white matter integrity could influence connectivity between brain regions (O'Sullivan et al., 2001; Pfefferbaum et al., 2005; Salat et al., 2005; Andrews-Hanna et al., 2007; Bennett et al., 2011; Lu et al., 2011; Borghesani et al., 2013). While similar biological processes were altered across regions with age, very few genes were similarly affected across regions. Dissimilarities may relate to regional differences in vulnerability. Furthermore, vulnerability to aging is influenced by environment and lifestyle such that food restriction to promote operant behavior as well as the training procedure, may have differentially modified aging processes across regions (Lee et al., 2000; Zeier et al., 2011). Nevertheless, RNA-seq profiles confirmed that brain aging is associated with biological processes that involve increased expression of immune/defense response genes and decreased mitochondria and neuronal/synaptic genes (Prolla, 2002; Blalock et al., 2003; Lu et al., 2004; Verbitsky et al., 2004; Bordner et al., 2011; VanGuilder et al., 2011; Zeier et al., 2011; Cribbs et al., 2012; Berchtold et al., 2013; Primiani et al., 2014).
Age-related differences in cortical transcription are observed across species possibly due to evolutionary constraints on aging, examination of disparate brain regions, or differences in the age range examined (Erraji-Benchekroun et al., 2005; Fraser et al., 2005; Loerch et al., 2008; Cribbs et al., 2012; Berchtold et al., 2013). The mPFC of rats is thought to be functionally related to the dorsolateral PFC in humans and monkeys (Dias et al., 1997; Birrell and Brown, 2000; Kesner and Churchwell, 2011). Similar to the PFC of humans, we found that genes linked to excitatory and inhibitory transmitter systems decline with age in the mPFC. Significantly, age-related changes in gene expression did not predict cognition.
An important contribution of the current research was the specificity of transcriptional changes that correlate with a cognitive function that depends on the mPFC. The mPFC contributes to cognitive flexibility and spatial learning (Churchwell et al., 2010); however, no correlation was observed between set shifting and acquisition of a spatial search strategy (Barense et al., 2002; Beas et al., 2013). Transcription in the mPFC reflected this distinction in that genes correlated with the discrimination index generally decreased expression (62%) in impaired animals, while a large proportion (73%) of genes associated with impaired set shifting exhibited increasing expression. Specificity of mPFC transcription was also reflected in the 3-fold increase in the number of genes that correlated with set shift performance relative to visual discrimination learning. Finally, many of the activity-related IEGs that increased in impaired animals, exhibit down regulation in the mPFC during aging; and down regulation of IEGs in other brain regions is associated with cognitive impairment and decreased responsiveness (Benloucif et al., 1997; Blalock et al., 2003; Rowe et al., 2007). The results indicate specificity of mPFC transcription with an age-related impairment of mPFC-dependent behavior.
Altered basal expression of PFC IEGs is observed across rat species that exhibit differences in executive function. Upregulation of neural activity and synaptic plasticity genes is observed for SHR rats, which exhibit impaired set shift behavior compared to Wistar-Kyoto rats (Qiu et al., 2010). Similarly, compared to successful aging of LOU/C/Jall rats, aging in Wistar rats is associated with working memory deficits and an increase in PFC expression of IEGs (Arc, Egr2, Fos, Junb, and Nr4a1; Paban et al., 2013). Differences in transcription across rodent strains may result from genetic polymorphisms. Thus, an important finding from the current study is that individual variability in cognition, within the same species, is associated with increased basal expression of mPFC IEGs indicating that increased expression is indicative of impaired cognitive flexibility during aging.
What mechanism could increase IEG expression in animals that delay set shift behavior? Differences in IEG expression could result from differences in epigenetic regulation of transcription during aging (Peleg et al., 2010; Hernandez et al., 2011). Indeed, set shifting behavior was associated with increased expression of genes involved in regulating histone deacetylase (Bcor, Ccdc101, Dnajb5, Kdm6b) and histone acetyltransferase (Ing3) activity. Alternatively, IEG expression is upregulated by increased neuronal activity (Ghosh et al., 1994; Guan et al., 2009; Rudenko et al., 2013). Impaired set shift behavior was correlated with increased expression of transcription factors of inhibitory neurons (Dlx1, Npas1) and genes linked to the strength of excitatory and inhibitory inputs (Arc, Npas4) suggesting altered synaptic plasticity and increased neural activity. Interestingly, an increase in frontal cortex neural activity is observed in older humans and may relate to performance of cognitive tasks (Rosano et al., 2005; Turner and Spreng, 2012; Maillet and Rajah, 2014).
One possible mechanism for impaired set shift performance and increased neuronal activity involves a decline in N-methyl-D-aspartate (NMDA) receptor function. The degree of mPFC NMDA receptor hypofunction is correlated with impaired attention starting in middle-age (Guidi et al., 2015). Furthermore, blockade of NMDA receptors in this region disrupts set shift behavior (Stefani and Moghaddam, 2005; Dalton et al., 2011). The mechanism appears to involve a shift in the balance of excitatory/inhibitory synaptic input since inhibition of NMDA receptors in the PFC decreases the activity of inhibitory interneurons and increases the discharge activity of pyramidal cells (Homayoun and Moghaddam, 2007). Thus, an age-related decline in NMDA receptor function may reduce inhibitory drive and increase expression of activity-related genes in pyramidal cells.
In summary, the results support the idea that aging is associated with an increase in expression of immune and defense response genes and a decline in synaptic and neural activity genes. Importantly, mPFC expression of IEGs related to neural activity and synaptic plasticity decline with age; however, expression is up regulated in aged animals that exhibit delayed set shift behavior. The mPFC transcriptional profile of impaired animals is in contrast to decreased IEG expression reported for the hippocampus and other brain regions during aging. The specificity of impairment on a mPFC-dependent task, associated with a particular mPFC transcriptional profile indicates that impaired executive function involves altered transcriptional regulation and neural activity/plasticity processes that are distinct from that described for impaired hippocampal function.
TF, AK, LI designed experiment, AR, LI, AK, and BB performed research, TF and LI analyzed data, constructed figures, and wrote paper
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Financial support by National Institutes of Aging Grants R01AG037984, R37AG036800, R01AG49711, and the Evelyn F. McKnight Brain Research Foundation is highly appreciated. Thanks to Dr. Jennifer Bizon for discussions about the behavioral tasks.
Aenlle, K. K., and Foster, T. C. (2010). Aging alters the expression of genes for neuroprotection and synaptic function following acute estradiol treatment. Hippocampus 20, 1047–1060. doi: 10.1002/hipo.20703
Aenlle, K. K., Kumar, A., Cui, L., Jackson, T. C., and Foster, T. C. (2009). Estrogen effects on cognition and hippocampal transcription in middle-aged mice. Neurobiol. Aging 30, 932–945. doi: 10.1016/j.neurobiolaging.2007.09.004
Andrews-Hanna, J. R., Snyder, A. Z., Vincent, J. L., Lustig, C., Head, D., Raichle, M. E., et al. (2007). Disruption of large-scale brain systems in advanced aging. Neuron 56, 924–935. doi: 10.1016/j.neuron.2007.10.038
Barense, M. D., Fox, M. T., and Baxter, M. G. (2002). Aged rats are impaired on an attentional set-shifting task sensitive to medial frontal cortex damage in young rats. Learn. Mem. 9, 191–201. doi: 10.1101/lm.48602
Benloucif, S., Masana, M. I., and Dubocovich, M. L. (1997). Light-induced phase shifts of circadian activity rhythms and immediate early gene expression in the suprachiasmatic nucleus are attenuated in old C3H/HeN mice. Brain Res. 747, 34–42. doi: 10.1016/S0006-8993(96)01182-1
Bennett, I. J., Madden, D. J., Vaidya, C. J., Howard, J. H. Jr., and Howard, D. V. (2011). White matter integrity correlates of implicit sequence learning in healthy aging. Neurobiol. Aging 32, 2317.e1-12. doi: 10.1016/j.neurobiolaging.2010.03.017
Berchtold, N. C., Coleman, P. D., Cribbs, D. H., Rogers, J., Gillen, D. L., and Cotman, C. W. (2013). Synaptic genes are extensively downregulated across multiple brain regions in normal human aging and Alzheimer's disease. Neurobiol. Aging 34, 1653–1661. doi: 10.1016/j.neurobiolaging.2012.11.024
Bizon, J. L., Foster, T. C., Alexander, G. E., and Glisky, E. L. (2012). Characterizing cognitive aging of working memory and executive function in animal models. Front. Aging Neurosci. 4:19. doi: 10.3389/fnagi.2012.00019
Blalock, E. M., Chen, K. C., Sharrow, K., Herman, J. P., Porter, N. M., Foster, T. C., et al. (2003). Gene microarrays in hippocampal aging: statistical profiling identifies novel processes correlated with cognitive impairment. J. Neurosci. 23, 3807–3819.
Bordner, K. A., Kitchen, R. R., Carlyle, B., George, E. D., Mahajan, M. C., Mane, S. M., et al. (2011). Parallel declines in cognition, motivation, and locomotion in aging mice: association with immune gene upregulation in the medial prefrontal cortex. Exp. Gerontol. 46, 643–659. doi: 10.1016/j.exger.2011.03.003
Borghesani, P. R., Madhyastha, T. M., Aylward, E. H., Reiter, M. A., Swarny, B. R., Schaie, K. W., et al. (2013). The association between higher order abilities, processing speed, and age are variably mediated by white matter integrity during typical aging. Neuropsychologia 51, 1435–1444. doi: 10.1016/j.neuropsychologia.2013.03.005
Burger, C., López, M. C., Feller, J. A., Baker, H. V., Muzyczka, N., and Mandel, R. J. (2007). Changes in transcription within the CA1 field of the hippocampus are associated with age-related spatial learning impairments. Neurobiol. Learn. Mem. 87, 21–41. doi: 10.1016/j.nlm.2006.05.003
Churchwell, J. C., Morris, A. M., Musso, N. D., and Kesner, R. P. (2010). Prefrontal and hippocampal contributions to encoding and retrieval of spatial memory. Neurobiol. Learn. Mem. 93, 415–421. doi: 10.1016/j.nlm.2009.12.008
Cribbs, D. H., Berchtold, N. C., Perreau, V., Coleman, P. D., Rogers, J., Tenner, A. J., et al. (2012). Extensive innate immune gene activation accompanies brain aging, increasing vulnerability to cognitive decline and neurodegeneration: a microarray study. J. Neuroinflammation 9:179. doi: 10.1186/1742-2094-9-179
Dalton, G. L., Ma, L. M., Phillips, A. G., and Floresco, S. B. (2011). Blockade of NMDA GluN2B receptors selectively impairs behavioral flexibility but not initial discrimination learning. Psychopharmacology 216, 525–535. doi: 10.1007/s00213-011-2246-z
De Bruin, N. M., Kiliaan, A. J., De Wilde, M. C., and Broersen, L. M. (2003). Combined uridine and choline administration improves cognitive deficits in spontaneously hypertensive rats. Neurobiol. Learn. Mem. 80, 63–79. doi: 10.1016/S1074-7427(03)00024-8
de Magalhães, J. P., Curado, J., and Church, G. M. (2009). Meta-analysis of age-related gene expression profiles identifies common signatures of aging. Bioinformatics 25, 875–881. doi: 10.1093/bioinformatics/btp073
Dennis, N. A., Hayes, S. M., Prince, S. E., Madden, D. J., Huettel, S. A., and Cabeza, R. (2008). Effects of aging on the neural correlates of successful item and source memory encoding. J. Exp. Psychol. Learn. Mem. Cogn. 34, 791–808. doi: 10.1037/0278-73184.108.40.2061
Dias, R., Robbins, T. W., and Roberts, A. C. (1997). Dissociable forms of inhibitory control within prefrontal cortex with an analog of the Wisconsin Card Sort Test: restriction to novel situations and independence from “on-line” processing. J. Neurosci. 17, 9285–9297.
Erraji-Benchekroun, L., Underwood, M. D., Arango, V., Galfalvy, H., Pavlidis, P., Smyrniotopoulos, P., et al. (2005). Molecular aging in human prefrontal cortex is selective and continuous throughout adult life. Biol. Psychiatry 57, 549–558. doi: 10.1016/j.biopsych.2004.10.034
Floresco, S. B., Block, A. E., and Tse, M. T. (2008). Inactivation of the medial prefrontal cortex of the rat impairs strategy set-shifting, but not reversal learning, using a novel, automated procedure. Behav. Brain Res. 190, 85–96. doi: 10.1016/j.bbr.2008.02.008
Foster, T. C. (2012). Dissecting the age-related decline on spatial learning and memory tasks in rodent models: N-methyl-D-aspartate receptors and voltage-dependent Ca2+ channels in senescent synaptic plasticity. Prog. Neurobiol. 96, 283–303. doi: 10.1016/j.pneurobio.2012.01.007
Grady, C. L., McIntosh, A. R., and Craik, F. I. (2005). Task-related activity in prefrontal cortex and its relation to recognition memory performance in young and old adults. Neuropsychologia 43, 1466–1481. doi: 10.1016/j.neuropsychologia.2004.12.016
Guan, J. S., Haggarty, S. J., Giacometti, E., Dannenberg, J. H., Joseph, N., Gao, J., et al. (2009). HDAC2 negatively regulates memory formation and synaptic plasticity. Nature 459, 55–60. doi: 10.1038/nature07925
Guidi, M., Kumar, A., and Foster, T. C. (2015). Impaired attention and synaptic senescence of the prefrontal cortex involves redox regulation of NMDA receptors. J. Neurosci. 35, 3966–3977. doi: 10.1523/JNEUROSCI.3523-14.2015
Guidi, M., Kumar, A., Rani, A., and Foster, T. C. (2014). Assessing the emergence and reliability of cognitive decline over the life span in Fisher 344 rats using the spatial water maze. Front. Aging Neurosci. 6:2. doi: 10.3389/fnagi.2014.00002
Haberman, R. P., Colantuoni, C., Stocker, A. M., Schmidt, A. C., Pedersen, J. T., and Gallagher, M. (2011). Prominent hippocampal CA3 gene expression profile in neurocognitive aging. Neurobiol. Aging 32, 1678–1692. doi: 10.1016/j.neurobiolaging.2009.10.005
Hernandez, D. G., Nalls, M. A., Gibbs, J. R., Arepalli, S., van der Brug, M., Chong, S., et al. (2011). Distinct DNA methylation changes highly correlated with chronological age in the human brain. Hum. Mol. Genet. 20, 1164–1172. doi: 10.1093/hmg/ddq561
Homayoun, H., and Moghaddam, B. (2007). NMDA receptor hypofunction produces opposite effects on prefrontal cortex interneurons and pyramidal neurons. J. Neurosci. 27, 11496–11500. doi: 10.1523/JNEUROSCI.2213-07.2007
Huang, D. W., Sherman, B. T., Tan, Q., Collins, J. R., Alvord, W. G., Roayaei, J., et al. (2007a). The DAVID gene functional classification tool: a novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol. 8:R183. doi: 10.1186/gb-2007-8-9-r183
Huang, D. W., Sherman, B. T., Tan, Q., Kir, J., Liu, D., Bryant, D., et al. (2007b). DAVID bioinformatics resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 35, W169–W175. doi: 10.1093/nar/gkm415
Jackson, T. C., Rani, A., Kumar, A., and Foster, T. C. (2009). Regional hippocampal differences in AKT survival signaling across the lifespan: implications for CA1 vulnerability with aging. Cell Death Differ. 16, 439–448. doi: 10.1038/cdd.2008.171
Kadish, I., Thibault, O., Blalock, E. M., Chen, K. C., Gant, J. C., Porter, N. M., et al. (2009). Hippocampal and cognitive aging across the lifespan: a bioenergetic shift precedes and increased cholesterol trafficking parallels memory impairment. J. Neurosci. 29, 1805–1816. doi: 10.1523/JNEUROSCI.4599-08.2009
Liao, Y., Smyth, G. K., and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930. doi: 10.1093/bioinformatics/btt656
Lipinski, M. M., Zheng, B., Lu, T., Yan, Z., Py, B. F., Ng, A., et al. (2010). Genome-wide analysis reveals mechanisms modulating autophagy in normal brain aging and in Alzheimer's disease. Proc. Natl. Acad. Sci. U.S.A. 107, 14164–14169. doi: 10.1073/pnas.1009485107
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Loerch, P. M., Lu, T., Dakin, K. A., Vann, J. M., Isaacs, A., Geula, C., et al. (2008). Evolution of the aging brain transcriptome and synaptic regulation. PLoS ONE 3:e3329. doi: 10.1371/journal.pone.0003329
Lu, P. H., Lee, G. J., Raven, E. P., Tingus, K., Khoo, T., Thompson, P. M., et al. (2011). Age-related slowing in cognitive processing speed is associated with myelin integrity in a very healthy elderly sample. J. Clin. Exp. Neuropsychol. 33, 1059–1068. doi: 10.1080/13803395.2011.595397
Maillet, D., and Rajah, M. N. (2014). Age-related differences in brain activity in the subsequent memory paradigm: a meta-analysis. Neurosci. Biobehav. Rev. 45, 246–257. doi: 10.1016/j.neubiorev.2014.06.006
Migo, E. M., O'Daly, O., Mitterschiffthaler, M., Antonova, E., Dawson, G. R., Dourish, C. T., et al. (2016). Investigating virtual reality navigation in amnestic mild cognitive impairment using fMRI. Neuropsychol. Dev. Cogn. B Aging Neuropsychol. Cogn. 23, 196–217. doi: 10.1080/13825585.2015.1073218
O'Sullivan, M., Jones, D. K., Summers, P. E., Morris, R. G., Williams, S. C., and Markus, H. S. (2001). Evidence for cortical “disconnection” as a mechanism of age-related cognitive decline. Neurology 57, 632–638. doi: 10.1212/WNL.57.4.632
Paban, V., Billard, J. M., Bouet, V., Freret, T., Boulouard, M., Chambon, C., et al. (2013). Genomic transcriptional profiling in LOU/C/Jall rats identifies genes for successful aging. Brain Struct. Funct. 218, 1501–1512. doi: 10.1007/s00429-012-0472-8
Peleg, S., Sananbenesi, F., Zovoilis, A., Burkhardt, S., Bahari-Javan, S., Agis-Balboa, R. C., et al. (2010). Altered histone acetylation is associated with age-dependent memory impairment in mice. Science 328, 753–756. doi: 10.1126/science.1186088
Persson, J., Nyberg, L., Lind, J., Larsson, A., Nilsson, L. G., Ingvar, M., et al. (2006). Structure-function correlates of cognitive decline in aging. Cereb. Cortex 16, 907–915. doi: 10.1093/cercor/bhj036
Pfefferbaum, A., Adalsteinsson, E., and Sullivan, E. V. (2005). Frontal circuitry degradation marks healthy adult aging: evidence from diffusion tensor imaging. Neuroimage 26, 891–899. doi: 10.1016/j.neuroimage.2005.02.034
Primiani, C. T., Ryan, V. H., Rao, J. S., Cam, M. C., Ahn, K., Modi, H. R., et al. (2014). Coordinated gene expression of neuroinflammatory and cell signaling markers in dorsolateral prefrontal cortex during human brain development and aging. PLoS ONE 9:e110972. doi: 10.1371/journal.pone.0110972
Qiu, J., Hong, Q., Chen, R. H., Tong, M. L., Zhang, M., Fei, L., et al. (2010). Gene expression profiles in the prefrontal cortex of SHR rats by cDNA microarrays. Mol. Biol. Rep. 37, 1733–1740. doi: 10.1007/s11033-009-9596-1
Rosano, C., Aizenstein, H. J., Cochran, J. L., Saxton, J. A., De Kosky, S. T., Newman, A. B., et al. (2005). Event-related functional magnetic resonance imaging investigation of executive control in very old individuals with mild cognitive impairment. Biol. Psychiatry 57, 761–767. doi: 10.1016/j.biopsych.2004.12.031
Rowe, W. B., Blalock, E. M., Chen, K. C., Kadish, I., Wang, D., Barrett, J. E., et al. (2007). Hippocampal expression analyses reveal selective association of immediate-early, neuroenergetic, and myelinogenic pathways with cognitive impairment in aged rats. J. Neurosci. 27, 3098–3110. doi: 10.1523/JNEUROSCI.4163-06.2007
Rudenko, A., Dawlaty, M. M., Seo, J., Cheng, A. W., Meng, J., Le, T., et al. (2013). Tet1 is critical for neuronal activity-regulated gene expression and memory extinction. Neuron 79, 1109–1122. doi: 10.1016/j.neuron.2013.08.003
Sagvolden, T., Russell, V. A., Aase, H., Johansen, E. B., and Farshbaf, M. (2005). Rodent models of attention-deficit/hyperactivity disorder. Biol. Psychiatry 57, 1239–1247. doi: 10.1016/j.biopsych.2005.02.002
Salat, D. H., Tuch, D. S., Hevelone, N. D., Fischl, B., Corkin, S., Rosas, H. D., et al. (2005). Age-related changes in prefrontal white matter measured by diffusion tensor imaging. Ann. N.Y. Acad. Sci. 1064, 37–49. doi: 10.1196/annals.1340.009
Stefani, M. R., and Moghaddam, B. (2005). Systemic and prefrontal cortical NMDA receptor blockade differentially affect discrimination learning and set-shift ability in rats. Behav. Neurosci. 119, 420–428. doi: 10.1037/0735-7044.119.2.420
Turner, G. R., and Spreng, R. N. (2012). Executive functions and neurocognitive aging: dissociable patterns of brain activity. Neurobiol. Aging 33, 826.e1-e13. doi: 10.1016/j.neurobiolaging.2011.06.005
Uylings, H. B., and van Eden, C. G. (1990). Qualitative and quantitative comparison of the prefrontal cortex in rat and in primates, including humans. Prog. Brain Res. 85, 31–62. doi: 10.1016/S0079-6123(08)62675-8
VanGuilder, H. D., Bixler, G. V., Brucklacher, R. M., Farley, J. A., Yan, H., Warrington, J. P., et al. (2011). Concurrent hippocampal induction of MHC II pathway components and glial activation with advanced aging is not correlated with cognitive impairment. J. Neuroinflammation 8:138. doi: 10.1186/1742-2094-8-138
Verbitsky, M., Yonan, A. L., Malleret, G., Kandel, E. R., Gilliam, T. C., and Pavlidis, P. (2004). Altered hippocampal transcript profile accompanies an age-related spatial memory deficit in mice. Learn. Mem. 11, 253–260. doi: 10.1101/lm.68204
Yuan, Y., Chen, Y. P., Boyd-Kirkup, J., Khaitovich, P., and Somel, M. (2012). Accelerated aging-related transcriptome changes in the female prefrontal cortex. Aging Cell 11, 894–901. doi: 10.1111/j.1474-9726.2012.00859.x
Keywords: aging, cognitive flexibility, prefrontal cortex, set shifting task, transcription
Citation: Ianov L, Rani A, Beas BS, Kumar A and Foster TC (2016) Transcription Profile of Aging and Cognition-Related Genes in the Medial Prefrontal Cortex. Front. Aging Neurosci. 8:113. doi: 10.3389/fnagi.2016.00113
Received: 08 March 2016; Accepted: 29 April 2016;
Published: 17 May 2016.
Edited by:P. Hemachandra Reddy, Texas Tech University, USA
Reviewed by:Corinna Burger, University of Wisconsin-Madison, USA
Willard Freeman, Oklahoma Health Sciences University, USA
Copyright © 2016 Ianov, Rani, Beas, Kumar and Foster. 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.