Identification of essential genes and immune cell infiltration in rheumatoid arthritis by bioinformatics analysis

Rheumatoid arthritis (RA) is a common autoimmune disease that can lead to severe joint damage and disability. And early diagnosis and treatment of RA can avert or substantially slow the progression of joint damage in up to 90% of patients, thereby preventing irreversible disability. Previous research indicated that 50% of the risk for the development of RA is attributable to genetic factors, but the pathogenesis is not well understood. Thus, it is urgent to identify biomarkers to arrest RA before joints are irreversibly damaged. Here, we first use the Robust Rank Aggregation method (RRA) to identify the differentially expressed genes (DEGs) between RA and normal samples by integrating four public RA patients’ mRNA expression data. Subsequently, these DEGs were used as the input for the weighted gene co-expression network analysis (WGCNA) approach to identify RA-related modules. The function enrichment analysis suggested that the RA-related modules were significantly enriched in immune-related actions. Then the hub genes were defined as the candidate genes. Our analysis showed that the expression levels of candidate genes were significantly associated with the RA immune microenvironment. And the results indicated that the expression of the candidate genes can use as predictors for RA. We hope that our method can provide a more convenient approach for the early diagnosis of RA.