Blood and skeletal muscle ageing determined by epigenetic clocks and their associations with physical activity and functioning
Clinical Epigenetics volume 13, Article number: 110 (2021)
The aim of this study was to investigate the correspondence of different biological ageing estimates (i.e. epigenetic age) in blood and muscle tissue and their associations with physical activity (PA), physical function and body composition. Two independent cohorts (N = 139 and N = 47) were included, whose age span covered adulthood (23–69 years). Whole blood and m. vastus lateralis samples were collected, and DNA methylation was analysed. Four different DNA methylation age (DNAmAge) estimates were calculated using genome-wide methylation data and publicly available online tools. A novel muscle-specific methylation age was estimated using the R-package ‘MEAT’. PA was measured with questionnaires and accelerometers. Several tests were conducted to estimate cardiorespiratory fitness and muscle strength. Body composition was estimated by dual-energy X-ray absorptiometry. DNAmAge estimates from blood and muscle were highly correlated with chronological age, but different age acceleration estimates were weakly associated with each other. The monozygotic twin within-pair similarity of ageing pace was higher in blood (r = 0.617–0.824) than in muscle (r = 0.523–0.585). Associations of age acceleration estimates with PA, physical function and body composition were weak in both tissues and mostly explained by smoking and sex. The muscle-specific epigenetic clock MEAT was developed to predict chronological age, which may explain why it did not associate with functional phenotypes. The Horvath’s clock and GrimAge were weakly associated with PA and related phenotypes, suggesting that higher PA would be linked to accelerated biological ageing in muscle. This may, however, be more reflective of the low capacity of epigenetic clock algorithms to measure functional muscle ageing than of actual age acceleration. Based on our results, the investigated epigenetic clocks have rather low value in estimating muscle ageing with respect to the physiological adaptations that typically occur due to ageing or PA. Thus, further development of methods is needed to gain insight into muscle tissue-specific ageing and the underlying biological pathways.
Human ageing is accompanied by a decline in physical activity (PA), partly resulting from a decline in muscle mass and function commonly referred to as sarcopenia . These age-related changes result in decreased mobility, loss of physical independence and increased morbidity and mortality. The precise underlying molecular mechanisms behind these physiological changes remain poorly characterised. It is clear that individual’s inherited genetic variants, epigenetic makeup, and the resultant regulation of gene expression, impacts of PA and PA-related phenotypes as well as the ageing process [2,3,4,5].
Epigenetic modifications are crucial hallmarks of ageing . As an interface between genome and external factors, epigenetics forms a fundamental link between genotype and environment in ageing. DNA methylation is the most widely studied epigenetic mechanism. It may modify gene function without changing the DNA sequence, in part by reacting to the internal and external environment. Some of the alterations in DNA methylation at specific CpG (a cytosine nucleotide followed by a guanine nucleotide) sites correlate with chronological age. This finding, together with technological innovations in the field of microarrays, development of infrastructures for data sharing and adaptation of statistical modelling methods for biological ageing studies, has offered promising tools for identifying epigenetic markers of biological age [7,8,9,10]. The novel epigenetic biomarkers, referred to as epigenetic clocks, are weighted sums of DNA methylation at specific CpG sites that can be used as indicators of biological age. Epigenetic clocks are based on penalty regression models in which chronological age or other ageing phenotypes are predicted using genome-wide DNA methylation data, and optimal CpG sites are selected for the age predictor by the algorithm .
First-generation clocks, Horvath’s and Hannum’s, were developed to predict chronological age [7, 8]. The epigenetic age estimate, DNA methylation age (DNAmAge), predicts differences in mortality better than chronological age after adjusting for several known risk factors . However, these first clocks were criticised for excluding CpGs whose methylation patterns may reflect deviation of epigenetic age from chronological age . Therefore, the next-generation epigenetic clocks were developed to capture CpGs with their methylation associated with the functional stage and chronological age of an individual. PhenoAge was developed using ‘phenotypic age measure’ instead of chronological age as a reference , and GrimAge was built using a combination of DNA methylation-based surrogate biomarkers for mortality, such as methylation-based estimates of health-related plasma proteins and smoking pack-years . These two novel epigenetic age estimates predict mortality risk better than Horvath’s and Hannum’s DNAmAge, while the latter ones perform better in predicting chronological age in large datasets [9, 10, 12].
The epigenetic clocks are developed using blood-based leukocyte DNA, but they seem to work in several tissues [7, 9]. There is some evidence suggesting that different tissues and organs in the same individual may vary in terms of their DNAmAge. The cerebellum shows slower ageing compared to other tissues in supercentenarians , while ageing of the liver, rather than blood, skeletal muscle or fat tissue, is accelerated in persons with obesity . These observations lend support to the assumption that tissue-specific analyses might yield further insights into age-related pathologies. Muscle is the key tissue for physical functioning and whole-body metabolism, affecting age-related decline and the onset or severity of several age-related diseases, such as type II diabetes and cardiovascular diseases and even mortality [15,16,17,18]. Understanding of DNAmAge in muscle is sparse. However, it could be a potential biomarker of the individual ageing process and may offer insights into the biological mechanisms that regulate muscle ageing.
Potential discrepancies in DNAmAge between tissues could also offer explanations for the relevance of different environmental exposures and their influence on health and morbidity. PA is a modifiable behaviour that has the potential to slow down the rate of cellular and molecular damage accumulation  and blunt the decline in physiological function with increasing age. The causal molecular and genetic links between PA and reduced age-related disease risk remain elusive, but it is clear that the effects of PA depend on its mode, intensity, duration and frequency and are tissue specific [20, 21].
Very recently, the first epigenetic clock for muscle tissue, ‘MEAT’ (Muscle Epigenetic Age Test), was developed to predict ageing in muscle . It has been suggested that this first muscle-specific epigenetic clock can predict chronological age with better accuracy than, for example, Horvath’s DNAmAge clock, which has been previously shown to have reasonable accuracy in many tissues . To the best of our knowledge, MEAT has not yet been investigated in relation to PA or physical function.
The aim of this study was to compare the biological ageing of blood and muscle tissue determined by several epigenetic clocks in two independent cohorts. In addition, we investigated whether ageing in these two tissues is associated with PA and physical function-related phenotypes.
The characteristics of the study participants, consisting of groups of younger [mean age 32.9 (SD 5.3) years] adults (Finnish twin cohort [FTC] twins), middle-aged [mean 51.8 (SD 2.0) years)] adults (Estrogenic Regulation of Muscle Apoptosis study [ERMA]) and older [mean age 63.7 (SD 4.2) years] adults (FTC twins) are presented in Table 1. The ERMA cohorts were all women, while 43% of the younger twins and 61% of the older twins were women. The younger and middle-aged adults were better educated, as 40–48% of them had a tertiary-level degree, while in older adults 20% had a tertiary-level degree. In all groups, alcohol use was rather low, and most were light drinkers. A quarter of the younger and 13% of the older twins were current smokers, while only 6% of the middle-aged adults were current smokers, which is due to the fact that smoking is more common among men.
All three groups were overweight on average (Table 1). The percentage of fat increased with age, being highest among the older FTC participants. The ERMA participants had the lowest amount of lean mass, which is logical, as this cohort included only women, who tend to be smaller than men.
Regarding PA, the sport and leisure indexes were lower and the work index was higher in older FTC participants than in younger ones. On average, the ERMA participants had daily 56 min of moderate-to-vigorous PA (MVPA) and 9400 steps per day, indicating that they were quite active. Maximal oxygen uptake was 35.2 ml/kg/min in younger FTC adults. The walking test, grip and knee extension strength and vertical jumping height were only measured in the ERMA (mean values are presented in Table 1).
Accuracy of different epigenetic clocks in estimating chronological age
In blood, all four clocks showed high correlations between the DNAmAge estimates and chronological age (r = 0.904–0.954, all p < 0.001, Fig. 1). The mean deviation between chronological age and predicted age was highest for Hannum’s DNAmAge (7.6 years) and PhenoAge (9.3 years). Lower deviations were observed for Horvath’s DNAmAge and GrimAge (5.6 and 3.6 years, respectively). Horvath’s clock and GrimAge tended to give lower age estimates to older participants than to younger ones, whereas the Hannum’s clock and PhenoAge DNAmAge estimates were consistently lower than the actual chronological age across all age groups (Fig. 1).
In muscle, DNAmAge estimated with MEAT showed a higher correlation with chronological age (r = 0.962) than DNAmAge estimated with Horvath’s pan tissue clock (r = 0.696, Fig. 2). However, GrimAge, although not developed for muscle tissue, gave age estimates that had an even slightly higher correlation with chronological age (r = 0.979). MEAT age estimates showed the highest accuracy (mean deviation 3.2 years). GrimAge age estimates differed from chronological age by 5.0 years, and Horvath’s age estimates differed by 9.5 years. In younger individuals, all clocks tended to give higher muscle age estimates compared with chronological age (Fig. 2).
Clock and tissue comparison in DNAmAge and age acceleration
In blood, the DNAmAge correlation coefficients were statistically significant between all clocks (r = 0.862–0.943) (Table 2). Blood DNAmAge estimates also correlated with all the estimates in muscle (r = 0.624–0.939). Similarly, DNAmAge estimates in muscle were significantly correlated between all clocks; a highest correlation was observed between MEAT and GrimAge (r = 0.940), while the correlations between Horvath’s clock and GrimAge and MEAT were lower (r = 0.651 and 0.714, respectively).
All DNAmAge acceleration estimates were also correlated in blood, but the associations were weak to modest (r = 0.244–0.566, Table 3). In muscle, the strongest correlation was found between Horvath’s and MEAT DNAmAge acceleration estimates (r = 0.365). Unlike the DNAmAge results, the blood DNAmAge acceleration values were at most very weakly (all r < 0.15) correlated with muscle age acceleration estimates.
Within-pair analysis of age acceleration for monozygotic (MZ) pairs
In blood, pairwise correlations were highest in Horvath’s (r = 0.807) and Hannum’s age acceleration (r = 0.824, Additional file 1, Fig. 1). Slightly lower within-pair correlations were observed with PhenoAge and GrimAge acceleration (r = 0.774 and r = 0.617, respectively). In muscle tissue, within-pair analysis showed similar moderate correlations for all three clocks (r = 0.523–0.585, Additional file 1, Fig. 2).
Associations between epigenetic clocks and physical activity
Associations between DNAmAge acceleration estimates and PA were analysed using both self-reported (FTC and ERMA) and accelerometer-measured activity (ERMA).
In blood, DNAmAge acceleration estimates had minor or non-existent associations with PA (Additional file 2, Table 1). In FTC, a weak association of borderline significance was observed between slower age acceleration GrimAge and a higher leisure index. In ERMA, slower age acceleration estimates in GrimAge were evident in the medium-activity group compared to the low-activity participants, and this association attenuated slightly after further adjustments for sex and smoking.
In muscle tissue, the muscle-specific MEAT clock was not associated with PA. Higher DNAmAge acceleration estimated by GrimAge was associated with a higher amount of MVPA and daily steps (Table 4), with consistent point estimates (although with larger confidence intervals from the two other clocks).
In addition to a separate analysis of the two cohorts, we conducted an analysis using both cohorts in the same model, aiming to study associations between DNAmAge estimates and moderate to high-intensity PA. The ERMA and FTC participants were classified into low-, medium- and high-activity groups based on tertials in the variables MVPA (min) and sport index, respectively. This analysis revealed accelerated ageing among high-activity participants compared to low-activity participants in muscle tissue with Horvath’s clock after adjusting for sex and smoking (Table 4).
Associations between epigenetic clocks and body composition
In the next chapters, for consistency, we only report results with regard to the adjusted models. Associations between blood and muscle DNAm age acceleration estimates and BMI and body composition (fat mass, lean mass and percentage of fat estimated with DXA) were investigated by combining the FTC and ERMA cohorts into the same models.
Associations between muscle DNAmAge acceleration estimates and BMI were not seen. Further, a higher Horvath’s and GrimAge acceleration in muscle was associated in the same direction with lower fat mass and fat percent, even after covariate adjustments (Table 5). Higher Horvaths’ age acceleration and MEAT age acceleration were not associated with body composition in adjusted models.
Associations between epigenetic clocks and physical functioning
Associations between blood and muscle DNAmAge acceleration estimates and maximal oxygen uptake were analysed using the younger FTC cohort. No associations were seen in blood or muscle after adjusting for sex and smoking (Table 4).
Other analyses of physical function (walking distance, muscle strength and power) were conducted in the ERMA cohort. Age acceleration estimates in blood and in muscle tissue were not associated with physical function phenotypes (Additional file 1, Table 2).
Because our study sample included 50 MZ and 4 DZ twin pairs with discordant BMI, it was important to determine whether the discordance would have affected the observed results. Therefore, we conducted several sensitivity analyses, excluding either the leaner or heavier twin in a pair from the analysis (Additional file 3). In general, the associations were not dependent on whether the analysis was conducted without the heavier or leaner twin in a pair. However, it must be noted that, after excluding the heavier twins, higher MEAT age acceleration was associated with lower maximal oxygen consumption (r = -0.366, p = 0.018) compared to a weaker association (r = -0.288) when excluding the leaner twin. In further linear regression analyses with both twins in the model, the association remained significant after adjusting for twin dependency, gender and smoking (β = -0.102, SE 0.046, p = 0.034).
Muscle tissue is important with respect to human health, well-being and quality of life. Its proper function is essential for locomotion and other important tasks related to energy metabolism, temperature regulation and recovery from adverse health events. Ageing challenges the optimal functioning of muscle tissue. Typical age-related changes include accelerated decline in muscle function and quantity, which results in problems in functional capacity and health . A method for estimating and predicting muscle ageing would be helpful in identifying persons who are at risk for age-related decline in functional capacity or morbidities, and it could also furnish mechanisms that are potential targets for treatments and interventions to improve muscle functioning among ageing populations. We investigated novel genome-wide epigenetic scores, developed to predict age, age span and healthy ageing, in blood and muscle tissue [7,8,9,10, 22]. Our results showed that these clocks are mainly suitable for predicting chronological age from both tissues, but variation in age estimates was not systematically associated with the variation in physical functional phenotypes, which suggests that the clocks lack value for research and health promotion purposes.
We found a strong linear relationship between chronological age and DNA methylation age estimates (i.e. Hannum’s clock, Horvath’s clock, GrimAge, PhenoAge and MEAT) in both blood and skeletal muscle tissues. This was expected, as most of these clocks have been developed to predict chronological age, and age is one component of the GrimAge algorithm. We also replicated the finding regarding the poorer calibration of Horvath’s clock in age estimation in muscle [7, 22]. The good overall performance of these clocks in predicting chronological age has been demonstrated repeatedly in previous studies [4, 10, 22, 24]. However, the age acceleration of blood was not associated with the age acceleration in muscle by any of the clocks, suggesting different ticking rates in blood and muscle. We also found age acceleration estimates of blood to be more similar within twin pairs than age acceleration estimates of muscle, suggesting either a larger familial component in blood-based age acceleration estimates or larger heterogeneity in muscle tissue composition. Because we only analysed MZ pairs, we cannot distinguish genetics from environmental effects (e.g. similar occupation or PA) that the twins share in the familial component estimate.
Regarding our hypothesis that PA will slow down age acceleration in muscle, we found that the age acceleration estimates from Horvath’s clock and GrimAge showed a minor but statistically significant positive association with PA measures. These results suggest that higher PA would be linked to accelerated biological ageing in muscle. As PA and ageing typically result in opposite structural and functional changes in muscle, it seems unlikely that PA would accelerate ageing processes in the muscle. Instead, we believe that the findings reflect the low capacity of the epigenetic clock algorithms to measure functional muscle ageing. Possible explanations for this include cell heterogeneity and differences in methylation profiles and structure between white blood cells and muscle tissue. Blood cells represent a constantly dividing cell population, while adult muscle tissue is more heterogeneous. It mainly includes multinucleated post-mitotic myofibers that do not divide, but also a smaller population of quiescent muscle stem cells, called satellite cells, which can re-enter the cell cycle to initiate proliferation and donate nuclei to muscle fibres. Furthermore, a small contribution to the muscle tissue DNA pool is provided by a network of capillaries formed by endothelial and epithelial cell layers as well as a variable, small amount of infiltrated leucocytes and fat and connective tissue cells. PA stimulates the formation of satellite cells, which are essential for muscle regeneration and are responsible for muscle growth . Skeletal muscle cells and satellite cells are differentially methylated, and PA may predispose them to different epigenetic changes depending on age . Satellite cells are important in the regeneration of muscle after cumulative injuries in exercisers, but their functionality declines with increasing age . In addition, ageing and PA typically lead to different structural changes in muscle. In particular, strength training increases the size of muscle cells, while ageing is associated with smaller and weaker muscle cells and changes in the cell-type composition of tissue. However, some structural changes caused by ageing and PA are shared, both resulting in increasing infiltration of fat inside the muscles . However, in the case of PA, infiltrated fat is mainly lipid droplets that do not contain DNA. Muscle tissue samples, although carefully collected and with all visible fat and connective tissue removed, may still contain traces of these other tissues. Whether this explains our finding that both the Horvath’s and GrimAge clock estimates were significantly associated with the amount of body fat in muscle tissue remains unanswered.
In our study, slower epigenetic age acceleration in blood was weakly or non-significantly associated with higher PA, depending on the epigenetic age estimate and PA measurement method. Substantial associations were observed only with the GrimAge epigenetic clock and with self-reported PA data. The results are in line with those of previous studies, which have reported null findings [4, 29] using Horvath’s clock and weak associations using GrimAge [30,31,32]. The lack of consistent findings regarding the association between epigenetic ageing and PA in our study is likely explained by different features of epigenetic age algorithms and PA data collection methods. We have previously shown that LTPA is associated with slower ageing, while occupational PA is associated with accelerated ageing . This may explain why device-based measures of PA representing whole day PA did not show associations with epigenetic ageing, although they are typically considered more reliable than self-reports in estimating PA.
In line with previous studies [33, 34], associations between DNAmAge estimates and different physical function phenotypes were non-existent or minor in blood. Within the same individual, different tissues may yield different age estimates [13, 14, 35]. Some previous studies have suggested that tissue and cell type-specific analyses might be able to track age-related pathologies in target tissues. Muscle is the key tissue affecting age-related decline in physical functioning and PA, and thus we hypothesised that muscle tissue might better reflect biological ageing relevant to functional capacity. We found that the muscle-specific epigenetic clock MEAT is not associated with physical functioning, PA or body composition. MEAT was constructed only to predict chronological age without using any other biological phenotype data and may therefore be independent from molecular alterations caused by changes in the modifiable behavioural phenotypes, such as PA. Muscle-specific clocks need to be further developed taking into account biological and functional parameters, as has been done with respect to the other epigenetic clocks GrimAge and PhenoAge [9, 10], which, however, were built for blood tissue only. We also tested the two multi-tissue epigenetic clocks, Horvath’s and GrimAge with muscle tissue. In addition to the previously mentioned association with body fat, the only significant association was between Horvath’s age acceleration and cardiorespiratory fitness. However, this weak association attenuated after adjusting for sex and smoking, suggesting that it originates mostly from biological sex differences and the smoking habits of the participants.
Strengths and limitations
To the best of our knowledge, this study is the first to examine the associations between epigenetic clocks and PA and physical functioning in muscle tissue, utilising the novel muscle clock MEAT developed by Voisin et al. . Our methods for estimating PA and function covered a wide variety of phenotypes measured in a research laboratory by trained professionals as well as measurement of PA with reliable device-based methods under free-living conditions. By using two independent datasets, we were able to cover a large age span of individuals representing the same genetic ancestry. Since we only included individuals with both blood and muscle data available, we acknowledge that our sample size (in total N = 186) might not be sufficient to capture the small effect sizes that some phenotypes may have with respect to accelerated ageing. However, the observed trends in the associations between blood GrimAge age acceleration and self-reported PA corresponds to previously reported results in a larger sample . Our sample also included twins with discordant BMI, which may result in larger within-pair differences in methylation profiles. However, we conducted comprehensive sensitivity analyses, excluding either heavier or leaner twins from pairs discordant for BMI from the analysis, and the results were very similar (Additional file 3). Our results should be considered preliminary findings that need to be replicated with other, larger cohorts.
Both multi-tissue and muscle-specific novel epigenetic clocks seem to perform well in estimating chronological age in muscle tissue. Based on current knowledge and our results, however, these epigenetic clocks have rather low value in estimating muscle ageing with respect to physiological adaptations that typically occur due to ageing or PA. Practical issues related to the sampling of muscle tissue have limited the development of epigenetic research and scoring methods in this area. As age-related functional limitations are expanding rapidly due to increasing life expectancy throughout the world, there is an urgent need to gain more insight into muscle tissue-specific ageing and the underlying biological pathways. This type of knowledge could help the development of medical and lifestyle treatments for age-related functional limitations and muscle pathologies.
The data used in the current analysis originate from two Finnish population-based cohort studies: the FTC and the ERMA study.
In brief, FTC includes three large cohort studies: (1) study of twins born before 1958, an older twin cohort; (2) in 1975–1979, Finntwin16 (FT16); and (3) in 1983–1987, Finntwin12 (FT12) [36,37,38]. All three twin cohorts have been studied in an intensive longitudinal manner, and several subprojects have been conducted in laboratory conditions. All twins who had taken part in clinical in-person sub-studies with sampling for both whole blood and muscle tissue, subsequent DNA methylation analyses and phenotype collection were included (N = 139) [39, 40]. The dataset comprises two age groups: young adults from 23 to 44 years of age (N = 93) and older adults from 57 to 69 years of age (N = 46). The dataset includes 54 twin pairs, who were discordance for body mass index (BMI; difference more than 3 kg/m2) and have participated in the sub-study including laboratory measurements for body composition and physical performance . Zygosity of the twins was confirmed by genetic analysis of multiple polymorphic markers. The majority (45 pairs) were MZ, and only four same sex dizygotic twin pairs were present.
The recruitment process, study design and methodology of the ERMA study has been described in detail elsewhere . In brief, the original ERMA cohort includes 47- to 55-year-old women living in the city of Jyväskylä or neighbouring municipalities. The study participants were randomly selected from the national registers of the Digital and Population Data Services Agency, Helsinki (dvv.fi). The ERMA study was designed to reveal how hormonal differences during the menopausal stages affect the physiological and psychological functioning of middle-aged women. The participants for DNA extraction and subsequent methylation analysis (N = 47) were randomly selected amongst ERMA participants who had both whole blood and muscle tissue samples available and who did not have missing or low-quality accelerometer data in PA measures.
The FTC data collection and analysis were approved by the ethics committee of the Helsinki University Central Hospital (Dnro 249/E5/01, 270/13/03/01/2008, 154/13/03/00/2011), and the ERMA study was approved by the ethics committee of the Central Finland health care district in 2014 (K-SSHP Dnro 8U/2014). Written informed consent was provided by the participants before the sample collection and PA and functioning measurements.
Tissue sampling and DNA extraction
Whole blood samples were collected in fasting conditions from the antecubital vein. Muscle biopsy samples were obtained from the mid-part of the m. vastus lateralis, defined as the midpoint between the greater trochanter and the lateral joint line of the knee, by percutaneous needle biopsy under local anaesthesia. Following the removal of all visible fat residues, connective tissue and blood, the biopsy samples were frozen in liquid nitrogen and stored at -150 °C until use. Blood white cell and muscle tissue DNA were extracted using a DNeasy Blood & Tissue Kit (Qiagen, Cat No./ID: 69,504).
Genome-wide methylation analysis
High molecular weight DNA was bisulfite-converted using an EZ-96 DNA Methylation-Gold Kit (Zymo Research, Irvine, CA, USA). Genome-wide DNA methylation was measured using Illumina’s Infinium MethylationEPIC BeadChip (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. The Illumina BeadChip measures single-CpG resolution DNA methylation levels across the human genome.
Data were preprocessed using the R package minfi. Detection p-values comparing the total signal for each probe to the background signal level were calculated to evaluate the quality of the samples . Poor-quality samples (mean detection p > 0.01) were excluded from further analysis. Beta values representing CpG methylation levels were calculated as the ratio of methylated intensities (M) to the overall intensities (Beta value = M/(M + U + 100), where U is unmethylated probe intensity. Data were normalised using the single-sample Noob normalisation and Beta Mixture Quantile normalisation methods [43, 44].
The resulting methylation data were uploaded to an online epigenetic age calculator (https://dnamage.genetics.ucla.edu/new), selecting the options ‘Normalise data’ and ‘Advanced analysis’ to generate the various clock estimates . For blood, the four most commonly used epigenetic age estimates, Horvath’s and Hannum’s DNA methylation age (DNAmAge), PhenoAge and GrimAge, were computed. For muscle tissue, only Horvath’s DNAmAge and GrimAge were used, as Hannum’s clock and PhenoAge were clearly not suitable for muscle tissue according to current knowledge and our preliminary analyses, which showed very low DNAmAge estimates for muscle tissue samples (data not shown). In addition to these two epigenetic clocks, the novel MEAT was computed using the open access R-package MEAT available on Bioconductor . The data were screened for potential outliers e.g. due to methodological bias, or mistakes in sampling or data processing steps by comparing epigenetic age estimate of each sample to the mean epigenetic age of the study population. No outliers were detected in ERMA or FTC data. DNAmAge age acceleration, which describes the rate of individual ageing and is unrelated to chronological age, was calculated separately for FTC and ERMA participants [Age acceleration = residual (linear model (predicted age ~ chronological age)].
The Baecke questionnaire was used to assess PA in FTC . The questionnaire consists of 16 items, from which three indexes were calculated: the work index refers to PA at work, the sport index refers to sports participation during leisure time and the leisure-time index refers to PA during leisure time excluding sport activities. All responses were given on a five-point scale, with the exception of questions on main occupation and the types of two main sports. In the original publication, the test–retest reliability of the work index, sport index and leisure-time indices were 0.88, 0.81 and 0.74, respectively . The Baecke questionnaire has been validated for cardiorespiratory fitness among Finnish twins .
In the ERMA study, PA was monitored using hip-worn GT3X + or wGT3X + ActiGraph accelerometers at 60 Hz (Pensacola, FL) for seven consecutive days [47, 48]. The raw acceleration data were converted into 60-s epoch counts and divided into the mean daily duration spent at different PA intensities using tri-axial vector magnitude cut-points of 450 counts per minute (cpm), 2690 cpm and 6166 cpm to separate sedentary from light PA, light from moderate PA and moderate from vigorous PA, respectively; 25,000 cpm was used as an upper limit. Mean time spent at MVPA intensity level and daily steps during waking hours were calculated for each participant. Self-reported PA was assessed with a seven-scale questionnaire [49, 50]. Due to the low number of participants in some categories, answers were re-categorised as low, medium and high activity.
Anthropometrics and body composition were measured after overnight fasting. Body mass index was calculated as a function of measured body weight and height (kg/m2). Percentage body fat, total body fat mass and lean body mass were measured through dual-energy X-ray absorptiometry (DXA, Lunar Prodigy).
Physical function phenotypes included measures for both cardiorespiratory fitness and upper and lower extremity muscle strength and muscle power. Maximal oxygen consumption was measured only in FTC and other phenotypes and only in the ERMA cohort.
Maximal oxygen consumption, which describes the level of cardiorespiratory fitness, was measured with a graded maximal exercise test using an electrically braked bicycle ergometer, as described earlier . The initial workload (40 W for women and 50 W for men) was increased by 40 W (women) or 50 W (men) every 3 min until the maximal exercise capacity was reached (i.e. until volitional exhaustion following generally approved limits) . Oxygen uptake (V̇O2), carbon dioxide production (VCO2), ventilation, breathing frequency and other standardised respiratory parameters were measured continuously breath by breath with a Vmax spiroergometer (Sensormedics, Yorba Linda, CA). A strict calibration protocol was followed to ensure data quality .
The 6-min walking test, in which the participants attempted to walk the longest possible distance in the allotted time, was performed on a 20-m indoor track. Participants were instructed to complete as many laps as possible within 6 min, and the distance they walked was recorded . The 6-min walking test is a performance-based measure of functional exercise capacity, and it can be used as an indicator of exercise tolerance or cardiorespiratory fitness. A standardised protocol and safety guidelines issued by the American Thoracic Society (ATS Committee on Proficiency Standards for Clinical Pulmonary Function Laboratories 2002) were followed.
Hand grip strength (N) and knee extension torque (Nm) were measured using an isometric dynamometer (Good Strength, Metitur Oy, Jyväskylä, Finland) . To test hand grip strength, the participant’s dominant arm was fixed to the armrest of the dynamometer chair with the elbow flexed at a 90° angle, while knee extension strength was measured at a knee angle of 60° from full extension. During the measurements, participants were encouraged verbally to squeeze the handle or extend their knee to produce maximal contraction. Each participant was allowed to have three to five trials with 1 min rest between trials. The best attempt with the highest 2–3 s lasting peak force value was used in the analysis. Handgrip strength measurement can be used as an indicator of general muscle strength among older people . Maximal Nm was calculated from the best attempt measured in Newtons (N) multiplied by lower leg length (the distance between the lateral joint line of the knee and the midline of a cuff around the ankle) .
Vertical jumping height (cm) was calculated from flight time during a maximal effort countermovement jump and measured on a contact mat . Three maximal efforts were conducted, and the best result was used in the analysis. Vertical jumping height can be used as an estimate of lower body muscle power (i.e. the ability of the neuromuscular system to produce the greatest possible force as rapidly as possible). Lower extremity muscle power declines earlier and to a greater extent than muscle strength during ageing . It also correlates strongly with mobility among older people .
Sociodemographic variables and life habits were collected with validated questionnaires. These questionnaires have been published previously on the FTC webpage (https://wiki.helsinki.fi/display/twinstudy) and regarding the ERMA study as supplementary digital contents . Level of education was classified into primary, secondary and tertiary education groups following the levels of the Finnish education system. Use of alcohol was calculated as the number of alcoholic drinks (1 drink is 12 g ethanol) consumed per week, with alcohol consumption computed based on quantity and frequency measures of alcohol beverage use [59, 60]. Smoking was classified as never, former (abstinence for at least 6 months) and current (daily or occasional) smoker.
Statistical analyses were carried out with Stata 16 software (StataCorp, Inc. College Station, TX, USA) and R studio (3.6.0). Data are shown as means and standard deviations unless otherwise stated. The mean absolute difference between the predicted DNAmAge and chronological age was calculated to see whether the estimated DNAmAge was younger or older compared to the chronological age. Correlation coefficients were computed to estimate the level of similarity in DNAmAge estimates in different tissues. DNAmAge age acceleration was calculated for all participants as the residuals from a linear regression model of the DNAmAge estimate on chronological age. Intraclass correlation coefficients (ICC) were computed for the MZ twin pairs to estimate the level of within-pair similarity. Associations between different DNAmAge acceleration estimates, PA, body composition and PA-related phenotypes were analysed using linear regression analysis, which were further adjusted for sex and smoking. In the models, the within-pair dependency of twin individuals was taken into account using the cluster option in the analysis. The level of significance was set at p ≤ 0.05.
Availability of data and materials
All twin data used in the current study are located within the Biobank of the National Institute for Health and Welfare, Finland. All the biobanked data are publicly available for use by qualified researchers following a standardised application procedure.
Body mass index, kg/m2
A cytosine nucleotide followed by a guanine nucleotide
DNA methylation age
Estrogenic Regulation of Muscle Apoptosis study
Finnish twin cohort
Intraclass correlation coefficients
Leisure time physical activity
Muscle epigenetic age test
A moderate-to-vigorous physical activity
Rosenberg IH. Sarcopenia: origins and clinical relevance. J Nutr. 1997;127(5 Suppl):990S–991S.
Aasdahl L, Nilsen TIL, Meisingset I, Nordstoga AL, Evensen KAI, Paulsen J, et al. Genetic variants related to physical activity or sedentary behaviour: a systematic review. Int J Behav Nutr Phys Act. 2021;18(1):15–25.
Contrepois K, Wu S, Moneghetti KJ, Hornburg D, Ahadi S, Tsai MS, et al. Molecular choreography of acute exercise. Cell. 2020;181(5):1112–1130.e16.
Sillanpää E, Ollikainen M, Kaprio J, Wang X, Leskinen T, Kujala UM, et al. Leisure-time physical activity and DNA methylation age-a twin study. Clin Epigenetics. 2019;11(1):12–5.
Garatachea N, Lucia A. Genes, physical fitness and ageing. Ageing Res Rev. 2013;12(1):90–102.
Lopez-Otin C, Blasco MA, Partridge L, Serrano M, Kroemer G. The hallmarks of aging. Cell. 2013;153(6):1194–217.
Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14(10):R115.
Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49(2):359–67.
Levine ME, Lu AT, Quach A, Chen BH, Assimes TL, Bandinelli S, et al. An epigenetic biomarker of aging for lifespan and healthspan. Aging (Albany NY). 2018;10(4):573–91.
Lu AT, Quach A, Wilson JG, Reiner AP, Aviv A, Raj K, et al. DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging (Albany NY). 2019;11(2):303–27.
Chen BH, Marioni RE, Colicino E, Peters MJ, Ward-Caviness CK, Tsai PC, et al. DNA methylation-based measures of biological age: meta-analysis predicting time to death. Aging (Albany NY). 2016;8(9):1844–65.
Horvath S, Raj K. DNA methylation-based biomarkers and the epigenetic clock theory of ageing. Nat Rev Genet. 2018;19(6):371–84.
Horvath S, Mah V, Lu AT, Woo JS, Choi OW, Jasinska AJ, et al. The cerebellum ages slowly according to the epigenetic clock. Aging (Albany NY). 2015;7(5):294–306.
Horvath S, Erhart W, Brosch M, Ammerpohl O, von Schonfels W, Ahrens M, et al. Obesity accelerates epigenetic aging of human liver. Proc Natl Acad Sci USA. 2014;111(43):15538–43.
Park SW, Goodpaster BH, Strotmeyer ES, de Rekeneire N, Harris TB, Schwartz AV, et al. Decreased muscle strength and quality in older adults with type 2 diabetes: the health, aging, and body composition study. Diabetes. 2006;55(6):1813–8.
Volaklis KA, Halle M, Koenig W, Oberhoffer R, Grill E, Peters A, et al. Association between muscular strength and inflammatory markers among elderly persons with cardiac disease: results from the KORA-Age study. Clin Res Cardiol. 2015;104(11):982–9.
Volaklis KA, Halle M, Meisinger C. Muscular strength as a strong predictor of mortality: a narrative review. Eur J Intern Med. 2015;26(5):303–10.
Rantanen T, Guralnik JM, Foley D, Masaki K, Leveille S, Curb JD, et al. Midlife hand grip strength as a predictor of old age disability. JAMA. 1999;281(6):558–60.
Bonaldo P, Sandri M. Cellular and molecular mechanisms of muscle atrophy. Dis Model Mech. 2013;6(1):25–39.
Sillanpää E, Laaksonen DE, Hakkinen A, Karavirta L, Jensen B, Kraemer WJ, et al. Body composition, fitness, and metabolic health during strength and endurance training and their combination in middle-aged and older women. Eur J Appl Physiol. 2009;106(2):285–96.
Sillanpää E, Hakkinen A, Nyman K, Mattila M, Cheng S, Karavirta L, et al. Body composition and fitness during strength and/or endurance training in older men. Med Sci Sports Exerc. 2008;40(5):950–8.
Voisin S, Harvey NR, Haupt LM, Griffiths LR, Ashton KJ, Coffey VG, et al. An epigenetic clock for human skeletal muscle. J Cachexia Sarcopenia Muscle. 2020;11(4):887–98.
Sillanpää E, Cheng S, Hakkinen K, Finni T, Walker S, Pesola A, et al. Body composition in 18- to 88-year-old adults–comparison of multifrequency bioimpedance and dual-energy X-ray absorptiometry. Obesity (Silver Spring). 2014;22(1):101–9.
Christiansen L, Lenart A, Tan Q, Vaupel JW, Aviv A, McGue M, et al. DNA methylation age is associated with mortality in a longitudinal Danish twin study. Aging Cell. 2016;15(1):149–54.
Alway SE, Myers MJ, Mohamed JS. Regulation of satellite cell function in sarcopenia. Front Aging Neurosci. 2014;6:246.
Turner DC, Gorski PP, Maasar MF, Seaborne RA, Baumert P, Brown AD, et al. DNA methylation across the genome in aged human skeletal muscle tissue and muscle-derived cells: the role of HOX genes and physical activity. Sci Rep. 2020:10(1):15360.
Henze H, Jung MJ, Ahrens HE, Steiner S, von Maltzahn J. Skeletal muscle aging—stem cells in the spotlight. Mech Ageing Dev. 2020;189:111283.
Li X, Li Z, Zhao M, Nie Y, Liu P, Zhu Y, et al. Skeletal muscle lipid droplets and the athlete’s paradox. Cells. 2019. https://doi.org/10.3390/cells8030249.
Quach A, Levine ME, Tanaka T, Lu AT, Chen BH, Ferrucci L, et al. Epigenetic clock analysis of diet, exercise, education, and lifestyle factors. Aging (Albany NY). 2017;9(2):419–46.
Kankaanpää A, Tolvanen A, Bollepalli S, Leskinen T, Kujala UM, Kaprio J, et al. Leisure-time and occupational physical activity associates differently with epigenetic aging. Med Sci Sports Exerc. 2021;53(3):487–95.
McCrory C, Fiorito G, Hernandez B, Polidoro S, O'Halloran AM, Hever A, et al. GrimAge outperforms other epigenetic clocks in the prediction of age-related clinical phenotypes and all-cause mortality. J Gerontol A Biol Sci Med Sci 2021;76(5):741–749.
Kresovich Jacob K, Taylor JA, Sandler DP, White AJ, Xu Z, Niehoff NM, et al. Body composition and physical activity associations with multiple measures of epigenetic age acceleration. Am J Epidemiol. 2020. https://doi.org/10.1093/aje/kwaa251.
Marioni RE, Shah S, McRae AF, Ritchie SJ, Muniz-Terrera G, Harris SE, et al. The epigenetic clock is correlated with physical and cognitive fitness in the Lothian Birth Cohort 1936. Int J Epidemiol. 2015;44(4):1388–96.
Sillanpää E, Laakkonen EK, Vaara E, Rantanen T, Kovanen V, Sipilä S, et al. Biological clocks and physical functioning in monozygotic female twins. BMC Geriatr. 2018;18(1):83–6.
Levine AJ, Quach A, Moore DJ, Achim CL, Soontornniyomkij V, Masliah E, et al. Accelerated epigenetic aging in brain is associated with pre-mortem HIV-associated neurocognitive disorders. J Neurovirol. 2016;22(3):366–75.
Kaidesoja M, Aaltonen S, Bogl LH, Heikkila K, Kaartinen S, Kujala UM, et al. FinnTwin16: a longitudinal study from age 16 of a population-based Finnish Twin Cohort. Twin Res Hum Genet. 2019;22(6):530–9.
Rose RJ, Salvatore JE, Aaltonen S, Barr PB, Bogl LH, Byers HA, et al. FinnTwin12 Cohort: an updated review. Twin Res Hum Genet. 2019;22(5):302–11.
Kaprio J, Bollepalli S, Buchwald J, Iso-Markku P, Korhonen T, Kovanen V, et al. The older Finnish Twin Cohort—45 years of follow-up. Twin Res Hum Genet. 2019;22(4):240–54.
Huang Y, Ollikainen M, Muniandy M, Zhang T, van Dongen J, Hao G, et al. Identification, heritability, and relation with gene expression of novel DNA methylation loci for blood pressure. Hypertension. 2020;76(1):195–205.
Naukkarinen J, Rissanen A, Kaprio J, Pietilainen KH. Causes and consequences of obesity: the contribution of recent twin studies. Int J Obes (Lond). 2012;36(8):1017–24.
Kovanen V, Aukee P, Kokko K, Finni T, Tarkka IM, Tammelin T, et al. Design and protocol of Estrogenic Regulation of Muscle Apoptosis (ERMA) study with 47 to 55-year-old women’s cohort: novel results show menopause-related differences in blood count. Menopause. 2018;25(9):1020–32.
Maksimovic J, Phipson B, Oshlack A. A cross-package bioconductor workflow for analysing methylation array data. F1000Res. 2016;5:1281.
Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, et al. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29(2):189–96.
Fortin JP, Triche TJ, Hansen KD. Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics. 2017;33(4):558–60.
Baecke JA, Burema J, Frijters JE. A short questionnaire for the measurement of habitual physical activity in epidemiological studies. Am J Clin Nutr. 1982;36(5):936–42.
Mustelin L, Latvala A, Pietilainen KH, Piirila P, Sovijarvi AR, Kujala UM, et al. Associations between sports participation, cardiorespiratory fitness, and adiposity in young adult twins. J Appl Physiol. 2011;110(3):681–6.
Laakkonen EK, Kulmala J, Aukee P, Hakonen H, Kujala UM, Lowe DA, et al. Female reproductive factors are associated with objectively measured physical activity in middle-aged women. PLoS ONE. 2017;12(2):e0172054.
Sasaki JE, John D, Freedson PS. Validation and comparison of ActiGraph activity monitors. J Sci Med Sport. 2011;14(5):411–6.
Hyvarinen M, Sipila S, Kulmala J, Hakonen H, Tammelin TH, Kujala UM, et al. Validity and reliability of a single question for leisure-time physical activity assessment in middle-aged women. J Aging Phys Act. 2019;28(2):231–41.
Hirvensalo M, Lampinen P, Rantanen T. Physical exercise in old age: an eight-year follow-up study on involvement, motives, and obstacles among persons age 65–84. J Aging Phys Act. 1998;6:157–68.
Mustelin L, Pietilainen KH, Rissanen A, Sovijarvi AR, Piirila P, Naukkarinen J, et al. Acquired obesity and poor physical fitness impair expression of genes of mitochondrial oxidative phosphorylation in monozygotic twins discordant for obesity. Am J Physiol Endocrinol Metab. 2008;295(1):148.
Enright PL. The six-minute walk test. Respir Care. 2003;48(8):783–5.
Ronkainen PH, Kovanen V, Alen M, Pollanen E, Palonen EM, Ankarberg-Lindgren C, et al. Postmenopausal hormone replacement therapy modifies skeletal muscle composition and function: a study with monozygotic twin pairs. J Appl Physiol. 2009;107(1):25–33.
Rantanen T, Era P, Heikkinen E. Maximal isometric strength and mobility among 75-year-old men and women. Age Ageing. 1994;23(2):132–7.
Bondarev D, Finni T, Kokko K, Kujala UM, Aukee P, Kovanen V, et al. Physical performance during the menopausal transition and the role of physical activity. J Gerontol A Biol Sci Med Sci. 2020;. https://doi.org/10.1093/gerona/glaa292.
Bosco C, Luhtanen P, Komi PV. A simple method for measurement of mechanical power in jumping. Eur J Appl Physiol Occup Physiol. 1983;50(2):273–82.
Bean JF, Leveille SG, Kiely DK, Bandinelli S, Guralnik JM, Ferrucci L. A comparison of leg power and leg strength within the InCHIANTI study: which influences mobility more? J Gerontol A Biol Sci Med Sci. 2003;58(8):728–33.
Runge M, Rittweger J, Russo CR, Schiessl H, Felsenberg D. Is muscle power output a key factor in the age-related decline in physical performance? A comparison of muscle cross section, chair-rising test and jumping power. Clin Physiol Funct Imaging. 2004;24(6):335–40.
Sipilä P, Rose RJ, Kaprio J. Drinking and mortality: long-term follow-up of drinking-discordant twin pairs. Addiction. 2016;111(2):245–54.
Stephenson M, Bollepalli S, Cazaly E, Salvatore JE, Barr P, Rose RJ, et al. Associations of alcohol consumption with epigenome-wide DNA methylation and epigenetic age acceleration: individual-level and co-twin comparison analyses. Alcohol Clin Exp Res. 2020;45(2):318–28.
ES was supported by the Academy of Finland (Grant No. 260001), the Juho Vainio foundation and the Yrjö Jahnsson foundation. VK was supported by the Academy of Finland (275323). KHP was funded by the Academy of Finland (Grant Nos. 314383 and 266286), the Academy of Finland Center of Excellence in Research on Mitochondria, Metabolism and Disease (FinMIT; Grant No. 272376), the Finnish Medical Foundation, the Gyllenberg Foundation, the Novo Nordisk Foundation (Grant Nos. NNF17OC0027232 and NNF10OC1013354), the Finnish Diabetes Research Foundation, the Finnish Foundation for Cardiovascular Research, Government Research Funds, the University of Helsinki and Helsinki University Hospital. JK was supported by the Academy of Finland (213506, 308248, 312073, 336823), EC FP5 GenomEUtwin, NIH NIH/NHLBI (grant HL104125), EC MC ITN Project EPITRAIN and Sigrid Juselius Foundation. MO was supported by the Academy of Finland (297908 and 251316), EC MC ITN Project EPITRAIN, University of Helsinki Research Funds and Sigrid Juselius Foundation. EKL was supported by the Academy of Finland (309504 and 314181) and the Juho Vainio Foundation. The Gerontology Research Center is a joint effort between the University of Jyväskylä and the University of Tampere.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: Within-pair correlations in age acceleration in blood and in muscle. Additional file 2: Associations between DNAmAge age acceleration estimates and body composition and physical activity in blood. Additional file 3: Sensitivity analyses related to twin pair discordance in body mass index.
About this article
Cite this article
Sillanpää, E., Heikkinen, A., Kankaanpää, A. et al. Blood and skeletal muscle ageing determined by epigenetic clocks and their associations with physical activity and functioning. Clin Epigenet 13, 110 (2021). https://doi.org/10.1186/s13148-021-01094-6