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.

Aspartate-β-hydroxylase and hypoxia marker expression in head and neck carcinomas: implications for HPV-associated tumors

Diverse Molecular Functions Of Aspartate Beta-Hydroxylase In Cancer

Noncanonical HPV carcinogenesis drives radiosensitization of head and neck tumors

Expression and prognostic significance of a panel of tissue hypoxia markers in head-and-neck squamous cell carcinomas

Clinical relevance and therapeutic predictive ability of hypoxia biomarkers in head and neck cancer tumour models

Aspartate-β-hydroxylase as a promising prognostic and therapeutic biomarker for HBV-positive hepatocellular carcinoma

Association between Tumor Microbiome and Hypoxia across Anatomic Subsites of Head and Neck Cancers

Tumor microenvironment in head and neck squamous cell carcinomas: Predictive value and clinical relevance of hypoxic markers. A review

Aspartate β-hydroxylase Regulates Expression of Ly6 Genes

A novel defined hypoxia-related gene signature to predict the prognosis of oral squamous cell carcinoma

The Interplay of HIF1A, DNA Methylation, and Prolyl Hydroxylases in Head and Neck Squamous Cell Carcinoma

Construction of a hypoxia-derived gene model to predict the prognosis and therapeutic response of head and neck squamous cell carcinoma

Clinical implications of activation of the LIMD1-VHL-HIF1α pathway during head-&-neck squamous cell carcinoma development

Aspartate β-hydroxylase as a target for cancer therapy

Bioinformatic Analysis Identifies HPV-related Tumor Microenvironment Remodeling Prognostic Biomarkers in Head and Neck Squamous Cell Carcinoma.

Hydroxylase Activity of ASPH Promotes Hepatocellular Carcinoma Metastasis Through Epithelial-to-Mesenchymal Transition Pathway.

Hypoxia-derived molecular subtype and gene signature characterize prognoses and therapeutic responses in head and neck squamous cell carcinoma

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

Diverse Molecular Functions of Aspartate Β-Hydroxylase in Cancer

Oxygen Deprivation Modulates EGFR and PD-L1 in Squamous Cell Carcinomas of the Head and Neck

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