Epigenetic predictor of age

Sven Bocklandt, Wen Lin, Mary E Sehl, Francisco J Sánchez, Janet S Sinsheimer, Steve Horvath, Eric Vilain, Sven Bocklandt, Wen Lin, Mary E Sehl, Francisco J Sánchez, Janet S Sinsheimer, Steve Horvath, Eric Vilain

Abstract

From the moment of conception, we begin to age. A decay of cellular structures, gene regulation, and DNA sequence ages cells and organisms. DNA methylation patterns change with increasing age and contribute to age related disease. Here we identify 88 sites in or near 80 genes for which the degree of cytosine methylation is significantly correlated with age in saliva of 34 male identical twin pairs between 21 and 55 years of age. Furthermore, we validated sites in the promoters of three genes and replicated our results in a general population sample of 31 males and 29 females between 18 and 70 years of age. The methylation of three sites--in the promoters of the EDARADD, TOM1L1, and NPTX2 genes--is linear with age over a range of five decades. Using just two cytosines from these loci, we built a regression model that explained 73% of the variance in age, and is able to predict the age of an individual with an average accuracy of 5.2 years. In forensic science, such a model could estimate the age of a person, based on a biological sample alone. Furthermore, a measurement of relevant sites in the genome could be a tool in routine medical screening to predict the risk of age-related diseases and to tailor interventions based on the epigenetic bio-age instead of the chronological age.

Conflict of interest statement

Competing Interests: The authors have declared that no competing interests exist.

Figures

Figure 1. Detection of gene co-methylation modules…
Figure 1. Detection of gene co-methylation modules in human saliva in twins.
(a) Branches of the hierarchical cluster tree define five co-methylation modules which are assigned a color as can be seen from the first color band underneath the tree. Probes that could not be clustered into one of these modules were coded grey. Every probe represents a line in the hierarchical cluster tree. Distance between two probes is shown as height on the y-axis. The second color band encodes the age relationships of each gene. Genes with positive age correlations are colored in blue. (b) Barplots showing age relationships of modules. Specifically, the y-axis shows the mean Student T-statistic testing whether the methylation status of a probe is correlated with age. Note that the green module is enriched for probes that have a significant positive correlation with age. A t-statistic value of 2 or higher indicates a significant correlation (p

Figure 2. Percentage methylation versus age for…

Figure 2. Percentage methylation versus age for three markers validated in three sample sets.

Original…

Figure 2. Percentage methylation versus age for three markers validated in three sample sets.
Original twin samples are blue, male control samples are red, female control samples green. Linear trendlines are shown in the colors of the individual sample sets a) Edaradd r = −0.81 (twins), r = −0.73 (male controls), r = −0.75 (female controls) b) NPTX2 r = 0.52 (twins), r = 0.79 (male controls), r = 0.03 (female controls) c) Tom1L1 r = −0.70 (twins), r = −0.49 (male controls), r = −0.24 (female controls).

Figure 3. Predicted versus observed age of…

Figure 3. Predicted versus observed age of all subjects using a leave-one-out model.

A multivariate…

Figure 3. Predicted versus observed age of all subjects using a leave-one-out model.
A multivariate regression model was fit on all but one sample and its predicted age (y-axis) was related to the truly observed age of the left out sample (x-axis). The predicted values are highly correlated with the observed ages (r = 0.83, p = 2.2×10−16, n = 66), and the average absolute difference between the predicted and the observed age is 5.2 years.
Figure 2. Percentage methylation versus age for…
Figure 2. Percentage methylation versus age for three markers validated in three sample sets.
Original twin samples are blue, male control samples are red, female control samples green. Linear trendlines are shown in the colors of the individual sample sets a) Edaradd r = −0.81 (twins), r = −0.73 (male controls), r = −0.75 (female controls) b) NPTX2 r = 0.52 (twins), r = 0.79 (male controls), r = 0.03 (female controls) c) Tom1L1 r = −0.70 (twins), r = −0.49 (male controls), r = −0.24 (female controls).
Figure 3. Predicted versus observed age of…
Figure 3. Predicted versus observed age of all subjects using a leave-one-out model.
A multivariate regression model was fit on all but one sample and its predicted age (y-axis) was related to the truly observed age of the left out sample (x-axis). The predicted values are highly correlated with the observed ages (r = 0.83, p = 2.2×10−16, n = 66), and the average absolute difference between the predicted and the observed age is 5.2 years.

References

    1. Goyns MH. Genes, telomeres and mammalian ageing. Mech Ageing Dev. 2002;123:791–799.
    1. Boks MP, Derks EM, Weisenberger DJ, Strengman E, Janson E, et al. The relationship of DNA methylation with age, gender and genotype in twins and healthy controls. PLoS One. 2009;4:e6767.
    1. Fraga MF, Ballestar E, Paz MF, Ropero S, Setien F, et al. Epigenetic differences arise during the lifetime of monozygotic twins. Proc Natl Acad Sci U S A. 2005;102:10604–10609.
    1. Bork S, Pfister S, Witt H, Horn P, Korn B, et al. DNA methylation pattern changes upon long-term culture and aging of human mesenchymal stromal cells. Aging Cell. 2009;9:54–63.
    1. Maegawa S, Hinkal G, Kim HS, Shen L, Zhang L, et al. Widespread and tissue specific age-related DNA methylation changes in mice. Genome Res. 2010;20:332–340.
    1. Teschendorff AE, Menon U, Gentry-Maharaj A, Ramus SJ, Weisenberger DJ, et al. Age-dependent DNA methylation of genes that are suppressed in stem cells is a hallmark of cancer. Genome Res. 2010;20:440–446.
    1. Gronniger E, Weber B, Heil O, Peters N, Stab F, et al. Aging and chronic sun exposure cause distinct epigenetic changes in human skin. PLoS Genet. 2010;6:e1000971.
    1. Rakyan VK, Down TA, Maslau S, Andrew T, Yang TP, et al. Human aging-associated DNA hypermethylation occurs preferentially at bivalent chromatin domains. Genome Res. 2010;20:434–439.
    1. Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8:118–127.
    1. Horvath S, Dong J. Geometric interpretation of gene coexpression network analysis. PLoS Comput Biol. 2008;4:e1000117.
    1. Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005;4:Article17.
    1. Storey JD, Tibshirani R. Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003;100:9440–9445.
    1. Tibshirani R. The lasso method for variable selection in the Cox model. Stat Med. 1997;16:385–395.
    1. Wu TT, Chen YF, Hastie T, Sobel E, Lange K. Genome-wide association analysis by lasso penalized logistic regression. Bioinformatics. 2009;25:714–721.
    1. Park JK, Ryu JK, Lee KH, Lee JK, Yoon WJ, et al. Quantitative analysis of NPTX2 hypermethylation is a promising molecular diagnostic marker for pancreatic cancer. Pancreas. 2007;35:e9–15.
    1. Moran LB, Hickey L, Michael GJ, Derkacs M, Christian LM, et al. Neuronal pentraxin II is highly upregulated in Parkinson's disease and a novel component of Lewy bodies. Acta Neuropathol. 2008;115:471–478.
    1. Yan M, Zhang Z, Brady JR, Schilbach S, Fairbrother WJ, et al. Identification of a novel death domain-containing adaptor molecule for ectodysplasin-A receptor that is mutated in crinkled mice. Curr Biol. 2002;12:409–413.
    1. Langton AK, Herrick SE, Headon DJ. An extended epidermal response heals cutaneous wounds in the absence of a hair follicle stem cell contribution. J Invest Dermatol. 2008;128:1311–1318.
    1. Wang SC, Oelze B, Schumacher A. Age-specific epigenetic drift in late-onset Alzheimer's disease. PLoS One. 2008;3:e2698.
    1. Marks PA, Breslow R. Dimethyl sulfoxide to vorinostat: development of this histone deacetylase inhibitor as an anticancer drug. Nat Biotechnol. 2007;25:84–90.
    1. Thiede C, Prange-Krex G, Freiberg-Richter J, Bornhauser M, Ehninger G. Buccal swabs but not mouthwash samples can be used to obtain pretransplant DNA fingerprints from recipients of allogeneic bone marrow transplants. Bone Marrow Transplant. 2000;25:575–577.
    1. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
    1. R Development Core Team. 2005. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL .

Source: PubMed

3
Iratkozz fel