Construction of Molecular Subtype and Prognosis Prediction Model of Osteosarcoma Based on Aging-Related Genes

Abstract

Background. Osteosarcoma (OS) is a rare form of malignant bone cancer that is usually detected in young adults and adolescents. This disease shows a poor prognosis owing to its metastatic status and resistance to chemotherapy. Hence, it is necessary to design a risk model that can successfully forecast the OS prognosis in patients. Methods. The researchers retrieved the RNA sequencing data and follow-up clinical data related to OS patients from the TARGET and GEO databases, respectively. The coxph function in R software was used for carrying out the Univariate Cox regression analysis for deriving the aging-based genes related sto the OS prognosis. The researchers conducted consistency clustering using the ConcensusClusterPlus R package. The R software package ESTIMATE, MCPcounter, and GSVA packages were used for assessing the immune scores of various subtypes using the ssGSEA technique, respectively. The Univariate Cox and Lasso regression analyses were used for screening and developing a risk model. The ROC curves were constructed, using the pROC package. The performance of their developed risk model and designed survival curve was conducted, with the help of the Survminer package. Results. The OS patients were classified into 2 categories, as per the aging-related genes. The results revealed that the Cluster 1 patients showed a better prognosis than the Cluster 2 patients. Both clusters showed different immune microenvironments. Additional screening of the prognosis-associated genes revealed the presence of 5 genes, i.e., ERCC4, GPX4, EPS8, TERT, and STAT5A, and these data were used for developing the risk model. This risk model categorized the training set samples into the high- and low-risk groups. The patients classified into the high-risk group showed a poor OS prognosis compared to the low-risk patients. The researchers verified the reliability and robustness of the designed 5-gene signature using the internal and external datasets. This risk model was able to effectively predict the prognosis even in the samples having differing clinical features. Compared with other models, the 5- gene model performs better in predicting the risk of osteosarcoma. Conclusion. The 5-gene signature developed by the researchers in this study could be effectively used for forecasting the OS prognosis in patients.