Development and Verification of a Hypoxic Gene Signature for Predicting Prognosis, Immune Microenvironment, and Chemosensitivity for Osteosarcoma

Objective: Hypoxic tumors contribute to local failure and distant metastases. Nevertheless, the molecular hallmarks of hypoxia remain ill-defined in osteosarcoma. Here, we developed a hypoxic gene signature in osteosarcoma prognoses. Methods: With the random survival forest algorithm, a prognostic hypoxia-related gene signature was constructed for osteosarcoma in the TARGET cohort. Overall survival (OS) analysis, receiver operating characteristic (ROC) curve, multivariate cox regression analysis, and subgroup analysis were utilized for assessing the predictive efficacy of this signature. Also, external validation was presented in the GSE21257 cohort. GSEA was applied for signaling pathways involved in the high- and low-risk samples. Correlation analyses between risk score and immune cells, stromal/immune score, immune checkpoints, and sensitivity of chemotherapy drugs were performed in osteosarcoma. Then, a nomogram was built by integrating risk score, age, and gender. Results: A five-hypoxic gene signature was developed for predicting survival outcomes of osteosarcoma patients. ROC curves confirmed that this signature possessed the well predictive performance on osteosarcoma prognosis. Furthermore, it could be independently predictive of prognosis. Metabolism of xenobiotics by cytochrome P450 and nitrogen metabolism were activated in the high-risk samples while cell adhesion molecules cams and intestinal immune network for IgA production were enriched in the low-risk samples. The low-risk samples were characterized by elevated immune cell infiltrations, stromal/immune scores, TNFRSF4 expression, and sensitivity to cisplatin. The nomogram accurately predicted 1-, 3-, and 5-years survival duration. Conclusion: These findings might offer an insight into the optimization of prognosis risk stratification and individualized therapy for osteosarcoma patients.