Impact of RA treatment strategies on lipids and vascular inflammation in rheumatoid arthritis: a secondary analysis of the TARGET randomized active comparator trial.

BACKGROUND: Treatments for rheumatoid arthritis (RA) are associated with complex changes in lipids and lipoproteins that may impact cardiovascular (CV) risk. The objective of this study was to examine lipid and lipoprotein changes associated with two common RA treatment strategies, triple therapy or tumor necrosis factor inhibitor (TNFi), and association with CV risk. METHODS: In this secondary data analysis of the TARGET trial, methotrexate (MTX) inadequate responders with RA were randomized to either add sulfasalazine and hydroxychloroquine (triple therapy), or TNFi for 24-weeks. The primary trial outcome was the change in arterial inflammation measured in the carotid arteries or aorta by FDG-PET/CT at baseline and 24-weeks; this change was described as the target-to-background ratio (TBR) in the most diseased segment (MDS). Routine lipids and advanced lipoproteins were measured at baseline and 24-weeks; subjects on statin therapy at baseline were excluded. Comparisons between baseline and follow-up lipid measurements were performed within and across treatment arms, as well as change in lipids and change in MDS-TBR. RESULTS: We studied 122 participants, 61 in each treatment arm, with median age 57 years, 76% female, and 1.5 year median RA disease duration. When comparing treatment arms, triple therapy had on average a larger reduction in triglycerides (15.9 mg/dL, p = 0.01), total cholesterol to HDL-C ratio (0.29, p-value = 0.01), and LDL particle number (111.2, p = 0.02) compared to TNFi. TNFi had on average a larger increase in HDL particle number (1.6umol/L, p = 0.006). We observed no correlation between change in lipid measurements and change in MDS-TBR within and across treatment arms. CONCLUSIONS: Both treatment strategies were associated with improved lipid profiles via changes in different lipids and lipoproteins. These effects had no correlation with change in CV risk as measured by vascular inflammation by FDG-PET/CT. TRIAL REGISTRATION: ClinicalTrials.gov ID NCT02374021.

The association of TNF inhibitor use with incident cardiovascular events in radiographic axial spondyloarthritis.

BACKGROUND: Whether tumor necrosis factor inhibitor (TNFi) use is cardioprotective among individuals with radiographic axial spondyloarthritis (r-axSpA), who have heightened cardiovascular (CV) risk, is unclear. We tested the association of TNFi use with incident CV outcomes in r-axSpA. METHODS: We identified a r-axSpA cohort within a Veterans Affairs database between 2002 and 2019 using novel phenotyping methods and secondarily using ICD codes. TNFi use was assessed as a time-varying exposure using pharmacy dispense records. The primary outcome was incident CV disease identified using ICD codes for coronary artery disease, myocardial infarction or stroke. We fit Cox models with inverse probability weights to estimate the risk of each outcome with TNFi use versus non-use. Analyses were performed in the overall cohort, and separately in two periods (2002-2010, 2011-2019) to account for secular trends. RESULTS: Using phenotyping we identified 26,928 individuals with an r-axSpA diagnosis (mean age 63.4 years, 94 % male); at baseline 3633 were TNFi users and 23,295 were non-users. During follow-up of a mean 3.3 ± 4.2 years, 674 (18.6 %) TNFi users had incident CVD versus 11,838 (50.8 %) non-users. In adjusted analyses, TNFi use versus non-use was associated with lower risk of incident CVD (HR 0.34, 95 % CI 0.29-0.40) in the cohort overall, and in the two time periods separately. CONCLUSION: In this r-axSpA cohort identified using phenotyping methods, TNFi use versus non-use had a lower risk of incident CVD. These findings provide reassurance regarding the CV safety of TNFi agents for r-axSpA treatment. Replication of these results in other cohorts is needed.

Interplay Between Systemic Inflammation, Myocardial Injury, and Coronary Microvascular Dysfunction in Rheumatoid Arthritis: Results From the LiiRA Study.

BACKGROUND: Coronary microvascular dysfunction as measured by myocardial flow reserve (MFR) is associated with increased cardiovascular risk in rheumatoid arthritis (RA). The objective of this study was to determine the association between reducing inflammation with MFR and other measures of cardiovascular risk. METHODS AND RESULTS: Patients with RA with active disease about to initiate a tumor necrosis factor inhibitor were enrolled (NCT02714881). All subjects underwent a cardiac perfusion positron emission tomography scan to quantify MFR at baseline before tumor necrosis factor inhibitor initiation, and after tumor necrosis factor inhibitor initiation at 24 weeks. MFR <2.5 in the absence of obstructive coronary artery disease was defined as coronary microvascular dysfunction. Blood samples at baseline and 24 weeks were measured for inflammatory markers (eg, high-sensitivity C-reactive protein [hsCRP], interleukin-1b, and high-sensitivity cardiac troponin T [hs-cTnT]). The primary outcome was mean MFR before and after tumor necrosis factor inhibitor initiation, with Δhs-cTnT as the secondary outcome. Secondary and exploratory analyses included the correlation between ΔhsCRP and other inflammatory markers with MFR and hs-cTnT. We studied 66 subjects, 82% of which were women, mean RA duration 7.4 years. The median atherosclerotic cardiovascular disease risk was 2.5%; 47% had coronary microvascular dysfunction and 23% had detectable hs-cTnT. We observed no change in mean MFR before (2.65) and after treatment (2.64, P=0.6) or hs-cTnT. A correlation was observed between a reduction in hsCRP and interleukin-1b with a reduction in hs-cTnT. CONCLUSIONS: In this RA cohort with low prevalence of cardiovascular risk factors, nearly 50% of subjects had coronary microvascular dysfunction at baseline. A reduction in inflammation was not associated with improved MFR. However, a modest reduction in interleukin-1b and no other inflammatory pathways was correlated with a reduction in subclinical myocardial injury. REGISTRATION: URL: https://www.clinicaltrials.gov; Unique identifier: NCT02714881.

Association of the multi-biomarker disease activity score with arterial 18-fluorodeoxyglucose uptake in rheumatoid arthritis.

OBJECTIVES: Rheumatoid arthritis (RA) and atherosclerosis share many common inflammatory pathways. We studied whether a multi-biomarker panel for RA disease activity (MBDA) would associate with changes in arterial inflammation in an interventional trial. METHODS: In the TARGET Trial, RA patients with active disease despite methotrexate were randomly assigned to the addition of either a TNF inhibitor or sulfasalazine+hydroxychloroquine (triple therapy). Baseline and 24-week follow-up 18F-fluorodeoxyglucose (FDG) positron emission tomography/computed tomography scans were assessed for change in arterial inflammation measured as the maximal arterial target-to-blood background ratio of FDG uptake in the most diseased segment of the carotid arteries or aorta (MDS-TBRmax). The MBDA test, measured at baseline and weeks 6, 18, and 24, was assessed for its association with the change in MDS-TBRmax. RESULTS: Interpretable scans were available at baseline and week 24 for n = 112 patients. The MBDA score at week 24 was significantly correlated with the change in MDR-TBRmax (Spearman's rho = 0.239; p= 0.011) and remained significantly associated after adjustment for relevant confounders. Those with low MBDA at week 24 had a statistically significant adjusted reduction in arterial inflammation of 0.35 units vs no significant reduction in those who did not achieve low MBDA. Neither DAS28-CRP nor CRP predicted change in arterial inflammation. The MBDA component with the strongest association with change in arterial inflammation was serum amyloid A (SAA). CONCLUSIONS: Among treated RA patients, achieved MBDA predicts of changes in arterial inflammation. Achieving low MBDA at 24 weeks was associated with clinically meaningful reductions in arterial inflammation, regardless of treatment.

Biomarkers of Cardiovascular Risk in Patients With Rheumatoid Arthritis: Results From the TARGET Trial.

Cardiovascular disease remains an important comorbidity in patients with rheumatoid arthritis (RA), but traditional models do not accurately predict cardiovascular risk in patients with RA. The addition of biomarkers could improve prediction. METHODS AND RESULTS: The TARGET (Treatments Against RA and Effect on FDG PET/CT) trial assessed whether different treatment strategies in RA differentially impact cardiovascular risk as measured by the change in arterial inflammation on arterial target to background ratio on fluorodeoxyglucose positron emission tomography/computed tomography scans conducted 24 weeks apart. A group of 24 candidate biomarkers supported by prior literature was assessed at baseline and 24 weeks later. Longitudinal analyses examined the association between baseline biomarker values, measured in plasma EDTA, and the change in arterial inflammation target to background ratio. Model fit was assessed for the candidate biomarkers only, clinical variables only, and models combining both. One hundred nine patients with median (interquartile range) age 58 years (53-65 years), RA duration 1.4 years (0.5-6.6 years), and 82% women had biomarkers assessed at baseline and follow-up. Because the main trial analyses demonstrated significant target to background ratio decreases with both treatment strategies but no difference across treatment groups, we analyzed all patients together. Baseline values of serum amyloid A, C-reactive protein, soluble tumor necrosis factor receptor 1, adiponectin, YKL-40, and osteoprotegerin were associated with significant change in target to background ratio. When selected candidate biomarkers were added to the clinical variables, the adjusted R2 improved from 0.20 to 0.33 (likelihood ratio P=0.0005). CONCLUSIONS: A candidate biomarker approach identified several promising biomarkers that associate with baseline and treatment-associated changes in arterial inflammation in patients with RA. These will now be tested in an external validation cohort.

Associations Between Rheumatoid Arthritis Clinical Factors and Synovial Cell Types and States.

OBJECTIVE: Recent studies have uncovered diverse cell types and states in the rheumatoid arthritis (RA) synovium; however, limited data exist correlating these findings with patient-level clinical information. Using the largest cohort to date with clinical and multicell data, we determined associations between RA clinical factors with cell types and states in the RA synovium. METHODS: The Accelerated Medicines Partnership Rheumatoid Arthritis study recruited patients with active RA who were not receiving disease-modifying antirheumatic drugs (DMARDs) or who had an inadequate response to methotrexate (MTX) or tumor necrosis factor inhibitors. RA clinical factors were systematically collected. Biopsies were performed on an inflamed joint, and tissue were disaggregated and processed with a cellular indexing of transcriptomes and epitopes sequencing pipeline from which the following cell type percentages and cell type abundance phenotypes (CTAPs) were derived: endothelial, fibroblast, and myeloid (EFM); fibroblasts; myeloid; T and B cells; T cells and fibroblasts (TF); and T and myeloid cells. Correlations were measured between RA clinical factors, cell type percentage, and CTAPs. RESULTS: We studied 72 patients (mean age 57 years, 75% women, 83% seropositive, mean RA duration 6.6 years, mean Disease Activity Score-28 C-reactive Protein 3 [DAS28-CRP3] score 4.8). Higher DAS28-CRP3 correlated with a higher T cell percentage (P < 0.01). Those receiving MTX and not a biologic DMARD (bDMARD) had a higher percentage of B cells versus those receiving no DMARDs (P < 0.01). Most of those receiving bDMARDs were categorized as EFM (57%), whereas none were TF. No significant difference was observed across CTAPs for age, sex, RA disease duration, or DAS28-CRP3. CONCLUSION: In this comprehensive screen of clinical factors, we observed differential associations between DMARDs and cell phenotypes, suggesting that RA therapies, more than other clinical factors, may impact cell type/state in the synovium and ultimately influence response to subsequent therapies.

Semi-supervised calibration of noisy event risk (SCANER) with electronic health records.

OBJECTIVE: Electronic health records (EHR), containing detailed longitudinal clinical information on a large number of patients and covering broad patient populations, open opportunities for comprehensive predictive modeling of disease progression and treatment response. However, since EHRs were originally constructed for administrative purposes not for research, in the EHR-linked studies, it is often not feasible to capture reliable information for analytical variables, especially in the survival setting, when both accurate event status and event times are needed for model building. For example, progression-free survival (PFS), a commonly used survival outcome for cancer patients, often involves complex information embedded in free-text clinical notes and cannot be extracted reliably. Proxies of PFS time such as time to the first mention of progression in the notes are at best good approximations to the true event time. This leads to difficulty in efficiently estimating event rates for an EHR patient cohort. Estimating survival rates based on error-prone outcome definitions can lead to biased results and hamper the power in the downstream analysis. On the other hand, extracting accurate event time information via manual annotation is time and resource intensive. The objective of this study is to develop a calibrated survival rate estimator using noisy outcomes from EHR data. MATERIALS AND METHODS: In this paper, we propose a two-stage semi-supervised calibration of noisy event rate (SCANER) estimator that can effectively overcome censoring induced dependency and attains more robust performance (i.e., not sensitive to misspecification of the imputation model) by fully utilizing both a small-labeled set of gold-standard survival outcomes annotated via manual chart review and a set of proxy features automatically captured via EHR in the unlabeled set. We validate the SCANER estimator by estimating the PFS rates for a virtual cohort of lung cancer patients from one large tertiary care center and the ICU-free survival rates for COVID patients from two large tertiary care centers. RESULTS: In terms of survival rate estimates, the SCANER had very similar point estimates compared to the complete-case Kaplan Meier estimator. On the other hand, other benchmark methods for comparison, which fail to account for the induced dependency between event time and the censoring time conditioning on surrogate outcomes, produced biased results across all three case studies. In terms of standard errors, the SCANER estimator was more efficient than the KM estimator, with up to 50% efficiency gain. CONCLUSION: The SCANER estimator achieves more efficient, robust, and accurate survival rate estimates compared to existing approaches. This promising new approach can also improve the resolution (i.e., granularity of event time) by using labels conditioning on multiple surrogates, particularly among less common or poorly coded conditions.

Generate Analysis-Ready Data for Real-world Evidence: Tutorial for Harnessing Electronic Health Records With Advanced Informatic Technologies.

Although randomized controlled trials (RCTs) are the gold standard for establishing the efficacy and safety of a medical treatment, real-world evidence (RWE) generated from real-world data has been vital in postapproval monitoring and is being promoted for the regulatory process of experimental therapies. An emerging source of real-world data is electronic health records (EHRs), which contain detailed information on patient care in both structured (eg, diagnosis codes) and unstructured (eg, clinical notes and images) forms. Despite the granularity of the data available in EHRs, the critical variables required to reliably assess the relationship between a treatment and clinical outcome are challenging to extract. To address this fundamental challenge and accelerate the reliable use of EHRs for RWE, we introduce an integrated data curation and modeling pipeline consisting of 4 modules that leverage recent advances in natural language processing, computational phenotyping, and causal modeling techniques with noisy data. Module 1 consists of techniques for data harmonization. We use natural language processing to recognize clinical variables from RCT design documents and map the extracted variables to EHR features with description matching and knowledge networks. Module 2 then develops techniques for cohort construction using advanced phenotyping algorithms to both identify patients with diseases of interest and define the treatment arms. Module 3 introduces methods for variable curation, including a list of existing tools to extract baseline variables from different sources (eg, codified, free text, and medical imaging) and end points of various types (eg, death, binary, temporal, and numerical). Finally, module 4 presents validation and robust modeling methods, and we propose a strategy to create gold-standard labels for EHR variables of interest to validate data curation quality and perform subsequent causal modeling for RWE. In addition to the workflow proposed in our pipeline, we also develop a reporting guideline for RWE that covers the necessary information to facilitate transparent reporting and reproducibility of results. Moreover, our pipeline is highly data driven, enhancing study data with a rich variety of publicly available information and knowledge sources. We also showcase our pipeline and provide guidance on the deployment of relevant tools by revisiting the emulation of the Clinical Outcomes of Surgical Therapy Study Group Trial on laparoscopy-assisted colectomy versus open colectomy in patients with early-stage colon cancer. We also draw on existing literature on EHR emulation of RCTs together with our own studies with the Mass General Brigham EHR.

Shared inflammatory pathways of rheumatoid arthritis and atherosclerotic cardiovascular disease.

The association between chronic inflammation and increased risk of cardiovascular disease in rheumatoid arthritis (RA) is well established. In the general population, inflammation is an established independent risk factor for cardiovascular disease, and much interest is placed on controlling inflammation to reduce cardiovascular events. As inflammation encompasses numerous pathways, the development of targeted therapies in RA provides an opportunity to understand the downstream effect of inhibiting specific pathways on cardiovascular risk. Data from these studies can inform cardiovascular risk management in patients with RA, and in the general population. This Review focuses on pro-inflammatory pathways targeted by existing therapies in RA and with mechanistic data from the general population on cardiovascular risk. Specifically, the discussions include the IL-1, IL-6 and TNF pathways, as well as the Janus kinase (JAK)-signal transducer and activator of transcription (STAT) signalling pathway, and the role of these pathways in RA pathogenesis in the joint alongside the development of atherosclerotic cardiovascular disease. Overall, some robust data support inhibition of IL-1 and IL-6 in decreasing the risk of cardiovascular disease, with growing data supporting IL-6 inhibition in both patients with RA and the general population to reduce the risk of cardiovascular disease.

Large-scale real-world data analyses of cancer risks among patients with rheumatoid arthritis.

Rheumatoid arthritis (RA) affects 24.5 million people worldwide and has been associated with increased cancer risks. However, the extent to which the observed risks are related to the pathophysiology of rheumatoid arthritis or its treatments is unknown. Leveraging nationwide health insurance claims data with 85.97 million enrollees across 8 years, we identified 92 864 patients without cancers at the time of rheumatoid arthritis diagnoses. We matched 68 415 of these patients with participants without rheumatoid arthritis by sex, race, age and inferred health and economic status and compared their risks of developing all cancer types. By 12 months after the diagnosis of rheumatoid arthritis, rheumatoid arthritis patients were 1.21 (95% confidence interval [CI] [1.14, 1.29]) times more likely to develop any cancer compared with matched enrollees without rheumatoid arthritis. In particular, the risk of developing lymphoma is 2.08 (95% CI [1.67, 2.58]) times higher in the rheumatoid arthritis group, and the risk of developing lung cancer is 1.69 (95% CI [1.32, 2.13]) times higher. We further identified the five most commonly used drugs in treating rheumatoid arthritis, and the log-rank test showed none of them is implicated with a significantly increased cancer risk compared with rheumatoid arthritis patients without that specific drug. Our study suggested that the pathophysiology of rheumatoid arthritis, rather than its treatments, is implicated in the development of subsequent cancers. Our method is extensible to investigating the connections among drugs, diseases and comorbidities at scale.

Identifying shared genetic architecture between rheumatoid arthritis and other conditions: a phenome-wide association study with genetic risk scores.

BACKGROUND: Rheumatoid arthritis (RA) shares genetic variants with other autoimmune conditions, but existing studies test the association between RA variants with a pre-defined set of phenotypes. The objective of this study was to perform a large-scale, systemic screen to determine phenotypes that share genetic architecture with RA to inform our understanding of shared pathways. METHODS: In the UK Biobank (UKB), we constructed RA genetic risk scores (GRS) incorporating human leukocyte antigen (HLA) and non-HLA risk alleles. Phenotypes were defined using groupings of International Classification of Diseases (ICD) codes. Patients with an RA code were excluded to mitigate the possibility of associations being driven by the diagnosis or management of RA. We performed a phenome-wide association study, testing the association between the RA GRS with phenotypes using multivariate generalized estimating equations that adjusted for age, sex, and first five principal components. Statistical significance was defined using Bonferroni correction. Results were replicated in an independent cohort and replicated phenotypes were validated using medical record review of patients. FINDINGS: We studied n = 316,166 subjects from UKB without evidence of RA and screened for association between the RA GRS and n = 1317 phenotypes. In the UKB, 20 phenotypes were significantly associated with the RA GRS, of which 13 (65%) were immune mediated conditions including polymyalgia rheumatica, granulomatosis with polyangiitis (GPA), type 1 diabetes, and multiple sclerosis. We further identified a novel association in Celiac disease where the HLA and non-HLA alleles had strong associations in opposite directions. Strikingly, we observed that the non-HLA GRS was exclusively associated with greater risk of the validated conditions, suggesting shared underlying pathways outside the HLA region. INTERPRETATION: This study replicated and identified novel autoimmune phenotypes verified by medical record review that share immune pathways with RA and may inform opportunities for shared treatment targets, as well as risk assessment for conditions with a paucity of genomic data, such as GPA. FUNDING: This research was funded by the US National Institutes of Health (P30AR072577, R21AR078339, R35GM142879, T32AR007530) and the Harold and DuVal Bowen Fund.

Autoantibodies against citrullinated and native proteins and prediction of rheumatoid arthritis-associated interstitial lung disease: A nested case-control study.

Background: To identify fine specificity anti-citrullinated protein antibodies (ACPA) associated with incident rheumatoid arthritis-associated interstitial lung disease (RA-ILD). Methods: This nested case-control study within the Brigham RA Sequential Study matched incident RA-ILD cases to RA-noILD controls on time of blood collection, age, sex, RA duration, and rheumatoid factor status. A multiplex assay measured ACPA and anti-native protein antibodies from stored serum prior to RA-ILD onset. Logistic regression models calculated odds ratios (OR) with 95% confidence intervals (CI) for RA-ILD, adjusting for prospectively-collected covariates. We estimated optimism-corrected area under the curves (AUC) using internal validation. Model coefficients generated a risk score for RA-ILD. Findings: We analyzed 84 incident RA-ILD cases (mean age 67 years, 77% female, 90% White) and 233 RA-noILD controls (mean age 66 years, 80% female, 94% White). We identified six fine specificity antibodies that were associated with RA-ILD. The antibody isotypes and targeted proteins were: IgA2 to citrullinated histone 4 (OR 0.08 per log-transformed unit, 95% CI 0.03-0.22), IgA2 to citrullinated histone 2A (OR 4.03, 95% CI 2.03-8.00), IgG to cyclic citrullinated filaggrin (OR 3.47, 95% CI 1.71-7.01), IgA2 to native cyclic histone 2A (OR 5.52, 95% CI 2.38-12.78), IgA2 to native histone 2A (OR 4.60, 95% CI 2.18-9.74), and IgG to native cyclic filaggrin (OR 2.53, 95% CI 1.47-4.34). These six antibodies predicted RA-ILD risk better than all clinical factors combined (optimism-corrected AUC=0·84 versus 0·73). We developed a risk score for RA-ILD combining these antibodies with the clinical factors (smoking, disease activity, glucocorticoid use, obesity). At 50% predicted RA-ILD probability, the risk scores both without (score=2·6) and with (score=5·9) biomarkers achieved specificity ≥93% for RA-ILD. Interpretation: Specific ACPA and anti-native protein antibodies improve RA-ILD prediction. These findings implicate synovial protein antibodies in the pathogenesis of RA-ILD and suggest clinical utility in predicting RA-ILD once validated in external studies. Funding: National Institutes of Health.

Reducing cardiovascular risk with immunomodulators: a randomised active comparator trial among patients with rheumatoid arthritis.

OBJECTIVE: Recent large-scale randomised trials demonstrate that immunomodulators reduce cardiovascular (CV) events among the general population. However, it is uncertain whether these effects apply to rheumatoid arthritis (RA) and if certain treatment strategies in RA reduce CV risk to a greater extent. METHODS: Patients with active RA despite use of methotrexate were randomly assigned to addition of a tumour necrosis factor (TNF) inhibitor (TNFi) or addition of sulfasalazine and hydroxychloroquine (triple therapy) for 24 weeks. Baseline and follow-up 18F-fluorodeoxyglucose-positron emission tomography/CT scans were assessed for change in arterial inflammation, an index of CV risk, measured as an arterial target-to-background ratio (TBR) in the carotid arteries and aorta. RESULTS: 115 patients completed the protocol. The two treatment groups were well balanced with a median age of 58 years, 71% women, 57% seropositive and a baseline disease activity score in 28 joints of 4.8 (IQR 4.0, 5.6). Baseline TBR was similar across the two groups. Significant TBR reductions were observed in both groups-ΔTNFi: -0.24 (SD=0.51), Δtriple therapy: -0.19 (SD=0.51)-without difference between groups (difference in Δs: -0.02, 95% CI -0.19 to 0.15, p=0.79). While disease activity was significantly reduced across both treatment groups, there was no association with change in TBR (ß=0.04, 95% CI -0.03 to 0.10). CONCLUSION: We found that addition of either a TNFi or triple therapy resulted in clinically important improvements in vascular inflammation. However, the addition of a TNFi did not reduce arterial inflammation more than triple therapy. TRIAL REGISTRATION NUMBER: NCT02374021.

A MUC5B Gene Polymorphism, rs35705950-T, Confers Protective Effects Against COVID-19 Hospitalization but Not Severe Disease or Mortality.

Rationale: A common MUC5B gene polymorphism, rs35705950-T, is associated with idiopathic pulmonary fibrosis (IPF), but its role in severe acute respiratory syndrome coronavirus 2 infection and disease severity is unclear. Objectives: To assess whether rs35705950-T confers differential risk for clinical outcomes associated with coronavirus disease (COVID-19) infection among participants in the Million Veteran Program (MVP). Methods: The MUC5B rs35705950-T allele was directly genotyped among MVP participants; clinical events and comorbidities were extracted from the electronic health records. Associations between the incidence or severity of COVID-19 and rs35705950-T were analyzed within each ancestry group in the MVP followed by transancestry meta-analysis. Replication and joint meta-analysis were conducted using summary statistics from the COVID-19 Host Genetics Initiative (HGI). Sensitivity analyses with adjustment for additional covariates (body mass index, Charlson comorbidity index, smoking, asbestosis, rheumatoid arthritis with interstitial lung disease, and IPF) and associations with post-COVID-19 pneumonia were performed in MVP subjects. Measurements and Main Results: The rs35705950-T allele was associated with fewer COVID-19 hospitalizations in transancestry meta-analyses within the MVP (Ncases = 4,325; Ncontrols = 507,640; OR = 0.89 [0.82-0.97]; P = 6.86 × 10-3) and joint meta-analyses with the HGI (Ncases = 13,320; Ncontrols = 1,508,841; OR, 0.90 [0.86-0.95]; P = 8.99 × 10-5). The rs35705950-T allele was not associated with reduced COVID-19 positivity in transancestry meta-analysis within the MVP (Ncases = 19,168/Ncontrols = 492,854; OR, 0.98 [0.95-1.01]; P = 0.06) but was nominally significant (P < 0.05) in the joint meta-analysis with the HGI (Ncases = 44,820; Ncontrols = 1,775,827; OR, 0.97 [0.95-1.00]; P = 0.03). Associations were not observed with severe outcomes or mortality. Among individuals of European ancestry in the MVP, rs35705950-T was associated with fewer post-COVID-19 pneumonia events (OR, 0.82 [0.72-0.93]; P = 0.001). Conclusions: The MUC5B variant rs35705950-T may confer protection in COVID-19 hospitalizations.

Temporal Trends in Clinical Evidence of 5-Year Survival Within Electronic Health Records Among Patients With Early-Stage Colon Cancer Managed With Laparoscopy-Assisted Colectomy vs Open Colectomy.

Importance: Temporal shifts in clinical knowledge and practice need to be adjusted for in treatment outcome assessment in clinical evidence. Objective: To use electronic health record (EHR) data to (1) assess the temporal trends in treatment decisions and patient outcomes and (2) emulate a randomized clinical trial (RCT) using EHR data with proper adjustment for temporal trends. Design, Setting, and Participants: The Clinical Outcomes of Surgical Therapy (COST) Study Group Trial assessing overall survival of patients with stages I to III early-stage colon cancer was chosen as the target trial. The RCT was emulated using EHR data of patients from a single health care system cohort who underwent colectomy for early-stage colon cancer from January 1, 2006, to December 31, 2017, and were followed up to January 1, 2020, from Mass General Brigham. Analyses were conducted from December 2, 2019, to January 24, 2022. Exposures: Laparoscopy-assisted colectomy (LAC) vs open colectomy (OC). Main Outcomes and Measures: The primary outcome was 5-year overall survival. To address confounding in the emulation, pretreatment variables were selected and adjusted. The temporal trends were adjusted by stratification of the calendar year when the colectomies were performed with cotraining across strata. Results: A total of 943 patients met key RCT eligibility criteria in the EHR emulation cohort, including 518 undergoing LAC (median age, 63 [range, 20-95] years; 268 [52%] women; 121 [23%] with stage I, 165 [32%] with stage II, and 232 [45%] with stage III cancer; 32 [6%] with colon adhesion; 278 [54%] with right-sided colon cancer; 18 [3%] with left-sided colon cancer; and 222 [43%] with sigmoid colon cancer) and 425 undergoing OC (median age, 65 [range, 28-99] years; 223 [52%] women; 61 [14%] with stage I, 153 [36%] with stage II, and 211 [50%] with stage III cancer; 39 [9%] with colon adhesion; 202 [47%] with right-sided colon cancer; 39 [9%] with left-sided colon cancer; and 201 [47%] with sigmoid colon cancer). Tests for temporal trends in treatment assignment (χ2 = 60.3; P < .001) and overall survival (χ2 = 137.2; P < .001) were significant. The adjusted EHR emulation reached the same conclusion as the RCT: LAC is not inferior to OC in overall survival rate with risk difference at 5 years of -0.007 (95% CI, -0.070 to 0.057). The results were consistent for stratified analysis within each temporal period. Conclusions and Relevance: These findings suggest that confounding bias from temporal trends should be considered when conducting clinical evidence studies with long time spans. Stratification of calendar time and cotraining of models is one solution. With proper adjustment, clinical evidence may supplement RCTs in the assessment of treatment outcome over time.

Real-world data analyses unveiled the immune-related adverse effects of immune checkpoint inhibitors across cancer types.

Immune checkpoint inhibitors have demonstrated significant survival benefits in treating many types of cancers. However, their immune-related adverse events (irAEs) have not been systematically evaluated across cancer types in large-scale real-world populations. To address this gap, we conducted real-world data analyses using nationwide insurance claims data with 85.97 million enrollees across 8 years. We identified a significantly increased risk of developing irAEs among patients receiving immunotherapy agents in all seven cancer types commonly treated with immune checkpoint inhibitors. By six months after treatment initialization, those receiving immunotherapy were 1.50-4.00 times (95% CI, lower bound from 1.15 to 2.16, upper bound from 1.69 to 20.36) more likely to develop irAEs in the first 6 months of treatment, compared to matched chemotherapy or targeted therapy groups, with a total of 92,858 patients. The risk of developing irAEs among patients using nivolumab is higher compared to those using pembrolizumab. These results confirmed the need for clinicians to assess irAEs among cancer patients undergoing immunotherapy as part of management. Our methods are extensible to characterizing the effectiveness and adverse effects of novel treatments in large populations in an efficient and economical fashion.

Testing the Effects of Disease-Modifying Antirheumatic Drugs on Vascular Inflammation in Rheumatoid Arthritis: Rationale and Design of the TARGET Trial.

Individuals with rheumatoid arthritis (RA) are at increased risk for atherosclerotic cardiovascular disease (ASCVD) events relative to the general population, potentially mediated by atherosclerotic plaques that are more inflamed and rupture prone. We sought to address whether RA immunomodulators reduce vascular inflammation, thereby reducing ASCVD risk, and whether such reduction depends on the type of immunomodulator. The TARGET (Treatments Against RA and Effect on 18-Fluorodeoxyglucose [18 F-FDG] Positron Emission Tomography [PET]/Computed Tomography [CT]) trial (NCT02374021) will enroll 150 patients with RA with active disease and an inadequate response to methotrexate. Participants will be randomized to add either a tumor necrosis factor (TNF) inhibitor (etanercept or adalimumab) or sulfasalazine and hydroxychloroquine to their background methotrexate. Participants will undergo full-body 18 F-FDG-labelled PET scanning at baseline and after 6 months. Efficacy and safety evaluations will occur every 6 weeks, with therapy modified in a treat-to-target approach. The primary outcome is the comparison of change in arterial inflammation in the wall of the aorta and carotid arteries between the randomized treatment groups, specifically, the change in the mean of the maximum target-to-background ratio of arterial 18 F-FDG uptake in the most diseased segment of either the aorta and carotid arteries. A secondary analysis will compare the effects of achieving low disease activity or remission with those of moderate to high disease activity on vascular inflammation. The TARGET trial will test, for the first time, whether RA treatments reduce arterial inflammation and whether such reduction differs according to treatment strategy with either TNF inhibitors or a combination of nonbiologic disease-modifying antirheumatic drugs.

Divergence of Cardiovascular Biomarkers of Lipids and Subclinical Myocardial Injury Among Rheumatoid Arthritis Patients With Increased Inflammation.

OBJECTIVE: Patients with rheumatoid arthritis (RA) are 1.5 times more likely to develop cardiovascular disease (CVD) attributed to chronic inflammation. A decrease in inflammation in patients with RA is associated with increased low-density lipoprotein (LDL) cholesterol. This study was undertaken to prospectively evaluate the changes in lipid levels among RA patients experiencing changes in inflammation and determine the association with concomitant temporal patterns in markers of myocardial injury. METHODS: A total of 196 patients were evaluated in a longitudinal RA cohort, with blood samples and high-sensitivity C-reactive protein (hsCRP) levels measured annually. Patients were stratified based on whether they experienced either a significant increase in inflammation (an increase in hsCRP of ≥10 mg/liter between any 2 time points 1 year apart; designated the increased inflammation cohort [n = 103]) or decrease in inflammation (a decrease in hsCRP of ≥10 mg/liter between any 2 time points 1 year apart; designated the decreased inflammation cohort [n = 93]). Routine and advanced lipids, markers of inflammation (interleukin-6, hsCRP, soluble tumor necrosis factor receptor II), and markers of subclinical myocardial injury (high-sensitivity cardiac troponin T [hs-cTnT], N-terminal pro-brain natriuretic peptide) were measured. RESULTS: Among the patients in the increased inflammation cohort, the mean age was 59 years, 81% were women, and the mean RA disease duration was 17.9 years. The average increase in hsCRP levels was 36 mg/liter, and this increase was associated with significant reductions in LDL cholesterol, triglycerides, total cholesterol, apolipoprotein (Apo B), and Apo A-I levels. In the increased inflammation cohort at baseline, 45.6% of patients (47 of 103) had detectable circulating hs-cTnT, which further increased during inflammation (P = 0.02). In the decreased inflammation cohort, hs-cTnT levels remained stable despite a reduction in inflammation over follow-up. In both cohorts, hs-cTnT levels were associated with the overall estimated risk of CVD. CONCLUSION: Among RA patients who experienced an increase in inflammation, a significant decrease in routinely measured lipids, including LDL cholesterol, and an increase in markers of subclinical myocardial injury were observed. These findings highlight the divergence in biomarkers of CVD risk and suggest a role in future studies examining the benefit of including hs-cTnT for CVD risk stratification in RA.

Patterns of Tumor Necrosis Factor Inhibitor (TNFi) Biosimilar Use Across United States Rheumatology Practices.

OBJECTIVE: It is unclear if biosimilars of biologics for inflammatory arthritis are realizing their promise to increase competition and improve accessibility. This study evaluates biosimilar tumor necrosis factor inhibitor (TNFi) utilization across rheumatology practices in the United States and compares whether patients initiating biosimilars remain on these treatments at least as long as new initiators of bio-originators. METHODS: We identified a cohort of patients initiating a TNFi biosimilar between January 2017 and September 2018 from an electronic health record registry containing data from 218 rheumatology practices and over 1 million rheumatology patients in the United States. We also identified a cohort of patients who initiated the bio-originator TNFi during the same period. We calculated the proportion of biosimilar prescriptions compared with other TNFi's and compared persistence on these therapies, adjusting for age, sex, diagnoses codes, and insurance type. RESULTS: We identified 909 patients prescribed the biosimilar infliximab-dyyb, the only biosimilar prescribed, and 4413 patients with a new prescription for the bio-originator infliximab. Biosimilar patients tended to be older, have a diagnosis code for rheumatoid arthritis, and covered by Medicare insurance. Over the study period, biosimilar prescriptions reached a maximum of 3.5% of all TNFi prescriptions. Patients persisted on the biosimilar at least as long as the bio-originator infliximab (hazard ratio [HR] 0.83, P = 0.07). CONCLUSION: The uptake of biosimilars in the United States remains low despite persistence on infliximab-dyyb being similar to the infliximab bio-originator. These results add to clinical studies that should provide greater confidence to patients and physicians regarding biosimilar use.