Utilizing biologic disease-modifying anti-rheumatic treatment sequences to subphenotype rheumatoid arthritis

Arthritis Research & Therapy

Research

Tianxi Cai¹  Priyam Das²  Katherine P. Liao³  Nancy A. Shadick⁴  Dana Weisenfeld⁴  Kumar Dahal⁴  Vivi Feathers⁴  Michael E. Weinblatt⁴  Jonathan S. Coblyn⁴  Debsurya De⁵ 

[1] Department of Biomedical Informatics, Harvard Medical School, Boston, MA, USA;Department of Biomedical Informatics, Harvard Medical School, Boston, MA, USA;Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA;Department of Biomedical Informatics, Harvard Medical School, Boston, MA, USA;Division of Rheumatology, Inflammation, and Immunity, Brigham and Women’s Hospital, 60 Fenwood Road, 02115, Boston, MA, USA;Division of Rheumatology, Inflammation, and Immunity, Brigham and Women’s Hospital, 60 Fenwood Road, 02115, Boston, MA, USA;Indian Statistical Institute, Kolkata, India;

关键词: Rheumatoid arthritis;    Medication prescriptions;    Biologic disease-modifying anti-rheumatic drugs;    Electronic health record;    Mixture model;    Markov chain;   

DOI  :  10.1186/s13075-023-03072-0

 received in 2022-11-23, accepted in 2023-05-20,  发布年份 2023

来源: Springer

【 摘 要 】

BackgroundMany patients with rheumatoid arthritis (RA) require a trial of multiple biologic disease-modifying anti-rheumatic drugs (bDMARDs) to control their disease. With the availability of several bDMARD options, the history of bDMARDs may provide an alternative approach to understanding subphenotypes of RA. The objective of this study was to determine whether there exist distinct clusters of RA patients based on bDMARD prescription history to subphenotype RA.MethodsWe studied patients from a validated electronic health record-based RA cohort with data from January 1, 2008, through July 31, 2019; all subjects prescribed ≥ 1 bDMARD or targeted synthetic (ts) DMARD were included. To determine whether subjects had similar b/tsDMARD sequences, the sequences were considered as a Markov chain over the state-space of 5 classes of b/tsDMARDs. The maximum likelihood estimator (MLE)-based approach was used to estimate the Markov chain parameters to determine the clusters. The EHR data of study subjects were further linked with a registry containing prospectively collected data for RA disease activity, i.e., clinical disease activity index (CDAI). As a proof of concept, we tested whether the clusters derived from b/tsDMARD sequences correlated with clinical measures, specifically differing trajectories of CDAI.ResultsWe studied 2172 RA subjects, mean age 52 years, RA duration 3.4 years, and 62% seropositive. We observed 550 unique b/tsDMARD sequences and identified 4 main clusters: (1) TNFi persisters (65.7%), (2) TNFi and abatacept therapy (8.0%), (3) on rituximab or multiple b/tsDMARDs (12.7%), (4) prescribed multiple therapies with tocilizumab predominant (13.6%). Compared to the other groups, TNFi persisters had the most favorable trajectory of CDAI over time.ConclusionWe observed that RA subjects can be clustered based on the sequence of b/tsDMARD prescriptions over time and that the clusters were correlated with differing trajectories of disease activity over time. This study highlights an alternative approach to consider subphenotyping of patients with RA for studies aimed at understanding treatment response.

【 授权许可】

CC BY   
© The Author(s) 2023

【 图 表 】

40517_2023_252_Article_IEq80.gif

40517_2023_252_Article_IEq69.gif