Perineural invasion is a prevalent pathological finding in head and neck squamous cell carcinoma and a risk factor for unfavorable survival. An adequate diagnosis of perineural invasion by pathologic examination is limited due to the availability of tumor samples from surgical resection, which can arise in cases of definitive nonsurgical treatment. To address this medical need, we established a random forest prediction model for the risk assessment of perineural invasion, including occult perineural invasion, and characterized distinct cellular and molecular features based on our new and extended classification. RNA sequencing data of head and neck squamous cell carcinoma from The Cancer Genome Atlas were used as a training cohort to identify differentially expressed genes that are associated with perineural invasion. A random forest classification model was established based on these differentially expressed genes and was validated by inspection of H&E-stained whole image slides. Differences in epigenetic regulation and the mutational landscape were detected by an integrative analysis of multiomics data and single-cell RNA-sequencing data were analyzed. We identified a 44-gene expression signature related to perineural invasion and enriched for genes mainly expressed in cancer cells according to single-cell RNA-sequencing data. A machine learning model was trained based on the expression pattern of the 44-gene set with the unique feature to predict occult perineural invasion. This extended classification model enabled a more accurate analysis of alterations in the mutational landscape and epigenetic regulation by DNA methylation as well as quantitative and qualitative differences in the cellular composition in the tumor microenvironment between head and neck squamous cell carcinoma with or without perineural invasion. In conclusion, the newly established model could not only complement histopathologic examination as an additional diagnostic tool but also guide the identification of new drug targets for therapeutic intervention in future clinical trials with head and neck squamous cell carcinoma patients at a higher risk for treatment failure due to perineural invasion.

Recent studies have highlighted neurons and their associated Schwann cells (SCs) as key regulators of cancer development. However, the mode of their interaction with tumor cells or other components of the tumor microenvironment (TME) remains elusive. We established an SC-related 43-gene set as a surrogate for peripheral nerves in the TME. Head and neck squamous cell carcinoma (HNSCC) from The Cancer Genome Atlas (TCGA) were classified into low, intermediate and high SC score groups based on the expression of this gene set. Perineural invasion (PNI) and TGF-β signaling were hallmarks of SC^high tumors, whereas SC^low tumors were enriched for HPV16-positive OPSCC and higher PI3K-MTOR activity. The latter activity was partially explained by a higher frequency of PTEN mutation and PIK3CA copy number gain. The inverse association between PI3K-MTOR activity and peripheral nerve abundance was context-dependent and influenced by the TP53 mutation status. An in silico drug screening approach highlighted the potential vulnerabilities of HNSCC with variable SC scores and predicted a higher sensitivity of SC^low tumors to DNA topoisomerase inhibitors. In conclusion, we have established a tool for assessing peripheral nerve abundance in the TME and provided new clinical and biological insights into their regulation. This knowledge may pave the way for new therapeutic strategies and impart proof of concept in appropriate preclinical models.

Purpose: Identification of molecularly-defined cancer subgroups and targeting tumor-specific vulnerabilities have a strong potential to improve treatment response and patient outcomes but remain an unmet challenge of high clinical relevance, especially in head and neck squamous cell carcinoma (HNSC). Experimental design: We established a UCHL1-related gene set to identify and molecularly characterize a UCHL1-related subgroup within TCGA-HNSC by integrative analysis of multi-omics data. An extreme gradient boosting model was trained on TCGA-HNSC based on GSVA scores for gene sets of the MSigDB to robustly predict UCHL1-related cancers in other solid tumors and cancer cell lines derived thereof. Potential vulnerabilities of UCHL1-related cancer cells were elucidated by an in-silico drug screening approach. Results: We established a 497-gene set, which stratified the TCGA-HNSC cohort into distinct subgroups with a UCHL1-related or other phenotype. UCHL1-related HNSC were characterized by higher frequencies of genomic alterations, which was also evident for UCHL1-related cancers of other solid tumors predicted by the classification model. These data indicated an impaired maintenance of genomic integrity and vulnerability for DNA-damaging treatment, which was supported by a favorable prognosis of UCHL1-related tumors after radiotherapy, and a higher sensitivity of UCHL1-related cancer cells to irradiation or DNA-damaging compounds (e.g., Oxaliplatin). Conclusion: Our study established UCHL1-related cancers as a novel subgroup across most solid tumor entities with a unique molecular phenotype and DNA-damaging treatment as a specific vulnerability, which requires further proof-of-concept in pre-clinical models and future clinical trials.