Disease activity drives transcriptomic heterogeneity in early untreated rheumatoid synovitis

OBJECTIVES: Transcriptomic profiling of synovial tissue from patients with early, untreated rheumatoid arthritis (RA) was used to explore the ability of unbiased, data-driven approaches to define clinically relevant subgroups. METHODS: RNASeq was performed on 74 samples, with disease activity data collected at inclusion. Principal components analysis (PCA) and unsupervised clustering were used to define patient clusters based on expression of the most variable genes, followed by pathway analysis and inference of relative abundance of immune cell subsets. Histological assessment and multiplex immunofluorescence (for CD45, CD68, CD206) were performed on paraffin sections. RESULTS: PCA on expression of the (n=894) most variable genes across this series did not divide samples into distinct groups, instead yielding a continuum correlated with baseline disease activity. Two patient clusters (PtC1, n=52; PtC2, n=22) were defined based on expression of these genes. PtC1, with significantly higher disease activity and probability of response to methotrexate therapy, showed upregulation of immune system genes; PtC2 showed upregulation of lipid metabolism genes, described to characterise tissue resident or M2-like macrophages. In keeping with these data, M2-like:M1-like macrophage ratios were inversely correlated with disease activity scores and were associated with lower synovial immune infiltration and the presence of thinner, M2-like macrophage-rich synovial lining layers. CONCLUSION: In this large series of early, untreated RA, we show that the synovial transcriptome closely mirrors clinical disease activity and correlates with synovial inflammation. Intriguingly, lower inflammation and disease activity are associated with higher ratios of M2:M1 macrophages, particularly striking in the synovial lining layer. This may point to a protective role for tissue resident macrophages in RA.