Increasing evidence has shown that hypoxia microenvironment relates to tumor initiation and progression. However, no studies focus on the application of hypoxia-associated genes in predicting osteosarcoma patients’ prognosis. This research aims to identify the hypoxia-associated genes related to osteosarcoma metastasis and construct a gene signature to predict osteosarcoma prognosis.
The differentially expressed messenger RNAs (DEmRNAs) related to osteosarcoma metastasis were identified from Therapeutically Applicable Research to Generate Effective Treatments (Target) database. Univariate and multivariate cox regression analyses were performed to develop the hypoxia-associated prognostic signature. The Kaplan–Meier (KM) survival analyses of patients with high and low hypoxia risk scores were conducted. The nomogram was constructed and the gene signature was validated in the external Gene Expression Omnibus (GEO) cohort. Single-sample gene set enrichment analysis (ssGSEA) was conducted to investigate the relationships between immune infiltration and gene signature.
Two genes, including decorin (DCN) and prolyl 4-hydroxylase subunit alpha 1 (P4HA1), were involved in the hypoxia-associated gene signature. In training and testing datasets, patients with high-risk scores showed lower survival rates and the gene signature was identified as the independent prognostic factor. Receiver operating characteristic (ROC) curves demonstrated the robustness of signature. Functional analyses of DEmRNAs among high- and low-risk groups revealed that immune-associated functions and pathways were significantly enriched. Furthermore, ssGSEA showed that five immune cells (DCs, macrophages, neutrophils, pDCs, and TIL) and three immune features (CCR, APC co inhibition, and Check-point) were down-regulated in the high-risk group.
The current study established and validated a novel hypoxia-associated gene signature in osteosarcoma. It could act as a prognostic biomarker and serve as therapeutic guidance in clinical applications.