Skip to main content
Download PDF

Genetic association of telomere length, obesity and tobacoo smoking with idiopathic pulmonary fibrosis risk

BMC Public Health volume 23, Article number: 868 (2023) Cite this article

Abstract

Background

Due to the inadequacy of published evidence, association of telomere length (TL), obesity and tobacco smoking with idiopathic pulmonary fibrosis (IPF) remains unclear. The aim of the study was to explore whether these exposures genetically affected the risk of the disease.

Methods

Genetic variants from genome-wide association studies for TL, body mass index (BMI), body fat percentage (BFP) and tobacco smoking (including maternal smoking) were used as instrumental variables. Inverse-variance weighted were mainly adopted to determine the genetic association of these exposures with IPF. All analyses were conducted by R-software (version 3.6.1).

Results

Firstly, longer TL was associated with the decreased risk of IPF (OR = 0.475 per SD increase in TL, 95%CI = 0.336 ~ 0.670, P<0.001). Secondly, higher levels of BMI and BFP were related to the increased risk of the disease (OR = 1.425 per SD increase in BMI level, 95%CI = 1.114 ~ 1.823, P = 0.005; OR = 1.702 per SD increase in BFP level, 95%CI = 1.202 ~ 2.409, P = 0.003). Thirdly, maternal smoking was implicated in the increased risk of the disease (OR = 13.183 per SD increase in the prevalence of maternal smoking, 95%CI = 1.820 ~ 95.484, P = 0.011).

Conclusion

TL should be a genetic risk factor for IPF. Obesity and exposure to tobacco smoking as a fetus might also contribute to the development of this fibrotic diseases. These findings should be verified by future studies.

Peer Review reports

Introduction

Idiopathic pulmonary fibrosis (IPF) is a chronic, progressive interstitial lung disease that tends to occur in middle-aged and elderly people [1]. According to statistics, the prevalence rate of the disease in the global population is about 2 to 29 per 100,000 each year [2]. In China, the number of IPF patients is conservatively estimated at least about 500,000 [3, 4]. More importantly, IFP is a fatal disorder with unsatisfactory prognosis [5,6,7]. So, exploring the risk factors and pathogenesis of the disease is of great significance for the control of the disease, and related research is being carried out.

Telomere length (TL) is a biomarker closely related to aging [8]. Previous observational studies suggested that shortening TL predicted the increased risk and poor prognosis of pulmonary fibrosis [9, 10]. One Mendelian randomization (MR) study which adopted genetic variants as proxies for exposure and outcome reported a casual relationship between TL and IPF for the first time [11]. However, this can only be regarded as a preliminary conclusion due to the small sample size of IPF and insufficient power of test, and needed to be validated by further research.

As we all know, a great number of factors have been shown to potentially affect TL. One recent systematic review and meta-analysis including 84 studies suggested that TL was shorter among ever-smokers compared to never-smokers, which may imply mechanisms linking tobacco smoking to ageing-related disease [12]. One published MR study revealed that higher genetically-predicted alcohol use disorder was associated with shorter TL, and non-linear analyses indicated a potential threshold relationship between alcohol and TL [13]. Two meta-analyses also revealed the possibility of an effect of paediatric obesity on telomerase length, although the current evidence had not been sufficient yet [14, 15].

Importantly, an increasing number of observational studies focused on the potential effects of these telomere-associated factors (i.e. tobacco smoking, alcohol drinking and obesity) on the risk of IPF, and provided some interesting findings [9, 10]. Considering the limitations of observational studies, these findings need to be validated at the genetic level by MR analysis.

Taken together, using the most updated and available genome-wide association study (GWAS) data for IPF, TL and several telomere-associated factors (i.e. tobacco smoking, alcohol drinking and obesity), the present MR study aimed to explore and verify the genetic effects of TL and these telomere-associated factors on the risk of IPF, so as to provide new knowledge about the pathogenesis of this fibrotic disease. In this process, we adopted heavy smoking, smoking initiation, age of smoking initiation, smoking cessation and maternal smoking to interpret different stages and levels of tobacco exposure, and also used body mass index (BMI) and body fat percentage (BFP) to more comprehensively describe obesity and fat accumulation.

Methods

Diagram of MR framework in the study was showed in Fig. 1.

Fig. 1
figure 1

Diagram of Mendelian randomization framework in the study

Full size image

Data sources for the exposures

The summary data for TL was obtained from a published GWAS which characterized the genetic architecture of naturally occurring variation in leukocyte TL using quantitative PCR assay and identified causal links between TL and biomedical phenotypes in 472,174 well-characterized UK Biobank participants [16]. An increase or decrease in TL indicated an increase or decrease of one or more standard deviations (SDs) in the mean base pairs.

The data for heavy smoking, smoking initiation, age of smoking initiation, smoking cessation, heavy drinking was collected from a GWAS and Sequencing Consortium of Alcohol and Nicotine use (GSCAN) study which explored the genetic etiology of tobacco and alcohol use in up to 1.2 million European participants [17]. Heavy smoking and heavy drinking were defined as daily tobacco use and weekly alcohol use respectively. Smoking initiation was a binary phenotype indicating whether an individual had ever smoked regularly, while smoking cessation was also a binary variable contrasting current versus former smokers [17].

The data for maternal smoking and BFP was included from studies of MRC Integrative Epidemiology Unit (MRC-IEU) using 397,732 to 454,633 UK Biobank participants from IEU Open GWAS Project (https://gwas.mrcieu.ac.uk/), and the IDs were separately ukb-b-17,685 and ukb-b-8909. Maternal smoking was defined as maternal tobacco smoking around birth. The data for BMI was collected from a meta-analysis of Genetic Investigation of ANthropometric Traits (GIANT) studies for the trait in up to 700,000 participants of European ancestry [18].

Characteristics of these exposures were shown in Table 1.

Table 1 Sources and characteristics of GWAS summary data in the study
Full size table

Data source for the outcome

The summary data for IPF was obtained from an available and updated FinnGen database (DF8) including 340,596 European participants with 1,812 cases and 338,784 controls [19]. The outcome was searched from hospital discharge records or death records using International Classification of Diseases 10th edition (ICD−10) codes (ICD−10 code: J84.1). Characteristics of the outcome were shown in Table 1.

Selection of genetic instruments

Single nucleotide polymorphisms (SNPs) are one of the most common types of genetic variants in humans and are used as instrumental variables to replace traits in MR studies [20]. The suitable SNPs were extracted from the above exposure-associated GWASs according to the following criteria: (1) A genome-wide significance-associated P value was less than 5 × 10− 8. If the number of available SNPs was insufficient, a more relaxed threshold (5 × 10− 6) was adopted. (2) A F-statistic should be higher than 10. (3) Remove SNPs with linkage disequilibrium [LD] (a clumping window of 10 MB and a r2 value less than 0.001). (4) All SNPs associated with the outcome and potential confounders were directly excluded using the PhenoScanner Pheno Scanner V2, and these confounders included air pollution, viral infection, autoimmune deficiency and gastroesophageal reflux disease [21].

Mendelian randomization analysis

Random-effect inverse-variance weighted (IVW) and MR-Egger were adopted to determine the causal association of TL and the associated factors with the risk of IPF. The IVW result is the slope of a weighted regression of SNP-outcome effect on SNP-exposure effect with an intercept equal to zero [22]. In MR studies, the intercept is the longitudinal coordinate of the line as it crosses the Y-axis. While the x-axis represents the SNP-exposure effect and the y-perimeter represents the SNP-outcome effect in a scatter plot, an intercept of zero means that when the SNP-exposure effect is zero, the SNP-outcome effect is also zero. It is widely accepted that the IVW is the most reliable method compared to the other methods. The MR-Egger complements the IVW, as the MR-Egger can be applied to a wider range of conditions and produce more conservative results, but is less efficient [23]. All instrumental variables in MR-Egger can have pleiotropy, but these pleiotropies cannot affect the association of instrumental variables with exposure. Here, pleiotropy refers primarily to horizontal pleiotropy, meaning that one genetic variant as an instrumental variable also affects outcome through other pathways independent of the exposure [24]. In the study, the correlation was considered statistically significant when the IVW P-value was less than 0.05 and the MR-Egger results were at least in the same direction as the IVW results [25]. In addition, scatter plots and forest plots were used to visualize the results obtained from the IVW and MR-Egger methods.

Four main methods were adopted for sensitivity analysis to control the quality of data and results. (1) Cochran’s Q test was adopted to determine underlying heterogeneity. (2) MR-Egger intercept test was used to measure horizontal pleiotropy. (3) Leave-one-out test evaluated whether the estimate was led by one single SNP, which was also related to pleiotropy. (4) Funnel plot assessed directional pleiotropy, which was similar to publication bias in meta-analysis. In MR studies, directional pleiotropy mainly refers to horizontal pleiotropy.

In addition, MR-Pleiotropy Residual Sum and Outlier method (MR-PRESSO) was also adopted to find outliers and correct horizontal pleiotropy. The method can also provide causal estimates including and removing outiers. MR-PRESSO is more accurate than the IVW and MR-Egger, when the number of genetic variants with horizontal pleiotropy is less than 10% of the total [24].

All analyses were conducted using TwoSampleMR package (version 0.4.25) and MRPRESSO package (version 1.0) in R software (version 3.6.1).

Results

Overview of the analysis

In Table 1, a total of 19 to 443 SNPs were separately extracted to genetically predict TL and its 8 associated factors. All SNPs met the inclusion criteria mentioned in Table 1, and were listed in supplemental Table 1 . The IVW and MR-Egger results as well as the Cochran’s Q and MR-Egger intercept results were showed in Table 2. Total R2, F-statistic and test power for each estimate was listed in supplemental Table 1 . MR-PRESSO results were showed in supplemental Table 3. The scatter plots, forest plots, Leave-one-out tests and funnel plots were showed in Figs. 2, 3 and 4 and supplemental Fig. 1 to 9 .

Fig. 2
figure 2

Mendelian randomization effect of telomere length on idiopathic pulmonary fibrosis risk (excluding the SNPs associated with smoking or obesity)

Full size image
Fig. 3
figure 3

Forest plots for effect of tobacco smoking and alcohol drinking on idiopathic pulmonary fibrosis risk (excluding the SNPs associated with telomere length or obesity in Figure 3C and 3E)

Full size image
Fig. 4
figure 4

Mendelian randomization effect of body mass index and body fat percentage on idiopathic pulmonary fibrosis risk (excluding the SNPs associated with telomere length or smoking)

Full size image
Table 2 Effect of telomere length and associated factors on idiopathic pulmonary fibrosis risk
Full size table

Effect of TL on the risk of IPF

Both IVW and MR-Egger results reported that longer TL was associated with the decreased risk of IPF (OR = 0.459 per SD increase in TL, 95%CI = 0.327 ~ 0.643, P<0.001; OR = 0.284 per SD increase in TL, 95%CI = 0.145 ~ 0.556, P<0.001). Pre-analyses showed that some markers of obesity and smoking can also affect the risk of IPF. Therefore, after excluding the obesity and smoking-related SNPs, the analysis was re-conducted. The updated IVW and MR-Egger results still reported that longer TL was associated with the decreased risk of the disease (OR = 0.475 per SD increase in TL, 95%CI = 0.336 ~ 0.670, P<0.001; OR = 0.298 per SD increase in TL, 95%CI = 0.152 ~ 0.586, P<0.001). These results were visualized in Fig. 2 and supplemental Fig. 1 .

For the updated results, P value for Cochran’s Q was 0.665 and P value for MR-Egger intercept was 0.121, which indicated that there was no significant heterogeneity and horizontal pleiotropy (Table 2). The “leave-one-SNP-out” analysis and funnel plot also found no horizontal pleiotropy and bias in the analysis (supplemental Fig. 1 ). The MR-PRESSO did not detect any outliers (P = 0.658), and its causal estimate was consistent with the results from the IVW and MR-Egger (P < 0.001) (supplemental Table 3 ).

These findings indicated that the risk of IPF was reduced by approximately 50–70% for each SD increase in TL.

Effect of tobacco smoking and alcohol drinking on the risk of IPF

Five variables (i.e. heavy smoking, smoking initiation, age of smoking initiation, smoking cessation and maternal smoking) were used to represent tobacco smoking in a multifaceted way. Heavy drinking was adopted to represent alcohol drinking. Their definitions had been mentioned above.

The IVW results revealed that older age of smoking initiation was associated with the decreased risk of IPF (OR = 0.136 per SD increase in age of smoking initiation, 95%CI = 0.040 ~ 0.464, P = 0.001). However, the MR-Egger results were not statistically significant and were in the opposite direction (OR = 102.692 per SD increase in age of smoking initiation, 95%CI = 0.049 ~ 2.136e + 05, P = 0.251). These results were showed in Table 2.

The IVW results revealed that maternal smoking was associated with the increased risk of the disease (OR = 15.487 per SD increase in the prevalence of maternal smoking, 95%CI = 2.679 ~ 89.543, P = 0.002). And, the MR-Egger results were in the same direction as the IVW results (OR = 2.995 per SD increase in the prevalence of maternal smoking, 95%CI = 0.001 ~ 2.446e + 04, P = 0.812). After removing the TL and obesity-related SNPs, the updated results from the IVW and MR-Egger did not change significantly (OR = 13.183 per SD increase in the prevalence of maternal smoking, 95%CI = 1.820 ~ 95.484, P = 0.011; OR = 1.228 per SD increase in the prevalence of maternal smoking, 95%CI = 6.682e−05 ~ 2.256e + 04, P = 0.968). These results were showed in Table 2, and were also visualized in Fig. 3, supplemental Fig. 6 .

For the updated results of the maternal smoking, the Cochran’s Q, MR-Egger intercept, “leave-one-SNP-out” analysis and funnel plot did not find any heterogeneity and horizontal pleiotropy (Table 2, supplemental Fig. 6 ). The MR-PRESSO reported no significant outliers (P = 0.451), and its causal estimates were consistent with the results from the IVW (P = 0.015) (supplemental Table 3 ).

In addition, the study did not find any association of heavy smoking, smoking initiation, smoking cessation and heavy drinking with the risk of IPF.

These findings indicated that for every SD increase in the prevalence of maternal smoking, the risk of IPF was likely to increase more than tenfold. Given the wide 95%CIs, these results should be further verified.

Effect of obesity on the risk of IPF

BMI and BFP were used as the markers of obesity in the study. After removing the TL and smoking-associated SNPs, the MR-Egger method reported that higher levels of BMI and BFP were associated with the increased risk of IPF (OR = 1.425 per SD increase in BMI level, 95%CI = 1.114 ~ 1.823, P = 0.005; OR = 1.702 per SD increase in BFP level, 95%CI = 1.202 ~ 2.409, P = 0.003). The MR-Egger reported similar results (OR = 3.200 per SD increase in BMI level, 95%CI = 1.586 ~ 6.455, P = 0.001; OR = 4.603 per SD increase in BFP level, 95%CI = 1.424 ~ 14.880, P = 0.011). The results were showed in Table 2; Fig. 4, supplemental Figs. 8 and 9 .

All P values for Cochran’s Q and MR-Egger intercept were larger than 0.05 (Table 2). The “leave-one-SNP-out” analyses and funnel plots did not find horizontal pleiotropy and bias in the analysis (supplemental Figs. 8 and 9 ). The MR-PRESSO did not detect any outliers (P = 0.200, P = 0.153), and their causal estimates also reported the causal association of BMI and BFP with the disease (P = 0.005, P = 0.003) (supplemental Table 3 ).

These findings indicated that the risk of IPF was increased by approximately 40–70% for each SD increase in BMI or BFP level.

Discussion

As mentioned above, one previous MR study had preliminarily reported a causal relationship between TL and the risk of IPF [11]. However, all current summary data on IPF (including the data used in the above study) were limited by small sample sizes and cannot be resolved in the short term. Therefore, iterative validation of this correlation was necessary. The present study adopted another set of summary data for the fibrotic disease, and provided a result that TL was associated with the risk of IPF (Detailed data was listed in the result section). This finding was consistent with that previous study, and further confirmed the role of telomere in the pathogenesis of IPF.

The study explored the causal relationship between several types of smoking history and the risk of IPF. Although most correlations were not statistically significant, there were still some interesting findings that may be useful for future studies. Specifically, at least in the IVW results, participants whose mothers had maternal smoking may be subject to a higher risk of the disease (Detailed data was listed in the result section). This finding was supported by the MR-PRESSO, and the MR-Egger result also had a same direction as the IVW result. These results were consistent with one important medical truism that many organs of the fetus (including the lungs) were at a developmental stage and were susceptible to damage from tobacco. In addition, genetic mechanisms were thought to play an important role in the pathogenesis of IPF, which partly supported the results in the study [26].

The study also found that obesity may significantly increase the risk of developing IPF (Detailed data was listed in the result section). A recent meta-analysis showed that BMI might be useful to predict mortality, disease progression and treatment-related toxicity in IPF [27]. Another MR study showed that gastro-esophageal reflux disease may increase the risk of IPF, but after adjustment for BMI, these associations were significantly attenuated. Therefore, the investigators concluded that treating obesity, but not gastro-esophageal reflux disease, was likely to have a protective effect against IPF [28]. So, these results were consistent with the present study.

In two-sample MR studies, sample overlap can lead to biased effect estimates [29]. However, in some cases, sample overlap is difficult to avoid. Researchers can only assess the degree of sample overlap by calculating the overlap percentage, and control it by selecting other data. Thus, sample overlap is still one of the main potential limitations in this type of MR studies. In the present study, the summary data of the exposures were mainly from UK Biobank and other institutions, while the summary data of the outcome (i.e. IPF) came from FinnGen. Since UK Biobank and FinnGen were independent in summary data, we considered the degree of sample overlap in the present study to be acceptable and not to have an significant impact on the results.

At present, the sample sizes of all the summary data on IPF were relatively small, which may be related to the insidious onset of the disease and the difficulty in diagnosing it. In the present study, the most updated summary data for IPF had been selected, but the small sample size remained a limitation. This was a reasonable explanation for the very wide 95%CIs from some of the results in this study [30]. The powers of test had been calculated and listed them in the supplemental Table 2 . As the sample size of IPF data was limited by objective conditions, researchers can adopt summary data from different sources for the same trait to conduct their analyses and validate the results obtained against each other, which may be a realistic solution. However, it must be noted that MR can be difficult to interpret. Its methodology was extremely difficult to implement for non-specialists. Hence, only the researchers with the appropriate co-workers could conduct their analysis and validate the results.

In MR studies, weak instrumental variable bias is an issue that cannot be avoided. A weak instrumental variable is a genetic variant that is less potent in explaining exposure [31]. This instrument is associated with the exposure, but the strength of this association is not very high, so it is fundamentally different from the null instrumental variable. In general, the main cause of weak instrumental variable bias is inadequate sample size. In the present study, we calculated F-statistic, and any instrumental variable with an F-statistic less than 10 was excluded, thus reducing the impact of this bias on the study.

The present study performed several sensitivity analyses to detect and control the horizontal pleiotropy and heterogeneity, such as Cochran’s Q, MR-Egger intercept, and so on. These may help to ensure the reliability of the conclusions in this study. In addition, with the increasing size and number of GWASs, individual SNP is increasingly found to associate with multiple traits. Vertical pleiotropy is one of the scenarios, and it means that exposure is also associated with other traits, and association with outcome can be established through these traits [32]. Vertical pleiotropy is very common. One mainstream view is that vertical pleiotropy does not affect the results, and does not need to be accounted for in MR studies. However, this view is also not fully substantiated. So, sensitivity analyses were also necessary in the present study.

Potential effect of binary exposures on causal estimates could be a potential limitation in MR studies. When a binary exposure is a dichotomization of a continuous risk factor, its genetic variants (instrument variables) might violate the exclusion restriction assumption, because the genetic variants are able to affect the outcome through the continuous risk factor, though the binary exposure does not change [33]. In the present study, the exposures with potential positive results, such as TL, BMI, BFP and maternal smoking, were not classic binary exposures. So, the limitation may not affect our main conclusions. However, some of the other exposures in the study were binary exposures, the results of which may be affected by this limitation and need to be validated by further research.

Telomere length can decrease progressively with age. So, when exploring the effect of telomere length on one disease, the exposure and outcome should be from the same age group population. In the present study, although both sets of summary data (TL and IPF) were from the adult population, there was still a significant difference in age, which was an important limitation of the present study. Therefore, further research was necessary to explore this issue as the study data allowed.

In conclusion, the present study confirmed the causal association of TL with the risk of IPF. The present study also provided some interesting genetic evidence to prove that obesity and exposure to tobacco smoking as a fetus might also contribute to the development of the fibrotic diseases. However, it was important to note that the evidence on tobacco smoking exposure was inadequate and should be verified by future studies.

Data Availability

Data was already available in the text and supplementary materials.

Abbreviations

IPF:

Idiopathic pulmonary fibrosis

TL:

Telomere length

MR:

Mendelian randomization

GWAS:

Genome-wide association study

BMI:

Body mass index

BFP:

Body fat percentage

GSCAN:

GWAS and Sequencing Consortium of Alcohol and Nicotine use

MRC-IEU:

MRC Integrative Epidemiology Unit

GIANT:

Genetic Investigation of ANthropometric Traits

SNP:

Single nucleotide polymorphisms

IVW:

Inverse-variance weighted

MR-PRESSO:

MR-Pleiotropy Residual Sum and Outlier method

References

  1. Richeldi L, Collard HR, Jones MG. Idiopathic pulmonary fibrosis. Lancet. 2017;389:1941–52.

    Article  PubMed  Google Scholar 

  2. Maher TM, Bendstrup E, Dron L, Langley J, Smith G, Khalid JM, Patel H, Kreuter M. Global incidence and prevalence of idiopathic pulmonary fibrosis. Respir Res. 2021;22:197.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Huang H, Peng X, Zhong C. Idiopathic pulmonary fibrosis: the current status of its epidemiology, diagnosis, and treatment in China. Intractable Rare Dis Res. 2013;2:88–93.

    PubMed  PubMed Central  Google Scholar 

  4. Xie B, Ren Y, Geng J, He X, Ban C, Wang S, Jiang D, Luo S, Chen Q, Liu M, Feng R, Zhao L, Dai H, Wang C. Idiopathic pulmonary fibrosis Registry China study (PORTRAY): protocol for a prospective, multicentre registry study. BMJ Open. 2020;10:e036809.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Sgalla G, Iovene B, Calvello M, Ori M, Varone F, Richeldi L. Idiopathic pulmonary fibrosis: pathogenesis and management. Respir Res. 2018;19:32.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Moss BJ, Ryter SW, Rosas IO. Pathogenic mechanisms underlying idiopathic pulmonary fibrosis. Annu Rev Pathol. 2022;17:515–46.

    Article  PubMed  Google Scholar 

  7. Wu W, Zhang G, Qiu L, Liu X, Zhou S, Wu J. Contribution of Adiponectin/Carnitine Palmityl Transferase 1A-Mediated Fatty Acid Metabolism during the Development of Idiopathic Pulmonary Fibrosis. Oxid Med Cell Longev. 2022 Aug 17;2022:5265616.

  8. Vaiserman A, Krasnienkov D. Telomere length as a marker of Biological Age: state-of-the-Art, Open Issues, and future perspectives. Front Genet. 2021;11:630186.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Bae W, Lee CH, Lee J, Kim YW, Han K, Choi SM. Impact of smoking on the development of idiopathic pulmonary fibrosis: results from a nationwide population-based cohort study. Thorax. 2022;77:470–6.

    Article  PubMed  Google Scholar 

  10. Andersson M, Blanc PD, Torén K, Järvholm B. Smoking, occupational exposures, and idiopathic pulmonary fibrosis among swedish construction workers. Am J Ind Med. 2021;64:251–7.

    Article  PubMed  Google Scholar 

  11. Duckworth A, Gibbons MA, Allen RJ, Almond H, Beaumont RN, Wood AR, Lunnon K, Lindsay MA, Wain LV, Tyrrell J, Scotton CJ. Telomere length and risk of idiopathic pulmonary fibrosis and chronic obstructive pulmonary disease: a mendelian randomisation study. Lancet Respir Med. 2021;9:285–94.

    Article  CAS  PubMed  Google Scholar 

  12. Astuti Y, Wardhana A, Watkins J, Wulaningsih W, PILAR Research Network. Cigarette smoking and telomere length: a systematic review of 84 studies and meta-analysis. Environ Res. 2017 Oct;158:480–9.

  13. Topiwala A, Taschler B, Ebmeier KP, Smith S, Zhou H, Levey DF, Codd V, Samani NJ, Gelernter J, Nichols TE, Burgess S. Alcohol consumption and telomere length: mendelian randomization clarifies alcohol’s effects. Mol Psychiatry. 2022 Oct;27(10):4001–8.

  14. Kahrizi MS, Patra I, Jalil AT, Achmad H, Alesaeidi S, Al-Gazally ME, Alesaeidi S. Leukocyte telomere length and obesity in children and adolescents: a systematic review and meta-analysis. Front Genet. 2022 Sep;7:13:861101.

  15. Lin L, Qin K, Chen D, Lu C, Chen W, Guo VY. Systematic review and meta-analysis of the association between paediatric obesity and telomere length. Acta Paediatr. 2021 Oct;110(10):2695–703.

  16. Codd V, Wang Q, Allara E, Musicha C, Kaptoge S, Stoma S, Jiang T, Hamby SE, Braund PS, Bountziouka V, Budgeon CA, Denniff M, Swinfield C, Papakonstantinou M, Sheth S, Nanus DE, Warner SC, Wang M, Khera AV, Eales J, Ouwehand WH, Thompson JR, Di Angelantonio E, Wood AM, Butterworth AS, Danesh JN, Nelson CP, Samani NJ. Polygenic basis and biomedical consequences of telomere length variation. Nat Genet. 2021 Oct;53(10):1425–33.

  17. Liu M, Jiang Y, Wedow R, Li Y, Brazel DM, Chen F, Datta G, Davila-Velderrain J, McGuire D, Tian C, 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 Feb;51(2):237–44.

  18. Yengo L, Sidorenko J, Kemper KE, Zheng Z, Wood AR, Weedon MN, Frayling TM, Hirschhorn J, Yang J, Visscher PM, GIANT Consortium. ;. Meta-analysis of genome-wide association studies for height and body mass index in ∼700000 individuals of European ancestry. Hum Mol Genet 2018 Oct 15;27(20):3641–9.

  19. Kurki MI, Karjalainen J, Palta P, Sipilä TP, Kristiansson K, Donner K, Reeve MP, Laivuori H, Aavikko M, Kaunisto MA et al. FinnGen: Unique genetic insights from combining isolated population and national health register data. medRxiv. 2022.

  20. Emdin CA, Khera AV, Kathiresan S. Mendelian Randomization JAMA. 2017;318:1925–6.

    Article  PubMed  Google Scholar 

  21. Richeldi L, Collard HR, Jones MG. Idiopathic pulmonary fibrosis. Lancet 2017 May 13;389(10082):1941–52.

  22. Birney E. Mendelian randomization. Cold Spring Harb Perspect Med. 2022 May;17(4):a041302.

  23. Burgess S, Thompson SG. Interpreting findings from mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017 May;32(5):377–89.

  24. Erbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and diseases. V Nat Genet. 2018 May;50(5):693–8.

  25. Ong JS, MacGregor S. Implementing MR-PRESSO and GCTA-GSMR for pleiotropy assessment in mendelian randomization studies from a practitioner’s perspective. Genet Epidemiol. 2019 Sep;43(6):609–16.

  26. Bonham CA, Noth I. Idiopathic pulmonary fibrosis genetic risk factors, function, and mechanisms? The Concepts are starting to gel. Am J Respir Crit Care Med. 2022 Nov;15(10):1197–9.

  27. Zinellu A, Carru C, Pirina P, Fois AG, Mangoni AA. A Systematic Review of the Prognostic Significance of the Body Mass Index in Idiopathic Pulmonary Fibrosis. J Clin Med. 2023 Jan 7;12(2):498.

  28. Cotton C, Alton P, Hughes DM, Zhao SS. Genetic liability to gastro-esophageal reflux disease, obesity, and risk of idiopathic pulmonary fibrosis. Respir Investig. 2023 Mar;16(3):335–8.

  29. Burgess S, Davies NM, Thompson SG. Bias due to participant overlap in two-sample mendelian randomization. Genet Epidemiol. 2016 Nov;40(7):597–608.

  30. Burgess S. Sample size and power calculations in mendelian randomization with a single instrumental variable and a binary outcome. Int J Epidemiol. 2014 Jun;43(3):922–9.

  31. Burgess S, Small DS, Thompson SG. A review of instrumental variable estimators for mendelian randomization. Stat Methods Med Res. 2017 Oct;26(5):2333–55.

  32. Yang Q, Sanderson E, Tilling K, Borges MC, Lawlor DA. Exploring and mitigating potential bias when genetic instrumental variables are associated with multiple non-exposure traits in mendelian randomization. Eur J Epidemiol. 2022 Jul;37(7):683–700.

  33. Burgess S, Labrecque JA. Mendelian randomization with a binary exposure variable: interpretation and presentation of causal estimates. Eur J Epidemiol. 2018 Oct;33(10):947–52.

Download references

Acknowledgements

Not applicable.

Funding

The project was supported by the Key Scientific and Technological Project of Henan Province (NO.222102310169).

Author information

Authors and Affiliations

  1. Department of Geriatrics Medicine, Henan Provincial People’s Hospital, Zhengzhou University, Zhengzhou, China

    Wenjuan Wu

  2. Stem cell program of clinical research center, People’s Hospital of Zhengzhou University and Henan Provincial People’s Hospital, Zhengzhou, China

    Chenghai Li

  3. Department of Thoracic Surgery, Henan Provincial People’s Hospital, Zhengzhou, China

    Xiaoming Zhu

  4. Department of Geriatrics Medicine, Henan Provincial People’s Hospital, Zhengzhou, China

    Xueya Liu

  5. Department of Respiratory and Critical Care Medicine, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China

    Ping Li

  6. Department of Geriatrics Medicine, Henan Provincial People’s Hospital, Zhengzhou, China

    Ruijie Wan

  7. Department of Traditional Chinese Medicine, Zhengzhou Shuqing Medical College, Zhengzhou, China

    Xinhui Wu

  8. Translational Research Institute, Henan Provincial People’s Hospital, Academy of Medical Science, Zhengzhou University, Zhengzhou, China

    Song Chen

Authors
  1. Wenjuan Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Chenghai Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Xiaoming Zhu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Xueya Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Ping Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Ruijie Wan
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Xinhui Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Song Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

WW, LC and CS conducted the study design. WW, LC, ZX, LX, LP, WR, WX and CS performed the literature research, data acquisition/collation and data analysis. WW, LC and CS preparated the manuscript.

Corresponding authors

Correspondence to Wenjuan Wu or Song Chen.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethics approval and consent to participate

The study was approved by the Ethics Committee of Henan Provincial People’s Hospital, and was performed in accordance with the Declaration of Helsinki. Informed consents were obtained from all subjects and/or their legal guardians.

Consent for publication

Not applicable.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, W., Li, C., Zhu, X. et al. Genetic association of telomere length, obesity and tobacoo smoking with idiopathic pulmonary fibrosis risk. BMC Public Health 23, 868 (2023). https://doi.org/10.1186/s12889-023-15733-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-023-15733-5

Keywords

BMC Public Health

ISSN: 1471-2458

Contact us