Objective

Rheumatoid arthritis (RA) fibroblast-like synoviocytes (FLS) derived from hip and knee have distinctive DNA methylation and transcriptome patterns in interleukin (IL)-6 signaling and Janus kinase (JAK)-signal transducers and activators of transcription (STAT) pathways. To determine the functional effects of these joint-specific signatures, we evaluated how RA hip and knee FLS differ in their response to IL-6.

Methods

Hip or knee RA FLS were obtained after arthroplasty. Previously published datasets on epigenetic landscape of FLS were mined to identify joint-specific IL-6-related epigenomic differences. RNA sequencing was performed on five RA hip and five knee FLS treated with or without IL-6. Differential gene expression was determined using edgeR software. STAT3 phosphorylation was measured using bead assays. Sensitivity to tofacitinib was evaluated by measuring CCL2 inhibition using quantitative polymerase chain reaction.

Results

Assay for Transposase-Accessible Chromatin sequencing and histone chromatin immunoprecipitation sequencing datasets from RA FLS were analyzed to identify epigenomic differences between hip and knee. Differential chromatin accessibility was associated with IL-6,IL-6R, and JAK1 genes. H3K27ac was also differentially marked at other JAK-STAT-related genes, including STAT3-STAT5A region. Principal component analysis of RNA sequencing data confirmed segregation between RA hip and knee FLS under basal conditions, that persisted following IL-6 treatment. STAT3 phosphorylation after IL-6 was significantly higher in knee than hip FLS and was highly correlated with JAK1 protein levels. Knee FLS were less sensitive to the JAK inhibitor tofacitinib than hip FLS.

Conclusion

RA hip and knee FLS have distinct transcriptomes, epigenetic marks, and STAT3 activation patterns in the IL-6 pathway. These joint-specific differences might contribute to a differential clinical response in individual joints to targeted therapies such as JAK inhibitors.

Rheumatoid arthritis (RA) is an immune-mediated disease affecting diarthrodial joints that remains an unmet medical need despite improved therapy. This limitation likely reflects the diversity of pathogenic pathways in RA, with individual patients demonstrating variable responses to targeted therapies. Better understanding of RA pathogenesis would be aided by a more complete characterization of the disease. To tackle this challenge, we develop and apply a systems biology approach to identify important transcription factors (TFs) in individual RA fibroblast-like synoviocyte (FLS) cell lines by integrating transcriptomic and epigenomic information. Based on the relative importance of the identified TFs, we stratify the RA FLS cell lines into two subtypes with distinct phenotypes and predicted active pathways. We biologically validate these predictions for the top subtype-specific TF RARα and demonstrate differential regulation of TGFβ signaling in the two subtypes. This study characterizes clusters of RA cell lines with distinctive TF biology by integrating transcriptomic and epigenomic data, which could pave the way towards a greater understanding of disease heterogeneity.