A whole lifespan mouse multi-tissue DNA methylation clock
Figures

Construction of the WLMT clock.
(A) Datasets used for the clock construction and initial validation. Radii of pie charts correspond to the number of samples in a group. The range of ages is shown above (in months). Performance of the whole-lifespan multi-tissue clock on training (B) and test (C) sets. MAE on training and test sets was 28.6 and 72.7 days, respectively. Tissues used are indicated on the right according to the color scheme.

Density plot for the number of CpG sites covered in all samples.
(A) and in samples with covered sites, which were used to train the WLMT clock (B).

Principle component analysis of samples sequenced in the study.
Library ‘M20’ was not used in the preparation of the clock.

Performance of the WLMT clock on male (blue) and female (pink) samples.
https://doi.org/10.7554/eLife.40675.006
Number of reads obtained for the sequenced samples.
https://doi.org/10.7554/eLife.40675.007
Overlap in CpG sites among mouse DNAm clocks.
(A) overlap in CpG sites. (B) overlap in genes containing the clock sites. * - significant with (two-tailed Chi-square test with Yates correction).

Distribution of WLMT clock CpG sites along the mouse genome.
Sites with positive weights (increasing methylation with age) are shown in green, and those with negative weights (decreasing methylation with age) in red.

Location of CpG sites forming the four epigenetic clocks in the mouse genome: WLMT sites are shown in blue, YOMT sites in green, blood clock sites in red and liver clock sites in orange.
https://doi.org/10.7554/eLife.40675.010
Frequency of the occurrence of WLMT clock sites in 100 clocks generated based on the robustness study.
100 clocks are trained on different subsets of the set of samples used in WLMT construction. (A) Density plot showing the overall distribution of frequencies of WLMT clock sites as found in the generated 100 clocks. Carpet plot shows the exact values of frequencies. (B) Relation of the absolute values of weight of sites in the WLMT clock and frequency of these sites in the 100 alternative clocks.

Age-related changes in DNA methylation of the ten WLMT clock sites with highest positive weights.
Tissues used in construction of the clock are shown.

Age-related changes in DNA methylation for the ten WLMT clock sites with lowest negative weights.
https://doi.org/10.7554/eLife.40675.013
Density plots showing distribution of weights of CpG clock sites in the four mouse clocks.
https://doi.org/10.7554/eLife.40675.014
Absolute errors of four mouse DNAm clocks applied to untreated wild-type samples.
(A) In this case, WLMT clock was applied to its test set, whereas for the other age predictors this set also included their training set samples. (B) Blood samples which were not used in any of the clocks. (C) Only liver samples. (D) Samples of all tissues, excluding those used in the YOMT clock. (E) Samples of all non-blood tissues. (F) All tissues, excluding liver.

Performance of the WLMT, blood, YOMT and liver clocks on the liver samples not included into the WLMT and liver training sets.
∆Age shows difference between DNAm age and chronological age. Chronological age of the included samples is 180 days. MAE is 244 days for the WLMT and 65 days for the liver clock.

Performance of the WLMT, blood, YOMT and liver clocks.
WLMT: only test samples are included. YOMT: test and training sets are shown. Blood samples have not been used in development of this clock. Blood clock: all blood samples belong to the training set, while all other tissues to the test set. Liver clock: only test samples are shown.

DNAm ages estimated for an independent dataset.
An additional dataset (Reizel et al., 2018) was tested with the four clocks. Solid line represents a perfect correspondence between DNAm age and chronological age. Dashed lines show slopes of linear regressions made for each DNAm clock based on the age values estimated with it. Slopes: 0.14 (YOMT), 0.16 (liver), 1.24 (WLMT), 0.45 (blood). p-values: (YOMT), (liver), (WLMT), (blood).

Difference in DNAm age estimated by four age predictors and chronological age based on WLMT, blood, YOMT and liver clocks.
Untreated and wild-type samples are shown in green, and intervention samples in orange. (A) Effect of caloric restriction in mouse strain C57BL/6, DNA from blood. (B) Effect of caloric restriction, mouse strain B6D2F1, DNA from blood. (C) Conversion of fibroblasts to iPSCs, DNA from cultured cells. (D) Growth hormone receptor knockout, DNA from blood. (E) Snell dwarf mice, DNA from blood. (F) Ames dwarf mice, DNA from liver. ΔAge was defined as the difference between the DNAm age and the chronological age. *-, **-, ***-.
Tables
Summary of available mouse DNAm age predictors.
https://doi.org/10.7554/eLife.40675.002Tissue type | Reported precision, R2 | Age range of training set, months | Number of CpG sites in the model | Publication |
---|---|---|---|---|
Blood | >0.9008 | 3.–35 | 90 | (Petkovich et al., 2017) |
Liver | 0.91 | 0.2–26 | 148* | (Wang et al., 2017) |
Multi-tissue | 0.7 | 0.2–9.5 | 329 | (Stubbs et al., 2017) |
Multi-tissue | 0.89 | 0.2–35 | 435 | This paper |
-
*Unlike other DNAm clocks, the liver clock is based on a combination of methylation levels on both strands and thus it corresponds to 296 genomic positions.
Additional files
-
Supplementary file 1
Sample information.
- https://doi.org/10.7554/eLife.40675.020
-
Supplementary file 2
Information on samples sequenced by RRBS in the study and included in the analysis.
- https://doi.org/10.7554/eLife.40675.021
-
Supplementary file 3
Mouse DNAm clocks.
- https://doi.org/10.7554/eLife.40675.022
-
Supplementary file 4
Genes overlapping in pairs of clocks.
- https://doi.org/10.7554/eLife.40675.023
-
Supplementary file 5
Distribution of samples across libraries.
- https://doi.org/10.7554/eLife.40675.024