Myeloid cells are a major heterogeneous cell population in the tumor immune microenvironment (TIME). Imbalance of myeloid response remains a major obstacle to a favorable prognosis and successful immune therapy. Therefore, we aimed to construct a risk model to evaluate the myeloid contexture, which may facilitate the prediction of prognosis and immune infiltration in patients with head and neck squamous cell carcinoma (HNSCC). In our study, six myeloid signature genes (including CCL13, CCR7, CD276, IL1B, LYVE1 and VEGFC) analyzed from 52 differentially expressed myeloid signature genes were finally pooled to establish a prognostic risk model, termed as myeloid gene score (MGS) in a training cohort and validated in a test cohort and an independent external cohort. Furthermore, based on the MGS subgroups, we were able to effectively identify patients with a poor prognosis, aggressive clinical parameters, immune cell infiltration status and immunotherapy response. Thus, MGS may serve as an effective prognostic signature and predictive indicator for immunotherapy response in patients with HNSCC.