Fig. 4

A Consensus matrixes for 770 HNSCC samples based on the DEGs among the three hypoxia clusters. Three subtypes were determined for all samples (k = 3). B OS differences among the three gene clusters. C Expression of hypoxia-related DEGs in different hypoxia clusters and gene clusters. D–E OS differences between HHSG and LHSG using a combination of TCGA-HNSCC and GSE65858, GSE41613 data. F Expression of different HRG in normal and tumor tissues (^nsP > 0.05, *P < 0.05, ** P < 0.01, *** P < 0.001). G Sankey diagram showing the distribution of patients with HPV16, hypoxia clusters, gene clusters, HS, and survival status. H Correlations between HS, immune-related functions, and the abundance of each immune cell in 770 HNSCC samples. I Differences in the HS among the three hypoxia clusters. J Differences in the HS among the three gene clusters. HNSCC head and neck squamous cell carcinoma, DEGs differentially expressed genes, OS overall survival, TCGA The Cancer Genome Atlas, HRG hypoxia-related gene, HHSG high hypoxia score group, LHSG low hypoxia score group, HS hypoxia score