Identifying hub circadian rhythm biomarkers and immune cell infiltration in rheumatoid arthritis

Rheumatoid arthritis (RA) is a chronic systemic autoimmune disease with symptoms characterized by typical circadian rhythmic changes. This study aimed to identify the hub circadian rhythm genes (CRGs) in RA and explore their association with immune cell infiltration and pathogenesis of RA.

The differentially expressed CRGs (DECRGs) between RA and normal control samples were screened from Datasets GSE12021 and GSE55235. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and Gene Set Enrichment Analysis were used to explore the potential functional mechanisms of DECRGs in RA. Weighted Gene Co-expression Network Analysis and Least Absolute Shrinkage and Selection Operator regression analysis were performed to identify hub CRGs of RA. CIBERSORT was conducted to compare the infiltration level of immune cells in RA and control synovial tissue and their relationship with hub genes. In addition, the diagnostic value of hub biomarkers was evaluated by the area under the receiver operator characteristic curve. Further, a nomogram prediction model was constructed and its significance for clinical decision-making was evaluated.

The green module was identified as the hub module associated with RA. Four hub CRGs (EGR1, FOSL2, GADD45B, and NFIL3) were identified and showed that they had the highest specificity and sensitivity for RA diagnosis, respectively. The expression levels and diagnostic values of these genes were externally validated in the dataset GSE55457. A nomogram prediction model based on the four hub CRGs was constructed and proved to have a certain clinical decision value. Additionally, the correlation analysis of immune cells with hub genes showed that all hub genes were significantly positively correlated with activated mast cells, resting memory CD4+ T cells, and monocytes. Whereas, all hub genes were negatively correlated with plasma cells, CD8+ T cells, and activated memory CD4+ T cells. Meanwhile, FOSL2 and GADD45B were negatively correlated with Tfh cells.

Four hub CRGs were identified and showed excellent diagnostic value for RA. These genes may be involved in the pathological process of RA by disrupting the rhythmic oscillations of cytokines through immune-related pathways and could be considered molecular targets for future chronotherapy against RA.