Identification of five hub immune genes and characterization of two immune subtypes of osteoarthritis

Background: Osteoarthritis (OA) is one of the most prevalent chronic diseases, leading to degeneration of joints, chronic pain, and disability in the elderly. Little is known about the role of immune-related genes (IRGs) and immune cells in OA.

Method: Hub IRGs of OA were identified by differential expression analysis and filtered by three machine learning strategies, including random forest (RF), least absolute shrinkage and selection operator (LASSO), and support vector machine (SVM). A diagnostic nomogram model was then constructed by using these hub IRGs, with receiver operating characteristic (ROC) curve, decision curve analysis (DCA), and clinical impact curve analysis (CICA) estimating its performance and clinical impact. Hierarchical clustering analysis was then conducted by setting the hub IRGs as input information. Differences in immune cell infiltration and activities of immune pathways were revealed between different immune subtypes.

Result: Five hub IRGs of OA were identified, including TNFSF11, SCD1, PGF, EDNRB, and IL1R1. Of them, TNFSF11 and SCD1 contributed the most to the diagnostic nomogram model with area under the curve (AUC) values of 0.904 and 0.864, respectively. Two immune subtypes were characterized. The immune over-activated subtype showed excessively activated cellular immunity with a higher proportion of activated B cells and activated CD8 T cells. The two phenotypes were also seen in two validation cohorts.

Conclusion: The present study comprehensively investigated the role of immune genes and immune cells in OA. Five hub IRGs and two immune subtypes were identified. These findings will provide novel insights into the diagnosis and treatment of OA.