Identification and verification of a BMPs-related gene signature for osteosarcoma prognosis prediction

Abstract Background This study aimed to get a deeper insight into new osteosarcoma (OS) signature based on bone morphogenetic proteins (BMPs)-related genes and to confirm the prognostic pattern to speculate on the overall survival among OS patients. Methods Firstly, pathway analyses using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) were managed to search for possible prognostic mechanisms attached to the OS-specific differentially expressed BMPs-related genes (DEBRGs). Secondly, univariate and multivariate Cox analysis was executed to filter the prognostic DEBRGs and establish the polygenic model for risk prediction in OS patients with the least absolute shrinkage and selection operator (LASSO) regression analysis. The receiver operating characteristic (ROC) curve weighed the model’s accuracy. Thirdly, the GEO database (GSE21257) was operated for independent validation. The nomogram was initiated using multivariable Cox regression. Immune infiltration of the OS sample was calculated. Finally, the three discovered hallmark genes’ mRNA and protein expressions were verified. Results A total of 46 DEBRGs were found in the OS and control samples, and three prognostic DEBRGs (DLX2, TERT, and EVX1) were screened under the LASSO regression analyses. Multivariate and univariate Cox regression analysis were devised to forge the OS risk model. Both the TARGET training and validation sets indicated that the prognostic biomarker-based risk score model performed well based on ROC curves. In high- and low-risk groups, immune cells, including memory B, activated mast, resting mast, plasma, and activated memory CD4 + T cells, and the immune, stromal, and ESTIMATE scores showed significant differences. The nomogram that predicts survival was established with good performance according to clinical features of OS patients and risk scores. Finally, the expression of three crucial BMP-related genes in OS cell lines was investigated using quantitative real-time polymerase chain reaction (qRT-PCR) and western blotting (WB). Conclusion The new BMP-related prognostic signature linked to OS can be a new tool to identify biomarkers to detect the disease early and a potential candidate to better treat OS in the future.