Data source
The FIRST registry is a multi-institutional cohort of RA patients treated with b/tsDMARDs, established by the University of Occupational and Environmental Health, Japan and its multiple affiliated hospitals. Detail of the cohort is available in other articles [5-9]. In this registry, all registered RA patients were enrolled in a long-term observational study at the point of new or switch-prescription of b/tsDMARDs. If a patient was treated with several b/tsDMARDs, each episode was treated as an independent episode.
By December 2020, 4,115 cases were enrolled in the registry. B/tsDMARDs with the following four different mechanisms of action (classes) were included:
TNFis: infliximab, etanercept, adalimumab, golimumab, and certolizumab,
Interleukin-6 receptor inhibitors (IL-6Ris): tocilizumab and sarilumab,
Cytotoxic T-lymphocyte–associated antigen-4 immunoglobulin (CTLA4-Ig): abatacept, and Janus kinase inhibitors (JAKis): tofacitinib, baricitinib, peficitinib, and upadacitinib. Biosimilar bDMARDs were also included. Rituximab was not included in this study because this drug is not approved as a treatment option for RA in Japan.
At the start of b/tsDMARD treatment, baseline data were collected, including demographics, disease duration, titres of anti–cyclic citrullinated protein antibody (ACPA), measures of disease activity, functional status, present and past treatments, serum creatinine levels, and coexistence of chronic lung diseases including chronic bronchitis, bronchial ectasia, interstitial pneumonia, old tuberculosis, and inflammatory lung nodule, and the names of coexisting diseases. Follow-up data on disease activity were collected at 2 weeks, 6 months, and 1 year after the start of the therapy. If a treatment was discontinued within a year due to severe adverse events (SAEs), data about the date and the reason of treatment cessation were also collected.
Patient selection and data collection
Eligibility criteria
As outcomes of the treatment may differ when treatment options are limited, this study included only cases that were enrolled in the FIRST registry after JAKis were first approved in Japan—from August 2013 to December 2020.
Exclusion criteria
To remove cases who administrated b/tsDMARDs as a part of treatment of comorbidities (e.g. interstitial lung diseases or vasculitis), cases have both significant comorbidities and PSL > 15 mg/day were excluded from the analysis. Information about the coexisting diseases was collected.
Definition of D2T-RA, very D2T-RA, and SAE
Based on the EULAR definition, cases that failed to achieve the treatment target with ≥ 2 classes of b/tsDMARDs were identified as D2T-RA. In addition to this definition, we categorised cases that failed ≥ 3 classes of b/tsDMARDs as very D2T-RA (vD2T-RA). Cases treated with a b/tsDMARD for the first time were assigned to the b/tsDMARD-naïve group.
Adverse reactions (e.g., allergic reactions, infections, malignancies, lymphoproliferative disorders, major adverse cardiovascular events, and abnormality in laboratory tests) which lead treatment discontinuation was recorded as SAEs.
Statistical analysis
Simple comparison of patient background
Backgrounds of the patients in each treatment subgroups were compared across the 4 classes of b/tsDMARDs using one-way analysis of variance (ANOVA) for numerical variables and chi-square test for categorical variables.
Latent class analysis
To identify different patterns of drug response, latent class analysis was conducted using the gsem suite of functions in Stata 16 (StataCorp, College Station, TX), which categorised patients into classes. We estimated 3 classes according to previous studies [6, 8]. The number and percentage of cases that fell in a particular latent class was calculated for each treatment type.
Panel analysis
Change in CDAI and HAQ-DI over the course of a year was compared using longitudinal panel data analysis. Regression analyses were conducted using the xt suite of functions in Stata16. A mixed-effects regression model was fitted with age, gender, disease duration, CDAI at week 0, positivity of ACPA, coexistence of pulmonary diseases, and serum creatinine as fixed effects and use of MTX and glucocorticoids as random effects.
Analyses using a propensity-based inverse-probability treatment weighted method
As the number of patients included for each group (D2T-RA and vD2T-RA groups) was limited, multivariate regression analysis tended to overfit the data, resulting in biased estimation. Therefore, propensity-based inverse-probability treatment weighted (IPTW) method was also performed for sensitivity analysis. For CDAI and HAQ-DI, delta (D) CDAI and DHAQ-DI were calculated, respectively, as follows and then used for outcomes.
DValue = (Value at 1year) − (Value at day 0)
A regression model was used to adjust for potential confounders. Variables that showed a significant correlation with failures to ≥ 2 classes of b/tsDMARDs were included as covariates. Missing data including loss to follow-up were managed as blank and no imputation was conducted because of homogeneity of the data.
Analysis of hazards of SAE
Development of SAEs is another major factor that causes D2T-RA status. Therefore, Cox regression analysis controlling for age, gender, dose of MTX and glucocorticoid at day 0 was conducted to assess the hazards of SAE over the course of a year among each treatment subgroup. As the risk of SAEs may differ depending on past failures of b/tsDMARDs, the risk among the b/tsDMARD-naïve group was also analysed. Nelson-Aalen cumulative hazard model was also used to show the time trend of hazard accumulation visually.