Abstract:Purpose This study aimed to provide a comprehensive understanding of the genome-wide expression patterns in the synovial tissue samples of patients with rheumatoid arthritis (RA) to investigate the potential mechanisms regulating RA occurrence and development. Methods Transcription profiles of the synovial tissue samples from nine patients with RA and 15 patients with osteoarthritis (OA) (control) from the East Asian population were generated using RNA sequencing (RNA-seq). Gene set enrichment analysis (GSEA) was used to analyze all the detected genes and the differentially expressed genes (DEGs) were identified using DESeq. To further analyze the DEGs, the Gene Ontology (GO) functional enrichment and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed. The protein–protein interaction (PPI) network of the DEGs was constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and the hub genes were identified by topology clustering with the Molecular Complex Detection (MCODE)-Cytoscape. The most important hub genes were validated using quantitative real-time PCR (qRT-PCR). Results Of the 17,736 genes detected, 851 genes were identified as the DEGs (474 upregulated and 377 downregulated genes) using the false discovery rate (FDR) approach. GSEA revealed that the significantly enriched gene sets that positively correlated with RA were CD40 signaling overactivation, Th1 cytotoxic module, overactivation of the immune response, adaptive immune response, effective vs. memory CD8+ T cells (upregulated), and naïve vs. effective CD8+ T cells (downregulated). Biological process enrichment analysis showed that the DEGs were significantly enriched for signal transduction ( P = 3.01 × 10 −6 ), immune response ( P = 1.65 × 10 −24 ), and inflammatory response ( P = 5.76 × 10 −10 ). Molecule function enrichment analysis revealed that the DEGs were enriched in calcium ion binding ( P = 1.26 × 10 −5 ), receptor binding ( P = 1.26 × 10 −5 ), and cytokine activity ( P = 2.01 × 10 −3 ). Cellular component enrichment analysis revealed that the DEGs were significantly enriched in the plasma membrane ( P = 1.91 × 10 −31 ), an integral component of the membrane ( P = 7.39 × 10 −13 ), and extracellular region ( P = 7.63 × 10 −11 ). The KEGG pathway analysis showed that the DEGs were enriched in the cytokine–cytokine receptor interaction ( P = 3.05 × 10 −17 ), chemokine signaling ( P = 3.50 × 10 −7 ), T-cell receptor signaling ( P = 5.17 × 10 −4 ), and RA ( P = 5.17 × 10 −4 ) pathways. We confirmed that RA was correlated with the upregulation of the PPI network hub genes, such as CXCL13, CXCL6, CCR5, CXCR5, CCR2, CXCL3 , and CXCL10 , and the downregulation of the PPI network hub gene such as SSTR1 . Conclusion This study identified and validated the DEGs in the synovial tissue samples of patients with RA, which highlighted the activity of a subset of chemokine genes, thereby providing novel insights into the molecular mechanisms of RA pathogenesis and identifying potential diagnostic and therapeutic targets for RA.

Laser capture microscopy (LCM)‐RNAseq for topological mapping of synovial pathology during rheumatoid arthritis

Cellular deconstruction of inflamed synovium defines diverse inflammatory phenotypes in rheumatoid arthritis

Clonal associations between lymphocyte subsets and functional states in rheumatoid arthritis synovium

A Single-Cell RNA-Sequencing Analysis of Distinct Subsets of Synovial Macrophages in Rheumatoid Arthritis

Deconstruction of rheumatoid arthritis synovium defines inflammatory subtypes

Automated multi-scale computational pathotyping (AMSCP) of inflamed synovial tissue

Defining inflammatory cell states in rheumatoid arthritis joint synovial tissues by integrating single-cell transcriptomics and mass cytometry

New Insights From Single-Cell Sequencing Data: Synovial Fibroblasts and Synovial Macrophages in Rheumatoid Arthritis

Synovial fibroblast gene expression is associated with sensory nerve growth and pain in rheumatoid arthritis

Three-dimensional spatial transcriptomics uncovers cell type localizations in the human rheumatoid arthritis synovium

Identification and preliminary validation of synovial tissue-specific genes and their-mediated biological mechanisms in rheumatoid arthritis

Cell-Specific Gene Networks and Drivers in Rheumatoid Arthritis Synovial Tissues

Identification of Three Rheumatoid Arthritis Disease Subtypes by Machine Learning Integration of Synovial Histologic Features and RNA Sequencing Data

Single-cell Sequencing of Immune Cells from Anticitrullinated Peptide Antibody Positive and Negative Rheumatoid Arthritis.

RNA-Seq and Network Analysis Reveal Unique Chemokine Activity Gene Expression Signatures in Synovial Tissue From Rheumatoid Arthritis

Transcriptomic profiling and immuno‐histological characterization of synovial biopsies in juvenile idiopathic arthritis (JIA) suggest two broad categories, with significant overlap with rheumatoid arthritis (RA)

The role of lactylation in plasma cells and its impact on rheumatoid arthritis pathogenesis: insights from single-cell RNA sequencing and machine learning

RNA-seq and Network Analysis Reveal Unique Chemokine Activity Signatures in the Synovial Tissue of Patients With Rheumatoid Arthritis

Applications of single-cell RNA sequencing in rheumatoid arthritis

Synovial Transcriptome Profiling for Predicting Biological Treatment Response in Rheumatoid Arthritis: A Feasibility study