Machine learning and bioinformatics analysis to identify and validate diagnostic model associated with immune infiltration in rheumatoid arthritis.

BACKGROUND

Rheumatoid arthritis (RA) is a chronic autoinflammatory condition that can result in significant disability. This study focuses on identifying immune infiltration-related diagnostic biomarkers of RA patients.

METHOD: Publicly available datasets from the Gene Expression Omnibus (GEO) were analyzed using ssGSEA and CIBERSORT algorithms to measure immune cell subset infiltration. Functional enrichment analyses, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), were performed. Additionally, least absolute shrinkage and selection operator (LASSO) regression and machine learning methods, such as random forest, were employed to identify key immune infiltration-related genes. Differential expression of these hub genes between subgroups was compared, and their diagnostic potential was evaluated through receiver operating characteristic (ROC) analysis, validated using GSE93777 and GSE205962 datasets.

RESULTS

Analysis of mRNA expression from GSE93272 revealed two distinct clusters: immunity_low (38 samples) and immunity_high (194 samples). A total of 320 differentially expressed genes (DEGs) were identified by intersecting DEGs from these clusters with those from RA and healthy controls (HC). Five hub genes (BMX, BTLA, CENPK, CMPK2, GBP3) were selected using LASSO and machine learning approaches, forming the basis of a diagnostic risk model. This five-gene model demonstrated strong diagnostic performance for distinguishing immune infiltration statuses (AUC = 0.977) and identifying RA patients (AUC = 0.942). External validation with GSE93777 (AUC = 0.807) and GSE205962 (AUC = 0.938) datasets confirmed its reliability.

CONCLUSION

Five key genes associated with immune infiltration were identified, enabling the construction of a diagnostic model for RA. This model shows potential to improve RA diagnosis and facilitate the development of personalized therapeutic strategies. Key Points •RA patients were stratified into two distinct immune subtypes (Immunity_H and Immunity_L) based on ssGSEA analysis of 29 immune gene sets, revealing marked differences in immune activity and HLA gene expression. •Five hub genes including BMX, BTLA, CENPK, CMPK2, and GBP3, were identified through LASSO and Random Forest algorithms, forming a robust risk model that accurately distinguishes RA patients and their immune subtypes. •The predictive model was validated in two independent external cohorts, confirming its diagnostic reliability and generalizability across RA datasets.