Single-Cell Transcriptome Analysis Reveals the Importance of IRF1/FSTL1 in Synovial Fibroblast Subsets for the Development of Rheumatoid Arthritis

Abstract

Objectives. This study aimed to investigate the potential role of synovial fibroblasts (SFs) in the development of rheumatoid arthritis (RA) to identify potential molecular targets and provide a theoretical basis for the treatment of RA. Methods. GSE109449, a fibroblast transcriptome dataset of synovial tissue from RA and osteoarthritis (OA), were obtained from the GEO database. After standard cell quality control, this single-cell transcriptome data was used to perform routine single-cell analysis processes. After completing dimensionality reduction, clustering, and cell subset identification of fibroblasts, the SCENIC analysis helped calculate the significant gene regulatory networks in fibroblasts and their subsets. From these computed gene regulatory networks, the regulon in which follistatin-like protein 1 (FSTL1) resides was extracted and used to analyze the transcriptional regulatory status of fibroblasts. Finally, the gene set enrichment analysis (GSEA) was used to calculate the respective enriched gene sets of IRF1 and FSTL1. Results. Three SF subgroups were identified from the single-cell transcriptome analysis; SF subset 3 was more abundant in RA than in OA ( $p<0.001$ ). From the SCENIC analysis, we obtained 269 regulons and the corresponding gene regulatory networks in SF from the RA datasets. Next, we screened and obtained a regulon-containing FSTL1, where IRF1 was the major transcription factor. The top five regulons in SF subset 3 were TWIST1, MECOM, KLF6, MAFB, and RUNX1. Among the 3 SF subsets, IRF1 regulon was ranked the highest in SF subset 3. Differential analysis of pseudobulk RNA-seq showed that IRF1 was up-regulated in RA compared to OA. Between the three SF subgroups, IRF1 and FSTL1 expression was more up-regulated in SF subset 3 compared to the other two subgroups. Conclusions. IRF1 was found to regulate the invasiveness of SFs by regulating FSTL1, which may influence the disease progression of RA.