UC Davis

BACKGROUND: Human papillomavirus (HPV)-positive head and neck squamous cell carcinoma (HNSCC) is an increasingly common malignancy. We aimed to explore the immune heterogeneity of natural killer (NK) cells in HPV-positive HNSCC. METHODS: Single-cell RNA-sequencing (scRNA-seq) and bulk RNA-sequencing datasets of HPV-positive HNSCC data were obtained from the Gene Expression Omnibus (GEO) database. Seurat, harmony, and SingleR were used to perform the scRNA-seq analysis. Subsequently, the cellphonedb package was used for the cell crosstalk analysis, and the clusterProfiler package was used for the hallmark pathway enrichment analysis. Finally, the gene set variation analysis (GSVA) package was used for the immune cell infiltration, Tumor Immune Dysfunction and Exclusion (TIDE), and risk-score analyses. RESULTS: A total of 30,562 cells were classified into 9 cell clusters that comprised 6 main cell types [i.e., T cells, natural killer T (NKT) cells, NK cells, B cells, plasma cells, and macrophages]. The NK cells were then further clustered into 3 tissue-resident NK (trNK0-2) and 2 tumor-associated NK (taNK0-1) cell types. The trNK0 cell type, which exhibited inhibitory cancer hallmark activity, appeared to exert potential anti-tumor effects via trNK0-macrophage crosstalk. The trNK score could serve as an independent and valuable prognostic classifier, as the patients with high-trNK scores had better outcomes, immune-infiltration levels, and immunotherapy effects. CONCLUSIONS: Using an scRNA-seq analysis, we identified a specific type of tissue-resident NK cell (i.e., trNK-0) that was involved in the anti-tumor and immunotherapy response in HPV-positive HNSCC.