Establishment of a Plasticity-Associated Risk Model Based on a SOX2- and SOX9-Related Gene Set in Head and Neck Squamous Cell Carcinoma

Abstract

Recent studies highlighted SOX2 and SOX9 as key determinants for cancer-cell plasticity and demonstrated that cisplatin-induced adaptation in oral squamous cell carcinoma (SCC) is acquired by an inverse regulation of both transcription factors. However, the association between SOX2/SOX9-related genetic programs with risk factors and genetic or epigenetic alterations in primary head and neck SCC (HNSCC), and their prognostic value is largely unknown. Here, we identified differentially-expressed genes (DEG) related to SOX2 and SOX9 transcription in The Cancer Genome Atlas (TCGA)-HNSC, which enable clustering of patients into groups with distinct clinical features and survival. A prognostic risk model was established by LASSO Cox regression based on expression patterns of DEGs in TCGA-HNSC (training cohort), and was confirmed in independent HNSCC validation cohorts as well as other cancer cohorts from TCGA. Differences in the mutational landscape among risk groups of TCGA-HNSC demonstrated an enrichment of truncating NSD1 mutations for the low-risk group and elucidated DNA methylation as modulator of SOX2 expression. Gene set variation analysis (GSVA) revealed differences in several oncogenic pathways among risk groups, including upregulation of gene sets related to oncogenic KRAS signaling for the high-risk group. Finally, in silico drug screen analysis revealed numerous compounds targeting EGFR signaling with significantly lower efficacy for cancer cell lines with a higher risk phenotype, but also indicated potential vulnerabilities. Implications: The established risk model identifies patients with primary HNSCC, but also other cancers at a higher risk for treatment failure, who might benefit from a therapy targeting SOX2/SOX9-related gene regulatory and signaling networks.