Skip to main content
Download PDF

Diabetes mellitus and idiopathic pulmonary fibrosis: a Mendelian randomization study

BMC Pulmonary Medicine volume 24, Article number: 142 (2024) Cite this article

Abstract

Background

The question as to whether or not diabetes mellitus increases the risk of idiopathic pulmonary fibrosis (IPF) remains controversial. This study aimed to investigate the causal association between type 1 diabetes (T1D), type 2 diabetes (T2D), and IPF using Mendelian randomization (MR) analysis.

Methods

We used two-sample univariate and multivariate MR (MVMR) analyses to investigate the causal relationship between T1D or T2D and IPF. We obtained genome-wide association study (GWAS) data for T1D and T2D from the European Bioinformatics Institute, comprising 29,652 T1D samples and 101,101 T1D single nucleotide polymorphisms (SNPs) and 655,666 T2D samples and 5,030,727 T2D SNPs. We also used IPF GWAS data from the FinnGen Biobank comprising 198,014 IPF samples and 16,380,413 IPF SNPs. All cases and controls in these datasets were derived exclusively from European populations. In the univariate MR analysis, we employed inverse variance-weighted (IVW), weighted median (WM), and MR-Egger regression methods. For the MVMR analysis, we used the multivariate IVW method primarily, and supplemented it with multivariate MR-Egger and multivariate MR- least absolute shrinkage and selection operator methods. Heterogeneity tests were conducted using the MR-IVW and MR-Egger regression methods, whereas pleiotropic effects were assessed using the MR-Egger intercept. The results of MR and sensitivity analyses were visualized using forest, scatter, leave-one-out, and funnel plots.

Results

Univariate MR revealed a significant causal relationship between T1D and IPF (OR = 1.118, 95% CI = 1.021–1.225, P = 0.016); however, no significant causal relationship was found between T2D and IPF (OR = 0.911, 95% CI = 0.796–1.043, P = 0.178). MVMR analysis further confirmed a causal association between T1D and IPF (OR = 1.133, 95% CI = 1.011–1.270, P = 0.032), but no causal relationship between T2D and IPF (OR = 1.009, 95% CI = 0.790–1.288, P = 0.950). Sensitivity analysis results validated the stability and reliability of our findings.

Conclusion

Univariate and multivariate analyses demonstrated a causal relationship between T1D and IPF, whereas no evidence was found to support a causal relationship between T2D and IPF. Therefore, in clinical practice, patients with T1D should undergo lung imaging for early detection of IPF.

Peer Review reports

Introduction

Idiopathic pulmonary fibrosis (IPF) is a chronic fibrotic interstitial pneumonia of unknown etiology characterized by cough, dyspnea, and a progressive decline in lung function that primarily affects the elderly population [1]. The radiological and histological features are predominantly consistent with usual interstitial pneumonia [1]. According to a statistical analysis of 12 countries, the incidence and prevalence rates of IPF range from 0.09–1.30/10,000 and 0.33–4.51/10,000, respectively, exhibiting a sharp increase with advancing age [2]. Globally, the number of patients with IPF is increasing, potentially due to factors such as an aging demographic, heightened disease awareness, and advancements in diagnostic technologies [2]. Patients diagnosed with IPF may have a poorer prognosis than those with many cancers that affect similar populations [3], as evidenced by a median survival rate of only 3–5 years in the absence of treatment [4]. Age, sex, and comorbidities significantly affect the clinical outcomes of patients with IPF, with mortality associated with congestive heart failure (CHF), diabetes mellitus (DM), and cancer, whereas hospitalization is linked to CHF and chronic obstructive pulmonary disease (COPD) [5].

DM is a cluster of systemic metabolic disorders characterized by persistent hyperglycemia, often resulting in chronic damage to the blood vessels, kidneys, retina, and nervous system [6]. DM presents a significant global health challenge, with the number of individuals affected worldwide reaching 424.9 million in 2017 and projected to increase by 48% to 628.6 million by 2045, resulting in substantial economic and social burdens [7].

DM is classified into type 1 diabetes (T1D), type 2 diabetes (T2D), gestational diabetes, and other specific types of diabetes caused by various etiologies [8]. Among these, T1D and T2D account for over 95% of the total diabetic population [8]. Patients with T2D exhibit varying degrees of relative insulin deficiency, which is often accompanied by insulin resistance. In contrast, T1D arises from autoimmune destruction of pancreatic β cells, resulting in an absolute reliance on exogenous insulin for blood glucose regulation [8]. A meta-analysis of observational studies examining the relationship between DM and IPF revealed that individuals with IPF had a 1.54-fold higher likelihood of developing DM than those without IPF; nevertheless, whether DM increases the risk of IPF remains controversial [9]. However, as DM was not classified in this study, causality between the different subtypes of DM and IPF could not be established.

Mendelian randomization (MR) studies have observational designs that employ randomly-assigned genetic variants as phenotypic instrumental variables (IVs) to establish reliable causal inferences regarding exposures and outcomes. Compared to traditional observational studies, MR is less susceptible to confounding or reverse causation [10].

Therefore, in this study, we aimed to examine the causal association between T1D and IPF, and between T2D and IPF using two-sample univariate and multivariate MR (MVMR) methods.

Methods

Study design

The causal relationships between T1D and IPF, as well as between T2D and IPF in the European population, were initially investigated using a two-sample univariate MR analysis. A sensitivity analysis was conducted to ensure data reliability and validity. Subsequently, MVMR analysis was conducted to examine the significance of the causal relationship between diabetes and IPF established by the univariate MR analysis. We applied a set of selection criteria to identify single nucleotide polymorphisms (SNPs) strongly associated with T1D and T2D, which were subsequently used as IVs in the MR analysis: (I) strong association between IVs and the exposure of interest; (II) no direct relationship between IVs and outcome, with their effect only indicated through exposure [11]; and (III) independence of IVs from any confounding factors affecting both exposure and outcome. The current research design is depicted in Fig. 1. This study adhered to the STROBE-MR reporting guidelines (Supplementary Table 4) [12].

Fig. 1
figure 1

Study design. (A) Univariable MR to identify the causal association between T1D or T2D and IPF, which involves IVs and requires three assumptions: (I) strong association between IVs and the exposure of interest; (II) no direct relationship between IVs and outcome, with their effect only indicated through exposure; (III) independence of IVs from any confounding factors affecting both exposure and outcome. (B) Multivariable MR to evaluate if T1D remains causal given T2D. MR, Mendelian randomization; T1D, type 1 diabetes; IPF, idiopathic pulmonary fibrosis; T2D, type 2 diabetes; IVs, instrumental variables

Full size image

Genome-wide association study (GWAS) data source and instrument selection

We utilized GWAS datasets from various databases for exposure and outcome to mitigate potential bias in causal effect estimation due to sample overlap. The IVs were chosen based on SNPs that were associated with the exposure, as determined by our selection criteria. The T1D and T2D datasets were obtained from the European Bioinformatics Institute. The T1D dataset (ID number: “ebi-a-GCST005536”) comprises 29,652 samples and 101,101 SNPs (https://gwas.mrcieu.ac.uk/datasets/ ebi-a-GCST005536/) [13]. The T2D dataset (ID number: “ebi-a-GCST006867”) comprises 655,666 samples and 5,030,727 SNPs (https://gwas.mrcieu.ac.uk/datasets/ ebi-a-GCST006867/) [14]. The GWAS dataset associated with IPF was established by FinnGen Biobank (ID number: “finn-b-IPF”), and it comprises 198,014 samples and 16,380,413 SNPs (https://gwas.mrcieu.ac.uk/datasets/finn-b-IPF/). The diagnosis of IPF and DM was determined according to the International Classification of Diseases 10th Revision (ICD-10) codes. The detailed characteristics of the exposure and outcome datasets are presented in Table 1.

Table 1 Characteristics of exposures and outcome
Full size table

A genome-wide search was conducted for SNPs (P < 5 × 10− 8) associated with T1D or T2D exposure, which were used as IVs for genetic prediction. Non-independent SNPs based on linkage disequilibrium (r2 < 0.001 within a 10,000 kb aggregation window) were excluded using the European population reference. To ensure independence and exclusivity, these SNPs were searched using Phenoscanner (University of Cambridge, Cambridge, UK; http://www.phenoscanner.medschl.cam.ac.uk/) and any SNPs related to the outcome or confounding factors were excluded. Confounding factors included risk factors other than exposure that may contribute to the outcome. The risk factors associated with IPF include smoking, gastroesophageal reflux, obstructive sleep apnea, herpes virus infection, and certain occupational interstitial lung diseases [15]. To ensure a strong relationship between the IVs and exposure, it is necessary to calculate F-statistics. The F-statistic formula used in this study was F = β2/se2 (β, effect size(exposure); se, standard error (exposure)). An F-statistic greater than 10 indicates a robust association between the IVs and exposure, thereby avoiding weak instrument bias.

Mendelian randomization analysis

We conducted a two-sample univariate MR analysis using T1D and T2D as exposures and IPF as the outcome. In particular, we extracted 36 and 118 strongly-associated SNPs from T1D and T2D GWAS datasets, respectively. Among the 36 SNPs that were strongly associated with T1D in the IPF dataset, rs3184504 was excluded because of its correlation with smoking status (Supplementary Table 1). Among the 118 SNPs strongly associated with T2D in the IPF dataset, four palindromic SNPs were excluded: rs13234269, rs1758632, rs2058913, and rs6494307. Additionally, rs2867125 was excluded due to its association with smoking (Supplementary Table 2).

After retrieving SNPs that were strongly correlated with T1D and T2D in Phenoscanner (University of Cambridge), we identified SNPs that were closely associated with both types of diabetes. To mitigate the potential bias of mutually-confounding factors in the MR analysis, we conducted an MVMR analysis. Fifty-nine SNPs strongly associated with either T1D or T2D were extracted from the IPF dataset, and the palindromic SNP, rs17411031, was excluded during reconciliation. Finally, in the MVMR analysis, 26 and 34 effective IVs were employed to estimate the causal effects of T1D and T2D, respectively, on IPF.

In the univariate MR analysis, we employed inverse variance-weighted (IVW), weighted median (WM), and MR-Egger regression methods to examine the causal association between diabetes and IPF. The IVW method was used as the primary analytical approach. To ensure the reliability of our research findings, we employed MR-IVW and MR-Egger regression methods to assess heterogeneity and quantified the results using Cochran’s Q statistics. We employed the MR-Egger intercept test to examine the presence of horizontal pleiotropy while utilizing a scatter plot of effective IVs for estimation of the causal effect of T1D and T2D on IPF to visually assess whether the outcome effect is zero when the IV effect is zero. A non-zero cutoff indicated the existence of horizontal pleiotropy. In addition, to assess the robustness of the findings, a leave-one-out approach was employed to remove SNPs individually and visualize the results using a leave-one-out plot in the MR analysis. The primary method employed for MVMR analysis was the extended multivariate IVW approach, complemented by secondary methods, such as multivariate MR-Egger and multivariate MR-least absolute shrinkage and selection operator (LASSO). The extended IVW and MR-Egger methods were employed to assess heterogeneity, while the MR-Egger intercept was used to evaluate pleiotropic effects. The statistical analyses in this study were performed using the “TwoSampleMR” and “MendelianRandomization” packages in R v4.3.0 (R Foundation for Statistical Computing, Vienna, Austria).

Ethics

This study relied solely on published GWAS datasets, and all original studies obtained appropriate ethical approval; therefore, separate ethical approval was not required for this study.

Results

Univariable MR estimates

In the two-sample univariate MR analysis, we included 35 and 113 SNPs associated with T1D and T2D, respectively, to estimate their effects on IPF (Supplementary Tables 1 and 2). The F-statistics of the IVs exceeded 10, indicating their robustness and absence of instrumental variable bias [16]. Table 2; Fig. 2 show the results of the analyses of the causal relationship between T1D or T2D and IPF according to the IVW, MR-Egger, and WM methods. The IVW analysis revealed that T1D was associated with an increased risk of IPF (OR = 1.118, 95% CI = 1.021–1.225, P = 0.016), with a mean increase of 11.8% of the likelihood of developing IPF. The WM and MR-Egger methods did not yield statistically significant results (WM: OR = 1.113, 95% CI = 0.983–1.259, P = 0.091; MR-Egger: OR = 1.068, 95% CI = 0.909–1.255, P = 0.427); however, there was a tendency towards a causal effect of T1D on IPF. On the other hand, investigations examining the impact of T2D on IPF consistently indicated an absence of a causal association between T2D and IPF, as evidenced by P-values exceeding 0.05 for IVW, WM, and MR-Egger regression analysis (Table 2; Fig. 2).

Table 2 Two-sample univariable MR analysis of the relationship between IPF and type 1 or 2 diabetes
Full size table
Fig. 2
figure 2

Forest plots of causal effect estimates in univariable Mendelian analysis. T1D, type 1 diabetes; T2D, type 2 diabetes

Full size image

Multivariable MR estimates

We conducted an MVMR analysis to confirm the persistence of a causal relationship between T1D and IPF even after adjusting for T2D. A total of 58 effective IVs were retained for MVMR analysis, of which 26 SNPs were included in the MR analysis of T1D and IPF and 34 SNPs were included in the MR analysis of T2D and IPF (Supplementary Table 3). When T1D was adjusted for T2D, the results of the multivariate analysis were consistent with the results of the univariate MR, suggesting that T1D was significantly associated with an increased risk of IPF (IVW: OR = 1.133, 95% CI = 1.011–1.270, P = 0.032, MR-Lasso: OR = 1.114, 95% CI = 1.004–1.236, P = 0.042, Table 3; Fig. 3). The results of the MR-Egger and median methods were not statistically significant, but tended towards a causal effect consistent with that of the IVW and MR-LASSO analyses, indicating that T1D increases the risk of IPF. Additionally, multivariate analysis of the effect of T2D on IPF yielded results consistent with those of the univariate analysis, indicating no causal relationship. Furthermore, all P-values from IVW, MR-Egger, MR-LASSO, and median analyses were greater than 0.05 (Table 3; Fig. 3).

Table 3 Two-sample multivariable Mendelian randomization analysis of the relationship for T1D or T2D on IPF
Full size table
Fig. 3
figure 3

Forest plots of causal effect estimates in multivariable Mendelian analysis. T1D, type 1 diabetes; T2D, type 2 diabetes; MVMR, multivariate Mendelian randomization; IVW, inverse variance weighted

Full size image

Sensitivity analysis

In the sensitivity analysis, the effective IVs used to estimate the effect of T1D on IPF showed no heterogeneity (IVW: I2 = 0.14, P = 0.232; MR-Egger: I2 = 0.15, P = 0.217, Table 4). The effective IVs used to estimate the effect of T2D on IPF also showed no heterogeneity (I2 = 0, P > 0.05, Table 4). The MR-Egger intercepts for the effects of T1D and T2D on IPF were close to 0, with P-values of 0.464 and 0.777, respectively, suggesting no substantial horizontal pleiotropic bias (Fig. 4). In addition, the leave-one-out analysis showed that no single SNP in T1D or T2D on IPF affected causal effect estimates (Figs. 5 and 6), whereas effect SNPs in the univariate MR analysis were approximately symmetrical in a funnel plot (Fig. 4). The results of the MVMR heterogeneity test showed no significant heterogeneity and no pleiotropic bias (MR-IVW: P = 0.125, MR-Egger: P = 0.109; MR-Egger intercept estimate = 0.003, P = 0.758). These sensitivity analysis results demonstrated the reliability and statistical power of our univariate and MVMR analyses.

Table 4 Sensitivity analyses of MR
Full size table
Fig. 4
figure 4

(A) Scatter plots of causal effect for T1D on IPF, (B) T2D on IPF, (C) Funnel plots to visualize overall heterogeneity of MR estimates for the effect of T1D on IPF, (D) T2D on IPF.T1D, type 1 diabetes; IPF, idiopathic pulmonary fibrosis; T2D, type 2 diabetes; IVW, inverse-variance weighted; MR, Mendelian randomization

Full size image
Fig. 5
figure 5

Forest plot of included SNPs in the MR analysis for causal effect of T1D on IPF. (A) Meta-analysis of SNPs’ effects on IPF based on IVW method with random effect model. (B) Sensitivity analysis by omitting every SNP. T1D, type 1 diabetes; IPF, idiopathic pulmonary fibrosis; IVW, inverse variance weighted

Full size image
Fig. 6
figure 6

Forest plot of included SNPs in the MR analysis for causal effect of T2D on IPF. (A) Meta-analysis of SNPs’ effects on IPF based on IVW method with random effect model. (B) Sensitivity analysis by omitting every SNP. T2D, type 2 diabetes; IPF, idiopathic pulmonary fibrosis; IVW, inverse variance weighted

Full size image

Discussion

In this study, we assessed the causal effects of T1D and T2D on IPF using a two-sample univariate MR approach. We identified a causal relationship between T1D and IPF, but no such association was observed between T2D and IPF. Furthermore, the findings of the MVMR analysis were consistent with those of the univariate analysis: namely, that only T1D was associated with an increased risk of IPF. The polygenic genetic predisposition to T2D is certainly greater than that to T1D; therefore, the current study’s finding of an association between T1D and IPF, rather than T2D, looks even more interesting.

Our univariate and multivariate MR analyses indicated that T1D is associated with an average increase in IPF risk of 11.8% and 13.3%, respectively, and genetically-predicted T1D increases the risk of developing IPF. Previous case-control studies have demonstrated a higher prevalence of diabetes among patients with IPF than among healthy individuals or those with other lung diseases [17,18,19,20,21]; however, these studies did not differentiate between specific types of diabetes or only examined patients with T2D. Owing to potential confounding factors in observational studies, recall bias in retrospective cohort studies, and limited sample sizes in some studies, it is not feasible to establish an accurate causal relationship between diabetes and IPF based solely on these types of observational cohort studies.

A large community-based study conducted in the United States demonstrated that patients with T1D exhibit a higher susceptibility to respiratory disease than those with T2D [22], Therefore, we postulated that T1D has a mechanism that impacts the lungs that is unrelated to hyperglycemia. In a UK population-based case-control study including 920 patients with IPF and 3593 controls, significant associations were observed between IPF and diabetes-related exposures, with the strongest association found with insulin use (OR = 2.36, 95% CI = 1.46–3.83, P < 0.001), and this association was much stronger than that of IPF and oral hypoglycemic drugs [19].

This study revealed the distinct effects of T1D and T2D on IPF. Patients with T2D do not typically require insulin therapy at disease onset or throughout their lifetime, whereas patients with T1D usually require insulin therapy from the outset because of absolute insulin deficiency [8]. Thus, we deduced that the stronger correlation between IPF and insulin use was attributable to a higher proportion of patients with T1D. However, since the diabetes type was not stratified in the original study and the causality direction remained uncertain, we did not make causal inferences. Future studies should conduct more in-depth analyses of the T1D and T2D subtypes to explore their distinct associations with IPF. Differentiating between the subtypes may provide valuable insights into the specific causal mechanisms underlying the relationship between diabetes and IPF.

Our study revealed no causal relationship between T2D and IPF; however, the results of previous observational studies on the association between these two conditions have been controversial. In a retrospective cohort study using the US Centers for Disease Control and Prevention Multi Cause of Death database, patients with diabetes had a lower overall probability of developing IPF compared to healthy controls (OR = 0.81, 95% CI = 0.79–0.82, P < 0.05); moreover, when data were stratified by age, race, and sex, patients with diabetes remained at a reduced risk of developing IPF [23]. This may be attributed to either the occurrence of cardiovascular mortality prior to the development of IPF in patients with diabetes or the prolonged survival of patients with coexisting diabetes and IPF [23]. Another population-based case-control study reported that the incidence of IPF was higher in individuals with T2D compared to non-T2D population (incidence rate ratio = 1.40, 95% CI = 0.76–2.44, P = 0.22), although this difference did not reach statistical significance [24]. Notably, the level of glycemic control was not mentioned among participants with T2D in the original study.

Exposure to genetically-determined T2D does not confer an increased risk of IPF; however, the opposite results have been reported in case-control studies, possibly due to the involvement of persistent hyperglycemia in the development of pulmonary fibrosis caused by unstable glycemic control in certain patients with diabetes. Since our study suggests no causal association between T2D and IPF, future research should focus on investigating the characteristics and pathophysiological processes in patients with T2D to elucidate why they may have a relatively lower risk of developing IPF.

IPF is a progressive interstitial lung disease characterized by autoreactive CD4 T cells, and abnormalities in multiple pathways involved in wound healing and inflammation lead to IPF [25, 26]. T1D is an organ-specific autoimmune disease caused by autoimmune-mediated destruction of insulin-producing β-cells in the pancreas [27], and its pathogenesis is intricate, involving a multitude of factors in its development. Characteristics of T1D such as autoimmune processes, autonomic nerve injury, and systemic vascular injury may increase the risk of respiratory disease [22], but the specific mechanisms by which T1D increases the risk of IPF are not clear. A study that included 120 patients with IPF found that the presence of specific anti-modifier protein antibodies (AMPA) in the serum of a substantial proportion of patients with IPF, suggesting that autoimmunity may characterize a subgroup of IPF and potentially exert an influence on disease progression [28]. KAT2B is a transcriptional co-activator that mediates anti-apoptotic effects under metabolic stress conditions and increases cellular resistance to cytotoxic compounds. High KAT2B expression has been identified as a predictive marker for disease progression to transplantation or mortality in patients with IPF, and was observed to be specifically and progressively upregulated in both individuals who progressed to T1D and those who developed islet cell autoantibodies [29]. Elevated KAT2B expression was observed in peripheral blood mononuclear cells (PBMC) of non-obese diabetic (NOD) mice of T1D, preceding the onset of diabetic hyperglycemia [29]. Collectively, this suggests that KAT2B may play a role in the shared pathogenesis of IPF and T1D during disease development.

Citrullination of proteins may contribute to the shared pathogenesis of both T1D and IPF. Citrullination was identified to modify the bioactivity of glucokinase and inhibit glucose-stimulated insulin secretion in the peripheral blood of T1D patients and NOD mice [30]. Autoreactive T cells and autoantibodies targeting citrullinated beta-cell antigens have been detected in the peripheral circulation and immune infiltrates of individuals with T1D, suggesting that citrullination of β-cell antigens may contribute to the pathogenesis of T1D [31,32,33]. Citrullinated proteins have been detected in bronchoalveolar lavage (BAL) cells from IPF patients, suggesting that the citrullinated pathway is upregulated in IPF [34]. Taken together, post-translational modifications of proteins, particularly citrullination, may contribute to the pathogenesis of IPF and activate immune responses through antibody reactivity [28]. However, limited research has been conducted on the shared pathogenesis of IPF and T1D, and further investigations are required to elucidate their underlying mechanisms.

Hyperglycemia is now known to exert a favorable effect in increasing the risk of IPF, diabetic damage to systemic multi-systems, such as that observed in diabetic nephropathy and diabetic retinopathy, is usually mediated by microvascular injury [35]. The pulmonary interstitium is composed of connective tissue, lymphatic vessels, nerve fibers, and blood vessels, suggesting that diabetes can lead to interstitial lung disease. Both T1D and T2D can affect microvessels through persistent hyperglycemia [35], leading to the development of IPF. However, further research is required to investigate the role of mechanisms other than elevated blood glucose levels in mediating the development of IPF.

In animal models of T1D and T2D, the simulation of clinical hyperglycemia has revealed an increase in DNA double-strand break signals that could not be repaired in a timely manner, resulting in sustained DNA damage [36]. Impaired DNA repair capacity leads to cell senescence, activates the senescence-associated secretory phenotype (including interleukins, inflammatory cytokines, and growth factors), and affects other cell tissues through paracrine signaling, which in turn leads to pulmonary fibrosis [37]. Future research should combine long-term follow-up and clinical case studies of individuals with T1D and IPF to help validate the causal associations and investigate potential biological mechanisms.

Our study presents novel evidence indicating that individuals with a genetic predisposition to T1D were associated with an increased risk of IPF, while no significant association was observed between T2D and IPF. Conflicting results from case-control study results may be due to hyperglycemia resulting from unstable blood glucose control during the development of pulmonary fibrosis. The utilization of diverse hypoglycemic agents could potentially ameliorate pulmonary fibrosis. Recently, several hypoglycemic agents, including metformin, liraglutide, rosiglitazone, and empagliflozin, have been shown to improve pulmonary fibrosis in animal models [38,39,40,41,42].

Currently, apart from metformin, the protective effects of hypoglycemic drugs against pulmonary fibrosis have only been demonstrated in animal and in vitro experiments. Further clinical research is necessary to gain a better understanding of the effects of different types of hypoglycemic drugs on IPF. Additionally, studying the effect of different types of antidiabetic medications on IPF could reveal their potential protective effects and therapeutic strategies.

This study has certain limitations. First, we examined the causal impact of both diabetes types on IPF without exploring the reverse causal effect. When we performed a reverse causality analysis, only six SNPs were screened from the IPF dataset, and these six SNPs did not correspond with the GWAS of T1D and T2D. However, our causal inferences about the two types of diabetes and IPF were reliable; the F-statistics of the IVs we selected were all greater than 10, and the sensitivity analysis did not suggest the existence of significant heterogeneity or pleiotropy. Furthermore, owing to the utilization of GWAS summary-level data instead of individual-level data, we were unable to investigate the impact of variables such as sex, age, and specific types of exposure on outcome or stratify the impact of different glycemic control conditions on IPF. The outcome population in our study consisted of the FinnGen population, although of European origin, represents a distinct and extant group, which may limit the generalization of our conclusions. When analyzing the association between T1D and IPF, researchers should consider exploring factors such as age, sex, specific types of diabetes, diabetes duration, and blood glucose control levels among patients with T1D, and their relationship with IPF development. This would enhance our understanding of the complex relationship between T1D and IPF, and guide clinical management and treatment strategies.

Conclusion

In summary, our MR study demonstrated that T1D significantly elevates the risk of IPF, whereas no significant association was observed between T2D and IPF. Future studies should consider using larger clinical cohorts or conducting clinical trials to validate our findings to ensure more robust and comprehensive conclusions. Moreover, DM should be considered in future studies on the pathogenesis of IPF, particularly focusing on factors other than hyperglycemia in T1D, and attention should be directed towards exploring the potential antifibrotic effects of hypoglycemic drugs. In clinical practice, patients with T1D should undergo lung imaging for early screening for IPF.

Data availability

The GWAS summary datasets for T1D (GWAS ID: ebi-a-GCST005536), T2D (GWAS ID: ebi-a-GCST006867), and IPF (GWAS ID: finn-b-IPF) are available through the ieu open gwas project (https://gwas.mrcieu.ac.uk/datasets).

Abbreviations

IPF:

Idiopathic pulmonary fibrosis

T1D:

Type 1 diabetes

T2D:

Type 2 diabetes

MR:

Mendelian randomization

MVMR:

Multivariate mendelian randomization

GWAS:

Genome-wide association study

SNPs:

Single nucleotide polymorphisms

IVW:

Inverse variance-weighted

WM:

Weighted median

CHF:

Congestive heart failure

DM:

Diabetes mellitus

COPD:

Chronic obstructive pulmonary disease

IVs:

Instrumental variables

LASSO:

Least absolute shrinkage and selection operator

FVC:

Forced vital capacity

References

  1. Raghu G, Remy-Jardin M, Richeldi L, Thomson CC, Inoue Y, Johkoh T, et al. Idiopathic pulmonary fibrosis (an update) and progressive pulmonary fibrosis in adults: an Official ATS/ERS/JRS/ALAT Clinical Practice Guideline. Am J Respir Crit Care Med. 2022;205(9):e18–47.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Maher TM, Bendstrup E, Dron L, Langley J, Smith G, Khalid JM, et al. Global incidence and prevalence of idiopathic pulmonary fibrosis. Respir Res. 2021;22(1):197.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Martinez FJ, Collard HR, Pardo A, Raghu G, Richeldi L, Selman M, et al. Idiopathic pulmonary fibrosis. Nat Reviews Disease Primers. 2017;3:17074.

    Article  PubMed  Google Scholar 

  4. Raghu G, Collard HR, Egan JJ, Martinez FJ, Behr J, Brown KK, et al. An official ATS/ERS/JRS/ALAT statement: idiopathic pulmonary fibrosis: evidence-based guidelines for diagnosis and management. Am J Respir Crit Care Med. 2011;183(6):788–824.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Caminati A, Madotto F, Conti S, Cesana G, Mantovani L, Harari S. The natural history of idiopathic pulmonary fibrosis in a large European population: the role of age, sex and comorbidities. Intern Emerg Med. 2021;16(7):1793–802.

    Article  PubMed  Google Scholar 

  6. ElSayed NA, Aleppo G, Aroda VR, Bannuru RR, Brown FM, Bruemmer D, et al. 2. Classification and diagnosis of diabetes: standards of Care in Diabetes-2023. Diabetes Care. 2023;46(Suppl 1):S19–40.

    Article  CAS  PubMed  Google Scholar 

  7. Standl E, Khunti K, Hansen TB, Schnell O. The global epidemics of diabetes in the 21st century: current situation and perspectives. Eur J Prev Cardiol. 2019;26(2suppl):7–14.

    Article  PubMed  Google Scholar 

  8. 2. Classification and diagnosis of diabetes: standards of Medical Care in Diabetes-2019. Diabetes Care. 2019;42(Suppl 1):S13–28.

    Google Scholar 

  9. Li C, Xiao Y, Hu J, Hu Z, Yan J, Zhou Z, et al. Associations between Diabetes and Idiopathic Pulmonary Fibrosis: a study-level pooled analysis of 26 million people. J Clin Endocrinol Metab. 2021;106(11):3367–80.

    Article  PubMed  Google Scholar 

  10. Davies NM, Holmes MV, Davey Smith G. Reading mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ (Clinical Res ed). 2018;362:k601.

    Article  Google Scholar 

  11. Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG. Using published data in mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30(7):543–52.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Skrivankova VW, Richmond RC, Woolf BAR, Davies NM, Swanson SA, VanderWeele TJ, et al. Strengthening the reporting of observational studies in epidemiology using mendelian randomisation (STROBE-MR): explanation and elaboration. BMJ (Clinical Res ed). 2021;375:n2233.

    Google Scholar 

  13. Onengut-Gumuscu S, Chen WM, Burren O, Cooper NJ, Quinlan AR, Mychaleckyj JC, et al. Fine mapping of type 1 diabetes susceptibility loci and evidence for colocalization of causal variants with lymphoid gene enhancers. Nat Genet. 2015;47(4):381–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Xue A, Wu Y, Zhu Z, Zhang F, Kemper KE, Zheng Z, et al. Genome-wide association analyses identify 143 risk variants and putative regulatory mechanisms for type 2 diabetes. Nat Commun. 2018;9(1):2941.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  15. Lederer DJ, Martinez FJ. Idiopathic pulmonary fibrosis. N Engl J Med. 2018;378(19):1811–23.

    Article  CAS  PubMed  Google Scholar 

  16. Pierce BL, Ahsan H, Vanderweele TJ. Power and instrument strength requirements for mendelian randomization studies using multiple genetic variants. Int J Epidemiol. 2011;40(3):740–52.

    Article  PubMed  Google Scholar 

  17. Enomoto T, Usuki J, Azuma A, Nakagawa T, Kudoh S. Diabetes mellitus may increase risk for idiopathic pulmonary fibrosis. Chest. 2003;123(6):2007–11.

    Article  PubMed  Google Scholar 

  18. García-Sancho Figueroa MC, Carrillo G, Pérez-Padilla R, Fernández-Plata MR, Buendía-Roldán I, Vargas MH, et al. Risk factors for idiopathic pulmonary fibrosis in a Mexican population. A case-control study. Respir Med. 2010;104(2):305–9.

    Article  PubMed  Google Scholar 

  19. Gribbin J, Hubbard R, Smith C. Role of diabetes mellitus and gastro-oesophageal reflux in the aetiology of idiopathic pulmonary fibrosis. Respir Med. 2009;103(6):927–31.

    Article  PubMed  Google Scholar 

  20. Collard HR, Chen SY, Yeh WS, Li Q, Lee YC, Wang A, et al. Health care utilization and costs of idiopathic pulmonary fibrosis in U.S. Medicare beneficiaries aged 65 years and older. Annals Am Thorac Soc. 2015;12(7):981–7.

    Article  Google Scholar 

  21. García-Sancho C, Buendía-Roldán I, Fernández-Plata MR, Navarro C, Pérez-Padilla R, Vargas MH, et al. Familial pulmonary fibrosis is the strongest risk factor for idiopathic pulmonary fibrosis. Respir Med. 2011;105(12):1902–7.

    Article  PubMed  Google Scholar 

  22. George C, Ducatman AM, Conway BN. Increased risk of respiratory diseases in adults with type 1 and type 2 diabetes. Diabetes Res Clin Pract. 2018;142:46–55.

    Article  PubMed  Google Scholar 

  23. Jeganathan N, Miot C, Sathananthan A, Sathananthan M. The association of pulmonary fibrosis with diabetes mellitus. ERJ open Res. 2020;6(4).

  24. Navaratnam V, Davis TME, Hubbard R, Davis WA. Incidence and predictors of idiopathic pulmonary fibrosis complicating type 2 diabetes: the Fremantle Diabetes Study Phase I. Intern Med J. 2021;51(2):276–9.

    Article  PubMed  Google Scholar 

  25. Shum AK, Alimohammadi M, Tan CL, Cheng MH, Metzger TC, Law CS, et al. BPIFB1 is a lung-specific autoantigen associated with interstitial lung disease. Sci Transl Med. 2013;5(206):206ra139.

    Article  PubMed  Google Scholar 

  26. Maher TM, Wells AU, Laurent GJ. Idiopathic pulmonary fibrosis: multiple causes and multiple mechanisms? Eur Respir J. 2007;30(5):835–9.

    Article  CAS  PubMed  Google Scholar 

  27. Taplin CE, Barker JM. Autoantibodies in type 1 diabetes. Autoimmunity. 2008;41(1):11–8.

    Article  CAS  PubMed  Google Scholar 

  28. Kalafatis D, Joshua V, Hansson M, Mathsson-Alm L, Hensvold A, Sköld M. Presence of anti-modified protein antibodies in idiopathic pulmonary fibrosis. Respirol (Carlton Vic). 2023;28(10):925–33.

    Article  Google Scholar 

  29. McKinney EF, Lee JC, Jayne DR, Lyons PA, Smith KG. T-cell exhaustion, co-stimulation and clinical outcome in autoimmunity and infection. Nature. 2015;523(7562):612–6.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  30. Yang ML, Horstman S, Gee R, Guyer P, Lam TT, Kanyo J, et al. Citrullination of glucokinase is linked to autoimmune diabetes. Nat Commun. 2022;13(1):1870.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  31. Buitinga M, Callebaut A, Marques Câmara Sodré F, Crèvecoeur I, Blahnik-Fagan G, Yang ML, et al. Inflammation-Induced Citrullinated glucose-regulated protein 78 elicits Immune responses in human type 1 diabetes. Diabetes. 2018;67(11):2337–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. James EA, Pietropaolo M, Mamula MJ. Immune Recognition of β-Cells: neoepitopes as key players in the loss of Tolerance. Diabetes. 2018;67(6):1035–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Babon JA, DeNicola ME, Blodgett DM, Crèvecoeur I, Buttrick TS, Maehr R, et al. Analysis of self-antigen specificity of islet-infiltrating T cells from human donors with type 1 diabetes. Nat Med. 2016;22(12):1482–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Samara KD, Trachalaki A, Tsitoura E, Koutsopoulos AV, Lagoudaki ED, Lasithiotaki I, et al. Upregulation of citrullination pathway: from autoimmune to idiopathic lung fibrosis. Respir Res. 2017;18(1):218.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Tsilibary EC. Microvascular basement membranes in diabetes mellitus. J Pathol. 2003;200(4):537–46.

    Article  CAS  PubMed  Google Scholar 

  36. Kumar V, Agrawal R, Pandey A, Kopf S, Hoeffgen M, Kaymak S, et al. Compromised DNA repair is responsible for diabetes-associated fibrosis. EMBO J. 2020;39(11):e103477.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Kumar V, Nawroth PP. Is the association between diabetes mellitus and pulmonary fibrosis real? Nat Reviews Endocrinol. 2021;17(12):703–4.

    Article  Google Scholar 

  38. Gamad N, Malik S, Suchal K, Vasisht S, Tomar A, Arava S, et al. Metformin alleviates bleomycin-induced pulmonary fibrosis in rats: pharmacological effects and molecular mechanisms. Volume 97. Biomedicine & pharmacotherapy = Biomedecine &; 2018. pp. 1544–53. pharmacotherapie.

  39. Rangarajan S, Bone NB, Zmijewska AA, Jiang S, Park DW, Bernard K, et al. Metformin reverses established lung fibrosis in a bleomycin model. Nat Med. 2018;24(8):1121–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Fandiño J, Toba L, González-Matías LC, Diz-Chaves Y, Mallo F. GLP-1 receptor agonist ameliorates experimental lung fibrosis. Sci Rep. 2020;10(1):18091.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  41. Yu W, Mi L, Long T. Efficacies of rosiglitazone and retinoin on bleomycin-induced pulmonary fibrosis in rats. Experimental Therapeutic Med. 2017;14(1):609–15.

    Article  CAS  Google Scholar 

  42. El-Horany HE, Atef MM, Abdel Ghafar MT, Fouda MH, Nasef NA, Hegab II et al. Empagliflozin ameliorates Bleomycin-Induced Pulmonary fibrosis in rats by modulating Sesn2/AMPK/Nrf2 signaling and targeting ferroptosis and autophagy. Int J Mol Sci. 2023;24(11).

Download references

Acknowledgements

Summary-level data for SNPs associated with T1D and T2D were extracted from the European Bioinformatics Institute. Summary-level data for SNPs associated with IPF were extracted from the FinnGen Biobank. We would like to thank all the participants for sharing the data.

Funding

We did not receive any funding for this study.

Author information

Authors and Affiliations

  1. Department of Pulmonary and Critical Care Medicine, The affiliated hospital of Qingdao University, Qingdao University, Qingdao, China

    Quou Kang, Jing Ren, Jinpeng Cong & Wencheng Yu

  2. Medical Department of Qingdao University, Qingdao, China

    Quou Kang & Jing Ren

Authors
  1. Quou Kang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Jing Ren
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jinpeng Cong
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Wencheng Yu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

QK: Study conception, Data analyses and interpretation, Methodology, Manuscript draft, Review and editing; JR: Study conception, Methodology, Review and editing; JC: Manuscript revision, Review and editing; WY: Manuscript revision, Review and editing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Wencheng Yu.

Ethics declarations

Ethics approval and consent to participate

Not applicable. This study relied solely on published GWAS datasets, and all original studies obtained appropriate ethical approval and consent to participate.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

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

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

Kang, Q., Ren, J., Cong, J. et al. Diabetes mellitus and idiopathic pulmonary fibrosis: a Mendelian randomization study. BMC Pulm Med 24, 142 (2024). https://doi.org/10.1186/s12890-024-02961-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12890-024-02961-7

Keywords

BMC Pulmonary Medicine

ISSN: 1471-2466

Contact us