Skip to main content

Unraveling the causality between chronic obstructive pulmonary disease and its common comorbidities using bidirectional Mendelian randomization

Abstract

Background

Chronic obstructive pulmonary disease (COPD) frequently coexists with various diseases, yet the causal relationship between COPD and these comorbidities remains ambiguous. As a result, the aim of our study is to elucidate the potential causality between COPD and its common comorbidities.

Methods

We employed the Mendelian randomization (MR) method to analyze single nucleotide polymorphism (SNP) data of common comorbidities with COPD from FinnGen and Integrative Epidemiology Unit (IEU) databases. Causality was primarily assessed using the inverse variance weighting (IVW) method. Multivariable Mendelian randomization (MVMR) analysis was also conducted to eliminate the interference of smoking-related phenotypes. Sensitivity analysis was conducted to ensure the reliability of our findings.

Results

Preliminary univariable MR revealed an increased risk of lung squamous cell carcinoma (LUSC) (IVW: OR = 1.757, 95% CI = 1.162–2.657, P = 0.008), chronic kidney disease (CKD) (IVW: OR = 1.193, 95% CI = 1.072–1.326, P < 0.001), chronic periodontitis (IVW: OR = 1.213, 95% CI = 1.038–1.417, P = 0.012), and heart failure (HF) (IVW: OR = 1.127, 95% CI = 1.043–1.218, P = 0.002). Additionally, the reverse MR analysis indicated that genetic susceptibility to HF (IVW: OR = 1.272, 95% CI = 1.084–1.493, P = 0.003), obesity (IVW: OR = 1.128, 95% CI = 1.056–1.205, P < 0.001), depression (IVW: OR = 1.491, 95% CI = 1.257–1.770, P < 0.001), and sleep apnea syndrome (IVW: OR = 1.209, 95% CI = 1.087–1.345, P < 0.001) could raise the risk of COPD. The MVMR analysis showed no causal effect of COPD on susceptibility to chronic periodontitis after adjusting for smoking.

Conclusions

Our study identified that COPD may elevate the risk of LUSC, HF, and CKD. Additionally, our analysis revealed that HF, sleep apnea symptoms, depression, and obesity might also increase the susceptibility to COPD. These findings revealed a potential causal relationship between COPD and several prevalent comorbidities, which may provide new insights for disease early prediction and prevention.

Background

COPD is one of the most prevalent chronic non-communicable diseases worldwide. Despite the availability of well-established treatment options, many of the symptoms of COPD patients remain uncontrolled, and the mortality is still alarmingly high [1]. Meanwhile, with the extension of human life expectancy, the elderly are confronted with a singular disease and a confluence of multiple systemic diseases referred to as multimorbidity. Patients with COPD also often coexist with other intrapulmonary or extrapulmonary clinical disorders [2, 3], potentially attributable to shared risk factors or the influence of systemic inflammation et al. [4, 5], which further increases the physical burden of patients and social medical expenditure [6, 7]. Strikingly, it was estimated that the final cause of death for more than two-thirds of COPD patients is non-respiratory disorders [2].

Due to time and financial resource constraints, the relationship between COPD and its comorbidities has not been adequately explored. Confounding variables such as aging, smoking, and environmental pollution pose challenges in large-scale cohort studies that provide higher levels of evidence for causation. Therefore, a cost-effective research method that minimizes confounding factors should be employed to establish a causal relationship between COPD and its comorbidities. The findings could lead to early prevention and intervention of multimorbidity within the context of COPD, ultimately improving patients’ quality of life, prolonging their lifespan, and reducing the burden on public health.

The MR analysis is an increasingly employed approach for investigating the causal relationship between observational risk factors and outcomes. Genetic variations, randomly assigned to individuals at birth, are utilized as instrumental variables (IVs) to maximally mitigate the impact of confounding factors mentioned above and reverse causation [8]. Previous studies have utilized MR methods to investigate the causal relationship between COPD and a few comorbidities, including osteoporosis [9], obesity [10], gastroesophageal reflux disease (GERD) [11], and iron deficiency anemia [12]. However, the bidirectional MR analysis has not yet been employed to thoroughly investigate the causal relationship between COPD and its remaining common comorbidities. Therefore, a comprehensive bidirectional MR analysis utilizing up-to-date data is warranted to explore the causal association between COPD and its commonly observed comorbidities.

Methods

Study design

This study is reported in accordance with the STROBE-MR guidelines (Additional file 1: Table S1). A schematic overview of the study design is presented in Fig. 1. In brief, we conducted a MR analysis using publicly available summary statistics from FinnGen and IEU datasets. Both exposure and outcome cohorts were limited to individuals of European ancestry to minimize bias arising from population stratification.

Fig. 1
figure 1

Overall design of our study

Data source for exposure and outcome

The GWAS summary data on COPD were obtained from the FinnGen Research public database, comprising 18,266 cases and 311,286 controls. Additionally, the GWAS dataset from FinnGen and IEU (https://gwas.mrcieu.ac.uk/) was utilized to derive information on the 24 common comorbidities associated with COPD. In the FinnGen study, diseases were defined using the International Classification of Disease (ICD) code. Specific information about related diseases is shown in Table 1.

Table 1 Details of studies included in the MR analyses for the association between COPD and its common comorbidities

Data source for smoking-related phenotypes

The variables chosen for MVMR analysis to represent smoking exposure were the lifetime smoking index and age of smoking initiation. Wootton et al. developed the comprehensive lifetime smoking index, which incorporates information on smoking intensity, duration, as well as initiation and cessation patterns [13]. Besides, the GWAS summary data on age of smoking initiation (id: ieu-b-24) were retrieved from IEU database.

Genetic instrumental variables selection

The fundamental requirements for IVs to fulfill the MR assumptions in this research are as follows: 1) The IVs must be associated with the exposure; 2) the IVs must not be associated with any confounders of the exposure-outcome association; and 3) the IVs should not influence the outcome, except possibly via its association with the exposure [14]. To fulfill the assumption 1 of MR, we exclusively included individuals of European ancestry in our study population and, whenever feasible, restricted the P value of the selected SNPs to 5 × 10–8; however, in cases where there were insufficient IVs for inclusion in the analysis, we could lower the threshold to 5 × 10–6 step by step [15,16,17] (Table 1). To ensure each SNP’s independence, we applied a stringent linkage disequilibrium (LD) correlation coefficient threshold of R2 < 0.001 and a clumping window width of 10000 kb [17]. R2 for each instrument variant: R2 = 2 × EAF × (1 − EAF) × β2, where EAF is the effect allele frequency [18] (Additional file 1: Table S2). Furthermore, we excluded SNPs, rs8040868 and rs16969968 that were associated with confounders or outcomes according to the comprehensive retrieved results of the Phenoscanner V2 database and GWAS catalog to fulfill the assumption 2 and 3 [19, 20]. Detailed traits of comprehensive SNPs retrieved in the Phenoscanner V2 database and GWAS catalog was shown in Additional file 1: Table S3. Additionally, to mitigate the impact of weak IVs on the experimental results, SNPs with an F statistic value less than ten were excluded based on the formula F = β2exposure/SE2exposure. To further verify our findings, we re-performed MR analysis using the threshold of R2 = 0.01 to select diverse eligible IVs [21].

Statistical analysis

UVMR and MVMR analysis

We employed UVMR to evaluate the causal relationship between COPD and its common comorbidities. Given the strong correlation between smoking and COPD, we also conducted the MVMR analysis to assess the independent causal effect of COPD on its comorbidities while adjusting for smoking, which itself has a causal impact on these comorbidities.

The primary approach for MR analysis in UVMR was the inverse-variance weighted (IVW) method, which is utilized to integrate the Wald ratio evaluations of each instrumental variable into a meta-analysis, and is equivalent to performing a weighted linear regression of the associations between the instrumental variables [22]. We also employed various other MR models such as MR-Egger regression, weighted median, and MR pleiotropy residual sum and outlier (MR-PRESSO) methods, were used to examine the aforementioned causality. The MR-Egger regression can be used even when all SNPs are invalid. It estimates the causal effect through a weighted linear regression of gene-outcome coefficients on gene-exposure coefficients [23]. The weighted median approach yields consistent effect estimates when at least half of the weighted variance attributable to horizontal pleiotropy is valid [24]. The MR-PRESSO method, which enables the detection and adjustment for horizontal pleiotropy through outlier removal, was employed to assess pleiotropy [25]. In cases where significant pleiotropy was observed, we utilized the MR-PRESSO method to eliminate outlier SNPs and reanalyze the data. Among these methods, the results obtained from IVW are more reliable compared to the other three approaches. For MVMR analysis, we utilized multivariable IVW (MV-IVW) as the main analysis and multivariable MR-Egger methods as the complementary analysis to estimate the effect of confounder factors on the outcome.

Sensitivity analysis

The sensitivity analysis encompassed assessments of heterogeneity, pleiotropy, and leave-one-out sensitivity tests. Heterogeneity was evaluated using the IVW method and MR-Egger regression, with Cochran's Q-test P > 0.05 means the absence of heterogeneity. If there existed heterogeneity, multiplicative random-effect IVW (IVW-MRE) was chosen as the primary MR analysis method [26]. Horizontal pleiotropy was detected through the intercept term of the MR-Egger method [23]. To address horizontal pleiotropy, we employed MR-PRESSO to identify and remove outliers that were corrected for horizontal pleiotropy (P < 0.05 for outlier detection), subsequently evaluating differences in estimates before and after outlier correction[25]. Additionally, a leave-one-out sensitivity test was conducted to evaluate the robustness of our MR findings [27]. For MVMR, the MVMR-Egger intercept test was performed to examine the horizontal pleiotropy [23]. Given that the exposure and outcome samples for the MR analysis were almost derived from the same database, it is inevitable to encounter sample overlap. To address this issue, we assessed the extent of sample overlap between exposure and outcome in the Finnish database and subsequently performed online calculations via a dedicated platform to determine Type I errors and potential bias (https://sb452.shinyapps.io/overlap/) [28].

The “TwosampleMR” package (version 0.5.7) of R software (version 4.2.2) was employed to analyze the causal-effect relationship between COPD and its common comorbidities. A conservative Bonferroni-corrected threshold (P < 0.05/24, 0.0021, because 24 comorbidities were evaluated for bidirectional analyses) was adopted to address multiple testing. Associations with P < 0.0021 were considered as significant evidence, while associations with P > 0.0021 and P < 0.05 were defined as suggestive evidence based on IVW method.

Results

Characteristics of included IVs and overall experimental design

Through the IVs screening process, the number of instrumental variables associated with 24 common comorbidities of COPD ranged from 3 to 187. The specific details regarding these instrumental variables can be found in Table 1 and Additional file 1: Table S3. However, given that the majority of samples were sourced from the FinnGen database and there was some degree of overlap between them, we searched the sample overlap rate in the FinnGen database and calculated the corresponding magnitude of biases and Type I error rate using available online tools (Additional file 1: Table S4). The potential bias caused by sample overlap was considered insignificant (bias estimate < 0.005). All SNPs involved exhibited F statistics greater than 20, indicating robust genetic instruments. Figure 1 provides an overview of the experimental design employed in this study.

Forward univariable MR analysis

For the preliminary forward MR analysis, it was discovered that COPD is potentially causally related to an increased risk of LUSC (IVW: OR = 1.757, 95% CI = 1.162–2.657, P = 0.008), HF (IVW: OR = 1.127, 95% CI = 1.043–1.218, P = 0.002), osteoporosis (IVW: OR = 1.197, 95% CI = 1.037–1.381, P = 0.014), chronic periodontitis (IVW: OR = 1.213, 95% CI = 1.038–1.417, P = 0.012) and CKD (IVW: OR = 1.193, 95% CI = 1.072–1.326, P < 0.001) (Fig. 2, Table 2 and Additional file 1: Table S5).

Fig. 2
figure 2

Preliminary UVMR assessments of the associations between genetic susceptibility to COPD and its prevalent comorbidities. Purple boxes denote positive correlations, while light blue boxes signify negative associations. The asterisk “*” represents P-values of MR estimates between the Bonferroni-adjusted threshold (P < 0.0021) and 0.05. The double asterisk “**” designates P < 0.0021

Table 2 The preliminary forward UVMR results of significant relationship between COPD and its common comorbidities mainly based on IVW method (P < 0.05)

Upon conducting a sensitivity analysis, heterogeneity was observed between studies of LUSC, osteoporosis and COPD as determined by Cochran’s Q statistics. Additionally, the MR-Egger intercept test revealed the presence of horizontal pleiotropy in the results of the analysis between COPD and osteoporosis. (Additional file 1: Table S6). For the leave-one-out analysis, our findings indicate that no individual genetic variant holds significant influence over the results (Additional file 1: Figure S1). Ultimately, due to the horizontal pleiotropy on the outcome involved between COPD and osteoporosis which cannot be adjusted by MR-PRSSO method, thus the preliminary result that liability to COPD could increase the risk of osteoporosis is unreliable.

Reverse univariable MR analysis

Preliminary reverse MR Analysis results suggest that genetic predisposition to specific comorbidities of COPD may also be associated with its outcome. From one perspective, genetic predisposition to IHD (IVW: OR = 0.867, 95% CI = 0.800–0.939, P < 0.001) may confer a protective effect on the development of COPD. Conversely, genetic susceptibility to HF (IVW: OR = 1.272, 95% CI = 1.084–1.493, P = 0.003), hypertension (IVW: OR = 1.067, 95% CI = 1.016–1.120, P = 0.010), obesity (IVW: OR = 1.128, 95% CI = 1.056–1.205, P < 0.001), depression (IVW: OR = 1.491, 95% CI = 1.257–1.770, P = 0.0005), anemias (IVW: OR = 1.227, 95% CI = 1.061–1.420, P = 0.006) and sleep apnea syndrome (IVW: OR = 1.209, 95% CI = 1.087–1.345, P < 0.001) may exert a detrimental effect on COPD (Fig. 3, Table 3 and Additional file 1: Table S7).

Fig. 3
figure 3

The preliminary MR assessments of the associations between genetic susceptibility to prevalent comorbidities and chronic obstructive pulmonary disease (COPD)

Table 3 The preliminary reverse UVMR results of significant relationship between COPD and its common comorbidities mainly based on IVW method (P < 0.05)

The sensitivity analysis revealed the presence of heterogeneity among the studies on IHD, hypertension, coronary atherosclerosis, anemias, depression, and COPD based on Cochran’s Q statistic. Furthermore, the MR-Egger intercept revealed the presence of horizontal pleiotropy in studies investigating the association among IHD, hypertension, coronary atherosclerosis, and COPD (Additional file 1: Table S8). No SNP was found to drive the above associations (Additional file 1: Figure S2). Finally, the outcomes that the association between IHD, hypertension, coronary atherosclerosis, anemias and COPD were unsteady due to the unadjustable horizontal pleiotropy.

Multivariable MR analysis

Considering the close relationship between smoking and COPD, we employed the multivariable MR to adjust the effect of smoking confounding factors on the outcome. The results showed that after adjusting the influence of lifetime smoking index and age of initiation smoking confounders through MVMR, there was no causal relationship between COPD and chronic periodontitis. However, individuals who are liable to COPD could be prone to suffer LUSC (OR = 1.998, 95%CI = 1.249–3.197, P = 0.004), HF (OR = 1.174, 95%CI = 1.018–1.354, P = 0.027) and CKD (OR = 1.293, 95%CI = 1.070–1.564, P = 0.008) was consistent with the preliminary results based on MVMR-IVW. No significant horizontal pleiotropy was detected in the sensitivity analysis (Additional file 1: Table S9).

The outcome of secondary analysis

The threshold of R2 for the selection of IVs was adjusted to 0.01 to investigate the potential impact of different IVs on the final outcome. Notably, our secondary MR analysis yielded consistent results with those obtained from the preliminary analysis (Additional file 1: Table S10–S14).

Summary results of bidirectional MR analysis

Collectively, after preliminary MR and subsequent sensitivity analysis, our study identified that liability to COPD may elevate the risk of LUSC, HF, and CKD. Additionally, our analysis revealed that genetic predisposition to HF, sleep apnea symptoms, depression, and obesity might also increase the susceptibility to COPD. Among them, there exists strong evidence linking CKD, sleep apnea syndrome, depression, obesity, and COPD. The comprehensive results of our research are depicted in Fig. 4.

Fig. 4
figure 4

The overall outcome of our study

Discussion

To the best of our knowledge, the current study represented the first comprehensive MR analysis to investigate the causal relationship between COPD and its common comorbidities using the UVMR and MVMR approaches. Our findings revealed that genetic liability to COPD caused a higher risk of LUSC, HF and CKD. Additionally, genetic predisposition towards HF, obesity, sleep apnea syndrome, and depression increases the risk of COPD after excluding studies of existing horizontal pleiotropy, which could influence the final accuracy of our outcome.

Lung cancer, known for its rapid progression and poor prognosis, is a challenging disease to manage. Smoking is a shared risk factor for both COPD and lung cancer, interfering with investigating their causal relationship. However, Huang et al. [29] examined 24 pooled case–control studies from the International Lung Cancer Consortium (ILCCO). They demonstrated that 86% of the increased risk of SCLC occurs in individuals with COPD independent of smoking status. Besides, Wang et al. [30] reported that individuals diagnosed with COPD, specifically those with an emphysema-predominant phenotype, exhibited a heightened risk profile for SCLC. Several different hypotheses have been proposed to explain the possibility that COPD increases the risk of lung cancer. One possible explanation is that chronic inflammation and the accumulating release of free radicals may facilitate genetic malignant transformation during tissue reparation [31, 32]. Additionally, the accumulation of detrimental substances resulting from impaired pulmonary ventilation function and airway cilia destruction in patients with COPD could also contribute to the development of lung cancer [33]. Our comprehensive MR analysis indicated that individuals liable to COPD could have a higher chance of developing LUSC, but is not associated with lung adenocarcinoma and SCLC development, indicating that COPD may not increase the risk of all types of lung cancer. Therefore, additional experiments are necessary to elucidate the relationship between COPD and lung malignancy.

Our preliminary UVMR analysis showed that genetically predisposition to COPD could elevate the risk of periodontitis, but this finding was not statistically significant after correcting for smoking confounding factors by MVMR analysis, suggesting that smoking exposure is primarily a contributing factor to periodontitis rather than COPD itself. Although a comprehensive longitudinal cohort study conducted by Shen et al. [34] revealed that individuals with COPD had a 1.19 times greater likelihood of developing periodontitis compared to those without COPD, the study did not consider the additional effect of smoking on periodontitis.

Our research provided significant evidence that individuals with a genetic predisposition to COPD may be more susceptible to developing CKD. A comprehensive case-cohort study revealed that the likelihood of CKD in COPD patients was 1.61 times greater than in those without COPD [35]. Our study results present a contrasting view to those obtained from a previous MR Analysis, where reduced renal function was found to be associated with lower lung function. At the same time, the reciprocal was not statistically significant. This inconsistency could be attributed to the fact that the MR Analysis only examined indicators of pulmonary, forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) and renal function, estimated glomerular filtration rate (eGFR) without considering the disease status [36]. The pathophysiological mechanisms underlying the development of CKD in patients with COPD are multifactorial and complex. However, previous studies have shown that the hypoxia and hypercapnia resulting from COPD can lead to renal vasoconstriction, thus compromising renal blood flow and sympathetic activation [37, 38]. Furthermore, the systematic inflammatory response triggered by COPD may also contribute to the development and progression of CKD.

It is worth noting that sleep apnea is a common condition that often coexists with COPD. The outcome of our reverse MR analysis provided significant evidence supporting the findings that individuals with a genetic predisposition for sleep apnea syndrome may be at an elevated risk of developing COPD. However, there have been varying findings between prior studies. Greenberg-Dotan et al. [39] discovered that sleep apnea can elevate the risk of COPD in a case–control study. In contrast, other community-based controlled studies of older adults have reported the opposite outcome [40]. The differences in the conclusions of the two studies may be due to the differences in the age and number of subjects enrolled in the study. Nonetheless, it has been well established that intermittent hypoxia-induced by sleep apnea syndrome can initiate the activation of hypoxia-driven and oxidative stress pathways, thereby triggering airway and systemic inflammatory responses [41]. To corroborate the results of our MR analysis, Additional file 1: studies with longer follow-up durations may be necessary.

Interestingly, our MR study revealed that genetic liability to depression could elevate the risk of developing COPD. A previous meta-analysis and systematic review, including sixteen relevant follow-up studies conducted by Atlantis et al. [42], demonstrated a link between depression and a higher risk of COPD outcomes. However, Martuccci [43] conducted a phenome-wide association study, which revealed no significant causal relationship between major depressive disorder and lung function at the genetic level, even after excluding the confounding factor of smoking. Nevertheless, it is essential to note that their study solely relied on a genome-wide association study (GWAS) of lung function, and their selection criteria exclusively encompassed patients with depression exhibiting severe symptoms. These methodological disparities may potentially contribute to discrepancies compared to our findings. The underlying mechanisms by which depression may contribute to the development of COPD remain unclear. However, several mechanisms have been speculated, such as depression leading to sympathetic activation and an increase in systemic inflammatory factors [27, 44]. Further studies are needed to establish the precise nature of the association between depression and COPD and the mechanisms that underlie this relationship.

Previous studies have established a significant relationship between COPD and obesity or Body Mass Index. Our latest analysis, which harnessed the power of the FinnGen database, corroborated these findings and suggested that obesity increases the susceptibility to COPD. The observed correlation may be attributed to the impact of obesity on chest volume, which can lead to physical repercussions. The inflammatory factors and hormonal changes associated with obesity are believed to exacerbate the progression of COPD [45, 46].

Our research has shown that liability to COPD are more likely to increase the risk of HF and vice versa. According to a prospective, multicenter, longitudinal cohort study conducted by LEKARZ et al. COPD was associated with a higher myocardial fibrosis burden and HF hospitalization compared to non-COPD patients [47]. Additionally, a recent experimental study has shown that elastase-induced COPD resulted in the development of diastolic cardiomyopathy in an animal model, independent of the confounding effects of cigarette smoke [48]. Given that HF is predominantly prevalent among elderly patients and the progression of COPD is a prolonged process, conducting a prospective cohort study to determine the causal relationship between HF and COPD would prove to be a formidable challenge. Notably, HF exists dysregulation of innate immunity and chronic inflammatory response. Furthermore, HF and COPD share common molecular and mechanistic pathways [49], which suggests a potential explanation for the increased vulnerability of HF patients to COPD development.

Our research has leveraged GWAS data as new as possible to scrutinize the causal connection between COPD and its common comorbidities. This approach has successfully addressed the challenge of establishing causality between diseases due to the impracticability of conducting a large-scale prospective study. Moreover, we have utilized various statistical techniques to ensure the precision of our outcomes and conducted a parallel sensitivity analysis to confirm the dependability of our findings.

Our study is subject to several limitations that warrant discussion. Firstly, to ensure homogeneity by race, we limited our study to populations of European ethnicity. Consequently, the generalizability of our findings is restricted to European populations, indicating the need for future studies that explore more diverse populations. Secondly, it is worth mentioning that majority of the study populations we analyzed were mainly derived from Finnish databases, which could potentially lead to some bias due to the inevitable overlap between exposure and outcome. In the future, the MR analysis based on GWAS databases from diverse sources is needed to verify our conclusions. Last but not least, although the MR analysis serves as an approach to infer the causal relationship between exposure and outcome using IVs, it is still crucial to confirm our findings with large-scale clinical studies or experiments.

Conclusions

To sum up, our study has effectively established a causal connection between COPD and prevalent comorbidities using MR. This discovery could offer fresh insights into preventing and managing concerned diseases in clinical settings. Nonetheless, extensive prospective studies and animal experiments should be conducted to verify the precision and applicability of our results in the future.

Availability of data and materials

The data used in this study were obtained from genome-wide association study summary statistics that were publicly released by genetic consortia. Interested parties may obtain the data by submitting a request to the corresponding author. Furthermore, all datasets generated for this study have been included in the article or additional files.

Abbreviations

MR:

Mendelian randomization

COPD:

Chronic obstructive pulmonary disease

HF:

Heart failure

LUSC:

Lung squamous cell carcinoma

CKD:

Chronic kidney disease

SCLC:

Small cell lung cancer

IEU:

Integrative Epidemiology Unit

SNP:

Single nucleotide polymorphism

IVW:

Inverse variance weighting

LD:

Linkage disequilibrium

OR:

Odds ratio

MR-PRESSO:

MR pleiotropy residual sum and outlier

GWAS:

Genome-wide association study

GERD:

Gastroesophageal reflux disease

FEV1:

Forced expiratory volume in 1 s

FVC:

Forced vital capacity

eGFR:

Estimated glomerular filtration rate

References

  1. Chen S, Kuhn M, Prettner K, Yu F, Yang T, Bärnighausen T, et al. The global economic burden of chronic obstructive pulmonary disease for 204 countries and territories in 2020–50: a health-augmented macroeconomic modelling study. Lancet Glob Health. 2023;11(8):e1183–93.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. McGarvey LP, John M, Anderson JA, Zvarich M, Wise RA. Ascertainment of cause-specific mortality in COPD: operations of the TORCH Clinical Endpoint Committee. Thorax. 2007;62(5):411–5.

    Article  PubMed  PubMed Central  Google Scholar 

  3. van Manen JG, Bindels PJ, van der Zee JS, Bottema BJ, Schadé E. Prevalence of comorbidity in patients with a chronic airway obstruction and controls over the age of 40. J Clin Epidemiol. 2001;54(3):287–93.

    Article  PubMed  Google Scholar 

  4. Venkatesan P. GOLD COPD report: 2023 update. Lancet Respir Med. 2023;11(1):18.

    Article  PubMed  Google Scholar 

  5. Miller J, Edwards LD, Agustí A, Bakke P, Calverley PMA, Celli B, et al. Comorbidity, systemic inflammation and outcomes in the ECLIPSE cohort. Respir Med. 2013;107(9):1376–84.

    Article  PubMed  Google Scholar 

  6. Divo M, Cote C, de Torres JP, Casanova C, Marin JM, Pinto-Plata V, et al. Comorbidities and risk of mortality in patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2012;186(2):155–61.

    Article  PubMed  Google Scholar 

  7. Mannino DM, Higuchi K, Yu T-C, Zhou H, Li Y, Tian H, et al. Economic Burden of COPD in the presence of comorbidities. Chest. 2015;148(1):138–50.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Emdin CA, Khera AV, Kathiresan S. Mendelian randomization. JAMA. 2017;318(19):1925–6.

    Article  PubMed  Google Scholar 

  9. Jiang R, Mou S, Luo F, Zhang Z. Causal relationship between chronic obstructive pulmonary disease and BMD at different sites: a bidirectional Mendelian randomization study. Medicine (Baltimore). 2023;102(41): e35495.

    Article  PubMed  CAS  Google Scholar 

  10. Yang W, Yang Y, Guo Y, Guo J, Ma M, Han B. Obesity and risk for respiratory diseases: a Mendelian randomization study. Front Endocrinol (Lausanne). 2023;14:1197730.

    Article  PubMed  Google Scholar 

  11. Zou M, Zhang W, Shen L, Xu Y, Zhu Y. Major depressive disorder plays a vital role in the pathway from gastroesophageal reflux disease to chronic obstructive pulmonary disease: a Mendelian randomization study. Front Genet. 2023;14:1198476.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Wei Y, Sun L, Liu C, Li L. Causal association between iron deficiency anemia and chronic obstructive pulmonary disease: a bidirectional two-sample Mendelian randomization study. Heart Lung. 2023;58:217–22.

    Article  PubMed  Google Scholar 

  13. Wootton RE, Richmond RC, Stuijfzand BG, Lawn RB, Sallis HM, Taylor GMJ, 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.

    Article  PubMed  Google Scholar 

  14. de Leeuw C, Savage J, Bucur IG, Heskes T, Posthuma D. Understanding the assumptions underlying Mendelian randomization. Eur J Hum Genet. 2022;30(6):653–60.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Grover S, Del Greco MF, Stein CM, Ziegler A. Mendelian randomization. Methods Mol Biol. 2017;1666:581–628.

    Article  PubMed  CAS  Google Scholar 

  16. Ruan X, Chen J, Sun Y, Zhang Y, Zhao J, Wang X, et al. Depression and 24 gastrointestinal diseases: a Mendelian randomization study. Transl Psychiatry. 2023;13(1):146.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Wang S, Jiang H, Qi H, Luo D, Qiu T, Hu M. Association between periodontitis and temporomandibular joint disorders. Arthritis Res Ther. 2023;25(1):143.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Meddens SFW, de Vlaming R, Bowers P, Burik CAP, Linnér RK, Lee C, et al. Genomic analysis of diet composition finds novel loci and associations with health and lifestyle. Mol Psychiatry. 2021;26(6):2056–69.

    Article  PubMed  CAS  Google Scholar 

  19. Kamat MA, Blackshaw JA, Young R, Surendran P, Burgess S, Danesh J, et al. PhenoScanner V2: an expanded tool for searching human genotype-phenotype associations. Bioinformatics. 2019;35(22):4851–3.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  20. Chen Y, Lu T, Pettersson-Kymmer U, Stewart ID, Butler-Laporte G, Nakanishi T, et al. Genomic atlas of the plasma metabolome prioritizes metabolites implicated in human diseases. Nat Genet. 2023;55(1):44–53.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Flatby HM, Ravi A, Damås JK, Solligård E, Rogne T. Circulating levels of micronutrients and risk of infections: a Mendelian randomization study. BMC Med. 2023;21(1):84.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37(7):658–65.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14.

    Article  PubMed  PubMed Central  Google Scholar 

  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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Wang X, Wang X, Wang H, Yang M, Dong W, Shao D. Association between psoriasis and lung cancer: two-sample Mendelian randomization analyses. BMC Pulm Med. 2023;23(1):4.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Laurin C, Moullec G, Bacon SL, Lavoie KL. Impact of anxiety and depression on chronic obstructive pulmonary disease exacerbation risk. Am J Respir Crit Care Med. 2012;185(9):918–23.

    Article  PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  29. Huang R, Wei Y, Hung RJ, Liu G, Su L, Zhang R, et al. Associated links among smoking, chronic obstructive pulmonary disease, and small cell lung cancer: a pooled analysis in the international lung cancer consortium. EBioMedicine. 2015;2(11):1677–85.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Wang W, Dou S, Dong W, Xie M, Cui L, Zheng C, et al. Impact of COPD on prognosis of lung cancer: from a perspective on disease heterogeneity. Int J Chron Obstruct Pulmon Dis. 2018;13:3767–76.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Ballaz S, Mulshine JL. The potential contributions of chronic inflammation to lung carcinogenesis. Clin Lung Cancer. 2003;5(1):46–62.

    Article  PubMed  CAS  Google Scholar 

  32. Peebles KA, Lee JM, Mao JT, Hazra S, Reckamp KL, Krysan K, et al. Inflammation and lung carcinogenesis: applying findings in prevention and treatment. Expert Rev Anticancer Ther. 2007;7(10):1405–21.

    Article  PubMed  CAS  Google Scholar 

  33. Zhai R, Yu X, Shafer A, Wain JC, Christiani DC. The impact of coexisting COPD on survival of patients with early-stage non-small cell lung cancer undergoing surgical resection. Chest. 2014;145(2):346–53.

    Article  PubMed  Google Scholar 

  34. Shen TC, Chang PY, Lin CL, Chen CH, Tu CY, Hsia TC, et al. Risk of periodontal diseases in patients with chronic obstructive pulmonary disease: a nationwide population-based cohort study. Medicine (Baltimore). 2015;94(46): e2047.

    Article  PubMed  Google Scholar 

  35. Chen CY, Liao KM. Chronic obstructive pulmonary disease is associated with risk of chronic kidney disease: a nationwide case-cohort study. Sci Rep. 2016;6:25855.

    Article  PubMed  PubMed Central  ADS  CAS  Google Scholar 

  36. Park S, Lee S, Kim Y, Cho S, Kim K, Kim YC, et al. Kidney function and obstructive lung disease: a bidirectional Mendelian randomisation study. Eur Respir J. 2021;58(6):2100848.

    Article  PubMed  CAS  Google Scholar 

  37. Chabot F, Mertes PM, Delorme N, Schrijen FV, Saunier CG, Polu JM. Effect of acute hypercapnia on alpha atrial natriuretic peptide, renin, angiotensin II, aldosterone, and vasopressin plasma levels in patients with COPD. Chest. 1995;107(3):780–6.

    Article  PubMed  CAS  Google Scholar 

  38. Hemlin M, Ljungman S, Carlson J, Maljukanovic S, Mobini R, Bech-Hanssen O, et al. The effects of hypoxia and hypercapnia on renal and heart function, haemodynamics and plasma hormone levels in stable COPD patients. Clin Respir J. 2007;1(2):80–90.

    Article  PubMed  CAS  Google Scholar 

  39. Greenberg-Dotan S, Reuveni H, Tal A, Oksenberg A, Cohen A, Shaya FT, et al. Increased prevalence of obstructive lung disease in patients with obstructive sleep apnea. Sleep Breath. 2014;18(1):69–75.

    Article  PubMed  Google Scholar 

  40. Zhao YY, Blackwell T, Ensrud KE, Stone KL, Omachi TA, Redline S. Sleep apnea and obstructive airway disease in older men: outcomes of sleep disorders in older men study. Sleep. 2016;39(7):1343–51.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Tuleta I, Stöckigt F, Juergens UR, Pizarro C, Schrickel JW, Kristiansen G, et al. Intermittent hypoxia contributes to the lung damage by increased oxidative stress, inflammation, and disbalance in protease/antiprotease system. Lung. 2016;194(6):1015–20.

    Article  PubMed  CAS  Google Scholar 

  42. Atlantis E, Fahey P, Cochrane B, Smith S. Bidirectional associations between clinically relevant depression or anxiety and COPD: a systematic review and meta-analysis. Chest. 2013;144(3):766–77.

    Article  PubMed  Google Scholar 

  43. Martucci VL, Richmond B, Davis LK, Blackwell TS, Cox NJ, Samuels D, et al. Fate or coincidence: do COPD and major depression share genetic risk factors? Hum Mol Genet. 2021;30(7):619–28.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Bremmer MA, Beekman ATF, Deeg DJH, Penninx BWJH, Dik MG, Hack CE, et al. Inflammatory markers in late-life depression: Results from a population-based study. J Affect Disord. 2008;106(3):249–55.

    Article  PubMed  CAS  Google Scholar 

  45. Cox AJ, West NP, Cripps AW. Obesity, inflammation, and the gut microbiota. Lancet Diabetes Endocrinol. 2015;3(3):207–15.

    Article  PubMed  CAS  Google Scholar 

  46. Kopelman PG. Hormones and obesity. Baillieres Clin Endocrinol Metab. 1994;8(3):549–75.

    Article  PubMed  CAS  Google Scholar 

  47. Lagan J, Schelbert EB, Naish JH, Vestbo J, Fortune C, Bradley J, et al. Mechanisms underlying the association of chronic obstructive pulmonary disease with heart failure. JACC Cardiovasc Imaging. 2021;14(10):1963–73.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Grillet PE, Desplanche E, Wynands Q, Gouzi F, Bideaux P, Fort A, et al. Diastolic cardiomyopathy secondary to experimentally induced exacerbated emphysema. Am J Respir Cell Mol Biol. 2023;69(2):230–41.

    Article  PubMed  CAS  Google Scholar 

  49. Rahman HH, Rashid MH, Miah NA, Israt S, Atiqullah S, Akbar MS. Correlation study between COPD and heart failure in elderly patient. Mymensingh Med J. 2022;31(2):498–505.

    PubMed  CAS  Google Scholar 

Download references

Acknowledgements

The authors wish to express their gratitude to the participants and researchers of the FinnGen and IEU studies. We sincerely thank for Biorender, a figure drawing tool, to be used to draw Figs. 1 and 4.

Funding

YC.S. is supported by the National Natural Science Foundation of China (81970041, 82170048).

Author information

Authors and Affiliations

Authors

Contributions

All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.

Corresponding author

Correspondence to Yongchang Sun.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

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:

Figure S1. Plots for "leave-one-out" analysis for the causal impact of COPD on potentially causal comorbidities. Figure S2. Plots for "leave-one-out" analysis for the causal impact ofpotentially causal comorbidities on COPD.

Additional file 2:

Table S1. STROBE-MR checklist of Mendelian randomization study. Table S2. The comprehensive retrieved results of traits corresponding to SNPs in Penoscanncer V2 and GWAS catalog. Table S3. Comprehensive details of incorporated SNPs for each disease investigated (R2<0.001). Table S4. Sample overlap rate between exposure and outcome and potential bias. Table S5. The UVMR results of COPD to its frequently coexisting comorbidities (R2<0.001). Table S6. Pleiotropy and heterogeneity analysis of concerned forward UVMR analysis (R2<0.001). Table S7. The reverse UVMR results of COPD to its frequently coexisting comorbidities (R2<0.001). Table S8. Pleiotropy and heterogeneity analysis of concerned reverse UVMR analysis (R2<0.001). Table S9. The MVMR results with adjusting the smoking exposure (R2<0.001). Table S10. Comprehensive details of incorporated SNPs for each disease investigated (R2<0.01). Table S11. The UVMR results of COPD to its frequently coexisting comorbidities (R2<0.01). Table S12. Pleiotropy and heterogeneity analysis of concerned forward UVMR analysis (R2<0.01). Table S13. The reverse UVMR results of COPD to its frequently coexisting comorbidities (R2<0.01). Table S14. Pleiotropy and heterogeneity analysis of concerned reverse UVMR analysis (R2<0.01).

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

Wang, Z., Sun, Y. Unraveling the causality between chronic obstructive pulmonary disease and its common comorbidities using bidirectional Mendelian randomization. Eur J Med Res 29, 143 (2024). https://doi.org/10.1186/s40001-024-01686-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40001-024-01686-x

Keywords