Identification of ferroptosis-related molecular clusters and genes for diabetic osteoporosis based on the machine learning

Abstract

Background: Diabetic osteoporosis exhibits heterogeneity at the molecular level. Ferroptosis, a controlled form of cell death brought on by a buildup of lipid peroxidation, contributes to the onset and development of several illnesses. The aim was to explore the molecular subtypes associated with ferroptosis in diabetic osteoporosis at the molecular level and to further elucidate the potential molecular mechanisms. Methods: Integrating the CTD, GeneCards, FerrDb databases, and the microarray data of GSE35958, we identified ferroptosis-related genes (FRGs) associated with diabetic osteoporosis. We applied unsupervised cluster analysis to divide the 42 osteoporosis samples from the GSE56814 microarray data into different subclusters based on FRGs. Subsequently, FRGs associated with two ferroptosis subclusters were obtained by combining database genes, module-related genes of WGCNA, and differentially expressed genes (DEGs). Eventually, the key genes from FRGs associated with diabetic osteoporosis were identified using the least absolute shrinkage and selection operator (LASSO), Boruta, support vector machine recursive feature elimination (SVM ­ RFE), and extreme gradient boosting (XGBoost) machine learning algorithms. Based on ROC curves of external datasets (GSE56815), the model’s efficiency was examined. Results: We identified 15 differentially expressed FRGs associated with diabetic osteoporosis. In osteoporosis, two distinct molecular clusters related to ferroptosis were found. The expression results and GSVA analysis indicated that 15 FRGs exhibited significantly different biological functions and pathway activities in the two ferroptosis subclusters. Therefore, we further identified 17 FRGs associated with diabetic osteoporosis between the two subclusters. The results of the comprehensive analysis of 17 FRGs demonstrated that these genes were heterogeneous and had a specific interaction between the two subclusters. Ultimately, the prediction model had a strong foundation and excellent AUC values (0.84 for LASSO, 0.84 for SVM ­ RFE, 0.82 for Boruta, and 0.81 for XGBoost). IDH1 is a common gene to all four algorithms thus being identified as a key gene with a high AUC value (AUC = 0.698). Conclusions: As a ferroptosis regulator, IDH1 is able to distinguish between distinct molecular subtypes of diabetic osteoporosis, which may offer fresh perspectives on the pathogenesis of the disease’s clinical symptoms and prognostic heterogeneity.