ObjectivesTo integrate published single-cell RNA sequencing (scRNA-seq) data and assess the contribution of synovial fibroblast (SF) subsets to synovial pathotypes and respective clinical characteristics in treatment-naïve early arthritis.MethodsIn this in silico study, we integrated scRNA-seq data from published studies with additional unpublished in-house data. Standard Seurat, Harmony and Liger workflow was performed for integration and differential gene expression analysis. We estimated single cell type proportions in bulk RNA-seq data (deconvolution) from synovial tissue from 87 treatment-naïve early arthritis patients in the Pathobiology of Early Arthritis Cohort using MuSiC. SF proportions across synovial pathotypes (fibroid, lymphoid and myeloid) and relationship of disease activity measurements across different synovial pathotypes were assessed.ResultsWe identified four SF clusters with respective marker genes: PRG4+ SF (CD55, MMP3, PRG4, THY1neg); CXCL12+ SF (CXCL12, CCL2, ADAMTS1, THY1low); POSTN+ SF (POSTN, collagen genes, THY1); CXCL14+ SF (CXCL14, C3, CD34, ASPN, THY1) that correspond to lining (PRG4+ SF) and sublining (CXCL12+ SF, POSTN+ + and CXCL14+ SF) SF subsets. CXCL12+ SF and POSTN+ + were most prominent in the fibroid while PRG4+ SF appeared highest in the myeloid pathotype. Corresponding, lining assessed by histology (assessed by Krenn-Score) was thicker in the myeloid, but also in the lymphoid pathotype + the fibroid pathotype. PRG4+ SF correlated positively with disease severity parameters in the fibroid, POSTN+ SF in the lymphoid pathotype whereas CXCL14+ SF showed negative association with disease severity in all pathotypes.ConclusionThis study shows a so far unexplored association between distinct synovial pathologies and SF subtypes defined by scRNA-seq. The knowledge of the diverse interplay of SF with immune cells will advance opportunities for tailored targeted treatments.

Background Genome-wide association studies have reported more than 100 risk loci for rheumatoid arthritis (RA). These loci are shown to be enriched in immune cell-specific enhancers, but the analysis so far has excluded stromal cells, such as synovial fibroblasts (FLS), despite their crucial involvement in the pathogenesis of RA. Here we integrate DNA architecture, 3D chromatin interactions, DNA accessibility, and gene expression in FLS, B cells, and T cells with genetic fine mapping of RA loci. Results We identify putative causal variants, enhancers, genes, and cell types for 30–60% of RA loci and demonstrate that FLS account for up to 24% of RA heritability. TNF stimulation of FLS alters the organization of topologically associating domains, chromatin state, and the expression of putative causal genes such as TNFAIP3 and IFNAR1. Several putative causal genes constitute RA-relevant functional networks in FLS with roles in cellular proliferation and activation. Finally, we demonstrate that risk variants can have joint-specific effects on target gene expression in RA FLS, which may contribute to the development of the characteristic pattern of joint involvement in RA. Conclusion Overall, our research provides the first direct evidence for a causal role of FLS in the genetic susceptibility for RA accounting for up to a quarter of RA heritability.

Glucocorticoids (GCs) exert potent anti-inflammatory effects in immune cells through the glucocorticoid receptor (GR). Dendritic cells (DCs), central actors for coordinating immune responses, acquire tolerogenic properties in response to GCs. Tolerogenic DCs (tolDCs) have emerged as a potential treatment for various inflammatory diseases. To date, the underlying cell type-specific regulatory mechanisms orchestrating GC-mediated acquisition of immunosuppressive properties remain poorly understood. In this study, we investigated the transcriptomic and epigenomic remodeling associated with differentiation to DCs in the presence of GCs. Our analysis demonstrates a major role of MAFB in this process, in synergy with GR. GR and MAFB both interact with methylcytosine dioxygenase TET2 and bind to genomic loci that undergo specific demethylation in tolDCs. We also show that the role of MAFB is more extensive, binding to thousands of genomic loci in tolDCs. Finally, MAFB knockdown erases the tolerogenic properties of tolDCs and reverts the specific DNA demethylation and gene upregulation. The preeminent role of MAFB is also demonstrated in vivo for myeloid cells from synovium in rheumatoid arthritis following GC treatment. Our results imply that, once directly activated by GR, MAFB plays a critical role in orchestrating the epigenomic and transcriptomic remodeling that define the tolerogenic phenotype.