Research Paper Volume 10, Issue 10 pp 2832—2854

A multi-tissue full lifespan epigenetic clock for mice

Michael J. Thompson1, *,, Karolina Chwiałkowska2, , Liudmilla Rubbi1, , Aldons J. Lusis3, , Richard C. Davis3, , Anuj Srivastava4, , Ron Korstanje5, , Gary A. Churchill5, , Steve Horvath6, *,, Matteo Pellegrini1, ,

  • 1 Molecular, Cell and Developmental Biology, University of California Los Angeles, Los Angeles, CA 90095, USA
  • 2 Centre for Bioinformatics and Data Analysis, Medical University of Bialystok, Bialystok, Poland
  • 3 Department of Microbiology, Immunology and Molecular Genetics, Department of Medicine, and Department of Human Genetics, University of California Los Angeles, Los Angeles, CA 90095, USA
  • 4 The Jackson Laboratory, Farmington, CT 06032, USA
  • 5 The Jackson Laboratory, Bar Harbor, Maine 04609, USA
  • 6 Department of Human Genetics and Biostatistics, University of California Los Angeles, Los Angeles, CA 90095, USA
* Equal contribution

Received: August 16, 2018       Accepted: October 5, 2018       Published: October 21, 2018      

https://doi.org/10.18632/aging.101590
How to Cite

Copyright: © 2018 Thompson et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Human DNA-methylation data have been used to develop highly accurate biomarkers of aging ("epigenetic clocks"). Recent studies demonstrate that similar epigenetic clocks for mice (Mus Musculus) can be slowed by gold standard anti-aging interventions such as calorie restriction and growth hormone receptor knock-outs. Using DNA methylation data from previous publications with data collected in house for a total 1189 samples spanning 193,651 CpG sites, we developed 4 novel epigenetic clocks by choosing different regression models (elastic net- versus ridge regression) and by considering different sets of CpGs (all CpGs vs highly conserved CpGs). We demonstrate that accurate age estimators can be built on the basis of highly conserved CpGs. However, the most accurate clock results from applying elastic net regression to all CpGs. While the anti-aging effect of calorie restriction could be detected with all types of epigenetic clocks, only ridge regression based clocks replicated the finding of slow epigenetic aging effects in dwarf mice. Overall, this study demonstrates that there are trade-offs when it comes to epigenetic clocks in mice. Highly accurate clocks might not be optimal for detecting the beneficial effects of anti-aging interventions.

Introduction

Our understanding of age-related epigenetic changes in DNA-methylation in humans has progressed rapidly with the technical advancement of genomic platforms [114]. For mammalian genomes, DNA methylation is a modification that regulates gene expression via its presence or absence at gene promoters and enhancers. During development, germline DNA methylation is erased, but re-established in tissue-specific patterns as tissue development programs unfold after implantation [15]. Age-based methylation changes accompany the functional decline of adult stem cells [1618], and even small changes can lead to loss of regulatory control of gene transcription, either directly or via additive effects [19].

The correlation between chronological age and DNA methylation over the course of an entire lifespan is strong [2023]. Recent studies have taken advantage of this relationship to accurately estimate chronological age based on the methylation levels of multiple CpG dinucleotides [10,13,24]. For example, the human multi-tissue epigenetic age estimation method combines the weighted average of DNA methylation levels of 353 CpGs into an age estimate that is referred to as DNAm age or epigenetic age [13]. Most importantly, we and others have shown that human epigenetic age relates to biological age, not just chronological age. This is demonstrated by the finding that the discrepancy between DNAm age and chronological age (what we term “epigenetic age acceleration”) is predictive of all-cause mortality even after adjusting for a variety of known risk factors [2529]. Epigenetic age acceleration is associated with lung cancer risk [30], cognitive and physical functioning [31], Alzheimer's disease [32], centenarian status [29,33], Down syndrome [34], Werner Syndrome [35], HIV infection [36], Huntington's disease [37], obesity [38], menopause [39], osteoarthritis [40], and Parkinson's disease [41]. Moreover, we have demonstrated that the human epigenetic clock applies without change to chimpanzees [13] but it loses utility for other animals as a results of evolutionary genome sequence divergence. Moving beyond primates into the broader mammalian arena, we recently constructed an epigenetic clock for canids using DNA-methylation data from Canis familiaris (domesticated dog) and Canis lupus (wolf) [42].

Recently, other groups constructed epigenetic clocks for mice and used these to evaluate gold standard longevity interventions [43]. Petkovich, et al., derived a clock from blood samples of approximately 250 mice in order to examine changes induced by diet treatments and changes associated with genetic backgrounds that produce dwarfism (and long-lived) phenotypes [44]. Similarly, Cole, et al. examined the effects of genetic background (dwarf genotypes) and diet interventions on longevity, and their data was utilized by Wang, et al. to construct a DNA-methylation clock [45,46]. Stubbs, et al, developed a clock for multiple tissue types. Application of their clock to samples from experimental interventions yielded biologically meaningful differences in epigenetic age [47]. Overall, these independent publications led to the important insight that epigenetic clocks for mice detect anti-epigenetic aging effects of gold standard interventions such as calorie restriction and growth hormone receptor knockouts.

Our current study addressed the following aims. First, to develop a multi-tissue DNA-methylation based estimator of chronological age across the entire lifespan based on new and existing reduced representation bisulfite sequencing (RRBS) data. Second, to evaluate the robustness of reported findings surrounding gold standard anti-aging interventions using the novel epigenetic clocks. Third, to assess whether one can develop an epigenetic clock based on roughly 1k CpGs in evolutionarily conserved genomic regions.

To address these aims, we combined hundreds of new DNA-methylation samples collected from several mouse tissues with publicly available data from previous studies of mouse DNA-methylation. These data include samples obtained with RRBS and whole-genome bisulfite sequencing (WGBS). We compared clocks built with different regression methods using hundreds of thousands of CpGs as input as well as a clock constructed from a limited set of mammalian-conserved CpGs. We evaluated the performance of these clocks across samples and tissues. We applied the most accurate clock to samples from previous longevity studies of mice to measure the effects of these interventions on epigenetic aging. And, finally, we performed a GWAS analysis using epigenetic age as a trait in a subset of age-matched mouse samples covering 88 strains.

Results

Data set

Based on calculations and criteria described in the Methods section, we constructed a matrix of high confidence methylation levels for 1189 mouse samples at 193,651 CpG sites. Of these 1189 mice, 893 were used as the training set for regression models described below. This was the largest matrix we could construct while minimizing missing values to 2% of total. The remaining 296 samples were held out entirely from the training so they could be used to investigate the effects of the experimental treatments (e.g. calorie restriction) and growth hormone receptor knockouts.

Four different epigenetic clocks

We considered 4 types of epigenetic clocks. The first two clocks are constructed on the basis of all 193,651 CpG sites (covariates). In particular, the "elastic net clock" used an elastic net model to regress chronological age (dependent variable) on all methylation levels. The second clock ("ridge regression clock") used a ridge regression model instead of an elastic net regression model. The two "conserved" clocks were constructed using elastic net regression and ridge regression, respectively, using 952 highly conserved CpGs, i.e. located in highly conserved stretches of DNA (Methods).

Accuracy with respect to chronological age

We compared the four different epigenetic clocks with respect to estimating chronological age at the time of DNA sample collection (Table 1). The training set estimates of accuracy are overly optimistic and should be ignored. To arrive at unbiased estimates of the age correlation R (defined as Pearson correlation between DNAm age and chronological age) and the median absolute error (mae), the table reports three types of cross-validation estimates: i) leave-one-batch-out estimate (row "batch" in Table 1), ii) leave-one-sample-out estimate (row "sample" in Table 1), and iii) a 10 fold cross validation estimate. The three different cross validation estimates lead to the same conclusion: elastic net regression outperforms the other clocks when it comes to CV estimates of age correlations and median error. For example, the elastic net clock leads to a (leave-one-batch-out) age correlation of R=0.82 and a median error of 2.5 months. Although the conserved clocks are clearly inferior to those based on all CpGs, their accuracy remains impressive. For example, the elastic net conserved clock leads to a (leave-one-batch-out) age correlation of R=0.68 and a median error of 3.8 months.

Table 1. Summary performance statistics of epigenetic aging models (“clocks”).

CpGsEstimateRegressionAge cor.maemean model sizemodel size std dev
All CpGs ClockTraining setRidge1.000.11936510
Elastic net0.990.75820
LO-Batch-OutRidge0.793.11936410
Elastic0.822.552981
LO-Sample-OutRidge0.852.11936510
Elastic0.891.844481
10-fold CVRidge0.880.31936410
Elastic0.891.2463134
Conserved CpGs ClockTrainingRidge0.852.79520
Elastic0.911.92740
LO-Batch-Out
Ridge0.644.09520
Elastic0.683.821439
LO-Sample-OutRidge0.753.39520
Elastic0.782.42366
10-fold CVRidge0.773.59520
Elastic0.802.524723
Accuracy of estimating chronological age for 4 different epigenetic clocks. The 4 clocks differ in terms of the CpGs that were used in their construction (first column) and in terms of the underlying regression model (third column). The second column describes the method for estimating the predictive accuracy. The training set estimates are overly optimistic and should be ignored. Leave-one-batch out estimates and leave-one-sample-out estimates provide accuracy estimates that are far less biased than those obtained in the training set. The mean model size refers to the number of CpGs selected by the penalized regression model. Since the ridge regression is based on all CpGs, the standard deviation is zero.

We find that these epigenetic clocks are multi-tissue clocks, i.e. they lead to accurate age estimates in all considered tissues: results for the ridge regression clock based on all CpGs can be found Fig. 1. Analogous results for elastic net clocks based on all CpGs or based on only conserved CpG clock can be found in Suppl. Fig. 1 and Suppl. Fig. 2, respectively. We also find that accurate age estimates are made for samples taken from time points from post-natal mice to mice of advanced age, as can be seen in these figures.

Accuracy of ridge regression epigenetic age predictions. DNA methylation age (y-axis) versus chronological age (x-axis) for all mouse samples. (a) Performance of ridge regression clock based on all 192K CpGs in all training samples. The training set estimates of the accuracy are overly optimistic and should be ignored. (b) Results by tissue type of cross-validated predictions obtained by iteratively withholding one “batch” (tissue x publication). For the batch cross-validation of this clock, the global Pearson correlation between predicted and chronological age was 0.79 (p c) Scatter plots by tissue type based on DNAm age estimates made with an iterative leave-one-sample-out cross-validation. The correlation between predicted and chronologic age was 0.85 (p

Figure 1. Accuracy of ridge regression epigenetic age predictions. DNA methylation age (y-axis) versus chronological age (x-axis) for all mouse samples. (a) Performance of ridge regression clock based on all 192K CpGs in all training samples. The training set estimates of the accuracy are overly optimistic and should be ignored. (b) Results by tissue type of cross-validated predictions obtained by iteratively withholding one “batch” (tissue x publication). For the batch cross-validation of this clock, the global Pearson correlation between predicted and chronological age was 0.79 (p < 2E-195) with a mae of 3.1 months. All models in these iterative cross-validations had the same size of 193,651 CpGs. (c) Scatter plots by tissue type based on DNAm age estimates made with an iterative leave-one-sample-out cross-validation. The correlation between predicted and chronologic age was 0.85 (p < 6E-258) with a mae of 2.1 months.

Statistically speaking, the construction of epigenetic clocks is highly degenerate. That is, there are many clocks that select different sites, use slightly different weights, and achieve similar performance. For this reason, we have not emphasized the specific sites used in our clocks, as we do not believe that they are unique, but rather one set among many that could be used to construct clocks. However, we have included two Supplementary Tables of information for two of the clocks. Suppl. Table 1 provides the CpGs used in the elastic net model derived from all available methylation data, along with the regression model coefficients, and the distance of each CpG to the transcription start site of the nearest gene(s). Suppl. Table 2 provides the same information but for the elastic net model derived from the subset of RRBS data corresponding to CpGs that are evolutionarily well-conserved in mammals. For the ridge regression clock based on all CpGs, this information is provided as a text file on the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/) under super-series accession number: GSE120137.

Table 2. Datasets.

ReferenceTissueStrainAge Dist.
Novel. Current studyAdipose (56)
Blood (72)
Kidney (56)
Liver (60)
Lung (60)
Muscle (60)
C57BL/6J (200)
BALB/cByJ (164)
N=364
mean = 10.9
median = 10
std dev = 7.4
min = 1.7
max=21.3
Stubbs (2017)Cortex (16)
Heart (15)
Liver (15)
Lung (15)
C57BL/6N=61
mean = 5.1
median = 6.21
std dev = 3.4
min = 0.2
max = 9.4
Cole (2017) [45]Liver (32)Ames Prop1 Dwarf (16)
UM-HET3 (16)
N = 32
mean = 13.5
median = 22
std dev = 9.8
min = 2
max = 22
Petkovich (2017) [44]Blood (231)C57BL/6 (161)
B6D2F1 (22)
GHRKO (26)
Snell ([DW/J x C3H/HEJ]/F2) (22)
N = 231
mean = 14.7
median = 9.5
min = 0.6
max = 32.2
Novel. Current study. JAX labKidney (190)Diversity Outbred (190)N = 190
mean = 12.1
median = 12
std dev = 4.9
min = 6
max = 18
Reizel (2015) [76]Cerebellum (8)
Liver (49)
Muscle (25)
Spleen (10)
C57BL/6N = 92
mean = 2.8
median = 4.6
std dev =1.9
min = 0.23
max = 4.6
Cannon (2016) [67]Heart (5)
Liver (22)
Muscle (5)
C57BL/6N=32
mean = 0.8
median = 0.6
std dev = 0.9
min = -0.6
max = 2.1
Cannon (2014) [52]Liver (40)C57BL/6N = 40
mean = 2.07
std dev = 0
min = 2.07
max = 2.07
Orozco (2014) [66]Liver (105)91 different strainsN = 105
mean = 4
median = 4
std dev = 0
min = 4
max = 4

Diet effects on epigenetic aging

We analyzed data from 3 calorie restriction experiments and 1 rapamycin diet treatment experiment. These "test data" had been left out of the training set used in the construction of our epigenetic clocks. The most significant results could be observed for the ridge regression clock based on all CpGs (Fig. 2): significantly delayed epigenetic aging effects can be observed in calorie restricted C57BL/6 mice (p=4.2E-6, Fig. 2a) and in B6D2F1 mice (p=0.041, Fig. 2b). A similar pattern could be observed for calorie restricted HET3 mice (Fig. 2c) but the results did not quite reach statistical significance (p=0.083), which might reflect the low sample size (4 CR vs 4 chow fed HET3 mice) or the fact that the latter data had been generated using a different platform (WGBS). However, the age estimates of the WGBS samples (HET3 strain) were consistent with those obtained for RRBS samples despite the absence of WGBS samples from the training set.

Age acceleration due to diet treatments. Results obtained from ridge regression clock. A meta-analysis p-value for the 3 calorie-restriction (CR) experiments is included. (a) Calorie restriction versus standard diet in the C57BL/J strain. (b) Calorie restriction versus standard chow diet in the B6D2F1 strain. (c) Calorie restriction versus standard diet for the HET3 strain. d) Rapamycin enriched diet versus standard diet for the HET3 strain.

Figure 2. Age acceleration due to diet treatments. Results obtained from ridge regression clock. A meta-analysis p-value for the 3 calorie-restriction (CR) experiments is included. (a) Calorie restriction versus standard diet in the C57BL/J strain. (b) Calorie restriction versus standard chow diet in the B6D2F1 strain. (c) Calorie restriction versus standard diet for the HET3 strain. d) Rapamycin enriched diet versus standard diet for the HET3 strain.

CR induced anti-epigenetic aging effects could also be observed with the 2 elastic net clocks (based on all CpGs and on highly conserved CpGs, respectively) but the results were less significant than the above mentioned effects observed for the ridge regression clock (Suppl. Fig. 3).

No significant effect for rapamycin

The single comparison of mice fed with a rapamycin-enriched diet to those fed a standard diet did not yield a significant difference in age acceleration irrespective of the underlying clock (e.g. p=0.15, Fig. 2d) which might reflect the low sample size (4 Chow vs 4 Rapamycin fed mice) or the fact that the latter data had been generated using a different platform (WGBS).

Delayed epigenetic aging in dwarf mice

A few transgenic strains of mice have maximum life spans substantially greater than that of most other strains. In particular, the Ames and Snell mice, which have mutations in pituitary transcription factors (and hence are deficient in growth hormones, luteinizing hormone, thyroid-stimulating hormone, and IGF1) have extensions in maximal lifespan of up to 65% [4850].

Using publicly available data, we aimed to replicate the findings from previous publications on delayed epigenetic aging effects in dwarf mice. Three different experiments within our composite dataset were designed to examine DNA methylation and dwarfism: i) growth hormone receptor knock out mice (GHR-KO) versus wild type mice ([C57BL/6J x BALB/cByJ]/F2), ii) Snell dwarf (SD) mice versus wildtype mice WT ([DW/J x C3H/HEJ]/F2), and iii) Ames dwarf mice versus WT where these were generated by mating either homozygous (df/df) or heterozygous (df/+) dwarf males with heterozygous females (df/+), respectively. The Ames dwarf mouse line carries a recessive mutation in the Prop1 gene and homozygous animals [Prop1(df)/Prop1(df)] show dwarfness and exhibit extended lifespan [51]. Heterozygous littermates [Prop1+/Prop1(df)] were generated by breeding heterozygous females with homozygous males are of normal size.

The ridge regression clock based on all CpGs managed to detect a delayed-epigenetic aging effect in all three types of dwarf mice (Fig. 3). Despite low sample sizes, the trend for homozygous dwarf strains to slow the epigenetic clock is clear and statistically significant with a meta-analysis p-value of 8E-7. However, these results are not statistically robust with respect to different epigenetic clocks: the association of dwarfism and slow epigenetic aging could not be detected with the same significance with the two elastic net regression clocks (Suppl. Fig. 4).

Age acceleration and Dwarfism in mice. Results obtained from ridge regression clock. A meta-analysis p-value for the 3 experiments is included. (a) Genetic knockout dwarf mice versus wild type. (b) Snell dwarf mice versus wild type. c) Ames Dwarf mice versus wild type.

Figure 3. Age acceleration and Dwarfism in mice. Results obtained from ridge regression clock. A meta-analysis p-value for the 3 experiments is included. (a) Genetic knockout dwarf mice versus wild type. (b) Snell dwarf mice versus wild type. c) Ames Dwarf mice versus wild type.

Maternal diet effects on epigenetic aging

Cannon, et al. investigated the potential influence of maternal diet on gene expression and DNA methylation in their offspring [52]. Our ridge regression clock reveals that the slowest epigenetic aging effects can be observed in low fat diet-fed offspring of low fat diet-fed mothers (Fig. 4).

Age acceleration and maternal diet. Results obtained from ridge regression clock. (a) Offspring of mothers fed a high fat diet (HFD) who were fed either a high fat or low fat diet (LFD). (b) Offspring of mothers fed a low fat diet who were fed either a high fat or low fat diet.

Figure 4. Age acceleration and maternal diet. Results obtained from ridge regression clock. (a) Offspring of mothers fed a high fat diet (HFD) who were fed either a high fat or low fat diet (LFD). (b) Offspring of mothers fed a low fat diet who were fed either a high fat or low fat diet.

Although, this finding is biologically plausible it is not statistically robust with respect to other epigenetic clocks. In particular, it cannot be observed for the two elastic net regression clocks (Suppl. Fig. 5)

GWAS of epigenetic age in mice

We used epigenetic age as trait in our Genome Wide Association Study (GWAS) in the Hybrid Mouse Diversity Panel (HMDP). We calculated epigenetic age via a cross validation approach in order to avoid overfitting. Specifically, epigenetic age was computed for 88 strains using the ridge regression based clock and leaving out from the training set the sample whose age was estimated. The mean calculated epigenetic age was 4.18 months (±0.95) in the range of 2.3-7.1 months with a median of 4.1 months. A set of 196,148 SNPs (MAF>5%) was used for association studies. We used linear mixed models to correct for population structure using the pyLMM software. As our cohort size was limited, we were not able to identify peaks whose significance was beyond the Bonferroni threshold in this analysis (Fig. 5a). Therefore the results presented in Table 3 with the top 10 SNPS are only suggestive of an association and will need to be confirmed in the future with a larger cohort.

Genome-wide association results for DNAm Age. (a) Manhattan plot presenting genome-wide association results for DNAm Age. Epigenetic age predictions were calculated using all CpGs clock with ridge regression and leave-one-sample-out estimates. GWAS analysis was based on linear mixed model and a set of 196,148 SNPs (MAF > 0.05) from HMDP mice strains. (b) This SNP as identified using GWAS analysis of epigenetic age predictions. It is located in an LD block on chromosome 6 and contains the genes Npy, Mpp6, Gsdme and Osbp13. A one-sided t-test of DNAm ages between the two allelic groups shown is statistically significant. (c) It is located in an LD block on chromosome 6 and contains the genes Npy, Mpp6, Gsdme and Osbp13. A one-sided t-test of DNAm ages between the two allelic groups shown is statistically significant.

Figure 5. Genome-wide association results for DNAm Age. (a) Manhattan plot presenting genome-wide association results for DNAm Age. Epigenetic age predictions were calculated using all CpGs clock with ridge regression and leave-one-sample-out estimates. GWAS analysis was based on linear mixed model and a set of 196,148 SNPs (MAF > 0.05) from HMDP mice strains. (b) This SNP as identified using GWAS analysis of epigenetic age predictions. It is located in an LD block on chromosome 6 and contains the genes Npy, Mpp6, Gsdme and Osbp13. A one-sided t-test of DNAm ages between the two allelic groups shown is statistically significant. (c) It is located in an LD block on chromosome 6 and contains the genes Npy, Mpp6, Gsdme and Osbp13. A one-sided t-test of DNAm ages between the two allelic groups shown is statistically significant.

Table 3. Top ten SNPs from GWAS analysis of DNAm age predictions corresponding to peaks connected to LD blocks in HDMP (~2 Mbp). P-values were computed with a linear mixed-model (LMM).

SNP IDChrPositionLMM
p-value
Nearest gene(s)
JAX001898826771044791,08E-04Ctnna2, Lrrtm1
JAX001411866551243511,59E-04Plekha8, Mturn, Znrf2, Nod1, Ggct, Gars, Crhr2, Inmt, Mindy4, Aqp1, Ghrhr, Adcyap1r1
JAX006138026736412781,67E-04Dnah6, Suclg1, 4931417E11Rik
JAX0065189871150700693,06E-04Calca, Calcb, Insc, Sox6
JAX0037352214146570814,48E-04Olfr720, Olfr31, Il3ra, Slc4a7, Nek10
JAX000499271494981344,83E-04Fhit
JAX001407996499763985,06E-04Npy, Mpp6, Gsdme, Osbp13
JAX001894884959907915,48E-04Fggy, Hook1, Cyp2j13, Cyp2j12, Cyp2j11, Cyp2j8
JAX0008797019257170225,66E-04Kank1, Dmrt1, Dmrt3, Dmrt2
JAX0037402014175878717,02E-04Thrb

We analysed the genes that were found proximal (500 kbp up- and downstream) to the identified peak SNP site. This list contains several genes that have previously been implicated in aging: Aqp1 (Aquaporin 1) [53], Npy (Neuropeptide Y) [54], Adcyap1r1 (Adenylate cyclase-activating polypeptide type I receptor 1) as a receptor for PACAP (Pituitary Adenylate Cyclase-Activating Peptide) [55,56]. The most notable is Npy, which encodes a hormone responsible for appetite control [57], regulation of fat metabolism [58], and plays a critical role in caloric restriction (CR) mediated lifespan extension [59]. As examples, we found statistically significant differences (p < 0.001) in average DNAm ages of mice strains having allele C or T within SNP JAX00140799 (Fig. 5b) and between mice strains having allele T or C within SNP JAX00373522 (Fig. 5c).

Discussion

Based on multiple tissue samples taken from previous studies and our own in-house collection we compiled a dataset of 1189 mouse DNA methylation measurements across hundreds of thousands of CpGs. These samples represent the most comprehensive dataset thus far of matched single base resolution methylomes in mice across multiple tissues and ages.

We demonstrate that these data enable construction of highly accurate multi-tissue age estimation methods (epigenetic clocks) for mice that apply to the entire life course (from birth to old age). We demonstrate that these clocks perform well on new tissues not included in the training of the clock by performing tissue exclusion cross-validation. This gives us confidence that these clocks will work on new samples from other tissue types as well. However, we cannot rule out that these clocks fail in specific cell types. Epigenetic age estimators that focus on specific tissues or cell types can have greater accuracy than pan tissue age estimators [60].

Our study leads to several novel insights. First, our first prototype of an age estimator based on fewer than 1000 highly conserved CpGs demonstrates that it will be feasible to build highly accurate DNAm age estimator on the basis of highly conserved CpGs.

Second, we find that epigenetic clocks that are optimal for estimating age (namely those based on elastic net regression) may be inferior to less accurate clocks (based on ridge regression) when it comes to gold standard anti-aging interventions. Only our ridge regression clock manages to corroborate most of the previously reported findings, e.g. only the ridge clock showed that dwarf strains show slower epigenetic aging relative to wild-type strains. The anti-epigenetic aging effects of calorie restriction are highly robust and could be observed with all clocks. Moreover, by utilizing epigenetic ages as phenotypic traits in a GWAS study of 88 strains of mice we found suggestive associations with several genes, including neuropeptide Y whose role in appetite control and calorie restriction mediated lifespan extension is well documented. However, none of our clocks managed to detect an anti-aging effect of rapamycin in a small data set which might reflect the low sample sizes or technical reasons including low coverage afforded by the measurement platform (WGBS).

All clocks were able to detect a slowing of the epigenetic clock in mice fed a calorie restricted diet, though with differing sensitivity, suggesting that the effects of calorie restriction are pervasive across the methylomes. In contrast, the slowing of epigenetic aging in mice fed a low fat diet for two generations were not detected by the clocks with fewer CpGs (elastic net), suggesting these effects are either more subtle or more localized in the methylome.

These results suggest that the multi-tissue ridge regression DNA-methylation clock is most useful in assessing “biological age” for a variety of treatments, experimental interventions, and genetic backgrounds. However, the elastic net clocks are better for assessing chronological age. We evaluated both ridge and lasso regression in previous studies with human data and found that lasso outperformed ridge not just in terms of accuracy but also in terms of interpretability (unreported findings). Therefore, it is a curious finding that ridge regression has some merits when it comes to mice.

We speculate that one reason that ridge regression works best in our context is that our dataset is more heterogeneous than those used in previous studies. Our dataset includes mice of different ages, strains and diverse tissues, all collected in different labs and resulting in a whole that is larger than any previous dataset. Because genetic diversity in mice is high, it is possible that lasso models that only use a limited number of sites are more prone to be influence by genetic variation (as DNA methylation is often associated with genetic variation). Thus, on the whole, it is possible that the ridge approach minimizes these effects by using all sites, and thus leads to the most robust overall performance.

The DNAm age estimates from our mouse clocks exhibit a correlation coefficient with chronological age that ranges from R=0.79 to 0.89 (Table 1). These correlations are only slightly weaker than those observed for human studies (R=0.96 for the pan tissue estimator from Horvath). We have no doubt that more accurate mouse clocks can be built by reducing technical variation and by employing even larger data sets.

We acknowledge several limitations. Our genetic study of epigenetic aging rates was under-powered. Large scale human studies have implicated genome-wide significant loci including the TERT gene [6163]. We did not assess the intra-assay variation of replicate samples in the current article but refer the interested reader to relevant articles [59,60].

Rigorous quantitative comparisons to previously published mouse clocks (e.g. comparing age correlation values) could not be made due to technical differences in the processing of sequence data, in the estimation and correction of methylation calls, in the limited tissue sampling of previous clocks, and, most importantly, in the simple absence of particular CpGs in the various datasets and clock models. The subset of CpGs with high coverage in one RRBS data set tends to exhibit poor overlap with a subset of CpGs from another RRBS data set. This poor overlap of CpGs makes it difficult to validate epigenetic clocks based on RRBS data. We are currently working on a custom methylation array platform that avoids these pitfalls.

Methods

Data sets

We generated reduced representation bisulfite sequencing (RRBS) methylation data for mouse adipose, blood, liver, and kidney, muscle, and lung tissue samples using the protocol below.

DNA methylation assay

Genomic DNA was isolated by standard phenol-chloroform extraction method and used as input to prepare Reduced Representation Bisulfite Sequencing (RRBS) libraries as described previously [64] with minor modifications. For each sample 50-100 ng of purified genomic DNA was digested with 20 U of MspI (NEB, cat # R0106L) at 37°C o/n in the presence of RNase Cocktail Mix (Ambion, cat # AM2286). End-repair and dA-tailing was performed by the addition of Klenow Fragment 3’->5’ exo- (NEB, cat # M0212L) in the presence of dATP, dGTP and d5mCTP (Fermentas). Adapter Ligation was performed by the addition of 0.3 µl of Illumina TruSeq methylated Adapters (Illumina, TruSeq Nano cat# FC-121-4001) and 2 µl of Illumina Ligation Mix 2 (Illumina, TruSeq Nano cat# FC-121-4001). Samples were pooled and purified using an equal volume of SPRI beads (Beckman Coulter, cat # B23318). Size-selection was performed using SPRI beads to enrich for fragments from 200 to 300 bp. Bisulfite treatment was performed using Epitect Bisulfite kit (QIAGEN, cat # 59104) according to manufacturer's protocol, except that two consecutive rounds of conversion are performed, for a total of 10 hr of incubation. Purified converted DNA was PCR amplified using MyTaq HS Mix (Bioline, cat# BIO-25045) and TruSeq PCR Primer Cocktail (Illumina, TruSeq Nano cat# FC-121-4001) according to the following protocol: initial denaturation at 98°C for 30s; 12 cycles of 98°C for 15s, 60°C for 30s, 72°C for 30s; final extension at 72°C for 5 min. Amplified libraries were purified twice with an equal amount of SPRI beads to remove primer and adapter dimers. Libraries were sequenced 100 bp single-end on an Illumina HiSeq4000. For the kidney data from JAX laboratories, the sequencing protocol was as follows. RRBS libraries were prepared using 100 ng DNA, the Ovation RRBS Methyl-Seq System 1–16 (NuGEN Technologies, San Carlo, CA) part number 0353, and the EpiTect Fast 96 Bisulfite Conversion kit (Qiagen, Hilden, Germany) part number 59720. The manufacturer’s protocols were followed except for the number of PCR cycles in the library amplification step, which was increased, from 12 to 13. Libraries were quantified using the Library Quantification Kit (Kapa Biosystems, Wilmington, MA) part number KK4835, normalized to 10nM, and pooled in groups of 12. Each pool was sequenced 1 x 100 bp on one lane of the HiSeq2500 (Illumina, San Diego, CA) at The New York Genome Center (New York, NY).

These datasets were integrated with RRBS data made available to the public via the GEO repository [65]. We included datasets from previous RRBS-based “epigenetic clock” studies [44,47] along with RRBS data from an EWAS study of metabolic traits [66], from a study of maternal diet effects on gene expression and DNA methylation [52], from a study of post-natal hepatocyte development [67], from a multi-tissue study of sex hormone effects on DNA methylation [68].

Kidney data from the Jackson Laboratory

Kidneys were collected from male and female Diversity Outbred mice at ages 6, 12 and 18 months. Mice were group housed in SPF condition and fed a standard lab chow diet (5K0G) with 6% calories from fat. Tissues were flash frozen in LN2, pulverized and mixed prior to DNA extraction. RRBS sequencing was carried out at the New York Genome Center. Although whole-genome bisulfite-sequencing data (WGBS) has a variety of characteristic differences from RRBS, we obtained the set of this data collected previously to examine longevity interventions in mice and build an epigenetic clock [45,46].

Data processing

Where possible, we downloaded the raw sequencing files from previous studies via (GEO) and performed alignments and methylation calling identically as for our in-house data using BS_Seeker2 with default parameters [69]. All mouse methylation data in this study utilized mouse genome mm10 coordinates. When technical replicates were available for a given sample, they were merged by summing the sequencing counts.

For the kidney data from JAX laboratories, a Bismark-based pipeline was initially used as follows [70]. All the samples were subjected to QC using the trim_galore module then trimming of the diversity adapters was performed by trimRRBSdiversityAdaptCustomers.py script from NuGen. High quality trimmed reads were aligned to all eight diversity outbred founder strains (A/J, C57BL/6J, 129S1/SvImJ, NOD/ShiLtJ, NZO/HlLtJ, CAST/EiJ, PWK/PhJ, and WSB/EiJ,) separately using Bismark at default parameters. The alignment to founder genomes except C57BL/6J were converted to reference genome coordinates (mm10) using g2gtools v1. Then, we selected the reads which were mapped to the same locus in multiple founders strains (assigned to founder strains with minimum edit distance) or mapped uniquely to one founder strains by custom in-house script. The bed file of estimated methylation proportion of each founder (except C57BL/6J ) was converted to reference genome coordinate by g2gtools v1 and, finally, combined to provide the methylation proportion in each diversity outbred animals.

For each CpG site in each sample we estimated the methylation frequency as the number of methylated mapped read counts over the total mapped read counts. Where available, the counts for the forward and reverse Cytosines of the CpG were pooled and treated as a single measurement. We then computed a 95% confidence interval with a Bayesian approach using a Beta distribution (0.5,0.5) (“Jeffrey’s Prior”) for all methylation values [71, 72]. For inclusion in our analysis, we required that each CpG site had confident methylation frequencies in at least 95% of samples. Confidence was defined as having a confidence interval smaller than 0.50. True missing values or measurements failing that confidence interval filter were imputed using k-nearest-neighbor approach with k=5. This select strategy resulted in CpGs whose mean methylation levels ranged from zero to 1 (Suppl. Fig. 6).

Sample exclusion

In order to maximize the number of samples and coverage of the methylomes, it was necessary to exclude a number of samples both from our new data and from previously published datasets. First, we removed samples with fewer than 500,000 measured CpGs. Next, after an initial matrix was constructed, we iteratively removed samples with the most missing values until we arrived at a matrix with ~2% total missing values.

Penalized regression models

Penalized regression models were created with glmnet [73]. We investigated models produced by both elastic net regression (alpha=0.5) and ridge regression (alpha=0). The optimal penalty parameters in all cases were determined automatically by using a 10 fold internal cross-validation (cv.glmnet) on the training set. By definition, the alpha value for the elastic net regression was set to 0.5 (midpoint between ridge and lasso type regression) and was not optimized for model performance. We omitted the results from lasso regression models (alpha=1) because the age estimates tended to be less accurate than those from elastic net regression.

The covariates in our data (methylation of CpGs) are known to have a high degree of multicollinearity. While lasso and elastic net regression allow regression coefficients to go to zero and thus yield a sort of “feature selection” which is desirable for interpretability, the correlations among methylation sites may contain subtle information that might be useful to retain (which supports the use of ridge regression).

Cross-validation estimates of accuracy

We performed three types of cross-validation schemes for arriving at unbiased (or at least less biased) estimates of the accuracy of the different DNAm based age estimators. One type consisted of leaving out a single sample (LOOCV) from the regression, predicting an age for that sample, and iterating over all samples. The second type (10-fold) was similar to the first except that 10% of samples were withheld per iteration. The third type consisted of iteratively leaving out all samples of a particular “batch” where batched was defined as combination of tissue type and publication of origin. For example, three batches resulted from a single publication if the underlying RRBS data were obtained from 3 distinct tissues. Samples from longevity intervention studies (CR, Rapamycin, dwarf mice) along with their respective control samples were excluded from all instances of training-set construction.

Conserved-CpG clock

Our team is currently developing a mammalian DNA-methylation array for measuring methylation levels in mammals. The primary selection criteria for CpGs for the probe-design of this bead-chip array was conservation of the local sequence context of the CpG across 60 sequenced mammalian genomes. We obtained the list of candidate CpGs in this effort and intersected it with CpGs in our high-quality RRBS data.

Epigenetic age acceleration

To investigate effects of biological interventions and genetic background on epigenetic aging, we employ a quantity termed “age acceleration”. In the simplest form, it is just the difference between the epigenetic age estimated by the clock and the chronological age. However, this measure can be age-dependent itself, causing difficulty in interpretation. Instead, age acceleration is computed as the residual, per sample, after fitting predicted ages to chronological ages. This fitting is done on a per-batch basis. P-values for age acceleration comparisons found in Figs 2, 3, 4, and Suppl. Figs 3, 4 and 5 were obtained using the non-parametric Kruskal-Wallis test.

Genome-wide association study

GWAS was performed using 88 strains from the Hybrid Mouse Diversity Panel (HMDP, listed in Suppl. Table 1), that have been extensively used as a resource for systems genetics analyses [66,74,75]. Of the total of 459 895 SNPs, we selected a set of 196 148 SNPs that exhibited minor allele frequency greater than 5%. DNAm ages, computed using all CpGs clock with ridge regression and leave-one-sample-out estimates, were treated as phenotypes in the association studies. All of the mice were at chronological age of 4 months. GWAS was conducted using the linear mixed model python package pyLMM to account for population structure and relatedness among the mouse strains. We selected the top SNPs in each of the peaks from the pyLMM analysis that were located a minimum of 2 Mbp apart, which is the average of LD block size for SNPs in the HMDP.

Data availability

Raw sequencing data and processed data for samples collected at UCLA and JAX, as well as re-processed data from previous studies have been made available at the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/) under super-series accession number: GSE120137

Conflicts of Interest

The authors declare there are no potential conflicts of interest.

Funding

SH, MP, and MT were supported by the Paul G. Allen Frontiers Group and by 1R21AG049400–01A1. MJT acknowledges support from a QCB Collaboratory Postdoctoral Fellowship, and the QCB Collaboratory community directed by Matteo Pellegrini. KC was a beneficiary of the ETIUDA scholarship funded by the Polish National Science Centre (NCN) under grant no. 2016/20/T/NZ2/00577. GAC and RK were supported by grant P30 AG038070-09 NIA.

References

  • 1. Christensen BC, Houseman EA, Marsit CJ, Zheng S, Wrensch MR, Wiemels JL, Nelson HH, Karagas MR, Padbury JF, Bueno R, Sugarbaker DJ, Yeh RF, Wiencke JK, Kelsey KT. Aging and environmental exposures alter tissue-specific DNA methylation dependent upon CpG island context. PLoS Genet. 2009; 5:e1000602. https://doi.org/10.1371/journal.pgen.1000602 [PubMed]
  • 2. Bollati V, Schwartz J, Wright R, Litonjua A, Tarantini L, Suh H, Sparrow D, Vokonas P, Baccarelli A. Decline in genomic DNA methylation through aging in a cohort of elderly subjects. Mech Ageing Dev. 2009; 130:234–39. https://doi.org/10.1016/j.mad.2008.12.003 [PubMed]
  • 3. Rakyan VK, Down TA, Maslau S, Andrew T, Yang TP, Beyan H, Whittaker P, McCann OT, Finer S, Valdes AM, Leslie RD, Deloukas P, Spector TD. Human aging-associated DNA hypermethylation occurs preferentially at bivalent chromatin domains. Genome Res. 2010; 20:434–39. https://doi.org/10.1101/gr.103101.109 [PubMed]
  • 4. Teschendorff AE, Menon U, Gentry-Maharaj A, Ramus SJ, Weisenberger DJ, Shen H, Campan M, Noushmehr H, Bell CG, Maxwell AP, Savage DA, Mueller-Holzner E, Marth C, et al. Age-dependent DNA methylation of genes that are suppressed in stem cells is a hallmark of cancer. Genome Res. 2010; 20:440–46. https://doi.org/10.1101/gr.103606.109 [PubMed]
  • 5. Horvath S, Zhang Y, Langfelder P, Kahn RS, Boks MP, van Eijk K, van den Berg LH, Ophoff RA. Aging effects on DNA methylation modules in human brain and blood tissue. Genome Biol. 2012; 13:R97. https://doi.org/10.1186/gb-2012-13-10-r97 [PubMed]
  • 6. Numata S, Ye T, Hyde TM, Guitart-Navarro X, Tao R, Wininger M, Colantuoni C, Weinberger DR, Kleinman JE, Lipska BK. DNA methylation signatures in development and aging of the human prefrontal cortex. Am J Hum Genet. 2012; 90:260–72. https://doi.org/10.1016/j.ajhg.2011.12.020 [PubMed]
  • 7. Alisch RS, Barwick BG, Chopra P, Myrick LK, Satten GA, Conneely KN, Warren ST. Age-associated DNA methylation in pediatric populations. Genome Res. 2012; 22:623–32. https://doi.org/10.1101/gr.125187.111 [PubMed]
  • 8. Johansson A, Enroth S, Gyllensten U. Continuous Aging of the Human DNA Methylome Throughout the Human Lifespan. PLoS One. 2013; 8:e67378. https://doi.org/10.1371/journal.pone.0067378 [PubMed]
  • 9. Day K, Waite LL, Thalacker-Mercer A, West A, Bamman MM, Brooks JD, Myers RM, Absher D. Differential DNA methylation with age displays both common and dynamic features across human tissues that are influenced by CpG landscape. Genome Biol. 2013; 14:R102. https://doi.org/10.1186/gb-2013-14-9-r102 [PubMed]
  • 10. Bocklandt S, Lin W, Sehl ME, Sánchez FJ, Sinsheimer JS, Horvath S, Vilain E. Epigenetic predictor of age. PLoS One. 2011; 6:e14821. https://doi.org/10.1371/journal.pone.0014821 [PubMed]
  • 11. Garagnani P, Bacalini MG, Pirazzini C, Gori D, Giuliani C, Mari D, Di Blasio AM, Gentilini D, Vitale G, Collino S, Rezzi S, Castellani G, Capri M, et al. Methylation of ELOVL2 gene as a new epigenetic marker of age. Aging Cell. 2012; 11:1132–34. https://doi.org/10.1111/acel.12005 [PubMed]
  • 12. Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, Klotzle B, Bibikova M, Fan JB, Gao Y, Deconde R, Chen M, Rajapakse I, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013; 49:359–67. https://doi.org/10.1016/j.molcel.2012.10.016 [PubMed]
  • 13. Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013; 14:R115. https://doi.org/10.1186/gb-2013-14-10-r115 [PubMed]
  • 14. Lin Q, Weidner CI, Costa IG, Marioni RE, Ferreira MR, Deary IJ, Wagner W. DNA methylation levels at individual age-associated CpG sites can be indicative for life expectancy. Aging (Albany NY). 2016; 8:394–401. https://doi.org/10.18632/aging.100908 [PubMed]
  • 15. Cedar H, Bergman Y. Programming of DNA methylation patterns. Annu Rev Biochem. 2012; 81:97–117. https://doi.org/10.1146/annurev-biochem-052610-091920 [PubMed]
  • 16. Przybilla J, Galle J, Rohlf T. Is adult stem cell aging driven by conflicting modes of chromatin remodeling? BioEssays. 2012; 34:841–48. https://doi.org/10.1002/bies.201100190 [PubMed]
  • 17. Przybilla J, Rohlf T, Loeffler M, Galle J. Understanding epigenetic changes in aging stem cells--a computational model approach. Aging Cell. 2014; 13:320–28. https://doi.org/10.1111/acel.12177 [PubMed]
  • 18. Beerman I, Bock C, Garrison BS, Smith ZD, Gu H, Meissner A, Rossi DJ. Proliferation-dependent alterations of the DNA methylation landscape underlie hematopoietic stem cell aging. Cell Stem Cell. 2013; 12:413–25. https://doi.org/10.1016/j.stem.2013.01.017 [PubMed]
  • 19. Beerman I, Rossi DJ. Epigenetic regulation of hematopoietic stem cell aging. Exp Cell Res. 2014; 329:192–99. https://doi.org/10.1016/j.yexcr.2014.09.013 [PubMed]
  • 20. Gibbs WW. Biomarkers and ageing: the clock-watcher. Nature. 2014; 508:168–70. https://doi.org/10.1038/508168a [PubMed]
  • 21. Jung M, Pfeifer GP. Aging and DNA methylation. BMC Biol. 2015; 13:7. https://doi.org/10.1186/s12915-015-0118-4 [PubMed]
  • 22. Benayoun BA, Pollina EA, Brunet A. Epigenetic regulation of ageing: linking environmental inputs to genomic stability. Nat Rev Mol Cell Biol. 2015; 16:593–610. https://doi.org/10.1038/nrm4048 [PubMed]
  • 23. Horvath S, Raj K. DNA methylation-based biomarkers and the epigenetic clock theory of ageing. Nat Rev Genet. 2018; 19:371–84. https://doi.org/10.1038/s41576-018-0004-3 [PubMed]
  • 24. Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, Klotzle B, Bibikova M, Fan JB, Gao Y, Deconde R, Chen M, Rajapakse I, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013; 49:359–67. https://doi.org/10.1016/j.molcel.2012.10.016 [PubMed]
  • 25. Marioni RE, Shah S, McRae AF, Chen BH, Colicino E, Harris SE, Gibson J, Henders AK, Redmond P, Cox SR, Pattie A, Corley J, Murphy L, et al. DNA methylation age of blood predicts all-cause mortality in later life. Genome Biol. 2015; 16:25. https://doi.org/10.1186/s13059-015-0584-6 [PubMed]
  • 26. Christiansen L, Lenart A, Tan Q, Vaupel JW, Aviv A, McGue M, Christensen K. DNA methylation age is associated with mortality in a longitudinal Danish twin study. Aging Cell. 2016; 15:149–54. https://doi.org/10.1111/acel.12421 [PubMed]
  • 27. Perna L, Zhang Y, Mons U, Holleczek B, Saum KU, Brenner H. Epigenetic age acceleration predicts cancer, cardiovascular, and all-cause mortality in a German case cohort. Clin Epigenetics. 2016; 8:64. https://doi.org/10.1186/s13148-016-0228-z [PubMed]
  • 28. Chen BH, Marioni RE, Colicino E, Peters MJ, Ward-Caviness CK, Tsai PC, Roetker NS, Just AC, Demerath EW, Guan W, Bressler J, Fornage M, Studenski S, et al. DNA methylation-based measures of biological age: meta-analysis predicting time to death. Aging (Albany NY). 2016; 8:1844–65. https://doi.org/10.18632/aging.101020 [PubMed]
  • 29. Horvath S, Pirazzini C, Bacalini MG, Gentilini D, Di Blasio AM, Delledonne M, Mari D, Arosio B, Monti D, Passarino G, De Rango F, D’Aquila P, Giuliani C, et al. Decreased epigenetic age of PBMCs from Italian semi-supercentenarians and their offspring. Aging (Albany NY). 2015; 7:1159–70. https://doi.org/10.18632/aging.100861 [PubMed]
  • 30. Levine ME, Hosgood HD, Chen B, Absher D, Assimes T, Horvath S. DNA methylation age of blood predicts future onset of lung cancer in the women’s health initiative. Aging (Albany NY). 2015; 7:690–700. https://doi.org/10.18632/aging.100809 [PubMed]
  • 31. Marioni RE, Shah S, McRae AF, Ritchie SJ, Muniz-Terrera G, Harris SE, Gibson J, Redmond P, Cox SR, Pattie A, Corley J, Taylor A, Murphy L, et al. The epigenetic clock is correlated with physical and cognitive fitness in the Lothian Birth Cohort 1936. Int J Epidemiol. 2015; 44:1388–96. https://doi.org/10.1093/ije/dyu277 [PubMed]
  • 32. Levine ME, Lu AT, Bennett DA, Horvath S. Epigenetic age of the pre-frontal cortex is associated with neuritic plaques, amyloid load, and Alzheimer’s disease related cognitive functioning. Aging (Albany NY). 2015; 7:1198–211. https://doi.org/10.18632/aging.100864 [PubMed]
  • 33. Horvath S, Mah V, Lu AT, Woo JS, Choi OW, Jasinska AJ, Riancho JA, Tung S, Coles NS, Braun J, Vinters HV, Coles LS. The cerebellum ages slowly according to the epigenetic clock. Aging (Albany NY). 2015; 7:294–306. https://doi.org/10.18632/aging.100742 [PubMed]
  • 34. Horvath S, Garagnani P, Bacalini MG, Pirazzini C, Salvioli S, Gentilini D, Di Blasio AM, Giuliani C, Tung S, Vinters HV, Franceschi C. Accelerated epigenetic aging in Down syndrome. Aging Cell. 2015; 14:491–95. https://doi.org/10.1111/acel.12325 [PubMed]
  • 35. Maierhofer A, Flunkert J, Oshima J, Martin GM, Haaf T, Horvath S. Accelerated epigenetic aging in Werner syndrome. Aging (Albany NY). 2017; 9:1143–52. https://doi.org/10.18632/aging.101217 [PubMed]
  • 36. Horvath S, Levine AJ. HIV-1 infection accelerates age according to the epigenetic clock. J Infect Dis. 2015; 212:1563–73. https://doi.org/10.1093/infdis/jiv277 [PubMed]
  • 37. Horvath S, Langfelder P, Kwak S, Aaronson J, Rosinski J, Vogt TF, Eszes M, Faull RL, Curtis MA, Waldvogel HJ, Choi OW, Tung S, Vinters HV, et al. Huntington’s disease accelerates epigenetic aging of human brain and disrupts DNA methylation levels. Aging (Albany NY). 2016; 8:1485–512. https://doi.org/10.18632/aging.101005 [PubMed]
  • 38. Horvath S, Erhart W, Brosch M, Ammerpohl O, von Schönfels W, Ahrens M, Heits N, Bell JT, Tsai PC, Spector TD, Deloukas P, Siebert R, Sipos B, et al. Obesity accelerates epigenetic aging of human liver. Proc Natl Acad Sci USA. 2014; 111:15538–43. https://doi.org/10.1073/pnas.1412759111 [PubMed]
  • 39. Levine ME, Lu AT, Chen BH, Hernandez DG, Singleton AB, Ferrucci L, Bandinelli S, Salfati E, Manson JE, Quach A, Kusters CD, Kuh D, Wong A, et al. Menopause accelerates biological aging. Proc Natl Acad Sci USA. 2016; 113:9327–32. https://doi.org/10.1073/pnas.1604558113 [PubMed]
  • 40. Vidal L, Lopez-Golan Y, Rego-Perez I, Horvath S, Blanco FJ, Riancho JA, Gomez-Reino JJ, Gonzalez A. Specific increase of methylation age in osteoarthritis cartilage. Osteoarthritis Cartilage. 2016; 24:S63. https://doi.org/10.1016/j.joca.2016.01.140
  • 41. Horvath S, Ritz BR. Increased epigenetic age and granulocyte counts in the blood of Parkinson’s disease patients. Aging (Albany NY). 2015; 7:1130–42. https://doi.org/10.18632/aging.100859 [PubMed]
  • 42. Thompson MJ, vonHoldt B, Horvath S, Pellegrini M. An epigenetic aging clock for dogs and wolves. Aging (Albany NY). 2017; 9:1055–68. https://doi.org/10.18632/aging.101211 [PubMed]
  • 43. Wagner W. Epigenetic aging clocks in mice and men. Genome Biol. 2017; 18:107. https://doi.org/10.1186/s13059-017-1245-8 [PubMed]
  • 44. Petkovich DA, Podolskiy DI, Lobanov AV, Lee SG, Miller RA, Gladyshev VN. Using DNA Methylation Profiling to Evaluate Biological Age and Longevity Interventions. Cell Metab. 2017; 25:954–960.e6. https://doi.org/10.1016/j.cmet.2017.03.016 [PubMed]
  • 45. Cole JJ, Robertson NA, Rather MI, Thomson JP, McBryan T, Sproul D, Wang T, Brock C, Clark W, Ideker T, Meehan RR, Miller RA, Brown-Borg HM, Adams PD. Diverse interventions that extend mouse lifespan suppress shared age-associated epigenetic changes at critical gene regulatory regions. Genome Biol. 2017; 18:58. https://doi.org/10.1186/s13059-017-1185-3 [PubMed]
  • 46. Wang T, Tsui B, Kreisberg JF, Robertson NA, Gross AM, Yu MK, Carter H, Brown-Borg HM, Adams PD, Ideker T. Epigenetic aging signatures in mice livers are slowed by dwarfism, calorie restriction and rapamycin treatment. Genome Biol. 2017; 18:57. https://doi.org/10.1186/s13059-017-1186-2 [PubMed]
  • 47. Stubbs TM, Bonder MJ, Stark AK, Krueger F, von Meyenn F, Stegle O, Reik W, Reik W, and BI Ageing Clock Team. Multi-tissue DNA methylation age predictor in mouse. Genome Biol. 2017; 18:68. https://doi.org/10.1186/s13059-017-1203-5 [PubMed]
  • 48. Brown-Borg HM. The somatotropic axis and longevity in mice. Am J Physiol Endocrinol Metab. 2015; 309:E503–10. https://doi.org/10.1152/ajpendo.00262.2015 [PubMed]
  • 49. Coschigano KT, Holland AN, Riders ME, List EO, Flyvbjerg A, Kopchick JJ. Deletion, but not antagonism, of the mouse growth hormone receptor results in severely decreased body weights, insulin, and insulin-like growth factor I levels and increased life span. Endocrinology. 2003; 144:3799–810. https://doi.org/10.1210/en.2003-0374 [PubMed]
  • 50. Flurkey K, Papaconstantinou J, Miller RA, Harrison DE. Lifespan extension and delayed immune and collagen aging in mutant mice with defects in growth hormone production. Proc Natl Acad Sci USA. 2001; 98:6736–41. https://doi.org/10.1073/pnas.111158898 [PubMed]
  • 51. Brown-Borg HM, Borg KE, Meliska CJ, Bartke A. Dwarf mice and the ageing process. Nature. 1996; 384:33. https://doi.org/10.1038/384033a0 [PubMed]
  • 52. Cannon MV, Buchner DA, Hester J, Miller H, Sehayek E, Nadeau JH, Serre D. Maternal nutrition induces pervasive gene expression changes but no detectable DNA methylation differences in the liver of adult offspring. PLoS One. 2014; 9:e90335. https://doi.org/10.1371/journal.pone.0090335 [PubMed]
  • 53. Dempsey KM, Ali HH. Identifying aging-related genes in mouse hippocampus using gateway nodes. BMC Syst Biol. 2014; 8:62. https://doi.org/10.1186/1752-0509-8-62 [PubMed]
  • 54. Botelho M, Cavadas C, Neuropeptide Y. Neuropeptide Y: An Anti-Aging Player? Trends Neurosci. 2015; 38:701–11. https://doi.org/10.1016/j.tins.2015.08.012 [PubMed]
  • 55. Banki E, Sosnowska D, Tucsek Z, Gautam T, Toth P, Tarantini S, Tamas A, Helyes Z, Reglodi D, Sonntag WE, Csiszar A, Ungvari Z. Age-related decline of autocrine pituitary adenylate cyclase-activating polypeptide impairs angiogenic capacity of rat cerebromicrovascular endothelial cells. J Gerontol A Biol Sci Med Sci. 2015; 70:665–74. https://doi.org/10.1093/gerona/glu116 [PubMed]
  • 56. Reglodi D, Jungling A, Longuespée R, Kriegsmann J, Casadonte R, Kriegsmann M, Juhasz T, Bardosi S, Tamas A, Fulop BD, Kovacs K, Nagy Z, Sparks J, et al. Accelerated pre-senile systemic amyloidosis in PACAP knockout mice - a protective role of PACAP in age-related degenerative processes. J Pathol. 2018; 245:478–90. https://doi.org/10.1002/path.5100 [PubMed]
  • 57. Beck B. [Feeding regulatory peptides: hopes and limits for the treatment of obesities]. J Soc Biol. 2006; 200:7–16. https://doi.org/10.1051/jbio:2006002 [PubMed]
  • 58. Park S, Komatsu T, Kim SE, Tanaka K, Hayashi H, Mori R, Shimokawa I. Neuropeptide Y resists excess loss of fat by lipolysis in calorie-restricted mice: a trait potential for the life-extending effect of calorie restriction. Aging Cell. 2017; 16:339–48. https://doi.org/10.1111/acel.12558 [PubMed]
  • 59. Chiba T, Tamashiro Y, Park D, Kusudo T, Fujie R, Komatsu T, Kim SE, Park S, Hayashi H, Mori R, Yamashita H, Chung HY, Shimokawa I. A key role for neuropeptide Y in lifespan extension and cancer suppression via dietary restriction. Sci Rep. 2014; 4:4517. https://doi.org/10.1038/srep04517 [PubMed]
  • 60. Horvath S, Oshima J, Martin G, Raj K, Matsuyama S. Epigenetic age estimator for skin and blood applied to Hutchinson Gilford Progeria. 2018. Aging (Albany NY). 2018; 10:1758–75. [PubMed]
  • 61. Lu AT, Xue L, Salfati EL, Chen BH, Ferrucci L, Levy D, Joehanes R, Murabito JM, Kiel DP, Tsai PC, Yet I, Bell JT, Mangino M, et al. GWAS of epigenetic aging rates in blood reveals a critical role for TERT. Nat Commun. 2018; 9:387. https://doi.org/10.1038/s41467-017-02697-5 [PubMed]
  • 62. Lu AT, Hannon E, Levine ME, Crimmins EM, Lunnon K, Mill J, Geschwind DH, Horvath S. Genetic architecture of epigenetic and neuronal ageing rates in human brain regions. Nat Commun. 2017; 8:15353. https://doi.org/10.1038/ncomms15353 [PubMed]
  • 63. Lu AT, Hannon E, Levine ME, Hao K, Crimmins EM, Lunnon K, Kozlenkov A, Mill J, Dracheva S, Horvath S. Genetic variants near MLST8 and DHX57 affect the epigenetic age of the cerebellum. Nat Commun. 2016; 7:10561. https://doi.org/10.1038/ncomms10561 [PubMed]
  • 64. Orozco LD, Farrell C, Hale C, Rubbi L, Rinaldi A, Civelek M, Pan C, Lam L, Montoya D, Edillor C, Seldin M, Boehnke M, Mohlke KL, et al. Epigenome-wide association in adipose tissue from the METSIM cohort. Hum Mol Genet. 2018; 27:2586–2586. https://doi.org/10.1093/hmg/ddy205 [PubMed]
  • 65. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, Yefanov A, Lee H, Zhang N, et al. NCBI GEO: archive for functional genomics data sets--update. Nucleic Acids Res. 2013; 41:D991–95. https://doi.org/10.1093/nar/gks1193 [PubMed]
  • 66. Orozco LD, Morselli M, Rubbi L, Guo W, Go J, Shi H, Lopez D, Furlotte NA, Bennett BJ, Farber CR, Ghazalpour A, Zhang MQ, Bahous R, et al. Epigenome-wide association of liver methylation patterns and complex metabolic traits in mice. Cell Metab. 2015; 21:905–17. https://doi.org/10.1016/j.cmet.2015.04.025 [PubMed]
  • 67. Cannon MV, Pilarowski G, Liu X, Serre D. Extensive Epigenetic Changes Accompany Terminal Differentiation of Mouse Hepatocytes After Birth. G3 (Bethesda). 2016; 6:3701–09. [PubMed]
  • 68. Reizel Y, Spiro A, Sabag O, Skversky Y, Hecht M, Keshet I, Berman BP, Cedar H. Gender-specific postnatal demethylation and establishment of epigenetic memory. Genes Dev. 2015; 29:923–33. https://doi.org/10.1101/gad.259309.115 [PubMed]
  • 69. Guo W, Fiziev P, Yan W, Cokus S, Sun X, Zhang MQ, Chen PY, Pellegrini M. BS-Seeker2: a versatile aligning pipeline for bisulfite sequencing data. BMC Genomics. 2013; 14:774. https://doi.org/10.1186/1471-2164-14-774 [PubMed]
  • 70. Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011; 27:1571–72. https://doi.org/10.1093/bioinformatics/btr167 [PubMed]
  • 71. Clopper CP, Pearson ES. E.S. The use of confidence or fiducial limits illustrated in the case of the binomial. Biometrika. 1934; 26:404–13. https://doi.org/10.1093/biomet/26.4.404
  • 72. Jeffreys H. An invariant form for the prior probability in estimation problems. Proc R Soc Lond A Math Phys Sci. 1946; 186:453–61. https://doi.org/10.1098/rspa.1946.0056 [PubMed]
  • 73. Friedman J, Hastie T, Tibshirani R. Regularization Paths for Generalized Linear Models via Coordinate Descent. J Stat Softw. 2010; 33:1–22. https://doi.org/10.18637/jss.v033.i01 [PubMed]
  • 74. Bennett BJ, Farber CR, Orozco L, Kang HM, Ghazalpour A, Siemers N, Neubauer M, Neuhaus I, Yordanova R, Guan B, Truong A, Yang WP, He A, et al. A high-resolution association mapping panel for the dissection of complex traits in mice. Genome Res. 2010; 20:281–90. https://doi.org/10.1101/gr.099234.109 [PubMed]
  • 75. Lusis AJ, Seldin MM, Allayee H, Bennett BJ, Civelek M, Davis RC, Eskin E, Farber CR, Hui S, Mehrabian M, Norheim F, Pan C, Parks B, et al. The Hybrid Mouse Diversity Panel: a resource for systems genetics analyses of metabolic and cardiovascular traits. J Lipid Res. 2016; 57:925–42. https://doi.org/10.1194/jlr.R066944 [PubMed]
  • 76. Reizel Y, Spiro A, Sabag O, Skversky Y, Hecht M, Keshet I, Berman BP, Cedar H. Gender-specific postnatal demethylation and establishment of epigenetic memory. Genes Dev. 2015; 29:923–33. https://doi.org/10.1101/gad.259309.115 [PubMed]