Furthermore, co-expression analysis in TCGA HPV − HNSCC cohort revealed that the correlation between IL6 and CCR2 exhibited the lowest value (R = 0.23) among the 5 DEGs identified from our scRNA-seq data analysis (Fig. 3E and F). This evidence concerns the gene IL6 and head and neck squamous cell carcinoma.