<title>Introduction</title> <p>A decline in cell function with advancing age is a ubiquitous characteristic of all organisms. In humans, the effects of aging become manifest on a variety of levels that extend from an accumulation of DNA mutations to lipid oxidation, protein modification, cell loss, and ultimately death that is primarily due to increased susceptibility to age-related diseases ##REF##11972154##[1]##. Apart from overt changes, such as declining muscle strength, extensive metabolic alterations also occur with aging, one of the most prominent being impaired glucose tolerance ##REF##9334293##[2]##. Two central evolutionary theories hypothesize that the detrimental effects of aging are due to an accumulation of mutations, or antagonistic pleiotropy, whereby genes with beneficial effects early in life become deleterious with age ##REF##12386342##[3]##. These are not necessarily exclusive, and there is at present relatively strong evidence for both in studies of model organisms and in natural populations ##REF##12386342##[3]##, ##REF##16618935##[4]##.</p> <p>Whatever the cause, there is value in charting the molecular sequelae of aging on as broad a scale as possible. Few other methodological approaches lend themselves as well to this as mRNA expression profiling. There have been a handful of studies that have attempted to catalogue how mRNA expression changes with age, the largest of which have been performed in kidney ##REF##15562319##[5]## and muscle samples ##REF##16789832##[6]##. An intriguing conclusion from the latter study is that there may be a common set of genes that change equivalently in different tissues. For example genes that make up the mitochondrial electron transport chain appear to decrease with age in different tissues, and this is supported in that decreases are also evident in mice and flies ##REF##16789832##[6]##. Importantly however, these studies remain relatively small in scale and few in number thus meriting larger studies in additional populations and tissues.</p> <p>The effects of aging are particularly pronounced in the human brain where characteristic changes in morphology include a reduction in both neuronal size and synaptic density ##REF##3606042##[7]##, ##REF##8423884##[8]##. On a behavioral level, decreases in motor and cognitive function are hallmarks of normal aging ##UREF##0##[9]##. Dementia is the most prevalent disorder of the human brain, affecting 20–25 million people worldwide and representing a tremendous burden in terms of years of suffering with disability ##REF##16360788##[10]##. To date, there has been only one study, involving 30 individuals, on the effects of aging on gene expression at a transcriptome-level in the normal human brain ##REF##15190254##[11]##. There the authors described a number of pathways affected by aging and also concluded that genes for which down-regulation was evident had an over-representation of mutation in gene promoters. A recent study examined the effects of gene polymorphism on gene expression in a relatively large set of human brain samples ##REF##17982457##[12]##. An attractive feature of that study was that individuals, prior to death, were free of neurological disease. Our primary focus in the present study has been to use the above sample to investigate the question of whether general changes in gene expression occur in the aging human brain. This has been complemented with an analysis of lymphocyte mRNA expression in order to explore for common molecular themes in different tissues as well as to enable an assessment of a broader age range than that typically available in post-mortem samples.</p>
<title>Materials and Methods</title> <title>Human Samples</title> <p>The primary sample used in this study consisted of 191 individuals with ages-at-death data ranging from 65–100 years for which brain autopsy tissue was obtained. Expression level data were obtained using the Illumina HumanRefSeq-8 expression BeadChip platform for a total of 14,078 transcripts in which expression was detected in 5 or more individuals. Detailed descriptions of the human samples as well as the expression profiling protocol are provided in the original publication for which this sample was presented ##REF##17982457##[12]##. To create a working dataset, expression phenotypes and covariate data were merged that included age-at-death, post-mortem interval (pmi), gender, brain region, and transcript detection rates. From this, we generated 2 additional covariates for 1) total average expression level from all transcripts and 2) total average expression level for all transcripts that were detected in all individuals (5269 transcripts fulfilled the latter criteria). All individual transcript levels, the global average transcript levels, and pmi were log<sub>10</sub> transformed prior to inclusion in analyses.</p> <p>The second sample consisted of 1240 individuals ranging in age from 15-94 in which fresh blood lymphocytes had been obtained and mRNA extracted. Expression phenotyping was conducted using the Illumina Sentrix Human Whole Genome (WG-6) Series I BeadChip platform. In total there were 19,648 individual transcript measures, representing approximately 18,500 genes. Details of the sample and expression protocols have been provided previously in the original publication describing this sample ##REF##17873875##[15]##. For our analysis, we used the normalized expression phenotypes without further transformation.</p> <title>Correlation of Age and Transcript Level</title> <p>We classified the extremes as outliers that are expected to be observed once or more in 1930 individual log transformed transcript estimates with assumption of independence and normal distribution of 193 measures for each transcript.where <italic>α′</italic> is the alpha level for each measure, <italic>α</italic> is 0.05 and <italic>k</italic> 193. Based on normal distribution assumption,where Φ is standard normal cumulative distribution fuction, z is z-score of each measure.</p> <p>The normality of log transformed observed values for each transcript was tested by means of a Shapiro-Wilk W test. To assess differences of transcript detection rates or global expression levels across brain regions, as well as to test differences in characteristics of genomic architecture between age-related gene sets, ANOVA was used. Contingency tables (for example comparing the proportions of housekeeping genes) were evaluated by means of a standard chi-square test. All statistical analyses not related to pathway based tests were performed using STATA se 9.0. For the brain sample, the dependency of expression level on age-at-death was tested by fitting a linear regression model for each transcript:where <italic>Y</italic> is the log<sub>10</sub> transformed expression level of the transcript, <italic>x</italic>
<sub>1</sub> age-at-death, <italic>x</italic>
<sub>2</sub> log<sub>10</sub> of global expression, <italic>β<sub>0</sub></italic> intercept, <italic>β<sub>1</sub></italic> and <italic>β<sub>2</sub></italic> slope of each independent variable, and <italic>ε</italic> random error. For the lymphocyte sample, the same test was applied but without the global expression covariate according to the following:where <italic>Y</italic> is the normalized expression level of each transcript, <italic>x</italic>
<sub>1</sub> age, <italic>β<sub>0</sub></italic> intercept, <italic>β<sub>1</sub></italic> slope of age, and <italic>ε</italic> random error. We also allowed for non-independence of expression measures given the relatedness of family members in the lymphocyte sample ##REF##17468176##[34]##.</p> <title>Gene Ontology and KEGG Pathway Analyses</title> <p>The Entrez Gene IDs, symbols, and descriptions of genes for all tested transcripts were attained by Entrez Programming Utilities (eUtils) using GI number with the aid of a Perl script. This facilitated a search for replaced sequence identifiers and extracted information of interest by scanning the output text files from Entrez eUtils on the transcripts. Using the attained Entrez Gene IDs as identifiers of genes, we obtained the total lists of both Gene Ontology terms and KEGG pathway descriptors, with which our selected set of genes was annotated and analyzed using the DAVID functional annotation tool ##REF##12734009##[14]## with the most recent data update (January, 2008). The EASE scores ##REF##14519205##[35]## from the DAVID tool were used in trimming to over-represented term lists, with Bonferroni corrections of the scores from the tool as an inclusion criteria of ##TAB##2##Tables 3##, ##TAB##3##4##, ##TAB##6##7## and ##TAB##7##8##. We scrutinized ancestor-descendant relationships in the Gene Ontology structure among the terms in the enrichment lists on the basis of <italic>is_a</italic>, <italic>part_of</italic>, and <italic>regulates</italic> relationship by scanning the master ontology file which was updated in April, 2008. Our goal with this was to determine if ancestor terms had emerged as significant primarily because of enrichment of one of their descendant terms. We detail our strategy for this in ##SUPPL##0##Figure S1##. Among the GO terms in the list, those which had no descendants were labeled as ‘end-terms’. For every term excluding end-terms, we created the artificial descendants. These were intended to represent the complement of the set of genes annotated with a descendant term in the list with respect to all genes in the ancestor. The artificial term for every descendant of every ancestor in the list was then tested to determine if it was over-represented by applying a hypergeometric test ##REF##12620386##[36]##. The ancestors for which all artificial descendant terms were over-represented at P&lt;0.005 were labeled ‘significant ancestor terms’. The terms whose artificial descendants for end-terms or significant ancestors were over-represented were added to the set of ‘significant ancestor terms’. After modifying the set, the terms in the set were checked if they fulfilled inclusion criteria (enriched complement at P&lt;0.005) and the set was updated iteratively until there was no change. The end-terms and significant ancestor terms are listed in ##SUPPL##1##Tables S1##, ##SUPPL##2##S2##, ##SUPPL##3##S3##, ##SUPPL##4##S4##. The remaining terms in the enrichment list follow in the same tables that are significantly over-represented mainly by a highly enriched descendant. P-values of the enrichment test for the artificial terms are also shown in the ##SUPPL##1##Tables S1##, ##SUPPL##2##S2##, ##SUPPL##3##S3##, ##SUPPL##4##S4##.</p> <title>Genomic Architecture of Age-related Gene Sets</title> <p>The genomic positions of start and end points of transcripts, their coding sequences, and their exons and introns in both brain and lymphocyte gene sets were collected from the UCSC genome browser using Genbank accession numbers (e.g. NM_018711) as identifiers. Information on some of the transcripts that were detected in both samples was not available from the UCSC browser due to the records being suppressed. Therefore, the number of transcripts in the genomic architecture analysis was reduced to 12,520 and 12,242 for brain and lymphocyte samples, respectively. Based on the positional data, lengths of pre-mRNA sequences, coding sequence, total intron, 5′UTR, 3′UTR, average intron and average exon of each transcript, as well as the number of exons and the ratio of intron per coding sequence length were calculated and transformed on a log<sub>10</sub> scale. For comparative analyses of gene characteristics, the transcripts were divided into 3 groups, up-regulated genes (positively correlated with age at a P&lt;0.05 threshold), down-regulated (negatively correlated with age at a P&lt;0.05 threshold), and a set of genes that were not significantly altered over the age spectrum. Differences between groups were assessed using ANOVA. Combined analyses were performed on both sets using second order factorial ANOVA with sample source as a covariate.</p> <title>URLs</title> <p>Entrez Gene: <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/sites/entrezdbgene">http://www.ncbi.nlm.nih.gov/sites/entrezdbgene</ext-link>; Entrez eUtils: <ext-link ext-link-type="uri" xlink:href="http://eutils.ncbi.nlm.nih.gov/entrez/eutils">http://eutils.ncbi.nlm.nih.gov/entrez/eutils</ext-link>; DAVID 2008: <ext-link ext-link-type="uri" xlink:href="http://david.abcc.ncifcrf.gov/">http://david.abcc.ncifcrf.gov/</ext-link>; Gene Ontology (GO): <ext-link ext-link-type="uri" xlink:href="http://www.geneontology.org/">http://www.geneontology.org/</ext-link>; UCSC genome browser: <ext-link ext-link-type="uri" xlink:href="http://genome.ucsc.edu/">http://genome.ucsc.edu/</ext-link>
</p>
<title>Results</title> <p>We began by performing several validation steps towards the initial goal of evaluating if age-related change in gene expression could be detected in the human brain. First, evidence was sought for systematic outlier effects among the 2,096,975 individual expression level estimates. We identified all expression values that were more than 3.4644 standard deviations from the mean expression value for that particular transcript (see <xref ref-type="sec" rid="s4">methods</xref> for an explanation for choosing this threshold). With this, we expected to see approximately 1,086 outliers by chance in the entire set assuming normal distributions of the log<sub>10</sub> transformed data. There were 7,079 outliers identified in this way, and these were eliminated from all further analyses. Second, the focus of the original study was primarily upon brain cortex, with 5 subsets being represented, 3 cortical, one group of 6 cerebellum samples and an additional group of samples for which no specific region could be assigned. We assessed whether these classifications differed systematically with regards to transcript detection rates and global expression levels using ANOVA (see <xref ref-type="sec" rid="s4">methods</xref>). Across the various brain region categories there were no significant differences for any these covariates. Simple linear regression was also performed for age-at-death and pmi versus global expression level. For this, there was marginal evidence that age-at-death correlated with global expression (P = 0.080, r<sup>2</sup> = 0.016). For pmi however, the trend was significant, but the direction of the curve suggested that expression was increasing with increasing pmi (P = 0.0039. r<sup>2</sup> = 0.046). Of note, linear regression of the 2 global expression covariates (see <xref ref-type="sec" rid="s4">methods</xref>) and detection rate metrics (one reflecting all 24,357 transcripts probed by the Illumina HumanRefseq-8 chip and the other reflecting the 14,078 transcripts detected in this brain sample) all showed strong correlation (r<sup>2</sup> in excess of 0.5 for all 6 comparisons). Based upon all of the above tests, we resolved to use the global average expression level that includes all 14,078 transcripts as a primary metric of quality as a covariate for initial analyses.</p> <p>For the brain sample, we regressed age-at-death on transcript levels one at a time for all 14,078 transcripts that had been detected in 5 or more individuals in the original study. The results of this are shown in ##TAB##0##Tables 1## and ##TAB##1##2##, where we document 54 transcripts that exhibited significant change after Bonferroni correction for 14,078 tests. For purposes of illustration, regression lines for the top 6 scoring genes (3 up-regulated and 3 down-regulated) are also displayed in ##FIG##0##Figure 1##. For this analysis, the most significant finding was for the <italic>SVOP</italic> gene (encoding a synaptic vesicle protein), and the single strongest effect in terms of variance explained was for the <italic>TAC3</italic> gene (encoding tachykinin 3). The smallest number of individuals for which a significant correlation with age was detected was 133 (also for the <italic>TAC3</italic> gene). We noted that 3 of the top scoring 54 genes were also present in a list of 532 putative “housekeeping” genes (see ##REF##12850439##[13]##), these being <italic>PIN1</italic>, <italic>GUK1</italic>, and <italic>ARPC2</italic>. This proportion is not significantly different from what would be expected by chance (3 vs. 51 compared to 532 vs. 6968). At this stage, we also scrutinized this list of 54 genes more closely for the potential impact of covariates. For this, multiple regression models were fitted as above, but for each transcript in turn we also included terms for pmi, brain region, or gender. All 54 genes remained highly significant when any of these were individually tested (not shown). Among the top 54 transcripts, all were detectable in more than 100 individuals, reinforcing the importance of power in analyses such as this. There were however a few cases where a small number of observations gave rise to apparently large effects. The gene that ranked as the 139<sup>th</sup> most significant finding (falling below the multiple test correction threshold), <italic>ZIC3</italic>, had only 14 expression level observations and an r<sup>2</sup> = 0.76 (P = 2.6×10<sup>−5</sup>). This highlighted the possibility that genes with lower maximum expression levels might be changing to a larger degree with age, but technologies with lower detection thresholds and/or larger sample sizes would be required to identify them.</p> <p>We contrasted the individual transcript results of this analysis with the first study to address the impact of aging on expression at the transcriptome level ##REF##15190254##[11]##. In that study, 30 individuals were included with an age range of 26–106. Our set was scrutinized for the 148 genes previously highlighted to be associated with aging (see specifically ##TAB##0##table 1## from ##REF##15190254##[11]##). Only one gene from our set that survived strict multiple testing correction was also present in their set, this being <italic>MRPL28</italic>. With a relaxed uncorrected significance threshold in our set of 0.005, we found 26 overlapping genes. This proportion is significant given that we observed 1,141 transcripts in our set that were significant at the 0.005 threshold, suggesting that some replication exists, albeit limited (26 vs. 122 compared to 1,141 vs. 12,937). In our set, the majority of significant genes exhibited a decrease in expression with advancing age. An intriguing aspect of this comparison was that among the 26 overlapping genes, there were 9 for which expression increased with age in the original study, and for all 9 these were increased with age in the present study. Due to the consistency in the direction of the effects, this might be considered a stronger case for replication. Seven of these 9 genes were previously considered as a class of stress response genes.</p> <p>We explored an expanded list of genes that exceeded a relaxed significance threshold (uncorrected P&lt;0.05) for enrichment or deficit based upon Gene Ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The entire set of genes was divided into two groups based upon whether expression was increasing or decreasing with age. We then searched for over- or under-represented terms and pathways in each list of genes using the DAVID web application ##REF##12734009##[14]##. For this we created a diminished non-redundant set of 13,216 base genes from the total set of 14,078 by eliminating multiple transcripts representing the same gene (see <xref ref-type="sec" rid="s4">methods</xref>). The enriched terms (Bonferroni corrected <italic>P</italic>&lt;0.05) are shown in ##TAB##2##Tables 3## and ##TAB##3##4##, excluding the ancestor terms that became significant mainly by an overrepresented descendant in the GO structure. A larger list of all overrepresented terms and pathways with a less stringent threshold (uncorrected <italic>P</italic>&lt;0.01) is shown in ##SUPPL##1##Tables S1## and ##SUPPL##2##S2##. In summarizing these results, the genes whose products are involved in DNA-binding were the most significantly overrepresented group among genes whose expression increases with age (##TAB##3##Table 4##). Genes with products involved in the regulation of DNA-dependent transcription and genes encoding proteins located in or a subcomponent of the nucleus were also significantly enriched. Among all genes negatively correlated with age, the most significantly enriched groups included genes involved in nervous system development, mitochondrial genes, and those whose products are constituents of the cytoplasm (##TAB##2##Table 3##). We also highlighted the degree of enrichment in terms of a fold change (the number of genes that give rise to this number are shown in ##SUPPL##1##Tables S1## and ##SUPPL##2##S2##). Thus, in quantitative terms the most enriched genes occurred in the set of negatively correlated transcripts, where dendrite genes and genes associated with neuronal projections were enriched more than 3 fold (##TAB##2##Table 3##). In contrast, although highly significant in some cases, genes whose expression is increasing with age are enriched a maximum of 1.80 fold (##TAB##3##Table 4##). As a comparison, the probabilities for the enriched terms to be underrepresented by chance within the opposite group of genes are also shown in the last two columns of ##TAB##2##Tables 3## and ##TAB##3##4##. For example, from ##TAB##3##Table 4##, we show that the term “nucleus” which is significantly enriched in the set of genes that are increasing with age is also in deficit in the set of genes that are decreasing with age.</p> <p>We extended the above analytical scheme to a second sample consisting of 1240 individuals with lymphocyte mRNA measures for 19,648 transcripts representing 18,519 genes. We note that a rigorous standardization procedure had been applied to this set previously ##REF##17873875##[15]## and we thus elected not to perform any additional transformations, nor to include any additional covariates (see <xref ref-type="sec" rid="s4">methods</xref>). Linear regression was performed as above on all 19,648 transcripts versus age for this sample. Given its size, and compared to the brain sample, there were many more transcripts that exceeded a strict Bonferroni significance threshold, and we chose to present only the top 50 (25 negatively and 25 positively correlated with age) which are shown in ##TAB##4##Tables 5## and ##TAB##5##6##, arbitrarily truncated from a total of 1080 (612 negative and 468 positive). The top 3 scoring negatively regulated genes and top 3 scoring positively regulated genes are shown in ##FIG##1##Figure 2##. The most significant finding in this set was for the <italic>LRRN3</italic> (encoding a membrane protein with unknown function), and this was also the gene upon which age exhibited the largest effect (28% variance explained). Interestingly, the second largest effect observed in this set was 14% variance explained, representing in our view an anomalous drop from the top ranking gene. As was done for the brain sample, we again divided this list into up- and down- regulated genes and performed term/pathway based analyses. An arbitrary significance level of 0.01 was chosen from the linear regression results to establish these 2 lists (this differed from the brain sample for which a threshold of 0.05 was set due to the number of genes exceeding significance in lymphocytes and a 2000 gene limit set by the DAVID web tool). The total set of 19,648 transcripts for which regression was performed, was reduced to 13,231 for which annotations could be obtained towards analysis using the DAVID web application. The results of this are shown in ##TAB##6##Tables 7## and ##TAB##7##8## where Bonferroni corrected (P&lt;0.05) enriched terms and pathways are documented. As was done for the brain sample, an expanded list of enriched terms and pathways according to relaxed uncorrected significance threshold of P&lt;0.01 is shown in ##SUPPL##3##Tables S3## and ##SUPPL##4##S4##. The main highlight in our view is the highly significant enrichment of mitochondrial genes in the negatively regulated group (##TAB##6##Table 7##), providing a replication of what was seen in the brain. In general, the magnitude of the fold-change in term/pathway enrichment was larger in the set of negatively regulated genes, also in agreement with what was seen in brain.</p> <p>Based upon the identified lists of age-related genes for both brain and lymphocyte samples, we also explored for differences in basic genomic architecture. Our guide for this was a detailed analysis recently presented on characteristics of housekeeping genes ##REF##12850439##[13]##. The results of this analysis are presented in ##TAB##8##Tables 9## and ##TAB##9##10##, where we document the differences between the up and down-regulated sets in comparison with the characteristics of the non-regulated genes. In summarizing these results, while we found some striking differences, for example with coding sequence length in brain, we note that this was not replicated in lymphocytes. Instead we highlight a single category that does appear to be replicating, namely the ratio of average intron to average exon sequence length. We used second-order factorial ANOVA models to perform combined analyses for all variables. For this, the intron/exon sequence ratio is highly significant (F<sub>2,24158</sub> = 41.57, P = 9.5×10<sup>−19</sup>) with no evidence of interaction between the two sets (P = 0.69). This can be taken in context with the same analysis for coding sequence length, where the tissue by group (group defined by up, down, or non-regulated genes) interaction term was highly significant (P&lt;0.0001). For all other combined analyses, the main group effect was either non-significant or the interaction term was highly significant (P&lt;0.0001).</p> <p>Finally, we used these two samples to pursue the question of whether variance in gene expression itself changes with age. Our hypothesis was that variance might increase with age as a consequence of accumulating somatic mutation ##REF##15190254##[11]## and/or increasingly heterogeneous environmental exposures. For the brain sample, an F-test for equality of variances was conducted on all 14,078 transcripts by dividing the sample in two groups according to the median age. We observed only 10 genes that exceeded a strict multiple test correction threshold and in each case there was evidence that the distributions of these genes deviated from normality. There was no overlap in this set for the highest scoring candidates with genes found to change with age in the linear regression analyses. There were 6 cases where an age-related gene also had a borderline significant F-test (P&lt;0.05). This proportion was not significantly different than what could be found in the entire set (not shown). For the lymphocyte set, we took a slightly different approach due to its size and rather than dividing the sample at the median chose to examine decile bins and applied a Levene's robust variance test to explore for differences across age groups. There was again no strong evidence that variance differed across these groupings, the most significant finding being for the <italic>NLP</italic> gene (P = 2.3×10<sup>−6</sup>), and with only 2 genes in total attaining significance after multiple test correction. As a final closing note, we enumerated the number of statistical tests used in this study in its entirety, arriving at an approximate number of 110,000. This may be used as a reference for any of the un-corrected P-values that are presented.</p>
<title>Discussion</title> <p>We consider the most important finding in this study to be that expression levels of genes involved in mitochondrial processes are decreasing with age, as originally proposed by Zahn et al. and also supported by Miller et al. ##REF##16789832##[6]##, ##REF##18256261##[16]##. In the former study, a key observation was that this is apparent in multiple tissues as well as in species other than humans. The present analysis provides a confirmation of this as a common characteristic of aging, in that this evident in both brain and lymphocytes. For these two samples there are also tissue specific themes that have emerged that we consider a validation of the quality of the expression phenotypes originally obtained. In the brain sample, consistent with previous results ##REF##15190254##[11]##, a decrease in genes involved in synaptic function was observed, which follows from the documented changes in synaptic function that occur with age ##REF##8423884##[8]##. In lymphocytes, there was evidence that genes involved in the immune response increase expression with age, which might be regarded as a reflection of chronic persisting viruses such as cytomegalovirus (CMV) ##REF##16098721##[17]##. While we acknowledge central differences in these samples that include ethnicity, the age ranges, sample size, and the cell types represented, in this discussion we highlight these and other pathways with a comparison of the two different tissues as a guide.</p> <p>A key distinction in contrasting the study by Zahn et al. was the specific focus on genes of the mitochondrial respiratory chain ##REF##16789832##[6]##. While many of these are included in our emergent lists of genes of the mitochondrion (see ##SUPPL##1##Tables S1## and ##SUPPL##3##S3##), results indicate that the effects of age on mitochondrial function may be broader. For example, in the brain the most significant finding across all pathways was for cytoplasm genes (##TAB##2##Table 3##). Although this categorization includes all genes of the mitochondrion, and indeed our results suggest it is the latter grouping that contributes to the strength of the statistic, there remain additional genes that may be enriched in this set (see ##SUPPL##1##Table S1##). Oxidative phosphorylation itself did emerge as a highly significant term associated with aging in the brain sample, but not in lymphocytes, indicating that there are likely to be tissue specific differences in the component genes related to mitochondrial function (##TAB##2##Tables 3## and ##TAB##6##7##). Nonetheless, these data together with previous studies provide an intriguing foundation to investigate if this is a primary event in the aging process or derivative of a general decrease in energy metabolism and activity that accompanies old age (e.g. ##REF##17289456##[18]##).</p> <p>That transcript levels are both decreasing and increasing is an important indication that the aging process does not lead to a unidirectional decline in expression. While this was evident in both brain and lymphocytes, in each case the relative fold-enrichment of gene categories was higher among negatively correlated transcripts. Thus, the up-regulated categories consisted of more unique transcripts in general than those that were down-regulated. In brain we noted several pathways that appear to have component genes that increase in expression with advancing age, the foremost among these encompassing genes that encode nuclear proteins. The most significant categorizations that emerged included a large contingent of genes related to DNA-binding and transcription, the most abundant being zinc-finger proteins. Although genes related to transcription were highlighted previously, results suggested a mix of both positive and negative regulation ##REF##15190254##[11]##. In contrast however, Miller et al. did observe an overlap of genes related to transcription that were increasing in relation to both ageing and to Alzheimer's disease (AD) ##REF##18256261##[16]##.</p> <p>In lymphocytes, immune response was the most significant category associated with increasing expression. However, there are several issues that are impossible to resolve for this latter result, especially regarding which cell types are represented in this sample (e.g. ##REF##18222733##[19]##). For example, this might relate to an accumulation of highly differentiated T cells due to acceleration by persistent CMV infection but this would require extensive further study ##REF##16098721##[17]##. There were also some interesting highlights among the many individual genes that were up-regulated with age in lymphocytes. In particular the most significant up-regulated gene <italic>IGFBP3</italic> is intriguing in the context that another member of the insulin-like growth factor binding protein family, <italic>IGFBP7</italic> has recently been implicated as an inducer of apoptosis in human melanoma cell lines ##REF##18267069##[20]##, ##REF##18555773##[21]##. <italic>IGFBP3</italic> itself has also been shown to be increased in senescent human fibroblasts ##REF##1719537##[22]##. Those findings together with results of the present study provide support for the role of insulin-like growth factor signaling in cellular senescence.</p> <p>For the brain sample, we were intrigued by the emergent pattern that expression of cytoplasm genes may be decreasing with age, while expression of nuclear genes may be increasing with age. While this might reflect the large number of mitochondrial and transcription genes in these particular sets, we cannot ignore the possibility that a basic change in cell morphology might be at play. One explanation may be that neuronal number remains stable with age, but synaptic vesicle density decreases ##REF##3606042##[7]##, ##REF##8423884##[8]##. This relative change might give rise to our observations here. In support of this, genes related to neuronal projections (which includes dendrite genes) and synapse genes are major categories that appear to decrease with age (see ##TAB##2##Table 3##). Also in support of this, the highest ranking individual gene, <italic>SVOP</italic>, encoding the SV2 related protein, is localized to the synaptic vesicle and appears to be an ion transporter. We had anticipated that more neuronal specific genes would exhibit a decreasing pattern, but noted that even some common neuronal reference genes, such as <italic>ENO2</italic> (ranking 7597<sup>th</sup>), were not associated with age. Thus, in contrast to the enrichment of neuronal projection genes mentioned above (and see ##TAB##2##Table 3##), this might support the concept that a decrease in neuronal number is not a major feature of aging ##REF##16880353##[23]##. A potentially important source of confounding in analyses of post-mortem brain samples is the mode of death, which we have not examined in more detail in this study. In particular, agonal state and pH have been highlighted as contributors to brain mRNA expression ##REF##16636682##[24]##–##REF##18061308##[26]##. The present study has strong similarities with the study by Li et al., (2004)##REF##14734628##[25]##. However, the relationship between pH and age was previously explored, but there was no evidence of a significant correlation ##REF##16636682##[24]##, ##REF##14559355##[27]##. Interestingly, although it was a small sample Vawter et al. also noted significant correlations of age with mitochondrial genes regardless of agonal state (see specifically supplementary ##TAB##3##table 4## from ##REF##16636682##[24]##). Nonetheless, one explanation might be that transcription factor activation and mitochondrial deactivation is a natural response to hypoxic stress which may be more common in elderly individuals (e.g. ##REF##14734628##[25]##, ##REF##17624786##[28]##). Another possibility might be that both aging and extended hypoxia dependent on the mode of death share similarities in terms of gene regulation.</p> <p>For the lymphocyte sample, the most significantly enriched terms were found in the class of negatively regulated genes, where the ribosomal and translation machinery appears to be strongly affected by age. This could potentially contribute to “immunosenescence”, whereby the immune system decays with time, rendering elderly individuals more susceptible to infectious diseases ##REF##18400689##[29]##. While changes at the expression level for this class of genes have to our knowledge not been previously reported on this scale, the concept that protein synthesis is involved in immunosenescence is not new ##REF##3542542##[30]##. This result is in our view the most striking difference between lymphocytes and brain. Although this suggests a more general decreases in metabolism in aging lymphocytes, we consider this an important reminder that tissue specificity can play a major role in gene expression profiling ##REF##16480787##[31]##.</p> <p>Another important feature of this study is the relative degree of stability in gene expression across the age spectrum. To put this in context, the strongest effects of age on any genes were 18% variance explained for <italic>TAC3</italic> in brain and 29% for <italic>LRRN3</italic> in lymphocytes, but in examining ##TAB##0##Tables 1##, ##TAB##1##2##, ##TAB##4##5## and ##TAB##5##6##, it is apparent that effect sizes drop off rapidly. For each of these top-ranked genes there was an approximately 2-fold change in average expression levels over the age ranges. This is comparable to well documented changes in hormone levels with age, such as for IGF1 ##REF##9334293##[2]##, ##REF##8491152##[32]##. Also of note, genes claimed to be changing with age previously exhibited differences between young and older groups in the range of 2–3 fold ##REF##15190254##[11]##. This again might suggest that cell-loss and/or changes in organ function/morphology are more important than broad changes in gene expression in the decrease in health that typically accompanies increasing age. The possibility however cannot be excluded that relatively small changes in gene expression have a large impact on cell function. As an example, disorders such as Parkinson's disease (PD) can be caused by simple gene dosage effects, as has now been shown for the <italic>SNCA</italic> gene which encodes alpha-synuclein ##REF##15760766##[33]##. Although rare, we noted several genes that may be changing to a larger degree, but which fell below a strict significance threshold, for example <italic>ZIC3</italic> in brain. However, there were no similar cases in the much more powerful lymphocyte sample.</p> <p>Finally, we also explored basic gene architecture characteristics in an attempt to provide some insight into why expression might vary with age. The only parameter that emerged as significant and equivalent in both brain and lymphocytes was the intron/exon length ratio, this representing a metric of gene “compactness” ##REF##12850439##[13]##. This might be interpreted as suggesting that non-compact genes are more susceptible to mutations that disrupt regulation, and thus lead to decreasing expression with advancing age. However, we cannot exclude the possibility that what we are seeing is simply a result of the gene pathways that are over-represented, these being enriched for compact or non-compact genes. In other words, if specific pathways are affected by age, the genes that represent those pathways may have similar features. We consider this nonetheless an intriguing finding in that it is replicating in different tissues.</p> <p>In summary, we validate previous findings that a decrease in mitochondrial gene expression appears to be a common theme in the aging process. Whether this is a primary event that causes a decline in health with advancing age or a result of a general decrease in metabolism in the elderly remains a topic for further investigation. We also highlight additional novel pathways that may be age dependent but with dramatic differences between tissues, in particular with genes related to transcription and to translation. These results may provide a valuable foundation for understanding the molecular consequences of aging and emphasize the development of catalogues of senescence-related genes in additional tissues.</p>
<p>Analyzed the data: MGH AJM PKM JAP. Wrote the paper: MGH AJM PKM JAP.</p> <title>Background</title> <p>Identifying biological pathways that vary across the age spectrum can provide insight into fundamental mechanisms that impact disease and frailty in the elderly. Few methodological approaches offer the means to explore this question on as broad a scale as gene expression profiling. Here, we have evaluated mRNA expression profiles as a function of age in two populations; one consisting of 191 individuals with ages-at-death ranging from 65–100 years and with post-mortem brain mRNA measurements of 13,216 genes and a second with 1240 individuals ages 15–94 and lymphocyte mRNA estimates for 18,519 genes.</p> <title>Principal Findings</title> <p>Among negatively correlated transcripts, an enrichment of mitochondrial genes was evident in both populations, providing a replication of previous studies indicating this as a common signature of aging. Sample differences were prominent, the most significant being a decrease in expression of genes involved in translation in lymphocytes and an increase in genes involved in transcription in brain, suggesting that apart from energy metabolism other basic cell processes are affected by age but in a tissue-specific manner. In assessing genomic architecture, intron/exon sequence length ratios were larger among negatively regulated genes in both samples, suggesting that a decrease in the expression of non-compact genes may also be a general effect of aging. Variance in gene expression itself has been theorized to change with age due to accumulation of somatic mutations and/or increasingly heterogeneous environmental exposures, but we found no evidence for such a trend here.</p> <title>Significance</title> <p>Results affirm that deteriorating mitochondrial gene expression is a common theme in senescence, but also highlight novel pathways and features of gene architecture that may be important for understanding the molecular consequences of aging.</p>
<title>Supporting Information</title>
<fig id="pone-0003024-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003024.g001</object-id><label>Figure 1</label><caption><title>Scatter plots and fitted linear regression lines of the most significantly associated down-regulated (a–c) and up-regulated (d–f) genes with age-at-death in brain samples.</title></caption></fig> <fig id="pone-0003024-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003024.g002</object-id><label>Figure 2</label><caption><title>Scatter plots and fitted linear regression lines of the most significantly associated down-regulated (a–c) and up-regulated (d–f) genes with age in lymphocyte samples.</title></caption></fig>
<table-wrap id="pone-0003024-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003024.t001</object-id><label>Table 1</label><caption><title>Genes whose expression decreases with age in the human brain.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Symbol</td><td align="left" rowspan="1" colspan="1">Description</td><td align="left" rowspan="1" colspan="1">Chr.Pos</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic> value<xref ref-type="table-fn" rid="nt101">a</xref>
</td><td align="left" rowspan="1" colspan="1">R<sup>2 b</sup>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">SVOP</td><td align="left" rowspan="1" colspan="1">SV2 related protein homolog (rat)</td><td align="left" rowspan="1" colspan="1">12q</td><td align="left" rowspan="1" colspan="1">4.3(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">0.14</td></tr><tr><td align="left" rowspan="1" colspan="1">HBQ1</td><td align="left" rowspan="1" colspan="1">hemoglobin, theta 1</td><td align="left" rowspan="1" colspan="1">16p</td><td align="left" rowspan="1" colspan="1">8.6(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">0.13</td></tr><tr><td align="left" rowspan="1" colspan="1">MRPL28</td><td align="left" rowspan="1" colspan="1">mitochondrial ribosomal protein L28</td><td align="left" rowspan="1" colspan="1">16p</td><td align="left" rowspan="1" colspan="1">1.0(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.14</td></tr><tr><td align="left" rowspan="1" colspan="1">TAC3</td><td align="left" rowspan="1" colspan="1">tachykinin 3 (neuromedin K, neurokinin beta)</td><td align="left" rowspan="1" colspan="1">12q</td><td align="left" rowspan="1" colspan="1">1.0(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.18</td></tr><tr><td align="left" rowspan="1" colspan="1">C6orf154</td><td align="left" rowspan="1" colspan="1">chromosome 6 open reading frame 154</td><td align="left" rowspan="1" colspan="1">6p</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.13</td></tr><tr><td align="left" rowspan="1" colspan="1">SLC25A6</td><td align="left" rowspan="1" colspan="1">solute carrier family 25, member A6</td><td align="left" rowspan="1" colspan="1">Xp/Yp</td><td align="left" rowspan="1" colspan="1">1.9(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.14</td></tr><tr><td align="left" rowspan="1" colspan="1">ST3GAL2</td><td align="left" rowspan="1" colspan="1">ST3 beta-galactoside alpha-2,3-sialyltransferase 2</td><td align="left" rowspan="1" colspan="1">16q</td><td align="left" rowspan="1" colspan="1">2.2(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.15</td></tr><tr><td align="left" rowspan="1" colspan="1">PIN1</td><td align="left" rowspan="1" colspan="1">peptidylprolyl cis/trans isomerase, NIMA-interacting 1</td><td align="left" rowspan="1" colspan="1">19p</td><td align="left" rowspan="1" colspan="1">3.2(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.13</td></tr><tr><td align="left" rowspan="1" colspan="1">ARPC2</td><td align="left" rowspan="1" colspan="1">actin related protein 2/3 complex, subunit 2, 34 kDa</td><td align="left" rowspan="1" colspan="1">2q</td><td align="left" rowspan="1" colspan="1">3.2(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.14</td></tr><tr><td align="left" rowspan="1" colspan="1">KCNF1</td><td align="left" rowspan="1" colspan="1">potassium voltage-gated channel, subfamily F, member 1</td><td align="left" rowspan="1" colspan="1">2p</td><td align="left" rowspan="1" colspan="1">4.1(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">GSS</td><td align="left" rowspan="1" colspan="1">glutathione synthetase</td><td align="left" rowspan="1" colspan="1">20q</td><td align="left" rowspan="1" colspan="1">4.3(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">LOC255849</td><td align="left" rowspan="1" colspan="1">hypothetical LOC255849</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">4.5(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">HMGB3</td><td align="left" rowspan="1" colspan="1">high-mobility group box 3</td><td align="left" rowspan="1" colspan="1">Xq</td><td align="left" rowspan="1" colspan="1">4.8(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.13</td></tr><tr><td align="left" rowspan="1" colspan="1">TMEM121</td><td align="left" rowspan="1" colspan="1">transmembrane protein 121</td><td align="left" rowspan="1" colspan="1">14q</td><td align="left" rowspan="1" colspan="1">6.1(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.14</td></tr><tr><td align="left" rowspan="1" colspan="1">CAMK2N1</td><td align="left" rowspan="1" colspan="1">calcium/calmodulin-dependent protein kinase II inhibitor 1</td><td align="left" rowspan="1" colspan="1">1p</td><td align="left" rowspan="1" colspan="1">6.4(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">OLFM1</td><td align="left" rowspan="1" colspan="1">olfactomedin 1</td><td align="left" rowspan="1" colspan="1">9q</td><td align="left" rowspan="1" colspan="1">7.4(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.13</td></tr><tr><td align="left" rowspan="1" colspan="1">KCNIP1</td><td align="left" rowspan="1" colspan="1">Kv channel interacting protein 1</td><td align="left" rowspan="1" colspan="1">5q</td><td align="left" rowspan="1" colspan="1">9.2(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">FABP3</td><td align="left" rowspan="1" colspan="1">fatty acid binding protein 3</td><td align="left" rowspan="1" colspan="1">1p</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.08</td></tr><tr><td align="left" rowspan="1" colspan="1">GUK1</td><td align="left" rowspan="1" colspan="1">guanylate kinase 1</td><td align="left" rowspan="1" colspan="1">1q</td><td align="left" rowspan="1" colspan="1">1.2(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.13</td></tr><tr><td align="left" rowspan="1" colspan="1">COPS7A</td><td align="left" rowspan="1" colspan="1">COP9 constitutive photomorphogenic homolog subunit 7A (Arabidopsis)</td><td align="left" rowspan="1" colspan="1">12p</td><td align="left" rowspan="1" colspan="1">1.3(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.07</td></tr><tr><td align="left" rowspan="1" colspan="1">VIP</td><td align="left" rowspan="1" colspan="1">vasoactive intestinal peptide</td><td align="left" rowspan="1" colspan="1">6q</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">PQLC1</td><td align="left" rowspan="1" colspan="1">PQ loop repeat containing 1</td><td align="left" rowspan="1" colspan="1">18q</td><td align="left" rowspan="1" colspan="1">1.6(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.13</td></tr><tr><td align="left" rowspan="1" colspan="1">FLJ34048</td><td align="left" rowspan="1" colspan="1">hypothetical transcript</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1.6(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">CYP46A1</td><td align="left" rowspan="1" colspan="1">cytochrome P450, family 46, subfamily A, polypeptide 1</td><td align="left" rowspan="1" colspan="1">14q</td><td align="left" rowspan="1" colspan="1">2.0(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">ATG7</td><td align="left" rowspan="1" colspan="1">ATG7 autophagy related 7 homolog (S. cerevisiae)</td><td align="left" rowspan="1" colspan="1">3p</td><td align="left" rowspan="1" colspan="1">2.1(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">CXCL14</td><td align="left" rowspan="1" colspan="1">chemokine (C-X-C motif) ligand 14</td><td align="left" rowspan="1" colspan="1">5q</td><td align="left" rowspan="1" colspan="1">2.3(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.13</td></tr><tr><td align="left" rowspan="1" colspan="1">NXPH1</td><td align="left" rowspan="1" colspan="1">neurexophilin 1</td><td align="left" rowspan="1" colspan="1">7p</td><td align="left" rowspan="1" colspan="1">2.6(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">C17orf76</td><td align="left" rowspan="1" colspan="1">chromosome 17 open reading frame 76</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">2.7(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.17</td></tr><tr><td align="left" rowspan="1" colspan="1">NPM3</td><td align="left" rowspan="1" colspan="1">nucleophosmin/nucleoplasmin, 3</td><td align="left" rowspan="1" colspan="1">10q</td><td align="left" rowspan="1" colspan="1">3.1(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.08</td></tr><tr><td align="left" rowspan="1" colspan="1">LHX6</td><td align="left" rowspan="1" colspan="1">LIM homeobox 6</td><td align="left" rowspan="1" colspan="1">9q</td><td align="left" rowspan="1" colspan="1">3.2(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.07</td></tr><tr><td align="left" rowspan="1" colspan="1">FRMPD2</td><td align="left" rowspan="1" colspan="1">FERM and PDZ domain containing 2</td><td align="left" rowspan="1" colspan="1">10q</td><td align="left" rowspan="1" colspan="1">3.4(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.07</td></tr><tr><td align="left" rowspan="1" colspan="1">HSD11B1L</td><td align="left" rowspan="1" colspan="1">hydroxysteroid (11-beta) dehydrogenase 1-like</td><td align="left" rowspan="1" colspan="1">19p</td><td align="left" rowspan="1" colspan="1">3.5(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">SMAD3</td><td align="left" rowspan="1" colspan="1">SMAD family member 3</td><td align="left" rowspan="1" colspan="1">15q</td><td align="left" rowspan="1" colspan="1">3.5(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.09</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003024-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003024.t002</object-id><label>Table 2</label><caption><title>Genes whose expression increases with age in the human brain.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Symbol</td><td align="left" rowspan="1" colspan="1">Description</td><td align="left" rowspan="1" colspan="1">Chr.Pos</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic> value<xref ref-type="table-fn" rid="nt102">a</xref>
</td><td align="left" rowspan="1" colspan="1">R<sup>2 b</sup>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">ELF2</td><td align="left" rowspan="1" colspan="1">E74-like factor 2 (ets domain transcription factor)</td><td align="left" rowspan="1" colspan="1">4q</td><td align="left" rowspan="1" colspan="1">5.3(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">0.16</td></tr><tr><td align="left" rowspan="1" colspan="1">C1orf162</td><td align="left" rowspan="1" colspan="1">chromosome 1 open reading frame 162</td><td align="left" rowspan="1" colspan="1">1p</td><td align="left" rowspan="1" colspan="1">6.0(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">0.13</td></tr><tr><td align="left" rowspan="1" colspan="1">LRCH4</td><td align="left" rowspan="1" colspan="1">leucine-rich repeats and calponin homology (CH) domain containing 4</td><td align="left" rowspan="1" colspan="1">7q</td><td align="left" rowspan="1" colspan="1">2.9(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">0.09</td></tr><tr><td align="left" rowspan="1" colspan="1">MTUS1</td><td align="left" rowspan="1" colspan="1">mitochondrial tumor suppressor 1</td><td align="left" rowspan="1" colspan="1">8p</td><td align="left" rowspan="1" colspan="1">6.8(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">RUFY1</td><td align="left" rowspan="1" colspan="1">RUN and FYVE domain containing 1</td><td align="left" rowspan="1" colspan="1">5q</td><td align="left" rowspan="1" colspan="1">7.0(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">RDH5</td><td align="left" rowspan="1" colspan="1">retinol dehydrogenase 5 (11-cis/9-cis)</td><td align="left" rowspan="1" colspan="1">12q</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.06</td></tr><tr><td align="left" rowspan="1" colspan="1">TYK2</td><td align="left" rowspan="1" colspan="1">tyrosine kinase 2</td><td align="left" rowspan="1" colspan="1">19p</td><td align="left" rowspan="1" colspan="1">1.8(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">CLK1</td><td align="left" rowspan="1" colspan="1">CDC-like kinase 1</td><td align="left" rowspan="1" colspan="1">2q</td><td align="left" rowspan="1" colspan="1">2.1(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.14</td></tr><tr><td align="left" rowspan="1" colspan="1">TXNIP</td><td align="left" rowspan="1" colspan="1">thioredoxin interacting protein</td><td align="left" rowspan="1" colspan="1">1q</td><td align="left" rowspan="1" colspan="1">2.1(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">SLC16A9</td><td align="left" rowspan="1" colspan="1">solute carrier family 16 (monocarboxylic acid transporters), member 9</td><td align="left" rowspan="1" colspan="1">10q</td><td align="left" rowspan="1" colspan="1">4.8(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.06</td></tr><tr><td align="left" rowspan="1" colspan="1">ADORA3</td><td align="left" rowspan="1" colspan="1">adenosine A3 receptor</td><td align="left" rowspan="1" colspan="1">1p</td><td align="left" rowspan="1" colspan="1">5.1(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">UCKL1</td><td align="left" rowspan="1" colspan="1">uridine-cytidine kinase 1-like 1</td><td align="left" rowspan="1" colspan="1">20q</td><td align="left" rowspan="1" colspan="1">7.4(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">CTDSP2</td><td align="left" rowspan="1" colspan="1">nuclear LIM interactor-interacting factor 2</td><td align="left" rowspan="1" colspan="1">12q</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.05</td></tr><tr><td align="left" rowspan="1" colspan="1">HLA-DPB1</td><td align="left" rowspan="1" colspan="1">major histocompatibility complex, class II, DP beta 1</td><td align="left" rowspan="1" colspan="1">6p</td><td align="left" rowspan="1" colspan="1">1.2(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.09</td></tr><tr><td align="left" rowspan="1" colspan="1">PATL1</td><td align="left" rowspan="1" colspan="1">protein associated with topoisomerase II homolog 1 (yeast)</td><td align="left" rowspan="1" colspan="1">11q</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">GDPD3</td><td align="left" rowspan="1" colspan="1">glycerophosphodiester phosphodiesterase domain containing 3</td><td align="left" rowspan="1" colspan="1">16p</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.07</td></tr><tr><td align="left" rowspan="1" colspan="1">BHLHB3</td><td align="left" rowspan="1" colspan="1">basic helix-loop-helix domain containing, class B, 3</td><td align="left" rowspan="1" colspan="1">12p</td><td align="left" rowspan="1" colspan="1">1.7(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.07</td></tr><tr><td align="left" rowspan="1" colspan="1">RNASE4</td><td align="left" rowspan="1" colspan="1">ribonuclease, RNase A family, 4</td><td align="left" rowspan="1" colspan="1">14q</td><td align="left" rowspan="1" colspan="1">1.9(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.06</td></tr><tr><td align="left" rowspan="1" colspan="1">PLEKHM1</td><td align="left" rowspan="1" colspan="1">pleckstrin homology domain containing, family M (with RUN domain) member 1</td><td align="left" rowspan="1" colspan="1">17q</td><td align="left" rowspan="1" colspan="1">2.5(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.04</td></tr><tr><td align="left" rowspan="1" colspan="1">FAM46A</td><td align="left" rowspan="1" colspan="1">family with sequence similarity 46, member A</td><td align="left" rowspan="1" colspan="1">6q</td><td align="left" rowspan="1" colspan="1">3.4(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.06</td></tr><tr><td align="left" rowspan="1" colspan="1">CALCOCO1</td><td align="left" rowspan="1" colspan="1">calcium binding and coiled-coil domain 1</td><td align="left" rowspan="1" colspan="1">12q</td><td align="left" rowspan="1" colspan="1">3.5(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">0.09</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003024-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003024.t003</object-id><label>Table 3</label><caption><title>Terms in the Gene Ontology and KEGG pathway databases over-represented among genes that decreased expression with advancing age in the human brain (1,450 genes in total 13,216).</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Category</td><td align="left" rowspan="1" colspan="1">Term</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic> value<xref ref-type="table-fn" rid="nt103">a</xref>
</td><td align="left" rowspan="1" colspan="1">Bonferroni</td><td align="left" rowspan="1" colspan="1">Fold<sup>b</sup>
</td><td colspan="2" align="left" rowspan="1">Increasing group<sup>c</sup>
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Fold<sup>b</sup>
</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic> value<sup>d</sup>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">GO Bio.Process</td><td align="left" rowspan="1" colspan="1">nervous system development</td><td align="left" rowspan="1" colspan="1">7.8(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">3.8(10<sup>−5</sup>)</td><td align="left" rowspan="1" colspan="1">1.73</td><td align="left" rowspan="1" colspan="1">0.76</td><td align="left" rowspan="1" colspan="1">6.7(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">synaptic transmission</td><td align="left" rowspan="1" colspan="1">7.4(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">3.5(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">2.07</td><td align="left" rowspan="1" colspan="1">0.28</td><td align="left" rowspan="1" colspan="1">3.1(10<sup>−7</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">oxidative phosphorylation</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">6.7(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">2.73</td><td align="left" rowspan="1" colspan="1">0.15</td><td align="left" rowspan="1" colspan="1">5.3(10<sup>−5</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">GO Cell.Component</td><td align="left" rowspan="1" colspan="1">cytoplasmic part</td><td align="left" rowspan="1" colspan="1">9.0(10<sup>−11</sup>)</td><td align="left" rowspan="1" colspan="1">7.2(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">1.29</td><td align="left" rowspan="1" colspan="1">0.76</td><td align="left" rowspan="1" colspan="1">4.2(10<sup>−11</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">mitochondrion</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−10</sup>)</td><td align="left" rowspan="1" colspan="1">1.2(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">1.65</td><td align="left" rowspan="1" colspan="1">0.54</td><td align="left" rowspan="1" colspan="1">7.7(10<sup>−10</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">neuron projection</td><td align="left" rowspan="1" colspan="1">2.3(10<sup>−10</sup>)</td><td align="left" rowspan="1" colspan="1">1.8(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">3.41</td><td align="left" rowspan="1" colspan="1">0.40</td><td align="left" rowspan="1" colspan="1">1.0(10<sup>−2</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">synapse</td><td align="left" rowspan="1" colspan="1">1.8(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">2.55</td><td align="left" rowspan="1" colspan="1">0.42</td><td align="left" rowspan="1" colspan="1">6.4(10<sup>−4</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">mitochondrial membrane part</td><td align="left" rowspan="1" colspan="1">8.7(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">6.9(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">2.95</td><td align="left" rowspan="1" colspan="1">0.14</td><td align="left" rowspan="1" colspan="1">1.8(10<sup>−5</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">mitochondrial inner membrane</td><td align="left" rowspan="1" colspan="1">3.1(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">2.5(10<sup>−5</sup>)</td><td align="left" rowspan="1" colspan="1">2.11</td><td align="left" rowspan="1" colspan="1">0.45</td><td align="left" rowspan="1" colspan="1">5.4(10<sup>−5</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">mitochondrial respiratory chain</td><td align="left" rowspan="1" colspan="1">4.9(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">3.9(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">2.99</td><td align="left" rowspan="1" colspan="1">0.11</td><td align="left" rowspan="1" colspan="1">4.3(10<sup>−4</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">dendrite</td><td align="left" rowspan="1" colspan="1">5.2(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">4.2(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">3.60</td><td align="left" rowspan="1" colspan="1">0.17</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−2</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">GO Mol.Function</td><td align="left" rowspan="1" colspan="1">hydrogen ion transmembrane transporter activity</td><td align="left" rowspan="1" colspan="1">1.8(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">4.7(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">2.70</td><td align="left" rowspan="1" colspan="1">0.58</td><td align="left" rowspan="1" colspan="1">5.1(10<sup>−2</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">KEGG Pathway</td><td align="left" rowspan="1" colspan="1">Oxidative phosphorylation</td><td align="left" rowspan="1" colspan="1">1.9(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">3.8(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">2.66</td><td align="left" rowspan="1" colspan="1">0.25</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−4</sup>)</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003024-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003024.t004</object-id><label>Table 4</label><caption><title>Terms in the Gene Ontology and KEGG pathway databases over-represented among genes that increased expression with advancing age in the human brain (1,943 genes in total 13,216).</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Category</td><td align="left" rowspan="1" colspan="1">Term</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic> value<xref ref-type="table-fn" rid="nt104">a</xref>
</td><td align="left" rowspan="1" colspan="1">Bonferroni</td><td align="left" rowspan="1" colspan="1">Fold<sup>b</sup>
</td><td colspan="2" align="left" rowspan="1">Decreasing group<sup>c</sup>
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Fold<sup>b</sup>
</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic> value<sup>d</sup>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">GO Bio.Process</td><td align="left" rowspan="1" colspan="1">regulation of transcription, DNA-dependent</td><td align="left" rowspan="1" colspan="1">2.1(10<sup>−13</sup>)</td><td align="left" rowspan="1" colspan="1">1.0(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">1.43</td><td align="left" rowspan="1" colspan="1">0.56</td><td align="left" rowspan="1" colspan="1">6.1(10<sup>−13</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">chromosome organization and biogenesis</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">6.8(10<sup>−4</sup>)</td><td align="left" rowspan="1" colspan="1">1.80</td><td align="left" rowspan="1" colspan="1">0.53</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">DNA packaging</td><td align="left" rowspan="1" colspan="1">4.7(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">2.2(10<sup>−2</sup>)</td><td align="left" rowspan="1" colspan="1">1.78</td><td align="left" rowspan="1" colspan="1">0.49</td><td align="left" rowspan="1" colspan="1">1.9(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">DNA metabolic process</td><td align="left" rowspan="1" colspan="1">6.4(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">3.0(10<sup>−2</sup>)</td><td align="left" rowspan="1" colspan="1">1.45</td><td align="left" rowspan="1" colspan="1">0.59</td><td align="left" rowspan="1" colspan="1">7.6(10<sup>−5</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">GO Cel.Component</td><td align="left" rowspan="1" colspan="1">nucleus</td><td align="left" rowspan="1" colspan="1">1.2(10<sup>−13</sup>)</td><td align="left" rowspan="1" colspan="1">9.6(10<sup>−11</sup>)</td><td align="left" rowspan="1" colspan="1">1.26</td><td align="left" rowspan="1" colspan="1">0.75</td><td align="left" rowspan="1" colspan="1">3.7(10<sup>−11</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">GO Mol.Function</td><td align="left" rowspan="1" colspan="1">DNA binding</td><td align="left" rowspan="1" colspan="1">6.8(10<sup>−14</sup>)</td><td align="left" rowspan="1" colspan="1">1.7(10<sup>−10</sup>)</td><td align="left" rowspan="1" colspan="1">1.44</td><td align="left" rowspan="1" colspan="1">0.52</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−14</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">zinc ion binding</td><td align="left" rowspan="1" colspan="1">1.7(10<sup>−11</sup>)</td><td align="left" rowspan="1" colspan="1">4.3(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">1.39</td><td align="left" rowspan="1" colspan="1">0.64</td><td align="left" rowspan="1" colspan="1">6.3(10<sup>−9</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">transcription regulator activity</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">2.7(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">1.38</td><td align="left" rowspan="1" colspan="1">0.65</td><td align="left" rowspan="1" colspan="1">1.8(10<sup>−5</sup>)</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003024-t005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003024.t005</object-id><label>Table 5</label><caption><title>Genes whose expression decreases with age in human lymphocytes.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Symbol</td><td align="left" rowspan="1" colspan="1">Description</td><td align="left" rowspan="1" colspan="1">Chr.Pos</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic> value<xref ref-type="table-fn" rid="nt105">a</xref>
</td><td align="left" rowspan="1" colspan="1">R<sup>2</sup>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">LRRN3</td><td align="left" rowspan="1" colspan="1">leucine rich repeat neuronal 3</td><td align="left" rowspan="1" colspan="1">7q</td><td align="left" rowspan="1" colspan="1">5.1(10<sup>−94</sup>)</td><td align="left" rowspan="1" colspan="1">0.29</td></tr><tr><td align="left" rowspan="1" colspan="1">FCGBP</td><td align="left" rowspan="1" colspan="1">Fc fragment of IgG binding protein</td><td align="left" rowspan="1" colspan="1">19q</td><td align="left" rowspan="1" colspan="1">7.5(10<sup>−44</sup>)</td><td align="left" rowspan="1" colspan="1">0.14</td></tr><tr><td align="left" rowspan="1" colspan="1">FBLN2</td><td align="left" rowspan="1" colspan="1">fibulin 2</td><td align="left" rowspan="1" colspan="1">3p</td><td align="left" rowspan="1" colspan="1">8.1(10<sup>−44</sup>)</td><td align="left" rowspan="1" colspan="1">0.14</td></tr><tr><td align="left" rowspan="1" colspan="1">NRCAM</td><td align="left" rowspan="1" colspan="1">neuronal cell adhesion molecule</td><td align="left" rowspan="1" colspan="1">7q</td><td align="left" rowspan="1" colspan="1">3.3(10<sup>−43</sup>)</td><td align="left" rowspan="1" colspan="1">0.14</td></tr><tr><td align="left" rowspan="1" colspan="1">ITM2C</td><td align="left" rowspan="1" colspan="1">integral membrane protein 2C</td><td align="left" rowspan="1" colspan="1">2q</td><td align="left" rowspan="1" colspan="1">1.8(10<sup>−39</sup>)</td><td align="left" rowspan="1" colspan="1">0.13</td></tr><tr><td align="left" rowspan="1" colspan="1">PDE9A</td><td align="left" rowspan="1" colspan="1">phosphodiesterase 9A</td><td align="left" rowspan="1" colspan="1">21q</td><td align="left" rowspan="1" colspan="1">4.9(10<sup>−37</sup>)</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">ZNF154</td><td align="left" rowspan="1" colspan="1">zinc finger protein 154</td><td align="left" rowspan="1" colspan="1">19q</td><td align="left" rowspan="1" colspan="1">5.3(10<sup>−36</sup>)</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">ZSCAN18</td><td align="left" rowspan="1" colspan="1">zinc finger and SCAN domain containing 18</td><td align="left" rowspan="1" colspan="1">19q</td><td align="left" rowspan="1" colspan="1">4.7(10<sup>−35</sup>)</td><td align="left" rowspan="1" colspan="1">0.12</td></tr><tr><td align="left" rowspan="1" colspan="1">SATB1</td><td align="left" rowspan="1" colspan="1">SATB homeobox 1</td><td align="left" rowspan="1" colspan="1">3p</td><td align="left" rowspan="1" colspan="1">9.8(10<sup>−35</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">FLNB</td><td align="left" rowspan="1" colspan="1">filamin B, beta (actin binding protein 278)</td><td align="left" rowspan="1" colspan="1">3p</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−34</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">FAM134B</td><td align="left" rowspan="1" colspan="1">family with sequence similarity 134, member B</td><td align="left" rowspan="1" colspan="1">5p</td><td align="left" rowspan="1" colspan="1">3.1(10<sup>−34</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">SCD</td><td align="left" rowspan="1" colspan="1">stearoyl-CoA desaturase (delta-9-desaturase)</td><td align="left" rowspan="1" colspan="1">10q</td><td align="left" rowspan="1" colspan="1">8.8(10<sup>−34</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">SREBF1</td><td align="left" rowspan="1" colspan="1">sterol regulatory element binding transcription factor 1</td><td align="left" rowspan="1" colspan="1">17p</td><td align="left" rowspan="1" colspan="1">4.5(10<sup>−33</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">CCR7</td><td align="left" rowspan="1" colspan="1">chemokine (C-C motif) receptor 7</td><td align="left" rowspan="1" colspan="1">17q</td><td align="left" rowspan="1" colspan="1">6.3(10<sup>−33</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">PHGDH</td><td align="left" rowspan="1" colspan="1">phosphoglycerate dehydrogenase</td><td align="left" rowspan="1" colspan="1">1p</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−32</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">LEF1</td><td align="left" rowspan="1" colspan="1">lymphoid enhancer-binding factor 1</td><td align="left" rowspan="1" colspan="1">4q</td><td align="left" rowspan="1" colspan="1">5.3(10<sup>−32</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">NPM3</td><td align="left" rowspan="1" colspan="1">nucleophosmin/nucleoplasmin, 3</td><td align="left" rowspan="1" colspan="1">10q</td><td align="left" rowspan="1" colspan="1">9.8(10<sup>−32</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">OXNAD1</td><td align="left" rowspan="1" colspan="1">oxidoreductase NAD-binding domain containing 1</td><td align="left" rowspan="1" colspan="1">3p</td><td align="left" rowspan="1" colspan="1">2.1(10<sup>−31</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">TNNT3</td><td align="left" rowspan="1" colspan="1">troponin T type 3 (skeletal, fast)</td><td align="left" rowspan="1" colspan="1">11p</td><td align="left" rowspan="1" colspan="1">3.4(10<sup>−30</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">PLEKHG4</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">16q</td><td align="left" rowspan="1" colspan="1">5.8(10<sup>−30</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">MGC9913</td><td align="left" rowspan="1" colspan="1">hypothetical protein MGC9913</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1.5(10<sup>−29</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">SLC7A6</td><td align="left" rowspan="1" colspan="1">solute carrier family 7 (cationic amino acid transporter, y+ system), member 6</td><td align="left" rowspan="1" colspan="1">16q</td><td align="left" rowspan="1" colspan="1">4.4(10<sup>−29</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">CD27</td><td align="left" rowspan="1" colspan="1">CD27 molecule</td><td align="left" rowspan="1" colspan="1">12p</td><td align="left" rowspan="1" colspan="1">6.1(10<sup>−29</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">AEBP1</td><td align="left" rowspan="1" colspan="1">AE binding protein 1</td><td align="left" rowspan="1" colspan="1">7p</td><td align="left" rowspan="1" colspan="1">1.0(10<sup>−28</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">MGC29506</td><td align="left" rowspan="1" colspan="1">hypothetical protein MGC29506</td><td align="left" rowspan="1" colspan="1">5q</td><td align="left" rowspan="1" colspan="1">1.2(10<sup>−28</sup>)</td><td align="left" rowspan="1" colspan="1">0.09</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003024-t006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003024.t006</object-id><label>Table 6</label><caption><title>Genes whose expression increases with age in human lymphocytes.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Symbol</td><td align="left" rowspan="1" colspan="1">Description</td><td align="left" rowspan="1" colspan="1">Chr.Pos</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic> value<xref ref-type="table-fn" rid="nt106">a</xref>
</td><td align="left" rowspan="1" colspan="1">R<sup>2</sup>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">IGFBP3</td><td align="left" rowspan="1" colspan="1">insulin-like growth factor binding protein 3</td><td align="left" rowspan="1" colspan="1">7p</td><td align="left" rowspan="1" colspan="1">1.7(10<sup>−70</sup>)</td><td align="left" rowspan="1" colspan="1">0.22</td></tr><tr><td align="left" rowspan="1" colspan="1">SYT11</td><td align="left" rowspan="1" colspan="1">synaptotagmin XI</td><td align="left" rowspan="1" colspan="1">1q</td><td align="left" rowspan="1" colspan="1">7.4(10<sup>−60</sup>)</td><td align="left" rowspan="1" colspan="1">0.19</td></tr><tr><td align="left" rowspan="1" colspan="1">GZMH</td><td align="left" rowspan="1" colspan="1">granzyme H (cathepsin G-like 2, protein h-CCPX)</td><td align="left" rowspan="1" colspan="1">14q</td><td align="left" rowspan="1" colspan="1">2.0(10<sup>−48</sup>)</td><td align="left" rowspan="1" colspan="1">0.16</td></tr><tr><td align="left" rowspan="1" colspan="1">JAKMIP1</td><td align="left" rowspan="1" colspan="1">janus kinase and microtubule interacting protein 1</td><td align="left" rowspan="1" colspan="1">4p</td><td align="left" rowspan="1" colspan="1">5.4(10<sup>−45</sup>)</td><td align="left" rowspan="1" colspan="1">0.15</td></tr><tr><td align="left" rowspan="1" colspan="1">RCAN2</td><td align="left" rowspan="1" colspan="1">regulator of calcineurin 2</td><td align="left" rowspan="1" colspan="1">6p</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−43</sup>)</td><td align="left" rowspan="1" colspan="1">0.14</td></tr><tr><td align="left" rowspan="1" colspan="1">CRIP1</td><td align="left" rowspan="1" colspan="1">cysteine-rich protein 1 (intestinal)</td><td align="left" rowspan="1" colspan="1">14q</td><td align="left" rowspan="1" colspan="1">1.9(10<sup>−39</sup>)</td><td align="left" rowspan="1" colspan="1">0.13</td></tr><tr><td align="left" rowspan="1" colspan="1">PATL2</td><td align="left" rowspan="1" colspan="1">protein associated with topoisomerase II homolog 2 (yeast)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1.2(10<sup>−38</sup>)</td><td align="left" rowspan="1" colspan="1">0.13</td></tr><tr><td align="left" rowspan="1" colspan="1">MSC</td><td align="left" rowspan="1" colspan="1">musculin (activated B-cell factor-1)</td><td align="left" rowspan="1" colspan="1">8q</td><td align="left" rowspan="1" colspan="1">5.9(10<sup>−32</sup>)</td><td align="left" rowspan="1" colspan="1">0.11</td></tr><tr><td align="left" rowspan="1" colspan="1">GDPD5</td><td align="left" rowspan="1" colspan="1">glycerophosphodiester phosphodiesterase domain containing 5</td><td align="left" rowspan="1" colspan="1">11q</td><td align="left" rowspan="1" colspan="1">3.1(10<sup>−30</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">APOBEC3H</td><td align="left" rowspan="1" colspan="1">apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3H</td><td align="left" rowspan="1" colspan="1">22q</td><td align="left" rowspan="1" colspan="1">2.3(10<sup>−29</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">CCL5</td><td align="left" rowspan="1" colspan="1">chemokine (C-C motif) ligand 5</td><td align="left" rowspan="1" colspan="1">17q</td><td align="left" rowspan="1" colspan="1">7.1(10<sup>−29</sup>)</td><td align="left" rowspan="1" colspan="1">0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">GFI1</td><td align="left" rowspan="1" colspan="1">growth factor independent 1 transcription repressor</td><td align="left" rowspan="1" colspan="1">1p</td><td align="left" rowspan="1" colspan="1">2.4(10<sup>−28</sup>)</td><td align="left" rowspan="1" colspan="1">0.09</td></tr><tr><td align="left" rowspan="1" colspan="1">MANEAL</td><td align="left" rowspan="1" colspan="1">mannosidase, endo-alpha-like</td><td align="left" rowspan="1" colspan="1">1p</td><td align="left" rowspan="1" colspan="1">2.9(10<sup>−28</sup>)</td><td align="left" rowspan="1" colspan="1">0.09</td></tr><tr><td align="left" rowspan="1" colspan="1">KIF21A</td><td align="left" rowspan="1" colspan="1">kinesin family member 21A</td><td align="left" rowspan="1" colspan="1">12q</td><td align="left" rowspan="1" colspan="1">3.1(10<sup>−27</sup>)</td><td align="left" rowspan="1" colspan="1">0.09</td></tr><tr><td align="left" rowspan="1" colspan="1">GPR137B</td><td align="left" rowspan="1" colspan="1">G protein-coupled receptor 137B</td><td align="left" rowspan="1" colspan="1">1q</td><td align="left" rowspan="1" colspan="1">5.4(10<sup>−27</sup>)</td><td align="left" rowspan="1" colspan="1">0.09</td></tr><tr><td align="left" rowspan="1" colspan="1">PDGFRB</td><td align="left" rowspan="1" colspan="1">platelet-derived growth factor receptor, beta polypeptide</td><td align="left" rowspan="1" colspan="1">5q</td><td align="left" rowspan="1" colspan="1">3.7(10<sup>−26</sup>)</td><td align="left" rowspan="1" colspan="1">0.09</td></tr><tr><td align="left" rowspan="1" colspan="1">PCBP4</td><td align="left" rowspan="1" colspan="1">poly(rC) binding protein 4</td><td align="left" rowspan="1" colspan="1">3p</td><td align="left" rowspan="1" colspan="1">8.1(10<sup>−26</sup>)</td><td align="left" rowspan="1" colspan="1">0.09</td></tr><tr><td align="left" rowspan="1" colspan="1">B3GAT1</td><td align="left" rowspan="1" colspan="1">beta-1,3-glucuronyltransferase 1 (glucuronosyltransferase P)</td><td align="left" rowspan="1" colspan="1">11q</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−25</sup>)</td><td align="left" rowspan="1" colspan="1">0.08</td></tr><tr><td align="left" rowspan="1" colspan="1">LLGL2</td><td align="left" rowspan="1" colspan="1">lethal giant larvae homolog 2 (Drosophila)</td><td align="left" rowspan="1" colspan="1">17q</td><td align="left" rowspan="1" colspan="1">4.7(10<sup>−25</sup>)</td><td align="left" rowspan="1" colspan="1">0.08</td></tr><tr><td align="left" rowspan="1" colspan="1">LAG3</td><td align="left" rowspan="1" colspan="1">lymphocyte-activation gene 3</td><td align="left" rowspan="1" colspan="1">12p</td><td align="left" rowspan="1" colspan="1">5.1(10<sup>−25</sup>)</td><td align="left" rowspan="1" colspan="1">0.08</td></tr><tr><td align="left" rowspan="1" colspan="1">PPP2R2B</td><td align="left" rowspan="1" colspan="1">beta isoform of regulatory subunit B55, protein phosphatase 2 isoform b</td><td align="left" rowspan="1" colspan="1">5q</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−24</sup>)</td><td align="left" rowspan="1" colspan="1">0.08</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">hypothetical gene supported by BC040060</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">9.1(10<sup>−24</sup>)</td><td align="left" rowspan="1" colspan="1">0.08</td></tr><tr><td align="left" rowspan="1" colspan="1">PRSS23</td><td align="left" rowspan="1" colspan="1">protease, serine, 23</td><td align="left" rowspan="1" colspan="1">11q</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−23</sup>)</td><td align="left" rowspan="1" colspan="1">0.08</td></tr><tr><td align="left" rowspan="1" colspan="1">B4GALT5</td><td align="left" rowspan="1" colspan="1">UDP-Gal:betaGlcNAc beta 1,4- galactosyltransferase, polypeptide 5</td><td align="left" rowspan="1" colspan="1">20q</td><td align="left" rowspan="1" colspan="1">1.9(10<sup>−22</sup>)</td><td align="left" rowspan="1" colspan="1">0.07</td></tr><tr><td align="left" rowspan="1" colspan="1">MXRA7</td><td align="left" rowspan="1" colspan="1">matrix-remodelling associated 7</td><td align="left" rowspan="1" colspan="1">17q</td><td align="left" rowspan="1" colspan="1">3.0(10<sup>−22</sup>)</td><td align="left" rowspan="1" colspan="1">0.07</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003024-t007" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003024.t007</object-id><label>Table 7</label><caption><title>Terms in the Gene Ontology and KEGG pathway databases over-represented among genes that decreased expression with advancing age in human lymphocytes (1,878 genes in total 13,232).</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Category</td><td align="left" rowspan="1" colspan="1">Term</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic> value<xref ref-type="table-fn" rid="nt107">a</xref>
</td><td align="left" rowspan="1" colspan="1">Bonferroni</td><td align="left" rowspan="1" colspan="1">Fold<sup>b</sup>
</td><td colspan="2" align="left" rowspan="1">Increasing group<sup>c</sup>
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Fold<sup>b</sup>
</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic> value<sup>d</sup>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">GO Bio.Process</td><td align="left" rowspan="1" colspan="1">translation</td><td align="left" rowspan="1" colspan="1">4.2(10<sup>−21</sup>)</td><td align="left" rowspan="1" colspan="1">2.0(10<sup>−17</sup>)</td><td align="left" rowspan="1" colspan="1">2.19</td><td align="left" rowspan="1" colspan="1">0.60</td><td align="left" rowspan="1" colspan="1">8.7(10<sup>−4</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">cellular biosynthetic process</td><td align="left" rowspan="1" colspan="1">6.7(10<sup>−21</sup>)</td><td align="left" rowspan="1" colspan="1">3.2(10<sup>−17</sup>)</td><td align="left" rowspan="1" colspan="1">1.83</td><td align="left" rowspan="1" colspan="1">0.79</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−2</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">gene expression</td><td align="left" rowspan="1" colspan="1">4.0(10<sup>−18</sup>)</td><td align="left" rowspan="1" colspan="1">1.9(10<sup>−14</sup>)</td><td align="left" rowspan="1" colspan="1">1.36</td><td align="left" rowspan="1" colspan="1">0.70</td><td align="left" rowspan="1" colspan="1">3.1(10<sup>−11</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">ribosome biogenesis and assembly</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">7.2(10<sup>−5</sup>)</td><td align="left" rowspan="1" colspan="1">2.90</td><td align="left" rowspan="1" colspan="1">0.23</td><td align="left" rowspan="1" colspan="1">5.9(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">tRNA metabolic process</td><td align="left" rowspan="1" colspan="1">2.9(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−4</sup>)</td><td align="left" rowspan="1" colspan="1">2.57</td><td align="left" rowspan="1" colspan="1">0.52</td><td align="left" rowspan="1" colspan="1">4.7(10<sup>−2</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">RNA metabolic process</td><td align="left" rowspan="1" colspan="1">2.2(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">1.0(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">1.24</td><td align="left" rowspan="1" colspan="1">0.71</td><td align="left" rowspan="1" colspan="1">6.9(10<sup>−9</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">RNA processing</td><td align="left" rowspan="1" colspan="1">8.9(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">4.2(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">1.61</td><td align="left" rowspan="1" colspan="1">0.23</td><td align="left" rowspan="1" colspan="1">2.5(10<sup>−11</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">tRNA processing</td><td align="left" rowspan="1" colspan="1">4.3(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">2.0(10<sup>−2</sup>)</td><td align="left" rowspan="1" colspan="1">2.82</td><td align="left" rowspan="1" colspan="1">0.63</td><td align="left" rowspan="1" colspan="1">2.2(10<sup>−1</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">rRNA processing</td><td align="left" rowspan="1" colspan="1">8.4(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">3.9(10<sup>−2</sup>)</td><td align="left" rowspan="1" colspan="1">2.80</td><td align="left" rowspan="1" colspan="1">0.16</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−2</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">GO Cell.Component</td><td align="left" rowspan="1" colspan="1">ribosome</td><td align="left" rowspan="1" colspan="1">2.0(10<sup>−23</sup>)</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−20</sup>)</td><td align="left" rowspan="1" colspan="1">3.02</td><td align="left" rowspan="1" colspan="1">0.09</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−8</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">ribosomal subunit</td><td align="left" rowspan="1" colspan="1">4.2(10<sup>−21</sup>)</td><td align="left" rowspan="1" colspan="1">3.3(10<sup>−18</sup>)</td><td align="left" rowspan="1" colspan="1">3.59</td><td align="left" rowspan="1" colspan="1">0.08</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−5</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">ribonucleoprotein complex</td><td align="left" rowspan="1" colspan="1">1.6(10<sup>−20</sup>)</td><td align="left" rowspan="1" colspan="1">1.3(10<sup>−17</sup>)</td><td align="left" rowspan="1" colspan="1">2.19</td><td align="left" rowspan="1" colspan="1">0.22</td><td align="left" rowspan="1" colspan="1">3.5(10<sup>−12</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">organelle lumen</td><td align="left" rowspan="1" colspan="1">9.0(10<sup>−16</sup>)</td><td align="left" rowspan="1" colspan="1">7.0(10<sup>−13</sup>)</td><td align="left" rowspan="1" colspan="1">1.71</td><td align="left" rowspan="1" colspan="1">0.55</td><td align="left" rowspan="1" colspan="1">2.2(10<sup>−7</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">mitochondrion</td><td align="left" rowspan="1" colspan="1">5.0(10<sup>−12</sup>)</td><td align="left" rowspan="1" colspan="1">3.9(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">1.62</td><td align="left" rowspan="1" colspan="1">0.60</td><td align="left" rowspan="1" colspan="1">7.4(10<sup>−6</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">small ribosomal subunit</td><td align="left" rowspan="1" colspan="1">2.7(10<sup>−11</sup>)</td><td align="left" rowspan="1" colspan="1">2.2(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">3.72</td><td align="left" rowspan="1" colspan="1">0.00</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">large ribosomal subunit</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−10</sup>)</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">3.45</td><td align="left" rowspan="1" colspan="1">0.14</td><td align="left" rowspan="1" colspan="1">4.6(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">cytosolic part</td><td align="left" rowspan="1" colspan="1">3.1(10<sup>−10</sup>)</td><td align="left" rowspan="1" colspan="1">2.5(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">2.60</td><td align="left" rowspan="1" colspan="1">0.63</td><td align="left" rowspan="1" colspan="1">8.0(10<sup>−2</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">nucleolus</td><td align="left" rowspan="1" colspan="1">5.6(10<sup>−10</sup>)</td><td align="left" rowspan="1" colspan="1">4.4(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">2.42</td><td align="left" rowspan="1" colspan="1">0.47</td><td align="left" rowspan="1" colspan="1">8.7(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">intracellular organelle part</td><td align="left" rowspan="1" colspan="1">4.9(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">3.8(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">1.24</td><td align="left" rowspan="1" colspan="1">0.83</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−4</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">mitochondrial part</td><td align="left" rowspan="1" colspan="1">6.4(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">5.0(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">1.72</td><td align="left" rowspan="1" colspan="1">0.63</td><td align="left" rowspan="1" colspan="1">1.9(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">mitochondrial matrix</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">1.2(10<sup>−5</sup>)</td><td align="left" rowspan="1" colspan="1">2.25</td><td align="left" rowspan="1" colspan="1">0.39</td><td align="left" rowspan="1" colspan="1">1.8(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">mitochondrial small ribosomal subunit</td><td align="left" rowspan="1" colspan="1">4.5(10<sup>−5</sup>)</td><td align="left" rowspan="1" colspan="1">3.4(10<sup>−2</sup>)</td><td align="left" rowspan="1" colspan="1">4.25</td><td align="left" rowspan="1" colspan="1">0.00</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−1</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">GO Mol.Function</td><td align="left" rowspan="1" colspan="1">structural constituent of ribosome</td><td align="left" rowspan="1" colspan="1">9.4(10<sup>−22</sup>)</td><td align="left" rowspan="1" colspan="1">2.4(10<sup>−18</sup>)</td><td align="left" rowspan="1" colspan="1">3.08</td><td align="left" rowspan="1" colspan="1">0.10</td><td align="left" rowspan="1" colspan="1">1.2(10<sup>−7</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">nucleic acid binding</td><td align="left" rowspan="1" colspan="1">2.4(10<sup>−19</sup>)</td><td align="left" rowspan="1" colspan="1">6.0(10<sup>−16</sup>)</td><td align="left" rowspan="1" colspan="1">1.40</td><td align="left" rowspan="1" colspan="1">0.64</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−14</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">RNA binding</td><td align="left" rowspan="1" colspan="1">2.6(10<sup>−19</sup>)</td><td align="left" rowspan="1" colspan="1">6.7(10<sup>−16</sup>)</td><td align="left" rowspan="1" colspan="1">1.95</td><td align="left" rowspan="1" colspan="1">0.34</td><td align="left" rowspan="1" colspan="1">1.3(10<sup>−11</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">methyltransferase activity</td><td align="left" rowspan="1" colspan="1">3.0(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">7.6(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">2.11</td><td align="left" rowspan="1" colspan="1">0.33</td><td align="left" rowspan="1" colspan="1">1.3(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">oxidoreductase activity (NAD or NADP)</td><td align="left" rowspan="1" colspan="1">5.2(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">1.3(10<sup>−2</sup>)</td><td align="left" rowspan="1" colspan="1">2.60</td><td align="left" rowspan="1" colspan="1">0.76</td><td align="left" rowspan="1" colspan="1">3.1(10<sup>−1</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">KEGG Pathway</td><td align="left" rowspan="1" colspan="1">Ribosome</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−17</sup>)</td><td align="left" rowspan="1" colspan="1">2.2(10<sup>−15</sup>)</td><td align="left" rowspan="1" colspan="1">3.97</td><td align="left" rowspan="1" colspan="1">0.00</td><td align="left" rowspan="1" colspan="1">3.4(10<sup>−5</sup>)</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003024-t008" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003024.t008</object-id><label>Table 8</label><caption><title>Terms in the Gene Ontology and KEGG pathway databases over-represented among genes that increased expression with advancing age in human lymphocytes (1,430 genes in total 13,232).</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Category</td><td align="left" rowspan="1" colspan="1">Term</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic> value<xref ref-type="table-fn" rid="nt108">a</xref>
</td><td align="left" rowspan="1" colspan="1">Bonferroni</td><td align="left" rowspan="1" colspan="1">Fold<sup>b</sup>
</td><td colspan="2" align="left" rowspan="1">Decreasing group<sup>d</sup>
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Fold<sup>b</sup>
</td><td align="left" rowspan="1" colspan="1">
<italic>P</italic> value<sup>e</sup>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">GO Bio.Process</td><td align="left" rowspan="1" colspan="1">signal transduction</td><td align="left" rowspan="1" colspan="1">2.1(10<sup>−10</sup>)</td><td align="left" rowspan="1" colspan="1">1.0(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">1.34</td><td align="left" rowspan="1" colspan="1">0.72</td><td align="left" rowspan="1" colspan="1">8.8(10<sup>−11</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">immune response</td><td align="left" rowspan="1" colspan="1">9.2(10<sup>−10</sup>)</td><td align="left" rowspan="1" colspan="1">4.3(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">1.84</td><td align="left" rowspan="1" colspan="1">1.02</td><td align="left" rowspan="1" colspan="1">6.0(10<sup>−1</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">defense response</td><td align="left" rowspan="1" colspan="1">7.0(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">3.3(10<sup>−5</sup>)</td><td align="left" rowspan="1" colspan="1">1.90</td><td align="left" rowspan="1" colspan="1">0.76</td><td align="left" rowspan="1" colspan="1">2.7(10<sup>−2</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">response to external stimulus</td><td align="left" rowspan="1" colspan="1">1.3(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">6.2(10<sup>−5</sup>)</td><td align="left" rowspan="1" colspan="1">1.84</td><td align="left" rowspan="1" colspan="1">0.66</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">cytoskeleton organization and biogenesis</td><td align="left" rowspan="1" colspan="1">6.8(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">3.2(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">1.78</td><td align="left" rowspan="1" colspan="1">0.58</td><td align="left" rowspan="1" colspan="1">2.4(10<sup>−4</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">positive regulation of cellular process</td><td align="left" rowspan="1" colspan="1">9.6(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">4.5(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">1.51</td><td align="left" rowspan="1" colspan="1">0.97</td><td align="left" rowspan="1" colspan="1">4.1(10<sup>−1</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">response to wounding</td><td align="left" rowspan="1" colspan="1">7.2(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">3.3(10<sup>−2</sup>)</td><td align="left" rowspan="1" colspan="1">1.80</td><td align="left" rowspan="1" colspan="1">0.62</td><td align="left" rowspan="1" colspan="1">3.0(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">cell adhesion</td><td align="left" rowspan="1" colspan="1">7.6(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">3.5(10<sup>−2</sup>)</td><td align="left" rowspan="1" colspan="1">1.64</td><td align="left" rowspan="1" colspan="1">0.71</td><td align="left" rowspan="1" colspan="1">6.0(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">GO Cel.Component</td><td align="left" rowspan="1" colspan="1">plasma membrane</td><td align="left" rowspan="1" colspan="1">7.0(10<sup>−11</sup>)</td><td align="left" rowspan="1" colspan="1">5.5(10<sup>−8</sup>)</td><td align="left" rowspan="1" colspan="1">1.42</td><td align="left" rowspan="1" colspan="1">0.62</td><td align="left" rowspan="1" colspan="1">1.0(10<sup>−13</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">membrane part</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">8.9(10<sup>−5</sup>)</td><td align="left" rowspan="1" colspan="1">1.21</td><td align="left" rowspan="1" colspan="1">0.79</td><td align="left" rowspan="1" colspan="1">3.5(10<sup>−10</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">cytoskeleton</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−4</sup>)</td><td align="left" rowspan="1" colspan="1">1.56</td><td align="left" rowspan="1" colspan="1">0.69</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−4</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">integral to plasma membrane</td><td align="left" rowspan="1" colspan="1">2.3(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">1.8(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">1.49</td><td align="left" rowspan="1" colspan="1">0.69</td><td align="left" rowspan="1" colspan="1">1.0(10<sup>−4</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">GO Mol.Function</td><td align="left" rowspan="1" colspan="1">protein binding</td><td align="left" rowspan="1" colspan="1">9.6(10<sup>−10</sup>)</td><td align="left" rowspan="1" colspan="1">2.4(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">1.16</td><td align="left" rowspan="1" colspan="1">0.93</td><td align="left" rowspan="1" colspan="1">3.3(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">actin binding</td><td align="left" rowspan="1" colspan="1">9.1(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">2.3(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">2.07</td><td align="left" rowspan="1" colspan="1">0.57</td><td align="left" rowspan="1" colspan="1">4.4(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">signal transducer activity</td><td align="left" rowspan="1" colspan="1">1.3(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">3.2(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">1.37</td><td align="left" rowspan="1" colspan="1">0.79</td><td align="left" rowspan="1" colspan="1">4.8(10<sup>−4</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">GTPase activity</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">3.5(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">2.24</td><td align="left" rowspan="1" colspan="1">0.62</td><td align="left" rowspan="1" colspan="1">2.7(10<sup>−2</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">receptor binding</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−6</sup>)</td><td align="left" rowspan="1" colspan="1">3.7(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">1.68</td><td align="left" rowspan="1" colspan="1">0.76</td><td align="left" rowspan="1" colspan="1">1.8(10<sup>−2</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">GTP binding</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−5</sup>)</td><td align="left" rowspan="1" colspan="1">3.9(10<sup>−2</sup>)</td><td align="left" rowspan="1" colspan="1">1.79</td><td align="left" rowspan="1" colspan="1">0.86</td><td align="left" rowspan="1" colspan="1">2.0(10<sup>−1</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">KEGG Pathway</td><td align="left" rowspan="1" colspan="1">Regulation of actin cytoskeleton</td><td align="left" rowspan="1" colspan="1">3.1(10<sup>−9</sup>)</td><td align="left" rowspan="1" colspan="1">6.0(10<sup>−7</sup>)</td><td align="left" rowspan="1" colspan="1">2.44</td><td align="left" rowspan="1" colspan="1">0.52</td><td align="left" rowspan="1" colspan="1">6.4(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Natural killer cell mediated cytotoxicity</td><td align="left" rowspan="1" colspan="1">1.3(10<sup>−5</sup>)</td><td align="left" rowspan="1" colspan="1">2.6(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">2.25</td><td align="left" rowspan="1" colspan="1">0.26</td><td align="left" rowspan="1" colspan="1">2.3(10<sup>−4</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Focal adhesion</td><td align="left" rowspan="1" colspan="1">3.2(10<sup>−5</sup>)</td><td align="left" rowspan="1" colspan="1">6.4(10<sup>−3</sup>)</td><td align="left" rowspan="1" colspan="1">2.02</td><td align="left" rowspan="1" colspan="1">0.55</td><td align="left" rowspan="1" colspan="1">1.2(10<sup>−2</sup>)</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003024-t009" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003024.t009</object-id><label>Table 9</label><caption><title>Genomic architecture comparisons for age associated genes in human brain samples.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Negative asso.<sup>a</sup>
</td><td align="left" rowspan="1" colspan="1">Positive asso.<sup>b</sup>
</td><td align="left" rowspan="1" colspan="1">Unregulated</td><td align="left" rowspan="1" colspan="1">P-value</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">(n = 1,369)</td><td align="left" rowspan="1" colspan="1">(n = 1,746)</td><td align="left" rowspan="1" colspan="1">(n = 9,405)</td><td align="left" rowspan="1" colspan="1"/></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">pre-mRNA length</td><td align="left" rowspan="1" colspan="1">75652±4560</td><td align="left" rowspan="1" colspan="1">53683±2018</td><td align="left" rowspan="1" colspan="1">62344±1186</td><td align="left" rowspan="1" colspan="1">1.9(10<sup>−1</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">Coding sequence length</td><td align="left" rowspan="1" colspan="1">1311±28</td><td align="left" rowspan="1" colspan="1">2057±44</td><td align="left" rowspan="1" colspan="1">1722±18</td><td align="left" rowspan="1" colspan="1">7.0(10<sup>−60</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">Number of exons</td><td align="left" rowspan="1" colspan="1">9.2±0.2</td><td align="left" rowspan="1" colspan="1">13.0±0.3</td><td align="left" rowspan="1" colspan="1">11.2±0.1</td><td align="left" rowspan="1" colspan="1">1.1(10<sup>−22</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">Total intron length</td><td align="left" rowspan="1" colspan="1">74662±4636</td><td align="left" rowspan="1" colspan="1">51237±2033</td><td align="left" rowspan="1" colspan="1">60622±1202</td><td align="left" rowspan="1" colspan="1">1.4(10<sup>−1</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">5′ UTR length</td><td align="left" rowspan="1" colspan="1">215±6</td><td align="left" rowspan="1" colspan="1">227±6</td><td align="left" rowspan="1" colspan="1">210±2</td><td align="left" rowspan="1" colspan="1">3.0(10<sup>−2</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">3′ UTR length</td><td align="left" rowspan="1" colspan="1">1155±34</td><td align="left" rowspan="1" colspan="1">1219±29</td><td align="left" rowspan="1" colspan="1">1202±13</td><td align="left" rowspan="1" colspan="1">3.6(10<sup>−1</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">Average intron of ea. transcript</td><td align="left" rowspan="1" colspan="1">8651±485</td><td align="left" rowspan="1" colspan="1">5312±249</td><td align="left" rowspan="1" colspan="1">6648±144</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−9</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">Average exon of ea. transcript</td><td align="left" rowspan="1" colspan="1">424±14</td><td align="left" rowspan="1" colspan="1">435±13</td><td align="left" rowspan="1" colspan="1">429±5</td><td align="left" rowspan="1" colspan="1">9.9(10<sup>−1</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">Average intron / average exon</td><td align="left" rowspan="1" colspan="1">28.7±1.5</td><td align="left" rowspan="1" colspan="1">16.4±0.7</td><td align="left" rowspan="1" colspan="1">21.9±0.5</td><td align="left" rowspan="1" colspan="1">5.7(10<sup>−10</sup>)</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003024-t010" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003024.t010</object-id><label>Table 10</label><caption><title>Genomic architecture comparisons for age associated genes in human lymphocyte samples.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">Negative asso.<sup>a</sup>
</td><td align="left" rowspan="1" colspan="1">Positive asso.<sup>b</sup>
</td><td align="left" rowspan="1" colspan="1">Unregulated</td><td align="left" rowspan="1" colspan="1">P-value</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">(n = 2,364)</td><td align="left" rowspan="1" colspan="1">(n = 1,980)</td><td align="left" rowspan="1" colspan="1">(n = 7,898)</td><td align="left" rowspan="1" colspan="1"/></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">pre-mRNA length</td><td align="left" rowspan="1" colspan="1">51908±1755</td><td align="left" rowspan="1" colspan="1">49489±1880</td><td align="left" rowspan="1" colspan="1">54879±1115</td><td align="left" rowspan="1" colspan="1">1.2(10<sup>−6</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">Coding sequence length</td><td align="left" rowspan="1" colspan="1">1724±52</td><td align="left" rowspan="1" colspan="1">1628±51</td><td align="left" rowspan="1" colspan="1">1710±18</td><td align="left" rowspan="1" colspan="1">6.9(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">Number of exons</td><td align="left" rowspan="1" colspan="1">11.5±0.2</td><td align="left" rowspan="1" colspan="1">10.7±0.2</td><td align="left" rowspan="1" colspan="1">11.1±0.1</td><td align="left" rowspan="1" colspan="1">2.1(10<sup>−5</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">Total intron length</td><td align="left" rowspan="1" colspan="1">49705±1764</td><td align="left" rowspan="1" colspan="1">47549±1893</td><td align="left" rowspan="1" colspan="1">53364±1135</td><td align="left" rowspan="1" colspan="1">7.5(10<sup>−9</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">5′ UTR length</td><td align="left" rowspan="1" colspan="1">195±4</td><td align="left" rowspan="1" colspan="1">211±5</td><td align="left" rowspan="1" colspan="1">208±2</td><td align="left" rowspan="1" colspan="1">1.0(10<sup>−3</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">3′ UTR length</td><td align="left" rowspan="1" colspan="1">1104±25</td><td align="left" rowspan="1" colspan="1">1134±28</td><td align="left" rowspan="1" colspan="1">1185±15</td><td align="left" rowspan="1" colspan="1">3.9(10<sup>−4</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">Average intron of ea. transcript</td><td align="left" rowspan="1" colspan="1">5497±209</td><td align="left" rowspan="1" colspan="1">5463±256</td><td align="left" rowspan="1" colspan="1">5651±110</td><td align="left" rowspan="1" colspan="1">1.5(10<sup>−6</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">Average exon of ea. transcript</td><td align="left" rowspan="1" colspan="1">381±10</td><td align="left" rowspan="1" colspan="1">420±11</td><td align="left" rowspan="1" colspan="1">435±6</td><td align="left" rowspan="1" colspan="1">7.7(10<sup>−11</sup>)</td></tr><tr><td align="left" rowspan="1" colspan="1">Average intron / average exon</td><td align="left" rowspan="1" colspan="1">20.5±1.0</td><td align="left" rowspan="1" colspan="1">17.7±0.7</td><td align="left" rowspan="1" colspan="1">19.2±0.4</td><td align="left" rowspan="1" colspan="1">7.1(10<sup>−10</sup>)</td></tr></tbody></table></alternatives></table-wrap>
<disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula> <disp-formula></disp-formula>
<supplementary-material content-type="local-data" id="pone.0003024.s001"><label>Figure S1</label><caption><p>Overrepresented GO term analysis</p><p>(0.06 MB PDF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0003024.s002"><label>Table S1</label><caption><p>Terms in the Gene Ontology and KEGG pathway databases enriched among genes that decreased expression with advancing age in brain (1450 genes in total 13,216)</p><p>(0.15 MB PDF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0003024.s003"><label>Table S2</label><caption><p>Terms in the Gene Ontology and KEGG pathway databases enriched among genes that increased expression with advancing age in brain (1943 genes in total 13,216)</p><p>(0.09 MB PDF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0003024.s004"><label>Table S3</label><caption><p>Terms in the Gene Ontology and KEGG pathway databases enriched among genes that decreased expression with advancing age in the human lymphocytes (1878 genes in total 13,232)</p><p>(0.13 MB PDF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0003024.s005"><label>Table S4</label><caption><p>Terms in the Gene Ontology and KEGG pathway databases enriched among genes that increased expression with advancing age in the human lymphocytes (1430 genes in total 13,232)</p><p>(0.12 MB PDF)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label>a</label><p>
<italic>P</italic>-value with the null hypothesis <italic>β</italic>
<sub>1</sub> = 0 in linear regression model <italic>Y</italic> = <italic>β</italic>
<sub>0</sub>+<italic>β</italic>
<sub>1</sub>
<italic>x</italic>
<sub>1</sub>+<italic>β</italic>
<sub>2</sub>
<italic>x</italic>
<sub>2</sub>+<italic>ε</italic>, where <italic>Y</italic> is transcript expression, <italic>x</italic>
<sub>1</sub> age-at-death, <italic>x</italic>
<sub>2</sub> global expression, and <italic>ε</italic> random error. <sup>b</sup>Coefficient of determination of linear regression model <italic>Y</italic> = <italic>β</italic>
<sub>0</sub>+<italic>β</italic>
<sub>1</sub>
<italic>x</italic>
<sub>1</sub>+<italic>ε</italic>, where each variable is as described above.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt102"><label>a</label><p>
<italic>P</italic>-value with the null hypothesis <italic>β</italic>
<sub>1</sub> = 0 in linear regression model <italic>Y</italic> = <italic>β</italic>
<sub>0</sub>+<italic>β</italic>
<sub>1</sub>
<italic>x</italic>
<sub>1</sub>+<italic>β</italic>
<sub>2</sub>
<italic>x</italic>
<sub>2</sub>+<italic>ε</italic>, where <italic>Y</italic> is transcript expression, <italic>x</italic>
<sub>1</sub> age-at-death, <italic>x</italic>
<sub>2</sub> global expression, and <italic>ε</italic> random error. <sup>b</sup>Coefficient of determination of linear regression model <italic>Y</italic> = <italic>β</italic>
<sub>0</sub>+<italic>β</italic>
<sub>1</sub>
<italic>x</italic>
<sub>1</sub>+<italic>ε</italic>, where each variable is as described above.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt103"><label>a</label><p>EASE-score, <italic>P-</italic>value of a modified Fisher's exact test for overrepresentation ##REF##14519205##[35]##. <sup>b</sup>Fold enrichment in each gene group compared to the base set. <sup>c</sup>Gene group with expression that increases with advancing age-at-death (same genes used in ##TAB##3##Table 4##, 1943 genes). <sup>d</sup>Hypergeometric test for underrepresentation using annotated genes from a total set of 13,216 genes as base population ##REF##12620386##[36]##.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt104"><label>a</label><p>EASE-score, <italic>P-</italic>value of a modified Fisher's exact test for overrepresentation ##REF##14519205##[35]##. <sup>b</sup>Fold enrichment in each gene group compared to the base set. <sup>c</sup>Gene group with expression that decreases with advancing age-at-death (same genes used in ##TAB##2##Table 3##, 1450 genes). <sup>d</sup>Hypergeometric test for underrepresentation using annotated genes from a total set of 13,216 genes as base population ##REF##12620386##[36]##.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt105"><label>a</label><p>
<italic>P</italic>-value with the null hypothesis <italic>β</italic>
<sub>1</sub> = 0 in linear regression model <italic>Y</italic> = <italic>β</italic>
<sub>0</sub>+<italic>β</italic>
<sub>1</sub>
<italic>x</italic>
<sub>1</sub>+<italic>ε</italic>, where <italic>Y</italic> is transcript expression, <italic>x</italic>
<sub>1</sub> age, and <italic>ε</italic> random error.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt106"><label>a</label><p>
<italic>P</italic>-value with the null hypothesis <italic>β</italic>
<sub>1</sub> = 0 in linear regression model <italic>Y</italic> = <italic>β</italic>
<sub>0</sub>+<italic>β</italic>
<sub>1</sub>
<italic>x</italic>
<sub>1</sub>+<italic>ε</italic>, where <italic>Y</italic> is transcript expression, <italic>x</italic>
<sub>1</sub> age, and <italic>ε</italic> random error.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt107"><label>a</label><p>EASE-score, <italic>P-</italic>value of a modified Fisher's exact test for overrepresentation ##REF##14519205##[35]##. <sup>b</sup>Fold enrichment in each gene group compared to the base set. <sup>c</sup>Gene group with expression that increases with advancing age (same genes used in ##TAB##7##Table 8##, 1430 genes). <sup>d</sup>Hypergeometric test for underrepresentation using annotated genes from a total set of 13,232 genes as a base population ##REF##12620386##[36]##.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt108"><label>a</label><p>EASE-score, <italic>P</italic>-value of a modified Fisher's exact test for overrepresentation ##REF##14519205##[35]##. <sup>b</sup>Fold enrichment in each gene group compared to the base set. <sup>c</sup>Gene group with expression that decreases with advancing age (same genes used in ##TAB##6##Table 7##, 1878 genes). <sup>e</sup>Hypergeometric test for underrepresentation using annotated genes from a total set of 13,232 genes as a base population ##REF##12620386##[36]##.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt109"><label/><p>All data are base pair ±SEM. P-values were calculated from log<sub>10</sub> transformed data using ANOVA. <sup>a</sup>Negatively associated genes were those whose expression decreases with age-at-death (P&lt;0.05). <sup>b</sup>Positively associated genes were those for which expression increased with age-at-death (P&lt;0.05). Unregulated genes were those not significantly correlated with age.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt110"><label/><p>All data are base pair ±SEM. P-values were calculated from log<sub>10</sub> transformed data using ANOVA. <sup>a</sup>Negatively associated genes were those whose expression decreases with age (P&lt;0.05). <sup>b</sup>Positively associated genes were those for which expression increased with age (P&lt;0.05). Unregulated genes were those not significantly correlated with age.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>We are grateful for financial support from The Swedish Medical Research Council, the Swedish Alzheimer Foundation (Alzheimerfonden), and from the Johnnie B Byrd Institute. For the brain sample, data and biomaterials were collected from several National Institute on Aging (NIA) and National Alzheimer's Coordinating Center (NACC, grant #U01 AG016976) funded sites. John A. Hardy, PhD, Marcelle Morrison-Bogorad, PhD, Tony Phelps, PhD and Walter Kukull PhD are thanked for helping to co-ordinate this collection. The directors, pathologist and technicians involved include: National Institute on Aging: Ruth Seemann, John Hopkins Alzheimer's Disease Research Center (NIA grant # AG05146): Juan C. Troncoso, MD, Dr. Olga Pletnikova, University of California, Los Angeles (NIA grant # P50 AG16570):Harry Vinters, MD, The Kathleen Price Bryan Brain Bank, Duke University Medical Center (NIA grant #AG05128, NINDS grant # NS39764, NIMH MH60451 also funded by Glaxo Smith Kline): Christine Hulette, MD, Director, John F. Ervin, Stanford University: Dikran Horoupian, MD, Ahmad Salehi, MD, PhD, New York Brain Bank, Taub Institute, Columbia University (NYBB): Jean Paul Vonsattel, MD, Katerina Mancevska, Massachusetts General Hospital: E. Tessa Hedley-Whyte, MD, Karlotta Fitch, University of Michigan (NIH grant P50-AG08671): Dr. Roger Albin, University of Kentucky (NIH #AG05144): William Markesbery, MD, Mayo Clinic, Jacksonville: Dennis W. Dickson, MD, University Southern California: Caroll A. Miller, MD, Washington University, St Louis Alzheimer's Disease Research Center (NIH #P50AG05681): Dan McKeel, MD, John C. Morris, MD, Eugene Johnson, Jr., PhD, Virginia Buckles, PhD, University of Washington, Seattle (NIH #P50 AG05136):Thomas Montine, MD, PhD, University of Pennsylvania School of Medicine, Alzheimer's Disease Research Center: John Q Trojanowski, MD, Virginia M Lee, MD, Vivianna Van Deerlin, MD, Boston University Alzheimer's Disease Research Center (NIH grant P30-AG13846): Ann C. McKee, Sun Health Research Institute, Arizona (NIA #P30 AG19610): Joseph Rogers, PhD, Thomas G. Beach, MD, PhD, Emory University: Bruce H. Wainer, MD, PhD, Marla Gearing, PhD, University of Texas, Southwestern Medical School: Charles L. White, III, M.D., Roger Rosenberg, University of California, Davis: William Ellis, MD, Rush University Medical Center, Rush Alzheimer's Disease Center (NIH #AG10161): David A. Bennett, M.D. Julie A. Schneider, MD, MS, University of Miami/NPF Brain Endowment Bank: Deborah C. Mash, MD, Margaret J Basile. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group>
<graphic xlink:href="pone.0003024.g001"/> <graphic id="pone-0003024-t001-1" xlink:href="pone.0003024.t001"/> <graphic id="pone-0003024-t002-2" xlink:href="pone.0003024.t002"/> <graphic id="pone-0003024-t003-3" xlink:href="pone.0003024.t003"/> <graphic id="pone-0003024-t004-4" xlink:href="pone.0003024.t004"/> <graphic xlink:href="pone.0003024.g002"/> <graphic id="pone-0003024-t005-5" xlink:href="pone.0003024.t005"/> <graphic id="pone-0003024-t006-6" xlink:href="pone.0003024.t006"/> <graphic id="pone-0003024-t007-7" xlink:href="pone.0003024.t007"/> <graphic id="pone-0003024-t008-8" xlink:href="pone.0003024.t008"/> <graphic id="pone-0003024-t009-9" xlink:href="pone.0003024.t009"/> <graphic id="pone-0003024-t010-10" xlink:href="pone.0003024.t010"/> <graphic xlink:href="pone.0003024.e001.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0003024.e002.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0003024.e003.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0003024.e004.jpg" mimetype="image" position="float"/>
<media xlink:href="pone.0003024.s001.pdf"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0003024.s002.pdf"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0003024.s003.pdf"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0003024.s004.pdf"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0003024.s005.pdf"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["9"], "element-citation": ["\n"], "surname": ["Kluger", "Gianutsos", "Golomb", "Ferris", "George"], "given-names": ["A", "JG", "J", "SH", "AE"], "year": ["1997"], "article-title": ["Patterns of motor impairement in normal aging, mild cognitive decline, and early Alzheimer's disease."], "source": ["J Gerontol B Psychol Sci Soc Sci"], "volume": ["52"], "fpage": ["P28"], "lpage": ["39"]}]
PLoS One. 2008 Aug 20; 3(8):e3024