Hostname: page-component-cd9895bd7-gxg78 Total loading time: 0 Render date: 2024-12-24T16:01:41.820Z Has data issue: false hasContentIssue false

Differences in DNA Methylation-Based Age Prediction Within Twin Pairs Discordant for Cancer

Published online by Cambridge University Press:  08 September 2022

Hannes F. Bode*
Affiliation:
Institute for Molecular Medicine Finland FIMM, HiLIFE, University of Helsinki, Helsinki, Finland
Aino Heikkinen
Affiliation:
Institute for Molecular Medicine Finland FIMM, HiLIFE, University of Helsinki, Helsinki, Finland
Sara Lundgren
Affiliation:
Institute for Molecular Medicine Finland FIMM, HiLIFE, University of Helsinki, Helsinki, Finland
Jaakko Kaprio
Affiliation:
Institute for Molecular Medicine Finland FIMM, HiLIFE, University of Helsinki, Helsinki, Finland
Miina Ollikainen
Affiliation:
Institute for Molecular Medicine Finland FIMM, HiLIFE, University of Helsinki, Helsinki, Finland
*
Author for correspondence: Hannes Bode, Email: hannes.bode@helsinki.fi

Abstract

DNA methylation-based age acceleration (DNAmAA) is associated with cancer, with both cancer tissue and blood showing increased DNAmAA. We aimed to investigate whether DNAmAA is associated with cancer risk within twin pairs discordant for cancer, and whether DNAmAA has the potential to serve as a biomarker for such. The study included 47 monozygotic and 48 same-sex-dizygotic cancer-discordant twin pairs from the Finnish Twin Cohort study with blood samples available between 17 and 31 years after the cancer diagnosis. We studied all cancers (95 pairs), then separately breast cancer (24 pairs) and all sites other than breast cancer (71 pairs). DNAmAA was calculated for seven models: Horvath, Horvath intrinsic epigenetic age acceleration, Hannum, Hannum intrinsic epigenetic age acceleration, Hannum extrinsic epigenetic age acceleration, PhenoAge and GrimAge. Within-pair differences in DNAmAA were analyzed by paired t tests and linear regression. Twin pairs sampled before cancer diagnosis did not differ significantly in DNAmAA. However, the within-pair differences in DNAmAA before cancer diagnosis increased significantly the closer the cancer diagnosis was, and this acceleration extended for years after the diagnosis. Pairs sampled after the diagnosis differed for DNAmAA with the Horvath models capturing cancer diagnosis-associated DNAmAA across all three cancer groupings. The results suggest that DNAmAA in blood is associated with cancer diagnosis. This may be due to epigenetic alterations in relation to cancer, its treatment or associated lifestyle changes. Based on the current study, the biomarker potential of DNAmAA in blood appears to be limited.

Type
Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2022. Published by Cambridge University Press on behalf of International Society for Twin Studies

Age and aging are two of the biggest risk factors for many diseases, including cancer (Aunan et al., Reference Aunan, Cho and Søreide2017; López-Otín et al., Reference López-Otín, Blasco, Partridge, Serrano and Kroemer2013). However, aging can be approached in two ways: chronological aging and biological aging. Chronological age is defined as the calendar time between birth and the time of observation. Biological age can be defined as the physiological state of a person in terms of cellular and organ functions (Horvath & Raj, Reference Horvath and Raj2018; López-Otín et al., Reference López-Otín, Blasco, Partridge, Serrano and Kroemer2013). Biological aging rates may differ from chronological aging rates and are often associated with health-related factors (Ahadi et al., Reference Ahadi, Zhou, Schüssler-Fiorenza Rose, Sailani, Contrepois, Avina, Ashland, Brunet and Snyder2020; Horvath & Raj, Reference Horvath and Raj2018). From this perspective, the difference between chronological and biological age may serve as a better biomarker for individual cancer risk than chronological age alone.

For the study of biological aging, DNA methylation (DNAm, addition of a methyl group to a cytosine base next to guanine) provides a valuable source of information. DNAm-based age prediction measures age-related changes in DNAm and can serve as a surrogate for biological age (Horvath & Raj, Reference Horvath and Raj2018). Such prediction tools enable a straightforward study of biological aging of each individual in large cohorts (Horvath & Raj, Reference Horvath and Raj2018). Increased DNAm-based age has been shown to be associated, for example, with breast cancer in the cancer tissue (Horvath, Reference Horvath2013) and in the adjacent healthy tissue (Hofstatter et al., Reference Hofstatter, Horvath, Dalela, Gupta, Chagpar, Wali, Bossuyt, Storniolo, Hatzis, Patwardhan, Von Wahlde, Butler, Epstein, Stavris, Sturrock, Au, Kwei and Pusztai2018) and in the blood of breast cancer patients (Kresovich et al., Reference Kresovich, Xu, O’Brien, Weinberg, Sandler and Taylor2019). Further, biological age captured by DNAm is associated with different risk factors for cancer (M. Chen et al., Reference Chen, Wong, Nguyen, Dite, Stone, Dugué, Giles, Southey, Milne, Hopper and Li2019; Xiao, Miller et al., Reference Xiao, Miller, Peng, Levine, Conneely, Zhao, Eldridge, Wommack, Jeon, Higgins, Shin, Saba, Smith, Burtness, Park, Irwin, Ferrucci, Ulrich, Qian and Bruner2021), suggesting that DNAm-based age could summarize exposure to such risk factors. Higher exposure would lead to higher cancer risk and increased DNAm-based age; thus, DNAm-based age could serve as a biomarker for cancer risk.

To calculate DNAm-based age, the DNAm is measured in each sample using Illumina BeadChips (Illumina, San Diego, CA, USA). Using mathematical models, individual age estimates (termed DNAm age) are predicted based on the methylation status of a predefined subset of DNAm sites (CpG sites) whose methylation correlate with age (Horvath & Raj, Reference Horvath and Raj2018; Yu et al., Reference Yu, Hazelton, Luebeck and Grady2020). A variant of this is DNAm-based age acceleration (DNAmAA), which is calculated by regressing the DNAm age on chronological age. The DNAmAA values are the residuals of this regression. Consequently, the DNAmAA is independent of chronological age and indicates a difference between the chronological age and the measured DNAm age (i.e., biological aging of an individual). Using DNAmAA instead of DNAm age allows for comparison of individuals in a cohort sampled at different chronological ages (B. H. Chen et al., Reference Chen, Marioni, Colicino, Peters, Ward-Caviness, Tsai, Roetker, Just, Demerath, Guan, Bressler, Fornage, Studenski, Vandiver, Moore, Tanaka, Kiel, Liang, Vokonas and Horvath2016; Horvath & Raj, Reference Horvath and Raj2018).

Two of the first and frequently used DNAmAA models are the Horvath model (Horvath, Reference Horvath2013), also referred to as the ‘epigenetic clock’ or the ‘multi-tissue DNAm age estimator’, and the Hannum model (Hannum et al., Reference Hannum, Guinney, Zhao, Zhang, Hughes, Sadda, Klotzle, Bibikova, Fan, Gao, Deconde, Chen, Rajapakse, Friend, Ideker and Zhang2013). Both of these models were developed by selecting CpG sites that correlate with chronological age (Hannum et al., Reference Hannum, Guinney, Zhao, Zhang, Hughes, Sadda, Klotzle, Bibikova, Fan, Gao, Deconde, Chen, Rajapakse, Friend, Ideker and Zhang2013; Horvath, Reference Horvath2013). DNAmAA measured by both of these tools are associated with cancer (Horvath, Reference Horvath2013; Kresovich et al., Reference Kresovich, Xu, O’Brien, Weinberg, Sandler and Taylor2019). Later, to enable better prediction of morbidity and mortality, the second-generation DNAmAA models were developed, such as the PhenoAge models (Levine et al., Reference Levine, Lu, Quach, Chen, Assimes, Bandinelli, Hou, Baccarelli, Stewart, Li, Whitsel, Wilson, Reiner, Aviv, Lohman, Liu, Ferrucci and Horvath2018) and GrimAge (Lu et al., Reference Lu, Quach, Wilson, Reiner, Aviv, Raj, Hou, Baccarelli, Li, Stewart, Whitsel, Assimes, Ferrucci and Horvath2019). In addition to chronological age, these two models were built by incorporating blood-derived biomarkers of aging (PhenoAge and GrimAge) and time to death (GrimAge) as a reference. Association between DNAmAA and cancer has been shown by the PhenoAge (Levine et al., Reference Levine, Lu, Quach, Chen, Assimes, Bandinelli, Hou, Baccarelli, Stewart, Li, Whitsel, Wilson, Reiner, Aviv, Lohman, Liu, Ferrucci and Horvath2018) and the GrimAge model (Lu et al., Reference Lu, Quach, Wilson, Reiner, Aviv, Raj, Hou, Baccarelli, Li, Stewart, Whitsel, Assimes, Ferrucci and Horvath2019).

Twin pairs discordant for a disease (here, pairs with a cancer diagnosis in one twin and no cancer diagnosis in the other twin during the study period) provide a powerful study design, since the twins are matched for many known factors (such as age, sex, genetic relatedness and family background) and potentially for many unknown factors. These factors shared by the co-twins of each pair may affect DNAm and DNAmAA (Czyz et al., Reference Czyz, Morahan, Ebers and Ramagopalan2012). Hence, for cancer-discordant twin pairs, within-pair differences in DNAmAA can be considered cancer diagnosis-related, and confounding factors that affect DNAmAA but not cancer risk are largely mitigated by matching the twins. It is known that DNAmAA is correlated between twins, with a considerable amount of heritability (Jylhävä et al., Reference Jylhävä, Hjelmborg, Soerensen, Munoz, Tan, Kuja-Halkola, Mengel-From, Christensen, Christiansen, Hägg, Pedersen and Reynolds2019; Reynolds et al., Reference Reynolds, Tan, Munoz, Jylhävä, Hjelmborg, Christiansen, Hägg and Pedersen2020). However, this correlation decreases over time, and the twins diverge from each other in their DNAmAA (Li et al., Reference Li, Nguyen, Wong, Dugué, Dite, Armstrong, Craig, Mather, Sachdev, Saffery, Sung, Tan, Thalamuthu, Milne, Giles, Southey and Hopper2020). To date, it is not known whether twins in a pair discordant for cancer differ by their DNAmAA, and whether these differences are associated with cancer diagnosis. Knowing this could be important in understanding whether DNAmAA serves as an appropriate biomarker for individual cancer risk. In addition, this would allow us to determine the extent to which a cancer diagnosis affects biological aging post diagnosis, reflected by DNAmAA.

To this end, we identified twin pairs discordant for incident cancer from the Finnish Twin Cohort. After blood-based DNAm measurement, DNAmAA was calculated. Within-pair differences of DNAmAA were first computed to determine whether cancer diagnosis is associated with increased DNAmAA. These analyses were performed separately for pairs where DNA was sampled prior to cancer diagnosis, and for those where DNA was sampled after a diagnosis. This division allowed for investigation of DNAmAA, both independent of and associated with diagnosis and treatment as possible factors affecting DNAmAA. Furthermore, these within-pair differences in DNAmAA were analyzed for associations between the time interval from DNA sampling to cancer diagnosis. All the analyses were performed for all cancers together and then separately for breast cancer and nonbreast cancer.

Materials and Methods

Twin Cohort

The data used in this study originate from the older Finnish Twin Cohort, including same-sex twin pairs born before 1958, and followed up with lifestyle, health and behavioral surveys in four waves (years 1975, 1981, 1990 and 1999–2017). Cancer incidence and cause of death were updated from the Cancer Registry and Statistics Finland, with the latest update in 2018 (Kaprio et al., Reference Kaprio, Bollepalli, Buchwald, Iso-Markku, Korhonen, Kovanen, Kujala, Laakkonen, Latvala, Leskinen, Lindgren, Ollikainen, Piirtola, Rantanen, Rinne, Rose, Sillanpää, Silventoinen, Sipilä and Waller2019). The extracted cancer diagnoses of the cohort were provided in the form of ‘International Classification of Diseases 10th revision’ (ICD10) codes (Harris et al., Reference Harris, Hjelmborg, Adami, Czene, Mucci and Kaprio2019; Skytthe et al., Reference Skytthe, Harris, Czene, Mucci, Adami, Christensen, Hjelmborg, Holm, Nilsen, Kaprio and Pukkala2019). The ICD10 codes were studied using the ICD10 code browser provided by the WHO (World Health Organization, 2019) to identify the different sites of diagnosis (Supplementary Table S1 and Table 1). Epigenetic analyses of DNA samples in the older cohort have been undertaken for a number of specific projects, which provided the DNAm data (Kaprio et al., Reference Kaprio, Bollepalli, Buchwald, Iso-Markku, Korhonen, Kovanen, Kujala, Laakkonen, Latvala, Leskinen, Lindgren, Ollikainen, Piirtola, Rantanen, Rinne, Rose, Sillanpää, Silventoinen, Sipilä and Waller2019). A total of 95 twin pairs discordant for any cancer with available blood DNAm data were included in the current study. The methylation data of the co-twins in each pair is from the same time point, and all twin pairs are same-sex pairs, thereby age and sex are accounted for in the within-pair analyses. For all individuals, smoking status at the time of DNA collection was coded categorically: 77.9% (n = 148) were nonsmokers, 14.7% (n = 28) were former smokers and 7.4% (n = 14) were current smokers, with 21 pairs discordant for smoking behavior (Table 2).

Table 1. Number of the Finnish Twin Cohort twin pairs discordant for cancer shown for the different cancer groups

Note:

a other cancers than breast cancer.

Table 2. Description of the cancer-discordant Finnish Twin Cohort twin pairs

Note:

a ages and time to diagnosis are in years;

b per individual, at the time point of sampling.

DNAm Data

High-molecular-weight blood DNA was bisulphite-converted using the EZ-96 DNA Methylation-Gold Kit (Zymo Research, Irvine, CA, USA) according to the manufacturer’s instructions and hybridized on the Infinium HumanMethylation450 (‘450k’) or EPIC BeadChip (‘EPIC’) (Illumina, San Diego, CA, USA; 21 and 74 twin pairs, respectively) for DNAm profiling following the standard protocols. Samples of both twins in a pair were always processed at the same time.

Sample quality control (QC) was performed using the R package MethylAid (van Iterson et al., Reference van Iterson, Tobi, Slieker, den Hollander, Luijk, Slagboom and Heijmans2014) and minfi (Aryee et al., Reference Aryee, Jaffe, Corrada-Bravo, Ladd-Acosta, Feinberg, Hansen and Irizarry2014) to assess quality based on five control probe metrics with default thresholds, as follows. A sample passed the QC when the log2 intensity of converted red/green channel bisulphite Type I control was higher than 12.75, median methylated versus unmethylated log2 intensity was higher than 10.5, log2 intensity of sample-dependent control probes was higher than 11.75, log2 intensity of sample-independent control probes was higher than 13.25 and when less than 5% of the control probes in a sample showed only background signal. All samples passed the sample QC. Finally, data generated on the 450K and EPIC platforms were combined using minfi’s combineArrays function.

Probe QC was performed using the R packages minfi (Fortin et al., Reference Fortin, Triche and Hansen2017) and wateRmelon (Pidsley et al., Reference Pidsley, Wong, Volta, Lunnon, Mill and Schalkwyk2013). Probes with a bead count < 3, detection p value > 0.01 or intensity value of exactly zero were removed, and probes with a call rate of ≥ 95% across all samples were retained. In total, 1058 probes were set to missing. Afterward, ambiguously mapping and poor-quality probes according to Zhou and colleagues (negative selection criteria: single-nucleotide polymorphisms with minor allele frequency > 1% near target site, color-channel-switching single-nucleotide polymorphisms, off-target hybridization events and incorrect mapping on the GRCh38 genome; Zhou et al., Reference Zhou, Laird and Shen2017) were removed. Additionally, probes binding to sex chromosomes were set to missing.

Functional normalization was applied using ssnoob from the R package minfi (Fortin et al., Reference Fortin, Triche and Hansen2017) with default settings to remove unwanted between-array variability (Zhou et al., Reference Zhou, Laird and Shen2017). Afterward, beta mixture quantile normalization was applied to the normalized data using the R package wateRmelon (Pidsley et al., Reference Pidsley, Wong, Volta, Lunnon, Mill and Schalkwyk2013) to further correct the data for probe design bias.

DNAm-Based Age Prediction

The DNAm-based age was calculated on the preprocessed DNAm data using the online calculator (http://dnamage.genetics.ucla.edu/). DNAm-based age prediction was obtained for seven models: Horvath (Horvath, Reference Horvath2013) and Horvath intrinsic epigenetic age acceleration (Horvath IEAA; B. H. Chen et al., Reference Chen, Marioni, Colicino, Peters, Ward-Caviness, Tsai, Roetker, Just, Demerath, Guan, Bressler, Fornage, Studenski, Vandiver, Moore, Tanaka, Kiel, Liang, Vokonas and Horvath2016), Hannum (Hannum et al., Reference Hannum, Guinney, Zhao, Zhang, Hughes, Sadda, Klotzle, Bibikova, Fan, Gao, Deconde, Chen, Rajapakse, Friend, Ideker and Zhang2013), Hannum IEAA (B. H. Chen et al., Reference Chen, Marioni, Colicino, Peters, Ward-Caviness, Tsai, Roetker, Just, Demerath, Guan, Bressler, Fornage, Studenski, Vandiver, Moore, Tanaka, Kiel, Liang, Vokonas and Horvath2016), Hannum extrinsic epigenetic age acceleration (Hannum EEAA; B. H. Chen et al., Reference Chen, Marioni, Colicino, Peters, Ward-Caviness, Tsai, Roetker, Just, Demerath, Guan, Bressler, Fornage, Studenski, Vandiver, Moore, Tanaka, Kiel, Liang, Vokonas and Horvath2016), PhenoAge (Levine et al., Reference Levine, Lu, Quach, Chen, Assimes, Bandinelli, Hou, Baccarelli, Stewart, Li, Whitsel, Wilson, Reiner, Aviv, Lohman, Liu, Ferrucci and Horvath2018) and GrimAge (Lu et al., Reference Lu, Quach, Wilson, Reiner, Aviv, Raj, Hou, Baccarelli, Li, Stewart, Whitsel, Assimes, Ferrucci and Horvath2019). DNAmAA was calculated by linearly regressing DNAm age over chronological age, with the residuals serving as the DNAmAA measure that was used in further analyses. Detailed information on the epigenetic age acceleration models is presented in the supplements (Supplementary materials and methods on DNAm-based age prediction and Supplementary Table S2 available at the Cambridge Core website (https://www.cambridge.org/core/).

Statistical Analysis

A paired t test was used to examine within-pair differences in DNAmAA. This was done separately for the pairs sampled before and after diagnosis. The association between the time to diagnosis (the time interval between the blood sampling for the DNA and diagnosis) and the within-pair difference in DNAmAA (dDNAmAA) was examined using linear regression (1). A negative time to diagnosis indicates that the sample was collected prior to cancer diagnosis, and a positive time to diagnosis that the sample was collected after diagnosis. A negative dDNAmAA indicates that the cancer-diagnosed twin in a pair is predicted to be younger, and a positive dDNAmAA that the cancer-diagnosed twin in a pair is predicted to be older than the healthy co-twin. A categorical correction term was included in the regression model to correct for confounding by differences in smoking behavior between twins in a pair. This correction term describes the smoking status of the pair by combining the smoking status of each twin in a pair into a single term (‘nonsmoker/nonsmoker’, ‘former smoker/nonsmoker’, ‘current smoker/nonsmoker’, etc.). F-statistics’ p values and adjusted R2 (adR2) values were reported for the regression models.

  1. (1) dDNAmAA (Pair X) ∼ Intercept Time Point of Diagnosis + Slope * Time to Diagnosis (Pair X) + Smoking (Pair X)

A post hoc power analysis was performed for the paired t test and the linear regression using the R package pwr with a significance threshold of p ≤ .05 (Champely, Reference Champely2020; Selya et al., Reference Selya, Rose, Dierker, Hedeker and Mermelstein2012). Only significant results with sufficient power (≥0.8) were considered meaningful, and the conclusions of this study were drawn from such results exclusively. Effects and measures of DNAmAA and dDNAmAA are reported in years.

Results

Cohort Description

A total of 95 same-sex twin pairs discordant for any cancer were identified (Table 1). Of these, 47 were monozygotic (MZ) and 48 were same-sex dizygotic (ssDZ) pairs. Within this cohort, 24 pairs were discordant for breast cancer and 71 pairs were discordant for nonbreast cancer (Table 1 and Supplementary Table S1). Nine individuals had more than one diagnosis, and in such cases, the first diagnosis was considered as the primary site. Age at first cancer diagnosis ranged from 33 to 87 years. The blood samples were collected at least 6 months before the cancer diagnosis in 46 pairs and at least 6 months after the diagnosis in the remaining 49 pairs. The interval between blood collection and the time of first cancer diagnosis ranged from 17 years before the diagnosis to 31 years after the diagnosis (Table 2).

Pan-Cancer

We first examined whether cancer diagnosis is associated with DNAmAA in a co-twin control design. To this end, we analyzed within-pair differences in DNAmAA using all 47 MZ and 48 ssDZ pairs discordant for any cancer as a pan-cancer analysis (Table 3 and Supplementary Figure S1). The twins with cancer from pairs whose blood samples were collected before diagnosis did not differ from their healthy co-twins for DNAmAA, the point estimates being small, and both negative and positive. However, among the twin pairs whose blood samples were collected after the diagnosis, all the DNAmAA models showed between 1.12 and 2.89 years within-pair difference in DNAmAA. The Horvath model and Hannum EEAA model showed slightly larger within-pair differences in DNAmAA after the cancer diagnosis compared with the other models (Table 3 and Supplementary Figure S1).

Table 3. Within-pair differences in DNAmAA in cancer-discordant twin pairs (N = 95). The presented units are in years

Note:

a t-test p value;

b F-statistics p value;

c adjusted coefficient of determination (R 2), proportion of variance in the data explained by the linear model.

The within-pair differences measured by Horvath (2.48 years) and Horvath IEAA (2.18 years) predictions were highly similar, whereas the differences by the Hannum model (1.98 years) deviated from both the Hannum IEAA (1.12 years) and Hannum EEAA (2.89 years) models. These results may suggest that within-pair variations in predicted blood cell-type proportions did not have a large effect on the DNAmAA predictions of the Horvath model, whereas the Hannum models were affected by differences in blood cell-type proportions within twin pairs.

Next, we examined whether differences in DNAmAA within twin pairs depended on the time to cancer diagnosis by regressing within-pair differences in DNAmAA on the time to cancer diagnosis. According to all DNAmAA models investigated, statistically significant or not, in general, the closer to cancer diagnosis a twin pair was sampled and extending for years after the diagnosis, the higher the difference in DNAmAA within the twin pair was. This suggests that twins in a pair diverged in their DNAmAA toward the cancer diagnosis and afterward, where twins with a cancer diagnosis had higher DNAmAA. Rates for this divergence varied from 0.05 years to 0.17 years (dDNAmAA) per calendar year in this cohort. However, only the Horvath and Horvath IEAA models had statistically significant (p ≤ .05) p values, with estimates of about two months (0.16 years and 0.14 years, dDNAmAA) per calendar year. The Hannum, Hannum EEAA and GrimAge models had low p values as well (p ≤ .11), with estimates of one month (0.05 years to 0.11 years, dDNAmAA) per calendar year. Together, this could mean that the divergence of DNAmAA between the twins in a pair is captured differently by the different models. Further, for all DNAmAA models, a rather low fit of the models suggests that there were other factors besides cancer diagnosis that contribute to the observed within-pair differences in DNAmAA (Table 3 and Supplementary Figure S2).

Breast Cancer

As a next step, we analyzed separately, as the largest subgroup, the 24 twin pairs discordant for breast cancer. For the 10 twin pairs with their blood samples collected before breast cancer diagnosis, only the GrimAge model showed within-pair difference, where the twin with cancer is 1.22 years younger than their healthy co-twin (p = .053, Table 4). None of the other DNAmAA models showed that the twins in a pair differed for their DNAmAA. For the 14 twin pairs sampled after breast cancer diagnosis, the DNAmAA models predicted the cancer-diagnosed twin to be on average 1.88 years to 5.83 years older than their healthy co-twin (Table 4 and Supplementary Figure S3). The PhenoAge model for breast cancer resulted in the highest within-pair difference in DNAmAA (5.83 years, 95% CI [0.61 years, 11.04 years]), with the difference being much higher than in the pan-cancer analysis (1.66 years, 95% CI [0.69 years, 4.01 years). Thus, the morbidity captured by the PhenoAge model might be particularly relevant for breast cancer (Table 4 and Supplementary Figure S3).

Table 4. Within-pair differences in DNAmAA in breast cancer-discordant twin pairs (N = 24). The presented units are in years

Note:

a t-test p value;

b F-statistics p value;

c adjusted coefficient of determination (R 2), proportion of variance in the data explained by the linear model.

As in the pan-cancer analysis, the breast cancer analysis also showed similar effects with respect to the Horvath and Hannum models and their adjustments for predicted blood cell-type proportions. The predictions by the Horvath models were highly similar in their effect sizes, while the predictions by the Hannum models were affected by the blood cell-type proportion adjustments (Table 4).

In the regression of within-pair differences in DNAmAA over time to diagnosis, the Horvath, Horvath IEAA and GrimAge models showed a divergence of the twins with the cancer-diagnosed twin on average getting older faster than their healthy co-twin (p ≤ .05). For these, the rate for this within-pair divergence in DNAmAA in breast cancer was two months (0.14 years to 0.17 years, dDNAmAA) per calendar year toward the diagnosis and afterward. However, the low fits of the regression models indicated the presence of additional unknown factors contributing to the within-pair difference in DNAmAA (Table 4 and Supplementary Figure S4).

Nonbreast Cancers

In the third analysis, we separately analyzed the 71 twin pairs from the cohort discordant for cancers other than breast cancer (nonbreast cancers). Here again, the 36 twin pairs sampled before diagnosis showed no within-pair difference in DNAmAA. For the 35 twin pairs sampled after diagnosis, only the Horvath, Horvath IEAA and GrimAge models showed that the twin diagnosed with cancer is 0.97 years to 2.15 years older than their healthy co-twin. The other models yielded point estimates of the same magnitude with the exception of PhenoAge (−0.01 years; Table 5 and Supplementary Figure S5). However, the observed mean differences were smaller compared to the breast cancer analysis (Table 4). In the case of PhenoAge, the lack of any difference (−0.01 years, 95% CI [−2.51 years, 2.49 years]) contrasts with the large effect size in breast cancer (5.83 years, 95% CI [0.61 years, 11.04 years]).

Table 5. Within-pair differences in DNAmAA in nonbreast cancer-discordant twin pairs (N = 71). The presented units are in years

Note:

a t-test p value;

b F-statistics p value;

c adjusted coefficient of determination (R 2), proportion of variance in the data explained by the linear model.

For the nonbreast cancer analysis, when regressing the within-pair differences in DNAmAA over time to diagnosis, only the Horvath and Horvath IEAA models showed substantive effects between the time to diagnosis and the within-pair difference in DNAmAA. This means that the co-twins diverged over time across the cancer types summarized here, where the twin with the cancer diagnosis aged faster. The rate of within-pair divergence of DNAmAA measured by the Horvath and Horvath IEAA models was two months (0.14 years and 0.13 years, dDNAmAA) per calendar year. While a good fit of the regression models was also not observed in the breast cancer group (Table 4), the fit was even lower in the nonbreast cancer group of discordant pairs, indicating an even greater impact of unknown factors on DNAmAA (Table 5 and Supplementary Figure S6).

Discussion

To our knowledge, this is the first study to report within-pair differences in DNAmAA in twin pairs discordant for cancer. The observed within-pair differences increased the closer the sample collection was to the time of cancer diagnosis and afterward, with the cancer-diagnosed twin having aged faster. However, significant and substantial within-pair differences in DNAmAA were observed mainly in twin pairs sampled after cancer diagnosis. Additionally, the observed within-pair differences for breast cancer-discordant twin pairs were notably high.

DNAmAA estimated by GrimAge was lower in the twins diagnosed with breast cancer compared with their healthy co-twins, while for the other DNAmAA estimates and cancer groups (pan-cancer and nonbreast cancer), the twins of pairs sampled before the cancer diagnosis did not differ by their DNAmAA. The within-pair differences, however, increased the closer to cancer diagnosis the pair was sampled, which may relate to the subclinical development of cancer prior to diagnosis. Nevertheless, the current study provided only minimal evidence for substantial within-pair differences in DNAmAA prior to cancer diagnosis. This observation is consistent with previous prospective studies on unrelated individuals with multiple cancer types, including breast cancer (Ambatipudi et al., Reference Ambatipudi, Horvath, Perrier, Cuenin, Hernandez-Vargas, Le Calvez-Kelm, Durand, Byrnes, Ferrari, Bouaoun, Sklias, Chajes, Overvad, Severi, Baglietto, Clavel-Chapelon, Kaaks, Barrdahl, Boeing and Herceg2017; Durso et al., Reference Durso, Bacalini, Sala, Pirazzini, Marasco, Bonafé, do Valle, Gentilini, Castellani, Faria, Franceschi, Garagnani and Nardini2017; Levine et al., Reference Levine, Hosgood, Chen, Absher, Assimes and Horvath2015; Zheng et al., Reference Zheng, Joyce, Colicino, Liu, Zhang, Dai, Shrubsole, Kibbe, Gao, Zhang, Jafari, Vokonas, Schwartz, Baccarelli, Baccarelli and Hou2016), with comparable follow-up time to our cohort (≤20 years vs. ≤17 years) and with varying cohort sizes (43–451 cases and 424–2029 unrelated controls). Our findings of nonsubstantially elevated DNAmAA before diagnosis argue against the suitability of DNAmAA alone as a biomarker of cancer risk. Cancer risk predictors that are built directly on DNAm seem to be more promising cancer biomarkers (Terry et al., Reference Terry, McDonald, Wu, Eng and Santella2016). However, when the PhenoAge DNAmAA estimate was integrated into a model along with four other DNAm estimators (including estimates for blood cell types and age acceleration) and 19 breast cancer risk-associated CpGs, the resulting breast cancer risk score, mBCRS, was able to predict breast cancer risk with similar strength to polygenic risk scores (Kresovich et al., Reference Kresovich, Xu, O’Brien, Shi, Weinberg, Sandler and Taylor2022). The mBCRS combined with the genetic risk data may enhance breast cancer prediction (Kresovich et al., Reference Kresovich, Xu, O’Brien, Shi, Weinberg, Sandler and Taylor2022).

In general, larger within-pair differences in DNAmAA in twin pairs discordant for any cancer were observed in the pairs sampled after diagnosis compared with the pairs sampled before the cancer diagnosis. The observed elevated DNAmAA in the blood of the cancer-diagnosed twins sampled after the diagnosis may be due to the cancer or its treatment. In line with this, recent studies on breast (Sehl et al., Reference Sehl, Carroll, Horvath and Bower2020) and head and neck cancers (Xiao, Beitler et al., Reference Xiao, Beitler, Peng, Levine, Conneely, Zhao, Felger, Wommack, Chico, Jeon, Higgins, Shin, Saba, Burtness, Bruner and Miller2021) have shown that DNAmAA in blood increases significantly after cancer treatment. It has been further shown that the survivors of esophageal cancer have blood DNAmAA close to zero, whereas those who died within the 3-year follow-up had highly elevated DNAmAA (Beynon et al., Reference Beynon, Ingle, Langdon, May, Ness, Martin, Suderman, Ingarfield, Marioni, McCartney, Waterboer, Pawlita, Relton, Smith and Richmond2022). Additionally, adult survivors of childhood cancers have higher DNAmAA than adults with no childhood cancer diagnosis (Qin et al., Reference Qin, Li, Song, Wilson, Easton, Mulder, Plyler, Neale, Walker, Zhou, Pan, Hudson, Yasui, Robison, Zhang, Ness and Wang2021). The increased DNAmAA after cancer diagnosis could be a response to treatment-related mutational stress, as many cancer treatments have mutagenic effects (Venkatesan et al., Reference Venkatesan, Swanton, Taylor and Costello2017). It has been hypothesized that epigenetic maintenance programs are hyperactivated to ensure genome stability as a reaction to mutational stress; thereby, the increased DNAmAA would mark an accumulated record of the mutational stress (Horvath & Raj, Reference Horvath and Raj2018; Levine et al., Reference Levine, Lu, Quach, Chen, Assimes, Bandinelli, Hou, Baccarelli, Stewart, Li, Whitsel, Wilson, Reiner, Aviv, Lohman, Liu, Ferrucci and Horvath2018). This could have led to the increased within-pair differences we observed after the cancer diagnosis. In our study, increased DNAmAA after diagnosis was associated especially with breast cancer. This may be due to hormonal treatment strategies common in breast cancer therapy (Barzaman et al., Reference Barzaman, Karami, Zarei, Hosseinzadeh, Kazemi, Moradi-Kalbolandi, Safari and Farahmand2020). The association between female sex hormone exposures and DNAmAA is widely described (M. Chen et al., Reference Chen, Wong, Nguyen, Dite, Stone, Dugué, Giles, Southey, Milne, Hopper and Li2019; Horvath, Reference Horvath2013; Levine et al., Reference Levine, Lu, Chen, Hernandez, Singleton, Ferrucci, Bandinelli, Assimes and Horvath2016; Sehl et al., Reference Sehl, Henry, Storniolo, Horvath and Ganz2021). Therefore, hormone-targeted treatment of breast cancer may have resulted in the observed high within-pair differences in DNAmAA in the breast cancer-discordant pairs postdiagnosis.

Another cause of increased DNAmAA could be trauma and posttraumatic stress related to cancer diagnosis. Several studies have shown that trauma (Boks et al., Reference Boks, Mierlo, Rutten, Radstake, De Witte, Geuze, Horvath, Schalkwyk, Vinkers, Broen and Vermetten2015; Wolf, Maniates et al., Reference Wolf, Maniates, Nugent, Maihofer, Armstrong, Ratanatharathorn, Ashley-Koch, Garrett, Kimbrel, Lori, Aiello, Baker, Beckham, Boks, Galea, Geuze, Hauser, Kessler and Logue2018) and posttraumatic stress disorder (Wang et al., Reference Wang, Hui, Goldberg, Smith, Kaseer, Murrah, Levantsevych, Shallenberger, Diggers, Bremner, Vaccarino and Sun2022; Wolf, Logue et al., Reference Wolf, Logue, Stoop, Schichman, Stone, Sadeh, Hayes and Miller2018; Wolf, Maniates et al., Reference Wolf, Maniates, Nugent, Maihofer, Armstrong, Ratanatharathorn, Ashley-Koch, Garrett, Kimbrel, Lori, Aiello, Baker, Beckham, Boks, Galea, Geuze, Hauser, Kessler and Logue2018; Yang et al., Reference Yang, Wu, Verhoeven, Gautam, Reus, Kang, Flory, Abu-Amara, Systems Biology Consortium, Hood, Doyle, Yehuda, Marmar, Jett, Hammamieh, Mellon and Wolkowitz2020) are associated with DNAmAA in blood. Importantly for this study, traumatic life events have been shown to be associated with increased risk of breast cancer (Ginsberg et al., Reference Ginsberg, Price, Ingram and Nottage1996; Lin et al., Reference Lin, Wang, Zhong, Huang, Peng, Shan, Wang and Sun2013; Santos et al., Reference Santos, Horta, do Amaral, Branco, Fernandes, Galvão and Fernandes2009), including a study from the Finnish Twin Cohort (Lillberg et al., Reference Lillberg, Verkasalo, Kaprio, Teppo, Helenius and Koskenvuo2003). A traumatic life event, such as the death of a close person or separation/divorce, significantly increases the risk of being diagnosed with breast cancer after such an event. It has been hypothesized that the traumatic event triggers changes in hormone metabolism (Lillberg et al., Reference Lillberg, Verkasalo, Kaprio, Teppo, Helenius and Koskenvuo2003; Lin et al., Reference Lin, Wang, Zhong, Huang, Peng, Shan, Wang and Sun2013; Santos et al., Reference Santos, Horta, do Amaral, Branco, Fernandes, Galvão and Fernandes2009) or the immune system (Lillberg et al., Reference Lillberg, Verkasalo, Kaprio, Teppo, Helenius and Koskenvuo2003), but a biological mechanism for this observation has not been demonstrated. Similar results have been reported for colorectal cancer (Kune et al., Reference Kune, Kune, Watson and Rahe1991) and for melanoma, lung cancer and hematopoietic malignancies (Levav et al., Reference Levav, Kohn, Iscovich, Abramson, Tsai and Vigdorovich2000). Overall, the traumatic event of a cancer diagnosis might increase DNAmAA in the twin diagnosed with cancer, or the twin with a cancer diagnosis might have been more likely to have been exposed to traumatic events prior to diagnosis, resulting in increased DNAmAA. Overall, the pathways from severe traumatic experiences to cancer are unclear and to which extent they involve biological aging processes picked up by epigenetic clocks are not established and outside the scope of the present work.

While the largest within-pair differences in DNAmAA were observed in breast cancer-discordant pairs, the Horvath and Horvath IEAA models captured significant associations between DNAmAA and time to diagnosis across all three cancer groups (pan-cancer, breast cancer and nonbreast cancer). Notably, the Horvath model was constructed using chronological age as a metric for aging, a common risk factor for most cancers. Furthermore, unlike the other models presented here, the Horvath model was developed for the prediction of DNAmAA in multiple tissues, including applicability to cancer tissues (Horvath, Reference Horvath2013). Both features may result in the Horvath models being able to robustly detect an increase in DNAmAA in the blood of individuals diagnosed with cancer. Our study further showed that the Hannum models were affected by estimated white blood cell percentages, as has been described before (Horvath & Raj, Reference Horvath and Raj2018). Hematologic changes with age are a well-known phenomenon (Groarke & Young, Reference Groarke and Young2019), and moreover, chemotherapy and radiotherapy frequently lead to hematologic changes (Iorio et al., Reference Iorio, Spieler, Ricardi and Dal Pra2021; Spivak et al., Reference Spivak, Gascón and Ludwig2009). This would suggest that cancer diagnosis-associated events, most likely the treatment, induce hematologic aging, which was captured by the Hannum models in the current study.

The main strength of the current study is the discordant twin study design. The fact that twins are highly matched on many factors — such as genetics, age, sex, early environment and lifestyle — that influence DNAm and aging leads to the best matched case–control setting available in humans (Ceribelli & Selmi, Reference Ceribelli and Selmi2020). For the current study, especially the fact that the twins in a pair are of the same age is an important advantage, which is difficult to set up to this degree in studies investigating sibling pairs or unrelated individuals. Although DNAmAA is corrected for chronological age, this correction may not be complete and residual effects of chronological age may still confound the associations. Pairwise comparison of the co-twins sampled at the same time point, and thus at the same age, overcomes this problem.

Our study also has some limitations worth acknowledging. Although our prospective twin cohort covered a period of 48 years around the cancer diagnosis, repeated measurement of DNAm for individual twin pairs was not available. Therefore, we cannot draw any conclusions on individual epigenetic aging trajectories in relation to cancer diagnosis. Such investigations would be highly relevant, since DNAmAA of twins in a pair become less correlated over time (Li et al., Reference Li, Nguyen, Wong, Dugué, Dite, Armstrong, Craig, Mather, Sachdev, Saffery, Sung, Tan, Thalamuthu, Milne, Giles, Southey and Hopper2020), and especially since the DNAmAA rate is not constant over time (Snir et al., Reference Snir, Farrell and Pellegrini2019). Resources for conducting such longitudinal studies in genetically informative samples such as twin pairs are, unfortunately, rare.

Although to the best of our knowledge this is the largest DNAm study on cancer-discordant twin pairs to date, we only had sufficient power to examine breast cancer as a separate cancer group. Therefore, we may have missed associations between DNAmAA and other cancer types. This study served as a pilot study, and our results on the blood DNAmAA in breast cancer warrant further investigations on the value of DNAmAA as a biomarker for breast cancer risk.

In conclusion, DNAmAA as a proxy for biological aging is associated with cancer diagnosis in a manner dependent on cancer type and time to cancer diagnosis. However, only weak associations of DNAmAA with cancer were found before cancer diagnosis. Therefore, DNAmAA on its own may not serve as a reliable biomarker for individual cancer risk; however, larger sample sizes on different cancer types and longitudinal designs should be included in future studies to draw firm conclusions. Importantly, DNAmAA is increased after cancer diagnosis, potentially resulting from the cancer itself, from potential lifestyle changes and from cancer treatments, which makes DNAmAA a valuable tool to study cancer diagnosis-related aging and morbidity.

Supplementary Material

To view supplementary material for this article, please visit https://doi.org/10.1017/thg.2022.32.

Acknowledgments

The authors thank the participants for their invaluable contribution to the study. The technical staff at the Finnish Twin Cohort study are acknowledged for their help in collecting the data. The authors also wish to thank Alyce M Whipp, PhD, from Language Services, University of Helsinki, for the valuable help with language editing.

Financial Support

This research was funded by the European Union’s Horizon 2020 Research and Innovation Programme, Marie Skłodowska-Curie (JK, grant number 859860). This project has received further funding from the Academy of Finland (MO, grant numbers 297908, 328685 & JK, grant 336823) and the Sigrid Juselius Foundation (MO and JK).

Conflicts of Interest

None.

Ethical Standards

The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.

References

Ahadi, S., Zhou, W., Schüssler-Fiorenza Rose, S. M., Sailani, M. R., Contrepois, K., Avina, M., Ashland, M., Brunet, A., & Snyder, M. (2020). Personal aging markers and ageotypes revealed by deep longitudinal profiling. Nature Medicine, 26, 8390.CrossRefGoogle ScholarPubMed
Ambatipudi, S., Horvath, S., Perrier, F., Cuenin, C., Hernandez-Vargas, H., Le Calvez-Kelm, F., Durand, G., Byrnes, G., Ferrari, P., Bouaoun, L., Sklias, A., Chajes, V., Overvad, K., Severi, G., Baglietto, L., Clavel-Chapelon, F., Kaaks, R., Barrdahl, M., Boeing, H., … Herceg, Z. (2017). DNA methylome analysis identifies accelerated epigenetic ageing associated with postmenopausal breast cancer susceptibility. European Journal of Cancer, 75, 299307.CrossRefGoogle ScholarPubMed
Aryee, M. J., Jaffe, A. E., Corrada-Bravo, H., Ladd-Acosta, C., Feinberg, A. P., Hansen, K. D., & Irizarry, R. A. (2014). Minfi: A flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics, 30, 13631369.CrossRefGoogle ScholarPubMed
Aunan, J. R., Cho, W. C., & Søreide, K. (2017). The biology of aging and cancer: A brief overview of shared and divergent molecular hallmarks. Aging and Disease, 8, 628642.CrossRefGoogle ScholarPubMed
Barzaman, K., Karami, J., Zarei, Z., Hosseinzadeh, A., Kazemi, M. H., Moradi-Kalbolandi, S., Safari, E., & Farahmand, L. (2020). Breast cancer: Biology, biomarkers, and treatments. International Immunopharmacology, 84, 106535.CrossRefGoogle ScholarPubMed
Beynon, R. A., Ingle, S. M., Langdon, R., May, M., Ness, A., Martin, R. M., Suderman, M., Ingarfield, K., Marioni, R. E., McCartney, D. L., Waterboer, T., Pawlita, M., Relton, C., Smith, G. D., & Richmond, R. C. (2022). Epigenetic biomarkers of ageing are predictive of mortality risk in a longitudinal clinical cohort of individuals diagnosed with oropharyngeal cancer. Clinical Epigenetics, 14, 113.CrossRefGoogle Scholar
Boks, M. P., Mierlo, H. C. va., Rutten, B. P. F., Radstake, T. R. D. J., De Witte, L., Geuze, E., Horvath, S., Schalkwyk, L. C., Vinkers, C. H., Broen, J. C. A., & Vermetten, E. (2015). Longitudinal changes of telomere length and epigenetic age related to traumatic stress and post-traumatic stress disorder. Psychoneuroendocrinology, 51, 506512.CrossRefGoogle ScholarPubMed
Ceribelli, A., & Selmi, C. (2020). Epigenetic methods and twin studies. Advances in Experimental Medicine and Biology, 1253, 95104.CrossRefGoogle ScholarPubMed
Champely, S. (2020). Basic Functions for Power Analysis [R package pwr version 1.3-0]. Comprehensive R Archive Network (CRAN). https://cran.r-project.org/package=pwr Google Scholar
Chen, B. H., Marioni, R. E., Colicino, E., Peters, M. J., Ward-Caviness, C. K., Tsai, P.-C., Roetker, N. S., Just, A. C., Demerath, E. W., Guan, W., Bressler, J., Fornage, M., Studenski, S., Vandiver, A. R., Moore, A. Z., Tanaka, T., Kiel, D. P., Liang, L., Vokonas, P., … Horvath, S. (2016). DNA methylation-based measures of biological age: meta-analysis predicting time to death. Aging (Albany NY), 8, 1844.CrossRefGoogle ScholarPubMed
Chen, M., Wong, E. M., Nguyen, T. L., Dite, G. S., Stone, J., Dugué, P.-A., Giles, G. G., Southey, M. C., Milne, R. L., Hopper, J. L., & Li, S. (2019). DNA methylation-based biological age, genome-wide average DNA methylation, and conventional breast cancer risk factors. Scientific Reports, 9, 15055.CrossRefGoogle ScholarPubMed
Czyz, W., Morahan, J. M., Ebers, G. C., & Ramagopalan, S. V. (2012). Genetic, environmental and stochastic factors in monozygotic twin discordance with a focus on epigenetic differences. BMC Medicine, 10, 112.CrossRefGoogle ScholarPubMed
Durso, D. F., Bacalini, M. G., Sala, C., Pirazzini, C., Marasco, E., Bonafé, M., do Valle, Í. F., Gentilini, D., Castellani, G., Faria, A. M. C., Franceschi, C., Garagnani, P., & Nardini, C. (2017). Acceleration of leukocytes’ epigenetic age as an early tumorand sex-specific marker of breast and colorectal cancer. Oncotarget, 8, 2323723245.CrossRefGoogle Scholar
Fortin, J., Triche, T. J. Jr., & Hansen, K. D. (2017). Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics, 33, 558560.Google ScholarPubMed
Ginsberg, A., Price, S., Ingram, D., & Nottage, E. (1996). Life events and the risk of breast cancer: a case-control study. European Journal of Cancer, 32, 20492052.CrossRefGoogle Scholar
Groarke, E. M., & Young, N. S. (2019). Aging and hematopoiesis. Clinics in Geriatric Medicine, 35, 285293.CrossRefGoogle ScholarPubMed
Hannum, G., Guinney, J., Zhao, L., Zhang, L., Hughes, G., Sadda, S. V., Klotzle, B., Bibikova, M., Fan, J. B., Gao, Y., Deconde, R., Chen, M., Rajapakse, I., Friend, S., Ideker, T., & Zhang, K. (2013). Genome-wide methylation profiles reveal quantitative views of human aging rates. Molecular Cell, 49, 359367.CrossRefGoogle ScholarPubMed
Harris, J. R., Hjelmborg, J., Adami, H.-O., Czene, K., Mucci, L., Kaprio, J., & Nordic Twin Study of Cancer (NorTwinCan) Collaboration. (2019). The Nordic Twin Study on Cancer — NorTwinCan. Twin Research and Human Genetics, 22, 817823.CrossRefGoogle ScholarPubMed
Hofstatter, E. W., Horvath, S., Dalela, D., Gupta, P., Chagpar, A. B., Wali, V. B., Bossuyt, V., Storniolo, A. M., Hatzis, C., Patwardhan, G., Von Wahlde, M. K., Butler, M., Epstein, L., Stavris, K., Sturrock, T., Au, A., Kwei, S., & Pusztai, L. (2018). Increased epigenetic age in normal breast tissue from luminal breast cancer patients. Clinical Epigenetics, 10, 111.CrossRefGoogle ScholarPubMed
Horvath, S. (2013). DNA methylation age of human tissues and cell types. Genome Biology, 14, 115.CrossRefGoogle ScholarPubMed
Horvath, S., & Raj, K. (2018). DNA methylation-based biomarkers and the epigenetic clock theory of ageing. Nature Reviews Genetics, 19, 371384.CrossRefGoogle ScholarPubMed
Iorio, G. C., Spieler, B. O., Ricardi, U., & Dal Pra, A. (2021). The impact of pelvic nodal radiotherapy on hematologic toxicity: A systematic review with focus on leukopenia, lymphopenia and future perspectives in prostate cancer treatment. Critical Reviews in Oncology/Hematology, 168, 103497.CrossRefGoogle ScholarPubMed
Jylhävä, J., Hjelmborg, J., Soerensen, M., Munoz, E., Tan, Q., Kuja-Halkola, R., Mengel-From, J., Christensen, K., Christiansen, L., Hägg, S., Pedersen, N. L., & Reynolds, C. A. (2019). Longitudinal changes in the genetic and environmental influences on the epigenetic clocks across old age: Evidence from two twin cohorts. EBioMedicine, 40, 710716.CrossRefGoogle ScholarPubMed
Kaprio, J., Bollepalli, S., Buchwald, J., Iso-Markku, P., Korhonen, T., Kovanen, V., Kujala, U., Laakkonen, E. K., Latvala, A., Leskinen, T., Lindgren, N., Ollikainen, M., Piirtola, M., Rantanen, T., Rinne, J., Rose, R. J., Sillanpää, E., Silventoinen, K., Sipilä, S., … Waller, K. (2019). The Older Finnish Twin Cohort — 45 years of follow-up. Twin Research and Human Genetics, 22, 240254.CrossRefGoogle ScholarPubMed
Kresovich, J. K., Xu, Z., O’Brien, K. M., Shi, M., Weinberg, C. R., Sandler, D. P., & Taylor, J. A. (2022). Blood DNA methylation profiles improve breast cancer prediction. Molecular Oncology, 16, 42.CrossRefGoogle ScholarPubMed
Kresovich, J. K., Xu, Z., O’Brien, K. M., Weinberg, C. R., Sandler, D. P., & Taylor, J. A. (2019). Methylation-based biological age and breast cancer risk. Journal of the National Cancer Institute, 111, 10511058.CrossRefGoogle ScholarPubMed
Kune, S., Kune, G. A., Watson, L. F., & Rahe, R. H. (1991). Recent life change and large bowel cancer. Data from the Melbourne Colorectal Cancer Study. Journal of Clinical Epidemiology, 44, 5768.CrossRefGoogle ScholarPubMed
Levav, I., Kohn, R., Iscovich, J., Abramson, J. H., Tsai, W. Y., & Vigdorovich, D. (2000). Cancer incidence and survival following bereavement. American Journal of Public Health, 90, 16011607.Google ScholarPubMed
Levine, M. E., Hosgood, H. D., Chen, B., Absher, D., Assimes, T., & Horvath, S. (2015). DNA methylation age of blood predicts future onset of lung cancer in the women’s health initiative. Aging (Albany NY), 7, 690.CrossRefGoogle ScholarPubMed
Levine, M. E., Lu, A. T., Chen, B. H., Hernandez, D. G., Singleton, A. B., Ferrucci, L., Bandinelli, S., Assimes, T. L., & Horvath, S. (2016). Menopause accelerates biological aging. Proceedings of the National Academy of Sciences of the United States of America, 113, 93279332.CrossRefGoogle ScholarPubMed
Levine, M. E., Lu, A. T., Quach, A., Chen, B. H., Assimes, T. L., Bandinelli, S., Hou, L., Baccarelli, A. A., Stewart, J. D., Li, Y., Whitsel, E. A., Wilson, J. G., Reiner, A. P., Aviv, A., Lohman, K., Liu, Y., Ferrucci, L., & Horvath, S. (2018). An epigenetic biomarker of aging for lifespan and healthspan. Aging, 10, 573591.CrossRefGoogle ScholarPubMed
Li, S., Nguyen, T. L., Wong, E. M., Dugué, P. A., Dite, G. S., Armstrong, N. J., Craig, J. M., Mather, K. A., Sachdev, P. S., Saffery, R., Sung, J., Tan, Q., Thalamuthu, A., Milne, R. L., Giles, G. G., Southey, M. C., & Hopper, J. L. (2020). Genetic and environmental causes of variation in epigenetic aging across the lifespan. Clinical Epigenetics, 12, 112.CrossRefGoogle ScholarPubMed
Lillberg, K., Verkasalo, P. K., Kaprio, J., Teppo, L., Helenius, H., & Koskenvuo, M. (2003). Stressful life events and risk of breast cancer in 10,808 women: A cohort study. American Journal of Epidemiology, 157, 415423.CrossRefGoogle ScholarPubMed
Lin, Y., Wang, C., Zhong, Y., Huang, X., Peng, L., Shan, G., Wang, K., & Sun, Q. (2013). Striking life events associated with primary breast cancer susceptibility in women: A meta-analysis study. Journal of Experimental & Clinical Cancer Research, 32, 53.CrossRefGoogle ScholarPubMed
López-Otín, C., Blasco, M. A., Partridge, L., Serrano, M., & Kroemer, G. (2013). The hallmarks of aging. Cell, 153, 11941217.CrossRefGoogle ScholarPubMed
Lu, A. T., Quach, A., Wilson, J. G., Reiner, A. P., Aviv, A., Raj, K., Hou, L., Baccarelli, A. A., Li, Y., Stewart, J. D., Whitsel, E. A., Assimes, T. L., Ferrucci, L., & Horvath, S. (2019). DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging, 11, 303327.CrossRefGoogle ScholarPubMed
Pidsley, R., Wong, C. C. Y., Volta, M., Lunnon, K., Mill, J., & Schalkwyk, L. C. (2013). A data-driven approach to preprocessing Illumina 450K methylation array data. BMC Genomics, 14, 110.CrossRefGoogle ScholarPubMed
Qin, N., Li, Z., Song, N., Wilson, C. L., Easton, J., Mulder, H., Plyler, E., Neale, G., Walker, E., Zhou, X., Pan, H., Hudson, M. M., Yasui, Y., Robison, L. L., Zhang, J., Ness, K. K., & Wang, Z. (2021). Epigenetic age acceleration and chronic health conditions among adult survivors of childhood cancer. JNCI Journal of the National Cancer Institute, 113, 597605.CrossRefGoogle ScholarPubMed
Reynolds, C. A., Tan, Q., Munoz, E., Jylhävä, J., Hjelmborg, J., Christiansen, L., Hägg, S., & Pedersen, N. L. (2020). A decade of epigenetic change in aging twins: Genetic and environmental contributions to longitudinal DNA methylation. Aging Cell, 19, e13197.CrossRefGoogle ScholarPubMed
Santos, M. C. L., Horta, B. L., do Amaral, J. J. F., Branco, P. F. C., Fernandes, C., Galvão, C. M., & Fernandes, A. F. C. (2009). Association between stress and breast cancer in women: A meta-analysis. Cadernos de Saúde Pública, 25, S453S463.CrossRefGoogle ScholarPubMed
Sehl, M. E., Carroll, J. E., Horvath, S., & Bower, J. E. (2020). The acute effects of adjuvant radiation and chemotherapy on peripheral blood epigenetic age in early stage breast cancer patients. NPJ Breast Cancer, 6, 15.CrossRefGoogle ScholarPubMed
Sehl, M. E., Henry, J. E., Storniolo, A. M., Horvath, S., & Ganz, P. A. (2021). The effects of lifetime estrogen exposure on breast epigenetic age. Cancer Epidemiology Biomarkers and Prevention, 30, 12411249.CrossRefGoogle ScholarPubMed
Selya, A. S., Rose, J. S., Dierker, L. C., Hedeker, D., & Mermelstein, R. J. (2012). A Practical guide to calculating Cohen’s f2, a measure of local effect size, from PROC MIXED. Frontiers in Psychology, 3, 111.CrossRefGoogle ScholarPubMed
Skytthe, A., Harris, J. R., Czene, K., Mucci, L., Adami, H.-O., Christensen, K., Hjelmborg, J., Holm, N. V., Nilsen, T. S., Kaprio, J., & Pukkala, E. (2019). Cancer incidence and mortality in 260,000 Nordic twins with 30,000 prospective cancers. Twin Research and Human Genetics, 22, 99107.CrossRefGoogle ScholarPubMed
Snir, S., Farrell, C., & Pellegrini, M. (2019). Human epigenetic ageing is logarithmic with time across the entire lifespan. Epigenetics, 14, 912926.CrossRefGoogle ScholarPubMed
Spivak, J. L., Gascón, P., & Ludwig, H. (2009). Anemia management in oncology and hematology. The Oncologist, 14, 4356 CrossRefGoogle ScholarPubMed
Terry, M. B., McDonald, J. A., Wu, H. C., Eng, S., & Santella, R. M. (2016). Epigenetic biomarkers of breast cancer risk: Across the breast cancer prevention continuum. Advances in Experimental Medicine and Biology, 882, 3368.CrossRefGoogle ScholarPubMed
van Iterson, M., Tobi, E. W., Slieker, R. C., den Hollander, W., Luijk, R., Slagboom, P. E., & Heijmans, B. T. (2014). MethylAid: Visual and interactive quality control of large Illumina 450k datasets. Bioinformatics, 30, 34353437.CrossRefGoogle ScholarPubMed
Venkatesan, S., Swanton, C., Taylor, B. S., & Costello, J. F. (2017). Treatment-induced mutagenesis and selective pressures sculpt cancer evolution. Cold Spring Harbor Perspectives in Medicine, 7, a026617.CrossRefGoogle ScholarPubMed
Wang, Z., Hui, Q., Goldberg, J., Smith, N., Kaseer, B., Murrah, N., Levantsevych, O. M., Shallenberger, L., Diggers, E., Bremner, J. D., Vaccarino, V., & Sun, Y. V. (2022). Association between posttraumatic stress disorder and epigenetic age acceleration in a sample of twins. Psychosomatic Medicine, 84, 151158.CrossRefGoogle Scholar
Wolf, E. J., Logue, M. W., Stoop, T. B., Schichman, S. A., Stone, A., Sadeh, N., Hayes, J. P., & Miller, M. W. (2018). Accelerated DNA methylation age: Associations with Posttraumatic stress disorder and mortality. Psychosomatic Medicine, 80, 4248.CrossRefGoogle ScholarPubMed
Wolf, E. J., Maniates, H., Nugent, N., Maihofer, A. X., Armstrong, D., Ratanatharathorn, A., Ashley-Koch, A. E., Garrett, M., Kimbrel, N. A., Lori, A., Va Mid-Atlantic Mirecc Workgroup, Aiello, A. E., Baker, D. G., Beckham, J. C., Boks, M. P., Galea, S., Geuze, E., Hauser, M. A., Kessler, R. C., … Logue, M. W. (2018). Traumatic stress and accelerated DNA methylation age: A meta-analysis. Psychoneuroendocrinology, 92, 123.CrossRefGoogle ScholarPubMed
World Health Organization. (2019). ICD-10 Version:2019. https://icd.who.int/browse10/2019/en Google Scholar
Xiao, C., Beitler, J. J., Peng, G., Levine, M. E., Conneely, K. N., Zhao, H., Felger, J. C., Wommack, E. C., Chico, C. E., Jeon, S., Higgins, K. A., Shin, D. M., Saba, N. F., Burtness, B. A., Bruner, D. W., & Miller, A. H. (2021). Epigenetic age acceleration, fatigue, and inflammation in patients undergoing radiation therapy for head and neck cancer: A longitudinal study. Cancer, 127, 33613371.CrossRefGoogle ScholarPubMed
Xiao, C., Miller, A. H., Peng, G., Levine, M. E., Conneely, K. N., Zhao, H., Eldridge, R. C., Wommack, E. C., Jeon, S., Higgins, K. A., Shin, D. M., Saba, N. F., Smith, A. K., Burtness, B., Park, H. S., Irwin, M. L., Ferrucci, L. M., Ulrich, B., Qian, D. C., … Bruner, D. W. (2021). Association of epigenetic age acceleration with risk factors, survival, and quality of life in patients with head and neck cancer. International Journal of Radiation Oncology, Biology, Physics, 111, 157167.CrossRefGoogle ScholarPubMed
Yang, R., Wu, G. W. Y., Verhoeven, J. E., Gautam, A., Reus, V. I., Kang, J. I., Flory, J. D., Abu-Amara, D., Systems Biology Consortium, PTSD, Hood, L., Doyle, F. J. 3rd, Yehuda, R., Marmar, C. R., Jett, M., Hammamieh, R., Mellon, S. H., & Wolkowitz, O. M. (2020). A DNA methylation clock associated with age-related illnesses and mortality is accelerated in men with combat PTSD. Molecular Psychiatry, 26, 49995009.CrossRefGoogle ScholarPubMed
Yu, M., Hazelton, W. D., Luebeck, G. E., & Grady, W. M. (2020). Epigenetic aging: More than just a clock when it comes to cancer. Cancer Research, 80, 367374.CrossRefGoogle ScholarPubMed
Zheng, Y., Joyce, B. T., Colicino, E., Liu, L., Zhang, W., Dai, Q., Shrubsole, M. J., Kibbe, W. A., Gao, T., Zhang, Z., Jafari, N., Vokonas, P., Schwartz, J., Baccarelli, A. A., Baccarelli, A. A., & Hou, L. (2016). Blood epigenetic age may predict cancer incidence and mortality. EBioMedicine, 5, 6873.CrossRefGoogle ScholarPubMed
Zhou, W., Laird, P. W., & Shen, H. (2017). Comprehensive characterization, annotation and innovative use of Infinium DNA methylation BeadChip probes. Nucleic Acids Research, 45, e22e22.Google ScholarPubMed
Figure 0

Table 1. Number of the Finnish Twin Cohort twin pairs discordant for cancer shown for the different cancer groups

Figure 1

Table 2. Description of the cancer-discordant Finnish Twin Cohort twin pairs

Figure 2

Table 3. Within-pair differences in DNAmAA in cancer-discordant twin pairs (N = 95). The presented units are in years

Figure 3

Table 4. Within-pair differences in DNAmAA in breast cancer-discordant twin pairs (N = 24). The presented units are in years

Figure 4

Table 5. Within-pair differences in DNAmAA in nonbreast cancer-discordant twin pairs (N = 71). The presented units are in years

Supplementary material: PDF

Bode et al. supplementary material

Bode et al. supplementary material
Download Bode et al. supplementary material(PDF)
PDF 2 MB