Integrated Analysis of Multi-Omics Data to Establish a Hypoxia-Related Prognostic Model in Osteosarcoma

Osteosarcoma (OS) is the most common malignant bone tumor in clinical practice, and currently, the ability to predict prognosis in the diagnosis of OS is limited. There is an urgent need to find new diagnostic methods and treatment strategies for OS.

We downloaded the multi-omics data for OS from the TARGET database. Prognosis-associated methylation sites were used to identify clustered subtypes of OS, and OS was classified into 3 subtypes (C1, C2, C3). Survival analysis showed significant differences between the C3 subtype and the other subtypes. Subsequently, differentially expressed genes (DEGs) across subtypes were screened and subjected to pathway enrichment analysis.

A total of 249 DEGs were screened from C3 subtype to other subtypes. Metabolic pathway enrichment analysis showed that DEGs were significantly enriched to the hypoxic pathway. Based on univariate and multivariate COX regression analysis, 12 genes from the hypoxia pathway were further screened and used to construct hypoxia-related prognostic model (HRPM). External validation of the HRPM was performed on the GSE21257 dataset. Finally, differences in survival and immune infiltration between high and low risk score groups were compared.

In summary, we proposed a hypoxia-associated risk model based on a 12-gene expression signature, which is potentially valuable for prognostic diagnosis of patients with OS.