A 13-gene expression-based radioresistance score highlights the heterogeneity in the response to radiation therapy across HPV-negative HNSCC molecular subtypes

Jean-Philippe Foy, Louis Bazire, Sandra Ortiz-Cuaran, Sophie Deneuve, Janice Kielbassa, Emilie Thomas, Alain Viari, Alain Puisieux, Patrick Goudot, Chloé Bertolus, Nicolas Foray, Youlia Kirova, Pierre Verrelle, Pierre Saintigny

Radiotherapy for head and neck squamous cell carcinomas (HNSCC) is associated with a substantial morbidity and inconsistent efficacy. Human papillomavirus (HPV)-positive status is recognized as a marker of increased radiosensitivity. Our goal was to identify molecular markers associated with benefit to radiotherapy in patients with HPV-negative disease. Gene expression profiles from public repositories were downloaded for data mining. Training sets included 421 HPV-negative HNSCC tumors from The Cancer Genome Atlas (TCGA) and 32 HNSCC cell lines with available radiosensitivity data (GSE79368). A radioresistance (RadR) score was computed using the single sample Gene Set Enrichment Analysis tool. The validation sets included two panels of cell lines (NCI-60 and GSE21644) and HPV-negative HNSCC tumor datasets, including 44 (GSE6631), 82 (GSE39366), and 179 (GSE65858) patients, respectively. We finally performed an integrated analysis of the RadR score with known recurrent genomic alterations in HNSCC, patterns of protein expression, biological hallmarks, and patterns of drug sensitivity using TCGA and the E-MTAB-3610 dataset (659 pancancer cell lines, 140 drugs). We identified 13 genes differentially expressed between tumor and normal head and neck mucosa that were associated with radioresistance in vitro and in patients. The 13-gene expression-based RadR score was associated with recurrence in patients treated with surgery and adjuvant radiotherapy but not with surgery alone. It was significantly different among different molecular subtypes of HPV-negative HNSCC and was significantly lower in the "atypical" molecular subtype. An integrated analysis of RadR score with genomic alterations, protein expression, biological hallmarks and patterns of drug sensitivity showed a significant association with CCND1 amplification, fibronectin expression, seven hallmarks (including epithelial-to-mesenchymal transition and unfolded protein response), and increased sensitivity to elesclomol, an HSP90 inhibitor. Our study highlights the clinical relevance of the molecular classification of HNSCC and the RadR score to refine radiation strategies in HPV-negative disease.