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 of Genome‐wide Association Study and Expression Quantitative Trait Loci Datasets Identified Various Immune Cell‐related Pathways for Rheumatoid Arthritis

Integrative analysis of genome-wide association studies and gene expression analysis identifies pathways associated with rheumatoid arthritis.

Integrative Analyses for Functional Mechanisms Underlying Associations for Rheumatoid Arthritis

Transcriptome-wide association study identifies susceptibility genes for rheumatoid arthritis

Integrated analysis of single-cell RNA-seq, bulk RNA-seq, Mendelian randomization, and eQTL reveals T cell-related nomogram model and subtype classification in rheumatoid arthritis

Integrative Analysis for Identification of Shared Markers from Various Functional Cells/tissues for Rheumatoid Arthritis.

Identification of new susceptibility loci associated with rheumatoid arthritis

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

Identifying Immune Cell Infiltration and Effective Diagnostic Biomarkers in Rheumatoid Arthritis by Bioinformatics Analysis

The shared biomarkers and immune landscape in psoriatic arthritis and rheumatoid arthritis: Findings based on bioinformatics, machine learning and single-cell analysis

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

Potential Biological Molecules and Immune Cell for Rheumatoid Arthritis Synovial Lesions

Autophagy-Related Genes Associated with Immune Cell Infiltration in Rheumatoid Arthritis Identified by Integrated Weighted Gene Co-Expression Network

Analyzation of the Peripheral Blood Mononuclear Cells Atlas and Cell Communication of Rheumatoid Arthritis Patients Based on Single-Cell RNA-Seq

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

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

Identification of Critical Biomarkers and Immune Infiltration in Rheumatoid Arthritis Based on WGCNA and LASSO Algorithm

The Commonness in Immune Infiltration of Rheumatoid Arthritis and Atherosclerosis: Screening for Central Targets Via Microarray Data Analysis.

Integrative analysis highlighted susceptibility genes for rheumatoid arthritis

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

Identification of Diagnostic Signatures and Immune Cell Infiltration Characteristics in Rheumatoid Arthritis by Integrating Bioinformatic Analysis and Machine-Learning Strategies