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

Abstract

Rheumatoid arthritis (RA) is an autoimmune disease characterized by erosive arthritis, which has not been thoroughly cured yet, and standardized treatment is helpful for alleviating clinical symptoms. Here, various bioinformatics analysis tools were comprehensively utilized, aiming to identify critical biomarkers and possible pathogenesis of RA. Three gene expression datasets profiled by microarray were obtained from GEO database. Dataset GSE55235 and GSE55457 were merged for subsequent analyses. We identified differentially expressed genes (DEGs) in RStudio with limma package, performing functional enrichment analysis based on GSEA software and clusterProfiler package. Next, protein-protein interaction (PPI) network was set up through STRING database and Cytoscape. Moreover, CIBERSORT website was used to assess the inflammatory state of RA. Finally, we validated the candidate hub genes with dataset GSE77298. As a result, we identified 106 DEGs (72 upregulated and 34 downregulated genes). Through GO, KEGG, and GSEA analysis, we found that DEGs were mainly involved in immune response and inflammatory signaling pathway. With the help of Cytoscape software and MCODE plug-in, the most prominent subnetwork was screened out, containing 14 genes and 45 edges. For ROC curve analysis, eight genes with AUC >0.80 were considered as hub genes of RA. In conclusion, compared with healthy controls, the DEGs and their closely related biological functions were analyzed, and we held that chemokines and immune cells infiltration promote the progression of rheumatoid arthritis. Targeting the eight biomarkers we identified may be useful for the diagnosis and treatment of rheumatoid arthritis.