Abstract:Background. Head and neck squamous cell carcinoma (HNSCC) is a significant health problem and related to poor long-term outcomes, indicating more research to be done to deeply understand the underlying pathways. Objective. This current study aimed in the assessment of the viral- (especially human papilloma virus [HPV]) and carcinogen-driven head and neck squamous cell carcinoma (HNSCC) microenvironment based on single-cell sequencing analysis. Methods. Data were downloaded from GEO database (GSE139324), including 131224 cells from 18 HP- HNSCC patients and 8 HPV+ HNSCC patients. Following data normalization, all highly variable genes in single cells were identified, and batch correction was applied. Differentially expressed genes were identified using Wilcoxon rank sum test. A gene enrichment analysis was performed in each cell cluster using KEGG analysis. Single-cell pseudotime trajectories were constructed with MONOCLE (version 2.6.4). Cell-cell interactions were analyzed with CellChat R package. Additionally, cell-cell communication patterns in key signal pathways were compared in different tissue groups. A hidden Markov model (HMM) was used to predict gene expression states (on or off) throughout pseudotime. Five-year overall survival outcomes were compared in both HPV+ and HPV- subsets. Results. 20,978 high-quality individual cells passed quality control. RNA-seq data were used from 522 HNSCC primary tumor samples. 1,137 differentially expressed genes between HPV+ and HPV- HNSCC patients were investigated. 96 differentially expressed genes were associated with overall survival and highly enriched in B cell associated biological process. Cell composition differed between types of samples. MHC-I, MHC-II, and MIF signaling pathways were found to be most relevant. Within these pathways, some cells were either signal receiver or signal sender, depending on sample type, respectively. Six genes were obtained, AREG and TGFBI (upregulation), CD27, CXCR3, MS4A1, and CD19 (downregulation), whose expression and HPV types were highly associated with worse overall survival. AREG and TGFBI were pDC marker genes, CXCR3 and CD27 were significantly expressed in T cell-related cells, while MS4A1 and CD19 were mainly expressed in B naïve cells. Conclusions. This study revealed dynamic changes in cell percentage and heterogeneity of cell subtypes of HNSCC. AREG, TGFBI, CD27, CXCR3, MS4A1, and CD19 were associated with worse overall survival in HPV-related HNSCC. Especially B-cell related pathways were revealed as particularly relevant in HPV-related HNSCC. These findings are a basis for the development of biomarkers and therapeutic targets in respective patients.

The tumor immune microenvironments of HPV+ and HPV- head and neck cancers

The immune microenvironment of HPV-positive and HPV-negative oropharyngeal squamous cell carcinoma: a multiparametric quantitative and spatial analysis unveils a rationale to target treatment-naïve tumors with immune checkpoint inhibitors

Noncanonical HPV carcinogenesis drives radiosensitization of head and neck tumors

Single-Cell Transcriptome Analysis Reveals Different Immune Signatures in HPV- and HPV + Driven Human Head and Neck Squamous Cell Carcinoma

Immune microenvironment of human papillomavirus-positive head and neck squamous cell carcinoma

Mutational profiles of head and neck squamous cell carcinomas based upon human papillomavirus status in the Veterans Affairs National Precision Oncology Program

Comparative analysis of the tumor microbiome, molecular profiles, and immune cell abundances by HPV status in mucosal head and neck cancers and their impact on survival

Deciphering the Tumor–Immune–Microbe Interactions in HPV-Negative Head and Neck Cancer

Heterogeneity of the Head and Neck Squamous Cell Carcinoma Immune Landscape and Its Impact on Immunotherapy

The Microenvironment of Head and Neck Cancers: Papillomavirus Involvement and Potential Impact of Immunomodulatory Treatments

Viral transcript and tumor immune microenvironment-based transcriptomic profiling of HPV-associated head and neck squamous cell carcinoma identifies subtypes associated with prognosis

Tumor microenvironment: an evil nexus promoting aggressive head and neck squamous cell carcinoma and avenue for targeted therapy

Human Papillomavirus Infection in Head and Neck Squamous Cell Carcinomas: Transcriptional Triggers and Changed Disease Patterns

Molecular Tumor Subtypes of HPV-Positive Head and Neck Cancers: Biological Characteristics and Implications for Clinical Outcomes

Immunotherapy Advances in Locally Advanced and Recurrent/Metastatic Head and Neck Squamous Cell Carcinoma and Its Relationship With Human Papillomavirus

Deciphering head and neck cancer microenvironment: Single‐cell and spatial transcriptomics reveals human papillomavirus‐associated differences

Deciphering the role of HPV-mediated metabolic regulation in shaping the tumor microenvironment and its implications for immunotherapy in HNSCC

HPV-positive status associated with inflamed immune microenvironment and improved response to anti-PD-1 therapy in head and neck squamous cell carcinoma

HPV infection related immune infiltration gene associated therapeutic strategy and clinical outcome in HNSCC

HPV Integration in HNSCC Correlates with Survival Outcomes, Immune Response Signatures, and Candidate Drivers

Immune evasion in HPV − head and neck precancer–cancer transition is driven by an aneuploid switch involving chromosome 9p loss