Abstract
Purpose: Tumor markers that are related to hypoxia, proliferation, DNA damage repair and stem cell-ness, have a prognostic value in advanced stage HNSCC patients when assessed individually. Here we aimed to evaluate and validate this in a multifactorial context and assess interrelation and the combined role of these biological factors in determining chemo-radiotherapy response in HPV-negative advanced HNSCC. Methods: RNA sequencing data of pre-treatment biopsy material from 197 HPV-negative advanced stage HNSCC patients treated with definitive chemoradiotherapy was analyzed. Biological parameter scores were assigned to patient samples using previously generated and described gene expression signatures. Locoregional control rates were used to assess the role of these biological parameters in radiation response and compared to distant metastasis data. Biological factors were ranked according to their clinical impact using bootstrapping methods and multivariate Cox regression analyses that included clinical variables. Multivariate Cox regression analyses comprising all biological variables were used to define their relative role among all factors when combined. Results: Only few biomarker scores correlate with each other, underscoring their independence. The different biological factors do not correlate or cluster, except for the two stem cell markers CD44 and SLC3A2 (r = 0.4, p < 0.001) and acute hypoxia prediction scores which correlated with T-cell infiltration score, CD8+ T cell abundance and proliferation scores (r = 0.52, 0.56, and 0.6, respectively with p < 0.001). Locoregional control association analyses revealed that chronic (Hazard Ratio (HR) = 3.9) and acute hypoxia (HR = 1.9), followed by stem cell-ness (CD44/SLC3A2; HR = 2.2/2.3), were the strongest and most robust determinants of radiation response. Furthermore, multivariable analysis, considering other biological and clinical factors, reveal a significant role for EGFR expression (HR = 2.9, p < 0.05) and T-cell infiltration (CD8+T-cells: HR = 2.2, p < 0.05; CD8+T-cells/Treg: HR = 2.6, p < 0.01) signatures in locoregional control of chemoradiotherapy-treated HNSCC. Conclusion: Tumor acute and chronic hypoxia, stem cell-ness, and CD8+ T-cell parameters are relevant and largely independent biological factors that together contribute to locoregional control. The combined analyses illustrate the additive value of multifactorial analyses and support a role for EGFR expression analysis and immune cell markers in addition to previously validated biomarkers. This external validation underscores the relevance of biological factors in determining chemoradiotherapy outcome in HNSCC.