In the HPV-negative subcohort (Fig. 5), the comparison of the two tumor clusters revealed a significant increase in the number of budding cases in the EMT-upregulated and myogenesis-upregulated clusters (P = 0.003 and P = 0.002, respectively) and as well as an increase in the number of NSD1-mutated cases in the EMT-downregulated and myogenesis-downregulated clusters (p = 0.003 and P = 0.03, respectively). This evidence concerns the gene NSD1 and neoplasm.