Identification of key genes and their correlation with immune infiltration in osteoarthritis using integrative bioinformatics approaches and machine-learning strategies

Osteoarthritis (OA) is a common degenerative joint disease and is closely associated with chronic, low-grade inflammation. Regulating ferroptosis by targeting ferroptosis-related genes may be a fast and effective way to delay the degeneration of OA. However, the molecular mechanisms and gene targets related to ferroptosis in OA are still unclear. Data of OA samples from 3 gene expression omnibus (GEO) datasets were combined to identify differentially expressed genes (DEGs). Ferroptosis-related genes (FRGs) retrieved by the Ferroptosis database were intersected with DEGs, and the intersected hub genes were used for functional enrichment analysis. The feature genes were obtained from the least absolute shrinkage and selection operator (LASSO) algorithm, support vector machine recursive feature elimination (SVM-RFE) algorithm, and random forest (RF) algorithm. Single sample gene set enrichment analysis (ssGSEA) was used to compare immune infiltration between OA patients and normal controls, and the correlation between feature genes and immune cells was analyzed. The expression levels of feature genes were confirmed by RT-PCR. In addition, to explore the applicability of these genes, we extended the bioinformatics analysis of these feature genes to cancer. Finally, 4 feature genes, GABARAPL1, TNFAIP3, ARNTL, and JUN, were confirmed in OA. Theirs expression level were validated by RT-PCR. ROC curves of the 4 genes exhibit excellent diagnostic efficiency for OA, suggesting that the 4 genes were associated with the pathogenesis of OA. Another GEO dataset validated this result. Further analysis revealed that the 4 feature genes were all closely related to the immune infiltration cells in OA. Additionally, results of prognosis analysis indicated that JUN might be a promising therapeutic target for cancer. GABARAPL1, TNFAIP3, ARNTL, and JUN may be predicted biomarkers for OA. The feature genes and association between feature genes and immune infiltration may provide potential biomarkers for OA prediction along with the better assessment of the disease.