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.

Multimodal assessment of high-risk human papillomavirus in sinonasal squamous cell carcinoma

The presence of high-risk human papillomavirus (HPV) E6/E7 mRNA transcripts in a subset of sinonasal carcinomas is evidence of involvement of HPV in its etiopathogenesis

Sinonasal squamous cell carcinoma in the United States: Temporal and geographic patterns associated with HPV testing and positivity

Prevalence of Human Papillomavirus in Sinonasal Squamous Cell Carcinoma with and Without Association of Inverted Papilloma in Eastern China

Molecular patterns and mechanisms of tumorigenesis in HPV-associated and HPV-independent sinonasal squamous cell carcinoma

Human papillomavirus-related carcinomas of the sinonasal tract

Diagnostic utility of high‐risk human papillomavirus mRNA in situ hybridisation in squamous cell carcinoma of unknown primary in the head and neck and implementing American Society of Clinical Oncology guideline recommendations

Noncanonical HPV carcinogenesis drives radiosensitization of head and neck tumors

HPV Positive Status Is a Favorable Prognostic Factor in Non-Nasopharyngeal Head and Neck Squamous Cell Carcinoma Patients: A Retrospective Study From the Surveillance, Epidemiology, and End Results Database

Low Prevalence of Human Papillomavirus Infection in Sinonasal Inverted Papilloma and Oncocytic Papilloma

Exploring the roles of HPV16 variants in head and neck squamous cell carcinoma: current challenges and opportunities

Prognostic factors in sinonasal cancers: A multicenter pooled analysis.

Uncommon and Challenging Phenotypes of High-Risk Human Papillomavirus-Associated Head and Neck Carcinomas Revealed by High-Throughput Studies

Prevalence of genotype-specific human papillomavirus in cytology specimens and cervical biopsies, and its implication in cervical cancer risk stratification: a retrospective study of 10647 cases

A Proportion of Primary Squamous Cell Carcinomas of the Parotid Gland Harbour High-Risk Human Papillomavirus.

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

Neutrophils infiltration and its correlation with human papillomavirus status in the oral squamous cell carcinoma.

Differences in the Prevalence of Human Papillomavirus (HPV) in Head and Neck Squamous Cell Cancers by Sex, Race, Anatomic Tumor Site, and HPV Detection Method

Survival Outcomes and Prognostic Factors of Squamous Cell Carcinomas Arising from Sinonasal Inverted Papillomas: a Retrospective Analysis of 120 Patients

Human Papillomavirus in Metastatic Lymph Nodes from Unknown Primary Head and Neck Squamous Cell Carcinoma

Correlation of Immediate Prevalence of Cervical Precancers and Cancers with HPV Genotype and Age in Women with ASC-US/hrHPV+: a Retrospective Analysis of 2292 Cases