- Research
- Open access
- Published:
Association between modifiable lifestyle factors and telomere length: a univariable and multivariable Mendelian randomization study
Journal of Translational Medicine volume 22, Article number: 160 (2024)
Abstract
Background
Telomere length has long been recognized as a valuable biomarker of aging and is inversely correlated with chronological age. Various lifestyle factors have been implicated in telomere shortening or preservation; however, the association between lifestyle factors and telomere length remains controversial. To address this issue, we conducted a Mendelian randomization (MR) analysis to investigate the potential causal associations between multiple lifestyle factors and telomere length.
Methods
Independent genetic variants strongly associated with lifestyle factors (tobacco smoking, sleep duration, insomnia, and physical activity) were selected as instrumental variables from corresponding genome-wide association studies (GWASs). Summary-level data for telomere length was obtained from a GWAS comprising 472,174 European ancestries. Univariable and multivariable MR analyses were performed to assess the relationships.
Results
The genetic liability to lifetime smoking was robustly associated with shorter telomere length (odd ratio [OR]: 0.882; 95% confidence interval [CI]: 0.847–0.918). Genetically predicted insomnia was also linked to shorter telomere length (OR: 0.972; 95% CI: 0.959–0.985), while no significant association was observed between sleep duration and telomere length. Furthermore, a suggestive association was found between moderate-to-vigorous physical activity and longer telomere length (OR: 1.680; 95% CI: 1.115–2.531). In multivariable MR analyses, adjusting for potential mediators such as body mass index, type 2 diabetes, alcohol consumption, and alcohol use disorder, the associations of lifetime smoking and insomnia with telomere length remained robust.
Conclusion
Our findings suggest that smoking and insomnia may contribute to telomere shortening, while physical activity may play a role in telomere length maintenance. These findings underscore the importance of managing positive risk factors and adopting a healthy lifestyle to promote telomere health.
Introduction
With medical and technological advances in recent decades, human life expectancy has increased worldwide and is expected to increase further in the coming years [1]. Population aging is accompanied by an increase in the incidence of age-related diseases, such as cardiovascular diseases [2], diabetes [3], and cancer [4]. For this reason, aging has become a critical public health issue in many countries. Currently, telomere length has been developed as a surrogate measure of biological age [5]. Telomeres are conserved repetitive DNA sequences, together with associated protective protein, located at the ends of cap-shaped linear chromosomes. With each cell division, telomeres naturally shorten because the DNA polymerases cannot wholly replicate the end of the lagging DNA strand. Thus, telomere shortening occurs with increasing cell age. Telomere length is partially determined by genetic factors and correlates with sex, race, and paternal age [6, 7]. It also correlates with environmental and lifestyle factors, including sleep, exercise, tobacco use and alcohol consumption [8, 9].
Lifestyle factors and telomere length have been investigated in observational studies with conflicting results. The largest cross-sectional study, conducted on 7,813 women from the Nurses’ Health Study, found that even moderate physical activity was associated with longer telomere length compared to less active individuals [10]. Similarly, a study involving 2,312 American Indian participants reported a positive association between higher levels of physical activity and longer telomeres [11]. However, a few studies with smaller sample sizes have shown no significant association between physical activity and telomere length [12,13,14]. The relationship between tobacco smoking and telomere length has also been investigated in observational studies with mixed findings. A dose-dependent adverse association was observed between the number of cigarettes smoked and telomere length in a study involving 8,074 participants [15]. Similarly, a longitudinal study of 5,624 participants from the Health and Retirement Study found an inverse association between the number of cigarettes smoked and telomere length, particularly in women [16]. However, a case–control study of 29 healthy smokers and 29 non-smokers reported no impact of tobacco smoking on telomere length [17]. In comparison to physical activity and smoking, few studies have specifically focused on sleep duration. Several studies have shown a significant association between shorter sleep duration and shorter telomere length [18,19,20]. Additionally, conflicting results have been reported regarding the effect of long sleep duration on telomere length, with some studies suggesting a U-shaped relationship or a linear effect [18, 20].
Conducting randomized controlled trials to examine the impact of lifestyle factors on telomere length is impractical due to ethical considerations. Moreover, the causal nature of the associations between the aforementioned lifestyle factors and telomere length remains uncertain in observational studies due to potential residual confounding and reverse causality issues. Consequently, alternative methods that facilitate causal inference can provide valuable insights into whether these factors represent potentially modifiable risk factors. Mendelian randomization (MR) is an analytical approach that utilizes single-nucleotide polymorphisms (SNPs) as reliable instrumental variables to assess the potential causal relationship between an exposure and an outcome. [21]. Conceptually, MR shares similarities with randomized controlled trials as genetic variants are randomly assigned during meiosis, thereby reducing concerns related to confounding and reverse causality.
In MR analyses, instrumental variables can be derived from summary statistics of extensive genome-wide association studies (GWASs) that have been conducted on lifestyle factors and telomere length. As of now, the causal relationship between the mentioned lifestyle factors and telomere length has yet to be established through MR analysis. Hence, in this study, we employed MR analysis to investigate the potential causal association between lifestyle factors (tobacco smoking, sleep duration, insomnia, and physical activity) and the telomere length.
Method
Study design
We conducted a MR analysis utilizing SNPs as instrumental variables for the exposures under investigation. The study design is illustrated in Fig. 1. To ensure the validity of causal estimates in MR studies [21], three key assumptions need to be satisfied: First, the selected genetic variants serving as instrumental variables should be strongly associated with exposures. Second, genetic variants should not be associated with any confounders that could potentially influence the relationship between the exposures and the outcome. Third, the genetic variants should not independently affect the outcome apart from their effects on the exposures. Genetic instrument selection based on publicly available GWAS for physical activity [22], insomnia [23], sleep duration [24], and lifetime smoking [25]. Summary data for telomere length were derived from the largest-to-date GWAS [26]. Additional file 1: Table S1 provides detailed information on the data sources used in this study. The GWAS data used in this study have been made publicly available, and all ethical approval was obtained in all original studies.
Data sources and instruments
Lifetime smoking
The measures of smoking used in our study were the lifetime smoking index. The SNPs associated with lifetime smoking index at the genome-wide significance level (p-value < 5 × 10–8) were extracted from a GWAS including 462,690 individuals from the UK Biobank. [25]. The lifetime smoking index aimed to capture various aspects of smoking, including smoking heaviness, duration, and initiation. It considered factors such as smoking status, age at initiation and cessation, and the number of cigarettes smoked per day. In the GWAS analysis, linear mixed models adjusted for sex and genotyping chip were conducted using BOLT-LMM. To ensure independence among the SNPs, clumping was performed with a linkage disequilibrium (LD) threshold of r2 = 0.01 and a distance of 10,000 kb using the TwoSampleMR package [27], utilizing the LD reference panel from the 1000 Genomes Project specifically for individuals of European descent [28]. Ultimately, 126 SNPs were identified as genetic instruments for lifetime smoking index (Additional file 1: Table S2). It is noteworthy that an increase of one standard deviation in the lifetime smoking score corresponds to a smoking history equivalent to either 20 cigarettes per day for 15 years with cessation 17 years ago or 60 cigarettes per day for 13 years with cessation 22 years ago [25].
Sleep duration and insomnia
The SNPs associated with sleep duration were obtained from a recent GWAS of individuals with European ancestry in the UK Biobank, which included a total of 446,118 participants [24]. Phenotypic data of sleep duration were self-reported by participants based on responses to questionnaires or interviews. Participants with sleep duration less than 3 h or more than 18 h, as well as those using sleep medication, were excluded. Sleep duration was considered both as a continuous variable and as categorical, with categories defined as short (6 h or less), normal (7 or 8 h), or long (9 h or more) sleep duration [24]. Genetic association analysis was performed using BOLT-LMM linear mixed models and an additive genetic model adjusted for age, sex, ten principal components of ancestry, and genetic correlation matrix. The GWAS identified 78 SNPs that reached genome-wide significance (p-value < 5 × 10–8), collectively accounting for 0.69% of the variance in sleep duration. Following clumping based on a LD threshold of r2 = 0.01 and a distance of 10 000 kb, a total of 74 independent SNPs were identified as instrumental variables for sleep duration (Additional file 1: Table S3). Additionally, separate GWAS for short (< 7 h) and long sleep duration (≥ 9 h) revealed 26 and 7 significant SNPs, respectively (Additional file 1: Tables S4 and S5).
We utilized SNPs associated with insomnia from a GWAS that involved 1,331,010 individuals of European ancestry (397,959 cases and 933,057 controls) [23]. Insomnia was defined by self-reported insomnia symptoms based on responses to questionnaires in a touchscreen device or an online survey. Genetic association analysis was conducted using PLINK logistic regression and an additive genetic model adjusted for age, sex, ten principal components, and genotyping array. The GWAS identified 248 significant SNPs, which collectively accounted for approximately 2.6% of the variance in insomnia. Among these, 207 independent SNPs were selected as instrumental variables following clumping based on the aforementioned strategy (Additional file 1: Tables S6). Sensitivity analysis excluding shift workers or individuals who reported using psychiatric medications or have psychiatric illnesses revealed a strong association with the primary GWAS [23].
Physical activity
We obtained SNPs associated with physical activity from a recent GWAS [22]. Moderate physical activity (MPA) and vigorous physical activity (VPA) were defined based on self-reported information collected through a touchscreen questionnaire. Moderate-to-vigorous physical activity (MVPA) was calculated as the combination of MPA and VPA, following the previously described methodology [29]. Participants with missing physical activity data or reporting more than 16Â h of physical activity were excluded from the analysis. To adjust the phenotype, the GWAS employed regression analysis, where the dependent physical activity variables were regressed on several independent variables, including age, sex, genotyping chip, the first ten genomic principal components, center, and season (month) at the center visit or during accelerometer wear. Within the MVPA category, the GWAS identified nine significantly associated SNPs using data from 377,234 individuals (Additional file 1: Table S7), explaining approximately 0.09% of the variance in exposure. Moreover, five loci showed significant associations with vigorous physical activity (VPA), and six loci were significantly associated with strenuous sports or other exercises (SSOE) (Additional file 1: Tables S8 and S9).
Telomere length
Summary-level data for the exposure-associated SNPs with telomere length came from the largest GWAS in 472,174 UK Biobank participants [26]. Participants between the ages of 45 and 69 were recruited and extracted DNA from peripheral blood leukocytes using an automated process or manual method [30]. Leukocyte telomere length is measured as the ratio of telomere repeats copy number relative to that of a single copy gene using the qPCR method. Samples were re-run later to assess the measurements’ stability and repeatability. GWAS was performed using BOTL-LMM adjusting for age, sex, genotype array, and the first ten principal components.
Potential confounders
The effects of genetic variants on body mass index were obtained from publicly available summary statistics [31]. The GWAS meta-analysis included 339,224 individuals from 125 studies. Body mass index was measured or self-reported weight in kg per height in meters squared. The summary statistics of alcohol consumption were derived from the largest GWAS, including 1.2 million individuals of European ancestry [32]. Alcohol consumption was measured as weekly drinks (n = 941,280 individuals). We obtained the genetic variants for alcohol use disorder from the IEU Open GWAS Project (https://gwas.mrcieu.ac.uk), which included 12,204 cases and 206,588 controls.
Testing instrument strength and statistical power
F-statistics were used to evaluate the instrument strength by using the formula: F = R2 × (N − 1)/(1 − R2) [33], where N indicates the GWAS sample size. R2 is calculated using the formula: R2 = 2 × EAF × (1 − EAF) × beta2, where EAF is the effect allele frequency, and beta is the estimated effect on exposure [34]. The risk of weak instrument bias in MR analysis is relatively low when the F-statistics ≥ 10 [35]. We used an online tool (https://shiny.cnsgenomics.com/mRnd/) to estimate the power of our study based on the GWAS sample size and the variance explained by genetic instruments for the exposure [36]. The post hoc power calculations for our main IVW analyses were presented in Additional file 1: Table S10.
Statistics analyses
For primary MR analysis, we used inverse variance-weighted (IVW) based on a multiplicative random-effects model as the primary method. The Wald ratio estimate was calculated by dividing the effect size (β per additional risk allele) for the SNP-telomere length association by the corresponding coefficient for the SNP-exposure association [37]. Subsequently, we combined the estimates for each SNP to derive a single causal estimate for each exposure.
To assess the robustness of our findings, we conducted various sensitivity analyses to evaluate the impact of potential pleiotropic effects on our causal estimates. In addition to the IVW method, we employed the weighted median method. The weighted median estimate of a causal effect remains consistent as long as at least 50% of the genetic variants are valid instrumental variables [38]. Cochran’s Q and I2 statistics were used to evaluate the heterogeneity produced by different instrumental variables in the analysis [39]. A p-value of Cochran’s Q < 0.05 indicated the presence of heterogeneity, and the I2 statistic reflected the level of heterogeneity (I2 < 25%: small heterogeneity; I2 = 25–75%: moderate heterogeneity; I2 > 75%: large heterogeneity). When heterogeneity was detected, a random-effects IVW analysis was more appropriate. Additionally, we performed the MR-Egger intercept test to detect the potential presence of horizontal pleiotropy, with a deviation from zero (p-value < 0.05) indicating evidence of directional pleiotropic bias [40]. In the presence of horizontal pleiotropy, the slope coefficient from the MR-Egger regression provided a reliable estimate of the causal effect. MR pleiotropy residual sum and outlier (MR-PRESSO) was used to identify and correct horizontal pleiotropic outliers [41]. If any outlying SNPs were identified, we reassessed the effect estimates by excluding them from the exposure instruments.
Given the possibility of confounding factors influencing the relationship between these exposures and telomere length, we performed multivariable MR analyses to assess the direct effect of exposures on telomere length whilst accounting for potential confounding effects of body mass index, type 2 diabetes, alcohol consumption, and alcohol use disorder. Additionally, considering the positive association between tobacco smoking and insomnia [23], we also conducted multivariable MR analyses to adjust for the lifetime smoking index and insomnia in relation to each other.
In our complementary analysis, we employed statistical methods similar to those described above to investigate the causal relationships between short sleep duration, long sleep duration, VPA, and SSOE with telomere length.
Statistical significance was defined as a p-value < 0.05. To account for multiple exposures, we adjusted the significance threshold using Bonferroni correction, resulting in a p-value of < 0.0125 (= 0.05/4 exposures) considered statistically significant. All analyses were conducted using R version 3.6.1 with the TwoSampleMR (version 0.5.5), MendelianRandomization (version 0.4.3), and MRPRESSO (version 1.0) packages.
Results
Lifetime smoking index
The genetic liability to lifetime smoking was strongly associated with shorter telomere length (OR: 0.882; 95% CI: 0.847–0.918; p = 9.22 × 10–10) (Fig. 2). Consistent findings were obtained from the weight median method, as shown in Additional file 1: Table S11. Evaluation of heterogeneity using Cochran’s Q and I2 statistics revealed moderate variability across SNPs (Additional file 1: Table S12). Furthermore, the intercept of MR-Egger intercept indicated no evidence of directional pleiotropy (0.001; p = 0.532) (Additional file 1: Table S12). Outlier detection using MR-PRESSO identified three SNPs, but upon removal of these variants, the causal estimate remained unchanged (Additional file 1: Table S13). In the multivariable analysis adjusting for body mass index, type 2 diabetes, alcohol consumption, alcohol use disorder, the associations between genetic susceptibility to lifetime smoking and telomere length did not diminish and remained statistically significant (Fig. 3). Considering the observed genetic association between smoking and insomnia, we additionally accounted for insomnia in the multivariable MR analysis. Notably, the inclusion of insomnia as a covariate did not alter the association between lifetime smoking and telomere length (Fig. 3).
Insomnia
Genetic predisposition to insomnia was significantly associated with shorter telomere length, as evident in the random-effects IVW analysis (OR: 0.972; 95% CI: 0.959–0.985; p = 1.76 × 10–5) and weight median analysis (OR: 0.980; 95% CI: 0.968–0.991; p = 5.34 × 10–4) (Fig. 2 and Additional file 1: Table S11). Cochran’s Q and I2 statistics revealed moderate heterogeneity across SNPs (Additional file 1: Table S12). The intercept of the MR-Egger was close to zero (-0.001; p = 0.626), suggesting the absence of directional pleiotropy (Additional file 1: Table S12). Through MR-PRESSO, nine SNPs were identified as potential outliers, yet exclusion of these variants did not alter the results (Additional file 1: Table S13). Moreover, even after adjusting for potential confounders such as body mass index, type 2 diabetes, alcohol consumption, and alcohol use disorder, the association between genetic liability to insomnia and telomere length remained statistically significant (Fig. 4). However, upon further adjustment for lifetime smoking, the effect sizes of the association slightly attenuated and did not reach statistical significance after Bonferroni correction (OR: 0.986; 95% CI: 0.975–0.998; p = 0.020) (Fig. 4).
Sleep duration
There was no significant association found between genetically predicted sleep duration and telomere length (OR: 1.059; 95% CI: 0.984–1.140; p = 0.123) (Fig. 2). The weighted median method yielded consistent estimates with the primary analysis (Additional file 1: Table S11). Notably, there was substantial heterogeneity across SNPs for sleep duration, as indicated by significant Cochran’s Q and I2 statistics (Additional file 1: Table S12). Moreover, no evidence of directional pleiotropy was observed (MR-Egger intercept: 0.001; p = 0.605). The estimate remained unchanged even after excluding five outlier SNPs identified by the MR-PRESSO method (Additional file 1: Table S13).
In the complementary analyses, we employed similar statistical methods to investigate the causal effects of short and long sleep duration on telomere length. The findings suggested a potential association between short sleep duration and telomere length (OR: 0.959; 95% CI: 0.923–0.996; p = 0.032). However, genetically predicted long sleep duration showed no association with telomere length (OR: 1.004; 95% CI: 0.959–1.051; p = 0.865).
Physical activity
The results revealed a suggested association between genetically predicted MVPA and longer telomere length (OR: 1.680; 95% CI: 1.115–2.531; p = 0.013) (Fig. 2). However, the estimate obtained from the weighted median method showed inconsistency with the primary analysis (Additional file 1: Table S11). The presence of high heterogeneity among these SNPs was assessed using Cochran’s Q and overall I2 statistics (Additional file 1: Table S12). The MR-Egger intercept test detected a possible directional pleiotropy (-0.033; p = 0.007). After accounting for pleiotropy in the MR-Egger analysis, the genetic liability to MVPA exhibited a positive association with telomere length (OR 11.934; 95% CI, 4.167- 34.180; p = 0.002) (Additional file 1: Table S11). By applying the MR-PRESSO method, six SNPs were identified as outliers, and even after excluding these SNPs, a suggested association between genetically predicted MVPA and longer telomere length persisted (OR: 2.482; 95% CI: 1.171–5.263; p = 0.018) (Additional file 1: Table S13).
Complementary analyses were conducted to investigate the association between VPA and SSOE with telomere length. In contrast to MVPA, no causal associations were observed between genetically predicted VPA (OR: 0.988; 95% CI: 0.945–1.034; p = 0.606), SSOE (OR: 1.019; 95% CI: 0.963–1.080; p = 0.511), and telomere length.
Discussion
Principal findings
In this MR study, we utilized large-scale human genetics data to investigate the potential causal relationships between physical activity, lifetime smoking, sleep duration, insomnia, and telomere length. By leveraging data from multiple large-scale GWASs, we were able to incorporate a larger number of cases involving lifestyle factors and telomere length compared to previous observational studies. Notably, we observed robust associations between lifetime smoking, insomnia, and shorter telomere length. While suggestive associations were found between MVPA and longer telomere length, no significant association was observed between sleep duration and telomere length. Importantly, the associations of lifetime smoking and insomnia with telomere length remained significant even after adjusting for confounding factors such as alcohol consumption, alcohol use disorder, body mass index, and type 2 diabetes.
Previous researches
The association between lifetime smoking and shorter telomere length was strong and likely to be causal, as indicated by consistent estimates across the employed MR methods. Furthermore, even after adjusting for potential confounding factors, the association between smoking and telomere length remained significant, providing further support for the plausibility of the hypothesis. Notably, findings from prospective observational studies yielded similar results. For instance, a comprehensive longitudinal study involving 5,624 participants over a 16-year period revealed a consistent link between smoking and reduced telomere length, albeit with variations by gender [16]. Another study [42] reported a negative association between smoking and telomere length in older adults, indicating a potential difference of 73 base pairs between smoking and non-smoking individuals. Additionally, a large longitudinal study encompassing 1,356 individuals aged 30–70 suggested that smoking status was associated with both telomere length attrition and baseline telomere length, further supporting the adverse impact of smoking on telomeres [43].
Deriving a causal relationship between sleep duration and telomere length is not so straightforward. In this MR study, no compelling evidence was found for a direct causal effect of continuous sleep duration on telomere length. However, when utilizing SNPs associated with short and long sleep durations to estimate causal effects, the results suggested a link between genetically predicted short sleep duration (< 7 h) and shorter telomere length, while no significant association was observed for long sleep duration (≥ 9 h). Observational studies investigating the association between sleep duration and telomere length have yielded conflicting findings. Zhao et al. [44] conducted a study that revealed a positive correlation between telomere length and daily sleeping time. In contrast, another study found an association between long sleep duration and shorter telomere length, even after adjusting for factors such as sex, age, and body mass index. However, two cross-sectional studies indicated that individuals with shorter sleep duration exhibited decreased telomere length [45, 46]. Insomnia, primarily a mental disorder characterized by poor sleep quality and shortened sleep duration, was also explored in our MR analyses. The results further supported a connection between genetically predicted insomnia and shorter telomere length. In line with our findings, an observational study involving 925 individuals reported a relationship between insomnia disorder (OR = 2.654, 95% CI = 1.025–6.873) and shortened telomere length [20].
Drawing definitive conclusions about the causal role of MVPA in telomere length is currently not possible. However, our MR analyses revealed a suggestive positive association between genetic predisposition to MVPA and telomere length. The relatively wide 95% confidence interval for the estimate may be attributed to the small variance explained (0.09%) by MVPA. Furthermore, our additional analyses did not indicate any associations between VPA, SSOE, and telomere length. Several large cross-sectional studies have reported a positive association between physical activity and telomere length [47,48,49]. Additionally, two longitudinal studies in elderly subjects demonstrated a positive relationship between physical activity and telomere length [50, 51]. However, the largest longitudinal study, involving 4,576 participants from the Copenhagen City Heart Study, found no association between changes in physical activity and telomere length over a ten-year follow-up period [52]. Results from intervention studies have also been mixed. Most trials with small sample sizes reported that participants in active physical activity interventions showed increased telomere length compared to control groups [53,54,55]. In contrast, a trial involving 439 women comparing telomere length after a 12-month physical activity intervention did not find any significant differences [56]. The optimal type and amount of physical activity required to maintain telomere length remain a topic of debate, indicating the need for further research.
Potential mechanisms
The precise pathophysiological mechanisms underlying the association between lifestyle factors and telomere length remain incompletely understood. However, oxidative stress and inflammation have been identified as potential explanations for this relationship [57, 58]. Both tobacco smoking and insomnia have been shown to increase the levels of inflammatory cytokines and pro-oxidants, thereby disrupting the balance between inflammatory and oxidative stress pathways [59, 60]. Conversely, physical activity has been found to reduce inflammation and oxidative stress by decreasing the production of reactive oxygen species, tumor necrosis factor α, C-reactive protein, and interleukin-6 [61, 62]. Moreover, exercise has been associated with increased plasma levels of irisin and insulin-like growth factor 1, which are directly linked to telomere length [63,64,65].
Strengths and limitations
This study possessed several strengths. The primary strength lies in the MR design, which mitigates confounding biases often present in observational studies. By utilizing genetic associations obtained from available GWAS as instrumental variables for different lifestyle factors, we were able to address potential horizontal pleiotropy. Additionally, we conducted numerous sensitivity analyses, including multivariate MR, to assess the robustness of our findings. Another strength is the inclusion of participants of European ancestry exclusively, thereby minimizing the influence of population stratification bias.
Nevertheless, several limitations were observed in this study. Firstly, we cannot entirely dismiss the possibility that the genetic variants associated with lifestyle factors also influence telomere length through alternative causal pathways. However, our robust associations persisted even after adjusting for major genetically related risk factors such as alcohol use and body mass index in multivariate MR analyses. Another limitation is our inability to explore potential nonlinear associations between physical activity and telomere length. Despite the large sample size for physical activity, the genetic instruments explaining the variation in exposure were limited, leading to wide 95% confidence intervals. Thirdly, the genetic variants for both exposures and outcomes were measured in the same individuals within the UK Biobank, which is known as one-sample MR. Although this ensures identical demographic characteristics (e.g., age, sex, socioeconomic background), it can inflate false positive rates (type 1 error). Fourthly, it is important to note that while MR studies provide valuable insights into the causal effects of lifetime exposures [66], they are limited in capturing variations attributed to non-genetic exposures that occur throughout an individual’s lifetime. Lastly, the generalizability of our findings may be limited to populations of European ancestry, and caution should be exercised when extrapolating the results to other racial groups. For instance, tobacco smoking is likely to be a potent risk factor across all racial groups.
Conclusion
The current MR study provided evidence supporting the association between lifetime smoking and insomnia with shorter telomere length, while no causal relationship was found between sleep duration and telomere length. However, further investigation is required to explore the potential causal role of physical activity in promoting longer telomere length, particularly through larger-scale MR analyses in the future. Consequently, adopting a healthy and active lifestyle may have potential health benefits.
Availability of data and materials
Our study used publicly available summary-level data of GWAS. The summary statistics for lifetime smoking is available at the website https://doi.org/https://doi.org/10.5523/bris.10i96zb8gm0j81yz0q6ztei23d. The summary statistics for sleep duration is available at the website https://sleep.hugeamp.org/downloads.html. The summary statistics for insomnia is available at the website https://ctg.cncr.nl/software/summary_statistics. The summary statistics for telomere length is available at the website https://figshare.com/s/caa99dc0f76d62990195 (Exclude 23and me). The summary statistics for alcohol use is available at the website https://genome.psych.umn.edu/index.php/GSCAN. The summary statistics of GWAS for type 2 diabetes, body mass index, and alcohol use disorder can be accessed at IEU OpenGWAS project (https://gwas.mrcieu.ac.uk/); Each IDs were presented in Additional file 1: Table S1. The code used in the current study can be obtained from the additional information.
References
Paneni F, Diaz Cañestro C, Libby P, Lüscher TF, Camici GG. The aging cardiovascular system: understanding it at the cellular and clinical levels. J Am Coll Cardiol. 2017;69(15):1952–67.
North BJ, Sinclair DA. The intersection between aging and cardiovascular disease. Circ Res. 2012;110(8):1097–108.
Ismail L, Materwala H, Al KJ. Association of risk factors with type 2 diabetes: a systematic review. Comput Struct Biotechnol J. 2021;19:1759–85.
Ewertz M, Christensen K, Engholm G, Kejs AM, Lund L, Matzen LE, Pfeiffer P, Storm HH, Herrstedt J. Trends in cancer in the elderly population in Denmark, 1980–2012. Acta Oncol (Stockholm, Sweden). 2016;55(Suppl 1):1–6.
Jylhävä J, Pedersen NL, Hägg S. Biological age predictors. EBioMedicine. 2017;21:29–36.
Barrett EL, Richardson DS. Sex differences in telomeres and lifespan. Aging Cell. 2011;10(6):913–21.
Diez Roux AV, Ranjit N, Jenny NS, Shea S, Cushman M, Fitzpatrick A, Seeman T. Race/ethnicity and telomere length in the Multi-Ethnic Study of Atherosclerosis. Aging Cell. 2009;8(3):251–7.
Topiwala A, Taschler B, Ebmeier KP, Smith S, Zhou H, Levey DF, Codd V, Samani NJ, Gelernter J, et al. Alcohol consumption and telomere length: Mendelian randomization clarifies alcohol’s effects. Mol Psychiatry. 2022;27(10):4001–8.
Barragán R, Ortega-AzorÃn C, Sorlà JV, Asensio EM, Coltell O, St-Onge MP, Portolés O, Corella D. Effect of physical activity, smoking, and sleep on telomere length: a systematic review of observational and intervention studies. J Clin Med. 2021;11(1):76.
Du M, Prescott J, Kraft P, Han J, Giovannucci E, Hankinson SE, De Vivo I. Physical activity, sedentary behavior, and leukocyte telomere length in women. Am J Epidemiol. 2012;175(5):414–22.
Fretts AM, Mete M, Howard BV, Best LG, Siscovick DS, Eilat-Adar S, Zhao J. Physical activity and telomere length in American Indians: the Strong Heart Study. Eur J Epidemiol. 2018;33(5):497–500.
Cassidy A, De Vivo I, Liu Y, Han J, Prescott J, Hunter DJ, Rimm EB. Associations between diet, lifestyle factors, and telomere length in women. Am J Clin Nutr. 2010;91(5):1273–80.
Ding LX, Zhang YH, Xu XZ, Zhang J, Ming SU, Liu D, Zhao ZY, Zhou Y, Zhang Q, et al. Association between physical activity and telomere length in a north Chinese population: a China suboptimal health cohort study. Biomed Environ Sci BES. 2018;31(5):394–8.
Steptoe A, Hamer M, Lin J, Blackburn EH, Erusalimsky JD. The longitudinal relationship between cortisol responses to mental stress and leukocyte telomere attrition. J Clin Endocrinol Metab. 2017;102(3):962–9.
Huzen J, Wong LS, van Veldhuisen DJ, Samani NJ, Zwinderman AH, Codd V, Cawthon RM, Benus GF, van der Horst IC, et al. Telomere length loss due to smoking and metabolic traits. J Intern Med. 2014;275(2):155–63.
Zhang C, Lauderdale DS, Pierce BL. Sex-specific and time-varying associations between cigarette smoking and telomere length among older adults. Am J Epidemiol. 2016;184(12):922–32.
Walters MS, De BP, Salit J, Buro-Auriemma LJ, Wilson T, Rogalski AM, Lief L, Hackett NR, Staudt MR, et al. Smoking accelerates aging of the small airway epithelium. Respir Res. 2014;15(1):94.
Cribbet MR, Carlisle M, Cawthon RM, Uchino BN, Williams PG, Smith TW, Gunn HE, Light KC. Cellular aging and restorative processes: subjective sleep quality and duration moderate the association between age and telomere length in a sample of middle-aged and older adults. Sleep. 2014;37(1):65–70.
Grieshober L, Wactawski-Wende J, Hageman Blair R, Mu L, Liu J, Nie J, Carty CL, Hale L, Kroenke CH, et al. A cross-sectional analysis of telomere length and sleep in the women’s health initiative. Am J Epidemiol. 2019;188(9):1616–26.
Tempaku P, Hirotsu C, Mazzotti D, Xavier G, Maurya P, Brietzke E, Belangero S, Poyares D, Bittencourt L, et al. Long sleep duration, insomnia, and insomnia with short objective sleep duration are independently associated with short telomere length. J Clin Sleep Med JCSM. 2018;14(12):2037–45.
Smith GD, Ebrahim S. ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22.
Klimentidis YC, Raichlen DA, Bea J, Garcia DO, Wineinger NE, Mandarino LJ, Alexander GE, Chen Z, Going SB. Genome-wide association study of habitual physical activity in over 377,000 UK Biobank participants identifies multiple variants including CADM2 and APOE. Int J Obes (London). 2018;42(6):1161–76.
Jansen PR, Watanabe K, Stringer S, Skene N, Bryois J, Hammerschlag AR, de Leeuw CA, Benjamins JS, Munoz-Manchado AB, et al. Genome-wide analysis of insomnia in 1,331,010 individuals identifies new risk loci and functional pathways. Nat Genet. 2019;51(3):394–403.
Dashti HS, Jones SE, Wood AR, Lane JM, van Hees VT, Wang H, Rhodes JA, Song Y, Patel K, et al. Genome-wide association study identifies genetic loci for self-reported habitual sleep duration supported by accelerometer-derived estimates. Nat Commun. 2019;10(1):1100.
Wootton RE, Richmond RC, Stuijfzand BG, Lawn RB, Sallis HM, Taylor GMJ, Hemani G, Jones HJ, Zammit S, et al. Evidence for causal effects of lifetime smoking on risk for depression and schizophrenia: a Mendelian randomisation study. Psychol Med. 2020;50(14):2435–43.
Codd V, Wang Q, Allara E, Musicha C, Kaptoge S, Stoma S, Jiang T, Hamby SE, Braund PS, et al. Polygenic basis and biomedical consequences of telomere length variation. Nat Genet. 2021;53(10):1425–33.
Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, Laurin C, Burgess S, Bowden J et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife 2018; 7.
Clarke L, Zheng-Bradley X, Smith R, Kulesha E, Xiao C, Toneva I, Vaughan B, Preuss D, Leinonen R, et al. The 1000 Genomes Project: data management and community access. Nat Methods. 2012;9(5):459–62.
Craig CL, Marshall AL, Sjöström M, Bauman AE, Booth ML, Ainsworth BE, Pratt M, Ekelund U, Yngve A, et al. International physical activity questionnaire: 12-country reliability and validity. Med Sci Sports Exerc. 2003;35(8):1381–95.
Codd V, Denniff M, Swinfield C, Warner SC, Papakonstantinou M, Sheth S, Nanus DE, Budgeon CA, Musicha C, et al. Measurement and initial characterization of leukocyte telomere length in 474,074 participants in UK Biobank. Nature Aging. 2022;2(2):170–9.
Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, Powell C, Vedantam S, Buchkovich ML, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015;518(7538):197–206.
Liu M, Jiang Y, Wedow R, Li Y, Brazel DM, Chen F, Datta G, Davila-Velderrain J, McGuire D, et al. Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Nat Genet. 2019;51(2):237–44.
Burgess S, Thompson SG. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. 2011;40(3):755–64.
Yao S, Zhang M, Dong SS, Wang JH, Zhang K, Guo J, Guo Y, Yang TL. Bidirectional two-sample Mendelian randomization analysis identifies causal associations between relative carbohydrate intake and depression. Nat Hum Behav. 2022;6:1569.
Palmer TM, Lawlor DA, Harbord RM, Sheehan NA, Tobias JH, Timpson NJ, Davey Smith G, Sterne JA. Using multiple genetic variants as instrumental variables for modifiable risk factors. Stat Methods Med Res. 2012;21(3):223–42.
Burgess S. Sample size and power calculations in Mendelian randomization with a single instrumental variable and a binary outcome. Int J Epidemiol. 2014;43(3):922–9.
Palmer TM, Sterne JA, Harbord RM, Lawlor DA, Sheehan NA, Meng S, Granell R, Smith GD, Didelez V. Instrumental variable estimation of causal risk ratios and causal odds ratios in Mendelian randomization analyses. Am J Epidemiol. 2011;173(12):1392–403.
Burgess S, Bowden J, Fall T, Ingelsson E, Thompson SG. Sensitivity analyses for robust causal inference from Mendelian randomization analyses with multiple genetic variants. Epidemiology. 2017;28(1):30–42.
Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan N, Thompson J. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med. 2017;36(11):1783–802.
Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25.
Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50(5):693–8.
Müezzinler A, Mons U, Dieffenbach AK, Butterbach K, Saum KU, Schick M, Stammer H, Boukamp P, Holleczek B, et al. Smoking habits and leukocyte telomere length dynamics among older adults: results from the ESTHER cohort. Exp Gerontol. 2015;70:18–25.
Bendix L, Thinggaard M, Fenger M, Kolvraa S, Avlund K, Linneberg A, Osler M. Longitudinal changes in leukocyte telomere length and mortality in humans. J Gerontol A Biol Sci Med Sci. 2014;69(2):231–9.
Zhao H, Han L, Chang D, Ye Y, Shen J, Daniel CR, Gu J, Chow WH, Wu X. Social-demographics, health behaviors, and telomere length in the Mexican American Mano a Mano Cohort. Oncotarget. 2017;8(57):96553–67.
Liang G, Schernhammer E, Qi L, Gao X, De Vivo I, Han J. Associations between rotating night shifts, sleep duration, and telomere length in women. PLoS ONE. 2011;6(8): e23462.
Jackowska M, Hamer M, Carvalho LA, Erusalimsky JD, Butcher L, Steptoe A. Short sleep duration is associated with shorter telomere length in healthy men: findings from the Whitehall II cohort study. PLoS ONE. 2012;7(10): e47292.
Ogawa EF, Leveille SG, Wright JA, Shi L, Camhi SM, You T. Physical activity domains/recommendations and leukocyte telomere length in U.S. adults. Med Sci Sports Exer. 2017;49(7):1375–82.
Tucker LA. Physical activity and telomere length in U.S. men and women: an NHANES investigation. Prev Med. 2017;100:145–51.
Loprinzi PD, Loenneke JP, Blackburn EH. Movement-based behaviors and leukocyte telomere length among US adults. Med Sci Sports Exerc. 2015;47(11):2347–52.
Soares-Miranda L, Imamura F, Siscovick D, Jenny NS, Fitzpatrick AL, Mozaffarian D. Physical activity, physical fitness, and leukocyte telomere length: the cardiovascular health study. Med Sci Sports Exerc. 2015;47(12):2525–34.
Stenbäck V, Mutt SJ, Leppäluoto J, Gagnon DD, Mäkelä KA, Jokelainen J, Keinänen-Kiukaanniemi S, Herzig KH. Association of physical activity with telomere length among elderly adults—the oulu cohort 1945. Front Physiol. 2019;10:444.
Weischer M, Bojesen SE, Nordestgaard BG. Telomere shortening unrelated to smoking, body weight, physical activity, and alcohol intake: 4,576 general population individuals with repeat measurements 10 years apart. PLoS Genet. 2014;10(3): e1004191.
Brandao CFC, Nonino CB, de Carvalho FG, Nicoletti CF, Noronha NY, San Martin R, de Freitas EC, Junqueira-Franco MVM, Marchini JS. The effects of short-term combined exercise training on telomere length in obese women: a prospective, interventional study. Sports Med Open. 2020;6(1):5.
Diman A, Boros J, Poulain F, Rodriguez J, Purnelle M, Episkopou H, Bertrand L, Francaux M, Deldicque L, et al. Nuclear respiratory factor 1 and endurance exercise promote human telomere transcription. Sci Adv. 2016;2(7): e1600031.
Dimauro I, Sgura A, Pittaluga M, Magi F, Fantini C, Mancinelli R, Sgadari A, Fulle S, Caporossi D. Regular exercise participation improves genomic stability in diabetic patients: an exploratory study to analyse telomere length and DNA damage. Sci Rep. 2017;7(1):4137.
Mason C, Risques RA, Xiao L, Duggan CR, Imayama I, Campbell KL, Kong A, Foster-Schubert KE, Wang CY, et al. Independent and combined effects of dietary weight loss and exercise on leukocyte telomere length in postmenopausal women. Obesity (Silver Spring). 2013;21(12):E549-554.
von Zglinicki T. Oxidative stress shortens telomeres. Trends Biochem Sci. 2002;27(7):339–44.
Jurk D, Wilson C, Passos JF, Oakley F, Correia-Melo C, Greaves L, Saretzki G, Fox C, Lawless C, et al. Chronic inflammation induces telomere dysfunction and accelerates ageing in mice. Nat Commun. 2014;2:4172.
van der Vaart H, Postma DS, Timens W, ten Hacken NH. Acute effects of cigarette smoke on inflammation and oxidative stress: a review. Thorax. 2004;59(8):713–21.
Irwin MR, Olmstead R, Carroll JE. Sleep disturbance, sleep duration, and inflammation: a systematic review and meta-analysis of cohort studies and experimental sleep deprivation. Biol Psychiat. 2016;80(1):40–52.
Nimmo MA, Leggate M, Viana JL, King JA. The effect of physical activity on mediators of inflammation. Diabetes Obes Metab. 2013;15(Suppl 3):51–60.
Teixeira-Lemos E, Nunes S, Teixeira F, Reis F. Regular physical exercise training assists in preventing type 2 diabetes development: focus on its antioxidant and anti-inflammatory properties. Cardiovasc Diabetol. 2011;10:12.
Daskalopoulou SS, Cooke AB, Gomez YH, Mutter AF, Filippaios A, Mesfum ET, Mantzoros CS. Plasma irisin levels progressively increase in response to increasing exercise workloads in young, healthy, active subjects. Eur J Endocrinol. 2014;171(3):343–52.
Rana KS, Arif M, Hill EJ, Aldred S, Nagel DA, Nevill A, Randeva HS, Bailey CJ, Bellary S, et al. Plasma irisin levels predict telomere length in healthy adults. Age (Dordr). 2014;36(2):995–1001.
Ludlow AT, Roth SM. Physical activity and telomere biology: exploring the link with aging-related disease prevention. J Aging Res. 2011;2011: 790378.
Soliman GA, Schooling CM. Causal association between mTOR-dependent EIF-4E and EIF-4A circulating protein levels and type 2 diabetes: a Mendelian randomization study. Sci Rep. 2020;10(1):15737.
Acknowledgements
None.
Funding
This study was supported by the Natural Science Foundation for Youth of China (82200269) and the Natural Science Foundation of China (82271812).
Author information
Authors and Affiliations
Contributions
LM conceived and designed the experiments. MC and ZW performed the experiments and analyzed the data. MC wrote the manuscript. HFX, PT, and WDL proofread the manuscript and data. All authors read and approved the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable. We used summary-level data from publicly available GWAS studies which have received ethical approval from their respective institutional review boards and informed consent from all participants. No administrative permissions were required to access the data.
Consent for publication
Not applicable.
Competing interests
All authors declare no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1. Table S1.
Descriptive information of the datasets included in the analyses. Table S2. Associations of single nucleotide polymorphisms for lifetime smoking index. Table S3. Associations of single nucleotide polymorphisms for continuant sleep duration. Table S4. Associations of single nucleotide polymorphisms for short sleep duration. Table S5. Associations of single nucleotide polymorphisms for long sleep duration.Table S6. Associations of single nucleotide polymorphisms for insomnia. Table S7. Associations of single nucleotide polymorphisms for moderate-to-vigorous physical activity. Table S8. Associations of single nucleotide polymorphisms for vigorous physical activity. Table S9. Associations of single nucleotide polymorphisms for strenuous sports or other exercises. Table S10. Post-hoc power calculations for our main IVW analyses on lifestyle factors and telomere length. Table S11. Estimates for the association between lifestyle factors and telomere length. Table S12. Heterogeneity and MR-Egger test for Horizontal pleiotropy. Table S13. Associations between genetic liability to lifestyle factors and telomere length following exclusion of outlier SNPs identified by MR-PRESSO.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Chen, M., Wang, Z., Xu, H. et al. Association between modifiable lifestyle factors and telomere length: a univariable and multivariable Mendelian randomization study. J Transl Med 22, 160 (2024). https://doi.org/10.1186/s12967-024-04956-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12967-024-04956-8