Epigenetic age-predictor for mice based on three CpG sites
- RWTH Aachen University Medical School, Germany
- University Hospital RWTH Aachen, Germany
- Ulm University, Germany
- University Medical Center Groningen, Netherlands
- Instituto Biodonostia, Spain
- Tecnun-University of Navarra, Spain
- Cincinnati Children's Hospital Burnet Campus, United States
Abstract
Epigenetic clocks for mice were generated based on deep-sequencing analysis of the methylome. Here, we demonstrate that site-specific analysis of DNA methylation levels by pyrosequencing at only three CG dinucleotides (CpGs) in the genes Prima1, Hsf4, and Kcns1 facilitates precise estimation of chronological age in murine blood samples, too. DBA/2 mice revealed accelerated epigenetic aging as compared to C57BL6 mice, which is in line with their shorter life-expectancy. The three-CpG-predictor provides a simple and cost-effective biomarker to determine biological age in large intervention studies with mice.
https://doi.org/10.7554/eLife.37462.001eLife digest
Epigenetic marks are chemical modifications found throughout the genome – the DNA within cells. By influencing the activity of nearby genes, the marks govern developmental processes and help cells to adapt to changes in their surroundings. Some epigenetic marks can be gained or lost with age. A lot of aging research focuses on one type of mark, called “DNA methylation”. By measuring the presence or absence of specific methyl groups, scientists can estimate biological age – which may differ from calendar age.
Recent studies have developed computer models called epigenetic aging clocks to predict the biological age of mouse cells. These clocks use epigenetic data collected from the entire genomes of mice, and are useful for understanding how the aging process is affected by genetic parameters, diet, or other environmental factors. Yet, the genome sequencing methods used to construct most existing epigenetic clocks are expensive, labor-intensive, and cannot be easily applied to large groups of mice.
Han et al. have developed a new way to predict biological aging in mice that needs methylation information from just three particular sections of the genome. Even though this approach is much faster and less expensive than other epigenetic approaches to measuring aging, it has a similar level of accuracy to existing models. Han et al. use the new method to show that cells from different strains of laboratory mice age at different rates. Furthermore, in a strain that has a shorter life expectancy, aging seems to be accelerated.
The new approach developed by Han et al. will make it easier to study how aging in mice is affected by different interventions. Further studies will also be needed to better understand how epigenetic marks relate to biological aging.
https://doi.org/10.7554/eLife.37462.002Introduction
Age-associated DNA methylation (DNAm) was first described for humans after Illumina Bead Chip microarray data became available to enable cross comparison of thousands of CpG loci (Bocklandt et al., 2011; Koch and Wagner, 2011). Many of these age-associated CpGs were then integrated into epigenetic age-predictors (Hannum et al., 2013; Horvath, 2013; Weidner et al., 2014). However, site-specific DNAm analysis at individual CpGs can also provide robust biomarkers for aging. For example, we have described that DNAm analysis at only three CpGs enables age-predictions for human blood samples with a mean absolute deviation (MAD) from chronological age of less than five years (Weidner et al., 2014). Such simplistic age-predictors for human specimen are widely used because they enable fast and cost-effective analysis in large cohorts.
Recently, epigenetic clocks were also published for mice by using either reduced representation bisulfite sequencing (RRBS) or whole genome bisulfite sequencing (WGBS) (Petkovich et al., 2017; Stubbs et al., 2017; Wang et al., 2017). For example, Petkovich et al. described a 90 CpG model for blood (Petkovich et al., 2017), and Stubbs and coworkers a 329 CpG model for various different tissues (Stubbs et al., 2017). Nutrition and genetic background seem to affect the epigenetic age of mice – and thereby possibly aging of the organism (Cole et al., 2017; Hahn et al., 2017; Maegawa et al., 2017). In analogy, epigenetic aging of humans is associated with life expectancy, indicating that it rather reflects biological age than chronological age (Lin et al., 2016; Marioni et al., 2015). However, DNAm profiling by deep sequencing technology is technically still challenging, relatively expensive, and not every sequencing-run covers all relevant CpG sites with enough reading depth.
Results
Therefore, we established pyrosequencing assays for nine genomic regions of previously published predictors (Petkovich et al., 2017; Stubbs et al., 2017). These regions were preselected to have multiple age-associated CpGs in close vicinity. DNAm was then analyzed in 24 blood samples of female C57BL/6 mice that covered a broad range of 12 different age groups (11 to 117 weeks old). The nine amplicons covered a total of 71 CpG sites (Supplementary file 1) and we used machine learning to identify the best fitted model for epigenetic age-predictions using cross-fold validation on the training set. The best results were observed for 15 CpGs from five different amplicons that provided an extremely high correlation with chronological age in the training set (R2 = 0.99; mean absolute deviation [MAD]=2.76 weeks; Supplementary file 2), albeit the training set might be too small for this approach. To make the method more easily applicable and more cost-effective, we wanted to focus on less CpGs. When we varied the regularization parameters for models with less CpGs, the precision declined significantly. For example the best model with three CpGs comprised the three CpGs of Hsf4 (CpGs# 3,4,5) that also revealed the overall highest Pearson correlations with chronological age (R2 = 0.95; MAD = 5.24 weeks). However, combination of different hypo- and hypermethylated amplicons might be advantageous to facilitate better assessment of plausibility of the results. Therefore, we alternatively selected those three CpGs that revealed the highest Pearson correlation with chronological age in different amplicons. These three CpGs were associated with the genes Proline rich membrane anchor 1 (Prima1: chr12:103214639; R2 = 0.71), Heat shock transcription factor 4 (Hsf4: chr8:105271000; R2 = 0.95) and Potassium voltage-gated channel modifier subfamily S member 1 (Kcns1: chr2:164168110; R2 = 0.83; Figure 1A–C; Figure 1—figure supplement 1). Notably, all three CpGs were derived from the epigenetic age-predictor for blood samples (Petkovich et al., 2017). A multivariable model for age-predictions was established for DNAm at the CpGs in Prima 1 (α), Hsf4 (β), and Kcns1 (γ):
Predicted ageC57BL/6 (in weeks) = −58.076 + 0.25788 α + 3.06845 β + 1.00879 γ
Age-predictions correlated very well with the chronological age of C57BL/6 mice in the training set (R2 = 0.96; MAD = 4.86 weeks; Figure 1D).
Figure 1 with 1 supplement see all
Three CpG epigenetic age-predictor for mice.
(a–c) DNA methylation (DNAm) of three CpGs in the genes Prima1, Hsf4 and Kcns1 was analyzed by pyrosequencing in 24 C57BL/6 mice (training set). Coefficient of determination (R2) of DNAm versus chronological age is indicated. (d) Based on these age-associated DNAm changes a multivariable model for age prediction was calculated. (e–g) Subsequently, two independent validation sets were analyzed: 21 C57BL/6 mice from the University of Ulm and 19 C57BL/6 mice from the University of Groningen (validation sets 1 and 2, respectively). (h) Age predictions with the three-CpG-model revealed a high correlation with chronological age in the independent validation sets (MAD = mean absolute deviation; MAE = median absolute error).
Our three CpG age-predictor was subsequently validated in a blinded manner for 21 C57BL/6J mice (7 to 104 weeks old) from the University of Ulm (validation set 1) and 19 C57BL/6J mice (14 to 109 weeks old) from the University of Groningen (validation set 2). The results of both validation sets revealed high correlations with chronological age (R2 = 0.95 and 0.91, respectively; Figure 1E–H) with relatively small MADs (6.9 and 7.1 weeks) and median absolute errors (MAE; 5.0 and 5.9 weeks). Thus, our age-predictions seem to have similar precision as previously described for multi-CpG predictors based on RRBS or WGBS data (Petkovich et al., 2017; Stubbs et al., 2017; Wang et al., 2017).
Gender did not have significant impact on our epigenetic age-predictions for mice (Figure 2), as described before (Maegawa et al., 2017; Petkovich et al., 2017; Stubbs et al., 2017). In contrast, the human epigenetic clock is clearly accelerated in male donors (Hannum et al., 2013; Horvath, 2013; Weidner et al., 2014). This coincides with shorter life expectancy in men than woman, whereas in mice there are no consistent sex differences in longevity (Goodrick, 1975).
Figure 2
Gender does not affect epigenetic age predictions in mice.
The deviations of predicted age by our three-CpG predictor versus chronological age did not reveal significant differences between female and male C57BL/6 mice (Mann–Whitney U test p=0.6).
To address the question if our three CpG signature was also applicable for other tissues than blood we analyzed the DNAm in skin, kidney, intestine, lung, liver, heart, brain, testis, and pancreas of 3 young (9.6 weeks old) and three old mice (56.9 weeks old). In all tissues tested the samples of old mice were predicted to be older using our three CpG signature. However, the different DNAm levels clearly demonstrate that the model needs to be retrained to be applied for these tissues (Figure 3).
Figure 3
Age-associated DNA methylation at the three CpG sites in different tissues.
Different tissues were isolated of three young (9.6 weeks) and three old mice (56.9 weeks) and DNAm was analyzed at the three relevant CpGs in (a) Prima1, (b) Hsf4, and (c) Kcns1. Epigenetic age-predictions using the 3 CpG model for blood demonstrated also significant differences between young and old mice in skin, intestine, brain, and testis (mean ± standard deviation; Student t-tests: *p<0.05; **p<0.01; ***p<0.001).
Subsequently, we analyzed epigenetic aging of DBA/2 mice that have a shorter life expectancy than C57BL/6 mice (Goodrick, 1975) (33 mice from Ulm and Groningen; 6 to 109 weeks old). The three CpGs in Prima1, Hsf4 and Kcns1 revealed high correlation with chronological age (R2 = 0.91, 0.88 and 0.83, respectively), albeit the offset in DNAm between DBA/2 and C57BL/6 mice indicated that the signature needs to be retrained for different mouse strains (Figure 4a–c). Notably, the slopes were higher in DBA/2 mice, particularly for the CpG in Prima1. Furthermore, DNAm of Hsf4 increased at a higher rate in young DBA/2 mice, indicating that it is more accurately modelled as a function of logarithmic age. This has also been described in human for many age-associated CpGs in pediatric cohorts (Alisch et al., 2012). In fact, epigenetic age-predictions in DBA/2 mice seemed to follow a logarithmic model of age (R2 = 0.89; Figure 4d) rather than a linear association (R2 = 0.86). These results provided evidence for accelerated epigenetic aging of DBA/2 mice.
Figure 4
Epigenetic aging is accelerated in DBA/2 mice as compared to C57BL/6 mice.
(a–c) Age-related DNA methylation (DNAm) determined by pyrosequencing assay for three candidate CpGs on 33 of DBA/2 blood samples (14 mice from the University of Ulm and 19 mice from the University of Groningen; red). For comparison we provided measurements of the C57BL/6 mice (only from validation sets; blue). (d) Epigenetic age-predictions using the three CpG multivariable model for the C57BL/6 mice (blue; linear regression) and DBA/2 mice (red, logarithmic regression). Age-predictions in DBA/2 mice rather followed a logarithmic regression (R = Pearson correlation); (e) Based on the DNAm measurements in DBA/2 we adjusted the multivariate regression model for age-predictions of this mouse strain as described in the text (DBA/2 predictor).
Either way, epigenetic age-predictions were overall significantly overestimated in the shorter-lived DBA/2 mice, suggesting that age-predictors need to be adjusted for different inbreed mice strains. To this end, we have retrained a multivariate model for DBA/2 mice:
Predicted ageDBA/2 (in weeks) = 87.54294–1.22221 α + 0.991558 β + 0.355444 γ
This adjusted model facilitated relatively precise age-predictions for DBA/2 mice (R2 = 0.95; MAD = 7.1 weeks; MAE = 5.3 weeks; Figure 4e).
Discussion
Generation of confined epigenetic signatures is always a tradeoff between integrating more CpGs for higher precision and higher costs for analysis (Wagner, 2017). It was somewhat unexpected that with only three CpGs our signature facilitated similar precision of epigenetic age-predictions as the previously published signatures based on more than 90 CpGs. This can be attributed to the higher precision of DNAm measurements at individual CpGs by bisulfite pyrosequencing, which is one of the most precise methods for determining DNAm at single CpG resolution (BLUEPRINT consortium, 2016). Particularly in RRBS data not all CpG sites are covered in all samples and a limited number of reads notoriously entails lower precision of DNAm levels at these genomic locations. Thus, genome wide deep sequencing approaches facilitate generation of robust large epigenetic age-predictors, while site specific analysis may compensate by higher precision of DNAm measurement at individual CpGs.
The ultimate goal of epigenetic age-predictors for mice is not to develop near perfect age predictors, but to provide a surrogate for biological aging that facilitates assessment of interventions on aging. In fact, using deep sequencing approaches (RRBS or WGBS) several groups already indicated that relevant parameters that affect aging of the organism - such as diet, genetic background, and drugs - do also impact on epigenetic aging (Cole et al., 2017; Hahn et al., 2017; Maegawa et al., 2017). It is yet unclear if epigenetic aging signatures can be specifically trained to either correlate with chronological age or biological age. For humans, recent studies indicate that this might be possible (Levine et al., 2018) and we have previously demonstrated that even individual age-associated CpGs can be indicative for life expectancy (Zhang et al., 2017). Further studies will be necessary to gain better understanding how epigenetic age predictions are related to the real state of biological aging, and how it is related to alternative approaches to quantify biological aging, such as telomere length (Belsky et al., 2018).
Our three CpG model has been trained for blood samples – a specimen that is commonly used in biochemical analysis and the small required volume can be taken without sacrificing the mice. However, epigenetic aging may occur at different rates in different tissues. It is difficult to address this question in humans because it is difficult to collect samples of various tissues in large aging cohorts, whereas this is feasible in mice. We demonstrate that age-associated DNAm changes occur in multiple tissues in our three CpGs albeit they were initially identified in blood (Petkovich et al., 2017). Furthermore, DNAm levels may vary between different hematopoietic subsets (Frobel et al., 2018; Houseman et al., 2014). In the future, sorted subsets should be analyzed to determine how the three CpG signature is affected by blood counts.
The results of our three CpG signature suggest that epigenetic aging is accelerated in DBA/2 mice. Notably, in elderly DBA/2 mice the epigenetic age predictions revealed higher ‘errors’ from chronological age, which might be attributed to the fact that the variation of lifespan is higher in DBA/2 than C57BL/6 mice (de Haan et al., 1998; Goodrick, 1975). It will be important to validate the association of the epigenetic age-predictions with biological age by additional correlative studies, including life expectancy in mice.
Taken together, we describe an easily applicable but quite precise approach to determine epigenetic age of mice. We believe that our assay will be instrumental to gain additional insight into mechanisms that regulate age-associated DNAm and for longevity intervention studies in mice.
Materials and methods
Mouse strains and blood collection
Request a detailed protocolBlood samples of C57BL/6J mice of the training set and of the validation set one were taken at the University of Ulm by submandibular bleeding (100–200 μl) of living mice or postmortem from the vena cava. C57BL/6J samples of the validation set two were taken at the University of Groningen from the cheek. DBA/2J samples were taken at the University of Ulm (n = 14) and Groningen (n = 19). All mice were accommodated under pathogen-free conditions. Experiments were approved by the Institutional Animal Care of the Ulm University as well as by Regierungspräsidium Tübingen and by the Institutional Animal Care and Use Committee of the University of Groningen (IACUC-RUG), respectively. To analyze age-associated changes in different tissues we used three young (9.6 weeks old) and three old mice (56.9 weeks old) C57BL/6J mice (JaxMice) in accordance with relevant Spanish and European guidelines after approval by the Biodonostia Animal Care Committee. These mice were sacrificed and dissected immediately. 25 mg of tissue (10 mg in the case of spleen) or 200 µl of blood were used for DNA extraction.
Genomic DNA isolation and bisulfite conversion
Request a detailed protocolGenomic DNA was isolated from 50 µl blood using the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany). Kidney and liver DNA extractions were digested with Ribonuclease A (100 mg/ml, Sigma R4875). DNA concentration was quantified by Nanodrop 2000 Spectrophotometers (Thermo Scientific, Wilmington, USA). 200 ng of genomic DNA was subsequently bisulfite-converted with the EZ DNA Methylation Kit (Zymo Research, Irvine, USA).
Pyrosequencing
Request a detailed protocolBisulfite converted DNA was subjected to PCR amplification. Primers were purchased at Metabion and the sequences are provided in Supplementary file 3. 20 µg PCR product was immobilized to 5 µl Streptavidin Sepharose High Performance Bead (GE Healthcare, Piscataway, NJ, USA), and then annealed to 1 µl sequencing primer (5 μM) for 2 min at 80°C. Amplicons were sequenced on PyroMark Q96 ID System (Qiagen, Hilden, Germany) and analyzed with PyroMark Q CpG software (Qiagen).
Alternative approaches to select CpGs for multivariable models
Request a detailed protocolWe used a penalized regression model from the R package glmnet on the training dataset to establish a predictor of mouse age based on CpG methylation. The alpha parameter of glmnet was set to 1 (lasso regression) and the lambda parameter was chosen by cross-fold validation of the training dataset (10-fold cross validation). Alternatively, we trained our multivariable model with preselected CpGs based on location in three different amplicons, high Pearson correlation (R) of DNAm with chronological age, and combination of hyper- and hypomethylated sites.
Statistical analysis
Request a detailed protocolLinear regressions, MAD and MAE were calculated with Excel. Statistical significance of the deviations between predicted and chronological age was estimated by Mann–Whitney U test or Student´s t-test as indicated.
Data availability
Raw data of pyrosequencing is provided as supplemental EXCEL table (Source data 1).
References
-
Age-associated DNA methylation in pediatric populationsGenome Research 22:623–632.https://doi.org/10.1101/gr.125187.111
-
Eleven telomere, epigenetic clock and biomarker-composite quantifications of biological aging: do they measure the same thing?American Journal of Epidemiology 187:1220–1230.https://doi.org/10.1093/aje/kwx346
-
Leukocyte counts based on DNA methylation at individual cytosinesClinical Chemistry 64:566–575.https://doi.org/10.1373/clinchem.2017.279935
-
Life-span and the inheritance of longevity of inbred miceJournal of Gerontology 30:257–263.https://doi.org/10.1093/geronj/30.3.257
-
Caloric restriction delays age-related methylation driftNature Communications 8:539.https://doi.org/10.1038/s41467-017-00607-3
-
Epigenetic aging clocks in mice and menGenome Biology 18:107.https://doi.org/10.1186/s13059-017-1245-8
Article and author information
Author details
Funding
Else Kröner-Fresenius-Stiftung (2014_A193)
- Wolfgang Wagner
Deutsche Forschungsgemeinschaft (WA 1706/8-1)
- Wolfgang Wagner
Bundesministerium für Bildung und Forschung (01KU1402B)
- Wolfgang Wagner
NIH Clinical Center (R01HL134617)
- Hartmut Geiger
Nederlandse Organisatie voor Wetenschappelijk Onderzoek
- Gerald de Haan
Deutsche Forschungsgemeinschaft (GRK 1789 CEMMA)
- Hartmut Geiger
Deutsche Forschungsgemeinschaft (GRK 2254 HEIST)
- Hartmut Geiger
Deutsche Forschungsgemeinschaft (SFBs 1074)
- Hartmut Geiger
Deutsche Forschungsgemeinschaft (SFBs 1149)
- Hartmut Geiger
Deutsche Forschungsgemeinschaft (SFBs 1275)
- Hartmut Geiger
NIH Clinical Center (R01DK104814)
- Hartmut Geiger
Bundesministerium für Bildung und Forschung (SyStarR)
- Hartmut Geiger
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by the Else Kröner-Fresenius-Stiftung (2014_A193; to WW), by the German Research Foundation (DFG; WA 1706/8-1 to WW; and GRK 1789 CEMMA, GRK 2254 HEIST and SFBs 1074, 1149 and 1275 to HG), by the German Ministry of Education and Research (BMBF; 01KU1402B to WW; and SyStarR to HG), and by the NIH (R01HL134617 and R01DK104814 to HG). The Groningen samples were obtained from the Mouse Clinic for Cancer and Ageing (http://www.mccanet.nl), which is supported by a grant from the Netherlands Organization for Scientific Research (NWO). The funding bodies were not involved in study design, data analysis, or writing of the manuscript.
Ethics
Animal experimentation: Experiments were approved by the Institutional Animal Care of the Ulm University as well as by Regierungspräsidium Tübingen and by the Institutional Animal Care and Use Committee of the University of Groningen (IACUC-RUG), respectively. To analyze age-associated changes in different tissues we used 3 young (67 days old) and 3 old (398 days old) C57BL/6J mice (JaxMice) in accordance with relevant Spanish and European guidelines after approval by the Biodonostia Animal Care Committee.
Version history
- Received: May 3, 2018
- Accepted: August 23, 2018
- Accepted Manuscript published: August 24, 2018 (version 1)
- Version of Record published: September 25, 2018 (version 2)
Copyright
© 2018, Han et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 5,283
- views
-
- 748
- downloads
-
- 49
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
A two-part list of links to download the article, or parts of the article, in various formats.
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Epigenetic age-predictor for mice based on three CpG sites
eLife 7:e37462.
https://doi.org/10.7554/eLife.37462
Further reading
-
- Developmental Biology
Inhibitory G alpha (GNAI or Gαi) proteins are critical for the polarized morphogenesis of sensory hair cells and for hearing. The extent and nature of their actual contributions remains unclear, however, as previous studies did not investigate all GNAI proteins and included non-physiological approaches. Pertussis toxin can downregulate functionally redundant GNAI1, GNAI2, GNAI3, and GNAO proteins, but may also induce unrelated defects. Here, we directly and systematically determine the role(s) of each individual GNAI protein in mouse auditory hair cells. GNAI2 and GNAI3 are similarly polarized at the hair cell apex with their binding partner G protein signaling modulator 2 (GPSM2), whereas GNAI1 and GNAO are not detected. In Gnai3 mutants, GNAI2 progressively fails to fully occupy the sub-cellular compartments where GNAI3 is missing. In contrast, GNAI3 can fully compensate for the loss of GNAI2 and is essential for hair bundle morphogenesis and auditory function. Simultaneous inactivation of Gnai2 and Gnai3 recapitulates for the first time two distinct types of defects only observed so far with pertussis toxin: (1) a delay or failure of the basal body to migrate off-center in prospective hair cells, and (2) a reversal in the orientation of some hair cell types. We conclude that GNAI proteins are critical for hair cells to break planar symmetry and to orient properly before GNAI2/3 regulate hair bundle morphogenesis with GPSM2.
-
- Computational and Systems Biology
- Developmental Biology
Organisms utilize gene regulatory networks (GRN) to make fate decisions, but the regulatory mechanisms of transcription factors (TF) in GRNs are exceedingly intricate. A longstanding question in this field is how these tangled interactions synergistically contribute to decision-making procedures. To comprehensively understand the role of regulatory logic in cell fate decisions, we constructed a logic-incorporated GRN model and examined its behavior under two distinct driving forces (noise-driven and signal-driven). Under the noise-driven mode, we distilled the relationship among fate bias, regulatory logic, and noise profile. Under the signal-driven mode, we bridged regulatory logic and progression-accuracy trade-off, and uncovered distinctive trajectories of reprogramming influenced by logic motifs. In differentiation, we characterized a special logic-dependent priming stage by the solution landscape. Finally, we applied our findings to decipher three biological instances: hematopoiesis, embryogenesis, and trans-differentiation. Orthogonal to the classical analysis of expression profile, we harnessed noise patterns to construct the GRN corresponding to fate transition. Our work presents a generalizable framework for top-down fate-decision studies and a practical approach to the taxonomy of cell fate decisions.
