Fig. 2

scRNA-seq data analysis identifies DEGs and enriched pathways within NK cells and tumor cells stratified by HPV infection status. Bioinformatic analysis of tumor cells based on the scRNA-seq data from GSE164690. A Louvain clustering on all analyzed cells identified 10 cell clusters. B Identification of cell types. C Visualization of identified tumor cells using UMAP, colored by sample ID. D Visualization of identified tumor cells using UMAP, colored by HPV status. E Enriched pathways in HPV − HNSCC cells vs HPV + HNSCC cells. F Five specific gene pairs identified between NK cells and HPV − HNSCC cells. G Visualization of the top 10 DEGs in (HPV − vs HPV +) in epithelial tumor cells