Fig. 1

NK cell enrichment in HPV − HNSCC by scRNA-seq analysis. A Schematic showing analysis from bioinformatics to experimental validation. (B-E) Bioinformatic analysis of immune cells based on the scRNA-seq data from GSE164690. B Visualization of all analyzed cells stratified by HPV infection status using UMAP, colored by sample ID. C Visualization of all analyzed cells using UMAP, colored by HPV status. D Louvain clustering on all analyzed cells identified 10 cell clusters. E Correlation heatmap of all identified cell types. F-I Bioinformatic analysis of immune cells based on the scRNA-seq data from GSE139324. F Visualization of all analyzed cells stratified by HPV infection status using UMAP, colored by sample ID. G Visualization of all analyzed cells using UMAP, colored by HPV status. H Louvain clustering on all analyzed cells identified 9 cell clusters. I Correlation heatmap of all identified cell types. J Enriched pathways derived from tumor infiltrating NK cells (HPV − tumors vs. HPV + tumors) using the scRNA-seq data from GSE164690 (upper panel) and GSE139324 (lower panel)