Development and Internal Validation of a Clinical and Genetic Risk Score for Rheumatoid Arthritis-Associated Interstitial Lung Disease.

OBJECTIVE: Although clinical and genetic risk factors have been identified for rheumatoid arthritis-associated interstitial lung disease (RA-ILD), there are no current tools allowing for risk stratification. We sought to develop and validate an ILD risk model in a large, multicentre, prospective RA cohort. METHODS: Participants in the Veterans Affairs RA (VARA) registry were genotyped for 12 single nucleotide polymorphisms (SNPs) associated with idiopathic pulmonary fibrosis. ILD was validated through systematic record review. A genetic risk score (GRS) was computed from minor alleles weighted by effect size with ILD, using backward selection. The GRS was combined with clinical risk factors within a logistic regression model. Internal validation was completed using bootstrapping, and model performance was assessed by the area under the receiver operating curve (AUC). RESULTS: Of 2,386 participants (89% male, mean age 69.5 years), 9.4% had ILD. Following backward selection, five SNPs contributed to the GRS. The GRS and clinical factors outperformed clinical factors alone in discriminating ILD (AUC 0.675 vs 0.635, p< 0.001). The shrinkage-corrected performance for combined and clinical-only models was 0.667 (95% CI 0.628, 0.712) and 0.623 (95% CI 0.584, 0.651), respectively. Twenty percent of the cohort had a combined risk score below a cut-point with >90% sensitivity. CONCLUSION: A clinical and genetic risk model discriminated ILD in a large, multicentre RA cohort better than a clinical-only model, excluding 20% of the cohort from low-yield testing. These results demonstrate the potential utility of a GRS in RA-ILD and support further investigation into individualized risk stratification and screening.

Extracting forced vital capacity from the electronic health record through natural language processing in rheumatoid arthritis-associated interstitial lung disease.

PURPOSE: To develop a natural language processing (NLP) tool to extract forced vital capacity (FVC) values from electronic health record (EHR) notes in patients with rheumatoid arthritis-interstitial lung disease (RA-ILD). METHODS: We selected RA-ILD patients (n = 7485) in the Veterans Health Administration (VA) between 2000 and 2020 using validated ICD-9/10 codes. We identified numeric values in proximity to FVC string patterns from clinical notes in the EHR. Subsequently, we performed processing steps to account for variability in note structure, related pulmonary function test (PFT) output, and values copied across notes, then assigned dates from linked administrative procedure records. NLP-derived FVC values were compared to values recorded directly from PFT equipment available on a subset of patients. RESULTS: We identified 5911 FVC values (n = 1844 patients) from PFT equipment and 15 383 values (n = 4982 patients) by NLP. Among 2610 date-matched FVC values from NLP and PFT equipment, 95.8% of values were within 5% predicted. The mean (SD) difference was 0.09% (5.9), and values strongly correlated (r = 0.94, p < 0.001), with a precision of 0.87 (95% CI 0.86, 0.88). NLP captured more patients with longitudinal FVC values (n = 3069 vs. n = 1164). Mean (SD) change in FVC %-predicted per year was similar between sources (-1.5 [30.0] NLP vs. -0.9 [16.6] PFT equipment; standardized response mean = 0.05 for both). CONCLUSIONS: NLP of EHR notes increases the capture of accurate, longitudinal FVC values by three-fold over PFT equipment. Use of this NLP tool can facilitate pharmacoepidemiologic research in RA-ILD and other lung diseases by capturing this critical measure of disease severity.

The impact of disease severity measures on survival in U.S. veterans with rheumatoid arthritis-associated interstitial lung disease.

OBJECTIVES: To determine whether RA and interstitial lung disease (ILD) severity measures are associated with survival in patients with RA-ILD. METHODS: We studied US veterans with RA-ILD participating in a multicentre, prospective RA cohort study. RA disease activity (28-joint DAS [DAS28-ESR]) and functional status (multidimensional HAQ [MDHAQ]) were collected longitudinally while pulmonary function tests (forced vital capacity [FVC], diffusing capacity for carbon monoxide) were obtained from medical records. Vital status and cause of death were determined from the National Death Index and administrative data. Predictors of death were assessed using multivariable Cox regression models adjusting for age, sex, smoking status, ILD duration, comorbidity burden and medications. RESULTS: We followed 227 RA-ILD participants (93% male and mean age of 69 years) over 1073 person-years. Median survival after RA-ILD diagnosis was 8.5 years. Respiratory diseases (28%) were the leading cause of death, with ILD accounting for 58% of respiratory deaths. Time-varying DAS28-ESR (adjusted hazard ratio [aHR] 1.21; 95% CI: 1.03, 1.41) and MDHAQ (aHR 1.85; 95% CI: 1.29, 2.65) were separately associated with mortality independent of FVC and other confounders. Modelled together, the presence of either uncontrolled disease activity (moderate/high DAS28-ESR) or FVC impairment (<80% predicted) was significantly associated with mortality risk. Those with a combination of moderate/high disease activity and FVC <80% predicted had the highest risk of death (aHR 4.43; 95% CI: 1.70, 11.55). CONCLUSION: Both RA and ILD disease severity measures are independent predictors of survival in RA-ILD. These findings demonstrate the prognostic value of monitoring the systemic features of RA-ILD.

Burden and trajectory of multimorbidity in rheumatoid arthritis: a matched cohort study from 2006 to 2015.

OBJECTIVES: To compare the onset and trajectory of multimorbidity between individuals with and without rheumatoid arthritis (RA). METHODS: A matched, retrospective cohort study was completed in a large, US commercial insurance database (MarketScan) from 2006 to 2015. Using validated algorithms, patients with RA (overall and incident) were age-matched and sex-matched to patients without RA. Diagnostic codes for 44 preidentified chronic conditions were selected to determine the presence (≥2 conditions) and burden (count) of multimorbidity. Cross-sectional comparisons were completed using the overall RA cohort and conditional logistic and negative binomial regression models. Trajectories of multimorbidity were assessed within the incident RA subcohort using generalised estimating equations. RESULTS: The overall cohort (n=277 782) and incident subcohort (n=61 124) were female predominant (76.5%, 74.1%) with a mean age of 55.6 years and 54.5 years, respectively. The cross-sectional prevalence (OR 2.29, 95% CI 2.25 to 2.34) and burden (ratio of conditions 1.68, 95% CI 1.66 to 1.70) of multimorbidity were significantly higher in RA than non-RA in the overall cohort. Within the incident RA cohort, patients with RA had more chronic conditions than non-RA (ß 1.13, 95% CI 1.10 to 1.17), and the rate of accruing chronic conditions was significantly higher in RA compared with non-RA (RA × follow-up year, ß 0.21, 95% CI 0.20 to 0.21, p<0.001). Results were similar when including the pre-RA period and in several sensitivity analyses. CONCLUSIONS: Multimorbidity is highly prevalent in RA and progresses more rapidly in patients with RA than in patients without RA during and immediately following RA onset. Therefore, multimorbidity should be aggressively identified and targeted early in the RA disease course.

Investigating changes in disease activity as a mediator of cardiovascular risk reduction with methotrexate use in rheumatoid arthritis.

OBJECTIVE: Examine the association of methotrexate (MTX) use with cardiovascular disease (CVD) in rheumatoid arthritis (RA) using marginal structural models (MSM) and determine if CVD risk is mediated through modification of disease activity. METHODS: We identified incident CVD events (coronary artery disease (CAD), stroke, heart failure (HF) hospitalisation, CVD death) within a multicentre, prospective cohort of US Veterans with RA. A 28-joint Disease Activity Score with C-reactive protein (DAS28-CRP) was collected at regular visits and medication exposures were determined by linking to pharmacy dispensing data. MSMs were used to estimate the treatment effect of MTX on risk of incident CVD, accounting for time-varying confounders between receiving MTX and CVD events. A mediation analysis was performed to estimate the indirect effects of methotrexate on CVD risk through modification of RA disease activity. RESULTS: Among 2044 RA patients (90% male, mean age 63.9 years, baseline DAS28-CRP 3.6), there were 378 incident CVD events. Using MSM, MTX use was associated with a 24% reduced risk of composite CVD events (HR 0.76, 95% CI 0.58 to 0.99) including a 57% reduction in HF hospitalisations (HR 0.43, 95% CI 0.24 to 0.77). Individual associations with CAD, stroke and CVD death were not statistically significant. In mediation analyses, there was no evidence of indirect effects of MTX on CVD risk through disease activity modification (HR 1.03, 95% CI 0.80 to 1.32). CONCLUSIONS: MTX use in RA was associated with a reduced risk of CVD events, particularly HF-related hospitalisations. These associations were not mediated through reductions in RA disease activity, suggesting alternative MTX-related mechanisms may modify CVD risk in this population.

Associations between Adiponectin and the Development of Diabetes in Rheumatoid Arthritis.

PURPOSE: We evaluated associations between adiponectin and the risk of diabetes among patients with rheumatoid arthritis (RA), a systemic inflammatory disease associated with metabolic disturbance. METHODS: This prospective cohort study included adults with RA from the Veteran's Affairs Rheumatoid Arthritis Registry. Adiponectin and inflammatory cytokines/chemokines were measured at enrollment on stored serum samples. Adiponectin levels were categorized and clinical variables were described across categories (<10 µg/mL; 10-40 µg/mL; > 40 µg/mL. Multivariable Cox proportional hazard models evaluated associations between adiponectin and incident diabetes adjusting for age, sex, race, smoking status, body mass index (BMI), disease-modifying therapy use, calendar year, and comorbidity. Testing for modification of effect in the context of elevated cytokines/chemokines was performed. RESULTS: Among 2595 patients included in the analysis, those with adiponectin levels >40 µg/mL (N = 379; 15%) were older and had lower BMI. There were 125 new cases of diabetes among 1,689 patients without prevalent disease at enrollment. There was an inverse association between adiponectin and incident diabetes, however, the association was positive among patients with adiponectin levels >40 µg/mL. Patients with levels >40 µg/mL were at higher risk compared to those with levels 10-40 µg/mL [HR: 1.70 (1.34,2.16) p < 0.001]. Those with adiponectin levels >40 µg/mL had significantly higher levels of inflammatory cytokines with evidence of a modified effect of adiponectin on diabetes risk in the setting of inflammation. CONCLUSIONS: The relationship between adiponectin and incident diabetes risk is U-shaped in RA. Patients with very high adiponectin levels have greater systemic inflammation and an altered relationship between adiponectin and diabetes risk.

The Risk of Lung Cancer in Rheumatoid Arthritis and Rheumatoid Arthritis-Associated Interstitial Lung Disease.

OBJECTIVE: To evaluate lung cancer risk in rheumatoid arthritis (RA) and RA-interstitial lung disease (ILD). METHODS: We performed a retrospective, matched cohort study of RA and RA-ILD within the Veterans Health Administration (VA) between 2000-2019. RA and RA-ILD patients were identified with validated administrative-based algorithms then matched (up to 1:10) on age, gender, and VA enrollment year to individuals without RA. Lung cancers were identified from a VA oncology database and the National Death Index. Conditional Cox regression models assessed lung cancer risk adjusting for race, ethnicity, smoking status, Agent Orange exposure, and comorbidity burden among matched individuals. Several sensitivity analyses were performed. RESULTS: We matched 72,795 RA to 633,937 non-RA persons (mean age 63 years; 88% male). Over 4,481,323 patient-years, 17,099 incident lung cancers occurred. RA was independently associated with an increased lung cancer risk (aHR 1.58 [1.52, 1.64]), which persisted in never-smokers (aHR 1.65 [1.22, 2.24]) and incident RA (aHR 1.54 [1.44, 1.65]). Compared to non-RA controls, prevalent RA-ILD (n=757) was more strongly associated with lung cancer risk (aHR 3.25 [2.13, 4.95]) than RA without ILD (aHR 1.57 (1.51, 1.64]). Analyses of both prevalent and incident RA-ILD produced similar results (RA-ILD vs. non-RA aHR 2.88 [2.45, 3.40]). CONCLUSION: RA was associated with a >50% increased risk of lung cancer, and RA-ILD represented a particularly high-risk group with an approximate 3-fold increased risk. Increased lung cancer surveillance in RA, and especially RA-ILD, may be a useful strategy for reducing the burden posed by the leading cause of cancer death.

Antibodies to Malondialdehyde-Acetaldehyde Adduct Are Associated With Prevalent and Incident Rheumatoid Arthritis-Associated Interstitial Lung Disease in US Veterans.

OBJECTIVE: The objective of this study is to determine the associations of protein-specific anti-malondialdehyde-acetaldehyde (MAA) antibodies with prevalent and incident rheumatoid arthritis-interstitial lung disease (RA-ILD). METHODS: Within a multicenter, prospective cohort of US veterans with RA, RA-ILD was validated by medical record review of clinical diagnoses, chest imaging, and pathology. Serum antibodies to MAA-albumin, MAA-collagen, MAA-fibrinogen, and MAA-vimentin (IgA, IgM, and IgG) were measured by a standardized enzyme-linked immunosorbent assay. Associations of anti-MAA antibodies with prevalent and incident RA-ILD were assessed using multivariable regression models adjusting for established RA-ILD risk factors. RESULTS: Among 2,739 participants with RA (88% male, mean age of 64 years), there were 114 with prevalent and 136 with incident RA-ILD (average time to diagnosis: 6.6 years). Higher IgM anti-MAA-collagen (per 1 SD: adjusted odds ratio [aOR] 1.28, 95% confidence interval [CI] 1.02-1.61), IgA anti-MAA-fibrinogen (aOR 1.48, 95% CI 1.14-1.92), and IgA (aOR 1.78, 95% CI 1.34-2.37) and IgG (aOR 1.48, 95% CI 1.14-1.92) anti-MAA-vimentin antibodies were associated with prevalent RA-ILD. In incident analyses, higher IgA (per one SD: adjusted hazards ratio [aHR] 1.40, 95% CI 1.11-1.76) and IgM (aHR 1.29, 95% CI 1.04-1.60) anti-MAA-albumin antibody concentrations were associated with increased ILD risk. Participants with IgA (aHR 2.13, 95% CI 1.16-3.90) or IgM (aHR 1.98, 95% CI 1.08-3.64) anti-MAA-albumin antibody concentrations in the highest quartile had an approximately two-fold increased risk of incident RA-ILD. Across all isotypes, anti-MAA-fibrinogen, anti-MAA-collagen, and anti-MAA-vimentin antibodies were not significantly associated with incident RA-ILD. CONCLUSION: Protein-specific anti-MAA antibodies to collagen, fibrinogen, and vimentin were associated with prevalent RA-ILD. IgA and IgM anti-MAA-albumin antibodies were associated with a higher risk of incident RA-ILD. These findings suggest that MAA modifications and resultant immune responses may contribute to RA-ILD pathogenesis.

Plasma Matrix Metalloproteinase Concentrations and Risk of Interstitial Lung Disease in a Prospective Rheumatoid Arthritis Cohort.

OBJECTIVE: To evaluate the associations of plasma matrix metalloproteinases (MMPs) with prevalent and incident interstitial lung disease (ILD) in people with rheumatoid arthritis (RA). METHODS: Within a multicenter, prospective cohort of US veterans with RA, we performed a cross-sectional study of prevalent ILD and cohort study of incident ILD. ILD diagnoses were validated by medical record review of provider diagnoses and chest imaging and/or pathology reports. MMP-1, 3, 7, and 9 concentrations were measured in plasma samples, then standardized and categorized into quartiles. The associations of MMPs with prevalent and incident ILD were assessed with logistic (prevalent) and Cox (incident) regression models adjusted for RA-ILD risk factors. RESULTS: Among 2,312 participants (88.9% male; mean age 63.8 years), 96 had prevalent ILD. Incident ILD developed in 130 participants over 17,378 person-years of follow-up (crude incidence rate 7.5/1,000 person-years). Participants with the highest quartile of MMP-7 concentrations had a nearly four-fold increased odds of prevalent ILD (adjusted odds ratio 3.78 [95% confidence interval (95% CI) 1.86-7.65]) and over two-fold increased risk of incident ILD (adjusted hazard ratio 2.33 [95% CI 1.35-4.02]). Higher MMP-9 concentrations were also associated with prevalent and incident ILD, as well as negatively correlated with forced vital capacity among those with prevalent ILD (r = -0.30, P = 0.005). CONCLUSION: MMP-7 and MMP-9 were strongly associated with both prevalent and incident ILD in this large, multicenter RA cohort after adjustment for other RA-ILD risk factors. These population-level findings further support a potential pathogenic role for MMPs in RA-ILD and suggest that their measurement could facilitate RA-ILD risk stratification.

Cause-Specific Mortality in Patients With Gout in the US Veterans Health Administration: A Matched Cohort Study.

OBJECTIVE: To compare all-cause and cause-specific mortality risk between patients with gout and patients without gout in the Veteran's Health Administration (VHA). METHODS: We performed a matched cohort study, identifying patients with gout in the VHA from January 1999 to September 2015 based on the presence of ≥2 International Classification of Diseases, Ninth Revision codes for gout (274.X). Gout patients were matched up to 1:10 on birth year, sex, and year of VHA enrollment with patients without gout and followed until death or end of study (December 2017). Cause of death was obtained from the National Death Index. Associations of gout with all-cause and cause-specific mortality were examined using multivariable Cox regression. RESULTS: Gout (n = 559,243) and matched non-gout controls (n = 5,428,760) had a mean age of 67 years and were 99% male. There were 246,291 deaths over 4,250,371 patient-years in gout patients and 2,000,000 deaths over 40,441,353 patient-years of follow-up in controls. After matching, gout patients had an increased risk of death (hazard ratio [HR] 1.09 [95% confidence interval (95% CI) 1.08-1.09]), which was no longer present after adjusting for comorbidities (HR 0.98 [95% CI 0.97-0.98]). The strongest association of gout with cause-specific mortality was observed with genitourinary conditions (HR 1.50 [95% CI 1.47-1.54]). Gout patients were at lower risk of death related to neurologic (e.g., Alzheimer's disease and Parkinson's disease) (HR 0.63 [95% CI 0.62-0.65]) and mental health (HR 0.66 [95% CI 0.65-0.68]) conditions. CONCLUSION: A higher risk of death among gout patients in the VHA was related to comorbidity burden. While deaths attributable to neurologic and mental health conditions were less frequent among gout patients, genitourinary conditions were the most overrepresented causes of death.

Updating and Validating the Rheumatic Disease Comorbidity Index to Incorporate ICD-10-CM Diagnostic Codes.

OBJECTIVE: To update and validate the Rheumatic Disease Comorbidity Index (RDCI) utilizing International Classification of Diseases, Tenth Revision, Clinical Modification (ICD-10-CM) codes. METHODS: We defined ICD-9-CM (n = 1,068) and ICD-10-CM (n = 1,425) era cohorts (n = 862 in both) spanning the ICD-9-CM to ICD-10-CM transition in a multicenter, prospective rheumatoid arthritis registry. Information regarding comorbidities was collected from linked administrative data over 2-year assessment periods. An ICD-10-CM code list was generated from crosswalks and clinical expertise. ICD-9- and ICD-10-derived RDCI scores were compared using intraclass correlation coefficients (ICC). The predictive ability of the RDCI for functional status and death during follow-up was assessed using multivariable regression models and goodness-of-fit statistics (Akaike's information criterion [AIC] and quasi information criterion [QIC]) in both cohorts. RESULTS: Mean ± SD RDCI scores were 2.93 ± 1.72 in the ICD-9-CM cohort and 2.92 ± 1.74 in the ICD-10-CM cohort. RDCI scores had substantial agreement in individuals who were in both cohorts (ICC 0.71 [95% confidence interval 0.68-0.74]). Prevalence of comorbidities was similar between cohorts with absolute differences <6%. Higher RDCI scores were associated with a greater risk of death and poorer functional status during follow-up in both cohorts. Similarly, in both cohorts, models including the RDCI score had the lowest QIC (functional status) and AIC (death) values, indicating better model performance. CONCLUSION: The newly proposed ICD-10-CM codes for the RDCI-generated comparable RDCI scores to those derived from ICD-9-CM codes and are highly predictive of functional status and death. The proposed ICD-10-CM codes for the RDCI can be used in rheumatic disease outcomes research spanning the ICD-10-CM era.

Identification of Multimorbidity Patterns in Rheumatoid Arthritis Through Machine Learning.

OBJECTIVE: Recognizing that the interrelationships between chronic conditions that complicate rheumatoid arthritis (RA) are poorly understood, we aimed to identify patterns of multimorbidity and to define their prevalence in RA through machine learning. METHODS: We constructed RA and age- and sex-matched (1:1) non-RA cohorts within a large commercial insurance database (MarketScan) and the Veterans Health Administration (VHA). Chronic conditions (n = 44) were identified from diagnosis codes from outpatient and inpatient encounters. Exploratory factor analysis was performed separately in both databases, stratified by RA diagnosis and sex, to identify multimorbidity patterns. The association of RA with different multimorbidity patterns was determined using conditional logistic regression. RESULTS: We studied 226,850 patients in MarketScan (76% female) and 120,780 patients in the VHA (89% male). The primary multimorbidity patterns identified were characterized by the presence of cardiopulmonary, cardiometabolic, and mental health and chronic pain disorders. Multimorbidity patterns were similar between RA and non-RA patients, female and male patients, and patients in MarketScan and the VHA. RA patients had higher odds of each multimorbidity pattern (odds ratios [ORs] 1.17-2.96), with mental health and chronic pain disorders being the multimorbidity pattern most strongly associated with RA (ORs 2.07-2.96). CONCLUSION: Cardiopulmonary, cardiometabolic, and mental health and chronic pain disorders represent predominant multimorbidity patterns, each of which is overrepresented in RA. The identification of multimorbidity patterns occurring more frequently in RA is an important first step in progressing toward a holistic approach to RA management and warrants assessment of their clinical and predictive utility.

Risk of Prostate Cancer in US Veterans With Rheumatoid Arthritis.

OBJECTIVE: Patients with rheumatoid arthritis (RA) have an increased risk of select cancers, including lymphoma and lung cancer. Whether RA influences prostate cancer risk is uncertain. We aimed to determine the risk of prostate cancer in patients with RA compared to patients without RA in the Veterans Health Administration (VA). METHODS: We performed a matched (up to 1:5) cohort study of male patients with and without RA in the VA from 2000 to 2018. RA status, as well as covariates, were obtained from national VA databases. Prostate cancer was identified through linked VA cancer databases and the National Death Index. Multivariable Cox models compared prostate cancer risk between patients with RA and patients without RA, including models that accounted for retention in the VA system. RESULTS: We included 56,514 veterans with RA and 227,284 veterans without RA. During 2,337,104 patient-years of follow-up, 6,550 prostate cancers occurred. Prostate cancer incidence (per 1,000 patient-years) was 3.50 (95% confidence interval [95% CI] 3.32-3.69) in patients with RA and 2.66 (95% CI 2.58-2.73) in patients without RA. After accounting for confounders and censoring for attrition of VA health care, RA was modestly associated with a higher prostate cancer risk (adjusted HR [HRadj ] 1.12 [95% CI 1.04-1.20]). There was no association between RA and prostate cancer mortality (HRadj 0.92 [95% CI 0.73-1.16]). CONCLUSION: RA was associated with a modestly increased risk of prostate cancer, but not prostate cancer mortality, after accounting for relevant confounders and several potential sources of bias. However, even minimal unmeasured confounding could explain these findings.

A Narrowing Mortality Gap: Temporal Trends of Cause-Specific Mortality in a National Matched Cohort Study in US Veterans With Rheumatoid Arthritis.

OBJECTIVE: To examine temporal trends in all-cause and cause-specific mortality in patients with rheumatoid arthritis (RA) in the Veterans Health Administration (VHA). METHODS: We conducted a matched cohort study in the VHA from January 1, 2000 to December 31, 2017. Incident RA patients were matched up to 1:10 on age, sex, and VHA enrollment year to non-RA patients, then followed until death or end of study period. Cause of death was obtained from the National Death Index. Multivariable Cox regression models stratified by RA diagnosis years were used to examine trends in RA-related risk of all-cause and cause-specific mortality. RESULTS: Among 29,779 incident RA patients (matched to 245,226 non-RA patients), 9,565 deaths occurred. RA patients were at increased risk of all-cause (adjusted hazard ratio [HRadj ] 1.23 [95% confidence interval (95% CI) 1.20-1.26]), cardiovascular (HRadj 1.19 [95% CI 1.14-1.23]), cancer (HRadj 1.19 [95% CI 1.14-1.24]), respiratory (HRadj 1.46 [95% CI 1.38-1.55]), and infection-related mortality (HRadj 1.59 [95% CI 1.41-1.80]). Interstitial lung disease was the cause of death most strongly associated with RA (HRadj 3.39 [95% CI 2.88-3.99]). Nearly 70% of excess deaths in RA were attributable to cardiopulmonary disease. All-cause mortality risk related to RA was lower among those diagnosed during 2012-2017 (HRadj 1.10 [95% CI 1.05-1.15]) compared to 2000-2005 (HRadj 1.31 [95% CI 1.26-1.36]), but still higher than for non-RA controls (P < 0.001). Cause-specific mortality trends were similar. CONCLUSION: Excess RA-related mortality was driven by cardiovascular, cancer, respiratory, and infectious causes, particularly cardiopulmonary diseases. Although our findings support that RA-related mortality risk is decreasing over time, a mortality gap remains for all-cause and cause-specific mortality in RA.

Multimorbidity Patterns and Rheumatoid Arthritis Disease Outcomes: Findings From a Multicenter, Prospective Cohort.

OBJECTIVE: To determine whether unique multimorbidity patterns are associated with long-term rheumatoid arthritis (RA) disease severity. METHODS: We conducted a cohort study within the Veterans Affairs Rheumatoid Arthritis registry. We applied previously derived multimorbidity patterns based on the presence of diagnostic codes for relevant conditions prior to enrollment using linked administrative data. Disease activity and functional status were assessed longitudinally up to 5 years after enrollment. The association of multimorbidity patterns with disease activity and functional status were assessed using generalized estimating equations models adjusting for relevant confounders. RESULTS: We studied 2,956 participants, of which 88.2% were male, 76.9% reported white race, and 79.3% had a smoking history. Mental health and substance abuse (ß 0.12 [95% confidence interval {CI} 0.00, 0.23]), cardiovascular (ß 0.25 [95% CI 0.12, 0.38]), and chronic pain (ß 0.21 [95% CI 0.11, 0.31]) multimorbidity were associated with higher Disease Activity Score in 28 joints (DAS28) scores. Mental health and substance abuse (ß 0.09 [0.03, 0.15]), cardiovascular (ß 0.11 [95% CI 0.04, 0.17]), and chronic pain multimorbidity (ß 0.15 [95% CI 0.10, 0.20]) were also associated with higher Multidimensional Health Assessment Questionnaire (MDHAQ) scores. The metabolic pattern of multimorbidity was not associated with DAS28 or MDHAQ. The number of multimorbidity patterns present was highly associated with DAS28 and MDHAQ (P trend < 0.001), and patients with all four multimorbidity patterns had the highest DAS28 (ß 0.59 [95% CI 0.36, 0.83]) and MDHAQ (ß 0.27 [95% CI 0.16, 0.39]) scores. CONCLUSION: Mental health and substance abuse, chronic pain, and cardiovascular multimorbidity patterns are associated with increased RA disease activity and poorer functional status. Identifying and addressing these multimorbidity patterns may facilitate achieving RA treatment targets.

Aortic Stenosis Risk in Rheumatoid Arthritis.

Importance: Although an increased risk of ischemic cardiovascular disease has been associated with rheumatoid arthritis (RA), the risk of aortic stenosis (AS) is unknown. Objective: To examine the risk of incident AS, aortic valve intervention, AS-related death, and risk factors for AS development in patients with RA. Design, Setting, and Participants: This cohort study linked data from the Veterans Health Administration (VHA) and Centers for Medicare & Medicaid Services from 2000 to 2019. Patients with RA were matched by age, sex, and VHA enrollment year with up to 10 patients without RA. The cohort was followed until incident AS, aortic valve intervention, or death. Data were analyzed from August 23, 2022, to March 3, 2023. Exposures: the primary exposure was the presence of RA, defined using validated RA algorithms. Main Outcomes and Measures: Aortic stenosis was defined as a composite of inpatient or outpatient diagnoses, surgical or transcatheter aortic valve replacement, or AS-related death using diagnostic and procedural codes. Risk of AS development was assessed with multivariable Cox proportional hazards models adjusted for race, ethnicity, smoking status, body mass index, rurality, comorbidities, and health care use. Results: The cohort included 73 070 patients with RA (64 008 [87.6%] males; mean [SD] age, 63.0 [11.9] years) matched with 639 268 patients without RA (554 182 [86.7%] males; mean [SD] age, 61.9 [11.7] years) and 16 109 composite AS outcomes that occurred over 6 223 150 person-years. The AS incidence rate was 3.97 (95% CI, 3.81-4.13) per 1000 person-years in patients with RA and 2.45 (95% CI, 2.41-2.49) per 1000 person-years in the control patients (absolute difference, 1.52 per 1000 person-years). Rheumatoid arthritis was associated with an increased risk of composite AS (adjusted hazard ratio [AHR], 1.48; 95% CI, 1.41-1.55), aortic valve intervention (AHR, 1.34; 95% CI, 1.22-1.48), and AS-related death (AHR, 1.26; 95% CI, 1.04-1.54). Conclusions and Relevance: In this cohort study, RA was associated with a higher risk of developing AS and the subsequent risks of undergoing aortic valve intervention and suffering from AS-related death. Future studies are needed to confirm whether valvular heart disease, specifically AS, may be an overlooked cardiovascular disease complication in RA.

Lowering Expectations: Glucocorticoid Tapering Among Veterans With Rheumatoid Arthritis Achieving Low Disease Activity on Stable Biologic Therapy.

OBJECTIVE: In the Steroid EliMination In Rheumatoid Arthritis (SEMIRA) trial, 65% of patients with rheumatoid arthritis (RA) in low disease activity (LDA) on stable biologic therapy successfully tapered glucocorticoids. We aimed to evaluate real-world rates of glucocorticoid tapering among similar patients in the Veterans Affairs Rheumatoid Arthritis registry. METHODS: Within a multicenter, prospective RA cohort, we used registry data and linked pharmacy claims from 2003 to 2021 to identify chronic prednisone users achieving LDA after initiating a new biologic or targeted synthetic disease-modifying antirheumatic drug (b/tsDMARD). We defined the index date as first LDA occurring 60 to 180 days after b/tsDMARD initiation. The primary outcome of successful tapering, assessed at day 180 after LDA, required a 30-day averaged prednisone dose both less than or equal to 5mg/day and at least 50% lower than at the index date. The secondary outcome was discontinuation, defined as a prednisone dose of 0 mg/day at days 180 through 210. We used univariate statistics to compare patient characteristics by fulfillment of the primary outcome. RESULTS: We evaluated 100 b/tsDMARD courses among 95 patients. Fifty-four courses resulted in successful tapering; 33 resulted in discontinuation. Positive rheumatoid factor, higher erythrocyte sedimentation rate, more background DMARDs, shorter time from b/tsDMARD initiation to LDA, and higher glucocorticoid dose 30 days before LDA were associated with greater likelihood of successful tapering. CONCLUSION: In a real-world RA cohort of chronic glucocorticoid users in LDA, half successfully tapered and a third discontinued prednisone within 6 months of initiating a new b/tsDMARD. Claims-based algorithms of glucocorticoid tapering and discontinuation may be useful to evaluate predictors of tapering in administrative data sets.

Enhancing the identification of rheumatoid arthritis-associated interstitial lung disease through text mining of chest computerized tomography reports.

OBJECTIVES: Algorithms have been developed to identify rheumatoid arthritis-interstitial lung disease (RA-ILD) in administrative data with positive predictive values (PPVs) between 70 and 80%. We hypothesized that including ILD-related terms identified within chest computed tomography (CT) reports through text mining would improve the PPV of these algorithms in this cross-sectional study. METHODS: We identified a derivation cohort of possible RA-ILD cases (n = 114) using electronic health record data from a large academic medical center and performed medical record review to validate diagnoses (reference standard). ILD-related terms (e.g., ground glass, honeycomb) were identified in chest CT reports by natural language processing. Administrative algorithms including diagnostic and procedural codes as well as specialty were applied to the cohort both with and without the requirement for ILD-related terms from CT reports. We subsequently analyzed similar algorithms in an external validation cohort of 536 participants with RA. RESULTS: The addition of ILD-related terms to RA-ILD administrative algorithms increased the PPV in both the derivation (improvement ranging from 3.6 to 11.7%) and validation cohorts (improvement 6.0 to 21.1%). This increase was greatest for less stringent algorithms. Administrative algorithms including ILD-related terms from CT reports exceeded a PPV of 90% (maximum 94.6% derivation cohort). Increases in PPV were accompanied by a decline in sensitivity (validation cohort -3.9 to -19.5%). CONCLUSIONS: The addition of ILD-related terms identified by text mining from chest CT reports led to improvements in the PPV of RA-ILD algorithms. With high PPVs, use of these algorithms in large data sets could facilitate epidemiologic and comparative effectiveness research in RA-ILD.

Comparison of Rates of Lower Extremity Amputation in Patients With and Without Gout in the US Department of Veterans Affairs Health System.

Importance: Cardiometabolic and other risk factors could render patients with gout more likely to undergo lower extremity amputation (LEA). Objective: To examine the rate of and factors associated with LEA in patients with gout. Design, Setting, and Participants: In this matched cohort study using national administrative data, multivariable Cox proportional hazards regression models were used to examine the associations of gout with LEA. In analyses limited to patients with gout, attributes of serum urate control and treatment with urate-lowering therapy were examined as factors associated with LEA. This study included patients who used US Department of Veterans Affairs services from January 1, 2000, to July 31, 2015. Patients with gout were identified using diagnostic codes and matched with up to 10 controls by age, sex, and year of benefit enrollment. Data analysis was performed from January 26, 2021, to September 3, 2021. Exposures: Gout classification served as the primary independent variable of interest. In analyses limited to patients with gout, factors associated with serum urate control and urate-lowering therapy were examined. Main Outcomes and Measures: Overall LEA, as well as toe, transmetatarsal, below-the-knee, and above-the-knee amputation. Results: This cohort study included 5 924 918 patients, 556 521 with gout (mean [SD] age, 67 [12] years; 550 963 (99.0%) male; 88 853 [16.0%] Black non-Hispanic; 16 981 [4.3%] Hispanic/Latinx; 345 818 [62.1%] White non-Hispanic; 80 929 [14.5%] with race and ethnicity data missing; and 23 940 [4.3%] classified as other) and 5 368 397 without gout (mean [SD] age, 67 [12] years; 5 314 344 [99.0%] male; 558 464 [10.4%] Black non-Hispanic; 204 291 [3.0%] Hispanic/Latinx; 3 188 504 [59.4%] White non-Hispanic; 1 257 739 [23.4%)] with race and ethnicity data missing; and 159 399 [3.0%] classified as other). Compared with patients without gout, patients with gout were more likely to undergo amputation, an increased rate that remained after adjustment (adjusted hazard ratio, 1.20; 95% CI, 1.16-1.24) and was highest for below-the-knee amputation (adjusted hazard ratio, 1.59; 95% CI, 1.39-1.81). In those with gout, poor serum urate control (mean >7 mg/dL during the preceding year) was associated with a 25% to 37% increase in the rate of amputation. In contrast, treatment with urate-lowering therapy was not associated with the LEA rate. Conclusions and Relevance: In this matched cohort study, patients with gout were more likely to undergo LEA. This increase was independent of other comorbidities that have been associated with amputation, including diabetes and peripheral vascular disease. Serum urate control was independently associated with the LEA rate, suggesting the possibility that lower extremity amputation may be preventable in some patients.

Associations between shortened telomeres and rheumatoid arthritis-associated interstitial lung disease among U.S. Veterans.

BACKGROUND: Shortened telomeres are associated with several different subtypes of interstitial lung disease (ILD), although studies of telomere length and ILD in rheumatoid arthritis (RA) are lacking. METHODS: Within the Veterans Affairs Rheumatoid Arthritis (VARA) registry, we performed cross-sectional and case-control studies of prevalent and incident ILD, respectively. We randomly selected a subset of RA patients with ILD and individually matched them to RA patients without ILD according to age, sex, and VARA enrollment date. Telomere length was measured on peripheral blood leukocytes collected at registry enrollment using quantitative PCR (T/S ratio). Short telomeres were defined as a T/S ratio in the lowest 10th percentile of the cohort. RESULTS: Our cross-sectional study cohort was comprised of 54 RA-ILD patients and 92 RA-non-ILD patients. T/S ratios significantly differed between patients with and without prevalent ILD (1.56 [IQR 1.30, 1.78] vs. 1.96 [IQR 1.65, 2.27], p < 0.001). Similarly, prevalence of ILD was significantly higher in patients with short vs. normal-length telomeres (73.3% vs. 32.8%, p = 0.002). Short telomeres were independently associated with an increased odds of prevalent ILD compared to normal-length telomeres (adjusted OR 6.60, 95% CI 1.78-24.51, p = 0.005). In our case-control analysis, comprised of 22 incident RA-ILD cases and 36 RA-non-ILD controls, short telomeres were not associated with incident RA-ILD (adjusted OR 0.90, 95% CI 0.06-13.4, p = 0.94). CONCLUSION: Short telomeres were strongly associated with prevalent but not incident ILD among patients with RA. Additional studies are needed to better understand telomere length dynamics among RA patients with and without ILD.