Transcriptome Analysis Identifies Novel Prognostic Genes in Osteosarcoma

Osteosarcoma (OS), a malignant primary bone tumor often seen in young adults, is highly aggressive. The improvements in high-throughput technologies have accelerated the identification of various prognostic biomarkers for cancer survival prediction. However, only few studies focus on the prediction of prognosis in OS patients using gene expression data due to small sample size and the lack of public datasets. In the present study, the RNA-seq data of 82 OS samples, along with their clinical information, were collected from the TARGET database. To identify the prognostic genes for the OS survival prediction, we selected the top 50 genes of contribution as the initial candidate genes of the prognostic risk model, which were ranked by random forest model, and found that the prognostic model with five predictors including CD180, MYC, PROSER2, DNAI1, and FATE1 was the optimal multivariable Cox regression model. Moreover, based on a multivariable Cox regression model, we also developed a scoring method and stratified the OS patients into groups of different risks. The stratification for OS patients in the validation set further demonstrated that our model has a robust performance. In addition, we also investigated the biological function of differentially expressed genes between two risk groups and found that those genes were mainly involved with biological pathways and processes regarding immunity. In summary, the identification of novel prognostic biomarkers in OS would greatly assist the prediction of OS survival and development of molecularly targeted therapies, which in turn benefit patients’ survival.