Differences in trajectory of disease activity according to biologic and targeted synthetic disease-modifying anti-rheumatic drug treatment in patients with rheumatoid arthritis
Arthritis Research & Therapy volume 24, Article number: 233 (2022)
The purpose of this study was to stratify patients with rheumatoid arthritis (RA) according to the trend of disease activity by trajectory-based clustering and to identify contributing factors for treatment response to biologic and targeted synthetic disease-modifying anti-rheumatic drugs (DMARDs) according to trajectory groups.
We analyzed the data from a nationwide RA cohort from the Korean College of Rheumatology Biologics and Targeted Therapy registry. Patients treated with second-line biologic and targeted synthetic DMARDs were included. Trajectory modeling for clustering was used to group the disease activity trend. The contributing factors using the machine learning model of SHAP (SHapley Additive exPlanations) values for each trajectory were investigated.
The trends in the disease activity of 688 RA patients were clustered into 4 groups: rapid decrease and stable disease activity (group 1, n = 319), rapid decrease followed by an increase (group 2, n = 36), slow and continued decrease (group 3, n = 290), and no decrease in disease activity (group 4, n = 43). SHAP plots indicated that the most important features of group 2 compared to group 1 were the baseline erythrocyte sedimentation rate (ESR), prednisolone dose, and disease activity score with 28-joint assessment (DAS28) (SHAP value 0.308, 0.157, and 0.103, respectively). The most important features of group 3 compared to group 1 were the baseline ESR, DAS28, and estimated glomerular filtration rate (eGFR) (SHAP value 0.175, 0.164, 0.042, respectively). The most important features of group 4 compared to group 1 were the baseline DAS28, ESR, and blood urea nitrogen (BUN) (SHAP value 0.387, 0.153, 0.144, respectively).
The trajectory-based approach was useful for clustering the treatment response of biologic and targeted synthetic DMARDs in patients with RA. In addition, baseline DAS28, ESR, prednisolone dose, eGFR, and BUN were important contributing factors for 4-year trajectories.
Patients with rheumatoid arthritis (RA) suffer from chronic inflammatory arthritis and a variety of extra-articular symptoms. Because RA patients develop long-term disabilities and decreased quality of life, they need intensive treatment early on, including treat-to-target strategies [1, 2]. Second-line disease-modifying anti-rheumatic drugs (DMARDs), such as biologic DMARDs (bDMARDs) or targeted synthetic DMARDs (tsDMARDs), are used in the treatment of RA patients who do not respond to treatment with conventional synthetic DMARDs (csDMARDs) [3,4,5]. However, because individual patients respond differently to bDMARDs or tsDMARDs, selecting the drug that will induce the best response in each patient remains challenging in clinical practice. Indeed, finding predictors of optimal treatment response to second-line DMARDs is necessary not only to improve the prognoses, but also to reduce patient suffering and medical costs [6, 7].
Although various predictors for treatment response have been proposed thus far [8,9,10,11], most studies have evaluated the treatment response at a specific time point (e.g., 3 or 6 months) without considering the changes in disease activity over time in the predictive models. Among patients treated with bDMARDs or tsDMARDs, some show improvements in disease activity early during treatment, while others respond slowly to treatment and their prescription is either stopped or switched to another drug. In addition, even if the initial response to treatment is good, some patients show aggravation of disease activity over time, and some patients do not improve disease activity even if the biological agent is maintained or switched. Given that changes can vary over a relatively long period of time, it is important to evaluate long-term treatment responses and to consider disease activity at multiple time points.
Trajectory-based clustering is a group-based approach that has been recently used in several studies on patients with RA [12,13,14,15,16,17,18,19,20]. This method can determine the clusters of individuals with a similar course of disease activity and progression over time and identify different patterns of response to treatment. Using predictive markers that are associated with each distinct disease activity pattern can increase the efficiency of the treatment response and maintain treatment effects for a long time.
The purpose of the present study was twofold: (1) to investigate trajectories by classifying long-term changes in treatment response to bDMARDs or tsDMARDs in patients with RA and (2) to find contributing factors in each trajectory using machine learning models.
Study population and data collection
This study used data from the KOBIO registry, a nationwide multicenter cohort in Korea that was established to evaluate the effectiveness and side effects of bDMARDs (abatacept, adalimumab, etanercept, golimumab, infliximab, rituximab, and tocilizumab) and tsDMARDs (tofacitinib and baricitinib) in patients with RA . Patients in the registry were recruited from 38 hospitals since 2012, and their demographics, medications, comorbidities, extra-articular manifestations, disease activities, radiographic findings, and laboratory findings have been recorded by investigators from each participating hospital. A total of 2,122 patients who were treated with second-line DMARDs including bDMARDs and tsDMARDs between December 2012 and June 2019, were registered. Patients who had been treated with bDMARDs or tsDMARDs before registration in the KOBIO registry were excluded from the study. The number of patients and number of follow-ups that each patient completed, were counted to find the most appropriate follow-up period for clustering. Patients who completed four follow-ups were selected. Patients who did not complete four follow-ups due to various reasons, including refusal of further registration, loss of follow-up, or less than four years of follow-up, were excluded from the study. Ethical approval of the KOBIO registry was obtained from the institutional review boards of all 38 participating institutions, including the Institutional Review Board of Inje University Seoul Paik Hospital (PAIK 2018–11-005).
Trajectory-based clustering model
Group-based trajectory modeling (GBTM) is a model for clustering sequential data and has often been used in clinical studies . The purpose of GBTM is to find groups with different characteristics by clustering response patterns over time and to compare the clustered groups. The model assumes that the dependent variable follows a certain probability distribution such as normal distribution or Poisson distribution and that the mean of the probability distribution takes the form of polynomials over time. Under the above assumptions, parameters and functional forms are fitted by the data using maximum likelihood estimation. In our study, the sequential treatment records using biologics were clustered by response pathway based on the model, and the probability distribution and parameters were determined by the medical rationale and data explanatory power. The probability of distribution of the dependent variable was chosen as the censored normal distribution. The number of clusters and the order of polynomials were determined to maximize the value of the model selection criterion, Bayesian Information Criterion, by performing several combinations. After clustering, we showed the individual response patterns to confirm the empirical results. Group-based trajectory modeling was performed using Stata software.
Machine learning model for finding contributing factors of trajectories
A machine learning model was built to find clustered groups and contributing factors. Clustered groups were used as outputs and potential important factors as inputs, and the contribution of each variable was evaluated. First, an XGBoost (eXtreme Gradient Boosting) model was employed to fit a prediction function whose output represented a long-term clinical path . We used an XGBoost as this is a tree-based ensemble model known for its high performance in many machine learning prediction tasks. To determine the degree of the added variable contribution to the prediction in the model, we used a SHAP (SHapley Additive exPlanations) value, an additive feature attribution method used for interpretability when building complex models for high prediction power such as ensemble models and deep learning . By using both XGBoost and the SHAP value, we were able to identify clinical parameters that contributed to the differentiation of poor and good response pathways and thus, had potential clinical implications. XGBoost and SHAP values were estimated using R software version 3.6.1 (R Foundation for Statistical Computing, Vienna, Austria).
For comparison of baseline characteristics among groups, the chi-squared and Kruskal–Wallis tests were used for categorical and continuous variables, respectively. Statistical analyses were performed using R software. All data are shown as mean (standard deviation [SD]) or percentage values.
Patient selection and grouping according to disease trajectory
Among the patients treated with tsDMARDs, only those treated with tofacitinib were included in the study. Of the 2122 patients, we excluded 353 patients treated with bDMARDs or tsDMARDs prior to registration in the KOBIO registry, and 66 patients with missing data. Among 1703, the second, third, fourth, fifth, sixth, and seventh follow-ups were 1307, 959, 688, 454, 217, and 61, respectively. Of the 1703 patients, 688 completed four follow-ups and therefore, were selected for clustering.
The 688 study patients were allocated to 4 distinct groups according to the trajectory of DAS28 (Fig. 1). In group 1 (n = 319), the disease activity rapidly decreased and remained below moderate; in group 2 (n = 36), the disease activity rapidly decreased but then increased; in group 3 (n = 290), the disease activity showed a continuous but slow decrease; and in group 4 (n = 43), the disease activity remained high. The trajectories of each individual patient are shown in Fig. 1B. The baseline characteristics, including the number of current smokers, BUN, cholesterol, and initial DAS28, differed significantly among the groups (Table 1).
Prescription patterns of bDMARDs and tsDMARDs according to trajectory groups
Table 2 shows the prescription patterns of bDMARDs and tsDMARDs according to trajectory groups. The proportion of patients who maintained their first prescribed bDMARDs or tsDMARDs during the follow-up period (never switching) in group 1 was 86.5%, which was numerically higher than those in groups 2, 3, and 4 (52.8%, 50.3%, and 25.6%, respectively). The mean number of bDMARDs or tsDMARDs prescribed was numerically higher in group 4 (2.3) than in groups 1 (1.1), 2 (1.7), and 3 (1.6). The proportion of patients who were prescribed non-tumor necrosis factor (TNF) inhibitor as initial treatment was numerically higher in group 1 (44.2%) than in groups 2, 3, and 4 (30.5%, 23.4%, and 18.6%, respectively).
Contributing factors for trajectory-based clustered groups
Figure 2 shows SHAP plots illustrating the contribution of each of the clinical features to the trajectory-based clustered groups. High SHAP values represent variables that have a high contribution to the trajectory group, either with a positive impact (right side of the zero point on the x-axis), or a negative impact (negative values on the x-axis); the color (yellow to purple) indicates the feature value, which varies from “low” (yellow) to “high” (purple). Figure 2A indicates the contributing features to the trajectory-based cluster group 1 compared to groups 2, 3, and 4. The baseline ESR was the most important feature (SHAP value 0.312) to the group 1 trajectory, indicating that the lower its feature value (yellow dots, right side), the more ESR related to group 1, and the higher its feature value (purple dots, left side), the less ESR related to group 1. The second important feature was the baseline DAS28 (SHAP value 0.293), which followed a similar pattern to baseline ESR: the higher the DAS28 feature value (purple dots), the less this feature related to group 1, and the lower the feature value (light purple and yellow dots), the more DAS28 was related to group 1. The third important feature was BUN (SHAP value 0.155): the higher it was (light purple dots), the more it was related to groups 2, 3, and 4 (right side).
Figure 2B shows the contributing features of group 2 compared to group 1. The first important feature was the baseline ESR (SHAP value 0.308). The SHAP plot indicates that the higher the ESR feature value, the more ESR was related to group 2. The second important feature was the prednisolone dose (SHAP value 0.157): the higher it was, the more it was related to group 2, although the tendency was not high. The third important feature was baseline DAS28 (SHAP value 0.103), and the median value was related to group 1. Figure 2C shows the important features of group 3 compared to group 1. The baseline ESR was the first important feature (SHAP value 0.175). As indicated in the SHAP plot, the higher the ESR feature value, the more it was related to group 3. The second important feature was the baseline DAS28 (SHAP value 0.164), with a similar pattern as baseline ESR’s contribution to group 3. A third important feature was eGFR (SHAP value 0.042): the higher it was, the more it was related to group 3. Figure 2D shows the important features of group 4 compared to group 1. The first important feature was the baseline DAS28 (SHAP value 0.387): the higher its feature value, the more it was related to group 4. ESR (SHAP value 0.167) and BUN (SHAP value 0.144) followed a similar contribution pattern to baseline DAS28: the higher their feature values, the more they related to group 4. Other variables, such as symptom duration and CRP (Fig. 2A), hemoglobin and AST (Fig. 2B), BUN and cholesterol (Fig. 2C), and hemoglobin and ALT (Fig. 2D) were the fourth and fifth important features.
In this study, we clustered the trajectories of the disease activity of RA patients during 4 years of treatment with bDMARDs or tsDMARDs into 4 patterns. From those trajectories, we found several distinguishing baseline features of 3 trajectories compared with the trajectory with good treatment response (i.e., group 1 with a rapid decrease in disease activity without subsequent increases). Among the baseline features, DAS28 and ESR were common important contributing factors associated with treatment response trajectory in patients with RA treated with bDMARDs or tsDMARDs. In addition, laboratory findings, such as BUN, CRP, hemoglobin, AST, eGFR, cholesterol, and ALT and clinical findings such as prednisolone dose and symptom duration were also high-ranked important features for trajectories.
Various trajectory studies have been conducted on RA patients. A study using the Swedish BARFOT cohort identified three trajectories: best (39.6%), moderate (41.5%), and worst (18.9%) outcomes . Trajectories using Australia’s cohort were similarly divided into good (43.8%), moderate (39.7%), and poor (16.5%) outcome groups and there was a significant difference in BMI and the proportion of ever smokers between the poor outcomes group and the good outcomes group. . A study using the CATCH cohort in Canada divided the cohort into 5 trajectories, which were similar to those of the BARFOT cohort . The strength of our study is that the treatment responses for bDMARDs or tsDMARDs were clustered into 4 trajectories from a 4-year longitudinal cohort data, which is a relatively long period compared with those of other studies.
The latent class mixed model by RA-MAP consortium from the TACERA longitudinal cohort study takes categorical latent variables and assumes each latent class is a function of covariates . On the other hand, group-based trajectory modeling in our study is suitable for classifying trajectories of the dependent variable by assuming that the dependent variable has a polynomial trajectory. Under the assumption that treatment response patterns can be expressed as changes in DAS28, we applied group-based trajectory modeling. In addition to estimating latent classes, important contributing factors of latent classes based on machine learning were provided for clinical implications.
In terms of the switching treatments (bDMARDs or tsDMARDs), as much as 86.5% of patients in group 1 maintained the initially prescribed drug without switching. Moreover, they responded more quickly to treatment with bDMARDs or tsDMARDs than those in other groups, and their disease activity remained low for nearly 4 years. Interestingly, patients in group 1 were more commonly prescribed non-TNF inhibitors as an initial treatment than in those in other groups. In groups 3 and 4, 50.3% and 25.6% of patients maintained one bDMARDs or tsDMARDs, respectively, and 10.0% and 41.9% of patients switched among the DMARDs twice or more, respectively. Differences in trajectories, such as delayed or poor treatment response despite frequent drug replacement, suggest the need for personalized drug selection for each RA patient [25,26,27]. Moreover, considering the multifactorial characteristics of RA, new drugs countering the various mechanisms of RA must be developed [28, 29].
SHAP value can be expressed as a plot to recognize the interaction of complex predictors . A method of finding predictors of the treatment response of bDMARDs or tsDMARDs that is affected by many variables can help us understand the trajectory of the therapeutic response. In this study, baseline ESR and DAS28 were the most important features in common for predicting trajectory. Interestingly, higher eGFR was associated with group 3 compared to group 1. That is, the higher the renal function, the slower the treatment response may be. In addition, median DAS28 value related to group 2 cannot be captured in a typical regression analysis.
Previous studies have suggested that hemoglobin is associated with disease activity [30, 31]. Our results indicate that hemoglobin may be a contributing factor related to long-term treatment response in some RA patients. Also, cholesterol was the fifth important feature among variables contributing to group 3 compared to group 1. Patients with high baseline cholesterol are unlikely to show a good treatment response .
Our study had some limitations. First, we did not consider changes in disease activity related to the type and dosage of DMARDs and individual adherence. Therefore, it may be difficult to generalize the results of this study to individual patients. Second, a small number of tsDMARD users were included in this study, which limits the generalizability of our results. Third, the number of patients and/or duration of the follow-up may have limited the number of trajectories. Indeed, we might discover more than four trajectories if increasing the number of patients and duration of follow-up. Fourth, the inclusion of only patients who completed the fourth follow-up in this study may have a selection bias in which patients who responded well to treatment were primarily selected. In this study, patients who due to a variety of reasons did not complete 4 follow-ups, were excluded. Fifth, important features of the baseline for the trajectory were identified. These important features are potential predictive factors, but statistical verification is required in various studies to confirm they are indeed risk factors.
In conclusion, we were able to cluster the trend of disease activity over 4 years in patients treated with bDMARDs or tsDMARDs into 4 distinct trajectories and found that DAS28, ESR, BUN, CRP, hemoglobin, AST, eGFR, cholesterol, ALT, prednisolone dose, and symptom duration may be important contributing factors for trajectories. Our study suggests that the trajectory-based clustering approach for disease activity may be useful in predicting treatment responses from longitudinal data in real-world practice and making decisions about treatment plans in patients with RA.
Availability of data and materials
Data are available from the Clinical Research Committee of KOBIO under the Korean College of Rheumatology for researchers who meet the criteria for access to confidential data. To request data, please contact Kichul Shin, MD, PhD, Director of the Korean College of Rheumatology Biologics Registry, Associate Professor of the Division of Rheumatology, Director of Logistics Planning at SMG-SNU, Boramae Medical Center, 20 Boramae-ro-5-gil, Dongjak-gu, Seoul, 07,061, Korea; Tel: + 82–2-870–3204; Fax: + 82–2-870–3866; Email: rk.ca.uns@1bedik.
Anti-cyclic citrullinated peptide antibody
Biologic disease-modifying anti-rheumatic drugs
Blood urea nitrogen
Disease activity score with 28-joint assessment
Disease activity scores in 28 joints using the erythrocyte sedimentation rate
Disease-modifying anti-rheumatic drugs
Estimated glomerular filtration rate
Erythrocyte sedimentation rate
Group-based trajectory modeling
Korean College of Rheumatology Biologics & Targeted therapy Registry
SHapley Additive exPlanations
Tumor necrosis factor
Targeted synthetic disease-modifying anti-rheumatic drugs
EXtreme Gradient Boosting
Grigor C, Capell H, Stirling A, McMahon AD, Lock P, Vallance R, Kincaid W, Porter D. Effect of a treatment strategy of tight control for rheumatoid arthritis (the TICORA study): a single-blind randomised controlled trial. Lancet. 2004;364(9430):263–9.
Kavanaugh A, van Vollenhoven RF, Fleischmann R, Emery P, Sainsbury I, Florentinus S, Chen S, Guerette B, Kupper H, Smolen JS. Testing treat-to-target outcomes with initial methotrexate monotherapy compared with initial tumour necrosis factor inhibitor (adalimumab) plus methotrexate in early rheumatoid arthritis. Ann Rheum Dis. 2018;77(2):289–92.
Lau CS, Chia F, Dans L, Harrison A, Hsieh TY, Jain R, Jung SM, Kishimoto M, Kumar A, Leong KP, et al. 2018 update of the APLAR recommendations for treatment of rheumatoid arthritis. Int J Rheum Dis. 2019;22(3):357–75.
Singh JA, Saag KG, Bridges SL Jr, Akl EA, Bannuru RR, Sullivan MC, Vaysbrot E, McNaughton C, Osani M, Shmerling RH, et al. 2015 American College of Rheumatology Guideline for the Treatment of Rheumatoid Arthritis. Arthritis Rheumatol. 2016;68(1):1–26.
Smolen JS, Landewe RBM, Bijlsma JWJ, Burmester GR, Dougados M, Kerschbaumer A, McInnes IB, Sepriano A, van Vollenhoven RF, de Wit M, 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.
Grabner M, Boytsov NN, Huang Q, Zhang X, Yan T, Curtis JR. Costs associated with failure to respond to treatment among patients with rheumatoid arthritis initiating TNFi therapy: a retrospective claims analysis. Arthritis Res Ther. 2017;19(1):92.
Navarro F, Martinez-Sesmero JM, Balsa A, Peral C, Montoro M, Valderrama M, Gomez S, de Andres-Nogales F, Casado MA, Oyaguez I. Cost-effectiveness analysis of treatment sequences containing tofacitinib for the treatment of rheumatoid arthritis in Spain. Clin Rheumatol. 2020;39(10):2919–30.
Lequerre T, Rottenberg P, Derambure C, Cosette P, Vittecoq O. Predictors of treatment response in rheumatoid arthritis. Joint Bone Spine. 2019;86(2):151–8.
Cuppen BV, Welsing PM, Sprengers JJ, Bijlsma JW, Marijnissen AC, van Laar JM, Lafeber FP, Nair SC. Personalized biological treatment for rheumatoid arthritis: a systematic review with a focus on clinical applicability. Rheumatology (Oxford). 2016;55(5):826–39.
Katchamart W, Johnson S, Lin HJ, Phumethum V, Salliot C, Bombardier C. Predictors for remission in rheumatoid arthritis patients: a systematic review. Arthritis Care Res (Hoboken). 2010;62(8):1128–43.
Wijbrandts CA, Tak PP. Prediction of response to targeted treatment in rheumatoid arthritis. Mayo Clin Proc. 2017;92(7):1129–43.
Barnabe C, Sun Y, Boire G, Hitchon CA, Haraoui B, Thorne JC, Tin D, van der Heijde D, Curtis JR, Jamal S, et al. Heterogeneous disease trajectories explain variable radiographic, function and quality of life outcomes in the Canadian Early Arthritis Cohort (CATCH). PLoS ONE. 2015;10(8):e0135327.
Leu Agelii M, Andersson M, Jones BL, Sjowall C, Kastbom A, Hafstrom I, Forslind K, Gjertsson I. Disease activity trajectories in rheumatoid arthritis: a tool for prediction of outcome. Scand J Rheumatol. 2021;50(1):1–10.
Norton S, Sacker A, Dixey J, Done J, Williams P, Young A, Early Rheumatoid Arthritis S. Trajectories of functional limitation in early rheumatoid arthritis and their association with mortality. Rheumatology (Oxford). 2013;52(11):2016–24.
Siemons L, Ten Klooster PM, Vonkeman HE, Glas CA, Van de Laar M. Distinct trajectories of disease activity over the first year in early rheumatoid arthritis patients following a treat-to-target strategy. Arthritis Care Res (Hoboken). 2014;66(4):625–30.
Wabe N, Wojciechowski J, Wechalekar MD, Cleland LG, McWilliams L, Lee A, Proudman S, Wiese MD. Disease activity trajectories in early rheumatoid arthritis following intensive DMARD therapy over 3 years: association with persistence to therapy. Int J Rheum Dis. 2017;20(10):1447–56.
Bykerk VP, Lee EB, van Vollenhoven R, Gruben DC, Fallon L, Woolcott JC, Keystone E. Identification of distinct disease activity trajectories in methotrexate-naive patients with rheumatoid arthritis receiving tofacitinib over twenty-four months. Arthritis Care Res (Hoboken). 2022;74(1):131–41.
Consortium R-M. Characterization of disease course and remission in early seropositive rheumatoid arthritis: results from the TACERA longitudinal cohort study. Ther Adv Musculoskelet Dis. 2021;13:1759720X211043977.
Norton S, Fu B, Scott DL, Deighton C, Symmons DP, Wailoo AJ, Tosh J, Lunt M, Davies R, Young A, et al. Health Assessment Questionnaire disability progression in early rheumatoid arthritis: systematic review and analysis of two inception cohorts. Semin Arthritis Rheum. 2014;44(2):131–44.
Pan Y, Norton S, Gwinnutt JM, Kearsley-Fleet L, Symmons DPM, Lunt M, Young A, Consortium B-RCC, Hyrich KL, Verstappen SMM. Not all moderate disease is the same - identification of disability trajectories among patients with rheumatoid arthritis and moderate disease activity. PLoS One. 2019;14(5):e0215999.
Shin K, Choi SJ, Kwak S-K, Park Y-B, Sung Y-K, Choi C-B, Jeon CH, Lee J, Nam EJ, Kim S-K, et al. AB1066 Korean Biologics Registry of Patients with Systemic Rheumatic Disease (KOBIO): a nationwide registry to assess adverse events associated with biologic treatment in Korea. Ann Rheum Dis. 2014;73:1153–4.
Nagin DS, Odgers CL. Group-based trajectory modeling in clinical research. Annu Rev Clin Psychol. 2010;6:109–38.
Chen T, Guestrin C. XGBoost a scalable tree boosting system. In: Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. San Francisco: Association for Computing Machinery; 2016. p. 785–94.
Lundberg SM, Lee S-I. A unified approach to interpreting model predictions: Proceedings of the 31st International Conference on Neural Information Processing Systems December 2017. Adv Neural Inform Process Syst. 2017;30:4768–77.
Isaacs JD, Ferraccioli G. The need for personalised medicine for rheumatoid arthritis. Ann Rheum Dis. 2011;70(1):4.
Tak PP. A personalized medicine approach to biologic treatment of rheumatoid arthritis: a preliminary treatment algorithm. Rheumatology (Oxford). 2012;51(4):600–9.
Koo BS, Eun S, Shin K, Yoon H, Hong C, Kim D-H, Hong S, Kim Y-G, Lee C-K, Yoo B, et al. Machine learning model for identifying important clinical features for predicting remission in patients with rheumatoid arthritis treated with biologics. Arthritis Res Ther. 2021;23(1):178.
Chaudhari K, Rizvi S, Syed BA. Rheumatoid arthritis: current and future trends. Nat Rev Drug Discovery. 2016;15(5):305–6.
Shams S, Martinez JM, Dawson JRD, Flores J, Gabriel M, Garcia G, Guevara A, Murray K, Pacifici N, Vargas MV, et al. The therapeutic landscape of rheumatoid arthritis: current state and future directions. Front Pharmacol. 2021;12(1233):680043.
Moller B, Scherer A, Forger F, Villiger PM, Finckh A, Swiss Clinical Quality Management Program for Rheumatic D. Anaemia may add information to standardised disease activity assessment to predict radiographic damage in rheumatoid arthritis: a prospective cohort study. Ann Rheum Dis. 2014;73(4):691–6.
Wilson A, Yu HT, Goodnough LT, Nissenson AR. Prevalence and outcomes of anemia in rheumatoid arthritis: a systematic review of the literature. Am J Med. 2004;116(Suppl 7A):50S-57S.
Park YJ, Cho CS, Emery P, Kim WU. LDL cholesterolemia as a novel risk factor for radiographic progression of rheumatoid arthritis: a single-center prospective study. PLoS ONE. 2013;8(7):e68975.
The registry was funded by the Korean College of Rheumatology, which had no involvement in the study design, collection, analysis, and interpretation of the data; in writing the manuscript; or in the decision to submit the manuscript for publication. We would like to thank all members of the Korean College of Rheumatology Biology (KOBIO) registry. We thank Dr. Joon Seo Lim from the Scientific Publications Team at Asan Medical Center for his editorial assistance in preparing this manuscript.
This work was supported by the following grants: the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (No. NRF-2021R1C1C1009815); a grant (2021IP0047) from the Asan Institute for Life Sciences, Asan Medical Center, Seoul, Korea; The Technology Innovation Program (20004927, Upgrade of CDM based Distributed Biohealth Data Platform and Development of Verification Technology) funded by the Ministry of Trade, Industry & Energy (MOTIE, Korea); and the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (HR21C0198).
Ethics approval and consent to participate
Ethical approval of the KOBIO registry was obtained from the institutional review boards of all 38 participating institutions, including the Institutional Review Board of Inje University Seoul Paik Hospital (PAIK 2018–11-005).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Koo, B.S., Eun, S., Shin, K. et al. Differences in trajectory of disease activity according to biologic and targeted synthetic disease-modifying anti-rheumatic drug treatment in patients with rheumatoid arthritis. Arthritis Res Ther 24, 233 (2022). https://doi.org/10.1186/s13075-022-02918-3
- Rheumatoid arthritis
- Trajectory clustering/trajectory modeling
- Treatment response