Abstract
Background/purpose: Previous epidemiological studies have associated interstitial lung disease (ILD) with rheumatoid arthritis (RA), yet the causality of this relationship remains uncertain. This study aimed to investigate the genetic causal link between ILD and RA.
Methods: Genome-wide association study (GWAS) statistics for ILD and RA were collected from public datasets. Relevant single-nucleotide polymorphisms (SNPs) were selected by executing quality control steps from the GWAS summary results. A two-sample bidirectional Mendelian randomization (MR) analysis was performed to assess the causal relationship between the two conditions. The MR analysis primarily used the inverse variance weighting (IVW), weighted median (WM), and MR-Egger regression methods. Sensitivity analyses, including MR-Egger, leave-one-out, and MR Pleiotropy RESidual Sum and Outlier (MR-PRESSO), were conducted to evaluate the heterogeneity and pleiotropy. Replication analyses using Asian datasets were also conducted to enhance the robustness of our findings.
Results: In the European population, RA was found to increase the risk of ILD by 9.6% (OR: 1.096, 95% CI: 1.023–1.174, p = 0.009). Conversely, ILD was associated with a 12.8% increased risk of RA (OR: 1.128, 95% CI: 1.013–1.256, p = 0.029). Replication analyses from Asian GWAS further supported these findings, particularly the increased risk of ILD attributable to RA (OR: 1.33, 95% CI: 1.18–1.49, p-value <0.001).
Conclusion: Our findings underscore the clinical importance of screening for ILD in RA patients and suggest that effective management of RA could significantly benefit ILD patients. The potential applicability of novel RA treatments to ILD warrants further exploration. Additionally, racial disparities in the manifestation of these diseases should not be overlooked, as they may offer new perspectives for targeted therapies in diverse populations.
1 Introduction
Rheumatoid arthritis (RA) is a chronic, systemic autoimmune disease whose main characteristic is persistent joint inflammation that can lead to lifelong damage and loss of function (). It affects approximately 1% of the world’s population and 0.32%–0.36% in China (). The exact etiology of RA is unknown, but it has been established that its chronic course could be correlated with genetic factors, and the heritability of RA was estimated to be 60% (Yarwood et al., 2016). So far, using genome-wide association studies (GWAS), research studies have revealed that several genetic susceptibility loci play a vital role in RA, such as the human leukocyte antigen D-related B1 gene (HLA-DRB1) and protein tyrosine phosphatase non-receptor type 22 (PTPN22) (Tanaka, 2020).
Interstitial lung disease (ILD) is a large heterogeneous group of diseases affecting the lung parenchyma through inflammation and fibrosis, characterized by symptoms of cough, shortness of breath, and hypoxia (Newton et al., 2016; Rivera-Ortega and Molina-Molina, 2019). Epidemiological surveys estimate the prevalence of IDL to be as high as 24.9% (). Genetic factors are considered important risk factors for ILD (; Rivera-Ortega and Molina-Molina, 2019). Currently, many alleles strongly associated with ILD risk have been identified by GWAS, such as MUC5B (Newton et al., 2018), HLA-DRB1 (), and the gene encoding surface active protein C (SFTPC) ().
The association between RA and ILD has been extensively discussed. In some studies, ILD was identified as a severe extra-articular manifestation of RA (; ). A prospective registry study noted that patients with a combination of rheumatoid factor (RF) or anti-citrullinated protein antibody (ACPA) seropositivity had a higher prevalence of ILD when compared to seronegative subjects (Natalini et al., 2021; ). In addition, studies have also shown that significant lung inflammation may be associated with higher local and systemic ACPA concentrations (Natalini et al., 2021). Gregory C McDermott et al. reported that a substantial proportion of patients usually develop ILD prior to the manifestation of arthritis, suggesting a role for the lungs in the progression of RA disease (). Similarly, studies have shown that 10%–17% of patients were diagnosed with ILD before they were diagnosed with RA (). However, the causal relationship between the two conditions still needs to be fully understood. Therefore, it is of great interest to establish their association with more basic evidence and elucidate the mechanisms underlying the association between RA and ILD.
Mendelian randomization (MR) is an epidemiological method that can infer causality between the exposure and outcome by using single-nucleotide polymorphisms (SNPs) as an instrumental variable (IV) (). This method is not prone to reverse causation since disease states usually do not change the germline DNA sequences. Moreover, MR can limit confounding since genotypes are randomly assorted at meiosis (). Though well-designed randomized controlled trials (RCTs) are usually the best approach to estimating a causal relationship between a risk factor and a disease, their implementation is limited by the small sample size, limited external validity, short duration of an intervention, and ethical concerns that limit the implementation of RCTs. Thus, MR has specific advantages that can complement another study (). This study took advantage of MR and investigated the causal relationship between RA and ILD using genetic determinants estimated from GWAS summary statistics of different populations.
2 Materials and methods
2.1 Study design
MR analysis was performed to investigate the causal relationship between RA and ILD. Relevant SNPs were selected as IVs by executing quality control steps from the GWAS summary results, which need to meet three basic assumptions (Figure 1): (i) SNPs were strongly associated with exposure (correlation hypothesis); (ii) SNPs were not associated with confounding factors, meaning that the results were not affected by confounding factors (independence assumption); and (iii) SNPs were not associated with the outcome or did not directly affect the outcome (exclusivity hypothesis). Based on the previous research, we have listed the STROBE-MR checklist of MR studies to ensure the integrity of our processes, as shown in Supplementary Table S1 (Skrivankova et al., 2021a; Skrivankova et al., 2021b).
FIGURE 1
2.2 Data sources
2.2.1 Genetic datasets associated with rheumatoid arthritis
The RA-related outcome dataset of Europeans, derived from a summary of the GWAS previously communicated by Okada et al., included 14,361 cases and 43,923 healthy controls of European descent, totaling 8,747,963 RA-related SNPs (Okada et al., 2014). Genetic datasets associated with RA in the East Asian population are derived from an OpenGWAS project previously compiled by Sakaue Scott et al. (2013), containing 1,046 cases and 176,974 healthy controls, and the number of SNPs is 12,454,608. All RA cases met the 1987 RA diagnostic criteria of the American College of Rheumatology () or were diagnosed with RA by a specialist rheumatologist.
2.2.2 Genetic datasets for interstitial lung disease
Genetic prediction of the ILD-related gene GWAS data was obtained from the available meta-analysis of published GWASs. The European dataset contains 1,969 cases and 196,986 healthy controls, covering 16,380,417 genotyped SNPs. Genetic datasets for ILD contain 806 cases and 211,647 healthy controls from Ishigaki K et al., covering 8,885,805 SNPs from East Asia (Sakaue et al., 2021). All cases were diagnosed according to the American Thoracic Society and the European Respiratory Society guidelines.
Since this MR study was conducted using publicly available GWAS summary data, ethical approval and informed consent given all subjects could be found in the original publications.
2.3 Selection of instrumental variables
Strict filtering was performed to control SNP quality before MR analysis. First, we extracted SNPs associated with the genome-wide significance threshold of exposure (p < 1 × 10−5). Second, for the condition of independence from exposure, we selected SNPs with independent inheritance (r2 < 0.001) and no linkage disequilibrium (LD, clump distance >10,000 kb) in the summary statistics. Then, to avoid possible statistical bias from the original GWAS, we excluded IVs with minor allele frequency (MAF) less than the threshold of 0.01. Finally, the strongly correlated IVs were obtained by the F statistic >10 (; Pierce and Burgess, 2013).
Furthermore, it is confirmed that BMI, smoking, and alcohol consumption were implemented with a varying degree of effects on RA and ILD in previous studies (; ; ; ; Poudel et al., 2020). Based on this, the PhenoScanner (http://phenoscanner.medschl.cam.ac.uk) search was applied to remove these variants, which are associated with confounding factors. At the same time, the SNPs selected from the exposure and outcome datasets were harmonized, removing the palindromic SNPs using intermediate allele frequencies. Finally, the remaining IVs were used as tools for further MR analysis.
2.4 Mendelian randomization analyses
The preliminary analysis of this study used inverse variance weighting (IVW), weighted median (WM), and MR-Egger regression to estimate the causal relationship between the exposure and outcome. The IVW approach calculated the Wald ratio for each SNP on the outcome and obtained a pooled causal estimate, and this method allows for overdispersion. WM estimates were used to complement the IVW approach for reliable estimation and accounted for at least 50% of the analysis weight provided by the validation IV. The MR-Egger method was used to test for bias in the results of analyses that result from influencing the results through pathways other than exposure. Of these methods, IVW was the most crucial method for estimating causality between the exposure and outcome, which provides an accurate estimate of the causal effect of the outcome risk. These MR methods are detailed in published studies (; ).
2.5 Sensitivity analysis
To assess the robustness of these results and prevent potential pleiotropy and heterogeneity, a series of sensitivity analyses were conducted. Cochran’s Q test was applied to quantify the heterogeneity between SNPs and was visualized with funnel plots. The intercept term of MR-Egger regression and the MR-PRESSO global test was used to assess whether horizontal pleiotropy affected the results of MR analysis (). In addition, leave-one-out analysis was performed to determine whether any single SNP drove the causal estimates. Results were examined as statistically significant at p < 0.05. All analyses mentioned above were implemented in R V4.1.0 using the “TwoSampleMR” package ().
3 Results
3.1 Causal effects of RA on ILD
We incorporated 104 SNPs with a p-value less than 1 × 10−5 as IV SNPs. All SNPs were consistent with the F statistic >10. Twelve SNPs related to RA were removed for being palindromic with intermediate allele frequencies, and the potential confounding factors found in the PhenoScanner including rs10774624 (associated with smoking), rs34695944 (associated with alcohol consumption), and rs4272, rs34431565, rs773125, rs10857135, rs2736337, and rs9275183 (associated with body mass index, BMI) were removed. Then, the remaining 84 SNPs were used as IVs for subsequent MR analysis. Detailed information of IVs for RA is presented in Supplementary Table S2.
When setting ILD as the outcome, RA was causally associated with ILD, as shown in Table 1 and Figures 2, 3A. According to the IVW analysis, the presence of RA may increase 9.6% risk of ILD (OR: 1.096, 95% CI: 1.023–1.174, p = 0.009). Sensitivity analysis shows that there was heterogeneity among the IVs (Q = 109.17, PQ = 0.029, Figure 3B), while the MR-Egger intercept analysis suggested no evidence of horizontal pleiotropy (p = 0.772). Furthermore, leave-one-out plots suggested that the causal estimates were unlikely to be influenced by certain SNPs. This suggested that the findings were stable and reliable (Figure 3C).
TABLE 1
| Exposure | Outcome | SNP | MR method | OR | 95% CI | p-val | Q_p-val |
|---|---|---|---|---|---|---|---|
| RA | ILD | 84 | IVW | 1.096 | (1.023, 1.174) | 0.009* | 0.029* |
| MR-Egger | 1.079 | (0.947, 1.228) | 0.258 | 0.025* | |||
| WM | 1.097 | (0.979, 1.229) | 0.100 | - | |||
| ILD | RA | 8 | IVW | 1.128 | (1.013, 1.256) | 0.029* | 0.170 |
| MR-Egger | 0.661 | (0.445, 0.982) | 0.086 | 0.812 | |||
| WM | 1.144 | (1.012, 1.293) | 0.040* | - |
Causal effects of RA on ILD in European population.
Odds ratio (OR) values, confidence interval (CI), p-val and Q_p-val of MR results were obtained by IVW, MR-Egger, and weighted median in the group of European population. *p < 0.05.
FIGURE 2
FIGURE 3
3.2 Causal effects of ILD on RA
There were 25 SNPs that were significantly associated with ILD (p < 1 × 10−5) and were incorporated as IVs. Similarly, the F-statistic values of these IVs were all more than 10. The palindromic alleles were then removed and synergized with RA-associated SNPs, and the remaining eight SNPs were used as IVs for subsequent MR analysis. Detailed information of IVs for ILD is presented in Supplementary Table S3.
When setting RA as the outcome, ILD was also causally associated with RA (Table 1; Figures 2, 4A). The existence of ILD may increase the risk of RA by 12.8% (OR: 1.128, 95% CI: 1.013–1.256, PQ = 0.029). Results of the weighted median method also support our findings (OR: 1.144, 95% CI: 1.012–1.293, p = 0.04). In the analysis of sensitivity, Cochran’s Q test (Q = 2.97, p = 0.812, Figure 4B) indicated no heterogeneity in the causal effect between ILD and RA. No significant directional horizontal pleiotropy between RA and ILD was presented in the MR-Egger regression analysis. In the leave-one-out sensitivity analysis, no single SNP significantly biased the causal effect of ILD on RA (Figure 4C).
FIGURE 4
3.3 Replication analyses in East Asian population
3.3.1 Causal effects of RA on ILD
There were 46 SNPs that were significantly associated with ILD without linkage disequilibrium (r2 < 0.001, p < 1 × 10−5). The confounding factors such as rs10821944, rs111335405, rs2069235, rs2647192, rs4728142, rs7732397, and rs2618476 (all associated with the body mass index) found in the PhenoScanner were removed. Then, the palindromic alleles were removed and synergized with ILD-associated SNPs. The remaining 36 SNPs were used as IVs for subsequent MR analysis. Details of IVs are presented in Supplementary Table S4.
The MR estimates of different methods are presented in Table 2 and Figures 5, 6A. Notably, the results showed a strong causal relationship between RA and ILD (OR: 1.33, 95% CI: 1.18–1.49, p < 0.001). The result of weighted median methods further endorsed the causal effect (OR: 1.27, 95% CI: 1.11–1.46, p < 0.001). As for sensitivity analysis, we found that there was heterogeneity among the IVs (Q = 49.27, PQ = 0.043, Figure 6B). Furthermore, the MR-Egger intercept analysis shows that there was no evidence of horizontal pleiotropy (p = 0.144). A leave-one-out analysis was conducted to avoid horizontal pleiotropy caused by a single SNP (Figure 6C).
TABLE 2
| Exposure | Outcome | SNP | MR method | OR | 95% CI | p-val | Q_p-val |
|---|---|---|---|---|---|---|---|
| RA | ILD | 36 | IVW | 1.328 | (1.182, 1.491) | <0.001* | 0.029* |
| MR-Egger | 1.162 | (0.943, 1.432) | 0.167 | 0.044* | |||
| WM | 1.272 | (1.107, 1.462) | 0.001* | - | |||
| ILD | RA | 12 | IVW | 1.040 | (0.988, 1.096) | 0.134 | 0.457 |
| MR-Egger | 1.231 | (1.061, 1.429) | 0.021* | 0.873 | |||
| WM | 1.017 | (0.945, 1.095) | 0.646 | - |
Causal effects of RA on ILD in East Asian population.
Odds ratio (OR) values, confidence interval (CI), p-val and Q_p-val of MR results were obtained by IVW, MR-Egger, and weighted median in the group of East Asian population. *p < 0.05.
FIGURE 5
FIGURE 6
3.3.2 Causal effects of ILD on RA
There were 14 SNPs that were significantly associated with ILD without linkage disequilibrium (r2 < 0.001, p < 1 × 10−5). The confounder rs3129960 (associated with the body mass index) found in the PhenoScanner was removed. Then, the remaining 12 SNPs were used as IVs for subsequent MR analysis. Details of IVs are presented in Supplementary Table S5. Horizontal pleiotropy was unlikely to skew the causality of ILD with RA, according to the results of the MR-PRESSO global test (p = 0.542). However, results of IVW showed that an increase in the risk of having ILD was not statistically related to an increased risk of having RA (RA: OR = 1.040, 95% CI: 0.997–1.10, p = 0.13). The detailed information is presented in Table 2 and Figures 5, 7A). Furthermore, no significant evidence of horizontal pleiotropy was observed for IVs (Q = 5.27, PQ = 0.873) (Figure 7B), and leave-one-out plots suggested that the causal estimates were unlikely to be influenced by a single SNP (Figure 7C).
FIGURE 7
4 Discussion
Our study suggested that the inflammatory process of RA had a positive causal effect on the presence of ILD and vice versa in the European population. However, in the Asian population, we only found that the presence of RA caused an increased risk of ILD. This is the first study to explore the bidirectional causal relationship between RA and ILD through a two-sample MR approach based on GWAS summary statistics, providing evidence of a genetic link between RA and ILD. Our findings indicated that these two common autoimmune diseases may share similar underlying pathophysiological mechanisms, but racial differences cannot be ignored.
Many studies have found that RA increases the risk of developing ILD. For example, in a longitudinal study, 51% of patients were diagnosed with RA–ILD more than 5 years after diagnosing RA (Mohning et al., 2021). In an incident cohort of RA patients in the UK, 4% developed clinically significant RA–ILD on high-resolution CT (HR-CT) imaging during 15 years of follow-up (). In a smaller study of RA–ILD patients in China, ILD was diagnosed after RA in 69% of cases, and the median time between RA and RA–ILD diagnosis was 60 months (). However, not every study has reached the same conclusion about the association between RA and ILD. In contrast, some studies have concluded that the lungs may be a mucosal (and potentially an initiating) site for generating RA-related autoimmunity (; ). A study predominantly about middle-aged women of European ancestry showed that pulmonary involvement was present early in the disease course in RA (Wilsher et al., 2012). Moreover, in a survey conducted on 10%–20% of RA cases, respiratory symptoms may precede the onset of articular symptoms (). Similarly, studies of RA–ILD patients from Denmark, the United States, and China showed that 10%–17% of patients were diagnosed with ILD before the articular diagnosis of RA (). Due to their mechanistic complexity and proximity, it can be seen that there is a debate about the cause and the effect of the relationship between RA and ILD. Our bi-directional MR study synthesized the two-sided stance, further complementing previous studies that the two conditions may share similar underlying pathophysiological mechanisms.
There are several possible reasons why RA can increase the risk of developing ILD. Ayodeji Adegunsoye et al. generalized that under the regulation of some harmful environmental factors such as smoking, diet, viral infection, and ambient air pollutants, individuals with genetic variation are more susceptible to pulmonary fibrosis (PF) through epigenetic modifications such as DNA methylation, which is referred to as the gene–environment interaction. Among them, the variations of the MUC5B promoter are one of the most common genetic risk factors of early and established PF (). Similarly, a seminal study by Juge et al. in France showed that the MUC5B promoter variant is associated with an increased risk of developing ILD in patients with RA (). Building on this finding, a Finnish study, by combining large-scale genotype data with clinical data from a national healthcare registry, found that among RA patients, the lifetime risk of developing ILD was 16.8% for MUC5B carriers and 6.1% for MUC5B non-carriers (Palomaki et al., 2021). In addition, the role of MUC5B SNPs in the developing RA–ILD may be associated with age factor. It significantly increased the risk of RA–ILD early in the RA course (before or within 2 years of RA diagnosis) and after 55 years of age (). In particular, in a large real-life European multicenter idiopathic pulmonary fibrosis (IPF) study, van der Vis JJ et al. showed that MUC5B minor allele carriers were significantly older at diagnosis (p = 0.001), and among patients ≥56 years of age at diagnosis, the 3-year cumulative incidence of death was lower among MUC5B minor allele carriers (39%) than among non-carriers (57%) (van der Vis et al., 2023).
In addition, several studies suggested that ILD exacerbation was due to the inflammatory process of RA. For example, high titers of rheumatoid factor [RF] and anti-cyclic citrullinated peptide [CCP] antibodies are specifically associated with the development of ILD in patients with RA (; ; ; ), in which the HLA-DRB1*15 allele may confer a specific susceptibility to the production of anti-citrullinated fibrinogen antibodies (). More specifically, the combination of matrix metalloproteinase-7 (MMP-7), pulmonary and activation-regulated chemokine, and surfactant protein D enhances the above association (). Among them, MMP-7 was identified as a potential biomarker for RA–ILD and was significantly elevated in the serum of RA patients with clinical and subclinical ILD (). A prospective cohort study evaluated the association of plasma matrix metalloproteinases (MMPs) with the incidence of interstitial lung disease in patients with rheumatoid arthritis in a large multicenter RA cohort, further supporting the potential pathogenic role of MMP-7 and MMP-9 for RA–ILD (). The association of RA-induced pulmonary fibrosis also involves MHC gene loci on chromosome 6, such as HLA-B54, HLA-DQ1B*0601, HLA-B40, and sites encoding a-1 protease inhibitors, which are associated with increased ILD risk in patients with RA (Spagnolo et al., 2014).
Although the exact mechanism of how ILD increases the risk of RA is unknown, several possible explanations exist for the results of reverse MR analysis. According to the “mucosal origins hypothesis,” RA-related autoimmunity disorder is initiated at a mucosal site, including the lungs, and then transitions to involve the synovial joints (; ). In a genetically susceptible individual, injury to the alveoli, airway epithelium, and mucosa caused by smoking, microbial dysbiosis, or other inhalant exposures can lead to citrullination of protein, production of neutrophil extracellular traps, generation of local pulmonary mucosal RA-associated autoantibodies, and establishing of systemic autoimmunity ultimately (; Valesini et al., 2015; Scher et al., 2016; ; Prisco et al., 2020). Specifically, in genetically susceptible individuals, the formation rate of the citrullinated protein is high in lung tissues, and citrullination can trigger an immune response leading to the production of ACPA, resulting in an increased risk of RA in patients (; Paulin et al., 2015). Another possible explanation is that the elevated cytokines and chemokines caused by the chronic inflammatory state in lung mucosal sites also affect lymphocyte subsets such as T cells, B cells, and macrophages, which recirculate and populate other sites and, in the process, influence the immune status in other tissues participating in the autoimmune response (Rangel-Moreno et al., 2006; ; ). However, findings about ILD increasing the risk of RA were not so consistent in patients with different races in our study. McFarlane et al. mentioned in their study that conflicting reports exist regarding the racial distribution of RA–ILD (McFarlane et al., 2019). Furthermore, disease activity of RA remained significantly different across the Asian groups versus European groups in the research of . These factors may have a substantial effect on the results, so further studies are needed to identify the genetic differences.
This study has several strengths. First, it is the first bidirectional two-sample MR analysis to reveal a causal relationship between RA and ILD. Second, the screening of instrumental variables set stringent conditions and therefore had high statistical power, allowing for a more robust analysis of possible causal relationships. Finally, in contrast to randomized controlled trials, MR uses an index of genetic variation to measure the causality of disease-related risk factors, overcoming the bias caused by confounding or reverse causality inherent in observational studies.
At the same time, this study has some limitations. First, the study population included in the MR analysis was of European and East Asian origin, and it remains to be verified whether the results represent the whole population. In addition, the minority differences associated with race cannot be ignored in the two diseases, which is worthy of further analysis. Second, though we have taken some steps to exclude heterogeneity and horizontal pleiotropy, the influence on the results is undetermined, and further analysis and studies are still warranted to understand the underlying mechanisms. Finally, it was impossible to determine whether there were overlapping participants in the GWAS data for the exposures and outcomes involved in this study. Fortunately, using strong instruments in this study (e.g., an F statistic much greater than 10) should minimize the potential bias of sample overlap.
In conclusion, our study suggested a causal relationship between RA and ILD, but there are minority racial differences. Based on our findings, it is reasonable to consider promoting a routine screening of ILD in RA patients. Proper management of inflammation of RA is essential for downregulating the risk of ILD. This result may facilitate future insights into the biological relationship between RA and ILD.
Statements
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.
Author contributions
RZ: data curation and writing–original draft. Y-WZ: data curation, visualization, and writing–original draft. J-CG: writing–original draft and formal analysis. JQ: writing–review and editing. SS: resources and writing–review and editing. T-TZ: writing–review and editing and project administration. H-YZ: writing–review and editing and validation. S-XZ: resources, writing–review and editing, and project administration.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. We acknowledge the Natural Science Foundation of Shanxi Province (No. 202203021221269) and the National Natural Science Foundation of China (No. 82001740) for invaluable funding.
Acknowledgments
The authors would like to thank Okada et al., Sun et al., and the genome-wide association study consortia, who made their summary statistics publicly available for this study, without whom this effort would not be possible, and they would also like to thank the study staff and participants for their significant contributions. They are also grateful to the Bolyn (Liaoning) Biotechnology Research Group for their expert assistance in data processing.
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/fgene.2024.1395315/full#supplementary-material
References
1
AdegunsoyeA.KropskiJ. A.BehrJ.BlackwellT. S.CorteT. J.CottinV.et al (2024). Genetics and genomics of pulmonary fibrosis: charting the molecular landscape and shaping precision medicine. Med. Am. J. Respir. Crit. Care Med.10.1164/rccm.202401-0238SO
2
AkiyamaM.KanekoY. (2022). Pathogenesis, clinical features, and treatment strategy for rheumatoid arthritis-associated interstitial lung disease. Autoimmun. Rev.21, 103056. 10.1016/j.autrev.2022.103056
3
ArnettF. C.EdworthyS. M.BlochD. A.McShaneD. J.FriesJ. F.CooperN. S.et al (1988). The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum.31, 315–324. 10.1002/art.1780310302
4
BakaZ.BuzásE.NagyG. (2009). Rheumatoid arthritis and smoking: putting the pieces together. Arthritis Res. Ther.11, 238. 10.1186/ar2751
5
BongartzT.NanniniC.Medina-VelasquezY. F.AchenbachS. J.CrowsonC. S.RyuJ. H.et al (2010). Incidence and mortality of interstitial lung disease in rheumatoid arthritis: a population-based study. a population-based study Arthritis Rheum.62, 1583–1591. 10.1002/art.27405
6
BurgessS.ButterworthA.ThompsonS. G. (2013). Mendelian randomization analysis with multiple genetic variants using summarized data. Genet. Epidemiol.37, 658–665. 10.1002/gepi.21758
7
BurgessS.ThompsonS. G.CRP CHD Genetics Collaboration (2011). Avoiding bias from weak instruments in Mendelian randomization studies. Int. J. Epidemiol.40, 755–764. 10.1093/ije/dyr036
8
CatrinaA. I.YtterbergA. J.ReynisdottirG.MalmströmV.KlareskogL. (2014). Lungs, joints and immunity against citrullinated proteins in rheumatoid arthritis. Nat. Rev. Rheumatol.10, 645–653. 10.1038/nrrheum.2014.115
9
ChangK.YangS. M.KimS. H.HanK. H.ParkS. J.ShinJ. I. (2014). Smoking and rheumatoid arthritis. Int. J. Mol. Sci.15, 22279–22295. 10.3390/ijms151222279
10
ChenJ.DoyleT. J.LiuY.AggarwalR.WangX.ShiY.et al (2015b). Biomarkers of rheumatoid arthritis-associated interstitial lung disease. Arthritis Rheumatol.67, 28–38. 10.1002/art.38904
11
ChenL.HuangZ.YangB.CaiB.SuZ.WangL. (2015a). Association of E26 transformation specific sequence 1 variants with rheumatoid arthritis in Chinese han population. PLoS One10, e0134875. 10.1371/journal.pone.0134875
12
ChenR. X.ZhaoL. D.XiaoX. Y.SongL.DuH. Y.XuZ. J.et al (2021). Distinctive clinical characteristics and outcome of ILD-onset rheumatoid arthritis and ACPA-positive ILD: a longitudinal cohort of 282 cases. Clin. Rev. Allergy Immunol.60, 46–54. 10.1007/s12016-020-08819-0
13
CojocaruM.CojocaruI. M.SilosiI.VrabieC. D.TanasescuR. (2010). Extra-articular manifestations in rheumatoid arthritis. Arthritis Maedica (Bucur)5, 286–291.
14
Davey SmithG.HemaniG. (2014). Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum. Mol. Genet.23, R89–R98. 10.1093/hmg/ddu328
15
DemoruelleM. K.WeismanM. H.SimonianP. L.LynchD. A.SachsP. B.PedrazaI. F.et al (2012). Brief report: airways abnormalities and rheumatoid arthritis-related autoantibodies in subjects without arthritis: early injury or initiating site of autoimmunity?Arthritis Rheum.64, 1756–1761. 10.1002/art.34344
16
DidelezV.SheehanN. (2007). Mendelian randomization as an instrumental variable approach to causal inference. Methods Med. Res.16, 309–330. 10.1177/0962280206077743
17
DoyleT. J.PatelA. S.HatabuH.NishinoM.WuG.OsorioJ. C.et al (2015). Detection of rheumatoid arthritis-interstitial lung disease is enhanced by serum biomarkers. Biomarkers Am. J. Respir. Crit. Care Med.191, 1403–1412. 10.1164/rccm.201411-1950OC
18
DuarteA. C.PorterJ. C.LeandroM. J. (2019). The lung in a cohort of rheumatoid arthritis patients-an overview of different types of involvement and treatment. Rheumatol. Oxf.58, 2031–2038. 10.1093/rheumatology/kez177
19
EggerM.Davey SmithG.SchneiderM.MinderC. (1997). Bias in meta-analysis detected by a simple, graphical test. Bmj315, 629–634. 10.1136/bmj.315.7109.629
20
EnglandB. R.HershbergerD. (2020). Management issues in rheumatoid arthritis-associated interstitial lung disease. Curr. Opin. Rheumatol.32, 255–263. 10.1097/bor.0000000000000703
21
EspositoA. J.ChuS. G.MadanR.DoyleT. J.DellaripaP. F. (2019). Thoracic manifestations of rheumatoid arthritis. Chest Med.40, 545–560. 10.1016/j.ccm.2019.05.003
22
FazeliM. S.KhaychukV.WittstockK.HanX.CrocketG.LinM.et al (2021). Rheumatoid arthritis-associated interstitial lung disease: epidemiology, risk/prognostic factors, and treatment landscape. Clin. Exp. Rheumatol.39, 1108–1118. 10.55563/clinexprheumatol/h9tc57
23
FriedlanderH. M.FordJ. A.ZaccardelliA.TerrioA. V.ChoM. H.SparksJ. A. (2020). Obstructive lung diseases and risk of rheumatoid arthritis. arthritis Expert Rev. Clin. Immunol.16, 37–50. 10.1080/1744666x.2019.1698293
24
García-SanchoC.Buendía-RoldánI.Fernández-PlataM. R.NavarroC.Pérez-PadillaR.VargasM. H.et al (2011). Familial pulmonary fibrosis is the strongest risk factor for idiopathic pulmonary fibrosis. Respir. Med.105, 1902–1907. 10.1016/j.rmed.2011.08.022
25
GeorgeM. D.BakerJ. F. (2016). The obesity epidemic and consequences for rheumatoid arthritis care. Curr. Rheumatol. Rep.18, 6. 10.1007/s11926-015-0550-z
26
GilesJ. T.DanoffS. K.SokoloveJ.WagnerC. A.WinchesterR.PappasD. A.et al (2014). Association of fine specificity and repertoire expansion of anticitrullinated peptide antibodies with rheumatoid arthritis associated interstitial lung disease. Ann. Rheum. Dis.73, 1487–1494. 10.1136/annrheumdis-2012-203160
27
GizinskiA. M.MascoloM.LoucksJ. L.KervitskyA.MeehanR. T.BrownK. K.et al (2009). Rheumatoid arthritis (RA)-specific autoantibodies in patients with interstitial lung disease and absence of clinically apparent articular RA. Clin. Rheumatol.28, 611–613. 10.1007/s10067-009-1128-9
28
GreenbergJ. D.SpruillT. M.ShanY.ReedG.KremerJ. M.PotterJ.et al (2013). Racial and ethnic disparities in disease activity in patients with rheumatoid arthritis. arthritis Am. J. Med.126, 1089–1098. 10.1016/j.amjmed.2013.09.002
29
GyetvaiA.SzekaneczZ.SoosL.SzaboZ.FeketeA.KapitanyA.et al (2010). New classification of the shared epitope in rheumatoid arthritis: impact on the production of various anti-citrullinated protein antibodies. Rheumatol. Oxf.49, 25–33. 10.1093/rheumatology/kep338
30
HemaniG.ZhengJ.ElsworthB.WadeK. H.HaberlandV.BairdD.et al (2018). The MR-Base platform supports systematic causal inference across the human phenome. Elife7, e34408. 10.7554/eLife.34408
31
HolersV. M.DemoruelleM. K.KuhnK. A.BucknerJ. H.RobinsonW. H.OkamotoY.et al (2018). Rheumatoid arthritis and the mucosal origins hypothesis: protection turns to destruction. Nat. Rev. Rheumatol.14, 542–557. 10.1038/s41584-018-0070-0
32
HolersV. M.KuhnK. A.DemoruelleM. K.NorrisJ. M.FiresteinG. S.JamesE. A.et al (2022). Mechanism-driven strategies for prevention of rheumatoid arthritis. arthritis Rheumatol. Autoimmun.2, 109–119. 10.1002/rai2.12043
33
HyldgaardC.HilbergO.PedersenA. B.UlrichsenS. P.LøkkeA.BendstrupE.et al (2017). A population-based cohort study of rheumatoid arthritis-associated interstitial lung disease: comorbidity and mortality. Ann. Rheum. Dis.76, 1700–1706. 10.1136/annrheumdis-2017-211138
34
InuiN.EnomotoN.SudaT.KageyamaY.WatanabeH.ChidaK. (2008). Anti-cyclic citrullinated peptide antibodies in lung diseases associated with rheumatoid arthritis. Clin. Biochem.41, 1074–1077. 10.1016/j.clinbiochem.2008.06.014
35
JugeP. A.LeeJ. S.EbsteinE.FurukawaH.DobrinskikhE.GazalS.et al (2018). MUC5B promoter variant and rheumatoid arthritis with interstitial lung disease. N. Engl. J. Med.379, 2209–2219. 10.1056/NEJMoa1801562
36
KaduraS.RaghuG. (2021). Rheumatoid arthritis-interstitial lung disease: manifestations and current concepts in pathogenesis and management. Eur. Respir. Rev.30, 210011. 10.1183/16000617.0011-2021
37
KellyC. A.SaravananV.NisarM.ArthanariS.WoodheadF. A.Price-ForbesA. N.et al (2014). Rheumatoid arthritis-related interstitial lung disease: associations, prognostic factors and physiological and radiological characteristics--a large multicentre UK study. Rheumatol. Oxf.53, 1676–1682. 10.1093/rheumatology/keu165
38
KropskiJ. A. (2020). Familial interstitial lung disease. Semin. Respir. Crit. Care Med.41, 229–237. 10.1055/s-0040-1708054
39
LiX.MengX.TimofeevaM.TzoulakiI.TsilidisK. K.IoannidisJ. P.et al (2017). Serum uric acid levels and multiple health outcomes: umbrella review of evidence from observational studies, randomised controlled trials, and Mendelian randomisation studies. Bmj357, j2376. 10.1136/bmj.j2376
40
LucchinoB.SpinelliF. R.IannuccelliC.GuzzoM. P.ContiF.Di FrancoM. (2019). Mucosa-environment interactions in the pathogenesis of rheumatoid arthritis. Arthritis cells.8, 700. 10.3390/cells8070700
41
LueddersB. A.WheelerA. M.AschermanD. P.BakerJ. F.DuryeeM. J.YangY.et al (2024). Plasma matrix metalloproteinase concentrations and risk of interstitial lung disease in a prospective rheumatoid arthritis cohort. Arthritis Rheumatol.10.1002/art.42812
42
ManousakiD.HarroudA.MitchellR. E.RossS.ForgettaV.TimpsonN. J.et al (2021). Vitamin D levels and risk of type 1 diabetes: a Mendelian randomization study. PLoS Med.18, e1003536. 10.1371/journal.pmed.1003536
43
MaxwellJ. R.GowersI. R.MooreD. J.WilsonA. G. (2010). Alcohol consumption is inversely associated with risk and severity of rheumatoid arthritis. arthritis Rheumatol. Oxf.49, 2140–2146. 10.1093/rheumatology/keq202
44
McDermottG.GillR.GagneS.ByrneS.HuangW.CuiJ.et al (2022). Associations of the MUC5B promoter variant with timing of interstitial lung disease and rheumatoid arthritis onset. Rheumatol. Oxf.61, 4915–4923. 10.1093/rheumatology/keac152
45
McDermottG. C.DoyleT. J.SparksJ. A. (2021). Interstitial lung disease throughout the rheumatoid arthritis disease course. Curr. Opin. Rheumatol.33, 284–291. 10.1097/bor.0000000000000787
46
McFarlaneI. M.ZhazS. Y.BhamraM. S.BurzaA.KollaS.AlvarezM. R.et al (2019). Assessment of interstitial lung disease among black rheumatoid arthritis patients. Clin. Rheumatol.38, 3413–3424. 10.1007/s10067-019-04760-6
47
MohningM. P.AmiguesI.DemoruelleM. K.Fernández PérezE. R.HuieT. J.KeithR. K.et al (2021). Duration of rheumatoid arthritis and the risk of developing interstitial lung disease. ERJ Open Res.7, 00633. 10.1183/23120541.00633-2020
48
NataliniJ. G.BakerJ. F.SinghN.MahajanT. D.RoulP.ThieleG. M.et al (2021). Autoantibody seropositivity and risk for interstitial lung disease in a prospective male-predominant rheumatoid arthritis cohort of U.S. Veterans. Ann. Am. Thorac. Soc.18, 598–605. 10.1513/AnnalsATS.202006-590OC
49
NewtonC. A.BatraK.TorrealbaJ.KozlitinaJ.GlazerC. S.AravenaC.et al (2016). Telomere-related lung fibrosis is diagnostically heterogeneous but uniformly progressive. Eur. Respir. J.48, 1710–1720. 10.1183/13993003.00308-2016
50
NewtonC. A.MolyneauxP. L.OldhamJ. M. (2018). Clinical genetics in interstitial lung disease. Front. Med. (Lausanne)5, 116. 10.3389/fmed.2018.00116
51
OkadaY.WuD.TrynkaG.RajT.TeraoC.IkariK.et al (2014). Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature506, 376–381. 10.1038/nature12873
52
PalomakiA.PalotieA.KoskelaJ.EklundK. K.PirinenM.et al (2021). Lifetime risk of rheumatoid arthritis-associated interstitial lung disease in MUC5B mutation carriers. Ann. Rheum. Dis.80, 1530–1536. 10.1136/annrheumdis-2021-220698
53
PaulinF.DoyleT. J.FletcherE. A.AschermanD. P.RosasI. O. (2015). Rheumatoid arthritis-associated interstitial lung disease and idiopathic pulmonary fibrosis: shared mechanistic and phenotypic traits suggest overlapping disease mechanisms. Rev. Invest. Clin.67, 280–286.
54
PierceB. L.BurgessS. (2013). Efficient design for Mendelian randomization studies: subsample and 2-sample instrumental variable estimators. Am. J. Epidemiol.178, 1177–1184. 10.1093/aje/kwt084
55
PoudelD.GeorgeM. D.BakerJ. F. (2020). The impact of obesity on disease activity and treatment response in rheumatoid arthritis. Arthritis Curr. Rheumatol. Rep.22, 56. 10.1007/s11926-020-00933-4
56
PriscoL. C.MartinL. W.SparksJ. A. (2020). Inhalants other than personal cigarette smoking and risk for developing rheumatoid arthritis. Curr. Opin. Rheumatol.32, 279–288. 10.1097/bor.0000000000000705
57
Rangel-MorenoJ.HartsonL.NavarroC.GaxiolaM.SelmanM.RandallT. D. (2006). Inducible bronchus-associated lymphoid tissue (iBALT) in patients with pulmonary complications of rheumatoid arthritis. J. Clin. Invest.116, 3183–3194. 10.1172/jci28756
58
Rivera-OrtegaP.Molina-MolinaM. (2019). Interstitial lung diseases in developing countries. Ann. Glob. Health85, 4. 10.5334/aogh.2414
59
SakaueS.KanaiM.TanigawaY.KarjalainenJ.KurkiM.KoshibaS.et al (2021). A cross-population atlas of genetic associations for 220 human phenotypes. Nat. Genet.53, 1415–1424. 10.1038/s41588-021-00931-x
60
ScherJ. U.JoshuaV.ArtachoA.Abdollahi-RoodsazS.ÖckingerJ.KullbergS.et al (2016). The lung microbiota in early rheumatoid arthritis and autoimmunity. Microbiome4, 60. 10.1186/s40168-016-0206-x
61
ScottI. C.TanR.StahlD.SteerS.LewisC. M.CopeA. P. (2013). The protective effect of alcohol on developing rheumatoid arthritis: a systematic review and meta-analysis. Rheumatol. Oxf.52, 856–867. 10.1093/rheumatology/kes376
62
SkrivankovaV. W.RichmondR. C.WoolfB. A. R.DaviesN. M.SwansonS. A.VanderWeeleT. J.et al (2021a). Strengthening the reporting of observational studies in epidemiology using mendelian randomisation (STROBE-MR): explanation and elaboration. BMJ375, n2233. 10.1136/bmj.n2233
63
SkrivankovaV. W.RichmondR. C.WoolfB. A. R.YarmolinskyJ.DaviesN. M.SwansonS. A.et al (2021b). Strengthening the reporting of observational studies in epidemiology using mendelian randomization: the STROBE-MR statement. JAMA326, 1614–1621. 10.1001/jama.2021.18236
64
SpagnoloP.GrunewaldJ.du BoisR. M. (2014). Genetic determinants of pulmonary fibrosis: evolving concepts. Lancet Respir. Med.2, 416–428. 10.1016/S2213-2600(14)70047-5
65
TanakaY. (2020). Rheumatoid arthritis. Rheum. arthritis Inflamm. Regen.40, 20. 10.1186/s41232-020-00133-8
66
ValesiniG.GerardiM. C.IannuccelliC.PacucciV. A.PendolinoM.ShoenfeldY. (2015). Citrullination and autoimmunity. Autoimmun. Rev.14, 490–497. 10.1016/j.autrev.2015.01.013
67
van der VisJ. J.PrasseA.RenzoniE. A.StockC. J. W.CaliskanC.MaherT. M.et al (2023). MUC5B rs35705950 minor allele associates with older age and better survival in idiopathic pulmonary fibrosis. Respirology28, 455–464. 10.1111/resp.14440
68
WilsherM.VoightL.MilneD.TehM.GoodN.KolbeJ.et al (2012). Prevalence of airway and parenchymal abnormalities in newly diagnosed rheumatoid arthritis. Respir. Med.106, 1441–1446. 10.1016/j.rmed.2012.06.020
69
YarwoodA.HuizingaT. W.WorthingtonJ. (2016). The genetics of rheumatoid arthritis: risk and protection in different stages of the evolution of RA. Rheumatol. Oxf.55, 199–209. 10.1093/rheumatology/keu323
Summary
Keywords
rheumatoid arthritis, interstitial lung disease, Mendelian randomization, causal relationship, gene
Citation
Zhao R, Zhang Y-W, Guo J-C, Qiao J, Song S, Zhang T-T, Zhang H-Y and Zhang S-X (2024) Genetic evidence reveals a causal relationship between rheumatoid arthritis and interstitial lung disease. Front. Genet. 15:1395315. doi: 10.3389/fgene.2024.1395315
Received
03 March 2024
Accepted
26 April 2024
Published
14 May 2024
Volume
15 - 2024
Edited by
Ali H. Hajeer, King Abdulaziz Medical City, Saudi Arabia
Reviewed by
Theodoros Karampitsakos, University of South Florida, United States
Jianan Zhao, Shanghai Guanghua Rheumatology Hospital, China
Updates
Copyright
© 2024 Zhao, Zhang, Guo, Qiao, Song, Zhang, Zhang and Zhang.
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: Sheng-Xiao Zhang, zhangshengxiao1@sxmu.edu.cn
† These authors have contributed equally to this work and share first authorship
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.