Abstract:Background: Rheumatoid arthritis (RA) is an autoimmune disease that affects individuals of all ages. The basic pathological manifestations are synovial inflammation, pannus formation, and erosion of articular cartilage, bone destruction will eventually lead to joint deformities and loss of function. However, the specific molecular mechanisms of synovitis tissue in RA are still unclear. Therefore, this study aimed to screen and explore the potential hub genes and immune cell infiltration in RA. Methods: Three microarray datasets (GSE12021, GSE55457, and GSE55235), from the Gene Expression Omnibus (GEO) database, have been analyzed to explore the potential hub genes and immune cell infiltration in RA. First, the LIMMA package was used to screen the differentially expression genes (DEGs) after removing the batch effect. Then the clusterProfiler package was used to perform functional enrichment analyses. Second, through weighted coexpression network analysis (WGCNA), the key module was identified in the coexpression network of the gene set. Third, the protein-protein interaction (PPI) network was constructed through STRING website and the module analysis was performed using Cytoscape software. Fourth, the CIBERSORT and ssGSEA algorithm were used to analyze the immune status of RA and healthy synovial tissue, and the associations between immune cell infiltration and RA-related diagnostic biomarkers were evaluated. Fifth, we used the quantitative reverse transcription-polymerase chain reaction (qRT-PCR) to validate the expression levels of the hub genes, and ROC curve analysis of hub genes for discriminating between RA and healthy tissue. Finally, the gene-drug interaction network was constructed using DrugCentral database, and identification of drug molecules based on hub genes using the Drug Signature Database (DSigDB) by Enrichr. Results: A total of 679 DEGs were identified, containing 270 downregulated genes and 409 upregulated genes. DEGs were primarily enriched in immune response and chemokine signaling pathways, according to functional enrichment analysis of DEGs. WGCNA explored the co-expression network of the gene set and identified key modules, the blue module was selected as the key module associated with RA. Seven hub genes are identified when PPI network and WGCNA core modules are intersected. Immune infiltration analysis using CIBERSORT and ssGSEA algorithms revealed that multiple types of immune infiltration were found to be upregulated in RA tissue compared to normal tissue. Furthermore, the levels of 7 hub genes were closely related to the relative proportions of multiple immune cells in RA. The results of the qRT-PCR demonstrated that the relative expression levels of 6 hub genes (CD27, LCK, CD2, GZMB, IL7R, and IL2RG) were up-regulated in RA synovial tissue, compared with normal tissue. Simultaneously, ROC curves indicated that the above 6 hub genes had strong biomarker potential for RA (AUC >0.8). Conclusions: Through bioinformatics analysis and qRT-PCR experiment, our study ultimately discovered 6 hub genes (CD27, LCK, CD2, GZMB, IL7R, and IL2RG) that closely related to RA. These findings may provide valuable direction for future RA clinical diagnosis, treatment, and associated research.

Integrative analysis highlighted susceptibility genes for rheumatoid arthritis

Transcriptome-wide association study identifies susceptibility genes for rheumatoid arthritis

Identification of new susceptibility loci associated with rheumatoid arthritis

Bioinformatics Analysis Identified the Hub Genes, Mrna–mirna–lncrna Axis, and Signaling Pathways Involved in Rheumatoid Arthritis Pathogenesis

Genetically Determined Circulating Levels of Cytokines and the Risk of Rheumatoid Arthritis

Exploring the causal relationship between inflammatory cytokines and inflammatory arthritis: A Mendelian randomization study

Integrative Analysis of Genome‐wide Association Study and Expression Quantitative Trait Loci Datasets Identified Various Immune Cell‐related Pathways for Rheumatoid Arthritis

Interleukins and rheumatoid arthritis: bi-directional Mendelian randomization investigation

Genome-wide Association Studies and Gene Expression Profiles of Rheumatoid Arthritis

Identification of shared disease marker genes and underlying mechanisms between rheumatoid arthritis and Crohn disease through bioinformatics analysis

Identification of hub genes in rheumatoid arthritis through an integrated bioinformatics approach

Comprehensive Bioinformatics Analysis Reveals Hub Genes and Inflammation State of Rheumatoid Arthritis

Genetic evidence reveals a causal relationship between rheumatoid arthritis and interstitial lung disease

Influence of MIF, CD40, and CD226 polymorphisms on risk of rheumatoid arthritis

Synovial transcriptome-wide association study implicates novel genes underlying rheumatoid arthritis risk

Genome-wide identification of RNA modification-related single nucleotide polymorphisms associated with rheumatoid arthritis

Assessment of genetic polymorphisms within nuclear factor-κB signaling pathway genes in rheumatoid arthritis: Evidence for replication and genetic interaction

Experimental Immunology the Association of Single Nucleotide Polymorphisms in TNFAIP3 with Rheumatoid Arthritis in the Chinese Population

Identification of Pathogenic Genes Related to Rheumatoid Arthritis Through Integrated Analysis of DNA Methylation and Gene Expression Profiling.

Assessing the causal role of immune traits in rheumatoid arthritis by bidirectional Mendelian randomization analysis

Screening and identification of potential hub genes and immune cell infiltration in the synovial tissue of rheumatoid arthritis by bioinformatic approach