ABSTRACT
Head and neck squamous cell carcinoma (HNSCC) is the sixth most common malignancy worldwide, and the development of novel therapeutic strategies for HNSCC requires a profound understanding of tumor cells and the tumor microenvironment (TME). Additionally, HNSCC has a poor prognosis, necessitating the use of genetic markers for predicting clinical outcomes in HNSCC. In this study, we performed single-cell sequencing analysis on tumor tissues from seven HNSCC patients, along with one adjacent normal tissue. Firstly, the analysis of epithelial cell clusters revealed two clusters of malignant epithelial cells, characterized by unique gene expression patterns and dysregulated signaling pathways compared to normal epithelial cells. Secondly, the examination of the TME unveiled extensive crosstalk between fibroblasts and malignant epithelial cells, potentially mediated through ligand-receptor interactions such as COL1A1-SDC1, COL1A1-CD44, and COL1A2-SDC1. Furthermore, transcriptional heterogeneity was observed in immune cells present in the TME, including macrophages and dendritic cells. Finally, leveraging the gene expression profiles of malignant epithelial cells, we developed a prognostic model comprising six genes, which we validated using two independent datasets. These findings shed light on the heterogeneity within HNSCC tumors and the intricate interplay between malignant cells and the TME. Importantly, the developed prognostic model demonstrates high efficacy in predicting the survival outcomes of HNSCC patients.
