Skip to main content

ORIGINAL RESEARCH article

Front. Med., 20 June 2024
Sec. Rheumatology

The causal relationship between rheumatoid arthritis and bronchiectasis: a bidirectional Mendelian randomization study

Yuanyuan Li
&#x;Yuanyuan Li1*Weina Wang&#x;Weina Wang1Dengfeng ZhouDengfeng Zhou1Lili LiLili Li2
  • 1Department of Respiratory and Critical Care Medicine, Wuhan Fourth Hospital, Wuhan, Hubei, China
  • 2Department of Gastroenterology, Wuhan Fourth Hospital, Wuhan, Hubei, China

Background: Several observational studies suggested an association between rheumatoid arthritis (RA) and bronchiectasis. Nevertheless, the presence of a causal relationship between these conditions is yet to be determined. This study aimed to investigate whether genetically predicted RA is associated with the risk of bronchiectasis and vice versa.

Methods: We obtained RA genome-wide association study (GWAS) data from FinnGen consortium, and bronchiectasis GWAS data from IEU Open GWAS project. Univariate Mendelian randomization (MR) analysis was performed using inverse variance weighted (IVW) estimation as the main method. Furthermore, bidirectional and replication MR analysis, multivariate MR (MVMR), Mediation analysis, and sensitivity analyses were conducted to validate the findings.

Results: In the UVMR analysis, the IVW results revealed that RA had an increased risk of bronchiectasis (OR = 1.18, 95% CI = 1.10–1.27; p = 2.34 × 10−6). In the reverse MR analysis, no evidence of a causal effect of bronchiectasis on the risk of RA was detected. Conversely, in the replication MR analysis, RA remained associated with an increased risk of bronchiectasis. Estimates remained consistent in MVMR analyses after adjusting for the prescription of non-steroidal anti-inflammatory drugs (NSAIDs) and glucocorticoids. Immunosuppressants were found to mediate 58% of the effect of the RA on bronchiectasis. Sensitivity analyses confirmed the stability of these associations.

Conclusion: This study demonstrated a positive causal relationship between RA and an increased risk of bronchiectasis, offering insights for the early prevention of bronchiectasis in RA patients and shedding new light on the potential role of immunosuppressants as mediators in promoting the effects of RA on bronchiectasis.

Introduction

Rheumatoid arthritis (RA) is a prevalent autoimmune and chronic inflammatory disease (1), affecting approximately 0.3 to 1% of the population in developed countries (2). RA is a chronic autoimmune disorder that primarily affects the joints but is also associated with several well-documented extra-articular manifestations, including ocular, cutaneous, gastrointestinal, cardiac, and pulmonary involvements (3). RA has the potential to impact various anatomical regions of the chest, encompassing the lung parenchyma, both large and small airways, the pleura, and, less commonly, the pulmonary vessels (4, 5). Bronchiectasis is a recognized extra-articular manifestation of RA, characterized by irreversible damage, widening, and thickening of the bronchi, which results in excessive mucus secretion (6, 7). The presence of both bronchiectasis and RA is associated with a diminished quality of life and an elevated risk of infection for patients with RA (8). Additionally, the coexistence of bronchiectasis and RA markedly impacts patients’ health and life expectancy, elevating their risk of mortality. Patients with both conditions face a mortality risk that is 7.3 times higher than that of the general population, five times greater than that of those with RA alone, and 2.4 times higher than that of those with bronchiectasis alone (9). Accordingly, understanding the interplay between bronchiectasis and RA and the underlying mechanisms that guide their relationship is crucial.

The coexistence of bronchiectasis in patients with RA has been recognized for several decades. The link between RA and bronchiectasis was initially delineated in 1955 (10). Since then, many observation studies have shown a close relationship between RA and bronchiectasis. A recent meta-analysis including 23 observational studies revealed that the prevalence of bronchiectasis in RA is 21.1% (95% CI: 15.0–28.9%), as estimated using a random-effects model. The pooled prevalence for clinically defined bronchiectasis was 2.69% (95% CI: 1.63–4.42), whereas that for radiologically identified disease was 24.9% (95% CI: 19.21–31.67). However, the majority of studies have been cross-sectional, enrolling a limited number of patients and lacking control groups (11). Consequently, the direction of the association between these two diseases remains ambiguous. Furthermore, these observational studies are prone to confounding factors and the possibility of reverse causation, leaving it uncertain whether having RA elevates the risk of bronchiectasis.

Mendelian randomization (MR) has garnered significant interest as a tool for elucidating causal relationships between risk factors and disease outcomes. MR capitalizes on random genetic variants, which arise from meiosis, to serve as instrumental variables (IVs) for examining the link between environmental exposures and disease development. The random assignment of genetic variants at conception, preceding the onset of disease, allows MR analyses to mitigate confounding and identify the causal determinants of specific outcomes (1214). Recent MR studies have confirmed the causal association of RA with pre-eclampsia (15), hypothyroidism (16), lung cancer, etc. (17). Therefore, this study collected published data and employed bidirectional MR analysis to ascertain whether a bidirectional causal relationship exists between RA and bronchiectasis.

Materials and methods

Study design

This study was conducted following the Strengthening the Reporting of Observational Studies in Epidemiology using MR (STROBE-MR) guidance (Supplementary Table S1) (18). This study did not necessitate a separate ethical review submission. Employing summary statistics from genome-wide association studies (GWAS), we performed two-sample MR analyses to estimate the causal impact of RA on bronchiectasis within European populations. IVs were selected based on three crucial criteria: they should exhibit a strong association with the exposure of interest, be independent of unmeasured confounders, and affect the results exclusively through their relationship with the exposure (19). In this MR study, RA was considered the exposure, while bronchiectasis was the outcome of interest. We initially conducted univariable MR (UVMR) to investigate the potential causal relationship between RA and bronchiectasis. Subsequently, we employed reverse and replication MR analyses to corroborate this relationship. Furthermore, we performed multivariable MR (MVMR) to evaluate the direct effect of RA on bronchiectasis, accounting for potential confounders such as the prescription of non-steroidal anti-inflammatory drugs (NSAIDs), glucocorticoids and immunosuppressants. A concise illustration of the MR study design is provided in Figure 1.

Figure 1
www.frontiersin.org

Figure 1. Assumptions and study design of the bidirectional MR study of the associations between RA and bronchiectasis. RA, rheumatoid arthritis; MR, Mendelian randomization; SNPs, single-nucleotide polymorphisms.

Data sources

SNPs associated with RA in European populations were identified from a comprehensive GWAS, which included 6,236 RA cases and 147,221 controls from the FinnGen consortium. For bronchiectasis, GWAS data from a European population were sourced from the IEU Open GWAS project, consisting of 2,888 cases and 440,263 controls (20).

Selection of IVs

In this MR study, single nucleotide polymorphisms (SNPs) were considered significant at the genome-wide level (p < 5.0 × 10−8). SNPs that exhibited no linkage disequilibrium (LD) with other SNPs (r2 < 0.001 within a 10 Mb clustering window) were selected as IVs. The F-statistic was computed to evaluate the extent of genetic variation, and SNPs with an F-statistic below 10 were excluded, as they suggested insufficient strength (21). For the missing SNPs in the resulting GWAS dataset, proxies were identified at the cut-off of R2 > 0.8.1 SNPs lacking suitable proxies were removed from the analysis. Harmonization processes were conducted to eliminate ambiguous and palindromic SNPs, as well as those exhibiting high correlation with known results (p < 5.0 × 10−8). Prior to each MR analysis, MR-pleiotropy residual sum and outlier (MR-PRESSO) tests were conducted. These tests are intended to identify horizontal pleiotropy in fewer than 50% of instruments, allowing for the detection and exclusion of potential outliers (22).

UVMR

We carried out an MR analysis to estimate the causal effect of RA on bronchiectasis risk. When a single SNP served as the IV, the Wald ratio was employed. In instances with more than two SNPs, the inverse variance weighted (IVW) estimation method was adopted as the primary analysis approach, which involved combining the individual SNPs with their corresponding Wald ratios to derive a pooled estimate of causality. This technique accommodates overdispersion (23). Additionally, complementary MR analyses, including MR-Egger regression and weighted median, were conducted to enhance the reliability of the estimates provided by the IVW method across a broader spectrum of scenarios. The MR-Egger regression is particularly adept at quantifying multinomial and substantial heterogeneity in exposure imbalances; however, it requires a larger sample size to detect comparable levels of variability in under-exposure when compared to alternative methods (24). When at least half of the weighted variance is due to horizontal pleiotropy and is valid, the weighted median method yields consistent effect estimates (25). The P < 0.05 (two-tailed) was statistically significant.

Bidirectional and replication MR analysis

To evaluate the bi-directional causation effects between RA and bronchiectasis, we used bronchiectasis as “exposure” and RA as “outcome.” In addition, to prove the reliability of our findings, we replicated the IVW analysis using another independent RA GWAS data from the IEU Open GWAS project, including 14,361 RA patients, 43,923 controls, as well as 13,108,512SNPs (26). Genetic variants associated with bronchiectasis from the FinnGen consortium, which contained 1,107 bronchiectasis patients, 186,723 controls, as well as 16,380,375 SNPs.

Similarly, we also validated our findings in an East Asian population. SNPs associated with RA were identified from a large GWAS involving 4,199 cases of East Asian descent and 208,254 controls. We further retrieved genetic data related to bronchiectasis in East Asian populations from the IEU Open GWAS project, which included 241 cases and 161,803 controls. The method for performing MR analyses in reverse and replication MR is the same as previously described.

MVMR

MVMR, using the MVMR-IVW method, was conducted to minimize the correlated pleiotropy introduced by potential confounding factors, including the prescription of NSAIDs, glucocorticoids, and immunosuppressants (27). In MVMR, we obtained summary data for NSAIDs (n = 164,520), glucocorticoids (n = 205,700), and immunosuppressants (n = 272,602) (20) from the IEU Open GWAS project.

Mediation analysis

We further conducted a mediation analysis using the two-step MR approach to investigate whether the causal effect of RA on bronchiectasis is mediated by the use of immunosuppressants. In the initial step, we determined the causal effect of RA on immunosuppressants (β1). In the second step, we estimated the causal effect of immunosuppressants on the development of bronchiectasis (β2). The mediation effect is calculated as β1*β2. Thus, the proportion mediated could be calculated as β1*β2/the overall effect (β).

Sensitivity analysis

To further assess the stability of the results, we also conducted Cochran’s Q statistic, MR-Egger intercept, funnel plot, and leave-one-out (LOO) analyses to detect the presence of heterogeneity and pleiotropy. Cochran’s Q-test (p < 0.05 indicates heterogeneity) and I2 statistic (I2-value > 50% indicates heterogeneity) were utilized to evaluate the heterogeneity among SNPs in the IVW estimates. A p-value of <0.05 was deemed indicative of significant heterogeneity, prompting the use of a random-effects model for subsequent analyses. The intercept term in the MR-Egger regression was employed to assess the presence of pleiotropy (25, 2830). All of the analyses were run by using the R package TwoSampleMR (version 0.5.6) in R (version 4.1.1).

Results

UVMR

In the UVMR analysis, 12 IVs were associated with the risk of RA following a series of quality control measures. All F-values for the inclusion of SNPs exceeded 10 (Supplementary Table S2). As depicted in Table 1, the primary results from the IVW method revealed a positive association between RA and the risk of bronchiectasis [odds ratio (OR) = 1.18, 95% confidence interval (CI) = 1.10–1.27, p = 2.34 × 10−6]. Meanwhile, similar risk estimates were obtained using the MR-Egger regression (OR = 1.15, 95% CI =1.03–1.29, p = 0.0353) and weighted median approaches (OR = 1.16, 95% CI =1.08–1.24, p = 2.06 × 10−5) (Figure 2). The Cochran Q-derived p-values (p > 0.05) suggested the absence of heterogeneity, and the MR-Egger regression intercept provided no significant evidence of horizontal pleiotropy. The funnel plot was symmetrical (Figure 3). Additionally, the LOO analysis demonstrated that the effect estimates remained unchanged despite the exclusion of any single variant (Figure 4).

Table 1
www.frontiersin.org

Table 1. Causal relationships between rheumatoid arthritis (RA) and bronchiectasis risk performed by MR.

Figure 2
www.frontiersin.org

Figure 2. Scatter plot for the causal effect of rheumatoid arthritis on the risk of bronchiectasis.

Figure 3
www.frontiersin.org

Figure 3. Funnel plot for the causal effect of rheumatoid arthritis on the risk of bronchiectasis.

Figure 4
www.frontiersin.org

Figure 4. Leave-one-out analysis for the causal effect of rheumatoid arthritis on the risk of bronchiectasis.

Bidirectional and replication MR analysis

In the bidirectional MR analysis, there was no evidence of a causal effect of bronchiectasis on RA risk (OR = 1.30, 95% CI = 0.99–1.70, p = 0.0562) (Table 1). To further validate our results, we performed replication MR analysis using GWAS databases from different sources. Unsurprisingly, we observed a similar trend in European (OR = 1.13, 95% CI = 1.05–1.21, p = 0.0011) and East Asian populations (OR = 1.30, 95% CI = 1.07–1.58; p = 0.0076) (Table 1). Sensitivity analysis indicated the results were stable.

MVMR and mediation analysis

Several factors may affect the association between RA and the risk of bronchiectasis. MVMR analysis was performed to examine whether RA has independent effects on bronchiectasis with adjustment of potential confounding factors. In the MVMR analysis, IVW estimates were consistent with the adjustments of NSAIDs (OR = 1.19, 95% CI = 1.10–1.29, p = 1.16 × 10−5), and glucocorticoids (OR = 1.18, 95% CI = 1.11–1.26, p = 7.34 × 10−8) (Table 2). However, no evidence was shown between RA and bronchiectasis after adjustment for immunosuppressants (OR = 0.99, 95% CI = 0.88–1.10, p = 0.7942), suggesting that immunosuppressants might play a role in the association between RA and bronchiectasis.

Table 2
www.frontiersin.org

Table 2. Multivariable inverse variance weighted estimates for adjusted associations with RAs.

We further performed mediation analyses to validate the mediating effect of immunosuppressants in bronchiectasis due to RA. The results revealed that immunosuppressants account for 58% of the association between RA and bronchiectasis (Supplementary Table S3).

Discussion

To our knowledge, this is the first MR study to explore the causal relationship between RA and bronchiectasis. In the present study, a two-sample MR study demonstrated a causal relationship between genetically predicted RA and increased risk of bronchiectasis. Contrary to this, we found no significant association between bronchiectasis and RA. Both the MVMR and sensitivity analyses confirmed the statistical robustness of these findings.

Previous observational studies have suggested an association between RA and bronchiectasis. For example, several retrospective studies have shown that the prevalence of bronchiectasis in patients with RA ranges from 0.6 to 52.0% (31, 32), whereas the prevalence observed in prospective studies ranges from 1.9 to 58% (33, 34). The variations in prevalence rates can be largely attributed to the utilization of diverse imaging modalities, such as clinical scans, computed tomography (CT), high-resolution CT (HRCT), and additional radiological methods. In addition, several cross-sectional studies also provide some supporting evidence (3439). In the cross-sectional study by Metafratzi et al., they found a prevalence of 58.1% bronchiectasis in early RA patients (without respiratory symptoms). Although the control group exhibited a similar frequency of bronchiectasis, the extent of the condition, as measured by HRCT scores, was substantially more severe among the RA patients compared to the controls (34). However, these study designs are limited to establishing causal inferences due to the potential for confounding bias. Randomized controlled trials (RCTs) are the gold standard for establishing causality; however, they are often limited by substantial time and financial constraints. MR can be considered a naturalistic approximation of RCTs. Therefore, our MR study provides robust evidence for a causal relationship between RA and bronchiectasis.

RA and bronchiectasis are both common diseases with a poor prognosis, and RA combined with bronchiectasis has an even worse prognosis, However, the mechanisms underlying the association of these conditions are unclear (40, 41). Current research suggests that three primary factors are likely involved in the development of the condition: inflammation, infection, and genetics. In a longitudinal population-based study by Chao et al., RA patients had a 2.12-fold increased risk of developing bronchiectasis compared with matched control patients, which remained significant after accounting for potential confounders. This elevated risk was particularly pronounced in patients with seropositive RA (SPRA), suggesting that rheumatic inflammation contributes significantly to the development of RA-bronchiectasis (8). Although NSAIDs and glucocorticoids are recognized for their proclivity to increase the risk of infections, the relationship between these medications and the heightened risk of bronchiectasis in patients with RA remains under-explored. Our MVMR analysis results suggest that these factors may not act as confounders; however, additional clinical studies are necessary to corroborate our findings. Standard RA treatments, such as corticosteroids, leflunomide, and biologic disease-modifying antirheumatic drugs (DMARDs), have been associated with an elevated risk of lower respiratory tract infections and mycobacterial lung infections, potentially leading to a marked increase in the incidence of bronchiectasis. The results of our MVMR analysis also suggest that immunosuppressants may play a role in RA leading to bronchiectasis. Mediation analyses further revealed that 58% of the effects of RA on bronchiectasis were mediated by immunosuppressants (4244). Besides, several studies have also reported 6 genetic risk factors associated with the development of bronchiectasis, and these studies have focused on CFTR and HLA variants (4548).

However, our study also has several limitations. First, the nature of GWAS datasets precludes the availability of comprehensive clinical details, such as the type of immunosuppressant used in RA patients, and therefore we were unable to perform subgroup analyses for different immunosuppressants. Second, due to the lack of detailed information on clinically defined or radiologically defined bronchiectasis in our GWAS data, we were unable to conduct subgroup analyses pertaining to these factors. Third, although the MR method provides us with a way of estimating causal relationships between genes and phenotypes, we cannot be completely certain of the causal relationship between RA and bronchiectasis due to the limitations of the effect size and the number of SNPs. Lastly, caution should be exercised when generalizing these findings to diverse racial or ethnic groups, given the study’s limited racial and ethnic scope.

In summary, RA could increase the risk of bronchiectasis and is mediated by immunosuppression, but there is no evidence that bronchiectasis increases the risk of RA. The findings of this study serve as a foundational point for future exploration and hold clinical relevance. Healthcare providers should be vigilant about the potential onset of bronchiectasis in patients with RA. The underlying pathophysiological mechanisms linking the development of bronchiectasis in RA patients warrant further comprehensive investigation.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

Ethical approval was not required for the studies involving humans because this study did not require separate ethical approval. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants' legal guardians/next of kin in accordance with the national legislation and institutional requirements because this study did not require separate ethical approval.

Author contributions

YL: Conceptualization, Data curation, Formal analysis, Investigation, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. WW: Data curation, Investigation, Writing – original draft. DZ: Investigation, Methodology, Writing – original draft. LL: Investigation, Methodology, Writing – original draft.

Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

Acknowledgments

We thank all the authors, investigators of the FinnGen study, and participants of GWASs.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2024.1403851/full#supplementary-material

Footnotes

References

1. Smolen, JS, Aletaha, D, and McInnes, IB. Rheumatoid arthritis. Lancet. (2016) 388:2023–38. doi: 10.1016/s0140-6736(16)30173-8

Crossref Full Text | Google Scholar

2. Smith, E, Hoy, DG, Cross, M, Vos, T, Naghavi, M, Buchbinder, R, et al. The global burden of other musculoskeletal disorders: estimates from the global burden of disease 2010 study. Ann Rheum Dis. (2014) 73:1462–9. doi: 10.1136/annrheumdis-2013-204680

PubMed Abstract | Crossref Full Text | Google Scholar

3. Cojocaru, M, Cojocaru, IM, Silosi, I, Vrabie, CD, and Tanasescu, R. Extra-articular manifestations in rheumatoid arthritis. Maedica. (2010) 5:286–91.

PubMed Abstract | Google Scholar

4. Brown, KK. Rheumatoid lung disease. Proc Am Thorac Soc. (2007) 4:443–8. doi: 10.1513/pats.200703-045MS

PubMed Abstract | Crossref Full Text | Google Scholar

5. Yunt, ZX, and Solomon, JJ. Lung disease in rheumatoid arthritis. Rheum Dis Clin N Am. (2015) 41:225–36. doi: 10.1016/j.rdc.2014.12.004

PubMed Abstract | Crossref Full Text | Google Scholar

6. McShane, PJ, Naureckas, ET, and Strek, ME. Bronchiectasis in a diverse US population: effects of ethnicity on etiology and sputum culture. Chest. (2012) 142:159–67. doi: 10.1378/chest.11-1024

PubMed Abstract | Crossref Full Text | Google Scholar

7. Weycker, D, Hansen, GL, and Seifer, FD. Prevalence and incidence of noncystic fibrosis bronchiectasis among US adults in 2013. Chron Respir Dis. (2017) 14:377–84. doi: 10.1177/1479972317709649

PubMed Abstract | Crossref Full Text | Google Scholar

8. Choi, H, Han, K, Jung, JH, Park, J, Kim, BG, Yang, B, et al. Impact of rheumatoid arthritis and seropositivity on the risk of non-cystic fibrosis bronchiectasis. Chest. (2024) 165:1330–40. doi: 10.1016/j.chest.2024.01.001

PubMed Abstract | Crossref Full Text | Google Scholar

9. Swinson, DR, Symmons, D, Suresh, U, Jones, M, and Booth, J. Decreased survival in patients with co-existent rheumatoid arthritis and bronchiectasis. Br J Rheumatol. (1997) 36:689–91. doi: 10.1093/rheumatology/36.6.689

PubMed Abstract | Crossref Full Text | Google Scholar

10. Aronoff, A, Bywaters, EG, and Fearnley, GR. Lung lesions in rheumatoid arthritis. Br Med J. (1955) 2:228–32. doi: 10.1136/bmj.2.4933.228

PubMed Abstract | Crossref Full Text | Google Scholar

11. Wiater, R, Håkansson, KEJ, and Ulrik, CS. A causal relationship between rheumatoid arthritis and bronchiectasis? A systematic review and meta-analysis. Chron Respir Dis. (2021) 18:1479973121994565. doi: 10.1177/1479973121994565

PubMed Abstract | Crossref Full Text | Google Scholar

12. Emdin, CA, Khera, AV, and Kathiresan, S. Mendelian randomization. JAMA. (2017) 318:1925–6. doi: 10.1001/jama.2017.17219

Crossref Full Text | Google Scholar

13. Davey Smith, G, and Hemani, G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. (2014) 23:R89–98. doi: 10.1093/hmg/ddu328

PubMed Abstract | Crossref Full Text | Google Scholar

14. Smith, GD, and Ebrahim, S. 'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. (2003) 32:1–22. doi: 10.1093/ije/dyg070

PubMed Abstract | Crossref Full Text | Google Scholar

15. Zhang, D, Hu, Y, Guo, W, Song, Y, Yang, L, Yang, S, et al. Mendelian randomization study reveals a causal relationship between rheumatoid arthritis and risk for pre-eclampsia. Front Immunol. (2022) 13:1080980. doi: 10.3389/fimmu.2022.1080980

PubMed Abstract | Crossref Full Text | Google Scholar

16. Duan, L, Chen, D, Shi, Y, Ye, S, Dou, S, and Feng, Y. Rheumatoid arthritis and hypothyroidism: a bidirectional Mendelian randomization study. Front Immunol. (2023) 14:1146261. doi: 10.3389/fimmu.2023.1146261

PubMed Abstract | Crossref Full Text | Google Scholar

17. Wu, X, Peng, H, Wen, Y, Cai, X, Li, C, Zhong, R, et al. Rheumatoid arthritis and risk of lung cancer: meta-analysis and Mendelian randomization study. Semin Arthritis Rheum. (2021) 51:565–75. doi: 10.1016/j.semarthrit.2021.03.015

Crossref Full Text | Google Scholar

18. 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. (2021) 375:n2233. doi: 10.1136/bmj.n2233

PubMed Abstract | Crossref Full Text | Google Scholar

19. Boef, AG, Dekkers, OM, and le Cessie, S. Mendelian randomization studies: a review of the approaches used and the quality of reporting. Int J Epidemiol. (2015) 44:496–511. doi: 10.1093/ije/dyv071

PubMed Abstract | Crossref Full Text | Google Scholar

20. Sakaue, S, Kanai, M, Tanigawa, Y, Karjalainen, J, Kurki, M, Koshiba, S, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet. (2021) 53:1415–24. doi: 10.1038/s41588-021-00931-x

PubMed Abstract | Crossref Full Text | Google Scholar

21. Pierce, BL, Ahsan, H, and Vanderweele, TJ. Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants. Int J Epidemiol. (2011) 40:740–52. doi: 10.1093/ije/dyq151

PubMed Abstract | Crossref Full Text | Google Scholar

22. Verbanck, M, Chen, CY, Neale, B, and Do, R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. (2018) 50:693–8. doi: 10.1038/s41588-018-0099-7

PubMed Abstract | Crossref Full Text | Google Scholar

23. Pierce, BL, and Burgess, S. Efficient design for Mendelian randomization studies: subsample and 2-sample instrumental variable estimators. Am J Epidemiol. (2013) 178:1177–84. doi: 10.1093/aje/kwt084

PubMed Abstract | Crossref Full Text | Google Scholar

24. Bowden, J, Davey Smith, G, and Burgess, S. Mendelian randomization with invalid instruments: effect estimation and bias detection through egger regression. Int J Epidemiol. (2015) 44:512–25. doi: 10.1093/ije/dyv080

PubMed Abstract | Crossref Full Text | Google Scholar

25. Bowden, J, Davey Smith, G, Haycock, PC, and Burgess, S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. (2016) 40:304–14. doi: 10.1002/gepi.21965

PubMed Abstract | Crossref Full Text | Google Scholar

26. Ha, E, Bae, SC, and Kim, K. Large-scale meta-analysis across east Asian and European populations updated genetic architecture and variant-driven biology of rheumatoid arthritis, identifying 11 novel susceptibility loci. Ann Rheum Dis. (2021) 80:558–65. doi: 10.1136/annrheumdis-2020-219065

Crossref Full Text | Google Scholar

27. Sanderson, E, Davey Smith, G, Windmeijer, F, and Bowden, J. An examination of multivariable Mendelian randomization in the single-sample and two-sample summary data settings. Int J Epidemiol. (2019) 48:713–27. doi: 10.1093/ije/dyy262

PubMed Abstract | Crossref Full Text | Google Scholar

28. Bowden, J, Del Greco, MF, Minelli, C, Davey Smith, G, Sheehan, N, and Thompson, J. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med. (2017) 36:1783–802. doi: 10.1002/sim.7221

PubMed Abstract | Crossref Full Text | Google Scholar

29. Greco, MF, Minelli, C, Sheehan, NA, and Thompson, JR. Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome. Stat Med. (2015) 34:2926–40. doi: 10.1002/sim.6522

PubMed Abstract | Crossref Full Text | Google Scholar

30. Higgins, JP, and Thompson, SG. Quantifying heterogeneity in a meta-analysis. Stat Med. (2002) 21:1539–58. doi: 10.1002/sim.1186

Crossref Full Text | Google Scholar

31. Shadick, NA, Fanta, CH, Weinblatt, ME, O'Donnell, W, and Coblyn, JS. Bronchiectasis: a late feature of severe rheumatoid arthritis. Medicine. (1994) 73:161–70. doi: 10.1097/00005792-199405000-00005

Crossref Full Text | Google Scholar

32. Akira, M, Sakatani, M, and Hara, H. Thin-section CT findings in rheumatoid arthritis-associated lung disease: CT patterns and their courses. J Comput Assist Tomogr. (1999) 23:941–8. doi: 10.1097/00004728-199911000-00021

Crossref Full Text | Google Scholar

33. Morrison, SC, Mody, GM, Benatar, SR, and Meyers, OL. The lungs in rheumatoid arthritis--a clinical, radiographic and pulmonary function study. S Afr Med J. (1996) 86:829–33.

PubMed Abstract | Google Scholar

34. Metafratzi, ZM, Georgiadis, AN, Ioannidou, CV, Alamanos, Y, Vassiliou, MP, Zikou, AK, et al. Pulmonary involvement in patients with early rheumatoid arthritis. Scand J Rheumatol. (2007) 36:338–44. doi: 10.1080/03009740701393957

Crossref Full Text | Google Scholar

35. Hassan, WU, Keaney, NP, Holland, CD, and Kelly, CA. High resolution computed tomography of the lung in lifelong non-smoking patients with rheumatoid arthritis. Ann Rheum Dis. (1995) 54:308–10. doi: 10.1136/ard.54.4.308

PubMed Abstract | Crossref Full Text | Google Scholar

36. Demir, R, Bodur, H, Tokoğlu, F, Olcay, I, Uçan, H, and Borman, P. High resolution computed tomography of the lungs in patients with rheumatoid arthritis. Rheumatol Int. (1999) 19:19–22. doi: 10.1007/s002960050093

Crossref Full Text | Google Scholar

37. Cortet, B, Perez, T, Roux, N, Flipo, RM, Duquesnoy, B, Delcambre, B, et al. Pulmonary function tests and high resolution computed tomography of the lungs in patients with rheumatoid arthritis. Ann Rheum Dis. (1997) 56:596–600. doi: 10.1136/ard.56.10.596

PubMed Abstract | Crossref Full Text | Google Scholar

38. Bilgici, A, Ulusoy, H, Kuru, O, Celenk, C, Unsal, M, and Danaci, M. Pulmonary involvement in rheumatoid arthritis. Rheumatol Int. (2005) 25:429–35. doi: 10.1007/s00296-004-0472-y

Crossref Full Text | Google Scholar

39. Attar, SM, Alamoudi, OS, and Aldabbag, AA. Prevalence and risk factors of asymptomatic bronchiectasis in patients with rheumatoid arthritis at a tertiary care center in Saudi Arabia. Ann Thorac Med. (2015) 10:176–80. doi: 10.4103/1817-1737.160836

PubMed Abstract | Crossref Full Text | Google Scholar

40. Wilczynska, MM, Condliffe, AM, and McKeon, DJ. Coexistence of bronchiectasis and rheumatoid arthritis: revisited. Respir Care. (2013) 58:694–701. doi: 10.4187/respcare.01857

PubMed Abstract | Crossref Full Text | Google Scholar

41. Puéchal, X, Génin, E, Bienvenu, T, Le Jeunne, C, and Dusser, DJ. Poor survival in rheumatoid arthritis associated with bronchiectasis: a family-based cohort study. PLoS One. (2014) 9:e110066. doi: 10.1371/journal.pone.0110066

PubMed Abstract | Crossref Full Text | Google Scholar

42. Geri, G, Dadoun, S, Bui, T, Del Castillo, PN, Paternotte, S, Dougados, M, et al. Risk of infections in bronchiectasis during disease-modifying treatment and biologics for rheumatic diseases. BMC Infect Dis. (2011) 11:304. doi: 10.1186/1471-2334-11-304

PubMed Abstract | Crossref Full Text | Google Scholar

43. Jenks, KA, Stamp, LK, O'Donnell, JL, Savage, RL, and Chapman, PT. Leflunomide-associated infections in rheumatoid arthritis. J Rheumatol. (2007) 34:2201–3.

PubMed Abstract | Google Scholar

44. Brode, SK, Jamieson, FB, Ng, R, Campitelli, MA, Kwong, JC, Paterson, JM, et al. Increased risk of mycobacterial infections associated with anti-rheumatic medications. Thorax. (2015) 70:677–82. doi: 10.1136/thoraxjnl-2014-206470

Crossref Full Text | Google Scholar

45. Puéchal, X, Bienvenu, T, Génin, E, Berthelot, JM, Sibilia, J, Gaudin, P, et al. Mutations of the cystic fibrosis gene in patients with bronchiectasis associated with rheumatoid arthritis. Ann Rheum Dis. (2011) 70:653–9. Epub 20101203. doi: 10.1136/ard.2010.142760

PubMed Abstract | Crossref Full Text | Google Scholar

46. Puéchal, X, Fajac, I, Bienvenu, T, Desmazes-Dufeu, N, Hubert, D, Kaplan, JC, et al. Increased frequency of cystic fibrosis deltaF508 mutation in bronchiectasis associated with rheumatoid arthritis. Eur Respir J. (1999) 13:1281–7. doi: 10.1183/09031936.99.13612889

PubMed Abstract | Crossref Full Text | Google Scholar

47. Hillarby, MC, McMahon, MJ, Grennan, DM, Cooper, RG, Clarkson, RW, Davies, EJ, et al. HLA associations in subjects with rheumatoid arthritis and bronchiectasis but not with other pulmonary complications of rheumatoid disease. Br J Rheumatol. (1993) 32:794–7. doi: 10.1093/rheumatology/32.9.794

Crossref Full Text | Google Scholar

48. Toussirot, E, Despaux, J, and Wendling, D. Increased frequency of HLA-DRB1*0401 in patients with RA and bronchiectasis. Ann Rheum Dis. (2000) 59:1001b–11001b. doi: 10.1136/ard.59.12.1001b

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: rheumatoid arthritis, bronchiectasis, Mendelian randomization, causal association, GWAS

Citation: Li Y, Wang W, Zhou D and Li L (2024) The causal relationship between rheumatoid arthritis and bronchiectasis: a bidirectional Mendelian randomization study. Front. Med. 11:1403851. doi: 10.3389/fmed.2024.1403851

Received: 20 March 2024; Accepted: 11 June 2024;
Published: 20 June 2024.

Edited by:

Cong-Qiu Chu, Oregon Health and Science University, United States

Reviewed by:

Carlos Montilla, University Hospital of Salamanca, Spain
Maria Leandro, University College London Hospitals NHS Foundation Trust, United Kingdom

Copyright © 2024 Li, Wang, Zhou and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yuanyuan Li, 18207190110@163.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.