Construction and Validation of a Prognostic Model for Osteosarcoma Patients Based on Autophagy-Related Genes

Background Osteosarcoma is the most frequent primary bone malignancy with a poor prognosis because of pulmonary metastasis. Autophagy is strongly associated with tumor metastasis, and it is valuable to construct an autophagy-related gene risk model for predicting the prognosis of osteosarcoma. Methods We obtained ARGs from the Human Autophagy Database and RNA-sequencing data of osteosarcoma patients from the Gene Expression Omnibus (GEO) database. Subsequently, univariate and multivariate cox regression analyses were performed to construct a three-gene prognostic model and its accuracy was further confirmed in the Therapeutic Applications Research to Generate Effective Treatments (TARGET) database. Afterward, we detected the expression levels and effects on osteosarcoma cells metastasis of MYC and MBTPS2, which were involved in the model. Results In both training and verification cohorts, patients with lower risk scores had longer OS, and the model was identified as an independent prognostic factor in osteosarcoma. Besides, the ROC curve demonstrated the reliability of the model. Furthermore, RT-qPCR, Western Blotting and IHC results indicated that MYC and MBTPS2 were differently expressed in osteosarcoma tissues and cell lines. MYC knockdown or MBTPS2 overexpression prevented the capacity of migration and invasion in osteosarcoma cell lines through inhibiting cellular autophagy. Conclusion The risk model based on three ARGs had a strong ability to predict the prognosis of osteosarcoma patients. Our findings also suggested that MYC and MBTPS2 were two major factors regulating autophagy in osteosarcoma, and could serve as potential therapeutic targets for osteosarcoma.