A comprehensive prognostic score for head and neck squamous cancer driver genes and phenotype traits

Abstract

Background Head and neck squamous cancer (HNSCC) presents variable phenotype and progression features. Clinically applicable, high-accuracy multifactorial prognostic models for HNSCC survival outcomes are warranted and an active area of research. This study aimed to construct a comprehensive prognostic tool for HNSCC overall survival by integrating cancer driver genes with tumor clinical and phenotype information. Methods Key overall survival-related cancer driver genes were screened from among main effector and reciprocal gene pairs using TCGA data using univariate Cox proportional hazard regression analysis. Independent validation was performed using the GSE41613 dataset. The main effector genes among these were selected using LASSO regression and transcriptome score modeling was performed using multivariate Cox regression followed by validation analysis of the prognostic score. Next, multivariate Cox regression analysis was performed using the transcriptome score combined with age, grade, gender, and stage. An ‘Accurate Prediction Model of HNSCC Overall Survival Score’ (APMHO) was computed and validated. Enriched functional pathways, gene mutational landscape, immune cell infiltration, and immunotherapy sensitivity markers associated with high and low APMHO scores were analyzed. Results Screening 107 overall survival-related cancer genes and 402 interacting gene pairs, 6 genes: CRLF2, HSP90AA1, MAP2K1, PAFAH1B2, MYCL and SET genes, were identified and a transcriptional score was obtained. Age, stage and transcriptional score were found to be significant predictors in Cox regression analysis and used to construct a final APMHO model showing an AUC > 0.65 and validated. Transcriptional score, age, pathologic_N, pathologic_T, stage, and TCGA_subtype were significantly different in distribution between high and low APMHO groups. High APMHO samples showed significantly higher mutation rate, enriched tumor-related pathways including Hypoxia, unfold_protein_response, Glycolysis, and mTORC1 signaling, along with differences in immune cell infiltration and immune checkpoint, interferon-γ pathway and m6A regulator expression patterns. Conclusion The APMHO score combining transcriptional and clinical variables showed good prognostic ability for HNSCC overall survival outcomes and was associated with different patterns of phenotypical features, immune and mutational landscape, and immunotherapy sensitivity marker expression. Future studies should validate this score in independent clinical cohorts.