Identification and Validation of Hub Genes for Predicting Treatment Targets and Immune Landscape in Rheumatoid...

Rheumatoid arthritis (RA) is recognized as a chronic inflammatory disease featured by pathological synovial inflammation. Currently, the underlying pathophysiological mechanisms of RA remain unclear. In the study, we attempted to explore the underlying mechanisms of RA and provide potential targets for the therapy of RA via bioinformatics analysis.

We downloaded four microarray datasets (GSE77298, GSE55235, GSE12021, and GSE55457) from the GEO database. Firstly, GSE77298 and GSE55457 were identified DEGs by the "limma" and "sva" packages of R software. Then, we performed GO, KEGG, and GSEA enrichment analyses to further analyze the function of DEGs. Hub genes were screened using LASSO analysis and SVM-RFE analysis. To further explore the differences of the expression of hub genes in healthy control and RA patient synovial tissues, we calculated the ROC curves and AUC. The expression levels of hub genes were verified in synovial tissues of normal and RA rats by qRT-PCR and western blot. Furthermore, the CIBERSORTx was implemented to assess the differences of infiltration in 22 immune cells between normal and RA synovial tissues. We explored the association between hub genes and infiltrating immune cells.

CRTAM, CXCL13, and LRRC15 were identified as RA's potential hub genes by machine learning and LASSO algorithms. In addition, we verified the expression levels of three hub genes in the synovial tissue of normal and RA rats by PCR and western blot. Moreover, immune cell infiltration analysis showed that plasma cells, T follicular helper cells, M0 macrophages, M1 macrophages, and gamma delta T cells may be engaged in the development and progression of RA.

In brief, our study identified and validated that three hub genes CRTAM, CXCL13, and LRRC15 might involve in the pathological development of RA, which could provide novel perspectives for the diagnosis and treatment with RA.