Assessing the causal role of epigenetic clocks in the development of multiple cancers: a Mendelian randomization study
Abstract
Background:
Epigenetic clocks have been associated with cancer risk in several observational studies. Nevertheless, it is unclear whether they play a causal role in cancer risk or if they act as a non-causal biomarker.
Methods:
We conducted a two-sample Mendelian randomization (MR) study to examine the genetically predicted effects of epigenetic age acceleration as measured by HannumAge (nine single-nucleotide polymorphisms (SNPs)), Horvath Intrinsic Age (24 SNPs), PhenoAge (11 SNPs), and GrimAge (4 SNPs) on multiple cancers (i.e. breast, prostate, colorectal, ovarian and lung cancer). We obtained genome-wide association data for biological ageing from a meta-analysis (N = 34,710), and for cancer from the UK Biobank (N cases = 2671–13,879; N controls = 173,493–372,016), FinnGen (N cases = 719–8401; N controls = 74,685–174,006) and several international cancer genetic consortia (N cases = 11,348–122,977; N controls = 15,861–105,974). Main analyses were performed using multiplicative random effects inverse variance weighted (IVW) MR. Individual study estimates were pooled using fixed effect meta-analysis. Sensitivity analyses included MR-Egger, weighted median, weighted mode and Causal Analysis using Summary Effect Estimates (CAUSE) methods, which are robust to some of the assumptions of the IVW approach.
Results:
Meta-analysed IVW MR findings suggested that higher GrimAge acceleration increased the risk of colorectal cancer (OR = 1.12 per year increase in GrimAge acceleration, 95% CI 1.04–1.20, p = 0.002). The direction of the genetically predicted effects was consistent across main and sensitivity MR analyses. Among subtypes, the genetically predicted effect of GrimAge acceleration was greater for colon cancer (IVW OR = 1.15, 95% CI 1.09–1.21, p = 0.006), than rectal cancer (IVW OR = 1.05, 95% CI 0.97–1.13, p = 0.24). Results were less consistent for associations between other epigenetic clocks and cancers.
Conclusions:
GrimAge acceleration may increase the risk of colorectal cancer. Findings for other clocks and cancers were inconsistent. Further work is required to investigate the potential mechanisms underlying the results.
Funding:
FMB was supported by a Wellcome Trust PhD studentship in Molecular, Genetic and Lifecourse Epidemiology (224982/Z/22/Z which is part of grant 218495/Z/19/Z). KKT was supported by a Cancer Research UK (C18281/A29019) programme grant (the Integrative Cancer Epidemiology Programme) and by the Hellenic Republic’s Operational Programme ‘Competitiveness, Entrepreneurship & Innovation’ (OΠΣ 5047228). PH was supported by Cancer Research UK (C18281/A29019). RMM was supported by the NIHR Biomedical Research Centre at University Hospitals Bristol and Weston NHS Foundation Trust and the University of Bristol and by a Cancer Research UK (C18281/A29019) programme grant (the Integrative Cancer Epidemiology Programme). RMM is a National Institute for Health Research Senior Investigator (NIHR202411). The views expressed are those of the author(s) and not necessarily those of the NIHR or the Department of Health and Social Care. GDS and CLR were supported by the Medical Research Council (MC_UU_00011/1 and MC_UU_00011/5, respectively) and by a Cancer Research UK (C18281/A29019) programme grant (the Integrative Cancer Epidemiology Programme). REM was supported by an Alzheimer’s Society project grant (AS-PG-19b-010) and NIH grant (U01 AG-18-018, PI: Steve Horvath). RCR is a de Pass Vice Chancellor’s Research Fellow at the University of Bristol.
Editor's evaluation
This paper is of broad interest to researchers seeking to disentangle the health impact of epigenetic age acceleration, and will provide a substantive empirical contribution to the literature. The authors were very meticulous in addressing all the concerns from the reviewers, which has further improved the paper.
https://doi.org/10.7554/eLife.75374.sa0eLife digest
Have you noticed that some people seem to get older faster than others? Scientists have previously found that a chemical tag on DNA known as DNA methylation can be used to predict an individual’s chronological age. However, age predicted using DNA methylation (also known as biological or epigenetic age) does not always perfectly correspond to chronological age. Indeed, some people’s biological age is higher than their years, while other people’s is lower.
When an individual’s biological age is higher than their chronological age, they are said to be experiencing ‘epigenetic age acceleration’. This type of accelerated ageing, which can be measured with ‘epigenetic clocks’ based on DNA methylation, has been associated with several adverse health outcomes, including cancer. This means that epigenetic clocks may improve our ability to predict cancer risk and detect cancer early. However, it is still unclear whether accelerated biological ageing causes cancer, or whether it simply correlates with the disease.
Morales-Berstein et al. wanted to investigate whether epigenetic age acceleration, as measured by epigenetic clocks, plays a role in the development of several cancers. To do so, they used an approach known as Mendelian randomization. Using genetic variants as natural experiments, they studied the effect of different measures of epigenetic age acceleration on cancer risk.
Their work focused on five types of cancer: breast, colorectal, prostate, ovarian and lung cancer. They used genetic association data from people of European ancestry to determine whether genetic variants that are strongly associated with accelerated ageing are also strongly associated with cancer. The results showed that one of the DNA methylation markers used as an estimate of biological ageing could be directly related to the risk of developing colorectal cancer.
This work provides new insights into the relationship between markers of biological ageing and cancer. Similar relationships should also be studied in other groups of people and for other cancer sites. The results suggest that reversing biological ageing by altering DNA methylation could prevent or delay the development of colorectal cancer.
Introduction
DNA methylation (DNAm) at specific cytosine-phosphate-guanine (CpG) sites has been found to be strongly correlated with chronological age. Biological age, as predicted by DNAm patterns at specific CpG sites, may differ from chronological age on an individual basis. Observational evidence suggests that epigenetic age acceleration (i.e. when an individual’s biological age is greater than their chronological age) may be associated with an increased risk of mortality and age-related diseases, including cancer (Fransquet et al., 2019).
Epigenetic clocks are heritable indicators of biological ageing derived from DNAm data. Each clock is based on DNAm levels measured at a different set of CpG sites, which capture distinctive features of epigenetic ageing (Liu et al., 2020). ‘First-generation’ epigenetic clocks, such as HannumAge (Hannum et al., 2013) and Intrinsic HorvathAge (Horvath, 2013), have been derived from DNAm levels at CpG sites found to be strongly associated with chronological age. HannumAge is trained on 71 age-related CpGs found in blood (Hannum et al., 2013), while Intrinsic HorvathAge is based on 353 age-related CpGs found in several human tissues and cell types, and is further adjusted for blood cell counts (Horvath, 2013). More recently, ‘second-generation’ epigenetic clocks, such as, PhenoAge (Levine et al., 2018) and GrimAge (Lu et al., 2019a), have been developed to predict age-related morbidity and mortality. PhenoAge incorporates data from 513 CpGs associated with mortality and nine clinical biomarkers (i.e. albumin, creatinine, serum glucose, C-reactive protein, lymphocyte percentage, mean corpuscular volume, red cell distribution width, alkaline phosphatase and leukocyte count) (Levine et al., 2018), and GrimAge includes data from 1,030 CpGs associated with smoking pack-years and seven plasma proteins (i.e. cystatin C, leptin, tissue inhibitor metalloproteinases 1, adrenomedullin, beta-2-microglobulin, growth differentiation factor 15, and plasminogen activation inhibitor 1 (PAI-1)) (Lu et al., 2019a). Due to differences in their composition, HannumAge and Intrinsic HorvathAge are better predictors of chronological age (Hannum et al., 2013; Horvath, 2013), while PhenoAge and GrimAge stand out for their ability to predict health and lifespan (Levine et al., 2018; Lu et al., 2019a; McCrory et al., 2021).
Several studies suggest that HannumAge, Intrinsic HorvathAge, PhenoAge and GrimAge acceleration are associated with cancer risk (Levine et al., 2018; Ambatipudi et al., 2017; Levine et al., 2015; Dugue et al., 2021; Kresovich et al., 2019b; Kresovich et al., 2019a; Zheng et al., 2016). In contrast, others indicate that evidence in support of this claim is weak or non existent (Dugué et al., 2018; Hillary et al., 2020; Durso et al., 2017; Wang et al., 2021). This lack of consensus could be explained by biases that often affect observational research, such as reverse causation (e.g. cancer influencing the epigenome and not the other way around) and residual confounding (e.g. unmeasured, or imprecisely measured confounders of the association between epigenetic age acceleration and cancer) (Relton and Davey Smith, 2012).
The strength of the associations between epigenetic age acceleration and different cancers has also been found to vary across epigenetic clocks. For instance, positive associations between epigenetic age acceleration and colorectal cancer seem to be much stronger when biological age is estimated using second-generation clocks (i.e. PhenoAge and GrimAge) (Dugue et al., 2021) rather than first-generation clocks (i.e. HannumAge and Intrinsic HorvathAge) (Dugué et al., 2018; Durso et al., 2017). Lack of consensus across epigenetic clocks could be explained by differences in their algorithms (which may reflect different mechanisms of biological ageing), as well as heterogeneity in study designs (Fransquet et al., 2019). Furthermore, even if there were a consensus, it would still be unclear whether age-related DNA methylation plays a causal role in cancer risk or if it merely acts as a non-causal prognostic biomarker.
Mendelian randomization (MR), a method that uses genetic variants as instrumental variables to infer causality between a modifiable exposure and an outcome, is less likely to be affected by residual confounding and reverse causation than traditional observational methods (Davey Smith and Ebrahim, 2003). A recent genome-wide association study (GWAS) meta-analysis has revealed 137 genetic loci associated with epigenetic age acceleration (as measured by six epigenetic biomarkers) that may be used within an MR framework (McCartney et al., 2021).
McCartney et al., 2021 used IVW MR, MR-Egger, weighted median and weighted mode methods to explore the genetically predicted effects of HannumAge, Intrinsic HorvathAge, PhenoAge and GrimAge acceleration on breast, ovarian, and lung cancer. Here, we extend this analysis to include colorectal and prostate cancer (two of the most common cancers worldwide Sung et al., 2021) and use additional methods and datasets to verify the robustness of our findings.
The aim of this two-sample MR study was to examine the genetically predicted effects of epigenetic age acceleration (as measured by HannumAge Hannum et al., 2013, Horvath Intrinsic Age Horvath, 2013, PhenoAge Levine et al., 2018 and GrimAge Lu et al., 2019a) on multiple cancers (i.e., breast, prostate, colorectal, ovarian and lung cancer) using summary genetic association data from (1) McCartney et al. (N = 34,710) (McCartney et al., 2021), (2) the UK Biobank (N cases = 2671–13,879; N controls = 173,493–372,016), (3) FinnGen (N cases = 719–8401; N controls = 74,685–174,006) and (4) several international cancer genetic consortia (N cases = 11,348–122,977; N controls = 15,861–105,974).
Materials and methods
Reporting guidelines
Request a detailed protocolThis study has been reported according to the STROBE-MR guidelines (Skrivankova et al., 2021; Supplementary file 2).
Genetic instruments for epigenetic age acceleration
Request a detailed protocolWe obtained summary genetic association estimates for epigenetic age acceleration measures of HannumAge (Hannum et al., 2013), Intrinsic HorvathAge (Horvath, 2013), PhenoAge (Levine et al., 2018), and GrimAge (Lu et al., 2019a) from a recent GWAS meta-analysis of biological ageing (McCartney et al., 2021), which included 34,710 participants of European ancestry. Across the 28 European ancestry studies considered in the analysis, 57.3% of participants were female. A detailed description of the methods that were used can be found in the publication by McCartney et al., 2021. In short, the Horvath epigenetic age calculator software (https://dnamage.genetics.ucla.edu) or standalone scripts were used to calculate age adjusted DNAm estimates. Outlier samples with clock methylation estimates of +/−5 s.d. from the mean were excluded from further analysis. SNPs were genotyped and imputed independently for each cohort included in the meta-analysis. Genotypes were imputed using either the HRC or the 1000 Genomes Project Phase 3 reference panels in all cohorts but the Sister Study (which did not have imputed data at the time of analysis) and the Genetics of Lipid Lowering Drugs and Diet Network Study (which used whole-genome sequencing data). GWAS summary statistics were obtained in each cohort using additive linear models adjusted for sex and genetic principal components, and they were later processed and harmonised using the ‘EasyQC’ R package. Fixed effect meta-analyses were performed using the METAL software (Willer et al., 2010).
We used the clump_data function in the ‘TwoSampleMR’ R package to select GWAS-significant SNPs (P < 5 × 10−8) for each epigenetic age acceleration measure and perform linkage disequilibrium (LD) clumping (r2 <0.001) using the European reference panel from the 1000 Genomes Project Phase 3 v5.
We identified 9 independent SNPs for HannumAge, 24 for Intrinsic HorvathAge, 11 for PhenoAge and 4 for GrimAge (Supplementary file 1 — Table 1). The proportions of trait variance explained by genetic instruments (R2) and instrument strength (F-statistic) were calculated using the following formulae: R2 = (2β2×MAF×(1-MAF))/(2β2×MAF×(1-MAF) + 2 N × MAF × (1-MAF)×SE2) and F = (R2×(N-2))/(1-R2) (where MAF = effect allele frequency, β = effect estimate of the SNP in the exposure GWAS, SE = standard error, N = sample size) (Palmer et al., 2012). The genetic instruments for HannumAge, Intrinsic HorvathAge, PhenoAge and GrimAge acceleration explained 1.48%, 4.41%, 1.86%, and 0.47% of the trait variance, respectively. All the selected SNPs had F-statistics greater than 10 (HannumAge median 38 and range 31–99, Intrinsic HorvathAge median 47 and range 31–240, PhenoAge median 45 and range 32–89, GrimAge median 36 and range 31–45).
Genetic Association Data sources for cancer outcomes
Request a detailed protocolWe obtained summary-level genetic association data for cancer outcomes from the UK Biobank, FinnGen and several international cancer genetic consortia: the Breast Cancer Association Consortium (BCAC), the Ovarian Cancer Association Consortium (OCAC), the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA), the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (PRACTICAL), the International Lung Cancer Consortium (ILCCO) and the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO) (Table 1). Further details of the studies and the data obtained are described in Appendix 1.
We extracted genetic association data for the selected SNPs from each cancer GWAS (for breast, prostate, colorectal, ovarian and lung cancers). LD proxies (r2 >0.8) were used when the SNPs of interest were missing from the cancer GWAS dataset. The proxies were located using the MR-Base platform, which calculates LD using the European subset of individuals from the 1000 Genomes Project reference panel as above (Hemani et al., 2018). The ‘LDlinkR’ R package version 1.1.2 was used to find proxies for cancer data that were not included in the MR-Base platform. The exposure and outcome datasets were then harmonised to ensure the genetic associations reflect the same effect allele. Palindromic SNPs with minor allele frequencies (MAF) <0.3 were aligned, while those with MAF ≥0.3 or mismatching strands were excluded.
Power calculations
Request a detailed protocolStatistical power was calculated using an online calculator for MR available at: https://shiny.cnsgenomics.com/mRnd/. Calculations were performed separately for each clock-cancer combination. They were based on a type one error rate of 0.05, the proportion of phenotypic variance explained by genetic variants (R2) for each measure of epigenetic age acceleration, and the total number of cases and controls included in the meta-analysis for each cancer. Across combinations of the four epigenetic clock acceleration and five cancer measures, we had 80% power to detect ORs as small as 1.04–1.39 (Supplementary file 1 — Table s2).
Statistical analysis
Request a detailed protocolWe estimated the genetically predicted effects of epigenetic age acceleration (as measured by HannumAge Hannum et al., 2013, Horvath Intrinsic Age Horvath, 2013, PhenoAge Levine et al., 2018 and GrimAge Lu et al., 2019a) on multiple cancers (i.e. breast, prostate, colorectal, ovarian, and lung cancer) using a two-sample MR framework (Figure 1).
Main analyses
Request a detailed protocolMain analyses were performed using multiplicative random effects inverse variance weighted (IVW) MR, a method that combines the genetically predicted effect of epigenetic age acceleration on cancer across genetic variants (Burgess et al., 2013). This is the default IVW MR method in the ‘TwoSampleMR’ R package, as it accounts for excess heterogeneity across SNP-specific estimates (as opposed to the fixed effect IVW method) and it does not affect the relative weights of individual SNP estimates (in contrast to the additive random effects IVW method) (Bowden et al., 2017).
We used fixed effect meta-analysis to pool results across studies (i.e. UK Biobank, FinnGen and international consortia). For colorectal cancer, we only pooled FinnGen and GECCO estimates, since UK Biobank participants were already included in GECCO. I2 statistics and their corresponding confidence intervals were used to estimate heterogeneity across study estimates (von Hippel, 2015). A Benjamini-Hochberg false discovery rate (FDR) < 5% was used to correct the pooled main IVW results for multiple testing (Benjamini and Hochberg, 1995). This correction was applied considering a total of 20 independent statistical tests (4 clocks x 5 cancers = 20).
Sensitivity analyses
Request a detailed protocolMR assumes genetic instruments for epigenetic age acceleration are (1) associated with epigenetic age acceleration (relevance assumption), (2) independent of confounders of the association between the instruments and cancer (independence assumption), and (3) only associated with cancer through their effect on epigenetic age acceleration (exclusion restriction assumption) (Didelez and Sheehan, 2007; Davies et al., 2018).
As a sensitivity analysis and to test for potential violations of the relevance assumption, we calculated F-statistics and the R2 for each measure of epigenetic age acceleration (Burgess and Thompson, 2011). Other sensitivity analyses included MR-Egger (Bowden et al., 2015), weighted median (Bowden et al., 2016) and weighted mode (Hartwig et al., 2017) methods, which are robust to some of the assumptions of the IVW approach (described in Appendix 1). These results were also pooled across studies, as explained above. Consistency across different MR methods would suggest that it is less likely that the independence and exclusion restriction assumptions are violated.
We further assessed the validity of the independence assumption by conducting MR analyses using negative control outcomes (i.e. skin colour, ease of skin tanning). Evidence of causality between our genetic instruments for epigenetic age acceleration and these negative control outcomes would suggest potential bias due to population stratification that has not been fully accounted for through adjustments in the GWAS (Sanderson et al., 2021). We also assessed the genetically predicted effect of epigenetic age acceleration on cancer risk factors (i.e. body mass index, waist circumference, pack years of smoking, time spent doing vigorous physical activity, age completed full time education, years of schooling, and alcohol intake frequency) to detect potential violations of the exclusion restriction assumption. GWAS data for negative control outcomes and cancer risk factors were obtained using the University of Bristol’s IEU OpenGWAS API (for more details, see Appendix 1).
Where associations between genetically predicted epigenetic age acceleration and cancer were identified, we additionally performed single-SNP two-sample MR analysis to assess whether the effects were likely to be driven by a single SNP. We used the METAL software (Willer et al., 2010) to conduct a GWAS meta-analysis of cancer genetic association data obtained from the UK Biobank, FinnGen and international cancer genetic consortia. We then used these meta-analysed summary statistics in two-sample MR analyses. Scatter plots showing the effects of genetic instruments on epigenetic clock acceleration against their effects on cancer were created using the ‘TwoSampleMR’ R package. Additionally, Cochran’s Q statistics were used to quantify global heterogeneity across SNP-specific MR estimates (Bowden et al., 2019) and MR-Egger intercept tests were performed to detect horizontal pleiotropy (Bowden et al., 2015).
We also used Causal Analysis using Summary Effect Estimates (CAUSE) (Morrison et al., 2020), a method that uses genome-wide summary statistics to disentangle causality (i.e. SNPs are associated with cancer through their effect on epigenetic age acceleration) from correlated horizontal pleiotropy (i.e. SNPs are associated with epigenetic age acceleration and cancer through a shared heritable factor), while taking into account uncorrelated horizontal pleiotropy (i.e. SNPs are associated with epigenetic age acceleration through separate mechanisms). It uses Bayesian modelling to assess whether the sharing model (i.e. model that fixes the causal effect at zero) fits the data at least as well as the causal model (i.e. model that allows a causal effect different from zero).
Secondary analyses
Request a detailed protocolAs a secondary analysis, we conducted two-sample MR of epigenetic age acceleration and cancer subtypes (i.e. breast cancer: ER+, ER-, triple negative, luminal B/HER2-negative-like, HER2-enriched-like, luminal A-like, luminal B-like, BRCA1 and BRCA2; ovarian cancer: high-grade serous, low-grade serous, invasive mucinous, clear cell, endometrioid, BRCA1 and BRCA2; prostate cancer: advanced, advanced [vs non-advanced], early onset, high risk [vs low risk], and high risk [vs low and intermediate risk]; lung cancer: adenocarcinoma and squamous cell; colorectal cancer: colon-specific, proximal colon-specific, distal colon-specific, rectal-specific, male and female) (Appendix 1).
We also performed two-sample MR analyses of epigenetic age acceleration and parental history of cancer in the UK Biobank for breast, prostate, lung and bowel cancer (Appendix 1). Data on parental history of ovarian cancer were not available in UK Biobank. Family history data correlate with combined hospital record and questionnaire data and it has been suggested that they provide better power to detect GWAS-significant associations for some phenotypes in the UK Biobank (DeBoever et al., 2020). Therefore, we expected these results to be consistent with those obtained in the main analyses.
MR results were reported as the odds ratio (OR) of site-specific cancer per one year increase in genetically predicted epigenetic age acceleration. These did not require any scale transformations, as the GWAS of biological ageing (McCartney et al., 2021) reported epigenetic age acceleration in years.
LD Score regression (Bulik-Sullivan et al., 2015) was used to identify genome-wide genetic correlations between epigenetic age acceleration and cancer. Genetic correlations were estimated using full GWAS summary statistics for the epigenetic clocks and cancer, as well as the 1000 Genomes Project European LD reference panel. Traits with mean heritability chi-square values < 1.02 were excluded from the analyses.
Finally, bidirectional MR analyses were conducted to assess the causality and directionality of the link between epigenetic clock acceleration and telomere length, another measure of biological ageing that has been shown to influence cancer risk in prior MR studies (Telomeres Mendelian Randomization Collaboration et al., 2017; Gao et al., 2020; Kuo et al., 2019). The MR Steiger test of directionality was used to confirm the assumption that the exposure causes the outcome is valid (Hemani et al., 2017). We also corroborated our findings by rerunning the analyses using data that had undergone Steiger filtering to remove SNPs that explained more variance in the outcome than in the risk factor. Genetic association data for measured telomere length were obtained from Codd et al., 2021, the largest GWAS of telomere length available through the OpenGWAS API at the time of analysis (N = 472,174, for more details, see Appendix 1).
All MR analyses were performed using R software version 4.0.2. Two sample MR analyses were conducted using the ‘TwoSampleMR’ package version 0.5.5. Meta-analyses of IVW results were performed using the ‘meta’ package version 4.18. GWAS meta-analyses used to perform single-SNP MR analyses were done using the METAL software (Willer et al., 2010). CAUSE analyses were conducted using the ‘cause’ package version 1.2.0. Forest plots were created using the ‘ggforestplot’ package version 0.1.0. LD Scores were computed using the ‘ldsc’ command line tool version 1.0.1. The code used in this study is available at: https://github.com/fernandam93/epiclocks_cancer.
Results
Breast cancer
We did not find strong evidence of causality between epigenetic age acceleration and breast cancer (GrimAge IVW OR = 0.98, 95% CI 0.95–1.00, p = 0.08; PhenoAge IVW OR = 0.99, 95% CI 0.98–1.01, p = 0.23; HannumAge IVW OR = 0.99, 95% CI 0.97–1.02, p = 0.63; and Intrinsic HorvathAge IVW OR = 0.99, 95% CI 0.98–1.00, p = 0.13) (Figure 2, Appendix 2—figure 1, Appendix 2—figure 2, Appendix 2—figure 3, Appendix 2—figure 4, Appendix 2—figure 5, Appendix 2—figure 6, Appendix 2—figure 7, Appendix 2—figure 8, Supplementary file 1 — Table s3, Supplementary file 1 — Table s4, Supplementary file 1 — Table s5, Supplementary file 1 — Table s6).
Ovarian cancer
There was also limited evidence of causality between epigenetic age acceleration and ovarian cancer (GrimAge IVW OR = 0.99, 95% CI 0.93–1.06, p = 0.78; PhenoAge IVW OR = 0.98, 95% CI 0.96–1.01, p = 0.24; HannumAge IVW OR = 1.00, 95% CI 0.96–1.04, p = 0.95; and Intrinsic HorvathAge IVW OR = 1.00, 95% CI 0.97–1.02, p = 0.89) (Figure 2, Appendix 2—figure 1, Appendix 2—figure 2, Appendix 2—figure 3, Appendix 2—figure 4, Appendix 2—figure 5, Appendix 2—figure 6, Appendix 2—figure 7, Appendix 2—figure 8, Supplementary file 1 — Table s3, Supplementary file 1 — Table s4, Supplementary file 1 — Table s5, Supplementary file 1 — Table s6).
Prostate cancer
Meta-analysed IVW MR findings suggested that genetically predicted GrimAge acceleration decreased the risk of prostate cancer (OR = 0.93 per year increase in GrimAge acceleration, 95% CI 0.87–0.99, p = 0.02) (Figure 2, Supplementary file 1 — Table s3, Supplementary file 1 — Table s4, Supplementary file 1 — Table s5, Supplementary file 1 — Table s6). Although the direction of the genetically predicted effect was consistent across main and sensitivity MR analyses (i.e. MR-Egger, weighted median and weighted mode) (Appendix 2—figure 1, Supplementary file 1 — Table s3, Supplementary file 1 — Table s4, Supplementary file 1 — Table s5, Supplementary file 1 — Table s6), the main IVW result for GrimAge and prostate cancer did not withstand multiple testing correction (FDR p = 0.16) (Supplementary file 1 — Table s6).
We did not find consistent evidence of causality between other measures of epigenetic age acceleration and prostate cancer (PhenoAge IVW OR = 1.01, 95% CI 0.99–1.03, p = 0.31; HannumAge IVW OR = 0.98, 95% CI 0.95–1.02, p = 0.39; and Intrinsic HorvathAge IVW OR = 0.99, 95% CI 0.98–1.01, p = 0.42) (Figure 2, Appendix 2—figure 1, Appendix 2—figure 2, Appendix 2—figure 3, Appendix 2—figure 4, Appendix 2—figure 5, Appendix 2—figure 6, Appendix 2—figure 7, Appendix 2—figure 8, Supplementary file 1 — Table s3, Supplementary file 1 — Table s4, Supplementary file 1 — Table s5, Supplementary file 1 — Table s6).
Lung cancer
Meta-analysed IVW MR findings suggested that genetically predicted Intrinsic HorvathAge acceleration decreased the risk of lung cancer (OR = 0.97 per year increase in Intrinsic HorvathAge acceleration, 95% CI 0.93–1.00, p = 0.03) (Figure 2, Supplementary file 1 — Table s3, Supplementary file 1 — Table s4, Supplementary file 1 — Table s5, Supplementary file 1 — Table s6). However, these results did not survive multiple testing correction (FDR p = 0.21) and were not strongly supported by sensitivity analyses (Appendix 2—figure 1, Supplementary file 1 — Table s3, Supplementary file 1 — Table s4, Supplementary file 1 — Table s5, Supplementary file 1 — Table s6).
We did not find evidence of causality between other measures of epigenetic age acceleration and lung cancer (GrimAge IVW OR = 1.00, 95% CI 0.91–1.09, p = 0.96; PhenoAge IVW OR = 0.97, 95% CI 0.94–1.00, p = 0.06; and HannumAge IVW OR = 0.99, 95% CI 0.95–1.04, p = 0.82) (Figure 2, Appendix 2—figure 1, Appendix 2—figure 2, Appendix 2—figure 3, Appendix 2—figure 4, Appendix 2—figure 5, Appendix 2—figure 6, Appendix 2—figure 7, Appendix 2—figure 8, Supplementary file 1 — Table s3, Supplementary file 1 — Table s4, Supplementary file 1 — Table s5, Supplementary file 1 — Table s6).
Colorectal cancer
Meta-analysed IVW MR findings suggested that genetically predicted GrimAge acceleration increased the risk of colorectal cancer (OR = 1.12 per year increase in GrimAge acceleration, 95% CI 1.04–1.20, p = 0.002) (Figure 2, Supplementary file 1 — Table s3, Supplementary file 1 — Table s5, Supplementary file 1 — Table s6). These results survived multiple testing correction (FDR p = 0.04) and there was little evidence of heterogeneity across FinnGen and GECCO estimates (I2 = 0%, 95% CI ‘NA’, p = 0.61). Additionally, the direction of the genetically predicted effect was consistent across main and sensitivity MR analyses (i.e. MR-Egger, weighted median, and weighted mode) (Figure 3, Supplementary file 1 — Table s3, Supplementary file 1 — Table s5, Supplementary file 1 — Table s6) and results were consistent when using UK Biobank data alone (IVW OR = 1.15, 95% CI 1.04–1.28, p = 0.007) (Figure 2, Supplementary file 1 — Table s4).
We did not find evidence of residual population stratification in MR analyses using negative control outcomes (Appendix 2—figure 9, Supplementary file 1 — Table s7), nor did we find evidence of horizontal pleiotropy via potential colorectal cancer risk factors (Appendix 2—figure 10, Supplementary file 1 — Table s8).
Single-SNP analysis revealed that the effect was not driven by a single SNP (Supplementary file 1 — Table s9). Figure 4 shows the effect of genetic instruments on GrimAge acceleration against their effect on colorectal cancer. Moreover, there was no detectable evidence of uncorrelated horizontal pleiotropy (MR-Egger intercept = –0.13, 95% CI –0.33–0.07, p = 0.33), or heterogeneity across individual SNP estimates (Cochran’s Q 7.12, p = 0.07) (Supplementary file 1 — Table s10). We further explored the genetically predicted effect of GrimAge on colorectal cancer using GECCO data only and found no evidence against bias due to correlated pleiotropy (CAUSE OR = 1.00, 95% credible intervals 0.96–1.04, p = 1.00; shared q = 4%, 95% credible intervals 0–24%) (Appendix 2—figure 11).
Among subtypes, we found strong evidence for a causal relationship between GrimAge acceleration and colon cancer (IVW OR = 1.15, 95% CI 1.09–1.21, p = 0.006). In contrast, we did not find such evidence for rectal cancer (IVW OR = 1.05, 95% CI 0.97–1.13, p = 0.24). After further stratification, the magnitude of the genetically predicted effect of GrimAge acceleration on colon cancer was the same for distal (IVW OR = 1.16, 95% CI 1.03–1.29, p = 0.01) and proximal colon cancer (IVW OR = 1.16, 95% CI 0.97–1.40, p = 0.11). Also, sex-stratified results suggest that GrimAge acceleration may influence colorectal cancer in both males (IVW OR = 1.12, 95% CI 1.00–1.25, p = 0.05) and females (IVW OR = 1.14, 95% CI 1.04–1.26, p = 0.008) (Figure 5, Supplementary file 1 Table s11).
These findings were further supported by evidence of a positive association between GrimAge acceleration and parental history of colorectal cancer (OR = 1.06, 95% CI 1.00–1.12, p = 0.03) (Figure 6, Supplementary file 1 — Table s12). Additionally, LD Score regression coefficients for GrimAge acceleration and colorectal cancer were also in the expected direction (GECCO rg = 0.28, p < 0.001; UK Biobank rg = 0.15, p = 0.21; FinnGen rg = 0.27, p = 0.29) (Appendix 2—figure 12, Supplementary file 1 — Table s13).
We did not find consistent evidence of causality between other measures of epigenetic age acceleration and colorectal cancer (PhenoAge IVW OR = 1.00, 95% CI 0.97–1.02, p = 0.73; HannumAge IVW OR = 1.02, 95% CI 0.97–1.08, p = 0.37; and Intrinsic HorvathAge IVW OR = 1.00, 95% CI 0.97–1.02, p = 0.79) (Figure 2, Appendix 2—figure 1, Appendix 2—figure 2, Appendix 2—figure 3, Appendix 2—figure 4, Appendix 2—figure 5, Appendix 2—figure 6, Appendix 2—figure 7, Appendix 2—figure 8, Supplementary file 1 — Table s3, Supplementary file 1 — Table s4, Supplementary file 1 — Table s5, Supplementary file 1 — Table s6).
Telomere length
In bidirectional MR analyses, we found evidence that genetically predicted GrimAge acceleration may be a cause of telomere shortening (IVW beta coefficient = −0.07 per year increase in GrimAge acceleration, 95% CI –0.09 to –0.05, p < 0.001) and that genetically predicted longer telomere length may increase Intrinsic HorvathAge acceleration (IVW beta coefficient = 0.57 per standard deviation increase in telomere length, 95% CI 0.39–0.77, p = 0.002) (Appendix 2—figure 13, Supplementary file 1 — Table s14).
Steiger filtering showed that all genetic instruments for GrimAge acceleration were stronger predictors of GrimAge acceleration than telomere length. In contrast, it identified 20 genetic instruments for telomere length that were better predictors of Intrinsic HorvathAge acceleration than telomere length (Supplementary file 1 — Table s15). After removing these SNPs from the analyses, the results were still suggestive of an effect of telomere length on Intrinsic HorvathAge acceleration (IVW beta coefficient = 0.71 per standard deviation increase in telomere length, 95% CI 0.57–0.85, p < 0.001) (Appendix 2—figure 14, Supplementary file 1 — Table s16).
There was little evidence of causality between other measures of epigenetic age acceleration and telomere length (Appendix 2—figure 13, Supplementary file 1 — Table s14).
Discussion
In this comprehensive two-sample MR study of epigenetic age acceleration and multiple cancers, we found evidence to suggest that genetically predicted GrimAge acceleration may increase the risk of colorectal cancer in both males and females. Among subtypes, effects appeared to be stronger in relation to colon than rectal cancer. Our MR results also suggested that genetically predicted GrimAge acceleration may decrease the risk of prostate cancer and that genetically predicted Intrinsic HorvathAge acceleration may be protective against lung cancer. Nevertheless, these did not pass multiple testing correction. Finally, we found no consistent evidence for other measures of epigenetic age acceleration and cancers.
Our MR estimates for the association between GrimAge and colorectal cancer were consistent with those reported in Dugue et al., 2021, an observational nested case-control study in the Melbourne Collaborative Cohort Study (RR = 1.04 per year increase in GrimAge acceleration, 95% CI 1.01–1.07, p = 0.02). However, our findings contrast with those highlighted in Hillary et al., 2020, an observational cohort study that used Generation Scotland data. The latter authors observed no evidence of an association between GrimAge acceleration and colorectal cancer after correcting for multiple testing. Nevertheless, it is possible that their analyses were underpowered, as their sample only included 63 colorectal cancer cases (0.66%). More importantly, the direction of the reported estimate is consistent with our findings and those presented in Dugue et al., 2021.
Observational evidence for the association between other measures of epigenetic ageing and cancer is inconclusive (the pre-existing evidence has been summarised in Supplementary file 1 — Table s17). For instance, epigenetic clock acceleration has been positively associated with breast (Ambatipudi et al., 2017; Kresovich et al., 2019b; Kresovich et al., 2019a) and lung cancer (Levine et al., 2018; Levine et al., 2015; Dugue et al., 2021) in some studies. However, (Durso et al., 2017, Hillary et al., 2020) and (Dugué et al., 2018) did not find strong evidence to support this. In some cases, observational evidence is stronger for some clocks than it is for others. For example, for colorectal cancer, evidence of a positive association is much stronger for second-generation clocks (Dugue et al., 2021) than for first-generation clocks (Dugué et al., 2018; Durso et al., 2017). In the case of prostate cancer, as in our study, apart from weak evidence of an inverse association with GrimAge, no other associations have been observed (Dugue et al., 2021; Dugué et al., 2018). To date, the association between epigenetic age acceleration and ovarian cancer has not been explored observationally. Although our findings were less susceptible to biases that often influence observational research, they still provide no compelling evidence of a causality between several measures of epigenetic clock acceleration and cancer.
This MR study had several strengths. For instance, we pooled results from multiple sources using fixed effect meta-analysis to improve the precision of the MR estimates presented in McCartney et al., 2021. We also conducted extra sensitivity analyses, such as MR of negative control outcomes, MR of cancer risk factors, single-SNP MR and CAUSE analyses, to assess the validity of the MR assumptions. Moreover, we performed subtype-specific MR analyses and sought to corroborate our results using UK Biobank GWAS data on parental history of cancer and LD Score regression. Additionally, our findings contribute to the identification of modifiable targets for future interventions aimed at reversing epigenetic ageing for the prevention of cancer. Compared to clinical trials, MR provides a cheaper, quicker, and ethical way of assessing the long-term impact of interventions on epigenetic ageing. This is especially relevant while attempts to develop interventions which reverse epigenetic ageing are still in early stages (Fahy et al., 2019; Fitzgerald et al., 2021; Gensous et al., 2020; Chen et al., 2019).
The findings from this study should be interpreted in light of its limitations. We only identified four genetic instruments for GrimAge acceleration, which explained 0.47% of the variance in the trait. This could lead to two issues: low statistical power and horizontal pleiotropy. First, our GrimAge analyses were underpowered to detect ORs < 1.20 for colorectal cancer. Therefore, it is possible that our findings do not reflect a true effect (we identified an OR = 1.12 for colorectal). Similarly, our study was underpowered to detect genetically predicted effects of GrimAge acceleration on cancer subtypes and cancers with smaller sample sizes (i.e. ovarian and lung cancer). Some of our sensitivity analyses, such as the MR-Egger intercept test used to detect uncorrelated horizontal pleiotropy, also had low power, resulting in imprecise estimates. The weighted mode method may also be misleading in this context, as its use is limited in the presence of very few SNPs. Although these limitations potentially undermine the validity of our results, it is reassuring that point estimates for the genetically predicted effect of GrimAge acceleration on colorectal cancer were consistent across MR methods and study populations. However, since CAUSE analyses did not provide evidence against confounding by correlated horizontal pleiotropy, it is possible that the genetically predicted effects identified are attributed to correlated pleiotropy (whereby SNPs are associated with epigenetic age acceleration and cancer through a shared heritable factor) rather than a causal effect of GrimAge on cancer risk.
One could argue that because the results for GrimAge acceleration were inconsistent with those obtained for other measures of epigenetic age acceleration, chance and horizontal pleiotropy are more likely explanations for our findings. However, inconsistencies across epigenetic ageing measures do not necessarily invalidate our results. They may simply reflect differences in how clocks were trained (i.e. they were trained on different outcomes, tissues, and populations). Different clocks may capture information on distinct underlying biological ageing mechanisms (Liu et al., 2020). For example, GrimAge was trained on mortality and smoking (factors which are closely related to cancer risk), which may explain why it outperforms other measures of epigenetic ageing in predicting time-to-cancer (Lu et al., 2019a).
Although little is known about the underlying mechanisms, GrimAge may plausibly influence cancer risk through hormonal, inflammatory and metabolic processes (Yu et al., 2020; Bottazzi et al., 2018; Lau and Robinson, 2021). In bidirectional MR analyses, we found evidence that genetically predicted GrimAge acceleration may be a cause of telomere shortening, another marker of biological ageing. Shorter telomeres have been shown to lower cancer risk in prior MR analyses (Telomeres Mendelian Randomization Collaboration et al., 2017; Gao et al., 2020; Kuo et al., 2019), so it is plausible that GrimAge acceleration decreases cancer risk, at least in part, via its effect on telomere length. GrimAge acceleration may still increase cancer risk via pathways other than those related to cellular division. The positive effect of GrimAge acceleration on cancer via these other pathways may counteract the negative effects mediated via telomere length, resulting in null MR results for GrimAge acceleration and breast, ovarian, prostate and lung cancer, and positive MR results for GrimAge acceleration and colorectal cancer. To better understand the biology of ageing, future studies should consider running MR analyses using data on DNAm-predicted telomere length, since DNAm telomere length is independent of telomerase activity and has been more strongly associated with age than measured telomere length (Lu et al., 2019b).
Although promising in terms of consistency and biological plausibility, further research is required to confirm our findings. For example, multivariable MR (Burgess and Thompson, 2015; Sanderson et al., 2019) could be used to disentangle the causal effects of GrimAge acceleration on cancer from shared heritable factors such as and blood cell composition. Additionally, our analyses could be replicated using other large independent cancer datasets. Although we conducted MR analyses on parental history of cancer, their effect estimates are not directly comparable to those obtained in the main analyses due to cases in the GWAS-by-proxy of parental endpoints being defined as either or both parents reportedly having a type of cancer. Furthermore, it would also be useful to replicate our analyses once a larger GWAS of epigenetic ageing with more genetic instruments for GrimAge acceleration is available. This would allow for a more rigorous assessment of horizontal pleiotropy and may be used to assess clustering of genetic variants to reveal distinct biological mechanisms underlying the effects (Foley et al., 2021). In spite of these suggestions, we acknowledge that it may be challenging to get access to suitable datasets for replication purposes in the short term.
The selection of ‘super controls’ (e.g. in UK Biobank, FinnGen and GECCO), with no other cancers, related lesions (i.e. benign, in situ, uncertain or unspecified behaviour neoplasms) or reported family history of cancer, could have inflated cancer GWAS effect sizes (and our MR estimates), because ‘super controls’ are healthier than the general population and are less likely to be genetically predisposed to develop cancer.
Another limitation is that we did not have access to individual level data. Therefore, we were unable to stratify the analyses by potential effect modifiers, such as sex, smoking, and menopausal status. Moreover, we did not have sex-specific instruments for sex-specific cancers. However, it is unlikely that the genetic architecture of epigenetic clock acceleration differs across sexes, as DNAm levels at individual clock CpGs are highly correlated between males and females (Grodstein et al., 2020; Tajuddin et al., 2019).
Finally, to reduce bias due to population stratification, this study was conducted using data from participants of European ancestry only. The GWAS data used for the analyses had been adjusted for the top genetic principal components for the same reason. Furthermore, our MR of negative control outcomes suggests that our MR results are unlikely to be biased by residual population stratification. Despite this, confounding due to population stratification, dynastic effects and assortative mating cannot be ruled out completely, as it is not possible to test the second MR assumption (i.e. independence assumption). Furthermore, more research is required to see if our results could translate to other ancestries.
From a public health perspective, our work provides potentially relevant findings. Observational and Mendelian randomization studies suggest that GrimAge acceleration may be influenced by several cancer risk factors, such as obesity and smoking (Lu et al., 2019a; McCartney et al., 2021). If GrimAge acceleration is a causal mediator between these risk factors and colorectal cancer, the GrimAge clock may be a treatable intermediary when targeting the underlying risk factors is not feasible or too difficult to accomplish. It could also be targeted in populations at high-risk of colorectal cancer. Nevertheless, we think it may be too early to make claims regarding the clinical utility of our findings. The GrimAge clock has only recently been created and very few studies have assessed its association with colorectal cancer. More research is required to corroborate our results and to evaluate whether GrimAge acceleration can be modified through lifestyle or clinical interventions.
In conclusion, our findings suggest that genetically predicted GrimAge acceleration may increase the risk of colorectal cancer. Findings were less consistent for other epigenetic clocks and cancers. Further work is required to investigate the potential mechanisms underlying the genetically predicted effects identified in this study.
Appendix 1
Additional Methods
Cancer datasets
UK Biobank
The UK Biobank is a large cohort study including around 500,000 individuals aged 40–69 years at the time of recruitment (2006–2010). The cohort has been described in detail in previous publications (Sudlow et al., 2015; Bycroft et al., 2018). In short, all participants provided written informed consent, after which baseline data were collected using sociodemographic, lifestyle and health-related questionnaires, physical and cognitive assessments, and biological samples. Participants’ data were linked to their health records for longitudinal follow-up. The study obtained ethical approval from the National Information Governance Board for Health and Social Care and the North-West Multicenter Research Ethics Committee (Ref: 11/NW/0382).
Cancer cases (diagnosed prior or after enrolment) were obtained from the UK Cancer Registry (updated to April 2019). They were then coded according to the ninth and tenth editions of the International Classification of Diseases (ICD-9 and ICD-10, respectively) as follows: breast (ICD-9: 174; ICD-10: C50), ovarian (ICD-9: 183; ICD-10: C56), prostate (ICD-9: 185; ICD-10: C61), lung (ICD-9: 162; ICD-10: C34) and colorectal cancer (ICD-9: 153; ICD-10: C18-C20). Controls excluded individuals with any type of cancer (self-reported and/or recorded in cancer registry), as well as those with benign, in situ, uncertain or unspecified behaviour neoplasms (ICD-9: 210–239; ICD-10: D00-D49).
Sample-level quality control (QC) involved removing any individuals who had non-white British genetic ancestry, sex chromosome aneuploidies, who withdrew consent from the UK Biobank study and who were closely related to other participants. Variant-level QC consisted in imputing SNPs using the Haplotype Reference Consortium (HRC) and restricting SNPs to a minor allele frequency (MAF) > 0.1%, a genotyping rate > 0.015 and a Hardy-Weinberg Equilibrium (HWE) P > 1 × 10–4. LD pruning was performed to an r2 cutoff of 0.1 using PLINK v2 (Mitchell et al., 2019). In order to reduce false positive signals, SNPs were removed when MAF was below our expectations (we would expect at least 25 minor alleles in cases), as recommended in http://www.nealelab.is/blog/2017/9/11/details-and-considerations-of-the-uk-biobank-gwas.
The GWAS analysis in the UK Biobank consisted of 13,879 cases and 198,523 controls for breast cancer, 1,218 cases and 198,523 controls for ovarian cancer, 9,132 cases and 173,493 controls for prostate cancer, 2,671 cases and 372,016 controls for lung cancer and 5,657 cases and 372,016 controls for colorectal cancer. It was performed using BOLT-LMM v2.3.5 (Loh et al., 2015; Elsworth et al., 2019), adjusting for sex and genotyping chip. BOLT-LMM uses a linear mixed model to account for population stratification and cryptic relatedness in the UK Biobank. Lung cancer associations were estimated twice, once adjusting for genotyping chip and once without. Since UKBiLEVE participants were genotyped using a different array and using adjusted lung cancer estimates may introduce collider bias, we only included MR results obtained using the unadjusted lung cancer GWAS estimates in the meta-analysis. For sex-specific cancers, analyses were limited to individuals of the pertinent sex (only females were used for breast and ovarian cancers, whereas only males were used for prostate cancer). Beta coefficients and their corresponding standard errors were finally transformed to log odds ratios (ORs) (Elsworth et al., 2019).
We also performed a GWAS analysis of parental history of cancer reported by UK Biobank participants (i.e. breast, prostate, lung and bowel cancer) using BOLT-LMM software v2.3.5 (Loh et al., 2015). Age and sex were included as covariates in the model as before. For sex-specific cancers, analyses were restricted to individuals of the relevant sex (i.e. maternal history only for breast cancer and paternal history only for prostate cancer). We obtained 35,356 breast cancer cases and 206,992 controls, in addition to 31,527 prostate cancer cases and 160,579 controls. For other cancers, we combined maternal and paternal history of cancer, thus obtaining a total of 51,073 lung cancer cases and 404,606 controls, as well as 45,213 bowel cancer cases and 412,429 controls. GWAS of these outcomes have previously provided strong concordance with those based on hospital records (DeBoever et al., 2020). They have also provided consistent results in MR (Richardson et al., 2021).
FinnGen
The FinnGen R5 release includes data on 218,792 individuals of Finnish ancestry, obtained from Finnish biobanks and digital health registry records (FinnGen, 2021). Complete study details are available elsewhere (https://www.finngen.fi/en). In brief, samples were excluded for the following reasons: ambiguous gender, genotype missingness > 5%, heterozygosity +–4 s.d. and non-Finnish ancestry. SNPs were genotyped using Illumina and Affymetrix arrays. Variants were excluded for the following reasons: missingness > 2%, HWE P < 1 × 10-6 and minor allele count < 3. Genotypes were imputed using the Finnish SISu v3 reference panel. The GWAS analysis was conducted using SAIGE v0.36.3.2, a mixed model logistic regression R/C ++ package. Sex, age, genotyping batch and the first 10 genetically derived principal components were included as covariates in the analysis. We used FinnGen R5 release data on breast (8,401 cases and 99,321 controls), ovarian (719 cases and 99,321 controls), prostate (6,311 cases and 74,685 controls), lung (1,681 cases and 173,933 controls) and colorectal cancer (3,022 cases and 174,006 controls). We used the “EXALLC” cancer variables, which excluded other cancers from controls.
Breast Cancer Association Consortium
The GWAS summary data for breast cancer were obtained from a Breast Cancer Association Consortium (BCAC) meta-analysis performed by Michailidou et al., 2017. This included 122,977 cases and 105,974 controls (69,501 cases of ER + and 21,468 of ER- breast cancer). All studies that contributed to this meta-analysis have been fully detailed in previous publications (Michailidou et al., 2017; Michailidou et al., 2013; Michailidou et al., 2015). In sum, samples were excluded if they had a low call rate ( < 95%), abnormally high or low heterozygosity (4.89 s.d. from the mean), < 80% European ancestry, probable duplicates and/or close relatives within and across studies. Genetic variants were genotyped using the Illumina OncoArray and iCOGS arrays and genotypes were imputed using the 1,000 Genomes Project Phase three reference panel. The GWAS analysis was performed using logistic regression models, adjusting for up to 10 principal components and either country or study. This was done using purpose-written software. OncoArray and iCOGS estimates were combined in a fixed-effect inverse variance weighted meta-analysis using the METAL software (Willer et al., 2010). Only SNPs with r2 ≥ 0.3 and MAF ≥ 0.005 were included in the meta-analysis.
We also obtained summary data for breast cancer subtypes from a BCAC GWAS meta-analysis by Zhang et al., 2020. The study comprised data on luminal A-like (7,325 cases), luminal B-like (1,682 cases), luminal B/HER2-negative-like (1,779 cases), HER2-enriched-like (718 cases) and triple-negative (2,006 cases) invasive breast cancer subtypes and 20,815 controls. The details of the study can be found in the publication. In brief, the analyses excluded cases of carcinoma in situ, cases missing data on tumour characteristics and cases for which there were no controls available in their respective countries. Participants were also excluded if age at diagnosis/enrolment was missing. Genotypes were obtained using OncoArray and iCOGS arrays. Imputation was performed using the 1,000 Genomes Project Phase three reference panel. OncoArray and iCOGS datasets were analysed separately using two-stage polytomous logistic regression analyses in R. Models were adjusted for age and the first 10 principal components. SNPs with r2 < 0.3 and MAF < 0.01 were excluded from the subtype analyses, as well as those in linkage disequilibrium (r2 ≥ 0.1) or within ± 500 kb of known susceptibility SNPs. GWAS results were then pooled using fixed-effect meta-analysis in METAL (Willer et al., 2010).
Ovarian Cancer Association Consortium
We used ovarian cancer genetic summary statistics from an Ovarian Cancer Association Consortium (OCAC) study by Phelan et al., 2017. This comprised 25,509 cases and 40,941 controls. Subtypes included high grade serous (13,037 cases), low grade serous (1,012 cases), invasive mucinous (1,417 cases), clear cell (1,366 cases) and endometrioid (2,810 cases) ovarian cancers. This study combined genotype data from OCAC and Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) genotyping projects. These have been fully described in the publication. In short, samples with > 27% non-European ancestry were excluded, as were those with a genotyping call rate < 95%, excessively low or high heterozygosity. Non-females and duplicates were also removed. SNPs were genotyped using several Illumina arrays (OncoArray, iSelect iCOGS, 550 k, HumanOmni 2.5 M, 610 Quad and 317 k). Imputations were performed separately for each genotyping project using the 1,000 Genomes Project v3 reference panel. GWAS analyses were conducted in custom-written software using logistic regression models adjusted for study and principal components. SNPs with r2 < 0.3 and MAF < 0.01 were excluded. GWAS estimates were pooled using fixed effect meta-analysis in METAL (Willer et al., 2010).
Consortium of Investigators of Modifiers of BRCA1/2
We also used CIMBA GWAS data for breast and ovarian cancers in BRCA1 and BRCA2 mutation carriers (Phelan et al., 2017; Milne et al., 2017). The genotyping and imputation procedures that were used have been described elsewhere (Phelan et al., 2017; Milne et al., 2017). In brief, samples were excluded if they were non-female, had discordant genotypes in known sample duplicates, had > 19% non-European ancestry, a genotyping call rate < 95% or extremely low or high heterozygosity (P < 1 × 10−6). SNPs were genotyped using Illumina’s Oncoarray and iSelect Collaborative Oncological Gene-Environment Study (iCOGS) arrays. Imputation was performed using the 1,000 Genomes Project Phase three reference panel. GWAS analyses were conducted separately for BRCA1 and BRCA2 mutation carriers and for iCOGS and OncoArray samples. Genetic association data was generated using a survival analysis framework, using a retrospective likelihood approach. Analyses were stratified by country of origin and Ashkenazi Jewish origin. Custom written functions in Fortran and Python were used to carry out the analyses and kinship-adjusted score test statistics were implemented in R software. OncoArray and iCOGS results were pooled using fixed-effect meta-analysis in METAL (Willer et al., 2010).
Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome
Prostate cancer GWAS summary data were acquired from a Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (PRACTICAL) study by Schumacher et al., 2018. This included 79,148 cases and 61,106 controls. It also comprised data on prostate cancer subtypes: 15,167 advanced cases vs. 58,308 healthy controls; 14,160 advanced cases vs. 62,421 non-advanced controls; 6,988 early-onset cases (age at diagnosis ≤ 55 years) vs. 44,256 healthy controls; 15,561 high aggressive cases vs. 9,739 low aggressive controls; and 20,658 high aggressive cases vs. 38,093 low/intermediate aggressive controls.
Prostate cancer aggressiveness was defined as follows:
Low aggressive: tumor stage ≤T1 AND Gleason score ≤6 AND prostate-specific antigen (PSA) <10 ng/mL.
Intermediate aggressive: tumor stage T2 OR Gleason score = 7 OR PSA 10–20 ng/mL.
High aggressive: tumor stage T3/T4, N1 or M1 OR Gleason score ≥8 OR PSA >20 ng/mL.
Advanced: metastatic disease OR Gleason score ≥ 8 OR PSA > 100 ng/mL OR death due to prostate cancer.
Study details are available in the publication. In brief, individuals were excluded if they presented a call rate < 95%, extreme heterozygosity ( > 4.9 s.d. from the mean), if they were duplicates or if they were related to other participants. Only men of European ancestry ( > 80%) were included in the GWAS. Studies were genotyped using Illumina (OncoArray, Human 610, 60 k, Infinium HumanHap 550, iSELECT, iCOGS and Human Omni 2.5) and Affymetrix GeneChip (500 k and 5.0 k) genotyping arrays and SNPs were imputed to the 1,000 Genomes Project Phase three reference panel. Genetic association data were obtained using logistic regression analysis. Models were adjusted for seven principal components and study-relevant covariates and stratified by country or study. Odds ratios were derived using either SNPTEST or a custom written C ++ software. GWAS estimates were combined using fixed-effect meta-analysis in METAL (Willer et al., 2010).
International Lung Cancer Consortium
For lung cancer, we used genetic summary data obtained from an International Lung Cancer Consortium (ILCCO) GWAS meta-analysis of 11,348 cases and 15,861 controls by Wang et al., 2014. We also used lung cancer subtype data including 3,275 squamous cell lung carcinoma cases and 15,038 controls, as well as 3,442 lung adenocarcinoma cases and 14,894 controls. Individual studies included in the meta-analysis have been explained in prior publications (Timofeeva et al., 2012; Amos et al., 2008; Wang et al., 2008; Hung et al., 2008). In summary, sample QC consisted in excluding any individuals of non-European ancestry, with low call rates ( < 90%), or abnormally high or low heterozygosity (P < 1 × 10–4). Duplicates and closely related individuals were also removed. Genotyping was performed using Illumina HumanHap 317 k, 317k + 240 S, 370Duo, 550 k, 610 k or 1 M arrays. SNPs were imputed from the 1,000 Genomes Project Phase 1 v3 reference panel. GWAS estimates were obtained by unconditional logistic regression in R v2.6, Stata v.10 and PLINK v1.06 software. Analyses were adjusted for principal components. Fixed-effect meta-analysis was used to pool estimates across studies.
Genetics and Epidemiology of Colorectal Cancer Consortium
Colorectal cancer GWAS summary statistics were retrieved from a Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO) GWAS meta-analysis by Huyghe et al., 2019. This comprised 58,131 cases (31,288 male and 26,843 female) and 67,347 controls (34,527 male and 32,820 female). Cases were defined as patients with colorectal cancer or advanced adenoma.
Data on colorectal cancer subtypes were obtained from another GECCO publication by Huyghe et al., 2021. This included 48,214 cases and 64,159 controls (32,002 colon, 15,706 proximal colon, 14,376 distal colon and 16,212 rectal cancer cases).
Colorectal cancer subtypes were defined as follows:
Proximal colon cancer: any primary tumour starting in the cecum, ascending colon, hepatic flexure, or transverse colon (ICD-9: 153.4, 153.6, 153.0, or 153.1, respectively).
Distal colon cancer: any primary tumour starting in the splenic flexure, descending colon, or sigmoid colon (ICD-9 codes: 153.7, 153.2, or 153.3, respectively)
Colon cancer: proximal and distal colon cancer cases, in addition to colon cancer cases with unspecified site.
Rectal cancer: any primary tumour starting in the rectum or rectosigmoid junction (ICD-9 codes: 154.1, or 154.0, respectively)
Controls excluded individuals with known history of cancer or reported family history of colorectal cancer. QC procedures have been explained in the publications (Huyghe et al., 2019; Huyghe et al., 2021). In brief, the studies excluded samples with evidence of DNA contamination, high missing genotype rates, unintentional duplicate pairs and sex discrepancies. Closely related individuals and those of non-European ancestry were also excluded. Genotyping was conducted using Illumina (300 k, Oncoarray, 1 M, 550 k, 610 k, OmniExpress, OmniExpressExome, 300/240 S and custom iSelect) and Affymetrix (Axiom and 500 k) arrays. Imputation was performed to the HRC reference panel. GWAS analyses were conducted for SNPs with an imputation accuracy r2 ≥ 0.3 and minor allele count ≥ 50 using logistic regression models adjusted for principal components, age, sex and study-specific covariates. METAL (Willer et al., 2010) was used to combine summary statistics across studies using fixed-effect meta-analysis.
Sensitivity analyses
MR-Egger assumes that the association between SNPs and epigenetic age acceleration is not correlated with SNPs that affect cancer via pleiotropic pathways (Instrument Strength Independent of Direct Effect—InSIDE assumption) (Bowden et al., 2015). The weighted median method assumes that at least half of the SNPs in the analysis are valid instruments. The weighted mode approach presupposes that the most frequent association estimate is not affected by pleiotropy, meaning it must correspond to the true causal effect (ZEro Modal Pleiotropy Assumption—ZEMPA) (Hartwig et al., 2017).
Data availability
Summary statistics for epigenetic age acceleration measures of HannumAge, Intrinsic HorvathAge, PhenoAge and GrimAge were downloaded from: https://datashareed.ac.uk/handle/10283/3645. Summary statistics for international cancer genetic consortiums were obtained from their respective data repositories. Colorectal cancer data were obtained following the submission of a written request to the GECCO committee, which may be contacted by email at kafdem@fredhutch.org/upeters@fredhutch.org. Breast, ovarian, prostate and lung cancer data were accessed via MR-Base (http://app.mrbase.org/), which holds complete GWAS summary data from BCAC, OCAC, PRACTICAL and ILCCO. Breast cancer subtype data were obtained from BCAC and can be downloaded from: http://bcacccgemedschlcam.ac.uk/bcacdata/oncoarray/oncoarray-and-combined-summary-result/gwas-summary-associations-breast-cancer-risk-2020/. Data on breast and ovarian cancer in BRCA1 and BRCA2 carriers were obtained from CIMBA and can be downloaded from: http://cimbaccgemedschlcam.ac.uk/oncoarray-complete-summary-results/. Prostate cancer subtype data are not publicly available through MR-Base but can be accessed upon request. These data are managed by the PRACTICAL committee, which may be contacted by email at practical@icr.ac.uk. FinnGen data is publicly available and can be accessed here: https://www.finngen.fi/en/access_results. UK Biobank data can be accessed through the MR-Base platform. Parental history of cancer data were obtained from the UK Biobank study under application #15825 and can be accessed via an approved application to the UK Biobank (https://www.ukbiobank.ac.uk/enable-your-research/apply-for-access). GWAS data for negative control outcomes and potential confounders were obtained via the MR-Base platform (GWAS IDs for negative control outcomes: "ukb-b-19560", "ukb-b-533"; GWAS IDs for confounders: "ukb-b-10831", "ukb-b-13702", "ukb-b-6134", "ieu-a-1239", "ukb-b-5779", "ieu-a-835", "ieu-a-61"). GWAS data for measured telomere length used in bidirectional MR analyses were also obtained via the MR-Base platform (GWAS ID: “ieu-b-4879”).
Code availability
The GWAS analysis for cancers in UK Biobank was performed using BOLT-LMM v2.3.5 (http://data.broadinstitute.org/alkesgroup/BOLT-LMM/). All MR analyses and visualisations were conducted using R software v4.0.2 (https://www.r-project.org/). For cancer datasets that were not obtained from the MR-Base platform, LD proxies were identified using the ‘LDlinkR’ v1.1.2 R package (https://github.com/CBIIT/LDlinkR; Myers, 2021). Two-sample MR analyses were conducted using the ‘TwoSampleMR’ v0.5.6 R package (https://github.com/MRCIEU/TwoSampleMR; Parker, 2021). Meta-analyses were performed using the ‘meta’ v4.18 R package (https://github.com/guido-s/meta/; Schwarzer, 2022). GWAS meta-analyses used to perform single-SNP MR analyses were done using the METAL software (https://genome.sph.umich.edu/wiki/METAL_Documentation). MR-CAUSE analyses were conducted using the ‘cause’ v1.2.0 R package (https://github.com/jean997/cause; Morrison, 2021). Plots were created using the ‘ggforestplot’ v0.1.0 R package (https://github.com/nightingalehealth/ggforestplot; Jagerroos, 2020). LD scores were computed using the ‘ldsc’ v1.0.1 command line tool (https://github.com/bulik/ldsc; Schorsch, 2020). The code used in this study is available at: https://github.com/fernandam93/epiclocks_cancer; Morales Berstein, 2021.
Appendix 2
Additional Figures
Data availability
Summary statistics for epigenetic age acceleration measures of HannumAge, Intrinsic HorvathAge, PhenoAge and GrimAge were downloaded from: https://datashare.ed.ac.uk/handle/10283/3645 (Datasets used: European-ancestries meta-analysis summary statistics: Hannum (645.4Mb), European-ancestries meta-analysis summary statistics: IEAA (645.7Mb), European-ancestries meta-analysis summary statistics: GrimAge (645.7Mb), European-ancestries meta-analysis summary statistics: PhenoAge (645.7Mb)). Summary statistics for international cancer genetic consortiums were obtained from their respective data repositories. Colorectal cancer data were obtained following the submission of a written request to the GECCO committee, which may be contacted by email at kafdem@fredhutch.org/upeters@fredhutch.org. Breast, ovarian, prostate and lung cancer data were accessed via MR-Base (http://app.mrbase.org/), which holds complete GWAS summary data from BCAC, OCAC, PRACTICAL and ILCCO. Breast cancer subtype data were obtained from BCAC and can be downloaded from: http://bcac.ccge.medschl.cam.ac.uk/bcacdata/oncoarray/oncoarray-and-combined-summary-result/gwas-summary-associations-breast-cancer-risk-2020/. Data on breast and ovarian cancer in BRCA1 and BRCA2 carriers were obtained from CIMBA and can be downloaded from: http://cimba.ccge.medschl.cam.ac.uk/oncoarray-complete-summary-results/. Prostate cancer subtype data are not publicly available through MR-Base but can be accessed upon request. These data are managed by the PRACTICAL committee, which may be contacted by email at practical@icr.ac.uk. FinnGen data is publicly available and can be accessed here: https://www.finngen.fi/en/access_results (Datasets used from release 5: Malignant neoplasm of breast (all cancers excluded), Malignant neoplasm of bronchus and lung (all cancers excluded), Colorectal cancer (all cancers excluded), Malignant neoplasm of ovary (all cancers excluded), Malignant neoplasm of prostate (all cancers excluded)). UK Biobank data can be accessed through the MR-Base platform. Parental history of cancer data were obtained from the UK Biobank study under application #15825 and can be accessed via an approved application to the UK Biobank (https://www.ukbiobank.ac.uk/enable-your-research/apply-for-access). GWAS data for negative control outcomes and potential confounders were obtained via the MR-Base platform (GWAS IDs for negative control outcomes: "ukb-b-19560", "ukb-b-533"; GWAS IDs for confounders: "ukb-b-10831", "ukb-b-13702", "ukb-b-6134", "ieu-a-1239", "ukb-b-5779", "ieu-a-835", "ieu-a-61"). GWAS data for measured telomere length used in bidirectional MR analyses were also obtained via the MR-Base platform (GWAS ID: "ieu-b-4879").
-
Edinburgh DataShareGenome-wide association studies identify 137 loci for DNA methylation biomarkers of aging.https://doi.org/10.7488/ds/2834
-
IEU OpenGWASID ieu-a-1126. Breast cancer (Combined Oncoarray; iCOGS; GWAS meta analysis).
-
IEU OpenGWASID ieu-a-1121. High grade serous ovarian cancer.
-
IEU OpenGWASID ieu-a-1122. Low grade serous ovarian cancer.
-
IEU OpenGWASID ieu-a-1123. Invasive mucinous ovarian cancer.
-
IEU OpenGWASID ieu-a-1124. Clear cell ovarian cancer.
-
IEU OpenGWASID ieu-a-1125. Endometrioid ovarian cancer.
-
IEU OpenGWASID ieu-a-1127. ER+ Breast cancer (Combined Oncoarray; iCOGS; GWAS meta analysis).
-
IEU OpenGWASID ieu-a-1128. ER- Breast cancer (Combined Oncoarray; iCOGS; GWAS meta analysis).
-
IEU OpenGWASID ukb-b-13702. Time spent doing vigorous physical activity.
-
IEU OpenGWASID ukb-b-6134. Age completed full time education.
-
IEU OpenGWASID ukb-b-5779. Alcohol intake frequency.
References
-
Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple TestingJournal of the Royal Statistical Society 57:289–300.https://doi.org/10.1111/j.2517-6161.1995.tb02031.x
-
Aging, inflammation and cancerSeminars in Immunology 40:74–82.https://doi.org/10.1016/j.smim.2018.10.011
-
Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regressionInternational Journal of Epidemiology 44:512–525.https://doi.org/10.1093/ije/dyv080
-
A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomizationStatistics in Medicine 36:1783–1802.https://doi.org/10.1002/sim.7221
-
Improving the accuracy of two-sample summary-data Mendelian randomization: moving beyond the NOME assumptionInternational Journal of Epidemiology 48:728–742.https://doi.org/10.1093/ije/dyy258
-
An atlas of genetic correlations across human diseases and traitsNature Genetics 47:1236–1241.https://doi.org/10.1038/ng.3406
-
Bias in causal estimates from Mendelian randomization studies with weak instrumentsStatistics in Medicine 30:1312–1323.https://doi.org/10.1002/sim.4197
-
Mendelian randomization analysis with multiple genetic variants using summarized dataGenetic Epidemiology 37:658–665.https://doi.org/10.1002/gepi.21758
-
Multivariable Mendelian randomization: the use of pleiotropic genetic variants to estimate causal effectsAmerican Journal of Epidemiology 181:251–260.https://doi.org/10.1093/aje/kwu283
-
Effects of Vitamin D3 Supplementation on Epigenetic Aging in Overweight and Obese African Americans With Suboptimal Vitamin D Status: A Randomized Clinical TrialThe Journals of Gerontology. Series A, Biological Sciences and Medical Sciences 74:91–98.https://doi.org/10.1093/gerona/gly223
-
‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of diseaseInternational Journal of Epidemiology 32:1–22.https://doi.org/10.1093/ije/dyg070
-
Reading Mendelian randomisation studies: a guide, glossary, and checklist for cliniciansBMJ (Clinical Research Ed.) 362:k601.https://doi.org/10.1136/bmj.k601
-
Assessing Digital Phenotyping to Enhance Genetic Studies of Human DiseasesAmerican Journal of Human Genetics 106:611–622.https://doi.org/10.1016/j.ajhg.2020.03.007
-
Mendelian randomization as an instrumental variable approach to causal inferenceStatistical Methods in Medical Research 16:309–330.https://doi.org/10.1177/0962280206077743
-
DNA methylation-based biological aging and cancer risk and survival: Pooled analysis of seven prospective studiesInternational Journal of Cancer 142:1611–1619.https://doi.org/10.1002/ijc.31189
-
MR-Clust: clustering of genetic variants in Mendelian randomization with similar causal estimatesBioinformatics (Oxford, England) 37:531–541.https://doi.org/10.1093/bioinformatics/btaa778
-
Characteristics of Epigenetic Clocks Across Blood and Brain Tissue in Older Women and MenFrontiers in Neuroscience 14:555307.https://doi.org/10.3389/fnins.2020.555307
-
Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumptionInternational Journal of Epidemiology 46:1985–1998.https://doi.org/10.1093/ije/dyx102
-
Methylation-Based Biological Age and Breast Cancer RiskJournal of the National Cancer Institute 111:1051–1058.https://doi.org/10.1093/jnci/djz020
-
DNA methylation age as a biomarker for cancerInternational Journal of Cancer 148:2652–2663.https://doi.org/10.1002/ijc.33451
-
GrimAge Outperforms Other Epigenetic Clocks in the Prediction of Age-Related Clinical Phenotypes and All-Cause MortalityThe Journals of Gerontology. Series A, Biological Sciences and Medical Sciences 76:741–749.https://doi.org/10.1093/gerona/glaa286
-
SoftwareUK Biobank Genetic Data: MRC-IEU Quality ControlVersion 2.
-
Using multiple genetic variants as instrumental variables for modifiable risk factorsStatistical Methods in Medical Research 21:223–242.https://doi.org/10.1177/0962280210394459
-
Two-step epigenetic Mendelian randomization: a strategy for establishing the causal role of epigenetic processes in pathways to diseaseInternational Journal of Epidemiology 41:161–176.https://doi.org/10.1093/ije/dyr233
-
An examination of multivariable Mendelian randomization in the single-sample and two-sample summary data settingsInternational Journal of Epidemiology 48:713–727.https://doi.org/10.1093/ije/dyy262
-
The use of negative control outcomes in Mendelian randomization to detect potential population stratificationInternational Journal of Epidemiology 50:1350–1361.https://doi.org/10.1093/ije/dyaa288
-
Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 CountriesCA: A Cancer Journal for Clinicians 71:209–249.https://doi.org/10.3322/caac.21660
-
Influence of common genetic variation on lung cancer risk: meta-analysis of 14 900 cases and 29 485 controlsHuman Molecular Genetics 21:4980–4995.https://doi.org/10.1093/hmg/dds334
-
The heterogeneity statistic I(2) can be biased in small meta-analysesBMC Medical Research Methodology 15:35.https://doi.org/10.1186/s12874-015-0024-z
-
Common 5p15.33 and 6p21.33 variants influence lung cancer riskNature Genetics 40:1407–1409.https://doi.org/10.1038/ng.273
-
METAL: fast and efficient meta-analysis of genomewide association scansBioinformatics (Oxford, England) 26:2190–2191.https://doi.org/10.1093/bioinformatics/btq340
Article and author information
Author details
Funding
Wellcome Trust (224982/Z/22/Z)
- Fernanda Morales Berstein
Cancer Research UK (C18281/A29019)
- Konstantinos K Tsilidis
- Philip C Haycock
- Richard M Martin
- Caroline L Relton
- George Davey Smith
Hellenic Republic's Operational Programme "Competitiveness, Entrepreneurship & Innovation" (OΠΣ 5047228)
- Konstantinos K Tsilidis
NIHR Biomedical Research Centre at University Hospitals Bristol
- Richard M Martin
Weston NHS Foundation Trust
- Richard M Martin
NIHR Senior Investigator (NIHR202411)
- Richard M Martin
Medical Research Council (MC_UU_00011/5 and MC_UU_00011/1)
- Caroline L Relton
- George Davey Smith
Alzheimer's Society (AS-PG-19b-010)
- Riccardo E Marioni
National Institutes of Health (U01 AG-18-018)
- Riccardo E Marioni
- Steve Horvath
de Pass Vice Chancellor's Research Fellow at the University of Bristol
- Rebecca C Richmond
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Richard Wilkinson for proofreading several versions of the manuscript and Dr Jean Morrison for helping us interpret the CAUSE analysis output. We acknowledge the participants and investigators of the FinnGen and UK Biobank studies. GWAS data on parental history of cancer were generated using the UK Biobank Resource under application number 15825. Finally, we thank the BCAC, OCAC, PRACTICAL, ILCCO, GECCO and CIMBA consortiums for their contributions.
CRUK and PRACTICAL consortium Funding Acknowledgements: This work was supported by the Canadian Institutes of Health Research, European Commission's Seventh Framework Programme grant agreement n° 223175 (HEALTH-F2-2009-223175), Cancer Research UK Grants C5047/A7357, C1287/A10118, C1287/A16563, C5047/A3354, C5047/A10692, C16913/A6135, and The National Institute of Health (NIH) Cancer Post-Cancer GWAS initiative grant: No. 1 U19 CA 148537-01 (the GAME-ON initiative).
We also thank the following for funding support: The Institute of Cancer Research and The Everyman Campaign, The Prostate Cancer Research Foundation, Prostate Research Campaign UK (now PCUK), The Orchid Cancer Appeal, Rosetrees Trust, The National Cancer Research Network UK, The National Cancer Research Institute (NCRI) UK. We are grateful for support of NIHR funding to the NIHR Biomedical Research Centre at The Institute of Cancer Research, The Royal Marsden NHS Foundation Trust, and Manchester NIHR Biomedical Research Centre. The Prostate Cancer Program of Cancer Council Victoria also acknowledge grant support from The National Health and Medical Research Council, Australia (126402, 209057, 251533, , 396414, 450104, 504700, 504702, 504715, 623204, 940394, 614296,), VicHealth, Cancer Council Victoria, The Prostate Cancer Foundation of Australia, The Whitten Foundation, PricewaterhouseCoopers, and Tattersall’s. EAO, DMK, and EMK acknowledge the Intramural Program of the National Human Genome Research Institute for their support.
Genotyping of the OncoArray was funded by the US National Institutes of Health (NIH) [U19 CA 148537 for ELucidating Loci Involved in Prostate cancer SuscEptibility (ELLIPSE) project and X01HG007492 to the Center for Inherited Disease Research (CIDR) under contract number HHSN268201200008I]. Additional analytic support was provided by NIH NCI U01 CA188392 (PI: Schumacher).
Research reported in this publication also received support from the National Cancer Institute of the National Institutes of Health under Award Numbers U10 CA37429 (CD Blanke), and UM1 CA182883 (CM Tangen/IM Thompson). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Funding for the iCOGS infrastructure came from: the European Community's Seventh Framework Programme under grant agreement n° 223175 (HEALTH-F2-2009-223175) (COGS), Cancer Research UK (C1287/A10118, C1287/A 10710, C12292/A11174, C1281/A12014, C5047/A8384, C5047/A15007, C5047/A10692, C8197/A16565), the National Institutes of Health (CA128978) and Post-Cancer GWAS initiative (1U19 CA148537, 1U19 CA148065 and 1U19 CA148112 - the GAME-ON initiative), the Department of Defence (W81XWH-10-1-0341), the Canadian Institutes of Health Research (CIHR) for the CIHR Team in Familial Risks of Breast Cancer, Komen Foundation for the Cure, the Breast Cancer Research Foundation, and the Ovarian Cancer Research Fund.
Ethics
Human subjects: This research did not require ethical approval as it used secondary, genome-wide association data from studies that obtained informed consent from all participants and ethical approval from review boards and/or ethics committees.
Copyright
© 2022, Morales Berstein 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,179
- views
-
- 720
- downloads
-
- 46
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
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)
Further reading
-
- Genetics and Genomics
Over the last three decades, human genetics has gone from dissecting high-penetrance Mendelian diseases to discovering the vast and complex genetic etiology of common human diseases. In tackling this complexity, scientists have discovered the importance of numerous genetic processes – most notably functional regulatory elements – in the development and progression of these diseases. Simultaneously, scientists have increasingly used multiplex assays of variant effect to systematically phenotype the cellular consequences of millions of genetic variants. In this article, we argue that the context of genetic variants – at all scales, from other genetic variants and gene regulation to cell biology to organismal environment – are critical components of how we can employ genomics to interpret these variants, and ultimately treat these diseases. We describe approaches to extend existing experimental assays and computational approaches to examine and quantify the importance of this context, including through causal analytic approaches. Having a unified understanding of the molecular, physiological, and environmental processes governing the interpretation of genetic variants is sorely needed for the field, and this perspective argues for feasible approaches by which the combined interpretation of cellular, animal, and epidemiological data can yield that knowledge.
-
- Developmental Biology
- Genetics and Genomics
Paternal obesity has been implicated in adult-onset metabolic disease in offspring. However, the molecular mechanisms driving these paternal effects and the developmental processes involved remain poorly understood. One underexplored possibility is the role of paternally-induced effects on placenta development and function. To address this, we investigated paternal high-fat diet-induced obesity in relation to sperm histone H3 lysine 4 tri-methylation signatures, the placenta transcriptome and cellular composition. C57BL6/J male mice were fed either a control or high-fat diet for 10 weeks beginning at 6 weeks of age. Males were timed-mated with control-fed C57BL6/J females to generate pregnancies, followed by collection of sperm, and placentas at embryonic day (E)14.5. Chromatin immunoprecipitation targeting histone H3 lysine 4 tri-methylation (H3K4me3) followed by sequencing (ChIP-seq) was performed on sperm to define obesity-associated changes in enrichment. Paternal obesity corresponded with altered sperm H3K4me3 at promoters of genes involved in metabolism and development. Notably, sperm altered H3K4me3 was also localized at placental enhancers. Bulk RNA-sequencing on placentas revealed paternal obesity-associated sex-specific changes in expression of genes involved in hypoxic processes such as angiogenesis, nutrient transport, and imprinted genes, with a subset of deregulated genes showing changes in H3K4me3 in sperm at corresponding promoters. Paternal obesity was also linked to impaired placenta development; specifically, a deconvolution analysis revealed altered trophoblast cell lineage specification. These findings implicate paternal obesity-effects on placenta development and function as one potential developmental route to offspring metabolic disease.