The invention relates generally to genetic algorithms, and more particularly to the identification of gene expression profile biomarkers and therapeutic targets for brain aging.
Brain aging processes are enormously complex phenomena that affect multiple systems, cell types and pathways, and result in cognitive decline and increased risk of Alzheimer's disease (AD). Landfield P W et al., J Neurobiol 23: 1247-1260 (1992). Although several biological mechanisms have been putatively linked to brain aging or Alzheimer's disease, including inflammation, oxidative stress, Ca2+ dyshomeostasis (Landfield, P W & Pitler T A, Science 226: 1089-1092 (1984); Landfield P W et al., J Neurobiol 23: 1247-1260 (1992)), mitochondrial dysfunction and chronic exposure to adrenal stress hormones (Landfield P W et al., Science 214: 581-584 (1981); Porter N M & Landfield P W, Nature Neurosci 1: 3-4 (1998)), the specific mechanisms and pathways, if any, through which they are linked to impaired brain function are not understood.
It is widely thought that gene expression changes contribute to many aspects of declining function with aging. Finch C E, Longevity, Senescence and the Genome, 37-42 (Univ. Chicago Press, Chicago, 1990). It is also thought that gene expression changes are important for processing and storage of memory. However, not all genes that change expression in the brain with aging are thought to be important for cognition.
Gene-expression changes that specifically contribute to age-related memory decline should selectively change with brain aging and should be correlated specifically with measures of age-associated cognitive decline; that is, a subset of the full set of aging-dependent genes should also correlate with age-related cognitive decline. See, Lockhart D J & Barlow C, Nat Rev Neurosci 2: 63-68 (2001) and Mimics K, Nat Rev Neurosci 2: 444-447 (2001).
If a subset of age-dependent genes also shows expression patterns directly correlated with age-related memory decline, then such a subset of “aging and cognition-related genes” (ACGs) would be extremely helpful as biological indexes (“biomarkers”) for assessing or diagnosing the degree of age-related cognitive impairment in individual subjects. In turn, the ability to measure aging-related cognitive impairment quantitatively is essential for discovering new therapeutic targets, and developing new strategies and pharmaceutical compounds for counteracting normal age-related cognitive decline and/or age-related neurodegenerative diseases, including Alzheimer's disease (AD) or Parkinson's disease (PD).
Identifying ACGs in any mammalian species therefore, might have great therapeutic usefulness. Moreover, because of the well-established homologies of most genes across mammalian species and because of the clear similarities in patterns of brain aging and cognitive decline across species, identification in any mammal would have human health implications. Furthermore, because the primary risk factor for Alzheimer's disease and Parkinson's disease is aging itself, therapeutic approaches developed for aging-related cognitive impairment should also help ameliorate cognitive decline from age-related neurodegenerative disease. Thus, there is a clear need for identifying ACGs but, to date, such genes have not been discovered for any mammal.
Gene microarray technology provides a powerful approach for unraveling the complex processes of aging. To date, however, its impact has been limited by statistical problems, small sample sizes, and difficulty in assessing functional relevance. Moreover, studies that have examined gene expression during brain aging using microarrays have not used sample sizes large enough to provide adequate statistical power for formal statistical testing. Lee C K et al., Nature Genetics 25: 294-297 (2000); Jiang C H et al., Proc Natl Acad Sci USA 98: 1930-1934 (2001) Therefore, even the genes they have reported to change with aging have not been validated by accepted statistical criteria.
The extremely large data sets generated by microarrays pose formidable bioinformatics and resource problems that have to date limited the impact of this powerful technology. Because of these difficulties, most microarray studies have relied on simple fold change comparisons in small samples. However, neither fold change analyses nor the small sample protocols widely used allow the direct estimates of variance necessary for defining type I error (false positives). In addition, fold change criteria, by definition, select for large changes. Therefore, they exhibit low detection sensitivity (high false negatives, or type II error), and are unable to identify the modest changes that often characterize functionally important (and, therefore, tightly regulated) genes. The inability to assign type I error is a particularly critical problem for microarray studies because the thousands of comparisons of gene expression in such analyses greatly increase the expected false positives. For example, even if group sizes were sufficient for formal statistical analyses, and 5000 gene transcripts were each tested by t-test for differences between two conditions at p≦0.05, the false positive rate is equal to the p-value and, consequently, 5% of the 5000 tested transcripts (250) would be expected to be found significant by chance alone.
Although microarray studies have some important offsetting advantages that improve statistical confidence (e.g., co-regulation of genes within a functional group), there is increasing recognition that microarray experiments should generally meet the same statistical standards as other biological experiments or, at least, should systematically estimate the degree of statistical uncertainty. Several strategies to improve statistical confidence have been developed for small-sample microarray studies, but these generally rely on indirect estimates of variance and/or greatly sacrifice sensitivity (i.e., stringent p-values).
Another highly important problem of microarray studies is that of determining which of the hundreds of expression changes that may be observed are likely to be functionally relevant. Correlation analysis is one quantitative approach to linking gene expression with function, although it also requires relatively large sets of independent samples. Expression-function correlations fulfill a key prediction of a causal relationship (i.e., that causally related variables should co-vary) and therefore, can serve as a valuable tool for the identification of candidate functionally relevant genes. Nonetheless, there have been few correlation studies attempting to link cognitive dysfunction with univariate gene expression patterns across individual subjects, much less using the massive amounts of data generated in microarray analyses.
The invention provides a statistical and functional correlation strategy to identify changes in cellular pathways specifically linked to impaired cognitive function with aging. The bioinformatics and functional correlation strategy improves the power of microarray analyses and provides the ability to test whether alterations in specific hippocampal pathways are correlated with aging-related cognitive impairment. The invention is useful for application in large, well-powered groups and for controlling type I error (false positives), enhancing detection sensitivity (reducing type II false negatives) and determining which aging changes in expression are most closely correlated with declining brain function.
Accordingly, the invention provides a method for identifying a biomarker for brain aging, where the biomarker is a polynucleotide or a polypeptide encoded by said polynucleotide. The method involves first obtaining a set of polynucleotides obtained from a set of brain samples (such as hippocampal samples), where the members of the set of brain samples were obtained from members of a set of mammals, wherein the set of mammals contains more than two members, with at least young, mid-aged and aged members, and then identifying the identity and amount of the members of the set of polynucleotides present in the brain samples. The method then involves the steps of deleting certain non-biomarker polynucleotides from the set of polynucleotides, testing by a conventional statistical method (such as) for a significant effect of aging across the young, mid-aged and aged members; and correlating the identity and amount of the members of the set of polynucleotides present in the brain samples with cognitive performance in behavioral tests.
By use of the methods of the invention, one skilled in the genomics art can identify multiple groups of related genes, many representing processes with previously unrecognized relationships to aging and/or cognitive dysfunction. Thus, the invention also provides compositions of matter comprising sets of genes, expressed sequence tags (ESTs), polynucleotides and polypeptides encoded by said polynucleotides identified as being involved in the aging processes. These sets usefully result in a statistically validated, comprehensive overview of mammalian, including human, functional brain aging. In particular, the set of genes can be used for the diagnosis of human age-related disease, such as an age-related neurodegenerative condition, including Alzheimer's disease or Parkinson disease.
The invention provides a set of biomarkers for brain aging, where (a) the set of biomarkers comprises at least two members; (b) the brain expression patterns of the members of the set are significantly altered with aging as determined by a conventional statistical method (such as ANOVA or student's t test), with p<0.05; (c) the brain expression patterns of the members of the set are correlated (using a conventional statistical correlation test, e.g., tested by Pearson's or Spearman's correlation test) across age groups with cognitive performance in behavioral tests, with a correlation of p<0.05 (or with a more stringent correlation of p<0.01 or p<0.001) between brain expression and cognitive performance; and (d) the cognitive performance in behavioral tests significantly altered with aging as determined by a conventional statistical method. The biomarkers may also correlate with a behavioral measure of functional impairment, such as an age-related neurodegenerative condition, including Alzheimer's disease or Parkinson's disease.
The invention also provides a set of at least two biomarkers for brain aging, where where the brain expression patterns of the members of the set are significantly altered with aging as measured by a conventional statistical correlation test at a significance level of p<0.01.
The invention further provides a set of at least two biomarkers for brain aging, where the brain expression patterns of the members of the set are significantly altered with aging as determined by a conventional statistical method, with p<0.05 (or a more stringent correlation, such as p<0.025, p<0.01 or p<0.001).
In one example of the invention, rats in three age groups (Young, Mid-Aged, Aged) were characterized on two memory tasks and each mammal's hippocampal CA1 region was analyzed by a microarray analysis for gene expression. These analyses identified multiple groups of genes, many representing pathways with previously unrecognized relationships to aging and/or cognitive decline. The analysis showed that for all groups, the aging changes in expression began by mid-life.
In one aspect of the invention, the known interactions of the identified processes suggest an integrative model of specific cellular cascades that begin in mid-life and eventually impair cognitive function and increase neuronal vulnerability. Initially decreased neuronal activity and/or oxidative metabolism trigger separate but parallel genomic cascades in neurons and glia. In neurons, the cascade results in reductions of immediate early gene signaling, biosynthesis, synaptogenesis and neurite remodeling. In contrast, glia undergo increased lipid metabolism and mediate a cycle of demyelination and remyelination that induces antigen presentation, inflammation, oxidative stress and extracellular restructuring. Intervention studies based on these findings can identify the cause and effect interactions among the complex processes of brain aging.
The concept of “biomarker” is well-known and useful concept for those of skill in the genomic art. In general, a biomarker is a measurable biological manifestation that is a quantitative indication of the presence and degree of an underlying biological process of interest.
We have devised a multi-stage method for the identification of biomarkers for brain aging, using gene expression microarrays and behavioral testing. The method of the invention allows one skilled in the genomics art to identify both “aging and cognition-related genes” (ACGs) and unique genes that change with brain aging alone, based on formal statistical testing.
As used in this specification, the word “cognitive” is defined as comprising the higher order intellectual/brain processes involved in learning, including attention, acquisition, short-term memory, long-term memory and memory retrieval, among others.
As used in this specification, across different mammalian species, age definitions are as follows: “Young” mammals are those at or beyond reproductive maturity for the species. “Mid-aged” is defined in two ways: at or around half the average lifespan for the species and at or around the midpoint between reproductive maturity and average lifespan. “Aged” mammals are those at or around average lifespan. Animals intermediate between two ages could be considered as part of the group to which they are most closely chronologically related (with the exception of young animals, for whom it would be inappropriate to include prepubescent individuals).
We used the bioinformatics and functional correlation strategy of the invention for microarray analyses. As a result, we were able to detect multiple groups of related genes that were altered by brain aging and also correlated with cognitive function across individual subjects. Most of the shifts in genomic regulation began by mid-life, well before the onset of measurable cognitive impairment, implying that cognitive function is not altered substantially without further progression and/or the cumulative effects of the initial changes in gene regulation.
This analysis depended on a novel combination of three approaches for microarray research: (a) the quantitative measurement of the dependent function of interest (cognitive performance), which provided a basis for large-scale expression-function correlation analyses; (b) the application of formal statistical analyses (ANOVA, Pearson's) to large groups of independent microarray samples, which conferred substantial statistical power and high detection sensitivity for even modest changes (low false negative type II error); and c) systematic estimates of the maximum probabilities of false positives in our data. Our results using the method of the invention provide a generally comprehensive overview of hippocampal genes/processes that are altered with brain aging and closely linked to brain functional decline.
To verify the method of the invention, we first tested young (3-4 months old), mid-aged (12-13 months old) and aged (24-25 months old) rats (n=9-10 per group) for performance on the Morris spatial water maze (SWM) and object memory task (OMT). Both behavioral tests clearly and reliably (statistically) revealed aging-related cognitive impairment (
We then anesthetized (for euthanasia) all animals and dissected out a region of the brain (CA1 region of the hippocampus) known to be important for memory. These brain tissues were then prepared for analyses of gene expression profiles (mRNA content) on Affymetrix GeneChip microarrays specific for the rat genome (RG-U34A arrays) (one array for each individual rat sample). The microarrays were then read and analyzed for expression profile data on an Affymetrix GeneChip System according to the manufacturer's instructions.
The behavioral and microarray methods that were used can reasonably be expected to apply as well to mice as to rats. Similar behavioral and microarray methods known to those of skill in the art can be used for testing of other mammals, including humans. The utility of the method of the invention for human testing is discussed below.
We then transferred the data into standard computer spreadsheets (e.g., Excel) for performing statistical analyses of the effects of aging. Using Analysis of Variance (ANOVA) we defined the set of genes whose degree of expression changed significantly with brain aging. We then used that set of “Aging Genes” and tested each gene's expression profile (across only the aged animals) for significant correlation with memory performance on the Object Memory Task (OMT) as well as the Morris spatial water maze (SWM). The “Aging Genes” whose expression patterns correlated significantly with cognitive performance were defined as the primary subset of “Aging and Cognition-Related Genes” (ACGs), and subcategorized as OMT-associated, SWM-associated, or both OMT and SWM-associated. We further included genes with no behavioral association that had an ANOVA p value≦0.001 since genes identified at this more stringent level are less subject to the error of multiple testing (
Based on those large-scale studies, we have developed a list of ACGs that appear to have considerable potential importance for assessing and generating new treatments for age-dependent functional decline (TABLES 1A and 1B).
These lists contain some genes that were identified previously as being linked to brain aging or neurodegeneration (e.g., inflammation or mitochondrial genes, Lee C K et al., Nature Genetics 25: 294-297 (2000)) but none has been previously shown to be specifically associated with both brain aging and aging-dependent cognitive impairment. Further, many genes on our list have not even been shown previously to be linked to brain aging alone or to cognition alone. Thus, our lists of ACGs are unique and useful biomarkers and therapeutic targets specifically for aging-dependent cognitive impairment. In addition, our list of all genes that change with brain aging contains many genes never before reported to change with brain aging, and therefore provides a useful and unique panel of gene biomarkers and therapeutic targets for study and treatment of brain aging.
In addition to these lists for identified genes, we have also performed the same analyses and compiled the same lists for unidentified expressed sequence tags (ESTs) that are on the same Affymetrix Chips (TABLE 2). These are valuable data, because once the ESTs are identified, they can provide therapeutic targets.
Using the method of the invention, we were able to identify a number of processes and pathways that previously have not been clearly associated with normal brain aging. The most unexpected findings included altered expression profiles suggestive of increased myelin and lipid turnover, as well as widespread changes indicating coordinated downregulation of oxidative metabolism, decreased neurite outgrowth and synaptogenesis. Other novel genes we identified appear to suggest alterations in general metabolic and biosynthetic chaperone functions. In addition, many of the identified groups confirmed previously described changes in expression for genes regulating several major processes (e.g., inflammation, glial reactivity, oxidative stress). However, our results also extend the earlier findings considerably by revealing the extent of the changes and the concurrent upregulation of potentially orchestrating transcription factors and cytokines that may provide important clues to pathogenic mechanisms.
In order to begin to develop an integrative overview of potential interactions among the multiple altered expression patterns observed here, we considered functional implications at the pathway level. Our interpretations rely on the functions that have been previously associated with many of the genes identified by those of skill in the genomics art. These are identified through PubMed literature searches, annotations provided by Affymetrix, entries in the SwissProtein database <http://www.expasy.ch/sprot/sprot-top.html> and associations reported in the Genome Ontology (GO; <http://www.geneontology.org>). We also rely on the general assumption held by those of skill in the genomics art that similar changes in the expression of multiple genes of a particular pathway imply like changes in the functions mediated by the encoded proteins of that pathway. Gene expression changes also can reflect compensatory negative feedback regulation (or other dissociations of gene expression and protein function), but the potential confound of dissociation is presumably less of a problem in microarray analyses in which multiple genes in a pathway are observed to change in the same direction. Some of the primary metabolic pathways and processes considered in the interpretations are depicted in TABLE 1.
Functional Groups. We found age-dependent upregulation of many ACGs involved in inflammatory/immune/stress responses and downregulation of many involved in energy metabolism. In addition, we found alterations of gene expression reflecting multiple categories/pathways not previously recognized to change with normal aging. These included upregulation of genes for myelin proteins, cholesterol biosynthesis and transport, amino acid metabolism, intracellular Ca2+ signaling, and protein processing, as strongly suggesting an ongoing cycle of remyelination and demyelination. We also found widespread downregulation of genes for biosynthesis, immediate early responses, and synaptic structural plasticity, suggestive of neuronal involution. Multiple transcriptional regulators and cytokines were also identified that may play orchestrating roles. Nearly all expression changes began by mid-life but cognition was not impaired until late life. Upregulated genes for inflammation and intracellular Ca2+ release were among those most closely correlated with impairment.
For TABLE 1A, “GenBank” is the gene accession number established at the web accessible GenBank database <http://www.ncbi.nlm.nih.gov/>, The “Description” includes a ‘common name’ (if applicable) as well as a brief description of the gene product. Values for Young, Mid-Aged, and Aged categories are the mean±SEM of expression values. Genes are put into functional categories (see, above) and grouped by their level of association with behavior (expression correlated significantly (Pearson's; p≦0.025) with both tasks, with the OMT, with the SWM, or with none of the tasks but highly significant across age (p≦0.001 on ANOVA across age, p>0.025 for correlation on both SWM and OMT). Within each level of association, genes are ranked by the significance of the age-dependent change in their expression level (ANOVA; p≦0.025). Asterisked (*) genes are those that also showed a significant behavioral correlation (Pearson; p≦0.025).
ACGs that were downregulated with aging (TABLE 1A) appeared primarily to represent metabolic and neuronal functions (
Metabolism. Multiple genes related to functions of the mitochondrial electron transport chain (e.g., glycerol 3-phosphate dehdrogenase, NADH dehydrogenase, Rieske's iron-sulpher protein) were downregulated with aging (TABLE 1A). Moreover, we found aging-dependent downregulation of several genes related to pathways important for glucogenic amino acid catabolism, including glutaminase and arginase (TABLE 1A).
Synaptic Structural Plasticity. One of the most prominent categories of identified genes showing decreased expression and behavioral correlation was that comprising genes involved in synaptic structural plasticity, including neurite outgrowth and synaptogenesis (e.g., decreased expression of genes encoding agrin, GAP-43, Homer 1a, Narp, Arc, etc.) (TABLE 1A). Many of these genes are activity-dependent in neurons and have been linked previously to synaptic plasticity, neurite remodeling or learning in univariate studies (e.g., Biewenga J E et al., Acta Biochim Pol 43: 327-38 (1996); Steward O et al., Neuron 21: 741-51 (1998), Mantych K B & Ferreira A, J Neurosci 21: 6802-9 (2001), Guzowski J F et al., J Neurosci 20: 3993-4001 (2000), Bezakova G, et al. Proc Natl Acad Sci USA 98 9924-9 (2001)), although Gap-43 is one of the few reported so far to change with aging. Similarly, many other neural activity-dependent genes, including IEGs in the Transcription Regulators and Signaling categories (e.g., Egr1, Egr 2, MAPKK, etc.), showed decreased expression with aging and were correlated with impaired cognition (TABLE 1A).
In addition, multiple genes important for general growth and biosynthetic mechanisms, chaperone functions and protein processing were also downregulated with aging (e.g., hsp60, histone H2AZ, proteasome subunit, DNA J-like homolog, etc.) as were specific neuronal signaling genes (e.g., GluR 5-2, the kainate receptor; and neuropeptide Y) (TABLE 1A). These widely downregulated biosynthetic and signaling genes appear to reflect a general involution of metabolic and neurite structural remodeling processes in neurons (e.g.,
The analyses for TABLE 1B are as described for TABLE 1A.
Upregulated Genes. Genes that were upregulated with aging and negatively correlated with behavior fit primarily into categories that appeared to reflect activated glial functions (
Lipid Metabolism. Multiple genes important for mitochondrial and cytosolic lipid β-oxidation (e.g., carnitine palmitoyltransferase, lecithin-cholesterol acyltransferase, etc.; TABLE 1B), the primary pathway for free fatty acid (FFA) catabolism, were upregulated.
Increased Myelin Synthesis, Cholesterol Biogenesis and Vesicle Transport. Importantly for identifying the trigger mechanism for elevated lipid catabolism, the expression of many genes encoding myelin-related proteins or myelin-related transcription factors on the microarray was increased with aging (and several also were correlated with cognitive impairment) (TABLE 1B). These observations strongly suggest that a major increase in myelin synthesis programs developed with aging. This interpretation is also supported by the upregulation of multiple genes important in lipogenesis for cholesterol biosynthesis (Add 1/SREBP1), and the packaging/transport of cholesterol esters and other complex lipids (ApoD, LCAT, etc.) (TABLE 1B). Recent studies have shown that stimulation of myelin synthesis programs in oligodendrocytes is associated with induction of genes for both myelin proteins and lipogenic pathways (Nagarajan R et al., Neuron 30 355-68 (2001)).
Cyloskeleton/Vesicles. Moreover, expression of genes related to actin assembly, transport or fusion of packaged vesicles (actin related complex, rsec8, tubulin, and syntaxin 7) was increased (TABLE 1B). These molecules are associated with vesicle transport and fusion in neurons. In addition, however actin assembly proteins are also known to play a major role in myelin vesicle transport in oligodendrocytes (Madison D L et al., J Neurochem 72: 988-98 (1999)). Given the upregulation of myelin programs and the downregulation of synaptic plasticity genes, therefore, the age-dependent upregulation of genes linked to vesicle transport capacity seems more likely to be associated with enhanced myelin transport in oligodendrocytes. Further support for the view that extensive oligodendrocyte activation and/or synthesis occurs in hippocampal aging is provided by the observation that many genes that were upregulated with aging are preferentially expressed in oligodendrocytes (e.g., myelin proteins, FAH, PGD-S, etc.) (e.g., Labelle Y et al., Biochim Biophys Acta 1180: 250-6 (1993)).
Myelin also is normally degraded to free fatty acids through the endosomal-lysosomal pathway. Consistent with elevation of myelin degradation, we also found increased expression of Cathepsin S and other genes encoding lysosomal enzymes (TABLE 1B). Cathepsin S is particularly important in the processing of antigenic myelin fragments.
Amino Acids. In contrast to enzymes for glucogenic amino acids (TABLE 1A), expression was upregulated for multiple genes encoding enzymes related to the metabolism of the ketogenic/glucogenic amino acids, tyrosine, phenylalanine and tryptophan (e.g., DHPR, KAT, FAH, see, TABLE 1B). Catabolism of ketogenic amino acids yields either acetoacetate or one of its precursors (e.g., acetyl CoA), which can be used either for energy metabolism or lipogenesis. Upregulation of DHPR, which catalyzes the formation of a critical cofactor (tetrahydrobiopterin) for tyrosine and monoamine synthesis, and concomitant upregulation of MAO-B (TABLE 3), together suggest elevated metabolism of tyrosine and tryptophan via greater monoamine turnover.
Inflammation/Defense/Immunity. There was massive upregulation of expression of genes encoding MHC class I antigen presenting molecules, and numerous other inflammatory/immune proteins (TABLE 1B). Genes in the inflammation category exhibited some of the most robust monotonic changes with aging seen in our results using the method of the invention (e.g., most were significant at the p<0.001 criterion with 0.025 FDR) (TABLE 1B). Moreover, most were inversely correlated with cognitive function (
Consistent with evidence of a role for oxidative stress in brain aging (Carney J M et al., Proc Natl Acad Sci USA 88: 3633-6 (1991), Hensley K et al., Ann NY Acad Sci 786: 120-34 (1996), Bickford P C et al., Brain Res 866: 211-7 (2000), Lee C K et al. Nat Genet 25: 294-7 (2000), Jiang C H et al., Proc Natl Acad Sci USA 98: 1930-4 (2001)), we also found increased expression for molecules important in defense against oxidative stress (GST, GSTa1) (TABLE 1B). One potentially key new finding here, as noted above, was that DHPR was upregulated with aging and correlated with cognitive decline (TABLE 1B). Its product, tetrahydrobiopterin, is also an essential cofactor for nitric oxide synthase (Boyhan A, et al. Biochem J 323 (Pt 1) 131-9 (1997)). Because oxyradicals formed from nitric oxide appear to play a major role in inflammatory neuronal damage (Bal-Price A & Brown G C, J Neurosci 21: 6480-91 (2001), Calingasan N Y & Gibson G E, Brain Res 885: 62-9 (2000)), this may be an important pathway through which the deleterious effects of inflammation are mediated in brain aging.
Glial Markers. Astrocyte reactivity and astrocyte markers are also well recognized to increase in the aged rodent and human hippocampus (Landfield P W et al., Science 214: 581-4 (1981), Landfield P W et al., J Neurobiol 23: 1247-60 (1992), Nichols N R et al., Neurobiol Aging 14: 421-9 (1993), Finch C E & Longo V D, Neuroinflammatory Mechanisms in Alzheimer's Disease: Basic and Clinical Research, 237-256 (2001)) and the present data confirm extensive upregulation of genes (Finch C E & Tanzi R E Science 278: 407-11 (1997)) for glial markers (e.g., vimentin, GFAP-cytoskeleton category, TABLE 1B). In addition, we extended those observations to show that genes for proteoglycans (TABLE 1B) and other extracellular proteins (e.g., fibronectin) that are components of astroglial scars also were upregulated. These changes may reflect astroglial-mediated reorganization of the extracellular matrix, a process known to be unfavorable for axonal remodeling.
Signal Transduction. Several genes in calcium regulating and G-protein-coupled signaling pathways were also identified (TABLE 1B). In particular, S100A1, which modulated Ca2+-induced Ca2+ release, and PI 4-kinase, which acts to produce IP3 were upregulated. Several other S100-related genes (e.g., S100A4 and P9K2; TABLE 3) were also upregulated with aging but failed to meet the strict criteria set forth herein (
Biosynthesis. Concomitantly, many ribosomal (growth) and protein processing genes were upregulated (TABLE 1B). The upregulated changes reflect increased protein synthesis, turnover and phagocytosis associated with strongly elevated biosynthetic processes in glial compartments (e.g., elevated myelin, MHC, proteoglycan synthesis).
Orchestrating Factors. Our data show that a number of transcriptional regulators and cytokines, including KZF-1, Roaz and members of the NFI family (TABLE 1B) were upregulated and therefore, may be strong candidates for coordinating factors. Under some conditions, several of these factors function as negative transcriptional regulators.
Relationship to Fold Change. The large majority of microarray analyses to date have used fold-change criteria to detect changes in expression. In addition to providing little basis for statistical assessment (e.g., Miller R A, et al. J Gerontol A Biol Sci Med Sci 56 B52-7 (2001)), however, fold-change criteria are relative insensitive. Among the 139 ACGs, most exhibited group mean fold changes between the Young and Aged groups of less than 1.5 (92), a few showed fold changes between 1.5 and 2.0 (26), and only a handful of genes exceeded 2-fold-change (20) (TABLES 1A and B). Thus, few of our results using the method of the invention would have been detected in the great majority of prior microarray studies, in which 1.7 to 2-fold change cutoffs are commonly used as minimum criteria for identifying differences, and many changes are reported in the 3-4 fold range. Further, the rank order correlation between group mean fold-change and p values on the ANOVA for all aging-significant genes, although significant, was modest according to Spearman's correlation test (Spearman's r=0.45, p<0.001). Armitage P & Berry G, Statistical Methods in Medical Research, 2nd Edn., 200-205 (1987). This indicates that fold-change accounted for only 20% of the variance (r2) in the degree of statistical significance on the ANOVA. Some of our results detected with the enhanced sensitivity of statistical analysis were extremely subtle (e.g., 1.1 fold for the L28 and L29 ribosomal proteins, TABLE 1B). Despite this enhanced sensitivity, however, numerous false negatives were still undoubtedly present in our data set.
Age Course of Gene Expression Changes. Using a design with three age groups enabled us to classify genes and categories according to their general patterns of age dependence of change (
Almost all categories comprising downregulated and cognitively correlated genes (TABLE 1A), exhibited their greatest change between the Young and Mid-Aged points, and many did not show much additional downregulation between the Mid-Aged and Aged groups (
Strongest Correlations of Pathways with Memory Performance. To determine which pathways were most closely correlated with memory performance, we calculated the percentage of genes in each of our categories that were correlated significantly (at p<0.025) with both memory tests. We reasoned that each test measures aspects of memory but each test also has its own error sources and confounding contributions from non-cognitive performance factors. Therefore, genes that correlated with both tasks seem more likely to be associated with cognitive processes.
Because memory performance changed most between the Mid-Aged and Aged groups (
Another way to examine closeness of correlation specifically with memory impairment is to correlate gene expression with performance only in the aged group. This correlation focuses on variation in the performance of aged animals and removes the overall age course pattern from contributing to the correlation with impairment. This correlation is independent of the ANOVA for aging effects and an FDR also can be calculated. Consequently, we tested each of the 139 primary aging- and behaviorally-related genes for correlation with 24 hr memory performance on the OMT in the aged group. The OMT was selected over the SWM for this test as it had the greater dispersion of performance needed for correlation analysis. The correlation tests in the aged group (n=10) of course had considerably less power than across all three groups (n=29) and the criterion for significance was set at p<0.025.
Only 3 (4.9%) of the downregulated ACGs, but 10 (12.2%) of the upregulated ACGs were correlated with Aged group performance on the OMT. The FDR for these genes was 0.28. Two of the 3 downregulated ACGs were accounted for by the Synaptic Structural Plasticity category (Fez-1, agrin). For upregulated genes, two of the 10 ACGs were from Inflammation (MHC and CD59 antigen), three from Cytoskeleton/Vesicle category (Vcp, rsec8 and p41-Arc), and three from Growth/Biosynthesis (2 ribosomal proteins and CD24 antigen). No other category had more than one, including Transcription (NF1-A) and Protein processing and trafficking (Skd3).
Thus, by the criterion of correlation on both tasks, the upregulated categories of Inflammation/Immune, signaling (particularly Ca2+ signaling), Cytoskeleton/Vesicle and Amino Acid Metabolism were ranked most highly. By the criterion of correlation in the aged group only, the upregulated categories of Cytoskeleton/Vesicle ({fraction (3/6)}), Biosynthesis ({fraction (3/12)}) and Inflammation ({fraction (2/20)}), and the downregulated category of Synaptic Plasticity ({fraction (2/7)}) were ranked most highly.
Benefits of the Invention. One of the major problems associated with developing treatments for aging-dependent functional decline is the lack of good genomic biomarkers or targets of brain aging needed for evaluating the efficacy of different treatments. Our ACGs, therefore, could serve as excellent biomarkers of cognitive aging. Using microarrays constructed to contain oligonucleotide sequences specific for hybridization with and measurement of mRNAs of the identified ACGs, laboratory animals could be assessed for degree of cognitive aging before, during and after treatment with a compound. Treatments that slowed or reversed the ACG profile during aging might be highly promising for development as new therapeutic approaches. Further, treatments that slowed or reversed expression profiles of particular genes in our panel of biomarkers might reveal which specific genes among the subset of ACGs are most critical for the age-dependent functional decline and, therefore, would suggest genes and gene products that should be targeted with high priority for development of therapeutic interventions. The same approach could be applied using our panel of unique brain aging genes that are not specifically clustered with cognition related genes, to evaluate and develop new therapies and compounds for treatment of brain aging in general.
The panel of ACGs identified here can be used on a microarray to perform diagnostic tests. Subjects suspected of having accelerated brain aging or early age-related neurodegenerative disease could provide a small brain biopsy sample for testing by microarray. This could then determine the subject's suitability for pharmacologic intervention.
Based on the gene lists described above, investigators can develop new drugs or treatments aimed at altering the activity of one or more genes in the lists, or products encoded by those genes, or targets of the products, with the goal of counteracting age-related cognitive impairment or brain aging in general.
A smaller subset of ACGs, specifically linked to some process or system (e.g., to inflammation, mitochondrial function, or lipid metabolism, etc.), could be used in a microarray to test efficacy of a new compound targeted to slowing or reversing aging and cognitive changes dependent on that set of genes or gene-impacted systems, either in experimental tests to develop new compounds, or as diagnostic or therapeutic guides.
Relevance to Human Brain Aging and Alzheimer's Disease. Normal human brain aging is associated with memory dysfunction and appears to set the stage for Alzheimer's disease and other age-related neurodegenerative conditions. It also shares many features with animal models of aging. Landfield P W et al., J Neurobiol 23: 1247-1260 (1992). Thus, many of the memory-correlated gene expression profiles seen here in rats may have implications for genomic mechanisms of human brain aging and/or Alzheimer's disease. This view is supported by several parallels between processes identified here and those seen in human aging or Alzheimer's disease. For example, myelin abnormalities are also found extensively in normal brain aging in humans (leukoaraiosis). These white matter changes in humans are also correlated with cognitive dysfunction and become more severe in disease states. Further, cerebral metabolism begins to decline by mid-life in humans, much as it apparently does in rats (
It is interesting, in view of the apparently altered lipid metabolism seen here, that activity of the cholesterol ester synthesizing enzyme acyl CoA:cholesterol acyltransferase (ACAT) is elevated in Alzheimer's disease and appears directly coupled to amyloid production. ACAT has lipogenic functions somewhat similar to those of LCAT, which was also upregulated here (TABLE 1A). Moreover, activity of glycerol-3-phosphate dehydrogenase (GPDH) is elevated in association with abnormal glucose metabolism in brains of patients with Down's syndrome. The gene encoding this glycolytic enzyme was also upregulated here (TABLE 1A). Other processes found in human aging or Alzheimer's disease brain for which we found corollaries in gene expression include, as noted, inflammation, oxidative stress and elevated KatII (kynurenine aminotransferase 2), among others. Thus, if these parallels depend, at least in part, on similar mechanisms, our our results show that widespread genomic regulatory changes would reasonably be expected to contribute to altered cerebral metabolism, lipid synthesis, neural activity and myelination in human brain aging as well.
Implications for a New Hypothesis of Brain Aging. Based on the functional implications of our results, as discussed above, we provide a new working model of brain aging (
The details of one or more embodiments of the invention are set forth in the accompanying description above. Although any methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present invention, the preferred methods and materials are now described. Other features, objects, and advantages of the invention will be apparent from the description and from the claims. In the specification and the appended claims, the singular forms include plural referents unless the context clearly dictates otherwise. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. All patents and publications cited in this specification are incorporated by reference.
The following EXAMPLES are presented in order to more fully illustrate the preferred embodiments of the invention. These examples should in no way be construed as limiting the scope of the invention, as defined by the appended claims.
Behavioral Results
Thirty animals in three age groups (n=10/group) were trained sequentially on two tasks, first in the Morris spatial water maze (SWM) and then in the object memory task (OMT). Male Fischer 344 rats aged 4 months (Young, n=10), 13 months (Mid-Aged, n=10) and 24 months (Aged, n=10) were used. Overall, the training/testing lasted seven days, and hippocampal tissue was collected 24 hr later. Training or testing occurred on each day except for the 2nd and 3rd days of the seven-day sequence.
Methods used here for cognition assessment in the Morris Spatial Water Maze (SWM), a task sensitive to both hippocampal function and aging, have been described previously Norris C M & Foster T C, Neurobiol Learn Mem 71, 194-206 (1999). Briefly, rats were trained in a black tank, 1.7 M in diameter, filled with water (27±2° C.). Behavioral data were acquired with a Columbus Instruments tracking system. After habituation to the pool, animals were given cue training with a visible platform (five blocks of three trials, maximum of 60 sec/trial, 20 sec intertrial interval and a 15 min interval between blocks). Rats remained in home cages under warm air after each block. Cue training was massed into a single day and the criterion for learning was finding the platform on 4 of the last 6 trials. For all animals that met this criterion, spatial discrimination training was initiated three days later in which the escape platform was hidden beneath the water but remained in the same location relative to the distal cues in the room. Fifteen min following the end of spatial training, a 1-min duration free-swim probe trial with the platform absent was administered, during which crossings over the former platform site (platform crossings) were recorded to test acquisition, followed by a refresher training block. Retention for platform location was again tested 24-hr later using a second 1-min free-swim probe trial.
During cue training in the SWM, all animals were able to locate the visible escape platform according to our criteria and therefore, were trained on the hidden platform spatial task. During acquisition, Aged animals performed more poorly (longer latencies) than Mid-Aged or Young. In addition, an aging-dependent decrease in 24 hr retention, measured by platform crossings (1-way ANOVA, p<0.05), was observed on the retention probe trial (
Methods used here for cognition assessment in the object memory task OMT) have been described previously by Ennaceur A & Delacour J, Behav Brain Res 31: 47-59. (1988). The object memory task (OMT) is also both sensitive to hippocampal function and affected by aging but is less dependent on physical strength and endurance. On the afternoon of the final spatial maze probe trial, animals were administered a habituation session (15-min) in the empty mesh cage to be used for the OMT (63.5 cm×63.5 cm). OMT training began 24 hr after habituation and consisted of a 15-min acquisition session during which two 3-dimensional objects were placed at opposite sides of the cage, followed by two 15-min retention test sessions at 1 and 24 hr posttraining. During the acquisition session, the cage contained two sample objects (A and B) and the time spent actively exploring each object was recorded. After 1 hr, the rat was reintroduced into the cage and the time spent exploring a novel object, C, relative to the familiar object, B, was recorded. On the 24 hr test, familiar object A was reintroduced and object B was replaced by a second novel object, D. Objects were randomized across individuals and timed measures of exploration were used to calculate a memory index (MI) as follows: MI=(N−F)/T, where N is time spent exploring the novel object, F is time spent exploring the familiar object, and Tis total time spent exploring the two objects. More time spent exploring the novel object (higher MI) is considered to reflect greater memory retention for the familiar object.
In the OMT, Aged animals performed as well as Young or Mid-Aged on the 1 hr retention test (not shown), but there was a significant age-related decline in recall (1-way ANOVA, p<0.001, for the main effect of age) on the 24 hr test (
Gene Microarray Chip Results
Microarray analyses were performed on hippocampal CA1 tissues from each of the same behaviorally characterized 30 animals (one chip per animal), but one chip was lost for technical reasons, leaving a data set of 29 microarrays (Young=9, Mid-Aged=10, Aged=10). For tissue preparation, twenty-four hours after completion of the OMT testing, animals were anesthetized with CO2 gas and decapitated. The brains were rapidly removed and immersed in ice-cold, oxygenated artificial cerebrospinal fluid consisting of (in mM): 124 NaCl, 2 KCl, 1.25 KH2PO4, 2 MgSO4, 0.5 CaCl2, 26 NaHCO3, and 10 dextrose. Hippocampi were removed and the CA1 region from one hippocampus per animal were dissected by hand under a stereomicroscope. The CA1 tissue block from each animal was placed in a microcentrifuge tube and flash frozen in dry ice for RNA isolation.
For RNA isolation, total RNA was isolated using the TRIzol reagent and following the manufacturer's RNA isolation protocol (Gibco BRL, #15596). One ml of TRIzol solution was added to each tube containing the frozen tissue block and the tissue was homogenized by 10 passages through an 18½ G syringe needle. After centrifugation, the RNA was precipitated from the aqueous layer, washed and redissolved in RNase-free water. RNA concentration and the integrity of RNA were assessed by spectophotometry and gel electrophoresis. The RNA samples were stored at −80° C.
Gene expression analyses were performed using the Affymetrix GeneChip System. The labeling of RNA samples, rat GeneChip (RG-U34A) hybridization and array scanning were carried out according to the Affymetrix GeneChip Expression Analysis Manual (r.4.0, 2000). Each animal's CA1 RNA was processed and run on a separate rat gene chip. Briefly, an average yield of 40 μg biotin-labeled cRNA target was obtained from 5 μg of total RNA from each CA1 sample, of which 20 ug cRNA was applied to one chip. The hybridization was run overnight in a rotating oven (Affymetrix) at 45° C. The chips were then washed and stained on a fluidics station (Affymetrix) and scanned at a resolution of 3 μm in a confocal scanner (Agilent Affymetrix GeneArray Scanner).
Each U34A rat chip (Affymetrix, Santa Clara, Calif.) contained 8,799 transcript probe sets (gene representations). Although the measured signal intensity for a transcript probe set (Methods) reflects mRNA content, it is referred to here as “gene expression”. However, it is well recognized that mRNA stability and other factors in addition to gene transcription can affect mRNA content.
We used the microarray suite (MAS 4.0) software (Affymetrix) to calculate the overall noise of the image (the amount of variation around the mean intensity, Qraw) for each array. Overall noise was highly similar across arrays in all 3 age groups (Young: 21.81±1.55; Mid Aged: 21.25±2.24; Aged: 20.66±2.06, N.S.). “All probe set scaling” was used to set overall intensities of different arrays to an arbitrary target central intensity of 1500. Thus, the average intensity of each array was adjusted to the 1500 value using a scaling factor (SF). There was no significant difference in SF across ages (Young: 1.58±0.14; Mid Aged: 1.46±0.20; Aged: 1.63±0.16, N.S.).
The algorithm used to determine Presence/Absence is listed in the Microarray Suite 4.0 Manual and is the basis upon which a particular transcript is determined to be reliably detectable by a given probe set. Average difference scores, the average of the difference in expression intensity (ADEI) of each probe pair within a probe set, formed the basis for determining expression (relative abundance) of transcripts, and throughout the text the term “expression level” refers to the ADEI score. When comparing across appropriately normalized arrays, the larger the ADEI score, the greater the relative expression for that particular message. However, ADEI scores are not comparable for relative expression levels among different messages on the same chip, as there are several other factors that can confound such an assessment (e.g., p. 356, Affymetrix Microarray Suite 4.0 User's Guide).
The Presence/Absence calls and Average Difference scores for all probe sets on all 29 arrays were then copied from the MAS pivot table to an Excel 9.0 (Microsoft, SR-1) workbook. From within Excel, the following data manipulations were performed.
Min-Max: For the purposes of filtering (
where x is ADEI score, {overscore (X)}min is the mean for the age group with the lowest ADEI score, and {overscore (X)}max is the mean for the age group with the largest ADEI score. Thus, normalized mean values varied between 0 (lowest) and 1 (highest) for each probe set.
Standardization (Z-score): For the purpose of obtaining the mathematical means within functional categories and graphing, the data was normalized using the Z-score method:
where {overscore (X)} is the mean, and SD(x) is the standard deviation of ADEI across all age groups for an individual probe set.
Statistical Analysis. All statistical tests were performed using a combination Excel (Microsoft, version 9, SR-1) and Sigma Stat (SPSS, version 2).
Multi-Step Gene Identification Algorithm
The analytic algorithm of the invention, which addresses the bioinformatics issues noted above, comprise three main steps aimed first, at reducing the number of comparisons (to manage type I error), second, at reliably detecting modest aging differences with global statistical analyses (by ANOVA), and third, at identifying aging-related expression changes that were quantitatively correlated with cognitive function (by Pearson's test; Armitage P & Berry G Statistical Methods in Medical Research, 2nd Edn., 200-205 (1987)) (see,
Multiple Comparison Reduction Step. The expected false positives in a series of multiple comparisons (false positive rate) are predicted to be a percentage of the total statistical comparisons to be made, as defined by the p-value (i.e., tests at p<0.05 will on average generate 5% false positives). Accordingly, the absolute numbers of expected false positives can be decreased simply by reducing the total transcript sets that are tested in a microarray analysis. This can be done by deleting all transcripts identified a priori as not likely to be relevant to the specific interests of the analysis.
Using this step of the method, we reduced the total transcripts to be tested in three phases. In the first phase, we deleted quality control oligonucleotide sets (“control”, n=60) and all gene transcripts (probe sets) rated “absent” by our criteria. As used in this specification, the term “quality control oligonucleotides” are those oligonucleotides and polypeptides used to test for the appropriate behavior of the technological system, rather than to measure expression levels of biological interest.
Of the original 8,799 sets, 4,118 gene transcript sets were removed at this stage, leaving 4,681 transcript sets that were called “present” for further consideration (
Group Statistical Testing Step (ANOVA). In this second main step of the algorithm, each of the remaining 1,985 transcript sets was tested by 1-way ANOVA for a significant effect of aging (at p≦0.025) across the 3 age groups (n=9-10/group). Of the 1,985 tested sets, 233 were found to change significantly with aging (observed total positives). As noted, at p<0.025, approximately 2.5% (50) of the 1,985 tested should be significant by chance alone (expected false positives). In order to estimate the proportion of false discoveries anticipated among our 233 observed positives (i.e., the fraction of observed positives expected to be false), we used the expected false positive value to calculate the false discovery rate (FDR) (Benjamini et al., Behav Brain Res 125: 279-284 (2001)). For any multiple comparison, the false discovery rate provides an empirical estimate of the anticipated chance error rate among all positives detected. It is partly analogous to the p value of statistical tests, in that the false discovery rate yields the probability that any positive found at the alpha level used (in this case p<0.025) is positive by chance alone.
For the ANOVA-positive results, the FDR was 50/233=0.21, indicating that up to 21% of the observed positives might be positive by chance alone or, that any one positive had a 21% chance of being a false positive.
In addition, we examined the FDR obtained using two other ANOVA p-value levels, p<0.01 and p<0.001. At the p<0.01, ˜20 genes should be found positive by chance alone among the 1,983 transcripts tested. A total of 145 total positives were observed, yielding an FDR of 20/145=0.14. At p<0.001 only 2 false positives are expected in 1,983 tests, and 70 total positives were found. This yields a FDR of 2/70=0.03. The latter, in particular, compares highly favorably with the 0.05 alpha level conventionally accepted for statistical significance in univariate analyses.
However, as noted, additional confidence and validation is gained in microarray analyses when similar patterns of regulation are found among multiple functionally similar genes (Prolla et al., J Gerontol A Biol Sci Med Sci 56: B327-330 (2001)). This is because such genes are not necessarily independent and their co-regulation can provide added cross-validation (e.g., Mimics et al., Neuron 28: 53-67 (2000); Prolla et al., J Gerontol A Biol Sci Med Sci 56: B327-330 (2001)). Consequently, in many cases, confidence advantages can be gained by relaxing p-value criteria in order to expand the numbers of genes included in functional categories. Mimics K, Nat Rev Neurosci 2: 444-447 (2001). Further, relaxing stringency of the p-value reduces the likelihood of type II error (false negatives). Based on these rationales, we used the set of 233 genes obtained at the less stringent p≦0.025 alpha level (rather than the set of 70 at p≦0.001) for the next main step of our algorithm, the behavioral correlation analysis (
Cognitive Performance Correlation Step (Pearson's Test). In this step we identify a specific subset of the 233 aging-significant (by ANOVA) genes that was also correlated with memory performance in both the OMT and SWM. We tested each of the 233 ANOVA-significant genes across animals for statistical correlation between that gene's expression value and behavioral scores (with Pearson's test).
The expression of 161 of the 233 ANOVA-significant genes was correlated significantly with behavioral performance on memory-dependent tasks (p≦0.025; ACGs). Of these, 84 were significantly correlated with both OMT and SWM performance, 40 were significantly correlated with OMT, and 37 were significantly correlated with SMW (
Because of the voluminous literature involved, many relevant citations are not included here. In addition to the “dual function” status of some genes, the functions of many are not completely understood, and therefore, the categorization here, while generally consistent with published reports, is not definitive.
ACGs that were downregulated with aging (TABLE 1A) appeared primarily to represent metabolic and neuronal functions. A substantial number of them fell into the category of oxidative metabolism (TABLE 1A). Many also fell into categories of synaptic/neuritic remodeling or other activity-dependent neuronal processes, e.g., immediate early genes (IEGs) (TABLE 1A). Conversely, ACGs that were upregulated with aging fit primarily into categories that appeared to reflect activated inflammatory response (TABLE 1B).
Additionally, among the main unexpected findings was a widespread upregulation in the expression of multiple genes encoding proteins for myelin synthesis (TABLE 1B) and lipid turnover (TABLE 1B). These various categories, overall, are consistent with a downward shift of oxidative metabolism in parallel with a major upregulation of lipid metabolism.
Genes Identified by the Method of the Invention
The following tables provide additional results from the tests performed above, and supplement the results presented in TABLES 1A and B.
Using the method of the invention, we have identified a set of genes and ESTs that changed with age by ANOVA (p≦0.05), but which are not ACGs. These include AA685974 (EST108806 cDNA); AA799396 (EST188893 cDNA); AA799479 (ESTs, Highly similar to NADH-ubiquinone oxidoreduct.); AA799481 (EST188978 cDNA); AA799529 (EST189026 cDNA); AA799599 (EST189096 cDNA); AA799636 (EST189133 cDNA); AA799680 (EST189177 cDNA); AA799724 (ESTs, Highly similar to DNA-directed RNA polymerase I); AA799773 (EST 189270 cDNA); AA799779 (acyl-CoA:dihydroxyacetonephosphate acyltransferase); AA799858 (EST189355 cDNA); AA800026 (EST189523 cDNA); AA800318 (EST189815 cDNA); AA800622 (EST190119 cDNA); AA800693 (EST190190 cDNA); AA800948 (Tuba4); AA801286 (Inositol (myo)-1(or 4)-monophosphatase 1); AA817887 (profilin); AA818025 (CD59 antigen); AA818240 (Nuclear pore complex protein); AA818487 (cyclophilin B); AA819500 (ESTs, Highly similar to AC12_HUMAN 37 kD subunit); AA819708 (Cox7a3); AA848831 (lysophosphatidic acid G-protein-coupled receptor, 2); AA852046 (EST194815 cDNA); AA859545 (ESTs, Weakly similar to hypothetical protein C09H6.3); AA859562 (UI-R-E0-bv-b-03-0-UI.s1 cDNA); AA859643 (UI-R-E0-bs-a-08-0-UI.s1 cDNA); AA859645 (attractin); AA859690 (UI-R-E0-bx-e-11-0-UI.s1 cDNA); AA859848 (UI-R-E0-cc-h-10-0-UI.s1 cDNA); AA859954 (Vacuole Membrane Protein 1); AA859980 (T-complex 1); AA860030 (UI-R-E0-bz-e-07-0-UI.s2 cDNA); AA866257 (ESTs); AA866299 (UI-R-A0-ac-f-12-0-UI.s3 cDNA); AA866299 (UI-R-A0-ac-f-12-0-UI.s3 cDNA); AA866432 (UI-R-E0-ch-e-06-0-UI.s1 cDNA); AA866477 (UI-R-E0-br-h-03-0-UI.s1 cDNA); AA874830 (UI-R-E0-cg-f-04-0-UI.s1 cDNA); AA874874 (ESTs, Highly similar to alcohol dehydrogenase class III); AA874969 (ESTs, Highly similar to c-Jun leucine zipper interactive); AA874995 (UI-R-E0-cf-d-08-0-UI.s1 cDNA); AA875004 (UI-R-E0-cb-b-07-0-UI.s1 cDNA); AA875019 (UI-R-E0-cb-f-08-0-UI.s1 cDNA); AA875032 (UI-R-E0-cb-h-09-0-UI.s1 cDNA); AA875037 (UI-R-E0-cb-a-03-0-UI.s1 cDNA); AA875129 (UI-R-E0-bu-e-01-0-UI.s2 cDNA); AA875257 (UI-R-E0-cq-d-12-0-UI.s1 cDNA); AA891037 (EST194840 cDNA); AA891041 (un B proto-oncogene); AA891221 (EST195024 cDNA); AA891476 (EST195279 cDNA); AA891537 (EST195340 cDNA); AA891690 (ESTs, Weakly similar to p-serine aminotransferase); AA891727 (EST195530 cDNA); AA891734 (EST195537 cDNA); AA891774 (EST195577 cDNA); AA891810 (EST195613 cDNA); AA891880 (Loc65042); AA891916 (membrane interacting protein of RGS16); AA891944 (EST195747 cDNA); AA891950 (EST195753 cDNA); AA892146 (EST195949 cDNA); AA892298 (EST196101 cDNA); AA892414 (EST196217 cDNA); AA892511 (EST196314 cDNA); AA892520 (EST196323 cDNA); AA892538 (EST196341 cDNA); AA892637 (EST196440 cDNA); AA892775 (Lysozyme); AA892796 (EST196599 cDNA); AA892813 (EST196616 cDNA); AA892986 (EST196789 cDNA); AA893082 (EST196885 cDNA); AA893185 (EST196988 cDNA); AA893199 (EST197002 cDNA); AA893224 (EST197027 cDNA); AA893320 (EST197123 cDNA); AA893584 (EST197387 cDNA); AA893690 (EST197493 cDNA); AA893708 (KIAA0560); AA893717 (EST197520 cDNA); AA893743 (EST197546 cDNA); AA893788 (ESTs, Highly similar to chromobox protein homolog 5); AA893946 (EST197749 cDNA); AA894305 (EST198108 cDNA); AA924925 (ER transmembrane protein Dri 42); AA942685 (cytosolic cysteine dioxygenase 1); AA955306 (ras-related protein rab10); AA955388 (Na+K+ transporting ATPase 2, beta polypeptide 2); AA957132 (N-acetylglucosaminyltransferase I); AB000778 (Phoshpolipase D gene 1); AB006451 (Tim23); AB008538 (HB2); AB016532 (period homolog 2 (Drosophila)); AF000899 (p58/p45, nucleolin); AF007554 (Mucin1); AF007758 (synuclein, alpha); AF007890 (resection-induced TPI (rs11)); AF008554 (implantation-associated protein (IAG2)); AF013144 (MAP-kinase phosphatase (cpg21)); AF016269 (kallikrein 6 (neurosin, zyme)); AF016296 (neuropilin); AF019974 (Chromogranin B, parathyroid secretory protein); AF020046 (integrin alpha E1, epithelial-associated); AF021935 (Ser-Thr protein kinase); AF023087 (Early growth response 1); AF030050 (replication factor C); AF030088 (RuvB-like protein 1); AF040954 (putative protein phosphatase 1 nuclear targeting subunit); AF051561 (solute carrier family 12, member 2); AF055477 (L-type voltage-dependent Ca2+ channel (?1D subunit)); AF074608 (MHC class I antigen (RT1.EC2) gene); AF076183 (cytosolic sorting protein PACS-1a (PACS-1)); AF095927 (protein phosphatase 2C); AI010110 (SH3-domain GRB2-like 1); AI012589 (glutathione S-transferase, pi 2); AI013627 (defender against cell death 1); AI013861 (3-hydroxyisobutyrate dehydrogenase); AI045249 (heat shock 70 kD protein 8); AI102031 (myc box dependent interacting protein 1); AI102299 (Bid3); AI102839 (cerebellar Ca-binding protein, spot 35 protein); AI102868 (EST212157 cDNA); AI104388 (heat shock 27 kD protein 1); AI136891 (zinc finger protein 36, C3H type-like 1); AI168942 (branched chain keto acid dehydrogenase E1); AI169265 (Atp6s1); AI171966 (ESTs, Highly similar to selenide, water dikinase 2); AI175973 (ESTs, Highly similar to NADH dehydrogenase); AI176491 (EST220076 cDNA); AI176595 (Cathepsin L); AI176621 (iron-responsive element-binding protein); AI177404 (EST221024 cDNA); AI178204 (EST221869 cDNA); AI178921 (Insulin degrading enzyme); AI228548 (ESTs, Highly similar to DKFZp586G0322.1); AI230247 (selenoprotein P, plasma, 1); AI230778 (ESTs, Highly similar to protein-tyrosine sulfotrans. 2); AI230914 (farnesyltransferase beta subunit); AI231807 (ferritin light chain 1); AI231807 (ferritin light chain 1); AI232268 (LDL receptor-related protein associated protein 1); AI235344 (geranylgeranyltransferase type I (GGTase-I)); AI237007 (ESTs, Highly similar to flavoprot.-ubiquin. Oxidoreduct.); AI638971 (mixed-tissue library cDNA clone rx04989 3); AI638997 (mixed-tissue library cDNA clone rx05048 3); AI639151 (mixed-tissue library cDNA clone rx02802 3); AI639209 (mixed-tissue library cDNA clone rx00680 3); AI639257 (mixed-tissue library cDNA clone rx01119 3); AI639477 (mixed-tissue library cDNA clone rx02351 3); D00569 (2,4-dienoyl CoA reductase 1, mitochondrial); D10262 (choline kinase); D10699 (ubiquitin carboxy-terminal hydrolase L1); D10854 (aldehyde reductase); D10874 (lysosomal vacuolar proton pump (16 kDa)); D16478 (mitochondrial long-chain enoyl-CoA hydratase); D17309 (delta 4-3-ketosteroid-5-beta-reductase); D28110 (myelin-associated oligodendrocytic basic protein); D28110 (myelin-associated oligodendrocytic basic protein); D28557 (cold shock domain protein A); D29766 (v-crk-associated tyrosine kinase substrate); D37951 (MIBP1 (c-myc intron binding protein 1)); D45247 (proteasome subunit RCX); D45249 (protease (prosome, macropain) 28 subunit, alpha); D78308 (calreticulin); D83948 (adult liver S1-1 protein); D88586 (eosinophil cationic protein); D89340 (dipeptidylpeptidase III); D89730 (Fibulin 3, fibulin-like extracellular matrix protein 1); D90211 (Lysosomal-associated membrane protein 2); E03229 (cytosolic cysteine dioxygenase 1); H33086 (EST108750 cDNA); H33725 (associated molecule with the SH3 domain of STAM); J02752 (acyl-coA oxidase); J02773 (heart fatty acid binding protein); J05022 (peptidylarginine deiminase); J05031 (Isovaleryl Coenzyme A dehydrogenase); J05132 (UDP-glucuronosyltransferase); K02248 (Somatostatin); L00191 (Fibronectin 1); L13202 (RATHFH2 HNF-3/fork-head homolog-2 (HFH-2)); L19998 (sulfotransferase family 1A, phenol-preferring, member 1); L24896 (glutathione peroxidase 4); L25605 (Dynamin 2); L26292 (Kruppel-like factor 4 (gut)); L29573 (neurotransmitter transporter, noradrenalin); L42855 (transcription elongation factor B (SIII) polypeptide 2); M10068 (NADPH-cytochrome P-450 oxidoreductase); M13100 (long interspersed repetitive DNA sequence LINE3); M13100 (long interspersed repetitive DNA sequence LINE3); M19936 (Prosaposin-sphingolipid hydrolase activator); M23601 (Monoamine oxidase B); M24104 (synaptobrevin 2); M24104 (Vesicle-associated membrane protein (synaptobrevin 2)); M24852 (Neuron specific protein PEP-19 (Purkinje cell protein 4)); M31174 (thyroid hormone receptor alpha); M36453 (Inhibin, alpha); M55015 (nucleolin); M57276 (Leukocyte antigen (Ox-44)); M58404 (thymosin, beta 10); M80550 (adenylyl cyclase); M83745 (Protein convertase subtilisin/kexin, type I); M89646 (ribosomal protein S24); M91234 (VL30 element); M93273 (somatostatin receptor subtype 2); M93669 (Secretogranin II); M99485 (Myelin oligodendrocyte glycoprotein); S53527 (S100 calcium-binding protein, beta (neural)); S61868 (Ryudocan/syndecan 4); S72594 (tissue inhibitor of metalloproteinase 2); S77492 (Bone morphogenetic protein 3); S77858 (non-muscle myosin alkali light chain); U04738 (Somatostatin receptor subtype 4); U07619 (Coagulation factor III (thromboplastin, tissue factor)); U08259 (Glutamate receptor, N-methyl D-aspartate 2C); U10357 (pyruvate dehydrogenase kinase 2 subunit p45 (PDK2)); U14950 (tumor suppressor homolog (synapse associ. protein)); U17254 (immediate early gene transcription factor NGFI-B); U17254 (immediate early gene transcription factor NGFI-B); U18771 (Ras-related protein Rab-26); U27518 (UDP-glucuronosyltransferase); U28938 (receptor-type protein tyrosine phosphatase D30); U38379 (Gamma-glutamyl hydrolase); U38801 (DNA polymerase beta); U67080 (r-MyT13); U67136 (A kinase (PRKA) anchor protein 5); U67137 (guanylate kinase associated protein); U72620 (Lot1); U75405 (procollagen, type I, alpha 1); U75917 (clathrin-associated protein 17); U77777 (interleukin 18); U78517 (cAMP-regulated guanine nucleotide exchange factor II); U89905 (alpha-methylacyl-CoA racemase); V01244 (Prolactin); X02904 (glutathione S-transferase P subunit); X05472 (2.4 kb repeat DNA right terminal region); X06769 (FBJ v-fos oncogene homolog); X06916 (S100 calcium-binding protein A4); X13905 (ras-related rab1B protein); X14323 (Fc receptor, IgG, alpha chain transporter); X16933 (RNA binding protein p45AUF1); X53427 (glycogen synthase kinase 3 alpha (EC 2.7.1.37)); X53504 (ribosomal protein L12); X54467 (cathepsin D); X55153 (ribosomal protein P2); X57281 (Glycine receptor alpha 2 subunit); X58294 (carbonic anhydrase 2); X59737 (ubiquitous mitochondrial creatine kinase); X60212 (ASI homolog of bacterial ribosomal subunit protein L22); X62950 (pBUS30 with repetitive elements); X67805 (Synaptonemal complex protein 1); X72757 (cox VIa gene (liver)); X74226 (LL5 protein); X76489 (CD9 cell surface glycoprotein); X76985 (latexin); X82445 (nuclear distribution gene C homolog (Aspergillus)); X84039 (lumican); X89696 (TPCR06 protein); X97443 (integral membrane protein Tmp21-I (p23)); X98399 (solute carrier family 14, member 1); Y17295 (thiol-specific antioxidant protein (1-Cys peroxiredoxin)); Z48225 (protein synthesis initiation factor eIF-2B delta subunit); Z49858 (plasmolipin).
Using the method of the invention, we have also identified a set of genes and ESTs that changed with age (p≦0.05), but which are correlated with cognitive performance in behavioral tests. These include L03294 (Lpl, lipoprotein lipase); M18416 (Egr1, Early growth response 1 (Krox-24)); S68245 (Ca4, carbonic anhydrase 4); M64780 (Agrn, Agrin); M27207 (Colla1, Procollagen-type I (alpha 1)); X16554 (Prps1, Phosphoribosyl pyrophosphate synthetase 1); M92433 (NGFI-C, Zinc-finger transcription factor (early response gene)); AA859975 (LOC64201, 2-oxoglutarate carrier); L08595 (Nuclear receptor subfamily 4, group A, member 2); M24542 (RISP, Rieske iron-sulfur protein); AI030089 (Nopp130, nucleolar phosphoprotein p130); AF104362 (Omd, Osteomodulin (osteoadherin)); L46873 (Slc15a1, Oligopeptide transporter); AI176689 (MAPKK 6, mitogen-activated protein kinase kinase 6); U66470 (rCGR11, Cell growth regulator); AF016387 (RXRG, retinoid X-receptor gamma); M18467 (Got2, glutamate oxaloacetate transaminase 2); X54793 (Hsp60, heat shock protein 60); X64401 (Cyp3a3, Cytochrome P450-subfamily 111A (polypeptide 3)); M37584 (H2afz, H2A histone family (member Z)); L21192 (GAP-43, membrane attached signal protein 2 (brain)); AA875047 (TCPZ, T-complex protein 1 (zeta subunit)); U90610 (Cxcr4, CXC chemokine receptor); AF003904 (CRH-binding protein); U83880 (GPDH-M, glycerol-3-phosphate dehydrogenase, mitochondrial); X89703 (TPCR19, Testis Polymerase Chain Reaction product 19); D63886 (MMP16, matrix metalloproteinase 16); J05499 (GLS, glutaminase (mitochondrial)); D21799 (Psmb2, Proteasome subunit (beta type 2)); AA800794 (HT2A, zinc-finger protein); U90887 (Arg2, arginase type II); S82649 (Narp, neuronal activity-regulated pentraxin); M74223 (VGF, neurosecretory protein); AA874794 (Bex3, brain expressed X-linked 3); M15191 (Tac1, Tachykinin); AA892506 (coronin, actin binding protein 1A); L04485 (MAPPK1, mitogen-activated protein kinase kinase 1); AA799641 (S164, Contains a PWI domain associated with RNA splicing); AA817892 (Gnb2, Guanine nucleotide binding protein (beta 2 subunit)); AA893939 (DSS1, deleted in split hand/split foot protein 1); AF000901 (P58/P45, Nucleoporin p58); AF087037 (Btg3, B-cell translocation gene 3); AB000280 (PHT1, peptide/histidine transporter); M87854 (Beta-ARK-1, beta adrenergic receptor kinase 1); U06099 (Prdx2, Peroxiredoxin 2); AF058795 (Gb2, GABA-B receptor); AA800517 (VAP1, vesicle associated protein); U63740 (Fez1, Protein kinase C-binding protein Zeta1); U53922 (Hsj2, DnaJ-like protein (RDJ1)); U78102 (Egr2, Early growth response 2); U44948 (SmLIM, smooth muscle cell LIM protein); U87627 (MCT3, putative monocarboxylate transporter); AB020504 (PMF31, highly homologus to mouse F-box-WD40 repeat protein 6); M21354 (Col3a1, collagen type III alpha-1); AA893664 (Temo, sertoli cell marker (KIAA0077 protein fragment)); AB010437 (CDH8, Cadherin-8); M22756 (Ndufv2, mitochondrial NADH dehydrogenase (24 kDa)); AA799389 (Rab3B, ras-related protein); AI172476 (Tieg-1, TGF-beta-inducible early growth response protein 1); AF091563 (Olfactory receptor); M64376 (Olfactory protein); J04488 (Ptgds, Prostaglandin D synthase); X71127 (c1qb, complement component 1-q (beta polypeptide)); J03752 (Microsomal GST-1, glutathione S-transferase); J03481 (Qdpr, Dihydropteridine reductase); L40362 (MHC class I RT1.C-type protein); M94918 (Hbb, beta hemoglobin); M55534 (Cryab, alpha crystallin polypeptide 2); U17919 (Aif1, allograft inflammatory factor 1); M15562 (MHC class II RT1.u-D-alpha chain); AA799645 (Phospholemman, FXYD domain-containing ion transport regulator 1); X13044 (Cd74, CD74 antigen); M24324 (RTS, MHC class I RT1 (RTS) (u haplotype)); U31866 (Nclone10); M32062 (Fcgr3, Fc IgG receptor III (low affinity)); AF095741 (Mg87); L03201 (Ctss, cathepsin S); M27905 (Rpl21, Ribosomal protein L21); D38380 (Tf, Transferrin); AA893493 (RPL26, Ribosomal protein L26); AJ222813 (I118, interleukin 18); E13541 (Cspg5, chondroitin sulfate proteoglycan 5); X54096 (Lcat, Lecithin-cholesterol acyltransferase); L40364 (RT1Aw2, RT1 class Ib); D28111 (MOBP, myelin-associated oligodendrocytic basic protein); M32016 (Lamp2, lysosomal-associated membrane protein 2); X13167 (NF1-A, nuclear factor 1 A); U26356 (S100A1, S100 protein (alpha chain)); AI231213 (Kangai 1, suppression of tumorigenicity 6); AI170268 (Ptgfr, Prostaglandin F receptor); X62952 (Vim, vimentin); AI014169 (Vdup1, vitamin D-upregulated); AA850219 (Anx3, Annexin A3); D84477 (Rhoa, ras-related homolog A2); X52477 (C3, Complement component 3); X52619 (Rpl28, Ribosomal protein L28); X06554 (S-MAG, myelin-associated glycoprotein C-term); Z50144 (Kat2, kynurenine aminotransferase II); X14181 (RPL18A, Ribosomal protein L18a); AA892333 (Tuba1, alpha-tubulin); U67082 (KZF-1, Kruppel associated box (KRAB) zinc finger 1); U11760 (Vcp, valosin-containing protein); AF048828 (VDAC1, voltage-dependent anion channel 1); M31076 (TNF-alpha, Transforming growth factor (alpha)); S83279 (HSDIV, 17-beta-hydroxysteroid dehydrogenase type IV); AI102103 (Pik4cb, Phosphatidylinositol 4-kinase); X56325 (Hba1, alpha 1 hemoglobin); X73371 (FCGR2, Low affinity immunoglobulin gamma Fc receptor II); X78848 (Gsta1, Glutathione-S-transferase (alpha type)); U92564 (Roaz, Olf-1/EBF associated Zn finger protein); AI171462 (Cd24, CD24 antigen); X83231 (PAIHC3, Pre-alpha-inhibitor, heavy chain 3); AF097593 (Ca4, cadherin 2-type 1 (neuronal)); X68283 (Rpl29, Ribosomal protein L29); S55427 (Pmp, peripheral myelin protein); AA818025 (Cd59, CD59 antigen); E01534 (Rps15, Ribosomal protein S15); U37138 (Sts, Steroid sulfatase); X55572 (Apod, Apolipoprotein D); AI028975 (AP-1, adaptor protein complex (beta 1)); L16995 (ADD1, adipocyte determination/differentiation-dependent factor 1); U07971 (Transamidinase, Glycine amidinotransferase, mitochondrial); L07736 (Cpt1a, Carnitine palmitoyltransferase 1 alpha (liver)); AI237535 (LitaF, LPS-induced TNF-alpha factor); AI175486 (Rps7, Ribosomal protein S7); U32498 (RSEC8, rat homolog of yeast sec8); X53504 (RPL12, Ribosomal protein L12); AF023621 (Sort1, sortilin); AF083269 (P41-Arc, actin-related protein complex 1b); AA891810 (GST, Glutathione S-transferase); M77694 (Fah, fimarylacetoacetate hydrolase); M22357 (MAG, myelin-associated glycoprotein); AI230712 (Pace4, Subtilisin-like endoprotease); AF008439 (NRAMP2, Natural resistance-associated macrophage protein 2); U77829 (Gas-5, growth arrest homolog); U92081 (Gp38, Glycoprotein 38); AA891445 (Skd3, suppressor of K+ transport defect 3); AI177161 (Nfe212, NF-E2-related factor 2); AF031430 (Stx7, Syntaxin 7); L35921 (Ggamma, GTP-binding protein (gamma subunit)); X62322 (Gm, Granulin); AF028784 (GFAP, glial fibrillary acidic protein); and AI234146 (Csrp1, Cysteine rich protein 1).
Using the method of the invention, we have further identified a set of genes and ESTs that changed with age (p≦0.01). These include AA891651 (rc_AA891651 EST195454 cDNA); AI070108 (rc_AI070108 UI-R-Y0-lu-a-09-0-UI.s1 cDNA); AI176689 (mitogen-activated protein kinase kinase 6); AI012051 (rc_AI012051 EST206502 cDNA); AI233365 (rc_AI233365 EST230053 cDNA); AA892532 (rc_AA892532 EST196335 cDNA); AA893185 (rc_AA893185 EST196988 cDNA); AA964320 (rc_AA964320 UI-R-C0-gu-e-09-0-UI.s1 cDNA); AA963449 (rc_AA963449 UI-R-E1-gj-e-08-0-UI.s1 cDNA); AA859632 (rc_AA859632 UI-R-E0-bs-h-08-0-UI.s1 cDNA); AI169265 (Atp6s1); AA850781 (rc_AA850781 EST193549 cDNA); AJ222813 (interleukin 18); D38380 (Transferrin); J03481 (dihydropteridine reductase); M24542 (Rieske iron-sulfur protein); L03294 (Lipoprotein lipase); L19998 (sulfotransferase family 1A, phenol-preferring, member 1); U53922 (DnaJ-like protein (RDJ1)); X54793 (liver heat shock protein (hsp60)); X62952 (vimentin); M55534 (Crystallin, alpha polypeptide 2); J03752 (microsomal glutathione S-transferase 1); X64401 (Cytochrome P450, subfamily 111A, polypeptide 3); X78848 (Gsta1); AF016387 (retinoid X receptor gamma); AF031430 (syntaxin 7); AF051561 (solute carrier family 12, member 2); AF076183 (cytosolic sorting protein PACS-1a (PACS-1)); AF095576 (adaptor protein with pleckstrin homology and src homology 2 domains); AF095741 (MG87); AF097593 (cadherin 2, type 1, N-cadherin (neuronal)); AF104362 (osteoadherin); D10699 (ubiquitin carboxy-terminal hydrolase L1); D28111 (myelin-associated oligodendrocytic basic protein); D37951 (MIBP1 (c-myc intron binding protein 1)); D84477 (RhoA); L13202 (RATHFH2 HNF-3/fork-head homolog-2 (HFH-2)); L26292 (Kruppel-like factor 4 (gut)); L46873 (solute carrier family 15 (oligopeptide transporter), member 1); M13100 (RATLIN3A long interspersed repetitive DNA sequence LINE3 (L1Rn)); M27207 (procollagen, type I, alpha 1); M92433 (Zinc-finger transcription factor NGFI-C (early response gene)); M94918 (Hemoglobin, beta); M94919 (Hemoglobin, beta); S55427 (Peripheral myelin protein); S68245 (carbonic anhydrase 4); S82649 (Narp=neuronal activity-regulated pentraxin); U10894 (allograft inflammatory factor 1); U26356 (RNSHUNA1 S100A1 gene); U75397 (RNKROX2 Krox-24); U75405 (procollagen, type I, alpha 1); U77829 (RNU77829 gas-5 growth arrest homolog non-translated sequence); U92081 (glycoprotein 38); X06554 (RNMAGSR myelin-associated glycoprotein (S-MAG) C-term); X13167 (Nuclear Factor IA); X14181 (RRRPL18A ribosomal protein L18a); X56325 (Hemoglobin, alpha 1); X60351 (Crystallin, alpha polypeptide 2); E13541 (chondroitin sulfate proteoglycan 5); M22357 (1B236/myelin-associated glycoprotein (MAG)); M24026 (RT1 class Ib gene); M24324 (MHC class I RT1 (RTS) (u haplotype)); J04488 (Prostaglandin D synthase); M15191 (Tachykinin (substance P, neurokinin A, neuropeptide K, neuropeptide gamma)); M74223 (VGF); U17254 (immediate early gene transcription factor NGFI-B); U08259 (Glutamate receptor, ionotropic, N-methyl D-aspartate 2C); U19866 (activity regulated cytoskeletal-associated protein); L40364 (RT1 class Ib gene); U17919 (allograft inflammatory factor 1); U78102 (early growth response 2); U67082 (KRAB-zinc finger protein KZF-1); U77777 (interleukin 18); D78018 (Nuclear Factor IA); U92564 (Olf-1/EBF associated Zn finger protein Roaz); AF008439 (Solute carrier family 11 member 2 (natural resistance-associated macrophage protein 2)); AB003726 (RuvB-like protein 1); M83561 (Glutamate receptor, ionotropic, kainate 1); AI639151 (mixed-tissue library cDNA clone rx02802 3); AI639247 (mixed-tissue library cDNA clone rx03939 3); AI639381 (mixed-tissue library cDNA clone rx01495 3); AI639532 (mixed-tissue library cDNA clone rx01030 3); AA799645 (FXYD domain-containing ion transport regulator 1); AA900516 (Pdi2); AI014169 (Vdup1); AI030089 (Nopp140); AI102299 (Bid3); AA818025 (CD59 antigen); AI170268 (Prostaglandin F receptor); AI171462 (CD24 antigen); AI171966 (ESTs, Highly similar to SPS2 MOUSE SELENIDE, WATER DIKINASE 2 [M. musculus]); AI76456 (ESTs, Weakly similar to ABP2_HUMAN ENDOTHELIAL ACTIN-BINDING PROTEIN [H. sapiens]); AI177161 (NF-E2-related factor 2); AI179576 (Hemoglobin, beta); AI230712 (Subtilisin-like endoprotease); AI230914 (farnesyltransferase beta subunit); AI231213 (kangai 1 (suppression of tumorigenicity 6), prostate); AI237731 (Lipoprotein lipase); M83745 (Protein convertase subtilisin/kexin, type I); M27905 (ribosomal protein L21); M32016 (Lysosomal-associated membrane protein 2); M11071 (RT1 class Ib gene); M15562 (MHC class II RT1.u-D-alpha chain); M15880 (Neuropeptide Y); L08595 (nuclear receptor subfamily 4, group A, member 2); M18416 (Early growth response 1); L40362 (MHC class I RT1.C-type protein); Z50144 (kynurenine/alpha-aminoadipate aminotransferase); X71127 (complement component 1, q subcomponent, beta polypeptide); U44948 (smooth muscle cell LIM protein (SmLIM)); AA850219 (Annexin A3); X73371 (FCGR2); X57281 (Glycine receptor alpha 2 subunit (glycine receptor, neonatal)); X83231 (pre-alpha-inhibitor); X52477 (Complement component 3); X16554 (phosphoribosyl pyrophosphate synthetase 1); X78605 ((Sprague Dawley) rab4b ras-homologous GTPase); X82445 (nuclear distribution gene C homolog (Aspergillus)); X52619 (ribosomal protein L28); X68283 (ribosomal protein L29); X13044 (CD74 antigen (invariant polpypeptide of major histocompatibility class II antigen-associated)); X54096 (Lecithin-cholesterol acyltransferase); U31866 (Nclone10); U72620 (Lot1); U66470 (rCGR11); M31018 (RT1 class Ib gene); U90887 (arginase type II); M18467 (Glutamate oxaloacetate transaminase 2, mitochondrial (aspartate aminotransferase 2)); M64780 (Agrin); U87627 (putative monocarboxylate transporter (MCT3)); AF019974 (Chromogranin B, parathyroid secretory protein); L03201 (cathepsin S); AB008538 (HB2); D89340 (dipeptidylpeptidase III); M77694 (fumarylacetoacetate hydrolase); M32062 (Fc-gamma receptor); L21192 (brain abundant, membrane attached signal protein 2); M37584 (H2afz); AA858588 (ESTs, Weakly similar to ODP2 RAT DIHYDROLIPOAMIDE ACETYLTRANSFERASE COMPONENT OF PYRUVATE DEHYDROGENASE COMPLEX [R. norvegicus]); AA858617 (rc_AA858617 UI-R-E0-bq-b-06-0-UI.s1 cDNA); AA859562 (rc_AA859562 UI-R-E0-by-b-03-0-UI.s1 cDNA); AA859626 (rc_AA859626 UI-R-E0-bs-h-02-0-UI.s1 cDNA); AA859690 (rc_AA859690 UI-R-E0-bx-e-11-0-UI.s1 cDNA); AA859777 (rc_AA859777 UI-R-E0-bu-e-10-0-UI.s1 cDNA); AA859975 (LOC64201); AA860030 (UI-R-E0-bz-e-07-0-UI.s2 cDNA); AA866291 (rc_AA866291 UI-R-A0-ac-e-12-0-UI.s3 cDNA); AA866409 (rc_AA866409 UI-R-E0-ch-a-03-0-UI.s1 cDNA); AA866411 (NDN); AA874794 (Bex3); A874887 (rc_AA874887 UI-R-E0-ci-g-10-0-UI.s1 cDNA); AA875004 (rc_AA875004 UI-R-E0-cb-b-07-0-UI.s1 cDNA); AA875037 (rc_AA875037 UI-R-E0-cb-a-03-0-UI.s1 cDNA); AA875047 (TCPZ); AA875059 (rc_AA875059 UI-R-E0-cb-f-04-0-UI.s1 cDNA); AA875129 (rc_AA875129 UI-R-E0-bu-e-01-0-UI.s2 cDNA); H31418 (rc_H31418 EST105434 cDNA); H31665 (rc_H31665 EST105952 cDNA); H32977 (rc_H32977 EST108553 cDNA); H33725 (associated molecule with the SH3 domain of STAM); AA891037 (rc_AA891037 EST194840 cDNA); AA891445 (Skd3); AA891690 (ESTs, Weakly similar to SERC_HUMAN PHOSPHOSERINE AMINOTRANSFERASE [H. sapiens]); AA891717 (USF1); AA891734 (rc_AA891734 EST195537 cDNA); AA891785 (rc_AA891785 EST195588 cDNA); AA891810 (ESTs, Highly similar to GTK1 RAT GLUTATHIONE S-TRANSFERASE, MITOCHONDRIAL [R. norvegicus]); AA891965 (rc_AA891965 EST195768 cDNA); AA892333 (Tuba1); AA892353 (ESTs, Moderately similar to JC5823 NADH dehydrogenase [H. sapiens]); AA892511 (rc_AA892511 EST196314 cDNA); AA892986 (rc_AA892986 EST196789 cDNA); AA893032 (ESTs, Moderately similar to CALX RAT CALNEXIN PRECURSOR [R. norvegicus]); AA893082 (rc_AA893082 EST196885 cDNA); AA893493 (RPL26); AA893607 (rc_AA893607 EST197410 cDNA); AA893708 (KIAA0560); AA893743 (rc_AA893743 EST197546 cDNA); AA894104 (rc_AA894104 EST197907 cDNA); AA799449 (EST, Weakly similar to UBP4 MOUSE UBIQUITIN CARBOXYL-TERMINAL HYDROLASE 4 [M. musculus]); AA799779 (acyl-CoA:dihydroxyacetonephosphate acyltransferase); AA799803 (ESTs, Weakly similar to K1CU RAT KERATIN, TYPE I CYTOSKELETAL 21 [R. norvegicus]); AA799854 (rc_AA799854 EST189351 cDNA); AA799996 (rc_AA799996 EST189493 cDNA); AA800693 (rc_AA800693 EST190190 cDNA); AA800708 (ESTs, Weakly similar to S28312 hypothetical protein F02A9.4-Caenorhabditis elegans [C. elegans]); AA800794 (HT2A); and AA800948 (Tuba4).
We have also identified age-related ESTs, including AA963449 (UI-R-E1-gj-e-08-0-UI.s1 cDNA); AA892532 (EST196335 cDNA); AA859626 (UI-R-E0-bs-h-02-0-UI.s1 cDNA); AA893743 (EST197546 cDNA); AI233365 (EST230053 cDNA); H31665 (EST105952 cDNA); AA892353 (ESTs, Moderately similar to JC5823 NADH dehydrogenase); AI639247 (mixed-tissue library cDNA clone rx03939 3); AA858617 (UI-R-E0-bq-b-06-0-UI.s1 cDNA); AI639429 (mixed-tissue library cDNA clone rx00973 3); AA858620 (UI-R-E0-bq-b-09-0-UI.s1 cDNA); AA866291 (UI-R-A0-ac-e-12-0-UI.s3 cDNA); AA894104 (EST197907 cDNA); AA799996 (EST189493 cDNA); AA892805 (EST196608 cDNA); AI639019 (mixed-tissue library cDNA clone rx01107 3); AA799538 (EST189035 cDNA); AI070108 (UI-R-Y0-lu-a-09-0-UI.s1 cDNA); AA866409 (UI-R-E0-ch-a-03-0-UI.s1 cDNA); AA859632 (UI-R-E0-bs-h-08-0-UI.s1 cDNA); AA891651 (EST195454 cDNA); AA893032 (ESTs, Moderately similar to CALX calnexin precursor); AA891965 (EST195768 cDNA); AA800708 (ESTs, Weakly similar to S28312 hypothetical protein F02A9.4); AA964320 (UI-R-C0-gu-e-09-0-UI.s1 cDNA); AA893173 (EST196976 cDNA); H32977 (EST108553 cDNA); AA874887 (UI-R-E0-ci-g-10-0-UI.s1 cDNA); AA850781 (EST193549 cDNA); AI176456 (ESTs, Weakly similar to endothelial actin-binding protein); H31418 (EST105434 cDNA); AA858588 (ESTs, Weakly similar to ODP2 dihydrolipoamide acetyl transferase); AA891785 (EST195588 cDNA); AA799803 (ESTs, Weakly similar to K1CU cytoskeletal keratin (type 1)); AA799449 (EST, Weakly similar to UBP4 ubiquitin carboxyl-terminal hydrolase 4); AA859777 (UL-R-E0-bu-e-10-0-UI.s1 cDNA); AI639532 (mixed-tissue library cDNA clone rx01030 3); AA875059 (UI-R-E0-cb-f-04-0-UI.s1 cDNA); AI012051 (EST206502 cDNA); AA800549 (EST190046 cDNA); AA799854 (EST189351 cDNA); and AA892520 (EST196323 cDNA).
Those of skill in the genomics art will understand that the identified genes and ESTs have utility as biomarkers of brain aging. Those of skill in the genomics art will understand that the mammalian homologues (including rat, mouse and human homologues) the identified genes and ESTs are also as biomarkers of brain aging. The easiest method for identifying mammalian homologues of the identified genes and ESTs is by identifying the homologues in the GenBank database, preferably, or in the SwissProtein and the Genome Ontology databases. Additional guidance as to homology can be obtained by using commercially available computer programs, such as DNA Strider and Wisconsin GCG, and following the instructions for the determination of the degree of homolgy between selected polynucleotides.
The foregoing description has been presented only for the purposes of illustration and is not intended to limit the invention to the precise form disclosed, but by the claims appended hereto.
This invention has been made in part with government support under grants AG04542, AG10836, AG18228 and AG14979 from the National Institute on Aging, and by MH59891. The government of the United States of America may have certain rights in this invention.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US02/25607 | 8/13/2002 | WO |
Number | Date | Country | |
---|---|---|---|
60311343 | Aug 2001 | US |