Objectives: Patients with rheumatoid arthritis (RA) display different trajectories towards improvement of disease. We aimed to disentangle the heterogeneity of RA disease trajectories from the first clinical visit onwards using graph-based pseudotime analysis. Methods: We studied early patients with RA over 1.5 years in 2 data sets: Leiden (Netherlands), n = 1237, with 5017 visits, and Towards a Cure for Early Rheumatoid Arthritis (TACERA) (United Kingdom), n = 243, with 750 visits. We created a pipeline for time-independent clustering of clinical and haematologic features to identify disease states. Sequence analyses of these states defined the trajectories. We studied the predictability of the trajectories with baseline features. Results: Clustering identified 8 disease states with localised inflammation (joints) and systemic inflammation (erythrocyte sedimentation rate [ESR] or leucocytes) as the main discriminating factors. The disease state sequences consisted of 4 trajectories, which we independently replicated in TACERA: A, high ESR; B, rapid progression from many inflamed joints towards remission; C, high leucocytes; and D, many inflamed joints with poor prognosis. Systemic vs local inflammation patterns showed moderate predictability at baseline (sensitivity of 71% and precision of 0.73 for trajectory A, although lower precision of 0.52 for trajectory B), while other trajectories were less predictable. Trajectories C and D had strong resemblance with B at baseline but deteriorated into less favourable trajectories. Patients in trajectory A were more often female and on average older. The trajectories were not explained by time till disease-modifying antirheumatic drug, baseline disease activity, or symptom duration. The suboptimal trajectories coincided with worse patient-reported outcomes, even when the inflammation was mainly systemic. Conclusions: We identified 4 distinct trajectories in early RA, differentiating RA into localised vs systemic inflammation. Our results highlight potential differences in disease pathology and opportunities for further targeted treatment. Inevitably, patterns without linkage to our selected features could not be detected.

Objective: To facilitate patient disease subset and risk factor identification by constructing a pipeline which is generalizable, provides easily interpretable results, and allows replication by overcoming electronic health records (EHRs) batch effects. Material and Methods: We used 1872 billing codes in EHRs of 102 880 patients from 12 healthcare systems. Using tools borrowed from single-cell omics, we mitigated center-specific batch effects and performed clustering to identify patients with highly similar medical history patterns across the various centers. Our visualization method (PheSpec) depicts the phenotypic profile of clusters, applies a novel filtering of noninformative codes (Ranked Scope Pervasion), and indicates the most distinguishing features. Results: We observed 114 clinically meaningful profiles, for example, linking prostate hyperplasia with cancer and diabetes with cardiovascular problems and grouping pediatric developmental disorders. Our framework identified disease subsets, exemplified by 6 "other headache"clusters, where phenotypic profiles suggested different underlying mechanisms: migraine, convulsion, injury, eye problems, joint pain, and pituitary gland disorders. Phenotypic patterns replicated well, with high correlations of ≥0.75 to an average of 6 (2-8) of the 12 different cohorts, demonstrating the consistency with which our method discovers disease history profiles. Discussion: Costly clinical research ventures should be based on solid hypotheses. We repurpose methods from single-cell omics to build these hypotheses from observational EHR data, distilling useful information from complex data. Conclusion: We establish a generalizable pipeline for the identification and replication of clinically meaningful (sub)phenotypes from widely available high-dimensional billing codes. This approach overcomes datatype problems and produces comprehensive visualizations of validation-ready phenotypes.