Abstract
BACKGROUND: osteoporosis is a skeletal disorder disease features low bone mass and poor bone architecture, which predisposes to increased risk of fracture. Copper death is a newly recognized form of cell death caused by excess copper ions, which presumably involve in various disease. Accordingly, we intended to investigate the molecular clusters related to the cuproptosis in osteoporosis and to construct a predictive model. METHODS: we investigated the expression patterns of cuproptosis regulators and immune signatures in osteoporosis based on the GSE56815 dataset. Through analysis of 40 osteoporosis samples, we investigated molecular clustering on the basis of cuproptosis--related genes, together with the associated immune cell infiltration. The WGCNA algorithm was applied to detect cluster-specific differentially expressed genes. Afterwards, the optimum machine model was selected by calculating the performance of the support vector machine model, random forest model, eXtreme Gradient Boosting and generalized linear model. Nomogram, decision curve analysis, calibration curves, and the GSE7158 dataset was utilizing to confirm the prediction efficiency. RESULTS: Differences between osteoporotic and non-osteoporotic controls confirm poorly adjusted copper death-related genes and triggered immune responses. In osteoporosis, two clusters of molecules in connection with copper death proliferation were outlined. The assessed levels of immune infiltration showed prominent heterogeneity between the different clusters. Cluster 2 was characterized by a raised immune score accompanied with relatively high levels of immune infiltration. The functional analysis we performed showed a close relationship between the different immune responses and specific differentially expressed genes in cluster 2. The random forest machine model showed the optimum discriminatory performance due to relatively low residuals and root mean square errors. Finally, a random forest model based on 5 genes was built, showing acceptable performance in an external validation dataset (AUC = 0.750). Calibration curve, Nomogram, and decision curve analyses also evinced fidelity in predicting subtypes of osteoporosis. CONCLUSION: Our study identifies the role of cuproptosis in OP and essentially illustrates the underlying molecular mechanisms that lead to OP heterogeneity.