Article Text
Abstract
Background Prior studies of SLE clusters based on autoantibodies have utilized cross-sectional data from single centers. We applied clustering techniques to longitudinal and comprehensive autoantibody data from a large multinational, multi-ethnic inception cohort of well characterized SLE patients to identify clusters associated with disease outcomes.
Methods We used demographic, clinical, and serological data at enrolment and follow-up visits years 3 and 5 from 805 patients who fulfilled the 1997 Updated ACR SLE criteria and were enrolled within 15 months of diagnosis. For each visit, ANA, dsDNA, Sm, U1-RNP, SSA/Ro60, SSB/La, Ro52/TRIM21, histones, ribosomal P, Jo-1, centromere B, PCNA, anti-DFS70, lupus anticoagulant (LAC), IgG and IgM for anticardiolipin, anti–β2GP1, and aPS/PT, and IgG anti-β2GP1 D1 were performed at a single lab (except LAC). K-means clustering algorithm on principal component analysis (10 dimensions) transformed longitudinal ANA/autoantibody profiles was used. We compared cluster demographic/clinical outcomes, including longitudinal disease activity (total and adjusted mean SLEDAI-2K), SLICC/ACR damage index and organ-specific domains, SLE therapies, and survival, using one-way ANOVA test and a Benjamini-Hochberg correction with false discovery rate alpha=0.05. Results were visualized using t-distributed stochastic neighbor embedding.
Results Four unique patient clusters were identified (table 1). Cluster 1, characterized by high frequency of anti-Sm and anti-RNP over time, was the youngest group at disease onset with a high proportion of subjects of Asian and African ancestry. At year 5, they had the highest disease activity, were more likely to have active hematologic and mucocutaneous involvement, and to be on/exposed to immunosuppressants/biologics. Cluster 2, the largest cluster, had low frequency of anti-dsDNA, were oldest at disease onset, and at year 5, had the lowest disease activity, and were least likely to have nephritis and be on/exposed to immunosuppressants/biologics. Cluster 3 had the highest frequency of antiphospholipid antibodies over time, were more likely to be of European ancestry, have an elevated BMI, be former smokers, and by year 5, to have nephritis, neuropsychiatric involvement, including strokes and seizures (SLICC/ACR damage index). Cluster 4 was characterized by anti-SSA/Ro60, SSB/La, Ro52/TRIM21, histone antibodies, and low complements at year 5. Overall, survival of the 805 subjects was 94% at 5 years, and none of the clusters predicted survival.
Conclusions Four SLE patient clusters associated with disease activity, organ involvement, and treatment were identified in this analysis of longitudinal ANA/autoantibody profiles in relation to SLE outcomes, suggesting these subsets might be identifiable based on extended autoantibody profiles early in disease and carry prognostic information.
Acknowledgments This study is presented on behalf of SLICC. We would also like to acknowledge the technical assistance of Ms. Haiyan Hou (MitogenDx Laboratory).
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.