Skip to main content

Management and treatment outcomes of rheumatoid arthritis in the era of biologic and targeted synthetic therapies: evaluation of 10-year data from the KURAMA cohort

Abstract

Background

Advances in rheumatoid arthritis (RA) treatment, highlighted by biological disease-modifying antirheumatic drugs (bDMARDs) and targeted synthetic DMARDs (tsDMARDs), have altered the paradigm of RA treatment in the last decade. Therefore, real-world clinical evidence is needed to understand how treatment strategies and outcomes have changed.

Methods

Using an observational cohort of RA from 2012 to 2021, we collected cross-sectional data of RA patients annually to analyze a trend in RA management. For patients who initiated b/tsDMRDs, we evaluated treatment outcomes between b/tsDMARDs. Mixed-effect models were applied to examine the statistical implications of changes over time in treatment outcomes with a background adjustment.

Results

We analyzed annual cross-sectional data from 5070 patients and longitudinal data from 1816 patients in whom b/tsDMARDs were initiated between 2012 and 2021. b/tsDMARD use increased, whereas glucocorticoid use decreased from 2012 to 2021. Disease activity and functional disability measures improved over time. The percentage of tsDMARD prescriptions considerably increased. All b/tsDMARDs showed clinical improvements in disease activity and functional disability. Statistically, TNFi showed better short-term improvements in b/tsDMARD-naïve patients, while IL6Ri demonstrated significant long-term benefits. IL6Ri had better retention rates in switched patients. After adjustment for patient characteristics, the annual change of RA disease activity and functional disability fared significantly better from 2012 to 2021.

Conclusions

With the development of new RA therapeutics, overall treatment outcomes advanced in the past decade.

Background

Treatment options for rheumatoid arthritis (RA) have progressed over the past two decades. Methotrexate (MTX), biologic disease-modifying antirheumatic drugs (bDMARDs), and targeted synthetic antirheumatic drugs (tsDMARDs) such as Janus kinase inhibitors (JAKis) have significantly improved disease activity, functional disability, and joint prognosis of RA patients [1, 2]. This clinical evidence is primarily based on well-designed randomized controlled trials; however, they cannot assess long-term efficacy and safety outside of the study period or advocate evidence in the groups of patients who met exclusion criteria due to age or comorbidities [3, 4]. Real-world data obtained from observational cohort studies could not only answer these questions but also describe how patient backgrounds and therapeutic methods affected treatment outcomes in the observational period [5].

The Kyoto University Rheumatoid Arthritis Management Alliance (KURAMA) cohort is an open-label, single-center, observational cohort study of RA and rheumatic diseases [6,7,8]. It was established in 2011, and a total of 4418 patients were registered until 2022. We recorded longitudinal treatment data at every clinical visit, including disease activity, functional disability, and adverse events. We also have conducted the “annual RA survey” that features radiographic examinations on joint destruction, testing for osteoporosis, and patient surveys regarding frailty, sarcopenia, mental status, and more every year since 2012 [9,10,11,12]. The KURAMA cohort also has a biobank of plasma, synovium, and other RA-related specimens [13,14,15].

Between 2012 and 2021, the Japanese Ministry of Health, Labour and Welfare approved two bDMARDs (certolizumab pegol, 2013; sarilumab, 2018) and five JAKis (tofacitinib, 2013; baricitinib, 2017; upadacitinib; peficitinib; and filgotinib, 2021) [16]. This study aimed to investigate how treatment strategies and outcomes evolved alongside advances in therapeutic strategies by reviewing the 10-year experience of the KURAMA cohort from 2012 to 2021.

Materials and methods

Patients

To observe the annual trends in RA management, we gathered clinical data from RA patients, including disease activity, functional disability, the usage of MTX, glucocorticoids (GCs), and b/tsDMARDs, along with demographic and anthropometric details, from the annual RA surveys conducted from 2012 to 2021. For patients treated with b/tsDMARDs, we recorded disease activity; functional disability before and 1, 2, 3, 6, and 12 months after initiation; and the date of initiation and discontinuation. Disease activity was assessed using the Disease Activity Score (DAS) 28-CRP and Clinical Disease Activity Index (CDAI), and functional disability was evaluated by the Health Assessment Questionnaire Disability Index (HAQ) [17, 18]. All patients met the 1987 American College of Rheumatology (ACR) classification criteria or the 2010 ACR/European League Against Rheumatism classification criteria for RA diagnosis. The KURAMA cohort study was approved by the Kyoto University Graduate School of Medicine Medical Ethics Committee (R0357), and informed consent was obtained from all patients.

Data analysis

Since the annual cross-sectional dataset had relatively low missing value rates, with the highest rate at 8.4% for DAS28-CRP, data with missing values were excluded from subsequent statistical analyses.

We employed propensity score matching to address potential confounding factors that could impact the effectiveness of b/tsDMARDs. The reference group consisted of patients treated with CTLA4-Ig (abatacept). The following variables were used to account for potential confounding factors: age, sex, CDAI, DAS28-CRP, HAQ, disease stage, disease class, methotrexate (MTX) and glucocorticoid (GC) dosages, rheumatoid factor (RF), and anti-cyclic citrullinated peptide antibody (ACPA) titers at baseline. One-to-one matching without replacement was performed using the nearest neighbor match on the logit of the propensity score with the caliper width set to the standard deviation of the logit of the propensity score. Student’s t-test and the chi-square test were used to analyze continuous variables and categorical variables, respectively.

For the time-series analysis, we fitted mixed-effects models to evaluate the trends in disease activity and functional disability measures [19]. All models were fitted with patient-specific random intercepts to account for inter-patient variations. Box-Cox transformation and square root transformation were applied to CDAI and HAQ, respectively, to fix skewness and gain a better fitting [20, 21]. Age, sex, BMI, RF, ACPA, Steinbrocker’s stage, class, dose of MTX, and dose of GCs were used as covariates.

Statistical analyses were performed using scipy v1.7.1, statsmodels v0.12.2, and scikit-learn v0.24.2. Propensity score matching was performed using psmpy v0.3.13. Kaplan‒Meier survival estimation was performed using lifelines v0.27.3 running on Python v3.9.7. We used ggalluvial v0.12.5 for the Sankey diagram, car v3.1–2 for Box-Cox transformation, and glmmTMB v1.1.7 for a mixed-effect model, running on R v4.2.0. Graphs were drawn using matplotlib v3.4.3, seaborn v0.11.2, and ggplot2 v3.4.2.

Results

Changes in RA patient background and RA treatment trends from 2012 to 2021

From 2012 to 2021, a total of 1156 patients and a cumulative total of 5070 patients participated in annual RA surveys. The mean age of RA patients increased from 62.9 in 2012 to 65.9 in 2021, and the percentage of patients younger than 60 years old was less than 20% in 2021 (Fig. 1A, B, and Supplementary Table S1). This trend may reflect the aging society of Japan. Regarding treatment, the percentage of patients treated with MTX was over 60% in all periods and slightly decreased between 2012 and 2021 from 70.8 to 64.3% (Fig. 1C). GC usage decreased (40.5–18.6%), whereas b/tsDMARD usage (29.5–53.2%) increased, indicating that more intensive treatment has been administered to RA patients in the past decade (Fig. 1C). The patients’ demographics and treatment methods during the observation period are summarized in Supplementary Table S1.

Fig. 1
figure 1

Change in patients’ demographics and treatment from 2012 to 2021. A Box plots showing the change in the mean patient age. B Change in the proportion of age groups. C Percentage of patients treated with MTX, GCs, and b/tsDMARDs

The transition of disease activity and functional disability in 10 years

We next investigated how disease activity and functional disability changed over 10 years. The average DAS28-CRP sharply decreased from 2012 to 2015, during which the number of patients treated with b/tsDMARDs increased (Fig. 1C) and then remained stable thereafter (Fig. 2A and Supplementary Table S2). The percentage of patients who achieved DAS28-CRP remission gradually increased and reached 79.7% in 2021 (Fig. 2B). Similarly, CDAI decreased, and the rate of CDAI remission also increased from 25.1 to 48.1% in 10 years (Fig. 2C, D, and Supplementary Table S2). Moreover, the median HAQ scores decreased over time from 0.69 to 0.25 (Fig. 2E and Supplementary Table S2).

Fig. 2
figure 2

Evolution of disease activity and functional disability over 10 years. A Change in DAS28-CRP. B Proportion of disease activity categories defined by DAS28-CRP. C Change in CDAI. D Proportion of disease activity categories defined by CDAI. E Change in HAQ

Treatment outcomes of b/tsDMARDs in the KURAMA cohort

Because an increasing number of patients in the KURAMA cohort were treated with b/tsDMARDs and achieved good clinical outcomes from 2012 to 2021, in which several new b/tsDMARDs were approved, we next analyzed the real-world treatment outcomes of b/tsDMARDs using longitudinal observational data on the treatment responses to b/tsDMARDs. We found 1816 new prescriptions (820 for b/tsDMARD-naïve patients and 996 for b/tsDMARD-switching patients) between 2012 and 2021. Baseline demographic data are shown in Supplementary Table S3; the patients in this dataset partially overlap with those in Supplementary Tables S1 and S2. Supplementary Fig. S1 shows the trends in the proportion of initiated b/tsDMARDs categorized by mode of action (MOA) from 2012 to 2021. The proportion of IL-6 receptor inhibitors (IL6Ri; tocilizumab and sarilumab) was almost consistent over 10 years, accounting for approximately 10% of naïve cases and 20% of switch cases. The initiation of TNF inhibitors (TNFi; infliximab, etanercept, adalimumab, golimumab, and certolizumab pegol, including their biosimilars) and CTLA4-Ig decreased in both types of cases from 2012 to 2021. The prescription of JAKis increased from 2014 in switch cases first, and JAKis were more recently initiated in over 20% of naïve patients (Supplementary Fig. S1). Concomitant usage of MTX in b/tsDMARD-naïve patients and switched patients was 74.5% and 63.3% in TNFi, 40.7% and 56.6% in IL6Ri, 47.3% and 48.5% in CTLA4-Ig, and 62.0% and 47.0% in JAKi, respectively (Supplementary Table S3).

We next compared the effectiveness of b/tsDMARDs, assessing the transition of DAS28-CRP, CDAI, and HAQ as depicted in Supplementary Fig. S2. Each MOA achieved low disease activity or remission according to the averaged DAS28-CRP and CDAI at 6 months (Supplementary Fig. S2). Similarly, in switched cases, all MOAs achieved CDAI levels below the threshold for low disease activity at 12 months (Supplementary Fig. S2). To address potential confounding factors, we applied propensity score matching, with the CTLA-4 group as the reference group, to evaluate the drug effectiveness in comparison to CTLA-4. Notably, we excluded JAKis from b/tsDMARDs-naïve cases, as the number of patients treated with JAKis as a first-line therapy was fewer than other MOAs (Supplementary Table S3). In b/tsDMARDs-naïve patients, TNFi showed a statistically better improvement in CDAI and DAS28-CRP at 1 and 2 months, and IL6Ri significantly ameliorated CDAI at 2, 6, and 12 months, and DAS28-CRP at 1, 2, 6, and 12 months compared to CTLA4-Ig (Table 1). All MOAs exhibited gradual improvements in HAQ, with no statistical differences (Table 1). In switched cases, all measures of disease activity and functional disability were comparable across all comparisons (Table 1).

Table 1 P-values computed to compare the effectiveness of b/tsDMARDs in the propensity score-matched groups

The retention of b/tsDMARDs is another indicator of successful treatment [22,23,24,25]. For Kaplan‒Meier analysis, we defined an event as when the prescription of b/tsDAMRDs ceased due to inefficacy and adverse events, and discontinuation due to remission, socioeconomic issues, or the patient’s decision was defined as a censored case. There was no significant difference in the retention rate in b/tsDMARD-naïve patients (Fig. 3A). Moreover, IL6Ri had a significantly better retention rate (p < 0.005) in switched patients (Fig. 3B), which agrees with a previous report [23]. The MOA and continuation/discontinuation status of the first b/tsDMARDs in naïve patients and the choice of MOA when switching are summarized in the Sankey diagram (Supplementary Fig. S3). Different MOAs tended to be chosen when switching from IL6Ri; however, the same MOA was selected in approximately half of the cases when switching from JAKi and TNFi (Supplementary Fig. S3). These results show treatment outcomes of b/tsDMARDs between different MOAs in a real-world setting.

Fig. 3
figure 3

Drug retention rate of b/tsDMARDs. A, B Kaplan–Meier survival curves showing drug retention comparing different MOAs. The X-axis shows the number of days after b/tsDMARDs initiation, and the Y-axis shows the retention probability. Statistics used: log-rank test. A Drug retention in b/tsDMARD-naïve patients. B Drug retention in switched patients

Statistical implications of the annual trends in disease activity and functional disability

To know the statistical implications of observed trends in DAS28-CRP, CDAI, and HAQ from 2012 to 2021 (Fig. 2A–D), we used a mixed-effect model for adjusting patient-specific characteristics [19]. The annual changes in CDAI and HAQ were statistically significant, but not DAS28-CRP (Fig. 4A–C and Table 2). Patients’ demographic data such as age, stage, and class significantly affected HAQ, and b/tsDMARDs, MTX, and GC use were significantly associated with DAS28-CRP and CDAI. b/tsDMARDs use was also associated with HAQ (Table 2).

Fig. 4
figure 4

Mixed-effect model fitted to the annual change in disease activity and functional disability. AC The dots denote the observed median values, and the lines represent the regression line fitted using mixed-effect models. A DAS28-CRP. B CDAI. C HAQ

Table 2 Results of the mixed-effect model

Discussions

RA treatment strategies have evolved over the past 20 years. Two bDMARDs and five tsDMARDs were launched in the last decade, and several drugs with novel MOAs are in clinical trials. Therefore, tracking the RA management of individual patients and analyzing real-world data is crucial for addressing whether the management of RA has improved in the real-world setting. In this study, we analyzed the trends and outcomes of RA treatments in the KURAMA cohort between 2012 and 2021; we showed that (1) the percentage of patients who were treated with GCs and b/tsDMARDs changed significantly, (2) treatment responses and drug retentions differed between different MOAs, and (3) metrics of disease activity and functional disability improved over time. Here, we demonstrated the disparity in b/tsDMARD effectiveness and yearly improvements in the overall disease activity and functional disability using propensity score matching and mixed-effects modeling to account for patient-specific factors.

We found that GC use consistently decreased while b/tsDMARD use increased over a decade, which is consistent with other studies that observed real-world trends of RA medications [5, 26]. However, some reported that the proportion of patients treated with b/tsDMARDs or GCs was stable over the observation period [16, 27]. Regarding treatment outcomes, studies have shown a stable decrease in disease activity and functional disability indices from the 1990s or early 2000s to the early 2010s [5, 27,28,29,30,31,32]. In contrast, some studies showed that the average DAS28 or DAS28 remission rate remained the same after the early 2010s, which is consistent with our data [5, 28]. We found that CDAI decreased, and the CDAI remission rate continued to increase from 2012 to 2021. This discordance could be attributed to the difference between DAS28-CRP and CDAI. CDAI criteria are reportedly more stringent than DAS28-CRP criteria, and DAS28-CRP does not correlate with CDAI when patients achieve DAS28 remission [33,34,35,36]. Aletaha et al. reported that the sensitivity of CDAI to subtle changes in disease activity (ACR < 20) is higher than DAS28 and Simplified Disease Activity Index (SDAI) [37]. In the KURAMA cohort, approximately 80% of patients achieved DAS28-CRP remission in 2021; thus, DAS28-CRP might not have correlated with CDAI in the current study.

Various treatment outcomes of b/tsDMARDs have been observed in randomized control trials and real-world studies [38,39,40,41,42,43]. In this study, we analyzed the treatment outcomes in b/tsDAMRD-naïve patients and switched patients and found the following: (1) in naïve patients, TNFi most effectively improved disease activity and functional disability after 3 months, though there were no significant differences in disease activity, functional disability, or drug retention rate between MOAs at 12 months, and (2) in switched patients, IL6Ri yielded significantly lower DAS28-CRP at 12 months and higher retention rate compared to reagents with other MOAs. These findings agree with a recent systematic review where non-TNFi drugs showed better retention rates, and there were no significant differences between bDMARDs in bDMARD-naïve patients [44]. We also observed lower concomitant usage of MTX in IL6Ri users. This trend may reflect the recommendations for RA treatment [45, 46]. The lower concomitant use of MTX and higher concomitant use of GCs in the CTLA4-Ig group may reflect the differences in patient characteristics such as age, comorbidities, and MTX intolerance [47].

Regarding drug persistence, IL6Ri had a better drug retention rate in the KURAMA cohort. Our findings are consistent with previous studies reported by Li et al., who assessed the treatment outcomes of over 8000 patients and observed the superior drug persistence of IL6Ri. Ebina et al. and Jinno et al. also found that IL6Ris had a better retention rate in the ANSWER cohort, a multicenter observational RA cohort where the KURMA cohort participates [23, 48, 49]. Finally, a systematic review by de Castro et al. showed greater persistence of non-TNFi over TNFi [44]. The improvement of disease activity and functional disability in the KURAMA cohort could be attributed to the evolution of RA treatment strategies during the past decade, in which the usage of b/tsDMARDs increased.

The emergence of b/tsDMARDs has altered RA treatment strategies. In this study, we aimed to assess whether the increased treatment options have indeed benefited RA patients in a real-world setting. To the best of our knowledge, this is the first study that revealed annual improvements in disease activity and functional disability using a mixed-effect model. Although it cannot be certain that the yearly trends in CDAI and HAQ in the mixed-effect models represent the increase in treatment options, some studies support that this may be the case. In the study period, RA refractory to multiple b/tsDMARDs was defined as difficult-to-treat RA [50,51,52,53]. Ochi et al. reported that JAKi is preferable for difficult-to-treat RA [54]. The emergence of JAKis and their usage optimization might positively affect RA management. Another factor that could account for the annual decline in disease activity and functional disability is the encouragement of exercise or physical therapy. Exercise is reported to improve functional disability in RA patients [55]. We previously reported that sarcopenia was associated with worse RA management [56]. Advances in treatment and care may have contributed to the annual improvements in disease activity and functional disability.

This study has the following limitations: First, individual patient circumstances, such as comorbidity, allergy history, and socioeconomic status, can affect the choice and outcome of RA treatment, but these factors were not considered here because this study analyzed a relatively large dataset of 5070 RA patients who participated in the annual RA survey and 1816 RA patients who initiated b/tsDMARDs. Given the nature of this real-world observational cohort study, the presence of other, possibly unknown, confounding factors may have influenced the results. Second, the annual RA survey does not include all RA patients who are regularly treated in our hospital.

Conclusions

In summary, we reviewed the real-world transition of therapeutic strategies and their outcomes in the 10-year history of the KURMA cohort. The disease activity and functional disability metrics of patients with RA improved over a decade with increased use of b/tsDMARDs.

Availability of data and materials

The datasets used during the current study are available from the corresponding author upon reasonable request.

Abbreviations

RA:

Rheumatoid arthritis

MTX:

Methotrexate

bDMARD:

Biologic disease-modifying antirheumatic drug

tsDMARD:

Targeted synthetic antirheumatic drug

JAKi:

Janus kinase inhibitor

KURAMA cohort:

Kyoto University Rheumatoid Arthritis Management Alliance cohort

GC:

Glucocorticoid

DAS:

Disease Activity Score

CDAI:

Clinical Disease Activity Index

HAQ:

Health Assessment Questionnaire Disability Index

ACR:

American College of Rheumatology

ANOVA:

Analysis of variance

MOA:

Mode of action

IL6Ri:

IL-6 receptor inhibitors

TNFi:

TNF inhibitors

References

  1. Aletaha D, Smolen JS. Diagnosis and management of rheumatoid arthritis: a review. JAMA. 2018;320(13):1360–72.

    Article  PubMed  Google Scholar 

  2. Smolen JS, Aletaha D, Barton A, Burmester GR, Emery P, Firestein GS, et al. Rheumatoid arthritis. Nat Rev Dis Primers. 2018;4:18001.

    Article  PubMed  Google Scholar 

  3. Kim HS, Lee S, Kim JH. Real-world evidence versus randomized controlled trial: clinical research based on electronic medical records. J Korean Med Sci. 2018;33(34):e213.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Black N. Why we need observational studies to evaluate the effectiveness of health care. BMJ. 1996;312(7040):1215–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Yamanaka H, Tanaka E, Nakajima A, Furuya T, Ikari K, Taniguchi A, et al. A large observational cohort study of rheumatoid arthritis, IORRA: providing context for today’s treatment options. Mod Rheumatol. 2020;30(1):1–6.

    Article  PubMed  Google Scholar 

  6. Ito Y, Hashimoto M, Hirota K, Ohkura N, Morikawa H, Nishikawa H, et al. Detection of T cell responses to a ubiquitous cellular protein in autoimmune disease. Science. 2014;346(6207):363–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  7. Tsuji H, Ikari K, Ohmura K, Yano K, Furu M, Hashimoto M, et al. Significant joint-destructive association of HLA-DRB1*04:05 independent of DAS28 in rheumatoid arthritis. Ann Rheum Dis. 2019;78(2):284–6.

    Article  PubMed  Google Scholar 

  8. Ishikawa Y, Ikari K, Hashimoto M, Ohmura K, Tanaka M, Ito H, et al. Shared epitope defines distinct associations of cigarette smoking with levels of anticitrullinated protein antibody and rheumatoid factor. Ann Rheum Dis. 2019;78(11):1480–7.

    Article  PubMed  CAS  Google Scholar 

  9. Itaya T, Torii M, Hashimoto M, Tanigawa K, Urai Y, Kinoshita A, et al. Prevalence of anxiety and depression in patients with rheumatoid arthritis before and during the COVID-19 pandemic. Rheumatology (Oxford). 2021;60(4):2023–4.

    Article  PubMed  CAS  Google Scholar 

  10. Minamino H, Katsushima M, Torii M, Yamamoto W, Fujita Y, Ikeda K, et al. Serum vitamin D status inversely associates with a prevalence of severe sarcopenia among female patients with rheumatoid arthritis. Sci Rep. 2021;11(1):20485.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Itaya T, Torii M, Hashimoto M, Jindai K, Yamamoto W, Tanigawa K, et al. Perceptions and behaviours related to COVID-19 in patients with rheumatoid arthritis: a cross-sectional study. Clin Rheumatol. 2022;41(1):45–51.

    Article  PubMed  Google Scholar 

  12. Onishi A, Torii M, Hidaka Y, Uozumi R, Oshima Y, Tanaka H, et al. Efficacy of personalized exercise program on physical function in elderly patients with rheumatoid arthritis at high risk for sarcopenia: study protocol for a randomized controlled trial. BMC Musculoskelet Disord. 2023;24(1):280.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Ishikawa M, Ito H, Furu M, Hashimoto M, Fujii T, Okahata A, et al. Plasma sLOX-1 is a potent biomarker of clinical remission and disease activity in patients with seropositive RA. Mod Rheumatol. 2016;26(5):696–701.

    Article  PubMed  CAS  Google Scholar 

  14. Iwasaki T, Watanabe R, Ito H, Fujii T, Okuma K, Oku T, et al. Dynamics of type I and type II interferon signature determines responsiveness to anti-TNF therapy in rheumatoid arthritis. Front Immunol. 2022;13:901437.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Hashimoto M, Yamazaki T, Hamaguchi M, Morimoto T, Yamori M, Asai K, et al. Periodontitis and Porphyromonas gingivalis in preclinical stage of arthritis patients. PLoS One. 2015;10(4):e0122121.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Fukuda W, Kadoya M, Omoto A, Yanagida T, Isoda Y, Sunaga A, et al. Treatment of rheumatoid arthritis and its outcomes in an aging society: a single-center cohort study in Japan from 2011 to 2020. Arthritis Res Ther. 2022;24(1):190.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. England BR, Tiong BK, Bergman MJ, Curtis JR, Kazi S, Mikuls TR, et al. 2019 update of the American College of Rheumatology recommended rheumatoid arthritis disease activity measures. Arthritis Care Res (Hoboken). 2019;71(12):1540–55.

    Article  PubMed  Google Scholar 

  18. Fries JF, Spitz P, Kraines RG, Holman HR. Measurement of patient outcome in arthritis. Arthritis Rheum. 1980;23(2):137–45.

    Article  PubMed  CAS  Google Scholar 

  19. Eldridge N, Wang Y, Metersky M, Eckenrode S, Mathew J, Sonnenfeld N, et al. Trends in adverse event rates in hospitalized patients, 2010–2019. JAMA. 2022;328(2):173–83.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Legendre P, Borcard D. Box-Cox-chord transformations for community composition data prior to beta diversity analysis. Ecography. 2018;41(11):1820–4.

    Article  Google Scholar 

  21. Yu H, Sang P, Huan T. Adaptive Box-Cox transformation: a highly flexible feature-specific data transformation to improve metabolomic data normality for better statistical analysis. Anal Chem. 2022;94(23):8267–76.

    Article  PubMed  CAS  Google Scholar 

  22. Ebina K, Hashimoto M, Yamamoto W, Ohnishi A, Kabata D, Hirano T, et al. Drug retention and discontinuation reasons between seven biologics in patients with rheumatoid arthritis - the ANSWER cohort study. PLoS One. 2018;13(3):e0194130.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Ebina K, Hirano T, Maeda Y, Yamamoto W, Hashimoto M, Murata K, et al. Drug retention of 7 biologics and tofacitinib in biologics-naive and biologics-switched patients with rheumatoid arthritis: the ANSWER cohort study. Arthritis Res Ther. 2020;22(1):142.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Neovius M, Arkema EV, Olsson H, Eriksson JK, Kristensen LE, Simard JF, et al. Drug survival on TNF inhibitors in patients with rheumatoid arthritis comparison of adalimumab, etanercept and infliximab. Ann Rheum Dis. 2015;74(2):354–60.

    Article  PubMed  CAS  Google Scholar 

  25. Simard JF, Arkema EV, Sundstrom A, Geborek P, Saxne T, Baecklund E, et al. Ten years with biologics: to whom do data on effectiveness and safety apply? Rheumatology (Oxford). 2011;50(1):204–13.

    Article  PubMed  Google Scholar 

  26. Crowson LP, Davis JM 3rd, Hanson AC, Myasoedova E, Kronzer VL, Makol A, et al. Time trends in glucocorticoid use in rheumatoid arthritis during the biologics era: 1999–2018. Semin Arthritis Rheum. 2023;61:152219.

    Article  PubMed  CAS  Google Scholar 

  27. Haugeberg G, Hansen IJ, Soldal DM, Sokka T. Ten years of change in clinical disease status and treatment in rheumatoid arthritis: results based on standardized monitoring of patients in an ordinary outpatient clinic in southern Norway. Arthritis Res Ther. 2015;17(1):219.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Gullick NJ, Ibrahim F, Scott IC, Vincent A, Cope AP, Garrood T, et al. Real world long-term impact of intensive treatment on disease activity, disability and health-related quality of life in rheumatoid arthritis. BMC Rheumatol. 2019;3:6.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Kievit W, Fransen J, de Waal Malefijt MC, den Broeder AA, van Riel PL. Treatment changes and improved outcomes in RA: an overview of a large inception cohort from 1989 to 2009. Rheumatology (Oxford). 2013;52(8):1500–8.

    Article  PubMed  Google Scholar 

  30. Mian AN, Ibrahim F, Scott IC, Bahadur S, Filkova M, Pollard L, et al. Changing clinical patterns in rheumatoid arthritis management over two decades: sequential observational studies. BMC Musculoskelet Disord. 2016;17:44.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Aga AB, Lie E, Uhlig T, Olsen IC, Wierod A, Kalstad S, et al. Time trends in disease activity, response and remission rates in rheumatoid arthritis during the past decade: results from the NOR-DMARD study 2000–2010. Ann Rheum Dis. 2015;74(2):381–8.

    Article  PubMed  CAS  Google Scholar 

  32. Minichiello E, Semerano L, Boissier MC. Time trends in the incidence, prevalence, and severity of rheumatoid arthritis: a systematic literature review. Joint Bone Spine. 2016;83(6):625–30.

    Article  PubMed  Google Scholar 

  33. Martins FM, da Silva JA, Santos MJ, Vieira-Sousa E, Duarte C, Santos H, et al. DAS28, CDAI and SDAI cut-offs do not translate the same information: results from the Rheumatic Diseases Portuguese Register Reuma.pt. Rheumatology (Oxford). 2015;54(2):286–91.

    Article  PubMed  Google Scholar 

  34. Aletaha D, Wang X, Zhong S, Florentinus S, Monastiriakos K, Smolen JS. Differences in disease activity measures in patients with rheumatoid arthritis who achieved DAS, SDAI, or CDAI remission but not Boolean remission. Semin Arthritis Rheum. 2020;50(2):276–84.

    Article  PubMed  Google Scholar 

  35. Rintelen B, Sautner J, Haindl PM, Andel I, Maktari A, Leeb BF. Comparison of three rheumatoid arthritis disease activity scores in clinical routine. Scand J Rheumatol. 2009;38(5):336–41.

    Article  PubMed  CAS  Google Scholar 

  36. Gaujoux-Viala C, Mouterde G, Baillet A, Claudepierre P, Fautrel B, Le Loet X, et al. Evaluating disease activity in rheumatoid arthritis: which composite index is best? A systematic literature analysis of studies comparing the psychometric properties of the DAS, DAS28, SDAI and CDAI. Joint Bone Spine. 2012;79(2):149–55.

    Article  PubMed  Google Scholar 

  37. Aletaha D, Nell VP, Stamm T, Uffmann M, Pflugbeil S, Machold K, et al. Acute phase reactants add little to composite disease activity indices for rheumatoid arthritis: validation of a clinical activity score. Arthritis Res Ther. 2005;7(4):R796-806.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Hetland ML, Christensen IJ, Tarp U, Dreyer L, Hansen A, Hansen IT, et al. Direct comparison of treatment responses, remission rates, and drug adherence in patients with rheumatoid arthritis treated with adalimumab, etanercept, or infliximab: results from eight years of surveillance of clinical practice in the nationwide Danish DANBIO registry. Arthritis Rheum. 2010;62(1):22–32.

    Article  PubMed  CAS  Google Scholar 

  39. Gartlehner G, Hansen RA, Jonas BL, Thieda P, Lohr KN. The comparative efficacy and safety of biologics for the treatment of rheumatoid arthritis: a systematic review and metaanalysis. J Rheumatol. 2006;33(12):2398–408.

    PubMed  CAS  Google Scholar 

  40. Gabay C, Emery P, van Vollenhoven R, Dikranian A, Alten R, Pavelka K, et al. Tocilizumab monotherapy versus adalimumab monotherapy for treatment of rheumatoid arthritis (ADACTA): a randomised, double-blind, controlled phase 4 trial. Lancet. 2013;381(9877):1541–50.

    Article  PubMed  CAS  Google Scholar 

  41. Burmester GR, Lin Y, Patel R, van Adelsberg J, Mangan EK, Graham NM, et al. Efficacy and safety of sarilumab monotherapy versus adalimumab monotherapy for the treatment of patients with active rheumatoid arthritis (MONARCH): a randomised, double-blind, parallel-group phase III trial. Ann Rheum Dis. 2017;76(5):840–7.

    Article  PubMed  CAS  Google Scholar 

  42. Weinblatt ME, Schiff M, Valente R, van der Heijde D, Citera G, Zhao C, et al. Head-to-head comparison of subcutaneous abatacept versus adalimumab for rheumatoid arthritis: findings of a phase IIIb, multinational, prospective, randomized study. Arthritis Rheum. 2013;65(1):28–38.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Schiff M, Keiserman M, Codding C, Songcharoen S, Berman A, Nayiager S, et al. Efficacy and safety of abatacept or infliximab vs placebo in ATTEST: a phase III, multi-centre, randomised, double-blind, placebo-controlled study in patients with rheumatoid arthritis and an inadequate response to methotrexate. Ann Rheum Dis. 2008;67(8):1096–103.

    Article  PubMed  CAS  Google Scholar 

  44. de Castro CT, de Queiroz MJ, Albuquerque FC, Brandao CC, Gerlack LF, Pereira DCR, et al. Real-world effectiveness of biological therapy in patients with rheumatoid arthritis: systematic review and meta-analysis. Front Pharmacol. 2022;13:927179.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Smolen JS, Landewe RBM, Bijlsma JWJ, Burmester GR, Dougados M, Kerschbaumer A, et al. EULAR recommendations for the management of rheumatoid arthritis with synthetic and biological disease-modifying antirheumatic drugs: 2019 update. Ann Rheum Dis. 2020;79(6):685–99.

    Article  PubMed  CAS  Google Scholar 

  46. Kawahito Y, Morinobu A, Kaneko Y, Kohno M, Hirata S, Kishimoto M, et al. Drug treatment algorithm and recommendations from the 2020 update of the Japan College of Rheumatology clinical practice guidelines for the management of rheumatoid arthritis-secondary publication. Mod Rheumatol. 2023;33(1):21–35.

    Article  PubMed  Google Scholar 

  47. Murata K, Ito H, Hashimoto M, Nishitani K, Murakami K, Tanaka M, et al. Elderly onset of early rheumatoid arthritis is a risk factor for bone erosions, refractory to treatment: KURAMA cohort. Int J Rheum Dis. 2019;22(6):1084–93.

    Article  PubMed  CAS  Google Scholar 

  48. Jinno S, Onishi A, Dubreuil M, Hashimoto M, Yamamoto W, Murata K, 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(1):116.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  49. Li KJ, Chang CL, Hsin CY, Tang CH. Switching and discontinuation pattern of biologic disease-modifying antirheumatic drugs and tofacitinib for patients with rheumatoid arthritis in Taiwan. Front Pharmacol. 2021;12:628548.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. de Hair MJH, Jacobs JWG, Schoneveld JLM, van Laar JM. Difficult-to-treat rheumatoid arthritis: an area of unmet clinical need. Rheumatology (Oxford). 2018;57(7):1135–44.

    PubMed  Google Scholar 

  51. Nagy G, Roodenrijs NMT, Welsing PM, Kedves M, Hamar A, van der Goes MC, et al. EULAR definition of difficult-to-treat rheumatoid arthritis. Ann Rheum Dis. 2021;80(1):31–5.

    Article  PubMed  Google Scholar 

  52. Nagy G, Roodenrijs NMT, Welsing PMJ, Kedves M, Hamar A, van der Goes MC, et al. EULAR points to consider for the management of difficult-to-treat rheumatoid arthritis. Ann Rheum Dis. 2022;81(1):20–33.

    Article  PubMed  CAS  Google Scholar 

  53. Watanabe R, Hashimoto M, Murata K, Murakami K, Tanaka M, Ohmura K, et al. Prevalence and predictive factors of difficult-to-treat rheumatoid arthritis: the KURAMA cohort. Immunol Med. 2022;45(1):35–44.

    Article  PubMed  CAS  Google Scholar 

  54. Ochi S, Sonomoto K, Nakayamada S, Tanaka Y. Preferable outcome of Janus kinase inhibitors for a group of difficult-to-treat rheumatoid arthritis patients: from the FIRST Registry. Arthritis Res Ther. 2022;24(1):61.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  55. Baillet A, Payraud E, Niderprim VA, Nissen MJ, Allenet B, Francois P, et al. A dynamic exercise programme to improve patients’ disability in rheumatoid arthritis: a prospective randomized controlled trial. Rheumatology (Oxford). 2009;48(4):410–5.

    Article  PubMed  Google Scholar 

  56. Torii M, Hashimoto M, Hanai A, Fujii T, Furu M, Ito H, et al. Prevalence and factors associated with sarcopenia in patients with rheumatoid arthritis. Mod Rheumatol. 2019;29(4):589–95.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors express our sincere gratitude to all people involved with the KURAMA cohort project. We thank Takashi Nakamura, the emeritus of the Department of Orthopeadic Surgery, Kyoto University, for supporting the establishment of the KURAMA cohort. We also thank Soichiro Masuda, Toshiki Fujisawa, and Koji Kawakami (Department of Pharmacoepidemiology, Graduate School of Medicine and Public Health, Kyoto University) for the advice and assistance with statistical analyses.

Funding

The Department of Advanced Medicine for Rheumatic Diseases, Kyoto University Graduate School of Medicine, is supported by Nagahama City, Shiga, Japan; Toyooka City, Hyogo, Japan; Asahi Kasei Pharma Corp.; and AYUMI Pharmaceutical Co. The KURAMA cohort study is supported by a grant from Daiichi Sankyo Co., Ltd. We would like to extend our heartfelt gratitude to Haruo Horii, Naohiro Ito, and Masatoshi Fujii for their generous financial support for the KURAMA cohort. The funders were not involved in the study design, data collection, analysis, interpretation of the data, writing of the manuscript, or in any decision to publish the results.

Author information

Authors and Affiliations

Authors

Contributions

Takayuki F. wrote the manuscript. Takayuki F, A.M., Koi. M., H.O., A.O., K.N., M.T., Kos. M, M.F., R.W., M.H., H.I, Takao F., and T.M collected and organized data. S.M. supervised the study.

Corresponding author

Correspondence to Takayuki Fujii.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Kyoto University Graduate School of Medicine Medical Ethics Committee (R0357), and informed consent was obtained from all patients.

Consent for publication

Not applicable.

Competing interests

The Department of Advanced Medicine for Rheumatic Diseases, Kyoto University Graduate School of Medicine, is supported by Nagahama City, Shiga, Japan; Toyooka City, Hyogo, Japan; Asahi Kasei Pharma Corp.; and AYUMI Pharmaceutical Co. The KURAMA cohort study is supported by a grant from Daiichi Sankyo Co. Ltd. Takayuki.F. speaker fees from AbbVie GK, Asahi Kasei Pharma Corp., Chugai Pharmaceutical Co., Ltd., Eisai Co., Ltd., and Janssen Pharmaceutical K.K. Kos.M. received speaker’s fees and/or consulting fees from Eisai Co. Ltd, Chugai Pharmaceutical Co. Ltd., Pfizer Inc., Bristol-Myers Squibb, Mitsubishi Tanabe Pharma Corporation, UCB Japan Co. Ltd, Daiichi Sankyo Co. Ltd., and Astellas Pharma Inc. H.O. has received speaker fees from AbbVie, Asahi Kasei, Astellas Pharma Inc., Eisai Co. Ltd., Janssen Pharmaceutical K.K., Mitsubishi Tanabe Pharma Corporation, and Daiichi Sankyo Co. Ltd. A.O. has received research grants and/or speaker fees from Pfizer Inc., Bristol-Myers Squibb., Advantest, Asahi Kasei Pharma Corp., Chugai Pharmaceutical Co. Ltd., Eli Lilly Japan K.K, Ono Pharmaceutical Co., UCB Japan Co., Mitsubishi Tanabe Pharma Co., Eisai Co. Ltd., Abbvie Inc., Takeda Pharmaceutical Co. Ltd., and Daiichi Sankyo Co. Ltd. K.N. received a speaking fee and/or consulting fees from Asahi Kasei Pharma Corp., Ayumi Pharmaceutical Co, Daiichi Sankyo Co. Ltd, Hisamitsu Pharmaceutical Co.,Inc, and Kyocera. M.H. received research grants and/or speaker fees from Abbvie, Asahi Kasei, Astellas, Ayumi, Brystol Meyers, Chugai, EA Pharma, Eisai, Daiichi Sankyo, Eli Lilly, Novartis Pharma, and Tanabe Mitsubishi. R.W. has received a research grant from AbbVie and speaker’s fee from Asahi Kasei, Chugai, Eli Lilly, and GSK. M.T. received research grants and/or speaker fees from AbbVie GK, Asahi Kasei Pharma Corp., Astellas Pharma Inc., Chugai Pharmaceutical Co., Ltd., Daiichi Sankyo Co., Ltd., Eisai Co., Ltd., Eli Lilly Japan K.K., Janssen Pharmaceutical K.K., Kyowa Kirin Co., Ltd., Pfizer Inc., Taisho Pharmaceutical Co., Ltd., Tanabe Mitsubishi Pharma Corp., Teijin Pharma, Ltd., UCB Japan Co., Ltd. Koi.M. received a speaking fee and/or consulting fees from Eisai Co. Ltd., Chugai Pharmaceutical Co. Ltd.; Asahi Kasei Pharma Corp., Bristol-Myers Squibb, Mitsubishi Tanabe Pharma Co., Janssen Pharmaceutical K.K. and Daiichi Sankyo Co. Ltd. S. M. received speaker fees from Pfizer Inc. A.M. has received honorarium from AbbVie G.K., Chugai Pharmaceutical Co. Ltd., Eli Lilly Japan K.K., Eisai Co. Ltd., Pfizer Inc., Bristol-Myers Squibb., Mitsubishi Tanabe Pharma Co., Astellas Pharma Inc., and Gilead Sciences Japan., and has received research grants from AbbVie G.K., Asahi Kasei Pharma Corp., Chugai Pharmaceutical Co. Ltd., Mitsubishi Tanabe Pharma Co. and Eisai Co. Ltd. Takao.F. has received lecture fees or honoraria from AbbVie GK; Asahi Kasei Pharma Corporation; Astellas Pharma Inc.; Chugai Pharmaceutical Co. Ltd.; Eisai Co. Ltd.; Eli Lilly Japan K.K.; Janssen Pharmaceutical K.K.; Mitsubishi Tanabe Pharma Corporation; Ono Pharmaceutical Co. Ltd.; and Pfizer Japan Inc.; has received grants from AbbVie GK; Asahi Kasei Pharma Corporation; Chugai Pharmaceutical Co. Ltd.; Eisai Co. Ltd.; and Mitsubishi Tanabe Pharma Corporation. T.M. received speaker’s fee from Asahi Kasei Pharma, AbbVie, Eisai, Jansen Mitsubishi Tanabe, Taisho, and Pfizer. H.I. has received a research grant and/or speaker fee from Bristol-Myers Squibb. The other authors declare that they have no conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: S. Table1.

Baseline characteristics of patients with RA who enrolled in the annual RA survey.

Additional file 2: S. Table 2.

Changes in disease activity and functional disability of the KURAMA cohort over the last decade.

Additional file 3: S. Table 3.

Baseline characteristics of b/tsDMARDs-initiated patients.

Additional file 4: Fig. S1.

Trend in b/tsDMARDs prescription. Bar charts showing the percentages of prescribed b/tsDMARDs in each year. Fig. S2. Treatment outcomes of b/tsDMARDs. A and B. Averaged treatment response comparing each MOA from baseline to 12 months. Left panel: DAS28-CRP, middle panel: CDAI, right panel: HAQ. Upper dotted line in DAS28-CRP: 2.7, lower dotted line in DAS28-CRP: 2.3, upper dotted line in CDAI: 10, lower dotted line in CDAI: 2.8, dotted line in HAQ: 0.5. A. b/tsDMARD-naïve patients. B. Switched patients. Fig. S3. Selection of the initial b/tsDMARDs for naïve patients and MOA selection pattern at the time of switching. Sankey diagram showing retention or switching of b/tsDMARDs initiated in KURAMA cohort for naïve patients, and the choice of MOA at the time of switching. The diagram shows the first through the fourth b/tsDMARDs. MOA: mode of action.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fujii, T., Murata, K., Onizawa, H. et al. Management and treatment outcomes of rheumatoid arthritis in the era of biologic and targeted synthetic therapies: evaluation of 10-year data from the KURAMA cohort. Arthritis Res Ther 26, 16 (2024). https://doi.org/10.1186/s13075-023-03251-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13075-023-03251-z

Keywords