RESUMEN
We analyzed transcriptional data from 104 HPV+ (Human papillomavirus) HNSCC (head and neck squamous cell carcinoma) tumors together with two publicly available sources to identify highly robust transcriptional programs (modules) which could be detected consistently despite heterogeneous sequencing and quantification methodologies. Among 22 modules identified, we found a single module that naturally subclassifies HPV+ HNSCC tumors based on a bimodal pattern of gene expression, clusters all atypical features of HPV+ HNSCC biology into a single subclass, and predicts patient outcome in four independent cohorts. The subclass-defining gene set was strongly correlated with Nuclear factor kappa B (NF-κB) target expression. Tumors with high expression of this NF-κB module were rarely associated with activating PIK3CA alterations or viral integration, and also expressed higher levels of HPHPV E2 and had decreased APOBEC mutagenesis. Alternatively, they harbored inactivating alterations of key regulators of NF-κB, TNF receptor associated factor 3 (TRAF3), and cylindromatosis (CYLD), as well as retinoblastoma protein (RB1). HPV+ HNSCC cells in culture with experimental depletion of TRAF3 or CYLD displayed increased expression of the subclass-defining genes, as well as robust radio-sensitization, thus recapitulating both the tumor transcriptional state and improved treatment response observed in patient data. Across all gene sets investigated, methylation to expression correlations were the strongest for the subclass-defining, NF-κB-related genes. Increased tumor-infiltrating CD4+ T cells and increased Estrogen receptors alpha (ERα) expression were identified in NF-κB active tumors. Based on the relatively high rates of cure in HPV+ HNSCC, deintensification of therapy to reduce treatment-related morbidity is being studied at many institutions. Tumor subclassification based on oncogenic subtypes may help guide the selection of therapeutic intensity or modality for patients with HPV+ HNSCC.