Abstract
Rheumatoid arthritis (RA) is a systemic autoimmune disease characterized by synovitis, joint damage and deformity. A newly described type of cell death, ferroptosis, has an important role in the pathogenesis of RA. However, the heterogeneity of ferroptosis and its association with the immune microenvironment in RA remain unknown. Synovial tissue samples from 154 RA patients and 32 healthy controls (HCs) were obtained from the Gene Expression Omnibus database. Twelve of twenty-six ferroptosis-related genes (FRGs) were differentially expressed between RA patients and HCs. Furthermore, the patterns of correlation among the FRGs were significantly different between the RA and HC groups. RA patients were classified into two distinct ferroptosis-related clusters, of which cluster 1 had a higher abundance of activated immune cells and a corresponding lower ferroptosis score. Enrichment analysis suggested that tumor necrosis factor-α signaling via nuclear factor-κB was upregulated in cluster 1. RA patients in cluster 1 responded better to anti-tumor necrosis factor (anti-TNF) therapy, which was verified by the GSE 198520 dataset. A diagnostic model to identify RA subtypes and immunity was constructed and verified, in which the area under the curve values in the training (70%) and validation (30%) cohorts were 0.849 and 0.810, respectively. This study demonstrated that there were two ferroptosis clusters in RA synovium that exhibited distinct immune profiles and ferroptosis sensitivity. Additionally, a gene scoring system was constructed to classify individual RA patients.