Bioinformatics-based study to identify immune infiltration and inflammatory-related hub genes as biomarkers for the treatment of rheumatoid arthritis.

Abstract

Rheumatoid arthritis (RA) is a systemic autoimmune disease whose principal pathological change is aggressive chronic synovial inflammation; however, the specific etiology and pathogenesis have not been fully elucidated. We downloaded the synovial tissue gene expression profiles of four human knees from the Gene Expression Omnibus database, analyzed the differentially expressed genes in the normal and RA groups, and assessed their enrichment in functions and pathways using bioinformatics methods and the STRING online database to establish protein-protein interaction networks. Cytoscape software was used to obtain 10 hub genes; receiver operating characteristic (ROC) curves were calculated for each hub gene and differential expression analysis of the two groups of hub genes. The CIBERSORT algorithm was used to impute immune infiltration. We identified the signaling pathways that play important roles in RA and 10 hub genes: Ccr1, Ccr2, Ccr5, Ccr7, Cxcl5, Cxcl6, Cxcl13, Ccl13, Adcy2, and Pnoc. The diagnostic value of these 10 hub genes for RA was confirmed using ROC curves and expression analysis. Adcy2, Cxcl13, and Ccr5 are strongly associated with RA development. The study also revealed that the differential infiltration profile of different inflammatory immune cells in the synovial tissue of RA is an extremely critical factor in RA progression. This study may contribute to the understanding of signaling pathways and biological processes associated with RA and the role of inflammatory immune infiltration in the pathogenesis of RA. In addition, this study shows that Adcy2, Cxcl13, and Ccr5 have the potential to be biomarkers for RA treatment.