Head and neck squamous cell carcinomas (HNSCC) are still one of the most common malignant tumors in China, with a high metastasis rate and poor prognosis. The tumor immune microenvironment can affect the occurrence, development and prognosis of tumors, but the underlying mechanism is still unclear. In this study, we tried to describe the correlation between the recurrence of HNSCC and the tumor microenvironment (TME). The expression data [estimate the level of tumor stromal and immune infiltration, expression data (ESTIMATE)] algorithm was used to identify and estimate highly correlated stromal cells, immune cells, and prognostic scores in 116 samples of head and neck cancer patients from The Cancer Genome Atlas (TCGA) dataset. The functional enrichment analysis and protein-protein interaction (PPI) networks of differential expressed genes (DEGs) were constructed. Subsequently, the abundance of various infiltrating immune cells was estimated with the tumor immune estimation resource (TIMER) and the infiltration pattern of immune cells were explored in HNSCC. A total of 407 immune-related genes were identified to involve in the TME. We found that CCR5, CD3E, CD4, and HLA -DRB1 were the most obvious DEGs and the dendritic cells (DCs) showed the highest abundance in the TME of HNSCC. In addition, the unsupervised cluster analysis determined 10 clusters of immune infiltration patterns, and indicated that immune infiltrated CD4 + T and B cells may be related to the prognosis of HNSCC. In conclusion, our research determined the list of immune genes and immune infiltrating cells related to the prognosis of HNSCC, and provided a perspective for HNSCC evolution, anti-tumor drugs selection, and drug resistance research.