Abstract The tumor microenvironment (TME) is of great clinical significance for predicting the therapeutic effect of tumors. Nonetheless, there was no systematic analysis of cellular interactions in the TME of head and neck cancer (HNSC). This study used gene expression data from 816 patients with HNSC to analyze the scores of 22 immune cells. On this basis, we have established a novel TMEscore-based prognostic risk model. The relationship between TMEscore and clinical and genomic characteristics was analyzed. The sample was divided into risk-H and risk-L groups based on the prognosis risk model of TMEscore, with significant differences in overall survival between the two groups (log rank p