Abstract
DNA methylation (DNAm) has been revealed to play a role in various diseases. Here we performed epigenome-wide screening and validation to identify mortality-related DNAm signatures in a general population-based cohort with up to 14 years follow-up. In the discovery panel in a case-cohort approach, 11,063 CpGs reach genome-wide significance (FDR<0.05). 58 CpGs, mapping to 38 well-known disease-related genes and 14 intergenic regions, are confirmed in a validation panel. A mortality risk score based on ten selected CpGs exhibits strong association with all-cause mortality, showing hazard ratios (95% CI) of 2.16 (1.10–4.24), 3.42 (1.81–6.46) and 7.36 (3.69–14.68), respectively, for participants with scores of 1, 2–5 and 5+ compared with a score of 0. These associations are confirmed in an independent cohort and are independent from the ‘epigenetic clock’. In conclusion, DNAm of multiple disease-related genes are strongly linked to mortality outcomes. The DNAm-based risk score might be informative for risk assessment and stratification.
Similar content being viewed by others
Introduction
DNA methylation (DNAm), as the most widely studied form of epigenetic programming, has been revealed to be modulated by lifestyle and environmental factors1,2 and to be involved in onset and progression of complex diseases, including various forms of malignant diseases, cardiovascular diseases (CVDs), metabolic diseases (for example, diabetes), neuropsychiatric disorders and autoimmune disorders3,4,5,6,7. Therefore, DNAm could plausibly be associated with the excess mortality from specific diseases and consequently with all-cause mortality. This was exemplified by the previous investigations on smoking-associated DNAm changes and their relationship with lung cancer incidence/mortality and mortality from any cause, cancer and CVD8,9,10.
In addition, evidence has accumulated that the recently established ‘epigenetic clock’ (also known as DNAm age) based on age-associated DNAm changes, which presumably reflects individuals’ biological age, is indicative for ageing-related outcomes and longevity11,12,13,14. Following the first study reporting an association of DNAm age with all-cause mortality by Marioni et al.13, the association was consistently demonstrated in various longitudinal studies15,16, for individual age-associated CpGs17 and also for newly identified age-associated CpGs18. On the other hand, several epigenome-wide association studies (EWASs) have pointed out that DNAm involved in ageing-related phenotypes are largely distinct from the established age-associated DNAm19,20,21.
To unravel the determinants of survival in the DNAm landscape, we performed an epigenome-wide screening and replication for mortality-related DNAm signatures in a general population-based cohort of older adults. Here we show that DNAm of 58 CpGs in baseline blood samples are associated with mortality from any causes during 14 years of follow-up. A mortality risk score based on ten selected CpGs strongly predicts all-cause, CVD and cancer mortality, also in an independent population-based cohort. The identified DNAm markers may thus bear implications in risk assessment and stratification in clinical practice.
Results
Study population
Table 1 presents the baseline characteristics of the ESTHER (Epidemiologische Studie zu Chancen der Verhütung, Früherkennung und optimierten Therapie chronischer Erkrankungen in der älteren Bevölkerung) study population. Of the 406 deaths in the case–cohort sample of the discovery panel, 90 were also included in the subcohort owing to random selection of subcohort at baseline. The time between blood sample collection and death ranged from 0.2 to 12.3 years (median (interquartile range (IQR), 7.4 (4.5–9.6) years) for these 406 participants. The corresponding figures for the 231 deaths in the validation panel were 0.2–13.8 years (range) and 8.6 (5.6–11.6) years (median (IQR)). The characteristics of the participants in the subcohort of the discovery panel are similar as those of the participants in the validation panel, except that the proportion of women was larger in the subcohort than in the validation cohort. In comparison with those two subgroups, the group of deceased participants in the discovery panel featured higher proportions of men, smokers, old (>70 years) and inactive participants, and participants with prevalence of hypertension, diabetes, CVD and cancer at baseline. The characteristics of the KORA (Kooperative Gesundheitsforschung in der Region Augsburg) study population are presented in Supplementary Table 1. The average age was similar in KORA and ESTHER participants (61 versus 62 years), but KORA participants had a much broader age range (31–82 years) than ESTHER participants (50–75 years).
Discovery and validation of mortality-related CpGs
In the discovery phase, a total of 11,063 CpGs passed the genome-wide significance threshold (false discovery rate (FDR) <0.05) (Supplementary Fig. 1). Associations with all-cause mortality were successfully replicated for 58 CpGs even after comprehensive confounder adjustment in the validation phase. Manhattan plots for the discovery and validation analyses are presented in Supplementary Fig. 2. Table 2 shows the results for the 58 CpGs. Methylation at the vast majority (49 of 58 CpGs) was inversely associated with mortality, with hazard ratios (HRs) and 95% confidence intervals (95% CIs) for a decrease in methylation by 1 s.d. ranging from 1.16 (1.04–1.28) to 1.95 (1.29–2.94). HRs (95% CI) for the other 9 CpGs showing positive associations with mortality ranged from 0.60 (0.47–0.77) to 0.83 (0.71–0.97) per s.d. decrease in methylation. The 58 loci are located at 38 genes and 14 intergenic regions across 19 chromosomes. In addition to three CpGs within AHRR, ten clusters within the identified sites were observed (Table 2), that is, 1p21.2 (2 CpGs), 2q37 (2 CpGs), 3q11/12 (2 CpGs), 6p21 (4 CpGs), 11p15 (3 CpGs), 11q13 (3 CpGs), 17q21 (2 CpGs), 17q25 (2 CpGs), 19p13 (3 CpGs) and 19q13 (7 CpGs). A literature search in PubMed for genes containing the identified CpGs found evidence that these genes or their methylation are involved in a variety of major diseases, including diabetes (for example, SARS, SQLE, NFE2L3, KCNQ1OT1 and SOCS3), CVD (for example, SARS, VCAM1, PLCL2, UTS2D, AHRR, 6p21.33, SQLE, KCNQ1OT1, SEMA7A, F2RL3, BCL3, PPP1R15A, PDE9A and MIR19A), various forms of cancers (for example, SOCS3, SLC1A5, MIR19A, MIR10A, CALR, ERCC1, BCL3, SQLE, RARA, LAPTM5, INPP5A, CSGALNACT1, KCNQ1OT1, CDC42BPB, PDE9A and MKL1), neuropsychiatric disorders (FOSL2, ATL3, SHANK2 and PPP1R15A) and HIV infection (for example, GPR15 and MIR10A) (Table 2 and Supplementary Table 2). Several of those genes, such as SQLE, KCNQ1OT1 and SOCS3, have been suggested to play roles in multiple types of diseases. Means and s.d. of the 58 CpGs at baseline among deceased participants and survivors are illustrated in Fig. 1.
(a) Mean and s.d. (error bar) of 22 mortality-related CpGs (also discovered to be associated with smoking in both current and previous studies) by vital status; (b) mean and s.d. (error bar) of 26 mortality-related CpGs (also discovered to be associated with smoking in the current study) by vital status; (c) mean and s.d. (error bar) of other 10 mortality-related CpGs by vital status.
Associations of risk factors with mortality-related CpGs
In the analyses of associations between the 58 CpGs and the covariates, differences in methylation levels with respect to age and sex were observed for 23 and 25 CpGs, respectively (Supplementary Table 3). However, none of the 58 CpGs overlapped with previously identified ageing-related sites11,12,18,22,23. Forty-eight of the 58 CpGs were differentially methylated according to smoking exposure and 22 of the 48 CpGs had also been found to be associated with smoking by previous EWASs2,24 (CpGs displayed in bold in Table 2). Five of the 48 smoking-associated CpGs and cg24397007 in FOSL2 were also associated with alcohol consumption (Supplementary Table 3). Four of the 48 smoking-associated CpGs and cg08362785 in MKL1 were also associated with prevalent diabetes; of these 5 sites, cg18181703 in SOCS3 was also recently identified to be associated with type 2 diabetes (T2D)5,25 and cg23190089 is located at SLC22A18AS, a locus near to known methylation-regulated genes implicated in T2D26,27. In addition, 4 of the 48 smoking-associated CpGs, including 2 diabetes-associated sites (cg18181703 in SOCS3 and cg26470501 in BLC3), were also associated with prevalent cancer. An illustration of the 48 CpGs is presented in Supplementary Fig. 3.
Mortality risk score and validation
Ten CpGs (cg01612140, cg05575921, cg06126421, cg08362785, cg10321156, cg14975410, cg19572487, cg23665802, cg24704287 and cg25983901) were selected by least absolute shrinkage and selection operator (LASSO) regression. Preliminary analyses in ESTHER samples showed that ⩾40% deaths occurred among participants with methylation levels in the highest quartile of cg08362785 (hypermethylated among deaths) or in the first quartile of the other 9 CpGs (demethylated among deaths) (Supplementary Fig. 4a). We therefore used the fourth quartile value of cg08362785 and first quartile values of other nine CpGs as the cutoff points, to define aberrant methylation for each CpG (the exact cutoff points are listed in Supplementary Table 4). Participants with aberrant methylation at 1–10 CpGs had a mortality score of 1–10, respectively, and participants without aberrant methylation at any of the 10 CpGs had score of 0. Table 3 shows the associations of score with all-cause mortality. Compared with participants with a score of 0, those who had a score of 1, 2–5 and 5+ had 2-, 3- and 7-fold risk of dying, controlling for all the potential confounding factors. Analyses restricted to only older participants (⩾60 years) yielded essentially the same risk estimates, for example, HRs (95% CI) were 2.14 (1.02–4.47), 3.38 (1.68–6.80) and 7.44 (3.50–15.84), respectively, for a score of 1, 2–5 and 5+ versus score=0. Similar patterns of distribution of deceased were also observed in KORA participants (Supplementary Fig. 4b). Using the cutoff points from the ESTHER cohort defining aberrant methylation of ten CpGs (Supplementary Table 4), replicated analyses in the KORA cohort showed consistent patterns and similar risk estimates (Table 3). Crude HRs (95% CI) for participants with score of 1, 2–5 and 5+ were 1.21 (0.37–3.97), 6.42 (2.55–16.18) and 19.29 (5.58–66.63), respectively, compared with score=0. In the fully adjusted model, three- and six-fold increases in mortality persisted for score levels of 2–5 and 5+, respectively. Using cutoff points (quartiles) of KORA itself defining aberrant methylation of the ten CpGs to build the mortality score, risk estimates were larger than those derived from using ESTHER’s cutoff points. For example, the HR (95% CI) in the fully adjusted model was 7.41 (1.61–34.07) for participants with score of 5+. In addition, a continuous risk score was computed through linear combination of LASSO regression coefficient weighted methylation values of the ten CpGs (the combination formula is presented in Supplementary Fig. 1). A similar trend that mortality monotonously increased with increasing continuous risk score was observed in both the ESTHER (risk score ranged from −3.92 to −0.72; median (IQR), −2.70 (−2.98 to −2.35)) and the KORA cohorts (risk score ranged from −4.40 to −1.51; median (IQR), −3.15 (−3.41 to −2.86)). Figure 2 shows the corresponding dose–response relationships derived from restricted cubic spline regression with adjustment for all the covariates again28.
Sex-specific analyses indicated the associations with all-cause mortality to be stronger among women than among men in both cohorts (Supplementary Table 5). Table 4 shows that the associations of score with CVD mortality were stronger than with cancer mortality in both cohorts. The corresponding survival curves in the ESTHER cohort are presented in Fig. 3. Similar survival curves were also obtained in the KORA cohort (Supplementary Fig. 5).
Table 5 presents the associations of score with all-cause and cause-specific mortality in the ESTHER cohort under consideration of the epigenetic age acceleration (determined by the algorithm of Hannum et al.11). The risk estimates of score for all three mortality outcomes were only very slightly attenuated by adjustment for the epigenetic age acceleration. On the contrary, HRs (95% CI) per 5 years of age acceleration dropped from 1.27 (1.10–1.46), 1.25 (0.98–1.59) and 1.34 (1.05–1.71), respectively, for all-cause, cancer and CVD mortality in the age- and sex-adjusted model to 1.08 (0.92–1.27), 1.15 (0.88–1.51) and 1.12 (0.85–1.48) in the full model. Similar results for the epigenetic age acceleration determined by the algorithm of Horvath et al.12 are presented in Supplementary Table 6.
Discussion
In this EWAS and subsequent validation based on approximately 1,900 older adults with up to 14 years of follow-up, we identified blood DNAm of 58 CpGs across 19 chromosomes to be associated with all-cause mortality. Although there is evidence that genes containing the identified CpGs are related to various types of common diseases, our study was the first to link DNAm of the vast majority of these genes to mortality in the general population. We additionally demonstrated that a risk score based on DNAm of ten identified CpGs was a very strong predictor for all-cause, CVD and cancer mortality, and we confirmed this finding in an independent cohort study. None of the newly identified CpGs overlapped with previously established ageing-related CpGs and the strong associations of score with mortality were also independent from the epigenetic clock.
Of the 58 identified CpGs, the top 1 locus showing the most significant association with mortality was cg05575921 in AHRR, followed by cg21161138 in AHRR, cg26963277 in KCNQ1OT1, cg19859270 in GPR15, cg03636183 in F2RL3, cg19572487 in RARA and cg06126421 in 6p21.33. All these CpGs (except cg26963277 in KCNQ1OT1) were also the top signals in previous EWASs on smoking2. In addition to the 22 CpGs identified to be associated with smoking in previous EWASs2,24, another 26 of the 58 CpGs were also smoking-associated in the current study. Furthermore, even though a few other CpGs were found to be associated with alcohol consumption, diabetes or cancer, such as cg18181703 in SOCS3 and cg26470501 in BCL3, most of them also showed associations with smoking exposure in our analyses (Supplementary Fig. 3). These findings suggest that tobacco smoking is the strongest factor leaving imprints on DNAm such that smoking rather than other common health risk factors accounts for the major burden of morbidity and mortality involving epigenetic programming. Regardless of the underlying mechanisms which remain to be elucidated in further research, it appears worthwhile pointing out that prevention of or intervention on smoking-related DNAm changes may provide major improvement in premature death prevention, given the reversibility of smoking-induced methylomic aberrations29,30.
The current study highlighted several genes or genetic regions as attractive targets for further investigation. The chromosome region 19q13.3 harbours six mortality-related CpGs mapped to BCL3, ERCC1, SLC1A5 and PPP1R15A. Although ERCC1 methylation has been previously reported in lung and breast cancer31,32, DNAm of BCL3, SLC1A5 and PPP1R15A were first linked to health-related outcomes in our study. In light of the known gene functions of BCL3 (pathogenesis of CVD and solid tumours)33,34,35, SLC1A5 (a glutamine transporter in various types of cancer development, progression and response to therapy)36 and PPP1R15A (neurological and CVD pathophysiology, as well as obesity and insulin resistance in animal models)37,38,39, it appears plausible that DNAm may play regulating roles in the development or progression of the respective diseases, which requires elucidation in future studies. This also applies to most of the other genes known to be related to specific diseases whose relationship to methylation-relevant outcomes were first disclosed in our study, such as DNAm of SARS, VCAM1, KCNQ1OT1, MIR19A, SEMA7A, BCL3, PPP1R15A and PDE9A for CVD, DNAm of SQLE, MIR19A, MIR10A, SOCS3, CALR, BCL3 and SLC1A5 for lung cancer, and DNAm of ATL3, SHANK2 and PPP1R15A for neurological diseases. In addition, it is known that the chromosomal region 11p15.5 contains clusters of epigenetically regulated genes, for example, KCNQ1 and KCNQ1OT1, which have been implicated in T2D26,27. We found two mortality-related CpGs (cg07123182 and cg26963277) in KCNQ1OT1 in 11p15.5. Of note, associations with prevalent T2D and mortality were also observed in the current study for cg23190089 in SLC22A18AS (on 11p15.4), a locus located ∼198 kb downstream of cg26963277 in KCNQ1OT1. The chromosome region 11p15.5/4, along with SLC1A5, SQLE and SOCS3 methylation that were suggested to be involved in T2D in both the current and previous studies5,25,40,41,42, therefore appear to be attractive targets for diabetes investigation, and even for CVD given the biological functions of these genes and their methylation in CVD40,43 and also the well-known causal relationship between diabetes and CVD. Similar to SQLE, KCNQ1OT1 and SOCS3, which are involved in diabetes, CVD and various cancers5,26,40,43,44,45,46, most identified genes are characterized by their relevance to multiple diseases, making them the most robust signals on an epigenome-wide scale, which may explain the extremely strong association of the risk score based on identified DNAm markers with all-cause mortality.
Compared with genetic variants related to longevity identified by GWAS, which typically show very small effect sizes of single SNPs, in particular in general population samples47,48, the effect size of even single CpGs identified in the current EWAS were substantial, with HRs ⩾1.17 or ⩽0.83 per s.d. increase of methylation, resulting in the strong overall prediction when combining these CpGs in a risk score. To our knowledge, no comparably strong prediction of mortality based on genetic data has been identified, suggesting that epigenetic data might be more informative for mortality prediction than genetic data.
The recently established epigenetic clock (DNAm age) has received growing attention as an increasing number of studies have uncovered it to be a proxy of biological ageing11,12,13,14,15 and thus potentially providing a measure for assessing health and mortality. Intriguingly, we targeted mortality-related DNAm changes and did not find any overlap with previously established CpGs that are used to determine the DNAm age11,12. Our findings are in line with evidence, suggesting that DNAm involved in ageing or health-related outcomes are mostly regulated by DNAm regions other than the established age-related DNAm19,20,21. The difference could also plausibly result from the fact that DNAm age was originally trained as precisely as possible to track chronological age and might thus be more indicative of natural ageing beyond the effect of disease, as exemplified by the much stronger association of DNAm age with mortality in oldest population (mean age 86.1 years)15 to whom common chronic diseases, such as CVD and cancer, might not continue to pose predominant risks49. Given characteristics of the identified genes and of our study population (mean age 62 years) that is at high-risk age for suffering from major diseases, the currently identified DNAm regions might be more indicative of disease-related outcomes and mortality. Only one previous study has also determined genome-wide methylomic mortality predictors, which were also distinct from the established ageing-associated sites11,12 but also different from signatures discovered in our study. A plausible explanation is that this study was conducted in a very old population (mean age 90+ years), in which causes of death might be distinct from those observed in our study49.
Lack of gene expression data hindered exploration of the roles of the identified DNAm sites in regulating the relevant gene expression. Diseases associated with the identified genes were determined based on a literature search. Whether and how DNAm of those genes are involved in development or progression of the described diseases needs to be elucidated by future multidisciplinary research. For example, genetic factors might potentially be involved in the observed methylation-related mortality and the interplay between genetic factors and these methylation markers warrants to be explored. In the analysis, we did not exclude probes that might be affected by known SNPs as annotated by ‘Infinium HD Methylation SNP List’ (http://support.illumina.com/array/array_kits/infinium_humanmethylation450_beadchip_kit/downloads.html). We later retrieved data of 32 such SNPs for 24 identified CpGs in 581 ESTHER participants of the validation set. Only one SNP-CpG pair (rs524-cg03707168) showed a significant association. However, no association was observed between rs524 and all-cause mortality irrespective of controlling for DNAm of cg03707168, whereas the strong association of cg03707168 with mortality did not change when controlling for rs524. In addition, no interaction was detected between rs524 and cg03707168 in relation to mortality. Nevertheless, potential genetic variants, that is, methylation quantitative trait loci for the identified candidates, should be systematically assessed in further studies. Despite the overall large size of the study population, sample size limitations restricted the list of identified sites, which should be extended in future larger longitudinal studies. In addition, the effect sizes (that is, average methylation difference between survivors and deaths) of most identified loci are relatively small as illustrated in Fig. 1. Plausible reasons are that methylation levels were measured on average 8.2 years before dying and presumably stronger methylation difference restricted to specific causes of death are expected to be diluted in an analysis of all-cause mortality. Another limitation is that DNAm was quantified in whole blood samples. Even though we controlled for the effect of potential cell shift by adjustment for leukocyte composition estimated according to an established and commonly applied algorithm50, residual confounding by leukocyte distribution cannot be ruled out. However, this would not diminish the value of the identified markers for mortality prediction, for which easy accessibility of blood samples is a major advantage. Finally, although we included a variety of covariates in the regression analyses, we cannot exclude the possibility that the observed associations between the identified methylation markers and mortality might be explained to some extent by incompletely controlled or uncontrolled confounding factors. For example, for smoking-related candidates, the observed associations might be partially confounded by imperfect controlling for smoking exposure or by potential confounders related to smoking. Despite its limitations, the prospective nature of the present study, the inclusion of large representative samples of participants from the general population, the long-term follow-up, the hypothesis-free approach with independent internal and external validation, as well as comprehensive adjustment for a variety of common risk factors in data analyses, are major strengths of the current study, which renders novel findings for future verification.
Our previous work using candidate gene approaches has demonstrated the potential use of mortality-related DNAm markers, such as F2RL3 and AHRR, for lung cancer and CVD risk prediction8,9. The clinical implications of other CpGs emerging from the present study for diagnosis, prognosis or even treatment of common diseases, in particular diabetes, CVD and cancer, warrant exploration by future studies. The methylation-based mortality risk score might be a useful tool for population stratification in disease screening and intervention, and its predictive value for ageing-related outcomes, such as frailty and dementia, is worthwhile investigating in future research.
Methods
Study population and data collection
The EWAS and subsequent validation were conducted in the ESTHER study, an ongoing population-based cohort study conducted in Saarland, Germany. The ESTHER cohort, as previously described in detail51, enroled 9,949 older adults (age 50–75 years) by their general practitioners during routine health check-ups between 2000 and 2002. The participants completed a standardized self-administered questionnaire and donated biological samples (blood, stool and urine) during baseline enrolment. Comprehensive medical data, such as the results of a physical assessment, medical diagnoses and drug prescriptions were additionally obtained from the general practitioner. Deaths during follow-up were identified through record linkage with population registries in Saarland. Information on the major cause of death was obtained from death certificates provided by the local health authorities and coded with ICD-10 codes. Deaths from CVD and malignant invasive cancers, respectively, were defined by ICD-10 codes I00-I99 and C00-C97 (excluding non-melanoma skin cancer (C44)).
Genome-wide DNAm measurements were performed in the baseline blood samples of two subsets of the ESTHER participants. Subset-I (discovery panel) consists of participants from a case–cohort study nested within 2,499 ESTHER participants who were consecutively recruited between October 2000 and March 2001, and had sufficient DNA available. Of the 2,499 participants, 406 participants who died during follow-up by March 2013 were the cases in the case–cohort design and 548 participants were randomly selected as the subcohort irrespective of death status during follow-up. The sampling fraction was thus 548/2,499=22%. Subset-II (validation panel) consists of 1,000 ESTHER participants who were recruited between July and October 2000, and who were non-overlapping with the case–cohort samples, among whom 231 deaths were ascertained during follow-up.
Replication in an independent cohort was performed in the KORA F4 study, a population-based cohort consisting of 3,080 participants (age 32–81 years) recruited between 2006 and 2008 from the region of Augsburg, Southern Germany52,53. The vital status of KORA participants was ascertained through population registries inside and outside the study area in December 2011. Causes of death were determined according to death certificates from the Regional Health Department and coded with ICD-9. A random baseline sample consisting of 1,727 participants were selected for methylation analysis, among whom 61 participants died.
All ESTHER and KORA F4 participants provided written informed consent. The ESTHER study was approved by the ethics committees of the University of Heidelberg and of the state medical board of Saarland, Germany. The KORA F4 study was approved by the Ethics Committee of the Bavarian Medical Association.
Methylation assessment
DNAm in whole blood was quantified using the Infinium HumanMethylation450K BeadChip (Illumina, Inc, San Diego, CA, USA) in both ESTHER and KORA F4. Details of methylation analysis in the ESTHER study have been reported previously8,54. According to the manufacturer’s protocol, data were normalized to internal controls provided by Illumina (Illumina normalization). In data pre-processing, probes with detection P-value>0.01, with missing values>10%, probes targeting the sex chromosomes, cross-reactive probes and polymorphic CpGs55 were excluded, leaving 430,363 CpGs for genome-wide screening. In the KORA study, data were pre-processed following the pipeline of Lehne et al.56, probes with detection P-value (1−P-value computed from the background model characterizing the probability that the target sequence signal was distinguishable from the negative controls) >0.01 and missing values >5% were removed, and quantile normalization was applied following stratification of the probe categories into six types, based on probe type and colour channel, using the R package limma57. Leukocyte composition was estimated using the algorithms of Houseman et al.50 in both studies.
Statistical analysis
Discovery and validation of mortality-related CpGs. The ESTHER study populations were described separately in the discovery and validation panel with respect to major sociodemographic characteristics, lifestyle factors and prevalent diseases at baseline. An epigenome-wide screening for mortality-related CpGs was first carried out in the case–cohort samples, using weighted Cox regression models that account for the case–cohort sampling design by Barlow weighting (the inverse of the subcohort sampling fraction, 1/(548/2499))58,59. The models with methylation β-values as explanatory variables were adjusted for age, sex and batch effects. After correcting for multiple testing using the Benjamini–Hochberg approach, CpGs that reached genome-wide significance (FDR <0.05) were entered into the validation phase, in which the associations with mortality were further analysed by multiple Cox regression adjusted for age, sex, batch effects, leukocyte composition50, smoking status (never, former and current smoker), body mass index (kg m−2), physical activity (inactive, low, medium/high), alcohol consumption (grams per day), systolic blood pressure (mmHg), total cholesterol level (mg dL−1), and prevalence of hypertension, CVD, diabetes and cancer. CpGs with FDR<0.05 in the validation panel were deemed as mortality-related loci. A flowchart of study design and data analysis is shown in Supplementary Fig. 1.
Associations of risk factors with mortality-related CpGs. To explore risk factors related to methylation associated with fatal endpoints, sociodemographic characteristics, lifestyle factors and prevalent diseases at baseline were assessed in relation to the methylation levels of the identified CpGs using mixed linear regression models in the validation panel, with batch as random effect, methylation β-value as the dependent variable and independent variables including age, sex, smoking status (never, former and current smoker), body mass index (underweight/normal weight, overweight and obesity), physical activity (inactive, low and medium/high), alcohol consumption (grams per day) and prevalent hypertension, diabetes, CVD and cancer, again controlling for leukocyte composition50. Multiple testing was again corrected for by the Benjamini–Hochberg approach (FDR<0.05).
Mortality risk score. To develop a DNAm-based mortality risk score, we applied the LASSO Cox regression60 with regularization parameter chosen by tenfold cross-validation following the ‘one standard error’ rule61,62, selecting candidates among the identified CpGs. The associations of the score with all-cause, CVD and cancer mortality were assessed first in the validation subset of the ESTHER cohort and then in the independent KORA cohort using multiple Cox regression models, adjusted for the covariates listed above (Supplementary Fig. 1). All analyses were then repeated in men and women separately. In addition, to compare the predictive value of score with that of recently established methylomic predictors of ‘epigenetic age acceleration’ (that is, Δage=DNAm age−chronological age), we assessed the associations of both score and Δage with all-cause mortality simultaneously. DNAm age was calculated according to two commonly applied algorithms introduced by Hannum et al.11 and Horvath et al.12.
The proportional hazards assumption was assessed by martingale-based residuals63. No violations were detected. The LASSO regression analyses were conducted using the R-package ‘glmnet’61. All other statistical analyses in the ESTHER study were carried out in SAS 9.4 (SAS Institute, Cary, NC) and the analyses in the KORA study were conducted in R (version 3.2.3).
Code availability
SAS codes for statistical analysis are available upon request.
Data availability
The data that support the findings of this study are available on reasonable request from the corresponding author (Y.Z.). The data are not publicly available due to restrictions of informed consent.
Additional information
How to cite this article: Zhang, Y. et al. DNA methylation signatures in peripheral blood strongly predict all-cause mortality. Nat. Commun. 8, 14617 doi: 10.1038/ncomms14617 (2017).
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
Dick, K. J. et al. DNA methylation and body-mass index: a genome-wide analysis. Lancet 383, 1990–1998 (2014).
Gao, X., Jia, M., Zhang, Y., Breitling, L. P. & Brenner, H. DNA methylation changes of whole blood cells in response to active smoking exposure in adults: a systematic review of DNA methylation studies. Clin. Epigenetics 7, 113 (2015).
Feinberg, A. P. Genome-scale approaches to the epigenetics of common human disease. Virchows Arch. 456, 13–21 (2010).
Zhong, J., Agha, G. & Baccarelli, A. A. The role of DNA methylation in cardiovascular risk and disease: methodological aspects, study design, and data analysis for epidemiological studies. Circ. Res. 118, 119–131 (2016).
Chambers, J. C. et al. Epigenome-wide association of DNA methylation markers in peripheral blood from Indian Asians and Europeans with incident type 2 diabetes: a nested case-control study. Lancet Diabetes Endocrinol. 3, 526–534 (2015).
Klengel, T., Pape, J., Binder, E. B. & Mehta, D. The role of DNA methylation in stress-related psychiatric disorders. Neuropharmacology 80, 115–132 (2014).
Dang, M. N., Buzzetti, R. & Pozzilli, P. Epigenetics in autoimmune diseases with focus on type 1 diabetes. Diabetes Metab. Res. Rev. 29, 8–18 (2013).
Zhang, Y. et al. Smoking-associated DNA methylation biomarkers and their predictive value for all-cause and cardiovascular mortality. Environ. Health Perspect. 124, 67–74 (2016).
Zhang, Y. et al. F2RL3 methylation, lung cancer incidence and mortality. Int. J. Cancer. 137, 1739–1748 (2015).
Zhang, Y. et al. F2RL3 methylation in blood DNA is a strong predictor of mortality. Int. J. Epidemiol. 43, 1215–1225 (2014).
Hannum, G. et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol. Cell 49, 359–367 (2013).
Horvath, S. DNA methylation age of human tissues and cell types. Genome Biol. 14, R115 (2013).
Marioni, R. E. et al. DNA methylation age of blood predicts all-cause mortality in later life. Genome Biol. 16, 25 (2015).
Breitling, L. P. et al. Frailty is associated with the epigenetic clock but not with telomere length in a German cohort. Clin. Epigenetics 8, 21 (2016).
Christiansen, L. et al. DNA methylation age is associated with mortality in a longitudinal Danish twin study. Aging Cell 15, 149–154 (2016).
Marioni, R. E. et al. The epigenetic clock and telomere length are independently associated with chronological age and mortality. Int. J. Epidemiol. 45, 424–432 (2016).
Lin, Q. et al. DNA methylation levels at individual age-associated CpG sites can be indicative for life expectancy. Aging (Albany NY) 8, 394–401 (2016).
Moore, A. Z. et al. Change in epigenome-wide DNA methylation over 9 years and subsequent mortality: results from the InCHIANTI Study. J. Gerontol. A Biol. Sci. Med. Sci. 71, 1029–1035 (2016).
Bell, J. T. et al. Epigenome-wide scans identify differentially methylated regions for age and age-related phenotypes in a healthy ageing population. PLoS Genet. 8, e1002629 (2012).
Marttila, S. et al. Ageing-associated changes in the human DNA methylome: genomic locations and effects on gene expression. BMC Genomics 16, 179 (2015).
Jung, M. & Pfeifer, G. P. Aging and DNA methylation. BMC Biol. 13, 7 (2015).
Florath, I., Butterbach, K., Müller, H., Bewerunge-Hudler, M. & Brenner, H. Cross-sectional and longitudinal changes in DNA methylation with age: an epigenome-wide analysis revealing over 60 novel age-associated CpG sites. Hum. Mol. Genet. 23, 1186–1201 (2014).
Weidner, C. I. et al. Aging of blood can be tracked by DNA methylation changes at just three CpG sites. Genome Biol. 15, R24 (2014).
Teschendorff, A. E. et al. Correlation of smoking-associated DNA methylation changes in buccal cells with DNA methylation changes in epithelial cancer. JAMA Oncol. 1, 476–485 (2015).
Al Muftah, W. A. et al. Epigenetic associations of type 2 diabetes and BMI in an Arab population. Clin. Epigenetics 8, 13 (2016).
Travers, M. E. et al. Insights into the molecular mechanism for type 2 diabetes susceptibility at the KCNQ1 locus from temporal changes in imprinting status in human islets. Diabetes 62, 987–992 (2013).
Nilsson, E. et al. Altered DNA methylation and differential expression of genes influencing metabolism and inflammation in adipose tissue from subjects with type 2 diabetes. Diabetes 63, 2962–2976 (2014).
Desquilbet, L. & Mariotti, F. Dose-response analyses using restricted cubic spline functions in public health research. Stat. Med. 29, 1037–1057 (2010).
Zhang, Y., Yang, R., Burwinkel, B., Breitling, L. P. & Brenner, H. F2RL3 methylation as a biomarker of current and lifetime smoking exposures. Environ. Health Perspect. 122, 131–137 (2014).
Guida, F. et al. Dynamics of smoking-induced genome-wide methylation changes with time since smoking cessation. Hum. Mol. Genet. 24, 2349–2359 (2015).
Do, H. et al. A critical re-assessment of DNA repair gene promoter methylation in non-small cell lung carcinoma. Sci. Rep. 4, 4186 (2014).
Singh, K. P., Treas, J., Tyagi, T. & Gao, W. DNA demethylation by 5-aza-2-deoxycytidine treatment abrogates 17 beta-estradiol-induced cell growth and restores expression of DNA repair genes in human breast cancer cells. Cancer Lett. 316, 62–69 (2012).
Marchetti, G. et al. An integrated genomic-transcriptomic approach supports a role for the proto-oncogene BCL3 in atherosclerosis. Thromb. Haemost. 113, 655–663 (2015).
Lu, W. et al. MCPIP1 selectively destabilizes transcripts associated with an anti-apoptotic gene expression program in breast cancer cells that can elicit complete tumor regression. Cancer Res. 76, 1429–1440 (2016).
Ahlqvist, K., Saamarthy, K., Syed Khaja, A. S., Bjartell, A. & Massoumi, R. Expression of Id proteins is regulated by the Bcl-3 proto-oncogene in prostate cancer. Oncogene 32, 1601–1608 (2013).
Bhutia, Y. D., Babu, E., Ramachandran, S. & Ganapathy, V. Amino Acid transporters in cancer and their relevance to “glutamine addiction”: novel targets for the design of a new class of anticancer drugs. Cancer Res. 75, 1782–1788 (2015).
Honjo, Y. et al. Increased GADD34 in oligodendrocytes in Alzheimer's disease. Neurosci. Lett. 602, 50–55 (2015).
Morton, E., Macrae, I. M., McCabe, C., Brown, S. M. & White, F. Identification of the growth arrest and DNA damage protein GADD34 in the normal human heart and demonstration of alterations in expression following myocardial ischaemia. Int. J. Cardiol. 107, 126–129 (2006).
Nishio, N. & Isobe, K. GADD34-deficient mice develop obesity, nonalcoholic fatty liver disease, hepatic carcinoma and insulin resistance. Sci. Rep. 5, 13519 (2015).
Ding, J. et al. Alterations of a cellular cholesterol metabolism network are a molecular feature of obesity-related type 2 diabetes and cardiovascular disease. Diabetes 64, 3464–3474 (2015).
Dayeh, T. et al. Genome-wide DNA methylation analysis of human pancreatic islets from type 2 diabetic and non-diabetic donors identifies candidate genes that influence insulin secretion. PLoS Genet. 10, e1004160 (2014).
Kriebel, J. et al. Association between DNA methylation in whole blood and measures of glucose metabolism: KORA F4 Study. PLoS ONE 11, e0152314 (2016).
Vausort, M., Wagner, D. R. & Devaux, Y. Long noncoding RNAs in patients with acute myocardial infarction. Circ. Res. 115, 668–677 (2014).
Parris, T. Z. et al. Frequent MYC coamplification and DNA hypomethylation of multiple genes on 8q in 8p11-p12-amplified breast carcinomas. Oncogenesis 3, e95 (2014).
Sunamura, N. et al. Regulation of functional KCNQ1OT1 lncRNA by beta-catenin. Sci. Rep. 6, 20690 (2016).
Chen, X. et al. Analysis of DNA methylation and gene expression in radiation-resistant head and neck tumors. Epigenetics 10, 545–561 (2015).
Walter, S. et al. A genome-wide association study of aging. Neurobiol. Aging 32, e28 (2011).
Hallden, S. et al. Smoking and obesity associated BDNF gene variance predicts total and cardiovascular mortality in smokers. Heart 99, 949–953 (2013).
Lee, S. J., Go, A. S., Lindquist, K., Bertenthal, D. & Covinsky, K. E. Chronic conditions and mortality among the oldest old. Am. J. Public Health 98, 1209–1214 (2008).
Houseman, E. A. et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics 13, 86 (2012).
Scho¨ttker, B. et al. Strong associations of 25-hydroxyvitamin D concentrations with all-cause, cardiovascular, cancer, and respiratory disease mortality in a large cohort study. Am. J. Clin. Nutr. 97, 782–793 (2013).
Holle, R., Happich, M., Lowel, H. & Wichmann, H. E. MONICA/KORA Study Group. KORA--a research platform for population based health research. Gesundheitswesen 67, S19–S25 (2005).
Zeilinger, S. et al. Tobacco smoking leads to extensive genome-wide changes in DNA methylation. PLoS ONE 8, e63812 (2013).
Zhang, Y., Florath, I., Saum, K. U. & Brenner, H. Self-reported smoking, serum cotinine, and blood DNA methylation. Environ. Res. 146, 395–403 (2016).
Chen, Y. A. et al. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics 8, 203–209 (2013).
Lehne, B. et al. A coherent approach for analysis of the Illumina HumanMethylation450 BeadChip improves data quality and performance in epigenome-wide association studies. Genome Biol. 16, 37 (2015).
Smyth, G. in Bioinformatics and Computational Biology Solutions Using R and Bioconductor (eds Gentleman, R. et al.) 397–420 (Springer, 2005).
Barlow, W. E. Robust variance estimation for the case-cohort design. Biometrics 50, 1064–1072 (1994).
Kulathinal, S., Karvanen, J., Saarela, O. & Kuulasmaa, K. Case-cohort design in practice - experiences from the MORGAM Project. Epidemiol. Perspect. Innov. 4, 15 (2007).
Tibshirani, R. The lasso method for variable selection in the Cox model. Stat. Med. 16, 385–395 (1997).
Friedman, J., Hastie, T. & Tibshirani, R. Regularization paths for generalized linear models via coordinate descent. J. Stat. Softw. 33, 1–22 (2010).
Hastie, T., Tibshirani, R. & Friedman, J. The Elements of Statistical Learning 2nd edn Springer (2009).
Lin, D., Wei, L. & Ying, Z. Checking the Cox model with cumulative sums of martingale-based residuals. Biometrika 80, 557–572 (1993).
Acknowledgements
The ESTHER study was supported by the Baden-Württemberg State Ministry of Science, Research and Arts (Stuttgart, Germany), the Federal Ministry of Education and Research (Berlin, Germany) and the Federal Ministry of Family Affairs, Senior Citizens, Women and Youth (Berlin, Germany). The KORA study was initiated and financed by the Helmholtz Zentrum München–German Research Center for Environmental Health, which is funded by the German Federal Ministry of Education and Research (BMBF) and by the State of Bavaria. The research leading to these results has received funding from the European Union Seventh Framework Programme (FP7/2007-2013) under grant agreement number 603288 (SysVasc). The sponsors had no role in the study design, in the collection, analysis and interpretation of data and preparation, review or approval of the manuscript.
Author information
Authors and Affiliations
Contributions
Y.Z. conceived the study, carried out the data analyses, interpreted the data and drafted the manuscript. R.W. conducted data analyses in the KORA study. J.H. computed leukocyte composition in the ESTHER study. L.P.B. contributed to the design of the study. K.S. was responsible for coordination of the ESTHER study. B.S. and B.H. were responsible for coordination of follow-up and work-up of follow-up data of the ESTHER study. M.W. and A.P. are the investigators of the KORA study. H.B. conducted the ESTHER study and contributed to all aspects of this work. All authors contributed to revision of the manuscript and approved the final version for submission.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Figures, Supplementary Tables, and Supplementary References (PDF 643 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Zhang, Y., Wilson, R., Heiss, J. et al. DNA methylation signatures in peripheral blood strongly predict all-cause mortality. Nat Commun 8, 14617 (2017). https://doi.org/10.1038/ncomms14617
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms14617