Abstract
Objective To determine the association between novel lifestyle factors on risk of rheumatoid arthritis (RA)-associated interstitial lung disease (ILD), define the threshold at which smoking increases RA-ILD risk, and calculate the degree to which known lifestyle and clinical factors predict RA-ILD.
Methods This nested case-control study matched incident RA-ILD cases to RA non-ILD controls on age, sex, RA duration, rheumatoid factor, and time from exposure assessment to RA-ILD. Exposures included education, BMI, smoking, anticyclic citrullinated peptide antibodies, race, joint erosions, rheumatoid nodules, C-reactive protein (CRP), disease activity score, functional status, disease-modifying antirheumatic drug use, and glucocorticoid use. OR for each exposure on risk of RA-ILD were obtained from logistic regression models. Area under the curve (AUC) was calculated based on all lifestyle and clinical exposures.
Results We identified 84 incident RA-ILD cases and 233 matched controls. After adjustment, obesity, high-positive CRP (≥ 10 mg/L), and poor functional status (multidimensional Health Assessment Questionnaire [MDHAQ] ≥ 1) were associated with increased risk of RA-ILD (OR 2.42, 95% CI 1.11–5.24 vs normal BMI; OR 2.61, 95% CI 1.21–5.64 vs CRP < 3 mg/L; OR 3.10, 95% CI 1.32–7.26 vs MDHAQ < 0.2). Smoking 30 pack-years or more was strongly associated with risk of RA-ILD compared to never smokers (OR 6.06, 95% CI 2.72–13.5). Together, lifestyle and clinical risk factors for RA-ILD had an AUC of 0.79 (95% CI 0.73–0.85).
Conclusion Obesity, CRP, functional status, and extensive smoking may be novel risk factors for RA-ILD that may be useful for RA-ILD risk assessment and prevention. The overall ability to predict RA-ILD remains modest.
Approximately 5–10% of patients with rheumatoid arthritis (RA) have clinically significant RA-associated interstitial lung disease (ILD), and an additional 20–30% may have subclinical RA-ILD1,2,3. Further, RA-ILD is associated with a median survival of fewer than 3 years after diagnosis2, with 2–10 times the odds of mortality compared to RA patients without ILD4,5. Identifying predictors of RA-ILD is crucial to better understand the underlying pathogenesis and to allow for earlier diagnosis and treatment, with the goal of preventing irreversible lung damage. In the future, ascertaining predictors of RA-ILD may even build the foundation for interventions that can improve its poor survival by early detection or altering the natural history prior to clinically severe RA-ILD.
There are several gaps in the RA-ILD literature. First, no studies have investigated the association between lifestyle factors such as education or obesity on RA-ILD risk, despite evidence that they are risk factors for RA6,7,8,9. Second, while many studies have shown smoking to be a risk factor for RA-ILD10–16, a few have not found an association17,18,19, possibly due to a threshold effect for smoking on RA-ILD risk. One study showed that smoking over 25 pack-years was associated with a dramatic increase in physiologic and radiologic abnormalities suggestive of RA-ILD compared to fewer than 25 pack-years, but this study did not explore additional smoking exposure thresholds or stratification groups, and it used pulmonary function abnormalities as the outcome instead of RA-ILD10. Third, existing studies of RA-ILD risk have typically lacked comprehensive adjustment for many potential confounders1,5,10–29, resulting in discordant results for certain risk factors including male sex, RA duration, rheumatoid factor (RF) positivity, disease activity, functional status, prednisone use, and rheumatoid nodules. Further, the degree to which known lifestyle and clinical risk factors together predict risk of RA-ILD remains unknown.
To address these gaps, our objectives were 3-fold: (1) To elucidate the association of lifestyle factors such as education and obesity with subsequent risk of RA-ILD; (2) to define the pack-years threshold at which smoking increases risk of RA-ILD; and (3) to determine the degree to which a comprehensive assessment of known risk factors for RA-ILD predicts disease development. We hypothesized that lower education, obesity, and smoking over a threshold of 20 pack-years would be associated with increased risk of RA-ILD but that all risk factors together would have only a modest ability to predict RA-ILD, pointing to the likely existence of unmeasured risk factors.
MATERIALS AND METHODS
Study design and participants. This nested case-control study was performed within the Brigham Rheumatoid Arthritis Sequential Study (BRASS), an ongoing large, single-center, prospective registry of patients with RA since 200330. With its low patient turnover, this population may have reasonable generalizability to other RA populations. For example, the mean age (57 yrs), and proportions of females (82%), new-onset RA (20%), and high functional status (82%) are nearly identical to RA patients in the Consortium of Rheumatology Researchers of North America (CORRONA) registry (59%, 75%, 20%, and 88%, respectively)30,31. All RA cases met 1987 American College of Rheumatology (ACR) or 2010 ACR/European League Against Rheumatism criteria as determined by a rheumatologist32,33, and those with verified RA were eligible to participate30. Blood samples, RA characteristics, and physician/patient-reported measures were collected at enrollment and semiannually30. This BRASS substudy received approval from the Partners HealthCare Institutional Review Board (2016P000226), followed STrengthening the Reporting of OBservational studies in Epidemiology (STROBE) guidelines for observational studies, and complied with the Declaration of Helsinki.
RA-ILD cases and RA non-ILD controls. High-resolution chest computed tomography (CT) scans of BRASS participants performed for clinical care up until April 14, 2016, were reviewed by 2 attending thoracic radiologists and 1 attending pulmonologist. RA-ILD was defined as presence of interstitial lung abnormalities on the research read of chest CT, as described previously28. In a previous study examining the same RA-ILD cases used in this analysis, the clinical significance of CT-defined ILD was found to be high in several ways. First, of 78 RA-ILD cases for whom clinical radiographic reports were available, 68 (87%) had features consistent with RA-ILD noted on the report. Second, 64 (75%) of the 85 RA-ILD cases had additional clinical evaluation and/or follow-up for the condition. In addition, 32 (38%) had changes to RA medications after CT scan. Finally, 32 (38%) of the RA-ILD cases died during study follow-up compared to only 126 (9%) in the previously analyzed RA population28. Controls were BRASS participants with absence of patient-reported ILD and absence of both billing codes and physician-reported history of asthma, chronic obstructive pulmonary disease, pulmonary fibrosis, bronchiectasis, bronchiolitis obliterans-organizing pneumonia, drug-induced pneumonitis, or tuberculosis. We did not require that all controls have chest CT performed, since doing so could introduce selection bias due to the clinical indication of the imaging. Each case was matched to up to 3 controls on age (within 5 yrs), sex, RA duration (within 5 yrs), RF positivity, and time from exposure assessment to index date (within 6 months). We matched by the above-known RA-ILD risk factors to improve power to investigate possible novel associations, knowing that this would slightly decrease the model’s ability to predict RA-ILD. Index date was the date of RA-ILD by the initial CT chest for the case and matched date for corresponding controls.
Exposures and covariates. To ensure that assessment of exposures preceded the date of ILD, we identified the earliest study visit date containing these data occurring before the index date. Exposure and covariate variables obtained from BRASS study visits included age, sex, RA duration, race/ethnicity (white, non-Hispanic vs other), education (lower than college degree vs higher), joint erosions (present vs absent on hand plain radiographs), rheumatoid nodules (present vs absent), BMI categories, smoking status (never, past, or current), smoking pack-years, disease activity score with 28-joint counts based on C-reactive protein (DAS28-CRP; remission [< 2.6], low [2.6 to < 3.2], moderate [3.2 to < 5.2], or high [≥ 5.1])34, functional status by multidimensional Health Assessment Questionnaire (MDHAQ; 0 to < 0.2, 0.2 < 1.0, or ≥ 1.0)35, and biologic disease-modifying antirheumatic drug (bDMARD), methotrexate (MTX), or prednisone use (never, past, current). RF positivity, anticyclic citrullinated peptide (anti-CCP) positivity, and CRP (normal [< 3 mg/L], low-positive [3 to < 10 mg/L], or high [≥ 10 mg/L]) came from study visit blood work.
Statistical analysis. We compared continuous variables using Wilcoxon rank-sum tests, and categorical variables using chi-square tests. For the association analyses, conditional logistic regression models of each exposure were fitted to obtain OR with 95% CI for RA-ILD case status, adjusting for anti-CCP, race/ethnicity, education, BMI (continuous), and smoking pack-years (continuous). Since disease activity and treatment closely interact, the main models of these factors did not adjust for each other. Exploratory sensitivity analyses added prednisone and MTX use to the disease activity model, and disease activity to each drug model, in order to investigate possible independent associations. For ordinal variables, we calculated P for trend by first calculating the median value within each category and then using those values as a continuous variable in a separate model. To determine the degree to which all risk factors predicted RA-ILD risk, we fit receiver-operating characteristic (ROC) curves and calculated the area under the curve (AUC). The threshold for statistical significance was a 2-sided P < 0.05. No participants were missing any data. Analyses were prespecified in a protocol and performed using SAS version 9.4 (SAS Institute Inc.).
RESULTS
Characteristics of RA-ILD cases and controls. Of the 1100 RA cases included in BRASS, we identified 84 cases of confirmed RA-ILD and 233 matched RA non-ILD controls. Among controls, 60 (26%) had at least 1 chest CT performed, and an additional 76 (33%) had at least 1 chest radiograph. Therefore, the majority (58%) of controls had chest imaging performed that did not clinically reveal RA-ILD. Characteristics of controls with CT compared to without CT were similar, with none of the factors associated with RA-ILD being associated with receipt of chest CT in controls, including education (Supplementary Table 1, available with the online version of this article). The median time between BRASS study visit where exposures/covariates were assessed and the index date of RA-ILD (or matched date) was 1.5 (IQR 0.6–2.5) years for RA-ILD cases and 1.9 (IQR 1.6–5.6) years for controls. At the time of exposure assessment, RA-ILD cases had more rheumatoid nodules, greater history of smoking, higher CRP, higher disease activity (DAS28-CRP), worse functional status (MDHAQ), less bDMARD and MTX use, and more prednisone use when compared to controls based on univariate analyses (Table 1).
Education and obesity. After further adjusting for anti-CCP, race, and smoking, we found that participants with lower education had a decreased risk of RA-ILD compared to participants with a college education or higher (Table 2). Obesity (BMI ≥ 30 kg/m2) was associated with a more than 2-fold increased risk of RA-ILD, with evidence of a dose-response relationship (P for trend 0.02; Table 2).
Smoking threshold. After adjustment, cumulative past smoking was not associated with RA-ILD (Table 3). However, current smoking was associated with a more than 3-fold increase in odds of RA-ILD compared to never smoking (Table 3). Further, smoking 30 pack-years or more was associated with a 6-fold increased risk compared to nonsmokers (P for trend < 0.001; Table 3). In contrast, smoking levels < 30 pack-years were not associated with increased risk of RA-ILD (Table 3).
Other risk factors for RA-ILD. Other factors that remained statistically significant predictors of RA-ILD after adjustment included high-positive CRP, moderate or high disease activity (DAS28-CRP 3.2 or higher), poor functional status (MDHAQ ≥ 1.0), and current prednisone use, while past or current MTX use was associated with decreased risk of RA-ILD (Table 4). A sensitivity analysis adjusting disease activity for prednisone and MTX use showed that moderate and high disease activity remained statistically significant predictors of RA-ILD (OR 2.22, 95% CI 1.04–4.72 and OR 5.09, 95% CI 1.58–16.4), respectively, compared to remission (P for trend 0.003). Similarly, the association between prednisone and RA-ILD remained significant even after additionally adjusting for disease activity (OR 4.08, 95% CI 1.23–13.6).
Ability to predict RA-ILD. Together, these known demographics (age, race, education), lifestyle (BMI, smoking pack-yrs), and RA-related (RA duration, RF positivity, anti-CCP positivity, joint erosions, rheumatoid nodules, bDMARD use, MTX use, prednisone use, disease activity, functional status) risk factors explain a modest amount of an individual’s risk for RA-ILD (Figure 1). The AUC for all factors combined was 0.79 (95% CI 0.73–0.85).
DISCUSSION
This case-control study discovered several novel predictors for incident RA-ILD, including obesity, higher education level, high-positive CRP, and poor functional status. We also made the discovery that smoking beyond a threshold of 30 pack-years had a strong association with RA-ILD, while there was no evidence of an association for lower levels of smoking. Third, we confirmed existing RA-ILD risk factors such as disease activity, suggesting disease control could reduce RA-ILD risk, and prednisone, which should be used with caution in patients with RA who are at risk for RA-ILD. Finally, we demonstrated that currently known clinical predictors of RA-ILD have only a modest ability to predict RA-ILD. These findings may shed light on the potential pathogenesis of RA-ILD, could improve patient-specific risk assessment and prevention efforts, and highlight the need for further studies on risk factors for RA-ILD.
We identified obesity and higher education as possible novel lifestyle risk factors for RA-ILD. Only 1 prior study reported data on the association between BMI and RA-ILD risk and found no association29, but it had only 29 RA-ILD cases. Thus, the association between obesity and RA-ILD is a novel finding to this study. Obesity has a modest but positive association with risk of RA, supporting the biological plausibility of such an association with RA-ILD8,9. It is possible that the atelectasis or poorer-quality CT scans that occur in obese persons could result in misclassification of RA-ILD in either direction. Nevertheless, this discovery has important clinical implications, as obesity is now a second potentially modifiable risk factor for RA-ILD on which clinicians can counsel their patients.
The association between higher education and RA-ILD is also novel to this study. The 1 prior study of RA-ILD that reported data on education showed no association, though it was limited by small sample size17. This association between higher education and RA-ILD was surprising since higher socioeconomic status has been shown to have no association with ILD36 or RA37,38. In fact, higher education has even been associated with a lower risk of RA6,7 and other autoimmune diseases like systemic lupus erythematosus39. Thus, high medical literacy, including willingness to undergo CT imaging, may explain this association rather than a true biologic effect. Future studies should confirm the association between education and RA-ILD and investigate the mechanism for it; this may provide clues to the pathogenesis of RA-ILD. Conversely, this association may be explained by health literacy as a construct of social determinants of health.
We also identified high-positive CRP and poor functional status as novel clinical risk factors for RA-ILD. Two prior studies showed no association between CRP and RA-ILD in multivariable models, but only included CRP as a continuous variable instead of separating out high CRP from those with low or normal CRP18,19. The validity of this association is supported by the association with ESR with RA-ILD10,14,16,17. Thus, high-positive CRP may be a useful biomarker for predicting increased risk of developing RA-ILD. Similarly, we found reduced functional status (in particular, MDHAQ ≥ 1) was associated with increased risk of RA-ILD. This is consistent with a prior study that showed a univariate association between poor functional status and RA-ILD10. Although 2 studies have shown no association between functional status and RA-ILD after adjusting for confounders, they were smaller than the current study11,17. Therefore, this study showed high-positive CRP (≥ 10 mg/L) and poor functional status (MDHAQ ≥ 1) to be risk factors for RA-ILD that can aid clinicians in prognosticating which patients are at increased risk for RA-ILD.
Another key finding from this study is that smoking may exhibit a threshold effect in its relationship with RA-ILD. We found that smoking 30 pack-years or more was associated with a 6-fold increase in RA-ILD risk, whereas smoking under this threshold was not associated with increased risk. Although most studies of RA-ILD risk show an association with smoking10,11,12,13,15,16, some others do not17,18,19. Based on our findings, a reason for this discrepancy could be that the smokers in the negative studies did not have sufficient (> 30 pack-yrs) smoking history. Alternatively, the lack of association could have resulted from lack of certain HLA-DR alleles, termed the “shared epitope,” as smoking may only increase RA-ILD risk in the presence of the shared epitope14. Interestingly, this 30 pack-year threshold for RA-ILD differs from RA, where the relationship between smoking and RA risk increases at a threshold of 5–10 pack-years and plateaus at 2040. Such dissimilarity suggests that the mechanism for the association of smoking with RA differs from its association with RA-ILD. Additionally, current smokers had increased risk for RA-ILD, confirming the findings from 2 previous small cross-sectional studies11,29. Together, these results underscore the importance of smoking cessation to prevent RA-ILD, especially in patients who have not yet reached the critical 30 pack-year threshold. Future studies of RA-ILD risk should consider stratifying by pack-years rather than ever/never smoking and investigate the pathophysiologic reasons for a threshold effect.
With its multivariable models, this study also clarified the association between several possible risk factors and RA-ILD. For example, some studies have shown a positive association for rheumatoid nodules10,14,16 and disease activity11,14,28, while other studies have shown no association between rheumatoid nodules11,17,28 or disease activity16,17,18 and risk of RA-ILD. We found that moderate to high disease activity was associated with increased risk of RA-ILD, suggesting that controlling disease activity can reduce risk. We did not find an association with rheumatoid nodules, though matching on RF may have reduced the association. Similarly, some studies have shown an association between prednisone use and RA-ILD14,19,21,28, whereas others have not11,13,18. In this study, even after accounting for disease activity, current prednisone use was associated with a 4-fold increased risk for RA-ILD. Thus, providers should use caution in prescribing prednisone in RA or even RA-ILD patients.
The lack of association with MTX use, while initially surprising given the reported association between MTX and lung disease41, has actually been shown repeatedly in prior studies of RA-ILD risk11,13,14,18,19,21,28. In fact, 3 prior studies have even shown MTX may have a protective effect on RA-ILD risk, similar to our present study10,16,20. While it is possible that MTX may protect against RA-ILD given its therapeutic effects in RA, this decreased association might also simply represent channeling bias away from MTX use in patients with lung disease. Even so, given the lack of treatments for RA-ILD, further exploration of these 2 possibilities may be warranted.
Finally, despite discovery of these additional risk factors for RA-ILD, we found that these RA-ILD risk factors had only a modest AUC of 0.79. This finding is consistent with prior studies, whose AUC ranged from 0.60 to 0.8513,14,29, and suggests a need to identify additional risk factors for RA-ILD. Fortunately, recent studies have shown promise for biomarkers like fine-specificity anticitrullinated peptide antibodies11,42,43,44,45, matrix metalloproteinase 713, surfactant protein D13,46, and pulmonary and activation-regulated chemokine13, as well as genetic markers such as the HLA shared epitope14 or MUC5B promoter variant47. However, these biomarkers are not usually available in routine clinical care. Further studies incorporating these biologic and genetic markers are needed.
Strengths of this study include its inclusion of numerous previously studied risk factors for RA-ILD and their prospective collection. Nevertheless, there are several important limitations. First, ILD can begin insidiously before CT is obtained, so it is possible that the predictors may have resulted from RA-ILD. Even if this were true, data were collected a median of 1.5 years before clinical detection, so these factors could still identify patients in an earlier phase in the course of RA-ILD. Second, selection bias is possible since not all BRASS participants had chest CT scans available for review, as these were performed for clinical purposes. For example, some controls might have had unrecognized, subclinical ILD since chest CT were performed on only 25% of controls, which would be expected to bias results toward the null. We mitigated this possibility by excluding controls with comorbidities that could be suggestive of ILD. Third, RA-ILD patients may be prone to differential follow-up, compared to patients with RA and no ILD, due to the serious nature of this condition. To mitigate any potential selection bias related to right/left censoring, our case-control study matched on time from exposure/covariate assessment to index/matched date. Fourth, matching controls on RA duration and RF positivity, while improving comparability of the groups in this study and power to detect novel associations, likely eliminated the association with anti-CCP and slightly lowered the calculated AUC since these factors were balanced between cases and controls by design. Finally, unmeasured factors exist and may be important for RA-ILD risk, such as biomarkers and genetic factors like the shared epitope and MUC5B variant.
In conclusion, we identified obesity, higher education, high-level CRP, and poor functional status as novel risk factors for RA-ILD. In addition, smoking 30 pack-years or greater was associated with a substantially increased risk of RA-ILD, whereas lower levels of smoking did not increase risk. These findings might improve our understanding of RA-ILD pathogenesis and prediction, and also suggest weight loss and smoking cessation might prevent RA-ILD. However, all known clinical risk factors for RA-ILD have only a modest ability to predict RA-ILD development. Given the poor prognosis and limited available treatment options for RA-ILD, this finding emphasizes the critical need for future studies of clinical, biological, and genetic risk factors for RA-ILD.
ACKNOWLEDGMENT
We thank Lily Martin for performing our technical review for this manuscript.
Footnotes
T.J. Doyle and J.A. Sparks contributed equally to this manuscript.
The authors declare no conflicts of interest.
Full Release Article. For details see Reprints and Permissions at jrheum.org.
This work was supported by the Rheumatology Research Foundation K Supplement Award; National Institutes of Health (NIH; grant numbers K23 AR069688, K23 HL119558, R03 AR075886, R03 HL148484, L30 AR066953, P30 AR070253, and P30 AR072577); and the R. Bruce and Joan M. Mickey Research Scholar Fund. The Brigham Rheumatoid Arthritis Sequential Study is funded by grants from Bristol-Myers Squibb, Crescendo Bioscience, and Sanofi. The funders had no role in study design, data collection, analysis, decision to publish, or preparation of the manuscript. The content is solely the responsibility of the authors and does not necessarily represent the official views of Harvard University, its affiliated academic health care centers, or the NIH.
- Accepted for publication October 23, 2020.
- Copyright © 2021 by the Journal of Rheumatology
Free online via JRheum Full Release option