Hannah-Ruth Engelbrecht1,2,3,4
Sarah M. Merrill1,2,3,4
Nicole Gladish5
Julie L. MacIsaac2,3,4David T. S. Lin2,3,4
Simone Ecker6
Christina A. Chrysohoou7
Giovanni M. Pes8Michael S. Kobor1,2,3,4*†David H. Rehkopf5*†- 1Edwin S. H. Leong Healthy Aging Program, Faculty of Medicine, University of British Columbia, Vancouver, BC, Canada
- 2Department of Medical Genetics, Faculty of Medicine, University of British Columbia, Vancouver, BC, Canada
- 3British Columbia Children’s Hospital Research Institute, University of British Columbia, Vancouver, BC, Canada
- 4Centre for Molecular Medicine and Therapeutics, University of British Columbia, Vancouver, BC, Canada
- 5Department of Epidemiology and Population Health, School of Medicine, Stanford University, Palo Alto, CA, United States
- 6UCL Cancer Institute, University College London, London, United Kingdom
- 7First Cardiology Clinic, School of Medicine, University of Athens, Athens, Greece
- 8Department of Clinical and Experimental Medicine, University of Sassari, Sassari, Italy
Sex differences in aging manifest in disparities in disease prevalence, physical health, and lifespan, where women tend to have greater longevity relative to men. However, in the Mediterranean Blue Zones of Sardinia (Italy) and Ikaria (Greece) are regions of centenarian abundance, male-female centenarian ratios are approximately one, diverging from the typical trend and making these useful regions in which to study sex differences of the oldest old. Additionally, these regions can be investigated as examples of healthy aging relative to other populations. DNA methylation (DNAm)-based predictors have been developed to assess various health biomarkers, including biological age, Pace of Aging, serum interleukin-6 (IL-6), and telomere length. Epigenetic clocks are biological age predictors whose deviation from chronological age has been indicative of relative health differences between individuals, making these useful tools for interrogating these differences in aging. We assessed sex differences between the Horvath, Hannum, GrimAge, PhenoAge, Skin and Blood, and Pace of Aging predictors from individuals in two Mediterranean Blue Zones and found that men displayed positive epigenetic age acceleration (EAA) compared to women according to all clocks, with significantly greater rates according to GrimAge (β = 3.55; p = 1.22 × 10−12), Horvath (β = 1.07; p = 0.00378) and the Pace of Aging (β = 0.0344; p = 1.77 × 10−08). Other DNAm-based biomarkers findings indicated that men had lower DNAm-predicted serum IL-6 scores (β = -0.00301, p = 2.84 × 10−12), while women displayed higher DNAm-predicted proportions of regulatory T cells than men from the Blue Zone (p = 0.0150, 95% Confidence Interval [0.00131, 0.0117], Cohen’s d = 0.517). All clocks showed better correlations with chronological age in women from the Blue Zones than men, but all clocks showed large mean absolute errors (MAE >30 years) in both sexes, except for PhenoAge (MAE <5 years). Thus, despite their equal survival to older ages in these Mediterranean Blue Zones, men in these regions remain biologically older by most measured DNAm-derived metrics than women, with the exception of the IL-6 score and proportion of regulatory T cells.
Introduction
Men and women seem to age differently—women tend to have greater morbidity and frailty in older ages but also live longer in most populations (Oksuzyan et al., 2008; Thinggaard et al., 2016; Hägg and Jylhävä, 2021). The diseases affecting aging men and women are also different (Hägg and Jylhävä, 2021). Chronic, non-lethal conditions and reduced health span disproportionately affect women more than men. Conversely, men display greater rates of heart disease, the leading cause of death globally, and higher incidence of non-reproductive cancers than women (Bots et al., 2017; Carmel, 2019; Meyer et al., 2020; Garmany et al., 2021). Given the lifespan and health span discrepancies between the sexes, studies of the oldest old (>80 years) provide useful opportunities to examine sex-related differences in healthy aging and biomarkers of the biological aging process, as the extended lifespans in these populations have allowed the accumulation of changes associated with aging (Puca et al., 2018).
While living to the age of 100 in most parts of the world is relatively unlikely, there exist four naturally occurring, validated, geographic regions of exceptional longevity, known as Blue Zones, with a high prevalence of centenarians (Buettner, 2012; Poulain et al., 2013): Sardinia (Italy), Ikaria (Greece), Nicoya Peninsula (Costa Rica) and Okinawa (Japan). These regions are geographically isolated, either being islands or separated from other regions by natural features such as mountains, and there is typically a traditional, “unmodernized” agricultural lifestyle observed (Poulain et al., 2011; Poulain and Mackowicz, 2021). A second unique feature in the Mediterranean Blue Zones, Ikaria and Sardinia (Poulain et al., 2011; Poulain and Mackowicz, 2021), is that equal proportions of men and women survive to extreme old age, in stark contrast to the preponderance of oldest old women commonly observed in most global populations (Passarino et al., 2002; Austad, 2011). While many aspects of lifestyle have been investigated in these locations, including psychological outlook, physical activity, and sociodemographic characterizations (Fastame et al., 2018; Hitchcott et al., 2018; Pes et al., 2018; Pes et al., 2020), there remains limited information regarding biological age differences between these men and women who survive to become centenarians more equally than the average global population. Biological age is a prediction of an individual’s chronological age that represents their functioning and health at a point in time, and has been used as a measure of well-being for many decades using a wide variety of health metrics (Furukawa et al., 1975; Rockwood et al., 2005; Jylhävä et al., 2017; Ji et al., 2021).
Epigenetic predictors have been developed to predict either chronological age or other biomarkers based on the DNA methylation (DNAm) levels in a particular combination of sites in the genome. DNAm, the addition of a methyl group to primarily cytosine-phosphate-guanine (CpG) sites of DNA, is a necessary part of the genomic regulatory suite that can change in response to external exposures and has both temporally stochastic changes and predictable changes at certain sites with increasing age (Horvath, 2013; Seale et al., 2022). Deviations between predicted epigenetic age and chronological age have been associated with differences in health, which may be related to behavioral and environmental exposures over the life course. With this conceptualization, lower predicted epigenetic age is typically indicative of better health, while higher predicted epigenetic age is indicative of poorer health (Jones et al., 2015). This has been demonstrated in many aging related conditions. For example, cancers and cardiovascular disease have been associated with increased epigenetic age (Perna et al., 2016; Soler-Botija et al., 2019). Epigenetic age predictors that are based on these consistent changes in DNAm are commonly termed epigenetic clocks. Many clocks exist, each using slightly different additive models of CpG site DNAm states to reach a prediction of epigenetic age, with some newer clocks incorporating other biological variables, such as DNAm-based predictors of serum protein levels (Bergsma and Rogaeva, 2020). In addition to epigenetic clocks, there are DNAm-based epigenetic predictors which estimate other age-associated biomarkers such as telomere length and cytokine levels, which have also been linked to health outcomes (Pusceddu et al., 2018; Villar-Fincheira et al., 2021).
Epigenetic clocks have been developed to detect different aspects of biological age using strict chronological age (Horvath, Hannum, Skin and Blood) (Hannum et al., 2013; Horvath, 2013; Horvath et al., 2018), mortality (GrimAge) (Lu A. T. et al., 2019), and general characteristics of aging such as frailty and increased heart disease risk (termed “phenotypic age”) (PhenoAge) (Levine et al., 2018). Furthermore, each epigenetic clock can be used to calculate a second measure, epigenetic age acceleration (EAA), the “ticking rate” of the clock. EAA is calculated by finding the residuals from the regression of the predicted epigenetic age on chronological age and can also be used to interrogate the rate at which aging is occurring. The aforementioned epigenetic clocks have been shown to discriminate between physical health and disease status (Quach et al., 2017; Hillary et al., 2020), and cognitive decline (Marioni et al., 2015), which also display sex differences in the aging process, positioning clocks as useful tools to investigate the sex differences in aging (Ahrenfeldt et al., 2019; Levine et al., 2021). According to a recently developed deep-learning based clock, epigenetic age in men accelerates sooner than in women (>55 years), supporting the ability of these epigenetic clocks to discriminate sex and age (Galkin et al., 2021). However, current epigenetic clocks were trained in predominantly middle-aged cohorts, with mean ages of approximately ∼50 years, and as such, their performance and behavior have not been thoroughly examined at the extreme upper end of life.
Additional DNAm-based biomarkers that complement epigenetic clocks have been developed, including a measure using DNAm to predict the annual rate of physiological decline (Pace of Aging) as a “cousin” measure to the clocks, a measure of telomere length (DNAmTL) (Lu A. T. et al., 2019), serum interleukin-6 (IL-6) (Stevenson et al., 2021), smoking status and score, and cell type proportions. The Pace of Aging was trained by using a longitudinal cohort of identically-aged individuals and monitoring the changes in DNAm at particular sites over time between 19 and 45 years of age. Thus, Pace of Aging represents a rate of physiological decline per annum (Belsky et al., 2022). This measure has been associated with health outcomes including dementia and frailty (Verschoor et al., 2021; Sugden et al., 2022), important features of the aging process, which are known to display sex differences (Anstey et al., 2021; Sindi et al., 2021). However, as the most recent version of the Pace of Aging has only been produced using individuals with a maximum age of 45 years, investigations into how this measure performs in a cohort of centenarians is needed to further characterize previously established associations and their relevance in the geroscience field.
Other biomarkers that have also been well-investigated in studies of both longevity and sex differences are telomere length and serum IL-6. Telomeres, the protective hexamer repeats that cap chromosomes, shorten over the lifespan, but are generally longer in women at all ages (Vaiserman and Krasnienkov, 2021). Studies of individuals from the Blue Zone in the Nicoya Peninsula in Costa Rica have indicated that Blue Zone residents have longer telomers than individuals from other regions in Costa Rica (Rehkopf et al., 2013). The longer telomeres observed in women have been postulated as both a mechanism and biomarker of longevity (Gardner et al., 2014; Öngel et al., 2021). However, telomere length has not been able to resolve aspects of healthy aging other than longevity (Arai et al., 2015). IL-6, an inflammatory cytokine, is known to increase with age, and is part of the chronic, low-grade, systemic inflammation that occurs with aging, known as inflammaging (Ferrucci and Fabbri, 2018). Centenarians have been found to possess lower inflammation measured by scores that include IL-6 as a main factor, potentially as a result of a counteracting anti-inflammaging response, and such findings have been associated with their better health (Arai et al., 2015; Minciullo et al., 2016). Smoking behavior has been associated with the risk of cardiovascular disease development in a sex-dependent manner (Vasiljevic et al., 2021), which has implications for the health of individuals. Cell-type differences predicted with DNAm have also been shown to differentiate amongst the oldest old in the Costa Rican Blue Zone, where a lower proportion of CD8+ memory T cells and higher proportion of naïve T cells were observed in Blue Zone residents compared to non-Blue Zone residents, indicating a younger immune system profile (McEwen et al., 2017). Similar to epigenetic clocks, the majority of these DNAm biomarker predictors have been developed or investigated predominantly in middle aged to young elderly cohorts (∼70 years), necessitating a characterization among the oldest old as well.
Focusing on the lack of characterization of epigenetic predictors among elderly populations, in these analyses we seek to characterize the performance of several predictors in a cohort of extremely elderly individuals from the Mediterranean Blue Zones, who remain understudied given their rarity. Initially, epigenetic clocks will be evaluated to assess at a molecular level whether long-surviving men from the Blue Zones display a comparable biological age to long-surviving women, as expected by their documented similar mortality rates. Subsequently, we use DNAm-based predictors of telomere length (DNAmTL), smoking score (EpiSmokER), and serum IL-6 (DNAm IL-6 score) which may indicate health differences in the oldest old from regions of exceptional longevity, reporting any differences related to sex (Lu A. T. et al., 2019; Bollepalli et al., 2019; Stevenson et al., 2021).
Methods
Cohort recruitment, sample collection, and data collection
A subset of participants from two separate studies, one in the Ikaria Blue Zone in Greece and one in the Sardinia Blue Zone in Italy, were selected for DNAm analysis (Pes et al., 2020; Foscolou et al., 2021) (Table 1). Given that these regions are in relatively close proximity geographically when compared to the other Blue Zones in Costa Rica and Japan, have relatively similar climates, diets, and cultures, and the same male survival phenomenon is observed, we combined the Sardinian and Ikarian Blue Zone residents into one cohort to represent the Mediterranean Blue Zones (Poulain et al., 2021; Pes et al., 2022). Venous blood was collected from participants with informed consent from both Blue Zone regions as part of each study. Age at time of blood collection and self-reported sex were used as the variables of interest in this study. DNA was extracted from the blood samples and shipped to the University of British Columbia (Vancouver, British Columbia, Canada). DNA was bisfulfite treated using the EZ-DNA methylation kit (Zymo Research, CA, United States), and DNAm data were measured using the MethylationEPIC BeadChip array (“EPIC”, Illumina, San Diego, CA, United States) according to manufacturer’s protocols. This study was approved by the Committee on Human Subjects at Stanford University.
TABLE 1. Age, age range, and number of participants by region and sex.
While matching Italian and Greek control participants from beyond the Blue Zone regions were recruited using the same process, the participants recruited from the Blue Zone regions were markedly older than the participants recruited from matching control regions, and there were far fewer control region participants than Blue Zone residents (Supplementary Figure S1). Given that health and DNAm can change quite dramatically in the extreme end of life, that age acts as a risk factor for most health-related outcomes (Niccoli and Partridge, 2012), and given that the sample size was underpowered for any statistical analysis, control participants could not be included in any statistical comparisons.
Data normalization, preprocessing
The Illumina BeadArray Data (IDAT) files from the array run were imported into R 4.0.3 (R Core Team, 2022). Poorly performing samples which were identified as outliers based on the ratio of the log2 unmethylated to methylated intensity of the control probes were removed (n = 3). Functional normalization using the preprocessFunnorm function from minfi was performed to adjust background and probe type intensity and (Fortin et al., 2014). Any samples identified by detectOutlier from the lumi package (Du et al., 2008), pcout from wateRmelon (Pidsley et al., 2013), locfdr from the locFDR package (Hannum et al., 2013), or pfilter from wateRmelon (Pidsley et al., 2013), were deemed outliers and removed (n = 4). DNAm-based sex prediction was performed using the getSex function from the minfi suite based on DNAm data on the X and Y chromosomes (Fortin et al., 2014). In addition to this measure, self-reported sex and DNAm-based log2 intensity of XY chromosome probes estimated sex were confirmed against each other, which led to the removal of 7 sex-mismatched samples. A second confirmation of sex was completed by using a k-means clustering approach of the DNAm from the 16,839 probes on the X and Y chromosomes, and no further mismatches were identified. The normalized cohort consisted of 112 individuals from the Blue Zone and control regions. Further processing involved the removal of a selection of probes from the EPIC array as is standard in the field, including SNP probes (n = 13,110), known polymorphic probes (n = 23,906), sex chromosome probes (n = 17,344), cross-hybridizing probes (n = 9,441) (Price et al., 2013), probes represented by three beads or less, failed detection in ≥1% of samples, or had a detection p-value > 1 × 10−16 (n = 59,154) (Pidsley et al., 2016; Fortin et al., 2017). The ComBat function from the sva package (Leek et al., 2019, 2012) was used for technical correction of plate, row position, and chip.
Epigenetic age prediction and DNAm-derived biological proxy measure calculation
After all processing steps were completed, the resulting DNAm beta matrix was uploaded to the DNA Methylation Age Calculator1 (Horvath, 2013) for epigenetic age prediction. When uploaded to the DNA Methylation Age Calculator, “Normalize Data” and “Advanced Analysis” were selected as options. When selecting the “Advanced Analysis” option, the DNA Methylation Age Calculator provides the Horvath, Hannum, Skin and Blood, PhenoAge, GrimAge, and DNAmTL predictions. This same DNAm data was used to calculate the Pace of Aging using the DunedinPACE predictor R package with the PoAmProjector function2 (Belsky et al., 2022). Epigenetic age acceleration (EAA) was calculated for the epigenetic clocks (Horvath, Hannum, Skin and Blood, PhenoAge and GrimAge) by obtaining the residuals from the linear model: Epigenetic age ∼ Chronological age. The features of the predictors used is provided in Supplementary Table S1.
In addition to the explicitly age-related DNAm-based biomarkers obtained above, two other DNAm-based biomarkers were calculated which predicts smoking and serum IL-6 levels. Smoking status and score provide an indication of smoking behavior and status and were calculated using EpiSmokEr in R (Bollepalli et al., 2019): Smoking status was calculated based on the DNAm of 121 CpGs identified by Bollepati et al., and smoking score was calculated based on DNAm of 187 previously identified CpGs (Elliott et al., 2014). The DNAm-based serum IL-6 score was calculated in R using a weighted score from the DNAm status of 12 CpGs (Stevenson et al., 2021). A summary of the epigenetic clocks, biomarkers, and their main features is provided in Supplementary Table S1. The data preprocessing and normalization steps selected for the prediction of epigenetic age were assessed based on the prediction concordance between chronological age-prediction clocks (Supplemental Figures 2 and 3).
Prediction and transformation of cell type proportions
DNAm is a key component of cell identity, and the source of DNAm data in this study is a heterogenous tissue, venous whole blood. To resolve and account for cell-based DNAm signal in downstream analyses, cell type deconvolution (a means of predicting cell type proportions in the tissue sample) was completed with constraint projections using the IDOL package in R for extended blood cell types (Koestler et al., 2016; Gervin et al., 2019; Salas et al., 2022). Prior to use in linear models, the cell type proportions were transformed using robust isometric principal components (PCs) (Filzmoser et al., 2009).
Assessment of epigenetic clock performance
Pearson’s correlation coefficients (r) between the chronological age of individuals and their predicted epigenetic age across the Horvath, Hannum, GrimAge, PhenoAge, and Skin and Blood epigenetic clocks were compared for men and women of the Mediterranean Blue Zones. As the remaining epigenetic predictors (Pace of Aging, DNAmTL, cell type proportions, IL-6 and smoking score) do not measure age in years, they were not included in these comparisons of year-based outcomes. The mean absolute error and maximum absolute error (MAE and MaxAE) were calculated for each epigenetic clock’s age prediction and chronological age using the mae and maxae functions in the R package mlr3measures (Lang et al., 2019).
Statistical model of biomarker comparisons
For the purposes of this study, the statistical significance threshold was defined as a p-value of 0.05 Bonferroni-corrected for nine predictor comparisons (p = 5.56 × 10–03). Linear regression was performed to compare differences between the sexes. In order to account for the low sample size and ensure the models were not overfitted, the covariates included: sex, age, and the first four PCs accounting for 80% of cell type, with female sex being considered the reference: Biomarker ∼ Sex + Age + Cell type PC1 + Cell type PC2 + Cell type PC3 + Cell type PC4.
Cell type proportion comparisons
The 12 cell type proportions were compared between men and women using t-tests of mean proportion per cell type. Tukey’s Honest Significant Difference test and resulting adjusted p-value were used to assess the significance of cell types with a p-value slightly <0.05.
Results
Epigenetic clocks perform differently between the sexes amongst the oldest old in the Blue Zones
As most epigenetic clocks are trained on approximately middle-aged individuals (Supplementary Table S1), we first assessed the performance of the Horvath, Hannum, Skin and Blood, GrimAge, and PhenoAge epigenetic clocks in this group of the oldest old. We identified that all epigenetic age predictors had higher correlations with chronological age in women (r = 0.52–0.92 in women and r = 0.46–0.87 in men across the tested clocks) (Table 2; Figure 1). PhenoAge demonstrated the lowest correlation with chronological age (r < 0.5 in men and women), while GrimAge outperformed the other clocks with strong correlations among both sexes (r > 0.8) (Table 2; Figure 1). However, the mean and maximum errors were consistently quite large for both sexes across all clocks (>30 years of error) except PhenoAge, which showed mean and maximum errors <5 years for both men and women, the smallest differences between predicted and chronological age (Table 2; Figure 2).
TABLE 2. Pearson’s correlations (r), Mean Absolute Error, and Maximum Absolute Error of the Horvath, Hannum, PhenoAge, GrimAge, and Skin and Blood epigenetic clocks in Men and Women from the Mediterranean Blue Zones.
FIGURE 1. Pearson’s correlation coefficients (r) between predicted epigenetic age and chronological age for several epigenetic clocks in (A) Female and (B) Male residents of the Mediterranean Blue Zones.
FIGURE 2. A Loess-smoothed curve of epigenetic predicted age with 95% confidence interval in grey against chronological age in years according to the PhenoAge, GrimAge, Hannum, Horvath, and Skin and Blood epigenetic age predictors in (A) female and (B) male residents of the Mediterranean Blue Zones.
Epigenetic age acceleration evident in oldest old men compared to oldest old women from the Blue Zones
To explore sex differences in biological age measures in this group of oldest old individuals from the Mediterranean Blue Zones, we examined the differences in the Pace of Aging measure and EAA across the Horvath, Hannum, GrimAge, PhenoAge, and Skin and Blood clocks. Men showed greater EAA than women across all epigenetic age predictors, with significant differences observed for the GrimAge (Adjusted R2 = 0.501, F (6, 87) = 16.6, p = 1.22 × 10–12, Cohen’s f2 = 1.00), Horvath (Adjusted R2 = 0.139, F (6,87) = 3.50, p = 0.00378, Cohen’s f2 = 0.161), and Pace of Aging predictors (Adjusted R2 = 0.383, F (6, 87) = 9.23, p = 1.77 × 10–08, Cohen’s f2 = 0.620) (Figure 3; Table 3).
FIGURE 3. Epigenetic age acceleration and Pace of Aging in male and female residents of the Mediterranean Blue Zones, according to the (A) Horvath, (B) Hannum, (C) PhenoAge, (D) GrimAge, (E) Skin and Blood, and (F) Pace of Aging predictors.
TABLE 3. Positive Epigenetic Age Acceleration (EAA) is evident in men from the Blue Zone relative to women from the Blue Zone. The coefficient estimate indicates the linear model β coefficient associated with the specified independent variable in brackets, and Cohen’s f2 effect size based on the Adjusted R2 of the linear model. Bolded p-values indicate significance (corrected p-value < 0.00556).
Regulatory T cell proportions higher in women from the Mediterranean Blue Zones
DNAm-predicted cell type proportions were used to make comparisons between the sexes. We observed that the oldest old women have significantly greater proportions of regulatory T cells (Tregs) than the oldest old men in the group (t (85.85) = 2.49, p = 0.0150, 95% CI [0.00131, 0.0117], Cohen’s d = 0.517) (Figure 4).
FIGURE 4. Epigenetically predicted cell type proportion by sex amongst the oldest old in the Mediterranean Blue Zones, including basophils (Bas), memory B cells (Bmem), naïve B cells (Bnv), memory CD4 T cells (CD4mem), memory CD8 T cells (CD8mem), naïve CD8 T cells (CD8nv), eosinophils (Eos), Monocytes (Mono), Neutrophils (Neu), Natural Killer cells (NK), and regulatory T cells (Treg).
Smoking score, predicted IL-6 score, and telomere length differ between the sexes in Mediterranean Blue Zones
We assessed smoking status and smoking score in the Mediterranean Blue Zone residents to determine whether the scores differed between the sexes of the oldest. Of the 94 Mediterranean Blue Zone residents, four were predicted to have never smoked, while 90 participants were classified as “former smokers” using DNAm-based detection methods in EpiSmokER. Men were not more likely to have smoked than women according to comparisons of smoking score. To examine potential signals of inflammaging, we assessed sex differences in predicted serum IL-6 scores between oldest old men and women from the Blue Zones. Men were predicted to have a lower serum IL-6 score than women (Adjusted R2 = 0.49, F (6, 87) = 15.94, p = 2.84 × 10–12, Cohen’s f2 = 0.96) in this group. The oldest old men also displayed shorter telomere length than women according to DNAmTL predictions (Adjusted R2 = 0.30, F (6, 87) = 7.66, p = 1.28 × 10–06, Cohen’s f2 = 0.04) (Figure 5).
FIGURE 5. Epigenetically predicted biomarkers by sex amongst the oldest old in the Mediterranean Blue Zones, (A) serum IL-6, (B) Smoking Score, and (C) leukocyte telomere length (DNAmTL).
Discussion
The aims of our study were to characterize the performance of the epigenetic clocks within a group of the oldest old, and secondly, to investigate markers of biological age in a unique longevous population, located in the Mediterranean Blue Zones, to develop a more comprehensive understanding of the sex differences that occur with extreme age. Men in these longevity regions, Sardinia and Ikaria, tend to survive into extreme old age at approximately the same rate as women, which suggest that they might have a health advantage relative to men from beyond these regions (Poulain et al., 2021). We assessed the performance of the epigenetic clocks in the Mediterranean Blue Zones, and all epigenetic clocks were found to have better correlations with chronological age in women compared to men, though a high error rate was found across all clocks except PhenoAge for both men and women. We then examined sex differences of EAA and observed that men had increased EAA, Pace of Aging, and smoking score relative to women from the Mediterranean Blue Zones (Horvath, 2013; Lu A. T. et al., 2019; Belsky et al., 2022). Though we observed minimal differences in cell type proportions, there was the exception in the proportion of Tregs where women were predicted to have higher levels than men. Additionally, we observed differences in DNAm-based biomarkers of telomere length and serum IL-6 where men were predicted to have shorter telomeres and lower levels of IL-6 relative to women. In general, DNAm-derived predictors of age and age-related biomarkers indicate men in these Mediterranean Blue Zones age biologically faster than women, despite their apparent survival advantage relative to men from other regions.
An examination of the performance of the epigenetic clocks indicated an overall poor performance of the epigenetic clocks in the Mediterranean Blue Zone group of the oldest old. The clocks had either poor correlations of epigenetic and chronological age, and/or large MAE and MaxAE rates according to every epigenetic age predictor assessed in both sexes. GrimAge, despite having a large error of approximately 30 years, along with the Horvath, Hannum, and Skin and Blood epigenetic clocks, also had the highest correlation with chronological age in men and women from the Mediterranean Blue Zones, which is comparable with epigenetic clock performance in other elderly cohorts (Ecker and Beck, 2019). We observed that all correlations of epigenetic and chronological age were stronger in women than men from the Mediterranean Blue Zones. Of all the tested clocks, GrimAge had the strongest correlation in women, indicating the aging health trajectory in this group of oldest old women is likely well-accounted for by this clock. Interestingly, all correlations between epigenetic age and chronological age were consistently lower in men from this region. The clocks tested here were not developed using populations with evidence of male longevity, and while each clock captures a unique aspect of the broader aging process, they currently do so in a manner that reflects female longevity. It may be that potential DNAm differences associated with the male longevity phenomenon of the Mediterranean Blue Zones are not captured by the current epigenetic clocks, which may indicate why increased male EAA, as is typical in most populations, is observed (Horvath et al., 2016).
The best performing clock by correlation was GrimAge, which is the only predictor tested here that incorporates a DNAm-based prediction of smoking pack-years (Lu A. T. et al., 2019). As observed via the EpiSmokER results, the vast majority of this cohort have smoked at some point in their lives, though this is very expected in this population, as smoking is quite common in the areas of Italy and Greece where Sardinia and Ikaria are located (Lugo et al., 2013; Christopoulou et al., 2021; Gallus et al., 2021). In fact, increased GrimAge, Horvath, and the Pace of Aging, which showed greater EAA and pace of aging in men, have all been associated with smoking (Klopack et al., 2022). GrimAge may display the highest correlation with chronological age because it is accounting for the molecular signature of well-established modifier of DNAm and health traits that may have sex-specific effects (Koo et al., 2021), even without sex-differential rates of smoking (Maas et al., 2020), and was developed using the oldest training cohort (Lu A. T. et al., 2019). It has also been observed that GrimAge EAA has more associations with other features of older biological age than PhenoAge, Horvath, and Hannum EAA, such as increased polypharmacy, lower grip strength, and slower walking speed, many of which display sex-differences themselves (Hägg and Jylhävä, 2021; McCrory et al., 2021). This suggests that GrimAge may be better suited to resolving differences in health status in the oldest old, rather than accurately predicting chronological age of the oldest old.
In contrast to the GrimAge EAA findings, we observed low correlations for PhenoAge and chronological age, yet also the lowest margins of error. PhenoAge clock training also included serum protein and cell counts in its training, but it has little similarity to those markers of aging used by GrimAge. Additionally, PhenoAge does not use DNAm-based predictions of these serum proteins, as it was trained in a two-stage process, which defined phenotypic age (based on serum proteins and white blood cell counts) in one cohort, and then trained to detect that phenotypic age using CpGs in a second cohort. The serum proteins used by these epigenetic clocks (GrimAge and PhenoAge) are distinct, with minimal similarity, barring C-reactive protein, which is used by both (Levine et al., 2018; Lu A. T. et al., 2019). The low error rates of PhenoAge in this cohort of oldest old individuals suggest that PhenoAge may be better suited to predicting oldest old age, and may be detecting general physiological trends that occur in the aging process. However, given the low correlation between predicted epigenetic and chronological age, PhenoAge is less able to resolve differences in health amongst the oldest old.
While women are, on average, less physically well than age-matched male counterparts, very few measures explain the lifespan gap that occurs. Epigenetic age has been well-associated with sex-specific features of aging, including frailty (Breitling et al., 2016; Maddock et al., 2019), cardiovascular disease (Roetker et al., 2018), and cognitive decline (Maddock et al., 2019; Vaccarino et al., 2021). Both telomere length and epigenetic age, however, better reflect the survival of women, although these molecular measures are not well correlated and cannot easily be compared to each other (Marioni et al., 2016). Women generally have longer telomeres than men at all ages (Gardner et al., 2014; Tucker, 2019), and men are found to have greater epigenetic age and EAA than women, with this effect becoming more pronounced at advanced ages (Horvath et al., 2016; Galkin et al., 2021; Kankaanpää et al., 2021). In the current Mediterranean Blue Zone population with equal male-female survival to advanced age, the commonly observed female aging advantage was also observed here by slower EAA in the most discriminating epigenetic clock, as well as longer predicted DNAm-derived telomere length. Additionally, the Pace of Aging predictor—one of the most current and well-validated predictors (Belsky et al., 2022, 2020)—produced higher scores among men relative to women, indicating that these molecular biomarkers of biological age may reflect the generally observed female survival to older ages. However, this is surprising given the hypothesized similarity of men’s and women’s biological aging, based on their similar survival rates in the Mediterranean Blue Zones. Therefore, it is reasonable to conjecture that these DNAm-derived measures do not capture the factors driving the protective effect evidenced by male survival to oldest old ages in the Mediterranean Blue Zones (Crimmins et al., 2021; Kankaanpää et al., 2021).
Alternatively, we did observe a lower DNAm-predicted serum IL-6 score in men than in women. Unlike EAA or DNAm-derived telomere length, this finding aligns with the predicted extended survival effect of men in these unique areas, as lower levels of serum IL-6 have been associated with healthy aging based on cognitive and physical ability (Puzianowska-Kuźnicka et al., 2016), thus potentially reflecting extended male survival. IL-6 is a cytokine secreted by several cell types during inflammation, which acts as a signaling molecule (Papanicolaou, 2000; Frasca and Blomberg, 2016; Kaur et al., 2020), and importantly, increases with age (Ershler, 1993). Previous research regarding the immune system of Sardinian centenarians has suggested a younger immune profile based on serum neopterin concentrations (a compound produced by macrophages in response to infection), which indicates lesser inflammation than expected for their age (Sotgia et al., 2017). The lower predicted IL-6 observed in this study in men from the Mediterranean Blue Zones further contributes to what may be an indicator of the survival advantage in men of this Blue Zone cohort, and suggests a modified inflammaging phenotype, specifically, rather than a general biological aging advantage.
Previous research interrogating immune system differences in the Nicoya Peninsula Blue Zone observed lower proportions of CD8+ memory T cells and higher proportions of naive T cells, suggesting a general younger immune profile (McEwen et al., 2017). In this study, men from the Mediterranean Blue Zones were predicted to have lower proportions of regulatory T cells (Tregs) than women. Tregs are usually involved in maintaining a level of homeostasis within the immune system by suppressing the immune response and promoting self-tolerance, but at older ages these cells tend to accumulate. Tregs play important protective roles in limiting autoimmune disorders (Tamosiuniene et al., 2018; Goodman et al., 2020), but functional Tregs are susceptible to both epigenetic alteration of the usually-expressed FOXP3 gene (Shu et al., 2017) and impaired estrogen signaling (Goodman et al., 2020, 2014), which may impair immunity and increase susceptibility to autoimmune diseases. The accumulation of dysregulated functional Tregs may contribute to immune system dysregulation by overly suppressing immune responses (Raynor et al., 2012; Klein and Flanagan, 2016; Müller et al., 2019; Churov et al., 2020; Rocamora-Reverte et al., 2021). Thus, the lower proportion of these Tregs in men from the Mediterranean Blue Zones compared to women may provide biological support for the similar male-female survival rate into old age and suggests an avenue worth investigating into male survival.
Sex differences in DNAm are known to exist (Grant et al., 2022), but the epigenetic clocks are trained to use age, aging phenotype, and mortality-related CpG sites as input, with the assumption that the age-related sites are likely the same in both sexes across the age spectrum. As found by McCrory et al. (2020), however, cardiovascular and metabolic conditions were more strongly associated with higher PhenoAge EAA in women, while Horvath EAA had greater associations with metabolic conditions in men. It is possible that in these understudied oldest old groups that age and age-related disease affected CpG sites become more sex-specific, or represent different aspects of the aging process that are sex-specific (McCrory et al., 2020). It is uncommon to have a very large cohort where oldest old men and women from a relatively homogenous environment can be compared, and ostensibly further studies of such cohorts may indicate that health phenotype differences between the sexes are sufficiently great that specific tools may be required to interrogate the complexity of different aspects of aging, either health trajectory or general biological age prediction. This is highlighted by the PhenoAge findings discussed above—that age-related CpG sites may remain similar, hence the low error metrics, but health-distinguishing sites may differ, and thus the low correlations appear.
While there are several limitations common to many studies of the oldest old, notably the small sample size, the lack of an age-matched control group, and a lack of health status data for this group of individuals that would further elucidate the findings of this study, we present suggestive data that add to the growing body of evidence that aging of the oldest old men and women cannot be analyzed adequately with extrapolations from tools designed in younger adults. Given the small number of control participants and their much younger age, it was not possible to make comparisons regarding possible sex differences within control regions and the Mediterranean Blue Zones. We were unable to measure any counteracting anti-inflammaging response in these long-lived individuals, and given that we only have DNA-predicted Treg proportion information, we cannot determine whether the difference is functionally relevant. Finally, the relatively poor performance of the clocks observed here may indicate that, as the epigenetic clocks assume a linear relationship between chronological age and epigenetic age, it may be that at the oldest old age range, the predicted relationship may not be linear making these clocks inaccurate in this age group (Snir et al., 2019). Despite the above limitations, this study has the strength of equal numbers of oldest old men and women from unique areas, the Mediterranean Blue Zones. In these regions where healthy aging, manifested as longevity, is prevalent, it is possible to investigate health span in order to garner insights into the aging process at the upper end of life.
This study analyzed the performance of increasingly popular biological aging measures, epigenetic clocks, and then compared these and other DNAm-derived health predictors in a rare population of equally long-surviving men and women from the Mediterranean Blue Zones—regions of centenarian abundance. We observed oldest old men had accelerated biological aging and shorter telomere length in comparison to oldest old women in these Blue Zones according to what we detected as the most precise predictors in this group, despite the equal longevity of the sexes in these regions. We also observed lower predicted IL-6 score in men and higher Treg proportion in women from the Mediterranean Blue Zone, which were both potential indicators of the survival advantage in men in these regions through specific immunological pathways perhaps not fully captured by DNAm-derived telomere length and epigenetic clocks. However, the relatively poor performance of all epigenetic clocks does indicate that these tools, currently trained using approximately middle-aged cohorts, performed inadequately, and epigenetic clocks with a focus on the oldest old should be developed to better assess health of the elderly and to capture sex-specific drivers of EAA. We found compelling evidence the epigenetic clocks support the female longevity phenomenon, even in regions of equal male longevity, but we also found evidence these predictors may fail to capture different aspects of the aging process amongst the oldest old. To this end, further research is required to understand the sex-dependent lifespan and health span gaps, which could be elucidated by an inflammation and immune-focused direction to interrogate the survival benefits of men in the Mediterranean Blue Zones.
Data availability statement
The datasets generated for this study are available on request to the corresponding authors.
Ethics statement
The studies involving human participants were reviewed and approved by the University of British Columbia C&W Research Ethics Board. The patients/participants provided their written informed consent to participate in this study.
Author contributions
H-RE completed data normalization, preparation and analyses, and wrote the manuscript. SM and NG were both involved in the preparation of the manuscript, and data analysis strategies. JM and DL performed the laboratory work to generate the DNA methylation data. SE provided commentary and edits to the manuscript. CC and GP contributed data from the participants. MK was involved in conceptualization of the analysis, provision of the facilities and training in which to generate the DNA methylation data and analyses, and manuscript review. DR was involved in obtaining participant data, conceptualizing the analysis, as well as manuscript review.
Funding
MK was supported by a grant from the Canadian Institutes of Health Research (EGM-141897) and is the Edwin S.H. Leong UBC Chair in Healthy Aging. DR was supported by a National Institutes of Health grant (MD011721).
Acknowledgments
The authors acknowledge and are grateful for the involvement of the participants in their respective studies. The authors are grateful for the support of Chaini Konwar and Alan Kerr for their review of this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fragi.2022.1007098/full#supplementary-material
Footnotes
1DNA Methylation Age Calculator: https://dnamage.genetics.ucla.edu/new
2DunedinPACE package: https://github.com/danbelsky/DunedinPACE
References
Ahrenfeldt, L. J., Möller, S., Thinggaard, M., Christensen, K., and Lindahl-Jacobsen, R. (2019). Sex differences in comorbidity and frailty in Europe. Int. J. Public Health 64, 1025–1036. doi:10.1007/s00038-019-01270-9
Anstey, K. J., Peters, R., Mortby, M. E., Kiely, K. M., Eramudugolla, R., Cherbuin, N., et al. (2021). Association of sex differences in dementia risk factors with sex differences in memory decline in a population-based cohort spanning 20–76 years. Sci. Rep. 11, 7710–10. doi:10.1038/s41598-021-86397-7
Arai, Y., Martin-Ruiz, C. M., Takayama, M., Abe, Y., Takebayashi, T., Koyasu, S., et al. (2015). Inflammation, but not telomere length, predicts successful ageing at extreme old age: A longitudinal study of semi-supercentenarians. EBioMedicine 2, 1549–1558. doi:10.1016/j.ebiom.2015.07.029
Austad, S. N. (2011). “Sex differences in longevity and aging,” in Handbook of the biology of aging. Handbooks of aging. Editors E. J. Masoro, and S. N. Austad. Seventh Edition (San Diego: Academic Press), 479–495. doi:10.1016/B978-0-12-378638-8.00023-3
Belsky, D. W., Caspi, A., Arseneault, L., Baccarelli, A., Corcoran, D. L., Gao, X., et al. (2020). Quantification of the pace of biological aging in humans through a blood test, the DunedinPoAm DNA methylation algorithm. eLife 9, e54870. doi:10.7554/eLife.54870
Belsky, D. W., Caspi, A., Corcoran, D. L., Sugden, K., Poulton, R., Arseneault, L., et al. (2022). Quantification of the pace of biological aging in humans through a blood test, the DunedinPoAm DNA methylation algorithm. eLife 11, e54870. doi:10.7554/eLife.54870
Bergsma, T., and Rogaeva, E. (2020). DNA methylation clocks and their predictive capacity for aging phenotypes and healthspan. Neurosci. Insights 15, 2633105520942221. doi:10.1177/2633105520942221
Bollepalli, S., Korhonen, T., Kaprio, J., Anders, S., and Ollikainen, M. (2019). EpiSmokEr: A robust classifier to determine smoking status from DNA methylation data. Epigenomics 11, 1469–1486. doi:10.2217/epi-2019-0206
Bots, S. H., Peters, S. A. E., and Woodward, M. (2017). Sex differences in coronary heart disease and stroke mortality: A global assessment of the effect of ageing between 1980 and 2010. BMJ Glob. Health 2, e000298. doi:10.1136/bmjgh-2017-000298
Breitling, L. P., Saum, K.-U., Perna, L., Schöttker, B., Holleczek, B., and Brenner, H. (2016). Frailty is associated with the epigenetic clock but not with telomere length in a German cohort. Clin. Epigenetics 8, 21. doi:10.1186/s13148-016-0186-5
Buettner, D. (2012). 9 lessons for living longer from the people who’ve lived the longest. Second Edition. Washington, DC: National Geographic.The Blue Zones
Carmel, S. (2019). Health and well-being in late life: Gender differences worldwide. Front. Med. 6, 218. doi:10.3389/fmed.2019.00218
Christopoulou, R., Mavropoulos, G., and Voucharas, G. (2021). The Greek smoking epidemic from a life-course perspective. J. Public Health, fdab342. fdab342. doi:10.1093/pubmed/fdab342
Churov, A. V., Mamashov, K. Y., and Novitskaia, A. V. (2020). Homeostasis and the functional roles of CD4+ Treg cells in aging. Immunol. Lett. 226, 83–89. doi:10.1016/j.imlet.2020.07.004
Crimmins, E. M., Thyagarajan, B., Levine, M. E., Weir, D. R., and Faul, J. (2021). Associations of age sex race/ethnicity and education with 13 epigenetic clocks in a nationally representative U.S. sample: The Health and Retirement Study. J. Gerontol. A Biol. Sci. Med. Sci. 76, 1117–1123. doi:10.1093/GERONA/GLAB016
Du, P., Kibbe, W. A., and Lin, S. M. (2008). lumi: a pipeline for processing Illumina microarray. Bioinformatics 24, 1547–1548. doi:10.1093/bioinformatics/btn224
Ecker, S., and Beck, S. (2019). The epigenetic clock: A molecular crystal ball for human aging? Aging 11, 833–835. doi:10.18632/aging.101712
Elliott, H. R., Tillin, T., McArdle, W. L., Ho, K., Duggirala, A., Frayling, T. M., et al. (2014). Differences in smoking associated DNA methylation patterns in South Asians and Europeans. Clin. Epigenetics 6, 4. doi:10.1186/1868-7083-6-4
Ershler, W. B. (1993). Interleukin-6: A cytokine for gerontologists. J. Am. Geriatr. Soc. 41, 176–181. doi:10.1111/j.1532-5415.1993.tb02054.x
Fastame, M., Hitchcott, P., Mulas, I., Ruiu, M., and Penna, M. (2018). Resilience in elders of the Sardinian Blue Zone: An explorative study. Behav. Sci. 8, 30. doi:10.3390/bs8030030
Ferrucci, L., and Fabbri, E. (2018). Inflammageing: Chronic inflammation in ageing, cardiovascular disease, and frailty. Nat. Rev. Cardiol. 15, 505–522. doi:10.1038/s41569-018-0064-2
Filzmoser, P., Hron, K., and Reimann, C. (2009). Principal component analysis for compositional data with outliers. Environmetrics 20, 621–632. doi:10.1002/env.966
Fortin, J.-P., Labbe, A., Lemire, M., Zanke, B. W., Hudson, T. J., Fertig, E. J., et al. (2014). Functional normalization of 450k methylation array data improves replication in large cancer studies. Genome Biol. 15, 503. doi:10.1186/s13059-014-0503-2
Fortin, J. P., Triche, T. J., and Hansen, K. D. (2017). Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics 33, 558–560. doi:10.1093/BIOINFORMATICS/BTW691
Foscolou, A., Chrysohoou, C., Dimitriadis, K., Masoura, K., Vogiatzi, G., Gkotzamanis, V., et al. (2021). The association of healthy aging with multimorbidity: IKARIA study. Nutrients 13, 1386. doi:10.3390/nu13041386
Frasca, D., and Blomberg, B. B. (2016). Inflammaging decreases adaptive and innate immune responses in mice and humans. Biogerontology 17, 7–19. doi:10.1007/s10522-015-9578-8
Furukawa, T., Inoue, M., Kajiya, F., Inada, H., Takasugi, S., Fukui, S., et al. (1975). Assessment of biological age by multiple regression analysis. J. Gerontol. 30, 422–434. doi:10.1093/geronj/30.4.422
Galkin, F., Mamoshina, P., Kochetov, K., Sidorenko, D., and Zhavoronkov, A. (2021). DeepMAge: A methylation aging clock developed with deep learning. Aging Dis. 12, 1252–1262. doi:10.14336/AD.2020.1202
Gallus, S., Lugo, A., Liu, X., Behrakis, P., Boffi, R., Bosetti, C., et al. (2021). Who smokes in europe? Data from 12 European countries in the TackSHS survey 2017–2018). J. Epidemiol. 31, 145–151. doi:10.2188/jea.JE20190344
Gardner, M., Bann, D., Wiley, L., Cooper, R., Hardy, R., Nitsch, D., et al. (2014). Gender and telomere length: Systematic review and meta-analysis. Exp. Gerontol. 51, 15–27. doi:10.1016/j.exger.2013.12.004
Garmany, A., Yamada, S., and Terzic, A. (2021). Longevity leap: Mind the healthspan gap. NPJ Regen. Med. 6, 57. doi:10.1038/s41536-021-00169-5
Gervin, K., Salas, L. A., Bakulski, K. M., Zelm, M. C., Koestler, D. C., Wiencke, J. K., et al. (2019). Systematic evaluation and validation of reference and library selection methods for deconvolution of cord blood DNA methylation data. Clin. Epigenetics 11, 125. doi:10.1186/s13148-019-0717-y
Goodman, W. A., Bedoyan, S. M., Havran, H. L., Richardson, B., Cameron, M. J., and Pizarro, T. T. (2020). Impaired estrogen signaling underlies regulatory T cell loss-of-function in the chronically inflamed intestine. Proc. Natl. Acad. Sci. U. S. A. 117, 17166–17176. doi:10.1073/pnas.2002266117
Goodman, W. A., Garg, R. R., Reuter, B. K., Mattioli, B., Rissman, E. F., and Pizarro, T. T. (2014). Loss of estrogen-mediated immunoprotection underlies female gender bias in experimental Crohn’s-like ileitis. Mucosal Immunol. 7, 1255–1265. doi:10.1038/mi.2014.15
Grant, O. A., Wang, Y., Kumari, M., Zabet, N. R., and Schalkwyk, L. (2022). Characterising sex differences of autosomal DNA methylation in whole blood using the Illumina EPIC array. Clin. Epigenetics 141 14, 62–16. doi:10.1186/S13148-022-01279-7
Hägg, S., and Jylhävä, J. (2021). Sex differences in biological aging with a focus on human studies. eLife 10, e63425–e63427. doi:10.7554/eLife.63425
Hannum, G., Guinney, J., Zhao, L., Zhang, L., Hughes, G., Sadda, S., et al. (2013). Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol. Cell 49, 359–367. doi:10.1016/j.molcel.2012.10.016
Hillary, R. F., Stevenson, A. J., McCartney, D. L., Campbell, A., Walker, R. M., Howard, D. M., et al. (2020). Epigenetic measures of ageing predict the prevalence and incidence of leading causes of death and disease burden. Clin. Epigenetics 12, 115. doi:10.1186/S13148-020-00905-6
Hitchcott, P. K., Fastame, M. C., and Penna, M. P. (2018). More to Blue Zones than long life: Positive psychological characteristics. Health Risk Soc. 20, 163–181. doi:10.1080/13698575.2018.1496233
Horvath, S. (2013). DNA methylation age of human tissues and cell types. Genome Biol. 14, 115. doi:10.1186/gb-2013-14-10-r115
Horvath, S., Gurven, M., Levine, M. E., Trumble, B. C., Kaplan, H., Allayee, H., et al. (2016). An epigenetic clock analysis of race/ethnicity, sex, and coronary heart disease. Genome Biol. 17, 171. doi:10.1186/s13059-016-1030-0
Horvath, S., Oshima, J., Martin, G. M., Lu, A. T., Quach, A., Cohen, H., et al. (2018). Epigenetic clock for skin and blood cells applied to Hutchinson Gilford Progeria Syndrome and ex vivo studies. Aging 10, 1758–1775. doi:10.18632/aging.101508
Ji, L., Jazwinski, S. M., and Kim, S. (2021). Frailty and biological age. Ann. Geriatr. Med. Res. 25, 141–149. doi:10.4235/agmr.21.0080
Jones, M. J., Goodman, S. J., and Kobor, M. S. (2015). DNA methylation and healthy human aging. Aging Cell 14, 924–932. doi:10.1111/acel.12349
Jylhävä, J., Pedersen, N. L., and Hägg, S. (2017). Biological age predictors. eBioMedicine 21, 29–36. doi:10.1016/J.EBIOM.2017.03.046
Kankaanpää, A., Tolvanen, A., Saikkonen, P., Heikkinen, A., Laakkonen, E. K., Kaprio, J., et al. (2021). Do epigenetic clocks provide explanations for sex differences in life span? A cross-sectional twin study. J. Gerontol. A Biol. Sci. Med. Sci. 77, 1898–1906. doi:10.1093/gerona/glab337
Kaur, S., Bansal, Y., Kumar, R., and Bansal, G. (2020). A panoramic review of IL-6: Structure, pathophysiological roles and inhibitors. Bioorg. Med. Chem. 28, 115327. doi:10.1016/j.bmc.2020.115327
Klein, S. L., and Flanagan, K. L. (2016). Sex differences in immune responses. Nat. Rev. Immunol. 16, 626–638. doi:10.1038/nri.2016.90
Klopack, E. T., Carroll, J. E., Cole, S. W., Seeman, T. E., and Crimmins, E. M. (2022). Lifetime exposure to smoking, epigenetic aging, and morbidity and mortality in older adults. Clin. Epigenetics 14, 72. doi:10.1186/s13148-022-01286-8
Koestler, D. C., Jones, M. J., Usset, J., Christensen, B. C., Butler, R. A., Kobor, M. S., et al. (2016). Improving cell mixture deconvolution by identifying optimal DNA methylation libraries (IDOL). BMC Bioinforma. 17, 120. doi:10.1186/s12859-016-0943-7
Koo, H.-K., Morrow, J., Kachroo, P., Tantisira, K., Weiss, S. T., Hersh, C. P., et al. (2021). Sex-specific associations with DNA methylation in lung tissue demonstrate smoking interactions. Epigenetics 16, 692–703. doi:10.1080/15592294.2020.1819662
Lang, M., Binder, M., Richter, J., Schratz, P., Pfisterer, F., Coors, S., et al. (2019). mlr3: A modern object-oriented machine learning framework in R. J. Open Source Softw. 4, 1903. doi:10.21105/joss.01903
Leek, J. T., Johnson, W. E., Parker, H. S., Fertig, E. J., Jaffe, A. E., Storey, J. D., et al. (2019). sva: surrogate variable analysis.
Leek, J. T., Johnson, W. E., Parker, H. S., Jaffe, A. E., and Storey, J. D. (2012). The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics 28 (6), 882–883. doi:10.1093/bioinformatics/bts034
Levine, D. A., Gross, A. L., Briceño, E. M., Tilton, N., Giordani, B. J., Sussman, J. B., et al. (2021). Sex differences in cognitive decline among US adults. JAMA Netw. Open 4, 210169. doi:10.1001/JAMANETWORKOPEN.2021.0169
Levine, M. E., Lu, A. T., Quach, A., Chen, B. H., Assimes, T. L., Bandinelli, S., et al. (2018). An epigenetic biomarker of aging for lifespan and healthspan. Aging 10, 573–591. doi:10.18632/aging.101414
Lu, A. T., Quach, A., Wilson, J. G., Reiner, A. P., Aviv, A., Raj, K., et al. (2019a). DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging 11, 303–327. doi:10.18632/aging.101684
Lu, A. T., Seeboth, A., Tsai, P. C., Sun, D., Quach, A., Reiner, A.P., et al. (2019b). DNA methylation-based estimator of telomere length. Aging (Albany. NY) 11, 5895–5923. doi:10.18632/AGING.102173
Lugo, A., Vecchia, C., Boccia, S., Murisic, B., and Gallus, S. (2013). Patterns of smoking prevalence among the elderly in Europe. Int. J. Environ. Res. Public Health 10, 4418–4431. doi:10.3390/ijerph10094418
Maas, S. C. E., Mens, M. M. J., Kühnel, B., van Meurs, J. B. J., Uitterlinden, A. G., Peters, A., et al. (2020). Smoking-related changes in DNA methylation and gene expression are associated with cardio-metabolic traits. Clin. Epigenetics 12, 157. doi:10.1186/s13148-020-00951-0
Maddock, J., Castillo-Fernandez, J., Wong, A., Cooper, R., Richards, M., Ong, K. K., et al. (2019). DNA methylation age and physical and cognitive ageing. J. Gerontol. A Biol. Sci. Med. Sci. glz246. doi:10.1093/gerona/glz246
Marioni, R. E., Harris, S. E., Shah, S., McRae, A. F., Zglinicki, T., Martin-Ruiz, C., et al. (2016). The epigenetic clock and telomere length are independently associated with chronological age and mortality. Int. J. Epidemiol. 45, 424–432. doi:10.1093/IJE/DYW041
Marioni, R. E., Shah, S., McRae, A. F., Ritchie, S. J., Muniz-Terrera, G., Harris, S. E., et al. (2015). The epigenetic clock is correlated with physical and cognitive fitness in the Lothian Birth Cohort 1936. Int. J. Epidemiol. 44, 1388–1396. doi:10.1093/ije/dyu277
McCrory, C., Fiorito, G., Hernandez, B., Polidoro, S., O’Halloran, A. M., Hever, A., et al. (2021). GrimAge outperforms other epigenetic clocks in the prediction of age-related clinical phenotypes and all-cause mortality. J. Gerontol. A Biol. Sci. Med. Sci. 76, 741–749. doi:10.1093/gerona/glaa286
McCrory, C., Fiorito, G., McLoughlin, S., Polidoro, S., Cheallaigh, C. N., Bourke, N., et al. (2020). Epigenetic clocks and allostatic load reveal potential sex-specific drivers of biological aging. J. Gerontol. A Biol. Sci. Med. Sci. 75, 495–503. doi:10.1093/gerona/glz241
McEwen, L. M., Morin, A. M., Edgar, R. D., MacIsaac, J. L., Jones, M. J., Dow, W. H., et al. (2017). Differential DNA methylation and lymphocyte proportions in a Costa Rican high longevity region. Epigenetics Chromatin 10, 21. doi:10.1186/s13072-017-0128-2
Meyer, A. C., Drefahl, S., Ahlbom, A., Lambe, M., and Modig, K. (2020). Trends in life expectancy: Did the gap between the healthy and the ill widen or close? BMC Med. 18, 41. doi:10.1186/s12916-020-01514-z
Minciullo, P. L., Catalano, A., Mandraffino, G., Casciaro, M., Crucitti, A., Maltese, G., et al. (2016). Inflammaging and anti-inflammaging: The role of cytokines in extreme longevity. Arch. Immunol. Ther. Exp. 64, 111–126. doi:10.1007/s00005-015-0377-3
Müller, L., Di Benedetto, S., and Pawelec, G. (2019). “The immune system and its dysregulation with aging,” in Biochemistry and cell biology of ageing: Part II clinical Science Subcellular biochemistry. Editors J. R. Harris, and V. I. Korolchuk (Singapore: Springer), 21–43. doi:10.1007/978-981-13-3681-2_2
Niccoli, T., and Partridge, L. (2012). Ageing as a risk factor for disease. Curr. Biol. 22, 741–752. doi:10.1016/J.CUB.2012.07.024
Oksuzyan, A., Juel, K., Vaupel, J. W., and Christensen, K. (2008). Men: Good health and high mortality. Sex differences in health and aging. Aging Clin. Exp. Res. 20, 91–102. doi:10.1007/BF03324754
Öngel, M. E., Yıldız, C., Akpınaroğlu, C., Yilmaz, B., and Özilgen, M. (2021). Why women may live longer than men do? A telomere-length regulated and diet-based entropic assessment. Clin. Nutr. 40, 1186–1191. doi:10.1016/j.clnu.2020.07.030
Papanicolaou, D. A., and Vgontzas, A. N. (2000). Interleukin-6: The endocrine cytokine. J. Clin. Endocrinol. Metab. 85, 1331–1333. doi:10.1210/jcem.85.3.6582
Passarino, G., Calignano, C., Vallone, A., Franceschi, C., Jeune, B., Robine, J., et al. (2002). Male/female ratio in centenarians: A possible role played by population genetic structure. Exp. Gerontol. 37, 1283–1289. doi:10.1016/S0531-5565(02)00140-7
Perna, L., Zhang, Y., Mons, U., Holleczek, B., Saum, K.-U., and Brenner, H. (2016). Epigenetic age acceleration predicts cancer, cardiovascular, and all-cause mortality in a German case cohort. Clin. Epigenetics 8, 64. doi:10.1186/s13148-016-0228-z
Pes, G. M., Dore, M. P., Errigo, A., and Poulain, M. (2018). Analysis of physical activity among free–living nonagenarians from a Sardinian longevous population. J. Aging Phys. Act. 26, 254–258. doi:10.1123/japa.2017-0088
Pes, G. M., Dore, M. P., Tsofliou, F., and Poulain, M. (2022). Diet and longevity in the Blue Zones: A set-and-forget issue? Maturitas 164, 31–37. doi:10.1016/j.maturitas.2022.06.004
Pes, G. M., Errigo, A., Tedde, P., and Dore, M. P. (2020). Sociodemographic, clinical and functional profile of nonagenarians from two areas of Sardinia characterized by distinct longevity levels. Rejuvenation Res. 23, 341–348. doi:10.1089/rej.2018.2129
Pidsley, R., Wong, C. Y., Volta, M., Lunnon, K., Mill, J., and Schalkwyk, L. C. (2013). A data-driven approach to preprocessing Illumina 450K methylation array data. BMC Genomics 14, 293. doi:10.1186/1471-2164-14-293
Pidsley, R., Zotenko, E., Peters, T. J., Lawrence, M. G., Risbridger, G. P., Molloy, P., et al. (2016). Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biol. 17, 208. doi:10.1186/s13059-016-1066-1
Poulain, M., Pes, G., and Salaris, L. (2011). A population where men live as long as women: Villagrande Strisaili, Sardinia. J. Aging Res. 2011, 153756. doi:10.4061/2011/153756
Poulain, M., Herm, A., and Pes, G. (2013). The Blue Zones: Areas of exceptional longevity around the world. Vienna Yearb. Popul. Res. Vol. 11, 87–108. doi:10.1553/populationyearbook2013s87
Poulain, M., Herm, A., Errigo, A., Chrysohoou, C., Legrand, R., Passarino, G., et al. (2021). Specific features of the oldest old from the longevity Blue Zones in Ikaria and Sardinia. Mech. Ageing Dev. 198, 111543. doi:10.1016/J.MAD.2021.111543
Poulain, M., and Mackowicz, J. (Editors) (2021). Positive ageing and learning from centenarians: Living longer and better (London: Routledge). doi:10.4324/9781003162216
Price, E. M., Cotton, A. M., Lam, L. L., Farré, P., Emberly, E., Brown, C. J., et al. (2013). Additional annotation enhances potential for biologically-relevant analysis of the Illumina Infinium HumanMethylation450 BeadChip array. Epigenetics Chromatin 6, 4. doi:10.1186/1756-8935-6-4
Puca, A. A., Spinelli, C., Accardi, G., Villa, F., and Caruso, C. (2018). Centenarians as a model to discover genetic and epigenetic signatures of healthy ageing. Mech. Ageing Dev. 174, 95–102. doi:10.1016/j.mad.2017.10.004
Pusceddu, I., Kleber, M., Delgado, G., Herrmann, W., März, W., and Herrmann, M. (2018). Telomere length and mortality in the ludwigshafen risk and cardiovascular health study. PLoS One 13, e0198373. doi:10.1371/JOURNAL.PONE.0198373
Puzianowska-Kuźnicka, M., Owczarz, M., Wieczorowska-Tobis, K., Nadrowski, P., Chudek, J., Slusarczyk, P., et al. (2016). Interleukin-6 and C-reactive protein, successful aging, and mortality: The PolSenior study. Immun. Ageing 13, 21. doi:10.1186/s12979-016-0076-x
Quach, A., Levine, M. E., Tanaka, T., Lu, A. T., Chen, B. H., Ferrucci, L., et al. (2017). Epigenetic clock analysis of diet, exercise, education, and lifestyle factors. Aging 9, 419–446. doi:10.18632/aging.101168
Raynor, J., Lages, C. S., Shehata, H., Hildeman, D. A., and Chougnet, C. A. (2012). Homeostasis and function of regulatory T cells in aging. Curr. Opin. Immunol. 24, 482–487. doi:10.1016/j.coi.2012.04.005
Rehkopf, D. H., Dow, W. H., Rosero-Bixby, L., Lin, J., Epel, E. S., and Blackburn, E. H. (2013). Longer leukocyte telomere length in Costa Rica’s Nicoya Peninsula: A population-based study. Exp. Gerontol. 48, 1266–1273. doi:10.1016/J.EXGER.2013.08.005
Rocamora-Reverte, L., Melzer, F. L., Würzner, R., and Weinberger, B. (2021). The complex role of regulatory T cells in immunity and aging. Front. Immunol. 11, 616949. doi:10.3389/fimmu.2020.616949
Rockwood, K., Song, X., MacKnight, C., Bergman, H., Hogan, D. B., McDowell, I., et al. (2005). A global clinical measure of fitness and frailty in elderly people. C. Can. Med. Assoc. J. 173, 489–495. doi:10.1503/CMAJ.050051
Roetker, N. S., Pankow, J. S., Bressler, J., Morrison, A. C., and Boerwinkle, E. (2018). Prospective study of epigenetic age acceleration and incidence of cardiovascular disease outcomes in the ARIC Study Atherosclerosis Risk in Communities). Circ. Genom. Precis. Med. 11, e001937. doi:10.1161/CIRCGEN.117.001937
Salas, L. A., Zhang, Z., Koestler, D. C., Butler, R. A., Hansen, H. M., Molinaro, A. M., et al. (2022). Enhanced cell deconvolution of peripheral blood using DNA methylation for high-resolution immune profiling. Nat. Commun. 13, 761. doi:10.1038/s41467-021-27864-7
Seale, K., Horvath, S., Teschendorff, A., Eynon, N., and Voisin, S. (2022). Making sense of the ageing methylome. Nat. Rev. Genet. 23, 585–605. doi:10.1038/s41576-022-00477-6
Shu, Y., Hu, Q., Long, H., Chang, C., Lu, Q., and Xiao, R. (2017). Epigenetic variability of CD4+CD25+ Tregs Contributes to the pathogenesis of autoimmune diseases. Clin. Rev. Allergy Immunol. 52, 260–272. doi:10.1007/s12016-016-8590-3
Sindi, S., Kåreholt, I., Ngandu, T., Rosenberg, A., Kulmala, J., Johansson, L., et al. (2021). Sex differences in dementia and response to a lifestyle intervention: Evidence from Nordic population-based studies and a prevention trial. Alzheimer’s Dement. 17, 1166–1178. doi:10.1002/ALZ.12279
Snir, S., Farrell, C., and Pellegrini, M. (2019). Human epigenetic ageing is logarithmic with time across the entire lifespan. Epigenetics 14, 912–926. doi:10.1080/15592294.2019.1623634
Soler-Botija, C., Gálvez-Montón, C., and Bayés-Genís, A. (2019). Epigenetic biomarkers in cardiovascular diseases. Front. Genet. 10, 950. doi:10.3389/fgene.2019.00950
Sotgia, S., Zinellu, A., Mangoni, A. A., Serra, R., Pintus, G., Caruso, C., et al. (2017). Cellular immune activation in Sardinian middle-aged, older adults and centenarians. Exp. Gerontol. 99, 133–137. doi:10.1016/j.exger.2017.10.005
Stevenson, A. J., Gadd, D. A., Hillary, R. F., McCartney, D. L., Campbell, A., Walker, R. M., et al. (2021). Creating and validating a DNA methylation-based proxy for Interleukin-6. J. Gerontol. A Biol. Sci. Med. Sci. 76, 2284–2292. doi:10.1093/gerona/glab046
Sugden, K., Caspi, A., Elliott, M. L., Bourassa, K. J., Chamarti, K., Corcoran, D. L., et al. (2022). Association of pace of aging measured by blood-based DNA methylation with age-related cognitive impairment and dementia. Neurology 99, e1402–e1413. doi:10.1212/WNL.000000000020089810.1212/WNL.0000000000200898
Tamosiuniene, R., Manouvakhova, O., Mesange, P., Saito, T., Qian, J., Sanyal, M., et al. (2018). Dominant role for regulatory T cells in protecting females against pulmonary hypertension. Circ. Res. 122, 1689–1702. doi:10.1161/CIRCRESAHA.117.312058
Thinggaard, M., McGue, M., Jeune, B., Osler, M., Vaupel, J. W., and Christensen, K. (2016). Survival prognosis in very old adults. J. Am. Geriatr. Soc. 64, 81–88. doi:10.1111/JGS.13838
Tucker, L. A. (2019). Serum and dietary folate and vitamin B12 levels account for differences in cellular aging: Evidence based on telomere findings in 5581 U.S. Adults. Oxid. Med. Cell. Longev. 2019, 4358717. doi:10.1155/2019/4358717
Vaccarino, V., Huang, M., Wang, Z., Hui, Q., Shah, A. J., Goldberg, J., et al. (2021). Epigenetic age acceleration and cognitive decline: A twin study. J. Gerontol. A Biol. Sci. Med. Sci. 76, 1854–1863. doi:10.1093/gerona/glab047
Vaiserman, A., and Krasnienkov, D. (2021). Telomere length as a marker of biological age: State-of-the-art open issues and future perspectives. Front. Genet. 11, 630186. doi:10.3389/fgene.2020.630186
Vasiljevic, Z., Scarpone, M., Bergami, M., Yoon, J., van der Schaar, M., Krljanac, G., et al. (2021). Smoking and sex differences in first manifestation of cardiovascular disease. Atherosclerosis 330, 43–51. doi:10.1016/J.ATHEROSCLEROSIS.2021.06.909
Verschoor, C. P., Lin, D. T. S., Kobor, M. S., Mian, O., Ma, J., Pare, G., et al. (2021). Epigenetic age is associated with baseline and 3-year change in frailty in the Canadian Longitudinal Study on Aging. Clin. Epigenetics 13, 163. doi:10.1186/s13148-021-01150-1
Keywords: blue zone, centenarian, epigenetic clock, biological age, telomere, biomarker, DNA methylation
Citation: Engelbrecht H-R, Merrill SM, Gladish N, MacIsaac JL, Lin DTS, Ecker S, Chrysohoou CA, Pes GM, Kobor MS and Rehkopf DH (2022) Sex differences in epigenetic age in Mediterranean high longevity regions. Front. Aging 3:1007098. doi: 10.3389/fragi.2022.1007098
Received: 29 July 2022; Accepted: 21 October 2022;
Published: 23 November 2022.
Edited by:
Victoria Belancio, Tulane University, United StatesReviewed by:
Yasumichi Arai, Keio University School of Medicine, JapanCalogero Caruso, University of Palermo, Italy
Copyright © 2022 Engelbrecht, Merrill, Gladish, MacIsaac, Lin, Ecker, Chrysohoou, Pes, Kobor and Rehkopf. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Michael S. Kobor, msk@bcchr.ca; David H. Rehkopf, drehkopf@stanford.edu
†These authors share senior authorship