Abstract:Background: Head and neck squamous cell carcinomas (HNSCCs) are derived from the mucosal linings of the upper aerodigestive tract, salivary glands, thyroid, oropharynx, larynx, and hypopharynx. The present study aimed to identify the novel genes and pathways underlying HNSCC. Despite the advances in HNSCC research, diagnosis, and treatment, its incidence continues to rise, and the mortality of advanced HNSCC is expected to increase by 50%. Therefore, there is an urgent need for effective biomarkers to predict HNSCC patients' prognosis and provide guidance to the personalized treatment. Methods: Both HNSCC clinical and gene expression data were abstracted from The Cancer Genome Atlas (TCGA) database. Intersecting analysis was adopted between the gene expression matrix of HNSCC patients from TCGA database to extract TME-related genes. Differential gene expression analysis between HNSCC tissue samples and normal tissue samples was performed by R software. Then, HNSCC patients were categorized into clusters 1 and 2 via NMF. Next, TME-related prognosis genes (p < 0.05) were analyzed by univariate Cox regression analysis, LASSO Cox regression analysis, and multivariate Cox regression analysis. Finally, nine genes were selected to construct a prognostic risk model and a prognostic gene signature. We also established a nomogram using relevant clinical parameters and a risk score. The Kaplan-Meier curve, survival analysis, time-dependent receiver operating characteristic (ROC) analysis, decision curve analysis (DCA), and the concordance index (C-index) were carried out to assess the accuracy of the prognostic risk model and nomogram. Potential molecular mechanisms were revealed by gene set enrichment analysis (GSEA). Additionally, gene correlation analysis and immune cell correlation analysis were conducted for further enriching our results. Results: A novel HNSCC prognostic model was established based on the nine genes (GTSE1, LRRN4CL, CRYAB, SHOX2, ASNS, KRT23, ANGPT2, HOXA9, and CARD11). The value of area under the ROC curves (AUCs) (0.769, 0.841, and 0.816) in TCGA whole set showed that the model effectively predicted the 1-, 3-, and 5-year overall survival (OS). Results of the Cox regression assessment confirmed the nine-gene signature as a reliable independent prognostic factor in HNSCC patients. The prognostic nomogram developed using multivariate Cox regression analysis showed a superior C-index over other clinical signatures. Also, the calibration curve had a high level of concordance between estimated OS and the observed OS. This showed that its clinical net can precisely estimate the one-, three-, and five-year OS in HNSCC patients. The gene set enrichment analysis (GSEA) to some extent revealed the immune- and tumor-linked cascades. Conclusion: In conclusion, the TME-related nine-gene signature and nomogram can effectively improve the estimation of prognosis in patients with HNSCC.

Identification of HOXB9 based on comprehensive bioinformatics analysis for predicting prognosis of head and neck squamous cell carcinoma

Pan-cancer Analysis of Homeobox B9 As a Predictor for Prognosis and Immunotherapy in Human Tumors.

Identification of HOXB9 to predict prognosis of endometrial cancer based on comprehensive bioinformatics analysis

Comprehensive analysis of a homeobox family gene signature in clear cell renal cell carcinoma with regard to prognosis and immune significance

HOXD9 is a Potential Prognostic Biomarker Involved in Immune Microenvironment of Glioma.

CD79A Work as a Potential Target For The Prognosis of Patients With HNSCC: Analysis of Immune Cell Infiltration In Head and Neck Squamous Cell Carcinoma Based on The CIBERSORTx Deconvolution Algorithm

The Effects of Differentially-Expressed Homeobox Family Genes on the Prognosis and HOXC6 on Immune Microenvironment Orchestration in Colorectal Cancer

Nine-gene signature and nomogram for predicting survival in patients with head and neck squamous cell carcinoma

Exploring the prognosis value, immune correlation, and drug responsiveness prediction of homeobox C6 (HOXC6) in lung adenocarcinoma

Data mining on identifying diagnosis and prognosis biomarkers in head and neck squamous carcinoma

Analysis of the Prognosis and Therapeutic Value of the CXC Chemokine Family in Head and Neck Squamous Cell Carcinoma

Multi-omic analysis and validation reveal ZBP1 as a potential prognostic and immunotherapy-related biomarker in head and neck squamous cell carcinoma

Bioinformatic Analysis Identifies HPV-related Tumor Microenvironment Remodeling Prognostic Biomarkers in Head and Neck Squamous Cell Carcinoma.

Pan-cancer exploration of oncogenic and clinical impacts revealed that HOXA9 is a diagnostic indicator of tumorigenesis

Diagnostic and prognostic value of FOXD1 expression in head and neck squamous cell carcinoma

Identification of HOXC Gene Family as Prognostic and Immune-Related Biomarkers in Breast Cancer Through mRNA Transcriptional Profile and Experimental Validation

Identification of Key Genes Associated with Head and Neck Squamous Cell Carcinoma in the Microenvironment

Mining of gene modules and identification of key genes in head and neck squamous cell carcinoma based on gene co-expression network analysis.

Prognostic value of novel immune-related genomic biomarkers identified in head and neck squamous cell carcinoma

Identification of M6a Methyltransferase-Related Genes Predicts Prognosis and Immune Infiltrates in Head and Neck Squamous Cell Carcinoma

Overexpression of HOXB9 Promotes Metastasis and Indicates Poor Prognosis in Colon Cancer.