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.

RNA-seq Analysis Reveals Different Gene Ontologies and Pathways in Rheumatoid Arthritis and Kashin-Beck Disease.

Differential Gene Expression in Articular Cartilage Between Rheumatoid Arthritis and Endemic Kashin–Beck Disease

Identification of the key genes and long non-coding RNAs in ankylosing spondylitis using RNA sequencing (vol 43, pg 1179, 2019)

Comprehensive comparative analysis of histopathology and gene expression in subchondral bone between kashin-beck disease and primary osteoarthritis

Identification of differentially expressed genes and pathways between primary osteoarthritis and endemic osteoarthritis (Kashin-Beck disease).

Comprehensive expression profiles of mRNAs, lncRNAs and miRNAs in Kashin-Beck disease identified by RNA-sequencing

Comparative Analysis of Signaling Pathways in Peripheral Blood from Patients with Kashin-Beck Disease and Osteoarthritis.

Whole-Transcriptome Sequencing of Knee Joint Cartilage from Kashin–Beck Disease and Osteoarthritis Patients

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

Genome-Wide Differentially Methylated Region Analysis to Reveal Epigenetic Differences of Articular Cartilage in Kashin–Beck Disease and Osteoarthritis

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

Comparison of microRNA expression profiles of Kashin-Beck disease, osteoarthritis and rheumatoid arthritis

Transcriptome‐wide RNA m6A methylation profiles in an endemic osteoarthropathy, Kashin‐Beck disease

Regulatory gene networks and signaling pathways from primary osteoarthritis and Kashin-Beck disease, an endemic osteoarthritis, identified by three analysis software.

Genome‐wide Study Identifies the Regulatory Gene Networks and Signaling Pathways from Chondrocyte and Peripheral Blood Monocyte of Kashin–Beck Disease

Difference in Apoptosis-Associated Genes Expression Profiling and Immunohistology Analysis Between Kashin-Beck Disease and Primary Osteoarthritis

ATAC-seq Reveals the Roles of Chromatin Accessibility in the Chondrocytes of Kashin–Beck Disease Compared with Primary Osteoarthritis

Identified molecular mechanism of interaction between environmental risk factors and differential expression genes in cartilage of Kashin-Beck disease.

Proteomic analysis of knee cartilage reveals potential signaling pathways in pathological mechanism of Kashin-Beck disease compared with osteoarthritis

Comparative Analysis of Gene Expression Profiles Between Primary Knee Osteoarthritis and an Osteoarthritis Endemic to Northwestern China, Kashin-Beck Disease

Identification of differently expressed genes and pathways in cartilage of Kashin-Beck disease using microarray