To date, knowledge regarding the epigenomic abnormalities related to transcriptional regulation has been accumulated from studies on cultured RASFs. Notably, most conventional studies involve RASFs purified by subculture; hence, the term “culture” has been emphasized.
DNA methylation
Changes in DNA methylation is the most widely studied epigenetic modification in RASFs. In principle, the expression of genes having CpG islands in the promoter region is suppressed upon methylation by DNA methyltransferase (DNMT). In a study comparing DNA methylation of unstimulated RASFs and SFs from osteoarthritis (OA) patients using methylation arrays, genes in the pathways associated with the pathogenesis of RA (e.g., cell adhesion, cell migration, and interaction between cells and ECM) were found to be specifically hypomethylated in RASFs [5]. It has been indicated that during RA development, epigenomic modifications occur selectively in SFs [6, 7]. Interestingly, differences in DNA methylation of SFs reportedly depend on the RA stage [8, 9]. The recent-onset RA has been found to harbor methylation abnormalities in the gene Shroom Family Member 1 (SHROOM1) that encodes a protein involved in microtubule reconstitution during neural tissue development and cell division. Such genes as SHROOM1 may therefore constitute biomarkers for the early diagnosis of RA [9]. Furthermore, advanced RA is more intractable to treatment and easily relapses after discontinuation of the drug compared with early RA. Disease relapse could occur from epigenetic modification of SFs caused by the chronic inflammatory environment in the joints. Inflammatory cytokines, represented by tumor necrosis factor (TNF)-α and IL-1β, have been reported to suppress DNMT expression and are involved in DNA demethylation.
Histone protein acetylation and methylation
Extensive details on abnormalities due to posttranslational modification of histone proteins in RASFs have emerged in recent years. Briefly, when histone proteins are acetylated, DNA binding is loosened to allow the binding of transcription factors, thereby facilitating transcription of the gene. Acetylation is regulated by a balance between histone acetyltransferase (HAT) and histone deacetylase (HDAC). Further, histone protein methylation is induced by histone methyltransferase (HMT). The state of transcriptional activation changes depends on the combination of these modifications and the specific histones or amino acids subjected to the modifications.
Reports state that unstimulated RASFs have transcriptionally active histone modifications (H3K4me3) in the genomic regions encoding MMPs that are directly involved in the degradation of the joints. MMP expression is suppressed by attenuating H3K4 methylation through the knockdown of WD Repeat Domain 5 (WDR5), a component of the HMT complex [10]. Another study reported that the promoter region of IL6 is subjected to activating histone modification (H3ac) in unstimulated RASFs, and HAT inhibitors suppress the expression of IL6 [11]. In addition, the promoter region of the T-Box Transcription Factor 5 (TBX5), which regulates the expression of chemokines (e.g., C-X-C Motif Chemokine Ligand [CXCL] 8, CXCL12, and C-C Motif Chemokine Ligand 20), is hypomethylated in RASFs compared with that in OASFs, accompanied by activating histone modifications (H3K4me3 and H3ac) [12]. Loh et al. reported that arthritis-related genes that are persistently expressed in TNF-α stimulate RASFs to show open chromatin structures with activating histone modifications (H3K27ac), and the binding motifs of specific transcription factors (nuclear factor-κB, interferon (IFN) regulatory factors, and activator protein 1) are enriched in their regulatory regions [13].
Some HDAC families are also expected to exert anti-inflammatory effects. For instance, the expression of HDAC5 is suppressed by inflammatory cytokines (e.g., TNF-α and IL-1β), with concurrent induction of the expression of inflammatory mediators such as chemokines, suggesting that HDAC5 plays an anti-inflammatory role [14]. In contrast, inhibition of HDAC3 suppresses the expression of inflammatory mediators induced by IL-1β, indicating a proinflammatory function of HDAC3 [15]. This effect is not observed with HDAC1/2, HDAC6, and HDAC8 inhibition. In other words, there is a functional difference within the HDAC family, and development of RA treatment by selectively inhibiting proinflammatory HDACs is anticipated.
Recent studies on integrated analysis using cultured RASFs
Several studies have reported the integrated analysis of the transcriptome and epigenome of cultured RASFs. Whitaker et al. identified differentially expressed genes (DEGs) and differentially methylated genes (DMGs) by comparing unstimulated RASFs and OASFs, and compared the DEGs and DMGs with candidate RA susceptibility genes obtained from different databases (cross-sectional analysis) in a GWAS [16]. Among the expected novel therapeutic targets based on these three databases (e.g., Engulfment And Cell Motility 1 [ELMO1], LBH Regulator Of WNT Signaling Pathway [LBH], and Protein Tyrosine Phosphatase Non-Receptor Type 11 [PTPN11]), the knockdown of ELMO1, a protein-coding gene related to cell phagocytosis and movement, suppresses the migration and invasion activities of RASFs. Further, the knockdown of LBH, a protein-coding gene involved in embryonic development, increases the growth activity upstream of the LBH gene, and the enhancer activity is significantly lower in the risk allele than in the alternative allele. However, the enhancer activity is suppressed by methylation of the enhancer region in either alleles, and therefore, the RA risk SNPs and epigenomic modifications (DNA methylation) are thought to cooperatively control the enhancer activity of LBH [17]. PTPN11, a gene encoding Src homology region 2 domain-containing phosphatase-2, is highly expressed in RASFs compared with OASFs and is known to be related to the invasion ability of SFs [18]. The enhancer region in which the glucocorticoid receptor-binding motif of PTPN11 exists is hypermethylated and is involved in enhanced cell sensitivity to glucocorticoid and PTPN11 expression [19].
Furthermore, Ai et al. built a database consisting of transcriptomic and epigenomic information (histone modifications [H3K27ac, H3K4me1, H3K4me3, H3K36me3, H3K27me3, and H3K9me3], open chromatin, and DNA methylation) of unstimulated RASFs and OASFs. As a result of this multilayer analysis, genes related to Huntington’s disease are reported to be activated in RASFs in addition to known pathways, and invasion of RASFs is found to be suppressed by knockdown of the gene encoding Huntingtin Interacting Protein 1 (HIP1) [20].
SF subpopulation identification by single-cell analysis of the fresh synovium
The recently developed single-cell RNA sequencing technology has revolutionized SF research. Mizoguchi et al. reported that single-cell RNA sequencing of RA synovium classifies SFs into at least three subpopulations [21]. It has been suggested that CD34−THY1+ SFs localized around the blood vessels in the sublining represent a pathological subpopulation that produces inflammatory cytokines. Following this study, the Accelerating Medicines Partnership (AMP), which has been active since 2014 at the National Institutes of Health (NIH) in the USA, has conducted single-cell RNA sequencing and mass cytometry of the RA synovium and reported that CD34−THY1+HLA−DRhigh SFs form a pathological subpopulation that highly expresses IL-6 [4]. Moreover, a study on a mouse model suggested that FAPα+THY1+ SFs localize to the sublining and are involved in synovial inflammation, while the FAPα+THY1− SFs in the lining are involved in bone destruction [22]. Based on these findings, Wei et al. reported that the lining and sublining fibroblasts exist along a gradient that corresponds to the anatomical localization of SFs in the synovium, regulated by endothelium-derived Notch3 signaling [23].