- Open Access
Factors associated with mortality in rheumatoid arthritis-associated interstitial lung disease: a systematic review and meta-analysis
Respiratory Research volume 22, Article number: 264 (2021)
Interstitial lung disease (ILD) is a common and potentially life-threatening complication for rheumatoid arthritis (RA) patients. However, there is a lack of clear prognostic factors in rheumatoid arthritis-associated interstitial lung disease (RA-ILD) patients. The purpose of this study was to complete a systematic review and meta-analysis of the factors associated with mortality in RA-ILD patients.
Medline, EMBASE and the Cochrane Library were searched up to September 1, 2020. The Newcastle–Ottawa Scale (NOS) was applied to assess the methodological quality of the eligible studies. Study characteristics and magnitude of effect sizes were extracted. Then, pooled hazard ratios (HRs) with the corresponding 95% confidence intervals (CIs) and pooled risk ratios (RRs) with 95% CIs were calculated to assess the factors associated with mortality in RA-ILD.
Twenty-three of 3463 articles were eligible, and ten factors associated with mortality for RA-ILD were evaluated in the meta-analysis. Older age (HRs = 1.04, 95% CI 1.03–1.05), male sex (HRs = 1.44, 95% CI 1.21–1.73), having a smoking history (HRs = 1.42, 95% CI 1.03–1.96), lower diffusing capacity of the lung for carbon monoxide (DLCO)% predicted (HRs = 0.98, 95% CI 0.97–1.00), forced vital capacity (FVC)% predicted (HRs = 0.99, 95% CI 0.98–1.00), composite physiological index (CPI) (HRs = 1.04, 95% CI 1.02–1.06), usual interstitial pneumonia (UIP) pattern on HRCT (HRs = 1.88, 95% CI 1.14–3.10 and RRs = 1.90, 95% CI 1.50–2.39), emphysema presence (HRs = 2.31, 95% CI 1.58–3.39), and acute exacerbation of ILD (HRs = 2.70, 95% CI 1.67–4.36) were associated with increased mortality in RA-ILD, whereas rheumatoid factor (RF) positive status was not associated.
Through this systematic review and meta-analysis, we found that older age, male sex, smoking history, higher CPI, lower DLCO% predicted, lower FVC% predicted, UIP pattern on HRCT, emphysema presence and acute exacerbation of ILD were associated with an increased risk of mortality in RA-ILD.
Rheumatoid arthritis (RA) is an autoimmune, inflammatory arthritis with extra-articular manifestations in up to 50% of patients . Among those manifestations, clinical and subclinical interstitial lung disease (ILD) are the most frequent extra-articular manifestations with increased morbidity and mortality, and it can be detected in 10–60% of RA patients [2, 3]. Rheumatoid arthritis-associated interstitial lung disease (RA-ILD) has a heterogeneous clinical presentation and disease course. Despite an increasing recognition of RA-ILD and growing knowledge of the related pathogenesis, ILD remains a major cause of mortality with a median survival of only 3–7 years after diagnosis [4,5,6]. Better means of identifying and treating RA-ILD are necessary to improve the prognosis. Multiple studies have previously investigated prognostic factors of RA-ILD. However, these findings were derived from retrospective studies with small sample sizes and the importance of their findings was unclear in isolation.
Identifying the factors consistently associated with mortality in RA-ILD would improve prognostication and provide for making better decision for clinicians in the management of these patients. Therefore, the aim of this systematic review and meta-analysis was to overcome the limitations of a primary study and identify the factors associated with all-cause mortality in RA-ILD.
We followed the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines throughout this review . We searched the Medline, EMBASE and Cochrane databases, and the Cochrane Library [the Cochrane Database of Systematic Reviews and the Cochrane Central Register of Controlled Trials (CENTRAL)] through September 1, 2020. Databases were searched and data abstracted by two authors working independently. We used subject headings and text words related to the study population to finish the search. The reference lists of eligible studies and relevant review articles were also hand-searched to find additional reports. The detailed search strategy is provided in Additional file 3: Table 1 in the online additional material.
Two authors independently evaluated each study for eligibility, sequentially reviewing the title, abstract and full text of each publication. Studies were included if they identified predictors of mortality in adults with RA-ILD, were published in English and enrolled at least 10 patients with RA-ILD. The primary outcome was all-cause mortality. Case reports, editorials, letters, review articles and conference proceedings were not considered. If RA-ILD patients overlapped between two studies, the study with the larger sample size was prioritized in the analysis. Any uncertainties or disagreements between two authors were resolved by discussion and consensus.
Data extraction and quality assessment
Two authors independently extracted data from the included articles based on a predefined data extraction form. Extracted data included the first author’s name, year of publication, study location, study design, sample size, demographic features of the subjects, outcomes, potential prognostic factors and their effect estimates. The methodological quality of each cohort or case–control study was assessed using the 9-point Newcastle–Ottawa scale (NOS). Studies were of high quality if the score was ≥ 7. Two investigators assessed the quality of the studies through consultations to reach consensus.
All of the statistical analyses were performed using STATA 13.0 software (Stata Corp. LP, College Station, TX, USA). Predictors of mortality in RA-ILD in more than one of the selected studies were quantified by calculating pooled hazard ratios (HRs), pooled risk ratios (RRs), and 95% confidence intervals (CIs) using Mantel–Haenszel fixed-effects or DerSimonian–Laird random-effects models. Only unadjusted effect estimates of potential prognostic factors were combined and the effect estimates derived from multivariate models were described qualitatively. The results from the fixed effects model were presented only when there was no heterogeneity between studies; otherwise, the results from the random-effects model were presented. Heterogeneity between studies was assessed using I2 statistics, and statistical significance was considered with a P value of < 0.05 . Forest plots were used to display the results from the individual studies and the pooled estimates. The potential for publication bias was evaluated by funnel plots and Egger’s test if 5 or more studies were available for meta-analysis . If publication bias was present, the trim and fill method was used to verify the asymmetric funnel plot . If combing data were deemed inappropriate, the results were reported qualitatively (due to a small number of studies or substantial clinical or methodological diversity). P-values < 0.05 were considered statistically significant.
A total of 3463 reports were identified through Medline, EMBASE, and the Cochrane Library. After excluding 471 duplicates, 264 non-English records, 1895 reports of ineligible types (consisting of 612 conference proceedings, 758 review articles, 432 case reports and 93 editorials or letters) and 715 irrelevant articles, the remaining 118 reports were screened as full-texts. Of these reports, 87 were excluded due to a lack of data on RA-ILD in 11 articles, no data on mortality due to RA-ILD in 32 articles, and no data on predictors of mortality in 44 articles. Thus, eight articles using the same study cohort were excluded [11,12,13,14,15,16,17,18]. Finally, 23 studies were eligible for this systematic review and meta-analysis and 16 studies were included in the quantitative synthesis [6, 19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40] (Fig. 1).
Study and subject characteristics
The characteristics of the included studies are shown in Table 1. The 23 studies included two prospective cohort studies [19, 22], 13 retrospective cohort studies [6, 23, 24, 27, 28, 30, 31, 33,34,35,36, 38, 40], and eight case–control studies [20, 21, 24, 26, 29, 32, 37, 39]. The studies were from nine different countries (US, UK, Japan, Republic of Korea, Mexico, Italy, Finland, Denmark and China). Of them, the majority of studies took place in the US (n = 6) and Japan (n = 6), followed by the Republic of Korea (n = 5), the UK (n = 4), China (n = 2), Italy (n = 1), Mexico (n = 1), Finland (n = 1), and Denmark (n = 1). The majority (62.5%) were single-centre studies. The 1987 American College of Rheumatology (ACR) diagnostic criteria of RA were used in 15 studies, and the 2010 ACR/EULAR diagnostic criteria of RA were used in eight studies [41, 42]. Variable criteria were used for the diagnosis of ILD, including findings on high-resolution computed tomography (HRCT) only; surgical lung biopsy only; HRCT and surgical lung biopsy; and a combination of clinical presentation, pulmonary function tests, and HRCT findings. The methodological quality of the included study is assessed by the NOS in Additional file 4: Table 2. The studies were generally of moderate-high quality. There were a total of 3062 RA-ILD patients in the included studies, with 932 deaths. The characteristics of the included study subjects are shown in Additional file 5: Table 3.
The patient-specific variables included age, sex, and smoking history in the meta-analysis. Fourteen studies evaluated the relationship between older age and mortality in RA-ILD. In the meta-analysis, older age was associated with increased mortality in the 11 studies reporting HRs (HRs = 1.04, 95% CI 1.03–1.04, P = 0.0001) [6, 22, 24,25,26, 28, 31, 35,36,37] (Fig. 2A). Heterogeneity across the studies was low. Visual inspection of the funnel plot showed that the funnel plot was symmetrical, suggesting that there was no publication bias (Additional file 2: Fig. 2A). In addition, the P-value of Egger’s test was 0.369 indicating that there was no evidence of publication bias. Two studies found no relationship between age and mortality in RA-ILD with odds ratio (OR) [19, 33]. However, the OR could not be pooled due to a lack of the same age threshold. Fourteen studies evaluated the impact of sex on RA-ILD mortality. There was an increased risk of mortality among men (HRs = 1.44, 95% CI 1.21–1.73, P = 0.001) [6, 22, 24,25,26,27, 30, 31, 35,36,37, 39] (Fig. 2B). Heterogeneity across the studies was low. The funnel plot was symmetrical and the P-value of Egger’s test was 0.093 (Additional file 2: Fig. 2B). There was no evidence of publication bias. Moreover, one study showed that female sex was associated with better survival (HR = 0.3, P = 0.008) in a multivariate model after adjusting for potential founders, including FVC% pre, DLCO% pre and UIP patterns . Eight studies evaluated the relationship between a smoking history and mortality in RA-ILD. Having ever smoked and/or being a current smoker were defined as having a smoking history in this meta-analysis. Having ever smoked was evaluated in five studies, while having ever smoked and current smoking were evaluated in three studies. Having a smoking history was significantly associated with increased mortality in the meta-analysis (HRs = 1.42, 95% CI 1.03–1.96, P = 0.03) (Fig. 2C). Given the heterogeneity in defining having a smoking history, subgroup analysis was conducted. Current smoking and having ever smoked were not statistically associated with mortality in RA-ILD (HRs = 1.28, 95% CI 0.87–1.89, P = 0.205). There was also no a significant difference between having ever smoked and mortality in RA-ILD (HRs = 1.53, 95% CI 0.91–2.57, P = 0.105). Heterogeneity across the studies was low or moderate. Visual inspection of the funnel plot showed that one study fell outside the funnel plot with another borderline study, together suggesting the possibility of publication bias (Additional file 2: Fig. 2C). However, using the Egger’s test for publication bias, the P-value was 0.108 indicating that there was no publication bias.
RA-specific variables included RF, ACPA, disease activity score (DAS28), C-reactive protein (CRP), ESR, health assessment questionnaire (HAQ) disability score and RA disease duration in the included studies. Only when these variables were collected when RA-ILD was first diagnosed, could the data be pooled. Two studies evaluated the association between positive RF status and mortality in RA-ILD [6, 27]. The included studies did not report a relevant correlation between positive RF status and mortality (HRs = 1.00, 95% CI 1.00–1.00, P = 0.009) (Additional file 1: Fig. 1). Heterogeneity across the studies was low (I2:0%, P = 0.659). The other three studies reported that there was no relationship between positive RF status and mortality in RA-ILD patients [22, 26, 37]. However, the data could not be pooled due to the unclear time of RF collection and different statistical methods. One study demonstrated no evidence of an association between positive ACPA status and mortality , whereas an ACPA titer ≥ 500 U/mL was associated with mortality (HR = 3.95, 1.12–13.93, P = 0.033) in another study . However, the HRs could not be pooled due to a lack of the same ACPA threshold. One study reported that DAS28 at RA-ILD diagnosis in a multivariable model was a predictor of mortality after adjusting for potential confounders, including age, sex, calendar year of entry, disease duration, HAQ score, steroid use and methotrexate use (HR = 1.43, 95% CI 1.11–1.85) . Although the DAS28 was administered at an unclear time in another study, it showed that the DAS28 was not associated with mortality in RA-ILD patients (HR = 1.01, 95% CI 1.00–1.03). Two studies collected CRP when RA-ILD was first diagnosed and explored the impact of CRP on mortality in a multivariable model [26, 30]. There was generally an increased risk of mortality for RA-ILD patients with high levels of CRP. Due to the lack of data in univariate analysis, the pooled effect could not be estimated. There was no evidence of an association between ESR and mortality in two studies (HR = 1.01, 95% CI 1–1.03, P = 0.084 and HR = 1.06, 95% CI 0.96–1.18, per 10 mm/h increase, P = 0.25) [22, 30]. Two studies showed that there could be an association between the HAQ and mortality (OR = 2.30, 95% CI 1.10–8.00, P = 0.03, and HR = 1.54, 95% CI 1.05–2.26) [19, 22]. However, it was not clear when HAQ was administered in two studies. The effect of RA disease duration at ILD diagnosis was heterogeneous across studies. Solomon et al. found that RA disease duration was not associated with mortality (HR = 1.00, 95% CI 0.97–1.02, P = 0.77) , whereas Zamora-Legoff et al. found that RA disease duration was associated with mortality (HR = 1.71, 95% CI 1.22–2.39, P = 0.002 per 10 years increase) . Due to the difference in the estimate effect unit, the pooled effect could not be evaluated.
We evaluated ILD-physiological variables including DLCO% predicted, FVC% predicted, and composite physiological index (CPI) in the meta-analysis. Only when these variables were collected when RA-ILD was first diagnosed, could the data be pooled. DLCO% predicted was associated with mortality in four studies (HRs = 0.98, 95% CI 0.97–1, P = 0.008) (Fig. 2D) [6, 26, 35, 36]. Heterogeneity across the studies was moderate (I2:42.4%, P = 0.157). Due to the limited number of studies, funnel plots could not be estimated. The pooled HRs (HRs = 0.99, 95% CI 0.98–1, P = 0.021) showed that FVC% predicted was significantly associated with mortality in five studies (Fig. 2E) [6, 25, 26, 35, 36]. Heterogeneity across the studies was low (I2:28%, P = 0.235). Visual inspection of the funnel plot was symmetrical (Additional file 2: Fig. 2D). Using the Egger’s test for publication bias, the P-value was 0.568 indicating that there was no publication bias. CPI was associated with mortality (HRs = 1.04, 95% CI 1.02–1.06, P = 0.0001) in two studies (Fig. 2F) [27, 36]. Heterogeneity across the studies was low (I2:22.5%, P = 0.256). Due to the limited number of studies, funnel plots could not be estimated.
ILD-radiological variables, including UIP pattern, emphysema presence, and acute exacerbation of ILD, were evaluated in the meta-analysis. Due to moderate heterogeneity across the studies (I2:63.4%, P = 0.008), the HRs were pooled in the random-effect model (HRs = 1.88, 95% CI 1.14–3.10, P = 0.014) and the results showed that the presence of the UIP pattern was associated with mortality in the meta-analysis (Fig. 3A) [6, 25, 27, 28, 30, 35, 37, 39]. There was no publication bias according to the funnel plot and by using the Egger’s test (P = 0.327) (Additional file 2: Fig. 2E). In another seven studies, the UIP pattern was also associated with almost twofold higher mortality in the meta-analysis (RRs = 1.9, 95% CI 1.5–2.39, P = 0.0001) [27,28,29,30, 33, 37, 40] (Fig. 3B). Visual inspection of the funnel plot was asymmetric and the P-value of Egger’s test was 0.027 (Additional file 2: Fig. 2F). However, after using the trim and fill method, there was no significant difference in RRs. Although there was low publication bias, the result was stable. One study of 54 biopsy-proven rheumatoid lung disease cases also showed that the survival rates were significantly poorer in UIP patients than in non-UIP RA patients (P = 0.0452) . Three studies found that emphysema presence was associated with mortality in RA-ILD patients (HRs = 2.31, 95% CI 1.58–3.39, P = 0.0001) (Fig. 3C) [31, 37]. In the included studies, the relationship between other ILD radiological variables and mortality was also evaluated, such as traction bronchiectasis, reticulation, honeycombing extent, fibrosis, ground glass score and emphysema extent. One study showed that the extent of reticulation (HR 1.144, P = 0.041) and traction bronchiectasis (HR 1.184, P = 0.030) were associated with decreased survival . A large study of RA-ILD in the UK showed that extensive disease (> 20% of lung involvement) was associated with increased mortality compared with limited disease (RR = 2.17, 95% CI 1.14–4.13, P = 0.018) . In another study in China, Cox hazard analysis revealed that extensive lung involvement on HRCT (> 30% of lung involved) was associated with worse survival (HR = 2.401, 95% CI 1.075–5.362, P = 0.033) . However, the quantification method of chest images was not equivalent in the included studies, and the results could not be pooled. Acute exacerbation of ILD in two studies was associated with twofold higher mortality in the meta-analysis (HRs = 2.70, 95% CI 1.67–4.36, P = 0.0001) (Fig. 3D) [25, 36]. Heterogeneity across studies was low (I2:0%, P = 0.619). Due to the limited number of studies, funnel plots could not be estimated.
In this systematic review and meta-analysis, we aimed to identify factors associated with mortality among RA-ILD patients. We identified nine statistically significant factors associated with increased mortality: older age, male sex, smoking history, lower DLCO% predicted, lower FVC% predicted, higher CPI, UIP pattern on HRCT, emphysema presence and acute exacerbation of ILD. Knowledge of prognostic factors, which consist of clinical information that is easily accessible in daily clinical practice, will be of great help in developing therapeutic strategies for this intractable disease and can be very informative to patients and families in facilitating their decision-making.
Currently, the exact aetiology of RA-associated ILD is not known. However, some factors, such as older age, male sex, smoking and other serum biomarkers have been reported to be involved in the pathogenesis of RA-ILD. Among the identified prognostic factors, patient-specific factors, older age and male sex are analogous to prognostic variables in other ILD subtypes . For example, older age in systemic sclerosis and male sex in idiopathic pulmonary fibrosis predict worse prognosis. The biological effect of cigarette smoking on the development of RA-ILD remains unclear. However, it is known that cigarette smoke components, such as nicotine, hydrocarbons and carbon monoxide could trigger an immune reaction that produces serum autoantibodies against multiple citrullinated proteins in the lung, in turn leading to inflammation and epithelial cell injury, ultimately resulting in ILD [44, 45].
In this systematic review and meta-analysis, we evaluated several RA-specific variables, including RF, ACPA, DAS28, CRP, ESR, HAQ disability score and RA disease duration at ILD diagnosis. If these variables were collected when RA-ILD was first diagnosed, the results could be pooled quantitatively using meta-analysis methods. As we all know, RF and ACPA can be expressed in two different ways: qualitative (positive/negative) and quantitative (titers). However, due to a lack of the same RF/ACPA threshold and different statistical methods of ACPA, the quantitative pooled results of RF/ACPA and the qualitative pooled results of ACPA could not be gotten. Although high titres of RF and ACPA are biomarkers associated with an increased risk of developing RA-ILD [46, 47], there was no relationship between RF positive status and mortality in RA-ILD patients in this meta-analysis. Moreover, ESR had no association with all-cause mortality. The ACPA, DAS28, CRP, HAQ disability score and RA disease duration were collected when RA-ILD was first diagnosed or at unclear times in these included studies. Therefore, the data could not be pooled to obtain results because of the differences in statistical time and methodology.
Pulmonary function testing (PFT) in RA-ILD patients often reveals restrictive ventilatory defects with decreases in gas exchange, and such abnormalities can be detected even in the absence of any clinical symptoms . In patients with idiopathic pulmonary fibrosis, FVC and DLCO are the most sensitive parameters for assessing the clinical course. Moreover, there is a strong correlation between progressively lower values of FVC and DLCO and clinical disease severity from ILD [48, 49]. In our meta-analysis, we identified that reduced DLCO% and FVC% predicted were associated with higher mortality in RA-ILD. The CPI is a “severity variable” constructed against the morphologic severity of disease to calibrate the quantification of pulmonary fibrosis using pulmonary function tests in isolation. The study shows that the CPI correlates with the extent of pulmonary fibrosis on CT more strongly and is linked to mortality more closely than individual pulmonary function in IPF patients [50, 51]. In our meta-analysis, the CPI was associated with mortality in RA-ILD patients. Therefore, the CPI provided a valuable prognostic assessment at disease presentation.
The UIP pattern is the predominant HRCT pattern in RA-ILD and is characterized by heterogenous honeycombing at the bases and periphery of the lungs, peripheral basilar predominant reticular abnormalities and architectural distortions [21, 52]. The UIP pattern seen in RA patients is very similar to that observed in IPF patients, and it predicts worse prognosis compared with a non-UIP pattern of RA-ILD [16, 24, 53]. In never-smokers and smokers, the presence of emphysema was independently associated with the presence of honeycombing on HRCT and an HRCT UIP pattern . Emphysema adversely affects prognosis by increasing the likelihood of developing pulmonary hypertension . Acute exacerbation reportedly occurs not only in patients with IPF but also in patients with other interstitial lung diseases, including RA-ILD . The in-hospital mortality associated with acute exacerbation in RA-ILD patients was demonstrated to be as high as that in IPF patients [25, 55]. Therefore, acute exacerbation has been suggested to be associated with mortality in RA-ILD patients as well as patients with IPF.
This systematic review and meta-analysis had several limitations. First, the studies included were heterogeneous and had variable methodological quality. Moreover, most of the studies were retrospective in design, with small sample sizes and conducted in single medical institutions. Second, the search criteria might have potentially excluded eligible studies not in English due to search being limited to the English language. Third, the diagnostic criteria of RA and the pattern of ILD in the included studies were also variable. The 1987 ACR diagnostic criteria of RA and the 2010 ACR/EULAR diagnostic criteria of RA were used in this study. The sensitivity of the 2010 ACR/EULAR criteria was greater than that of the 1987 ACR criteria, while the specificity was lower. However, the included studies were observational studies with follow-up periods, and there was little possibility of misdiagnosis. The diagnostic criteria for ILD included findings on HRCT only, surgical lung biopsy only, HRCT and surgical lung biopsy, and a combination of clinical presentation, pulmonary function tests, and HRCT findings. Currently, there are no uniform diagnostic criteria or classification guidelines for RA-ILD. Due to the risks associated with surgical lung biopsies, most patients are diagnosed and classified into the different subtypes based only on HRCT. In our study, we evaluated the correlation between the UIP pattern and mortality in RA-ILD. The UIP pattern was classified by HRCT findings in 10 of 11 included studies. UIP patterns were classified by HRCT and/or lung biopsy in one study. However, all patients with the histopathology of UIP exhibited the same HRCT pattern in this study . ILD is a highly heterogenous disorder, and the diagnosis of ILD must involve the integration of multidisciplinary team meetings, including a respiratory physician, a radiologist and a rheumatologist, as well as a pathologist when applicable . In addition, the difference between the factors assessed in the different studies could also limit the value of the meta-analysis. Heterogeneity across some studies was low or moderate. Sensitivity analysis was limited due to the small number of studies in this meta-analysis. Therefore, a prospective cohort study with high quality involving multiple institutions in different countries is imperative to draw a definitive conclusion.
In summary, our meta-analysis identified the older age, male sex, smoking history, lower DLCO% predicted, lower FVC% predicted, UIP pattern on HRCT, emphysema presence, and acute exacerbation of ILD were associated with mortality in RA-ILD. The data should be interpreted cautiously due to the limited number of high-quality studies, substantial heterogeneity, and the limitations inherent in observational studies. This analysis could provide for better decision making by clinicians to ascertain whether early interventions targeting these factors can reduce mortality in high-risk RA-ILD patients.
Availability of data and materials
The datasets used and/or analyzed in this study are available from the corresponding author upon reasonable request.
Interstitial lung disease
Rheumatoid arthritis-associated interstitial lung disease
Pooled hazard ratios
Pooled risk ratios
American College of Rheumatology/European League Against Rheumatism
Usual interstitial pneumonia
Forced vital capacity
Lower diffusing capacity of the lung for carbon monoxide
Composite physiological index
Anti-citrullinated peptide antibodies
Erythrocyte sedimentation rate
Ground glass opacity
High resolution computed tomography
Disease activity score
Health assessment questionnaire
Turesson C, O’Fallon WM, Crowson CS, Gabriel SE, Matteson EL. Extra-articular disease manifestations in rheumatoid arthritis: incidence trends and risk factors over 46 years. Ann Rheum Dis. 2003;62:722–7.
Olson AL, Swigris JJ, Sprunger DB, Fischer A, Fernandez-Perez ER, Solomon J, Murphy J, Cohen M, Raghu G, Brown KK. Rheumatoid arthritis-interstitial lung disease-associated mortality. Am J Respir Crit Care Med. 2011;183:372–8.
Gabbay E, Tarala R, Will R, Carroll G, Adler B, Cameron D, Lake FR. Interstitial lung disease in recent onset rheumatoid arthritis. Am J Respir Crit Care Med. 1997;156:528–35.
Bongartz T, Nannini C, Medina-Velasquez YF, Achenbach SJ, Crowson CS, Ryu JH, Vassallo R, Gabriel SE, Matteson EL. Incidence and mortality of interstitial lung disease in rheumatoid arthritis: a population-based study. Arthritis Rheum. 2010;62:1583–91.
Raimundo K, Solomon JJ, Olson AL, Kong AM, Cole AL, Fischer A, Swigris JJ. Rheumatoid arthritis-interstitial lung disease in the United States: prevalence, incidence, and healthcare costs and mortality. J Rheumatol. 2019;46:360–9.
Hyldgaard C, Ellingsen T, Hilberg O, Bendstrup E. Rheumatoid arthritis-associated interstitial lung disease: clinical characteristics and predictors of mortality. Respiration. 2019;98:455–60.
Moher D, Liberati A, Tetzlaff J, Altman DG, Group P. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. Ann Intern Med. 2009;151:264–9, w64. https://doi.org/10.7326/0003-4819-151-4-200908180-00135.
Cumpston M, Li T, Page MJ, Chandler J, Welch VA, Higgins JP, Thomas J. Updated guidance for trusted systematic reviews: a new edition of the Cochrane Handbook for Systematic Reviews of Interventions. Cochrane Database Syst Rev. 2019;10:Ed000142.
Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ. 1997;315:629–34.
Duval S, Tweedie R. Trim and fill: a simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis. Biometrics. 2000;56:455–63.
Lee HK, Kim DS, Yoo B, Seo JB, Rho JY, Colby TV, Kitaichi M. Histopathologic pattern and clinical features of rheumatoid arthritis-associated interstitial lung disease. Chest. 2005;127:2019–27.
Lee YS, Kim HC, Lee BY, Lee CK, Kim MY, Jang SJ, Lee HS, Moon J, Colby TV, Kim DS. The value of biomarkers as predictors of outcome in patients with rheumatoid arthritis-associated usual interstitial pneumonia. Sarcoidosis Vasc Diffuse Lung Dis. 2016;33:216–23.
Nurmi HM, Purokivi MK, Kärkkäinen MS, Kettunen HP, Selander TA, Kaarteenaho RL. Variable course of disease of rheumatoid arthritis-associated usual interstitial pneumonia compared to other subtypes. BMC Pulm Med. 2016;16:107.
Nurmi HM, Purokivi MK, Kärkkäinen MS, Kettunen HP, Selander TA, Kaarteenaho RL. Are risk predicting models useful for estimating survival of patients with rheumatoid arthritis-associated interstitial lung disease? BMC Pulm Med. 2017;17:16.
Fu Q, Wang L, Li L, Li Y, Liu R, Zheng Y. Risk factors for progression and prognosis of rheumatoid arthritis–associated interstitial lung disease: single center study with a large sample of Chinese population. Clin Rheumatol. 2019;38:1109–16.
Jacob J, Hirani N, van Moorsel CHM, Rajagopalan S, Murchison JT, van Es HW, Bartholmai BJ, van Beek FT, Struik MHL, Stewart GA, et al. Predicting outcomes in rheumatoid arthritis related interstitial lung disease. Eur Respir J. 2019;53:1800869.
Kim HC, Choi KH, Jacob J, Song JW. Prognostic role of blood KL-6 in rheumatoid arthritis–associated interstitial lung disease. PLoS ONE. 2020;15:e0229997.
Yamakawa H, Sato S, Nishizawa T, Kawabe R, Oba T, Kato A, Horikoshi M, Akasaka K, Amano M, Sasaki H, et al. Impact of radiological honeycombing in rheumatoid arthritis-associated interstitial lung disease. BMC Pulm Med. 2020;20:25.
Wolfe F, Caplan L, Michaud K. Rheumatoid arthritis treatment and the risk of severe interstitial lung disease. Scand J Rheumatol. 2007;36:172–8.
Dixon WG, Hyrich KL, Watson KD, Lunt M, Consortium BCC, Symmons DP, British Society for Rheumatology Biologics R. Influence of anti-TNF therapy on mortality in patients with rheumatoid arthritis-associated interstitial lung disease: results from the British Society for Rheumatology Biologics Register. Ann Rheum Dis. 2010;69:1086–91.
Kim EJ, Elicker BM, Maldonado F, Webb WR, Ryu JH, Van Uden JH, Lee JS, King TE Jr, Collard HR. Usual interstitial pneumonia in rheumatoid arthritis-associated interstitial lung disease. Eur Respir J. 2010;35:1322–8.
Koduri G, Norton S, Young A, Cox N, Davies P, Devlin J, Dixey J, Gough A, Prouse P, Winfield J, et al. Interstitial lung disease has a poor prognosis in rheumatoid arthritis: results from an inception cohort. Rheumatology. 2010;49:1483–9.
Nakamura Y, Suda T, Kaida Y, Kono M, Hozumi H, Hashimoto D, Enomoto N, Fujisawa T, Inui N, Imokawa S, et al. Rheumatoid lung disease: prognostic analysis of 54 biopsy-proven cases. Respir Med. 2012;106:1164–9.
Solomon JJ, Ryu JH, Tazelaar HD, Myers JL, Tuder R, Cool CD, Curran-Everett D, Fischer A, Swigris JJ, Brown KK. Fibrosing interstitial pneumonia predicts survival in patients with rheumatoid arthritis-associated interstitial lung disease (RA-ILD). Respir Med. 2013;107:1247–52.
Hozumi H, Nakamura Y, Johkoh T, Sumikawa H, Colby TV, Kono M, Hashimoto D, Enomoto N, Fujisawa T, Inui N, et al. Acute exacerbation in rheumatoid arthritis-associated interstitial lung disease: a retrospective case control study. BMJ Open. 2013;3:e003132.
Song JW, Lee HK, Lee CK, Chae EJ, Jang SJ, Colby TV, Kim DS. Clinical course and outcome of rheumatoid arthritis-related usual interstitial pneumonia. Sarcoidosis Vasc Diffuse Lung Dis. 2013;30:103–12.
Solomon JJ, Chung JH, Cosgrove GP, Demoruelle MK, Fernandez-Perez ER, Fischer A, Frankel SK, Hobbs SB, Huie TJ, Ketzer J, et al. Predictors of mortality in rheumatoid arthritis-associated interstitial lung disease. Eur Respir J. 2016;47:588–96.
Rojas-Serrano J, Herrera-Bringas D, Pérez-Román DI, Pérez-Dorame R, Mateos-Toledo H, Mejía M. Rheumatoid arthritis-related interstitial lung disease (RA-ILD): methotrexate and the severity of lung disease are associated to prognosis. Clin Rheumatol. 2017;36:1493–500.
Yang JA, Lee JS, Park JK, Lee EB, Song YW, Lee E. Clinical characteristics associated with occurrence and poor prognosis of interstitial lung disease in rheumatoid arthritis. Korean J Intern Med. 2017;34:434–41.
Zamora-Legoff JA, Krause ML, Crowson CS, Ryu JH, Matteson EL. Patterns of interstitial lung disease and mortality in rheumatoid arthritis. Rheumatology (Oxford). 2017;56:344–50. https://doi.org/10.1093/rheumatology/kew391.
Jacob J, Song JW, Yoon HY, Cross G, Barnett J, Woo WL, Adams F, Kokosi M, Devaraj A, Renzoni E, et al. Prevalence and effects of emphysema in never-smokers with rheumatoid arthritis interstitial lung disease. EBioMedicine. 2018;28:303–10.
Kelly CA, Saravanan V, Nisar M, Arthanari S, Woodhead FA, Price-Forbes AN, Dawson J, Sathi N, Ahmad Y, Koduri G, Young A. Rheumatoid arthritis-related interstitial lung disease: associations, prognostic factors and physiological and radiological characteristics—a large multicentre UK study. Rheumatology (United Kingdom). 2014;53:1676–82.
Koo BS, Hong S, Kim YJ, Kim YG, Lee CK, Yoo B. Mortality in patients with rheumatoid arthritis-associated interstitial lung disease treated with an anti-tumor necrosis factor agent. Korean J Intern Med. 2015;30:104–9.
Nurmi HM, Kettunen HP, Suoranta SK, Purokivi MK, Kärkkäinen MS, Selander TA, Kaarteenaho RL. Several high-resolution computed tomography findings associate with survival and clinical features in rheumatoid arthritis-associated interstitial lung disease. Respir Med. 2018;134:24–30.
Ito Y, Arita M, Kumagai S, Takei R, Noyama M, Tokioka F, Nishimura K, Koyama T, Tachibana H, Ishida T. Radiological fibrosis score is strongly associated with worse survival in rheumatoid arthritis-related interstitial lung disease. Mod Rheumatol. 2019;29:98–104.
Yamakawa H, Sato S, Tsumiyama E, Nishizawa T, Kawabe R, Oba T, Kamikawa T, Horikoshi M, Akasaka K, Amano M, et al. Predictive factors of mortality in rheumatoid arthritis-associated interstitial lung disease analysed by modified HRCT classification of idiopathic pulmonary fibrosis according to the 2018 ATS/ERS/JRS/ALAT criteria. J Thorac Dis. 2019;11:5247–57.
Kim HC, Lee JS, Lee EY, Ha YJ, Chae EJ, Han M, Cross G, Barnett J, Joseph J, Song JW. Risk prediction model in rheumatoid arthritis-associated interstitial lung disease. Respirology. 2020;25:1257–64.
Li L, Liu R, Zhang Y, Zhou J, Li Y, Xu Y, Gao S, Zheng Y. A retrospective study on the predictive implications of clinical characteristics and therapeutic management in patients with rheumatoid arthritis-associated interstitial lung disease. Clin Rheumatol. 2020;39:1457–70.
Wang N, Zhang Q, Jing X, Guo J, Huang H, Xu Z. The association between MUC5B mutations and clinical outcome in patients with rheumatoid arthritis-associated interstitial lung disease: a retrospective exploratory study in China. Med Sci Monit. 2020;26:e920137.
Tsuchiya Y, Takayanagi N, Sugiura H, Miyahara Y, Tokunaga D, Kawabata Y, Sugita Y. Lung diseases directly associated with rheumatoid arthritis and their relationship to outcome. Eur Respir J. 2011;37:1411–7.
Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988;31:315–24.
Aletaha D, Neogi T, Silman AJ, Funovits J, Felson DT, Bingham CO, Birnbaum NS, Burmester GR, Bykerk VP, Cohen MD, et al. 2010 rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Ann Rheum Dis. 2010;69:1580–8.
Winstone TA, Assayag D, Wilcox PG, Dunne JV, Hague CJ, Leipsic J, Collard HR, Ryerson CJ. Predictors of mortality and progression in scleroderma-associated interstitial lung disease: a systematic review. Chest. 2014;146:422–36.
Gochuico BR, Avila NA, Chow CK, Novero LJ, Wu HP, Ren P, MacDonald SD, Travis WD, Stylianou MP, Rosas IO. Progressive preclinical interstitial lung disease in rheumatoid arthritis. Arch Intern Med. 2008;168:159–66.
Albano SA, Santana-Sahagun E, Weisman MH. Cigarette smoking and rheumatoid arthritis. Semin Arthritis Rheum. 2001;31:146–59.
Mori S, Koga Y, Sugimoto M. Different risk factors between interstitial lung disease and airway disease in rheumatoid arthritis. Respir Med. 2012;106:1591–9.
Zhu J, Zhou Y, Chen X. A meta analysis of the increased risk of rheumatoid arthritis-related pulmonary disease as a result of serum anticitrullinated protein antibody positivity. J Rheumatol. 2014;41:1282–9.
Ley B, Collard HR, King TE. Clinical course and prediction of survival in idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2011;183:431–40.
Zamora-Legoff JA, Krause ML, Crowson CS, Ryu JH, Matteson E. Progressive decline of lung function in rheumatoid arthritis-associated interstitial lung disease. Arthritis Rheumatol. 2017;69:542–9.
Sharp C, Adamali HI, Millar AB. A comparison of published multidimensional indices to predict outcome in idiopathic pulmonary fibrosis. ERJ Open Res. 2017. https://doi.org/10.1183/23120541.00096-2016.
Wells AU, Desai SR, Rubens MB, Goh NS, Cramer D, Nicholson AG, Colby TV, du Bois RM, Hansell DM. Idiopathic pulmonary fibrosis: a composite physiologic index derived from disease extent observed by computed tomography. Am J Respir Crit Care Med. 2003;167:962–9.
Doyle TJ, Dellaripa PF. Lung Manifestations in the rheumatic diseases. Chest. 2017;152:1283–95.
Travis WD, Costabel U, Hansell DM, King TE Jr, Lynch DA, Nicholson AG, Ryerson CJ, Ryu JH, Selman M, Wells AU, et al. An official American Thoracic Society/European Respiratory Society statement: update of the international multidisciplinary classification of the idiopathic interstitial pneumonias. Am J Respir Crit Care Med. 2013;188:733–48.
Cottin V, Nunes H, Mouthon L, Gamondes D, Lazor R, Hachulla E, Revel D, Valeyre D, Cordier JF, Groupe d'Etudes et de Recherche sur les Maladies “Orphelines” P. Combined pulmonary fibrosis and emphysema syndrome in connective tissue disease. Arthritis Rheum. 2011;63:295–304.
Park IN, Kim DS, Shim TS, Lim CM, Lee SD, Koh Y, Kim WS, Kim WD, Jang SJ, Colby TV. Acute exacerbation of interstitial pneumonia other than idiopathic pulmonary fibrosis. Chest. 2007;132:214–20.
Walsh SLF, Wells AU, Desai SR, Poletti V, Piciucchi S, Dubini A, Nunes H, Valeyre D, Brillet PY, Kambouchner M, et al. Multicentre evaluation of multidisciplinary team meeting agreement on diagnosis in diffuse parenchymal lung disease: a case-cohort study. Lancet Respir Med. 2016;4:557–65.
Assayag D, Lubin M, Lee JS, King TE, Collard HR, Ryerson CJ. Predictors of mortality in rheumatoid arthritis-related interstitial lung disease. Respirology. 2014;19:493–550.
We really appreciate the efforts of all the researchers whose articles were included in this study. Furthermore, we also acknowledge previous systematic reviews on the same topic .
This work was supported by grants from the Science and Technology Program of Yantai City (No. 2020YD024).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Forrest plots of the meta-analysis for RF. RF rheumatoid factor.
Funnel plots of publication bias in age (A), male sex (B), smoking history (C), FVC% predicted (D), UIP with HRs (E), and UIP with RRs (F). HRs hazard ratios; RRs relative ratios; DLCO diffusing capacity of the lung for carbon monoxide; FVC forced vital capacity; UIP usual interstitial pneumonia.
Detailed search strategy and results.
The quality of the included studies in the meta-analysis based on the Newcastle–Ottawa scale.
About this article
Cite this article
Qiu, M., Jiang, J., Nian, X. et al. Factors associated with mortality in rheumatoid arthritis-associated interstitial lung disease: a systematic review and meta-analysis. Respir Res 22, 264 (2021). https://doi.org/10.1186/s12931-021-01856-z
- Interstitial lung disease
- Rheumatoid arthritis