In the present study, by integrating scRNA-seq data with TCGA HNSCC dataset, we identified five enriched pathways, including IL6/IL6R and CCL2/CCR2 signaling, within NK cells and tumor cells in the HPV − HNSCC cohort compared to the HPV + counterpart. The gene discussed is IL6; the disease is neoplasm.