- Research
- Open access
- Published:
Assessing the genetic relationship between gastroesophageal reflux disease and chronic respiratory diseases: a mendelian randomization study
BMC Pulmonary Medicine volume 23, Article number: 243 (2023)
Abstract
Background
Previous observational studies have found an association between gastroesophageal reflux disease (GERD) and chronic respiratory diseases, but it remains uncertain whether GERD causally influences these diseases. In this study, we aimed to estimate the causal associations between GERD and 5 chronic respiratory diseases.
Methods
88 GERD-associated single nucleotide polymorphisms (SNPs) identified by the latest genome-wide association study were included as instrumental variables. Individual-level genetic summary data of participants were obtained from corresponding studies and the FinnGen consortium. We applied the inverse-variance weighted method to estimate the causality between genetically predicted GERD and 5 chronic respiratory diseases. Furthermore, the associations between GERD and common risk factors were investigated, and mediation analyses were conducted using multivariable MR. Various sensitivity analyses were also performed to verify the robustness of the findings.
Results
Our study demonstrated that genetically predicted GERD was causally associated with an increased risk of asthma (OR 1.39, 95%CI 1.25–1.56, P < 0.001), idiopathic pulmonary fibrosis (IPF) (OR 1.43, 95%CI 1.05–1.95, P = 0.022), chronic obstructive disease (COPD) (OR 1.64, 95%CI 1.41–1.93, P < 0.001), chronic bronchitis (OR 1.77, 95%CI 1.15–2.74, P = 0.009), while no correlation was observed for bronchiectasis (OR 0.93, 95%CI 0.68–1.27, P = 0.645). Additionally, GERD was associated with 12 common risk factors for chronic respiratory diseases. Nevertheless, no significant mediators were discovered.
Conclusions
Our study suggested that GERD was a causal factor in the development of asthma, IPF, COPD and chronic bronchitis, indicating that GERD-associated micro-aspiration of gastric contents process might play a role in the development of pulmonary fibrosis in these diseases.
Background
Gastroesophageal reflux disease (GERD) is a multifaceted disorder with multiple symptoms, such as heartburn, regurgitation, chest or epigastric pain, nausea, bloating, and so on, which occurs when the contents of the stomach reflux into the esophagus or laryngopharynx [1]. The prevalence of GERD is around 15–20% in western countries and has increased from 18.1 to 27.8% in North America, from 8.8 to 25.9% in Europe, and from 2.5 to 7.8% in East Asia 2005 to 2011 [2, 3].
Chronic respiratory diseases, such as idiopathic pulmonary fibrosis (IPF), chronic obstructive disease (COPD), and asthma, share common features of airway inflammation that can be exacerbated by prolonged exposure to gastric fluid, indicating that GERD might contribute to their development [4]. Epidemiological evidence has been accumulating over the past decade linking GERD with an increased risk of chronic respiratory diseases. Several studies have reported that GERD was a significant risk factor for asthma [5, 6], COPD [7], and IPF [8, 9]. For example, compared with population controls, the odds of having GERD were higher in IPF cases [8]. A recent meta-analysis also found an association between GERD and increased risk of IPF, although the relationship was not significant after adjusting for smoking [9].
However, the results of observational studies may be influenced by potential confounders and reverse causality, leading to inconclusive conclusions. For instance, few studies have controlled for body mass index (BMI), which is known to be correlated with the risk of both GERD and chronic respiratory diseases. Additionally, other confounding factors such as smoking status and alcohol consumption may also interfere with the results [10, 11] .Therefore, the existence of a true causal relationship between GERD and chronic respiratory diseases remains uncertain.
Mendelian randomization (MR) analysis is an increasingly popular epidemiological approach for assessing the causation between an exposure and an outcome. By using genetic variants as instrumental variables (IVs), MR analysis is less susceptible to unobserved confounders and reverse causation [12]. Since genetic variations are randomly distributed at conception according to Mendel’s second law, they are generally independent of environmental risk factors and precede the development of diseases [13]. Moreover, with the availability of summary data from genome-wide association studies (GWASs) databases, it is possible to explore the potential relationship and causality between GERD and chronic respiratory diseases. In this study, we utilized published GWASs summary data and conducted MR analyses to comprehensively evaluate the causality [14, 15].
Materials and methods
Two-sample MR analyses were performed to investigate the causal relationship between GERD and the risk of 5 chronic respiratory diseases, as well as 12 related risk factors. In addition, multivariable MR analyses were performed to assess the direct causal impacts of GERD on chronic respiratory diseases by adjusting for each potential mediator, which was statistically significant. Our complete study design overview is depicted in Fig. 1.
Genetic instruments data for GERD
We obtained IVs for GERD from a meta-analysis of GWASs in European populations, which combined data from the UK Biobank study and the QSKIN study, including a total of 78,707 cases and 288,734 controls. In the meta-analysis, 88 independent lead single-nucleotide polymorphisms (SNPs) that surpassed the genome-wide significance threshold of p < 5 × 10− 8 were identified [16]. These IVs collectively explained 3.9% of the variance in GERD [17]. Subsequently, 14 out of the 88 SNPs were removed after conducting linkage disequilibrium (LD) analysis with a threshold of R2 < 0.001. Additionally, after harmonizing the genetic data of GERD and chronic respiratory diseases, three SNPs (rs2145318, rs2358016, rs957345) were excluded from the IV set, as they were prone to bias in two-sample MR analysis due to their palindromic nature and intermediate allele frequencies. Ultimately, a total of 71 SNPs were included in the final IV set as genetic variables. Detailed information on the IVs for GERD can be found in Supplementary Data Table S1.
GWASs summary data on 5 chronic respiratory diseases
Our research aimed to determine whether there is a causal link between GERD and several chronic respiratory diseases that commonly have high morbidity rates and wide-ranging impacts on both males and females. To achieve this, we included five chronic respiratory diseases as outcomes in our primary analyses, namely asthma, COPD, IPF, bronchiectasis, and chronic bronchitis. To ensure the accuracy of our results, we applied strict inclusion criteria, including using publicly available summary data on chronic respiratory diseases released within the past five years and studying patients of European ancestry. Additionally, the diagnostic approach of included patients based not only on ICD-coding, but also other codes such as examinations or drugs. The genetic participants’ sample for asthma was obtained from a large GWAS summary data, which included 56,167 cases and 352,255 controls [18]. Asthma cases from UK Biobank accepted a broad definition based on a diagnosis from UK Biobank Outcome Adjudication Group, hospital, primary care, and self-reported related records (data-field 20,002 in UK Biobank). Additionally, we included the summary data of IPF from the largest available GWASs to date, which was a meta-analysis including 4125 IPF cases and 20,464 controls [19]. IPF cases were identified based on primary or secondary ICD10 code J84.1 (Other interstitial pulmonary diseases with fibrosis). We are grateful to the FinnGen group for their hard work in providing us with summary statistics for COPD (16,410 cases and 283,589 controls), chronic bronchitis (1,035 cases and 283,589 controls), and bronchiectasis (1,967 cases and 283,589 controls) [20]. COPD cases were defined as those with primary or secondary ICD10 code J43 (Emphysema) or J44 (Other chronic obstructive pulmonary disease), and chronic bronchitis cases were defined as those with primary or secondary ICD10 code J42 (Unspecified chronic bronchitis) [20]. Detailed information on studies included in analyses is presented in Table 1. Additionally, we retrieved data from other groups in the UK Biobank and FinnGen group to replicate the casual association between genetically predicted GERD and chronic respiratory diseases (Supplementary Data Table S10 and S11).
GWASs summary data on 12 common risk factors of chronic respiratory diseases
We collected summary data for 12 common risk factors associated with chronic respiratory diseases. These included smoking status (previous and current) from Neale Lab (dataset ID: ukb-a-224 and ukb-a-225) [21], initiating smoking from the GWAS & Sequencing Consortium of Alcohol and Nicotine use (GSCAN; dataset ID: ieu-b-4877) [22], exposure to tobacco smoke at home and outside home from Neale Lab (dataset ID: ukb-a-19 and ukb-a-20), obesity class 1–3 from the Genetic Investigation of ANthropometric Traits (GIANT; dataset ID: ieu-a-90, ieu-a-91, ieu-a-92) [23], alcohol consumption from the UK Biobank (dataset ID: ieu-a-1283) [24], Type 2 diabetes from the open GWAS study (dataset ID: ebi-a-GCST006867) [25], Vitamin D deficiency from the FinnGen group (dataset ID: finn-b-E4_VIT_D_DEF) [20], and lung function (FEV1/FVC) from the open GWAS study (dataset ID: ebi-a-GCST007431) [26]. Detailed information on the studies mentioned above is presented in Supplementary Data Table S2.
Statistical analysis
Two-sample MR analyses between GERD and risk of 5 chronic respiratory diseases outcomes
We first applied a two-sample MR model to access the causal effect of GERD on 5 chronic respiratory diseases. The main MR analyses were performed using the inverse-variance weighted (IVW) method, which involves a meta-analysis of the Wald ratio of each single nucleotide polymorphism (SNP) between the exposure (GERD) and the outcome (chronic respiratory diseases) using a random-effects inverse-variance method. This approach weights each ratio based on its standard error, taking into account possible heterogeneity [27]. In this regard, given that phenotype of GERD (X) on the risk of outcomes (Y) and the genetic variants (G) as IVs, the IVW analyses were performed by combining genetic variants on GERD (bYG) with their standard errors (SEYG) to examine the causal relationship between genetically determined exposure factor (βGERD) and outcome phenotypes. Odds ratios (ORs) with 95% confidence intervals (CIs) were calculated and estimated for the five chronic respiratory diseases. We considered a causal relationship as significant if the Bonferroni-corrected threshold P-value was < 0.01 (0.05/5), indicating a stringent correction for multiple testing.
Two-sample MR analyses the associations of GERD and 12 common risk factors for chronic respiratory diseases
We further conducted Mendelian randomization (MR) analysis to investigate the potential pathogenesis from gastroesophageal reflux disease (GERD) to chronic respiratory diseases by examining whether genetic susceptibility for GERD could be associated with common risk factors of chronic respiratory diseases. Although the underlying pathophysiological processes of many of them has not been fully understood, previous studies have observed associations between potential genetic and environmental factors and exacerbation risk of these diseases [28,29,30,31,32]. In our analysis, we used the inverse-variance weighted (IVW) method to estimate the causal effects of genetically predicted GERD on 12 common risk factors for chronic respiratory diseases. We applied a Bonferroni correction with a significance threshold of P-value < 4.2 × 10− 3 (0.05/12) to account for multiple testing.
Two-sample MR analyses between genetically predicted risk factors and risk of five chronic respiratory diseases outcomes
Furthermore, we selected specific risk factors that showed a significant causal relationship with GERD to estimate their causal impact on individual chronic respiratory diseases, for which GERD was statistically significant. We obtained IVs that were independent genome-wide significant (R2 < 0.001) for smoking initiation, obesity class 1–3, Alcohol consumption, Type 2 diabetes, Vitamin D deficiency, and Lung function (FEV1/FVC) from the summary data of corresponding GWASs [20, 22,23,24,25,26]. IVs for smoking status, exposure to tobacco smoke, and vitamin D deficiency were obtained from the MR-base database (https://gwas.mrcieu.ac.uk/) (accessed on 25 December 2022) [21]. Additionally, we applied a Bonferroni-corrected P-value < 1.3 × 10− 2 (0.05/4; 4 outcomes that GERD was significantly associated with) as the significance threshold.
Assuming that causal relationships were observed in all three previously mentioned analyses, the risk factor being analyzed was considered as a potential pathway mediator linking GERD to this particular chronic respiratory disease outcome.
Multivariable MR analysis and mediation analysis
We further performed multivariable MR analysis to identify potential mediators between GERD and chronic respiratory diseases and investigate their influence on causal estimates. Multivariable MR analysis c allows for a more accurate estimation of the true causal effect by evaluating the effects of multiple genetic variations and multiple confounding factors in the causal relationship [33]. Specifically, multivariable MR analysis can simultaneously consider the effects of multiple sets of genetic variations and use multiple linear regression to control multiple confounding factors to eliminate their effects on the causal relationship. Multivariable MR analysis can also compare the causal effects between different genetic variations and outcome variables to further address potential collinearity and reverse causality problems. Through these methods, multivariable MR analysis can more accurately test causal hypotheses without the need for complex statistical methods and mathematical models [33]. If the analysis reveals that a specific risk factor significantly weakened the causal association between GERD and diseases of outcomes (P < 0.01 changed to P > 0.05), the risk factor could be identified as a mediator [33].
Sensitivity analysis
We took advantage of the online web tool (https://shiny.cnsgenomics.com/mRnd/) (accessed on 30 December 2022) to estimate the F statistics and power for the MR analysis [34]. In addition, we employed several other methods, including the Mendelian Randomization Pleiotropy Residual Sum and Outlier (MR-PRESSO) test, MR-Egger regression, and weighted-median estimator, to comprehensively evaluate our MR results [35]. Among them, the MR-PRESSO global test could detect the horizontal pleiotropy of overall IVs based on a regression framework. Moreover, its outlier test can produce a corrected estimate by performing outlier removal if desired, and the distortion test could compare whether there was a significant difference between the results before and after outlier correction [36]. Furthermore, MR-Egger regression, which assumed independent pleiotropic associations, was conducted to explore pleiotropy and estimate causality adjusted for pleiotropy [37]. We also conducted leave-one-out analyses to assess the potential bias from single sensitive SNPs in the IVW method. Additionally, we performed the MR-heterogeneity test using Cochran’s Q test on the IVW and MR-Egger estimate to identify SNPs contributing to heterogeneity in causal estimation [38].
All of our MR analyses were performed in R (version 4.2.1) using the package TwoSampleMR (version 0.5.6) and MRPRESSO (version 1).
Results
In our study, all 71 SNPs that we selected were chosen based on the genome-wide significance threshold of p < 5 × 10− 8, ensuring that they met the MR assumption. The F-statistics for the instrumental variables (IVs) were all above 10, indicating that we did not use weak genetic instruments in our study. Furthermore, the power was sufficient in all MR analyses, except for bronchiectasis, as shown in Supplementary Data Table S3.
Genetically predicted GERD and 5 chronic respiratory diseases
The result of the IVW method showed a positive association between genetically predicted GERD and asthma (OR 1.39, 95%CI 1.25–1.56, P < 0.001), COPD (OR 1.64, 95%CI 1.41–1.93, P < 0.001) and chronic bronchitis (OR 1.77, 95%CI 1.15–2.74, P = 0.009). Additionally, our study indicated a potential association between genetically predicted GERD and IPF (OR 1.43, 95%CI 1.05–1.95, P = 0.022). However, no significant association was observed between GERD and bronchiectasis (OR 0.93, 95%CI 0.68–1.27, P = 0.645) (Fig. 2). These results were consistent with MR-PRESSO analysis after excluding this outlier. Causal estimates for all five outcomes remained consistent in sensitivity analyses using the weighted median method. Nevertheless, we found that the result of the weighted median method and MR-Egger were not corresponding with them, while they offered less accurate estimates than conventional MR (IVW method). Detailed results of MR using MR-PRESSO, weighted median and MR-Egger method are presented in Table 2 and Supplementary Data Table S4.
Genetically predicted GERD and risk factors of chronic respiratory diseases
GERD were found to be positively associated with seven risk factors, including smoking initiation (OR 1.36, 95%CI 1.25–1.47, P < 0.001), exposure to tobacco smoke at home and outside home (OR 1.03, 95%CI 1.02–1.05, P < 0.001; and OR 1.08, 95%CI 1.06–1.10, P < 0.001 respectively), obesity class 1–3(OR 1.95, 95%CI 1.49–2.56, P < 0.001, OR 2.23, 95%CI 1.54–3.24, P < 0.001, and OR 2.55, 95%CI 1.47–4.44, P < 0.001, respectively), Type 2 diabetes (OR 1.86, 95%CI 1.53–2.26, P < 0.001), and negatively associated alcohol consumption (OR 0.93, 95%CI 0.89–0.96, P < 0.001). Additionally, evidence from the IVW method indicated potential associations between genetically predicted GERD and smoking status (previous and current; OR 1.03, 95%CI 1.01–1.05, P = 0.013, and OR 1.05, 95%CI 1.03–1.06, P = 0.033, respectively). Notably, our analyses showed no statistically significant associations between genetically predicted GERD and vitamin D deficiency (OR 1.12, 95%CI 0.40–3.13, P = 0.830) and lung function (FEV1/FVC; OR 1.02, 95%CI 0.97–1.07, P = 0.385). The MR results using the IVW model are presented in Fig. 3.
Genetically predicted risk factors and chronic respiratory diseases of outcomes
Further evaluation was conducted to determine whether smoking initiation, exposure to tobacco smoke at home and outside the home, obesity class 1–3, Type 2 diabetes, alcohol consumption, and smoking status (previous and current) were causally connected with four chronic respiratory diseases that GERD was positively associated with. The IVW method indicated that obesity class 1–2 (OR 1.06, 95%CI 1.03–1.09, P < 0.001, and OR 1.05, 95%CI 1.02–1.08, P < 0.001 respectively) and smoking initiation (OR 1.16, 95%CI 1.05–1.28, P = 0.003) were associated with an increased risk for asthma (Supplementary Data Table S5). Additionally, positive associations between smoking initiation (OR 1.72, 95%CI 1.48-2.00, P < 0.001), exposure to tobacco smoke outside home (OR 5.92, 95%CI 1.76–19.93, P = 0.004), current smoking status (OR 135.63, 95%CI 5.93-3100.91, P = 0.002) and COPD were found (Supplementary Data Table S5).
Multivariable MR analysis
We conducted further analysis to investigate the influence of potential mediators on the causal effects of GERD on asthma and COPD. The results from multivariable MR analyses showed that the causal estimates of GERD with outcomes remained robust even after adjusting for each potential mediator (Table 3).
Sensitivity analyses
We also performed sensitivity analyses to assess the robustness of our findings. Egger intercept analyses showed little evidence of pleiotropy for asthma, IPF, and COPD, with P-values ranging from > 0.01 to < 0.05 (Table 2 and Supplementary Data Table S6). MR-PRESSO tests identified two outliers for asthma and one outlier for COPD, but no significant statistical evidence of distortion was found through the distortion tests (Table 2 and Supplementary Data Table S7). Overall, we concluded from these analyses that horizontal pleiotropy has little or no interference with our causal estimates. Cochran Q statistics indicated the presence of heterogeneity in the correlation between asthma and COPD (Table 2 and Supplementary Data Table S8). Additionally, no single SNP strongly drove the overall effect of GERD on the five chronic respiratory diseases, as indicated by the leave-one-out sensitivity analysis (Supplementary Figure S1-S5).
Discussion
With the two-sample MR approach, our study illustrated that patients with GERD were causally associated with an increased risk of asthma, IPF, COPD, and chronic bronchitis, suggesting the process of GERD-associated micro-aspiration of gastric contents may play a role in the development of these respiratory diseases. However, no significant association was observed between GERD and bronchiectasis.
Our results were consistent with previous epidemiological evidence that has reported an increased prevalence of GERD in chronic respiratory diseases. A systematic review in 2010 found that asthma patients were more likely to have GERD compared to controls, with an average prevalence of 22.0% in asthma patients versus 4.8% in controls (pooled OR 5.6, 95% confidence interval [CI] 4.3–6.9) [6]. A recent systematic review and meta-analysis by Mallah et al. also demonstrated an association between GERD and asthma exacerbation (OR 1.27, 95% CI 1.18–1.35) in various study designs and populations, including European and non-European populations [5]. Furthermore, a prospective study found that 21% of patients with severe uncontrolled asthma had signs of GERD detected on bronchoscopy, suggesting that GERD may be a trigger for uncontrolled asthma [39]. The association between GERD and IPF has also been reported in previous studies. Tobin et al. demonstrated a high correlation between abnormal esophageal acid exposure and the incidence of IPF in a comparative study conducted 20 years ago [40]. Recent observational studies have also reported similar findings. A retrospective study by Lee et al. with a sample of 786 IPF patients suggested that the prevalence of GERD increased over the follow-up period, with eventually 107 patients (13.6%) being diagnosed with GERD [41]. A population-based, case-control study also found a higher incidence of GERD in IPF cases compared to population controls (OR 1.78, 95% CI 1.09–2.91, P = 0.02) [8]. In COPD, the prevalence of GERD has been reported to be high, and its association with adverse outcomes has been demonstrated in several studies. A Korean national cross-sectional cohort study found a prevalence of GERD of 28% in COPD patients, and the presence of GERD was associated with an increased risk of hospitalization (OR 1.54, 95% CI 1.50–1.58, P < 0.001) and frequent emergency room visits (OR 1.55, 95% CI 1.48–1.62, P < 0.001) [42]. Similar findings were reported by Huang et al. in a systematic review with meta-analysis, which showed an increased risk of COPD exacerbation associated with GERD (OR 5.37, 95% CI 2.71–10.64) [7]. Another study based on the COPDGene cohort reported a prevalence of GERD of 29.1% in COPD patients, who were more likely to have symptoms of chronic bronchitis and worse scores of life quality [43].
Notably, our MR analysis did not identify a causal relationship between bronchiectasis and GERD, which was somewhat different from previous studies. A recent meta study mentioned that the prevalence of GERD in bronchiectasis ranged from 34 to 74% according to self-reported symptoms and questionnaires and pepsin was detected in sputum samples in 26-70% of patients with mild to moderate bronchiectasis [44]. However, the relationship between bronchiectasis and GORD remains controversial actually. For example, a retrospective cohort study including 192 patients suggested that gastroesophageal reflux was not found to be an important risk factor for the development of bronchiectasis [45]. In addition, two case-control studies did not identify any association with decreased lung function or other markers of disease severity of bronchiectasis patients with GERD. Nevertheless, due to recruitment difficulties, the efficacy of these studies to detect such effects was clearly inadequate and a single time dimension may not be sufficient to accurately reflect the relationship between GERD and bronchiectasis [46, 47]. Consequently, more observational studies are warranted to support our results and investigate their potential connections.
Nevertheless, owing to the flaws of observational research design, there may be some heterogeneity in patient characteristics among the primary studies that could bias the results. For example, GERD could be defined from varied perspectives. According to the Montreal definition, GERD is a condition that develops when the reflux of stomach contents causes troublesome symptoms and/or complications [48]. However, some studies diagnosed GERD based solely on pH-metry, others on a symptom inventory, and only a few inferred from diagnostic code entries, which could generate multiple deviations. Therefore, the genuine relationship between GERD and them might also be biased owing to the inaccurate diagnosis of GERD. Additionally, the association estimates from case-control studies, as indicated by Bédard et al., may have little confidence [9]. Unfortunately, there is still no prospective large-scale longitudinal cohort study that has investigated the relationship between GERD and asthma, idiopathic pulmonary fibrosis (IPF), chronic obstructive pulmonary disease (COPD), bronchiectasis, and chronic bronchitis. The results of our MR analysis revealed that genetically predicted GERD might be associated with an increased risk of asthma, IPF, COPD, and chronic bronchitis, supporting a causal role of GERD in the development of them.
Current research findings have identified several risk factors associated with chronic respiratory diseases, such as Smoking, which is one of the most widely known risk factors for respiratory diseases. Additionally, a growing number of studies have suggested that tobacco smoke exposure may also contribute to the development of asthma, COPD, and chronic bronchitis [30,31,32]. A 2016 international report indicated that low FVC could be a potential risk factor for acute exacerbation of IPF [28]. Furthermore, higher body mass index and alcohol might also contribute to the exacerbation of IPF, COPD and chronic bronchitis [28, 31, 32]. Interestingly, our further detailed MR analysis confirmed that GERD may be causally associated with a broad range of them, including smoking, smoking status, exposure to tobacco, obesity, alcohol consumption, and type 2 diabetes, while no statistically significant associations were detected with vitamin D deficiency and lung function. Additionally, although smoking is a well-known risk factor for IPF or chronic bronchitis, it did not show statistical significance in this study, other two MR studies have reported similar results. One of them discussed the relationship between hypothyroidism and IPF, and the authors also did not find any statistically significant evidence for an association between smoking and IPF in the multivariable MR analysis based on BMA [49]. Another article explored the role of gastro-esophageal reflux disease in the development of IPF also discovered the same results in the MR analysis using identical data source [50]. Consequently, while smoking may be a risk factor for the development of IPF or chronic bronchitis from a macro epidemiological perspective, it may not be genetically associated with these conditions. Therefore, further GWAS studies with more precise assessment of smoking behavior are warranted. Furthermore, our study also revealed the causal relationship between obesity and asthma, smoking and asthma, tobacco smoke exposure outside home and COPD, smoking and COPD. Therefore, smoking and obesity may act as mediators between GERD and chronic respiratory diseases. After using MVMR analyses adjusted by smoking and obesity, none of the potential confounders mentioned above were statistically significant mediators and the direct causal estimates of genetically predicted GERD with the outcomes remained robust. Complementary multivariable MR analyses showed that none of the potential mediators had statistically significant effects on the GERD-chronic respiratory disease pathway.
In general, caution is required when interpreting the results of MR analysis due to the presence of pleiotropy. The third assumption in MR analysis is the exclusion restriction assumption, meaning that instrumental variables only affect the outcome through the exposure [51]. However, pleiotropy can occur, causing instrumental variables to affect both the outcome variable and other variables through different mechanisms, leading to biased results in MR analysis. Therefore, excluding pleiotropy is a critical prerequisite for MR analysis. While it is challenging to entirely eliminate pleiotropy in practical studies, various sensitivity analysis methods can be used to mitigate its effects and enhance the reliability of the results. As mentioned above, MR-Egger regression and MR-PRESSO method can be applied to evaluate whether the instrumental variable selection results in any potential bias that could substantially impact the results [36, 37]. In our present study, there is a little evidence of pleiotropy for asthma, IPF, and COPD, with P-values ranging from > 0.01 to < 0.05 using the Egger intercept analyses. Meanwhile, MR-PRESSO tests identified two outliers for asthma and one outlier for COPD, but no significant statistical evidence of distortion was found through the distortion tests. Overall, we can conclude from these analyses that horizontal pleiotropy has little or no interference with our causal estimates. It is essential to note that while these methods can alleviate the effects of pleiotropy, they cannot entirely eliminate the possibility of pleiotropic effects. Therefore, the results should be interpreted with caution.
Possible potential mechanisms that have been advanced to explain the causality between GERD and the chronic respiratory diseases we mentioned above. The reflux theory and the reflex15 theory as two main mechanisms have been suggested in the association between GERD and asthma. On the one hand, the reflux theory suggests that micro-aspiration of gastric contents from the proximal esophagus can cause mucosal reaction in the larynx and lungs, leading to respiratory symptoms [52]. On the other hand, the reflex theory describes an indirect mechanism in which distal oesophageal reflux leads to bronchoconstriction secondary to a vagal reflex [53]. Although the reflex theory requires further validation, it is not surprising considering that the esophagus and trachea share the same embryonic origin and innervation by the vagus nerve [54]. The mechanisms underlying the relationship between GERD and COPD and IPF were similar to those in asthma. Furthermore, an increasing number of supporting evidence were detected through the measurement of pepsin and bile salts in broncho alveolar lavage fluid (BALF). Pauwels et al. demonstrated that eight of 29 patients (28%) with asthma were found the presence of bile acids in sputum [55]. Similarly, Savarino E et al. reported higher bile acids and pepsin (P < 0.03) in BALF (62% and 67%, respectively) and saliva (61% and 68%, respectively) of IPF patients than that of non-IPF patients (25% and 25% in BALF, and 33% and 36%, respectively, in saliva) and controls (0% and 0% in BALF and saliva, respectively) [47]. Gastric fluid has been shown to induce the expression of matrix metalloproteinase-2 (MMP-2) and MMP-9, as well as activation of the nuclear factor-КB (NF-КB) signaling pathway in macrophages, which can subsequently induce pulmonary fibroblast differentiation [56]. Previous reports have suggested the potentially crucial role of IL-6 produced by macrophages in fibrotic progression and lung disease [57]. Elevated IL-6 expression has been observed in pulmonary fibroblast cells, suggesting that IL-6 may play a potential autocrine or paracrine role in gastric fluid-induced inflammation in the pulmonary microenvironment [58]. In the pathogenesis of asthma, vitamin D is believed to affect lung development, regulate immune responses, and remodel airway smooth muscle [59]. 1,25(OH)2D3, the active form of vitamin D, has been shown to down-regulate the expression of MMP-9 and improve airway remodeling [60]. In contrast, we did not find a causal association between genetically predicted GERD and vitamin D deficiency, suggesting that GERD may affect smooth muscle cells through remodeling the expression of MMP-2 instead of MMP-9 [57]. Further investigations are needed to confirm these hypotheses.
There are several strengths in our study. First, we conducted the first MR study to assess the causality between GERD and chronic respiratory diseases as well as their risk factors. Participants were grouped according to their randomly allocated genotype, and this procedure mimicked a randomized controlled trial, adding rigor to our findings. In addition, the MR design helps to mitigate issues such as reverse causation and potential confounding factors that are commonly present in conventional observational studies. Secondly, we obtained robustly associated instrumental variables (IVs) from studies with large samples, ensuring that our MR analyses had sufficient statistical power and could provide comparatively accurate estimates of causality.
Several shortcomings in our study cannot be neglected. Firstly, the choice of SNPs and related consortiums were of European origin, which raises concerns about the generalizability of our findings to other populations and regions. Furthermore, the replication group from the FinnGen study did not observe a positive association between GWRD and IPF, which suggests that the current causal results may not be generalizable to them. Additionally, we observed heterogeneity between asthma and COPD in our study. As mentioned previously, the GWASs from which we obtained the summary data of asthma and COPD accepted a wide range of case definitions, which could contribute to the observed heterogeneity. Additionally, we found some sample overlap (UK Biobank mainly) in the GWASs of GERD. Apart from that, comprehensive time-to-event analysis were not conducted to establish the temporal association between GERD and chronic respiratory diseases. Unfortunately, there are no established guidelines for dealing with time-to-event data within an MR context. One analysis using the second stage regression method discussed the relationship between BMI and breast cancer survival, in which survival was considered as the time-to-event outcome [61]. The second stage regression is substituted with either a Cox proportional hazard, or additive hazard regression, which required the outcome to make the Cox model valid [61]. However, the study from which we collected the summary statistics did not provide the necessary data for such analysis. We look forward to future research that will further investigate time-to-event traits and provide more guidance for conducting MR analyses in this context. Other than that, chronic pulmonary diseases have so complex clinical manifestations that its diagnosis should rely on not only used standard ICD codes, but also other examinations or drugs. The diagnosis criteria of our study were cited from corresponding GWAS, which might be inadequate to diagnosis. Future large scale GWAS with more specific inclusion criteria are wanted. Instrumental variables for 12 common risk factors including smoking status (previous and current), exposure to tobacco smoke at home and outside home, obesity class 1–3, alcohol consumption, Type 2 diabetes, Vitamin D deficiency and lung function (FEV1/FVC) were selected from the open MR-base platform ,which all passed the genome-wide significance threshold of p < 5 × 10− 8 and the LD analysis with a threshold of R2 < 0.001 [21]. Nevertheless, IVs that were not included in the platform need to be explored in future studies.
In conclusion, our study provides novel evidence suggesting that individuals with a genetic predisposition to GERD were associated with an increased risk of asthma, IPF, COPD, and chronic bronchitis, while no significant association was observed between GERD and bronchiectasis. The MR results indicated that the micro-aspiration of GERD might play a role in the development of chronic respiratory diseases. Additionally, further large-scale studies are needed to confirm our findings and explore the epidemiological and mechanistic interactions between GERD and chronic respiratory diseases.
Data Availability
The datasets used and analyzed during the current study available from the corresponding author on reasonable request. R code for the methods proposed in this paper can be found in the package mr.raps that is publicly available at https://github.com/chengbrick/mendelian-randomization-between-GERD-AND-Chronic-respiratory-diseases.git.
Abbreviations
- GERD:
-
Gastroesophageal reflux disease
- SNP:
-
Single-nucleotide polymorphism
- IPF:
-
Idiopathic pulmonary fibrosis
- COPD:
-
Chronic obstructive disease
- OR:
-
Odds ratio
- CI:
-
Confidence interval
- GWASs:
-
Genome-wide association studies
- IVW:
-
Inverse variance-weighted
- MR:
-
Mendelian randomization
- IVs:
-
Instrumental variables
- LD:
-
Linkage disequilibrium
References
Katz PO, Gerson LB, Vela MF. Guidelines for the diagnosis and management of gastroesophageal reflux disease. Am J Gastroenterol. 2013;108(3):308–28. quiz 29.
Dent J, El-Serag HB, Wallander MA, Johansson S. Epidemiology of gastro-oesophageal reflux disease: a systematic review. Gut. 2005;54(5):710–7.
El-Serag HB, Sweet S, Winchester CC, Dent J. Update on the epidemiology of gastro-oesophageal reflux disease: a systematic review. Gut. 2014;63(6):871–80.
Savarino E, Bazzica M, Zentilin P, Pohl D, Parodi A, Cittadini G, et al. Gastroesophageal reflux and pulmonary fibrosis in scleroderma: a study using pH-impedance monitoring. Am J Respir Crit Care Med. 2009;179(5):408–13.
Mallah N, Turner JM, González-Barcala FJ, Takkouche B. Gastroesophageal reflux disease and asthma exacerbation: a systematic review and meta-analysis. Pediatr allergy immunology: official publication Eur Soc Pediatr Allergy Immunol. 2022;33(1):e13655.
Thakkar K, Boatright RO, Gilger MA, El-Serag HB. Gastroesophageal reflux and asthma in children: a systematic review. Pediatrics. 2010;125(4):e925–30.
Huang C, Liu Y, Shi G. A systematic review with meta-analysis of gastroesophageal reflux disease and exacerbations of chronic obstructive pulmonary disease. BMC Pulm Med. 2020;20(1):2.
Baqir M, Vasirreddy A, Vu AN, Moua T, Chamberlain AM, Frank RD, et al. Idiopathic pulmonary fibrosis and gastroesophageal reflux disease: a population-based, case-control study. Respir Med. 2021;178:106309.
Bédard Méthot D, Leblanc É, Lacasse Y. Meta-analysis of gastroesophageal reflux Disease and Idiopathic Pulmonary Fibrosis. Chest. 2019;155(1):33–43.
Eusebi LH, Ratnakumaran R, Yuan Y, Solaymani-Dodaran M, Bazzoli F, Ford AC. Global prevalence of, and risk factors for, gastro-oesophageal reflux symptoms: a meta-analysis. Gut. 2018;67(3):430–40.
Raghu G, Ley B, Brown KK, Cottin V, Gibson KF, Kaner RJ, et al. Risk factors for disease progression in idiopathic pulmonary fibrosis. Thorax. 2020;75(1):78–80.
Sekula P, Del Greco MF, Pattaro C, Köttgen A. Mendelian randomization as an Approach to assess causality using Observational Data. J Am Soc Nephrology: JASN. 2016;27(11):3253–65.
Lawlor DA, Harbord RM, Sterne JA, Timpson N, Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. 2008;27(8):1133–63.
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.
Pierce BL, Burgess S. Efficient design for mendelian randomization studies: subsample and 2-sample instrumental variable estimators. Am J Epidemiol. 2013;178(7):1177–84.
Ong JS, An J, Han X, Law MH, Nandakumar P, Schumacher J, et al. Multitrait genetic association analysis identifies 50 new risk loci for gastro-oesophageal reflux, seven new loci for Barrett’s oesophagus and provides insights into clinical heterogeneity in reflux diagnosis. Gut. 2022;71(6):1053–61.
Sun X, Chen L, Zheng L. A mendelian randomization study to assess the genetic liability of gastroesophageal reflux Disease for Cardiovascular Diseases and Risk factors. Human molecular genetics; 2022.
Valette K, Li Z, Bon-Baret V, Chignon A, Bérubé JC, Eslami A, et al. Prioritization of candidate causal genes for asthma in susceptibility loci derived from UK Biobank. Commun biology. 2021;4(1):700.
Allen RJ, Stockwell A, Oldham JM, Guillen-Guio B, Schwartz DA, Maher TM, et al. Genome-wide association study across five cohorts identifies five novel loci associated with idiopathic pulmonary fibrosis. Thorax. 2022;77(8):829–33.
Kurki MI, Karjalainen J, Palta P, Sipilä TP, Kristiansson K, Donner K et al. FinnGen: unique genetic insights from combining isolated population and national health register data. 2022:2022.03.03.22271360.
Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D et al. The MR-Base platform supports systematic causal inference across the human phenome. eLife. 2018;7.
Liu M, Jiang Y, Wedow R, Li Y, Brazel DM, Chen F, et al. Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Nat Genet. 2019;51(2):237–44.
Berndt SI, Gustafsson S, Mägi R, Ganna A, Wheeler E, Feitosa MF, et al. Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture. Nat Genet. 2013;45(5):501–12.
Clarke TK, Adams MJ, Davies G, Howard DM, Hall LS, Padmanabhan S, et al. Genome-wide association study of alcohol consumption and genetic overlap with other health-related traits in UK Biobank (N = 112 117). Mol Psychiatry. 2017;22(10):1376–84.
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.
Shrine N, Guyatt AL, Erzurumluoglu AM, Jackson VE, Hobbs BD, Melbourne CA, et al. New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries. Nat Genet. 2019;51(3):481–93.
Emdin CA, Khera AV, Kathiresan S. Mendelian Randomization Jama. 2017;318(19):1925–6.
Collard HR, Ryerson CJ, Corte TJ, Jenkins G, Kondoh Y, Lederer DJ, et al. Acute Exacerbation of Idiopathic Pulmonary Fibrosis. An International Working Group Report. Am J Respir Crit Care Med. 2016;194(3):265–75.
Lederer DJ, Martinez FJ. Idiopathic pulmonary fibrosis. N Engl J Med. 2018;378(19):1811–23.
Stern J, Pier J, Litonjua AA. Asthma epidemiology and risk factors. Semin Immunopathol. 2020;42(1):5–15.
Wang G, Hallberg J, Um Bergström P, Janson C, Pershagen G, Gruzieva O et al. Assessment of chronic bronchitis and risk factors in young adults: results from BAMSE. Eur Respir J. 2021;57(3).
Yang IA, Jenkins CR, Salvi SS. Chronic obstructive pulmonary disease in never-smokers: risk factors, pathogenesis, and implications for prevention and treatment. The Lancet Respiratory medicine. 2022;10(5):497–511.
Burgess S, Thompson SG. Multivariable mendelian randomization: the use of pleiotropic genetic variants to estimate causal effects. Am J Epidemiol. 2015;181(4):251–60.
Brion MJ, Shakhbazov K, Visscher PM. Calculating statistical power in mendelian randomization studies. Int J Epidemiol. 2013;42(5):1497–501.
Burgess S, Bowden J, Fall T, Ingelsson E, Thompson SG. Sensitivity analyses for robust causal inference from mendelian randomization analyses with multiple genetic variants. Epidemiol (Cambridge Mass). 2017;28(1):30–42.
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.
Burgess S, Thompson SG. Interpreting findings from mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017;32(5):377–89.
Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan N, Thompson J. A framework for the investigation of pleiotropy in two-sample summary data mendelian randomization. Stat Med. 2017;36(11):1783–802.
CosÃo BG, Shafiek H, Mosteiro M, Iglesias A, Gómez C, Toledo-Pons N et al. Redefining the role of bronchoscopy in the workup of severe uncontrolled asthma in the era of Biologics: a prospective study. Chest. 2023.
Tobin RW, Pope CE 2nd, Pellegrini CA, Emond MJ, Sillery J, Raghu G. Increased prevalence of gastroesophageal reflux in patients with idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 1998;158(6):1804–8.
Lee CM, Lee DH, Ahn BK, Hwang JJ, Yoon H, Shin CM, et al. Protective effect of Proton Pump inhibitor for survival in patients with gastroesophageal reflux Disease and Idiopathic Pulmonary Fibrosis. J Neurogastroenterol Motil. 2016;22(3):444–51.
Kim J, Lee JH, Kim Y, Kim K, Oh YM, Yoo KH, et al. Association between chronic obstructive pulmonary disease and gastroesophageal reflux disease: a national cross-sectional cohort study. BMC Pulm Med. 2013;13:51.
Martinez CH, Okajima Y, Murray S, Washko GR, Martinez FJ, Silverman EK, et al. Impact of self-reported gastroesophageal reflux disease in subjects from COPDGene cohort. Respir Res. 2014;15(1):62.
McDonnell MJ, O’Toole D, Ward C, Pearson JP, Lordan JL, De Soyza A, et al. A qualitative synthesis of gastro-oesophageal reflux in bronchiectasis: current understanding and future risk. Respir Med. 2018;141:132–43.
Duncan DR, Cohen A, Golden C, Lurie M, Mitchell PD, Liu E, et al. Gastrointestinal factors associated with risk of bronchiectasis in children. Pediatr Pulmonol. 2023;58(3):899–907.
Lee AL, Button BM, Denehy L, Roberts S, Bamford T, Mu FT, et al. Exhaled Breath Condensate Pepsin: potential noninvasive test for gastroesophageal reflux in COPD and Bronchiectasis. Respir Care. 2015;60(2):244–50.
Lee AL, Button BM, Denehy L, Roberts SJ, Bamford TL, Ellis SJ, et al. Proximal and distal gastro-oesophageal reflux in chronic obstructive pulmonary disease and bronchiectasis. Respirol (Carlton Vic). 2014;19(2):211–7.
Vakil N, van Zanten SV, Kahrilas P, Dent J, Jones R. The Montreal definition and classification of gastroesophageal reflux disease: a global evidence-based consensus. Am J Gastroenterol. 2006;101(8):1900–20. quiz 43.
Zhang Y, Zhao M, Guo P, Wang Y, Liu L, Zhao J, et al. Mendelian randomisation highlights hypothyroidism as a causal determinant of idiopathic pulmonary fibrosis. EBioMedicine. 2021;73:103669.
Cotton C, Alton P, Hughes DM, Zhao SS. Genetic liability to gastro-esophageal reflux disease, obesity, and risk of idiopathic pulmonary fibrosis. Respiratory Invest. 2023;61(3):335–8.
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 research ed). 2021;375:n2233.
Raghavendran K, Nemzek J, Napolitano LM, Knight PR. Aspiration-induced lung injury. Crit Care Med. 2011;39(4):818–26.
Kasasbeh A, Kasasbeh E, Krishnaswamy G. Potential mechanisms connecting asthma, esophageal reflux, and obesity/sleep apnea complex–a hypothetical review. Sleep Med Rev. 2007;11(1):47–58.
Hom C, Vaezi MF. Extra-esophageal manifestations of gastroesophageal reflux disease: diagnosis and treatment. Drugs. 2013;73(12):1281–95.
Pauwels A, Decraene A, Blondeau K, Mertens V, Farre R, Proesmans M, et al. Bile acids in sputum and increased airway inflammation in patients with cystic fibrosis. Chest. 2012;141(6):1568–74.
Savarino E, Carbone R, Marabotto E, Furnari M, Sconfienza L, Ghio M, et al. Gastro-oesophageal reflux and gastric aspiration in idiopathic pulmonary fibrosis patients. Eur Respir J. 2013;42(5):1322–31.
Chiang CC, Chen CM, Suen JL, Su HH, Hsieh CC, Cheng CM. Stimulatory effect of gastroesophageal reflux disease (GERD) on pulmonary fibroblast differentiation. Digestive and liver disease: official journal of the italian society of Gastroenterology and the Italian Association for the study of the liver. 2020;52(9):988–94.
Camelo A, Dunmore R, Sleeman MA, Clarke DL. The epithelium in idiopathic pulmonary fibrosis: breaking the barrier. Front Pharmacol. 2014;4:173.
Salmanpour F, Kian N, Samieefar N, Khazeei Tabari MA, Rezaei N. Asthma and vitamin D Deficiency: occurrence, Immune Mechanisms, and New Perspectives. J Immunol Res. 2022;2022:6735900.
Berraies A, Hamzaoui K, Hamzaoui A. Link between vitamin D and airway remodeling. J asthma allergy. 2014;7:23–30.
Guo Q, Burgess S, Turman C, Bolla MK, Wang Q, Lush M, et al. Body mass index and breast cancer survival: a mendelian randomization analysis. Int J Epidemiol. 2017;46(6):1814–22.
Acknowledgements
The authors acknowledge the efforts of the consortia in providing high-quality GWAS resources for researchers.
Funding
We did not receive any funding for this study.
Author information
Authors and Affiliations
Contributions
Xiaoxue Cheng: Conceptualization, Data curation, Formal analysis, Methodology, Writing - original draft, Writing - review & editing; Jiang Shi: Conceptualization, Methodology, Writing - Review & Editing; Ding Zhang: Data curation, Writing - Review & Editing; Caichen Li: Writing - Review & Editing; Jianxing He: Supervision, Funding acquisition, Project administration, Writing - Review & Editing; Wenhua Liang: Supervision, Funding acquisition, Project administration, Writing - Review & Editing. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable, since the study is based on summary-level data. In all original studies, ethical approval and consent to participate were obtained.
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.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Cheng, X., Shi, J., Zhang, D. et al. Assessing the genetic relationship between gastroesophageal reflux disease and chronic respiratory diseases: a mendelian randomization study. BMC Pulm Med 23, 243 (2023). https://doi.org/10.1186/s12890-023-02502-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12890-023-02502-8