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.

Identification of the HNSC88 Molecular Signature for Predicting Subtypes of Head and Neck Cancer

Development and Validation of Radiomic Signatures of Head and Neck Squamous Cell Carcinoma Molecular Features and Subtypes.

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

Head and neck squamous cell carcinoma-specific prognostic signature and drug sensitive subtypes based on programmed cell death-related genes

Characterization of Molecular Subtypes in Head and Neck Squamous Cell Carcinoma With Distinct Prognosis and Treatment Responsiveness

Identification of Molecular Subtypes in Head and Neck Squamous Cell Carcinoma Based on Dysregulated Immune LncRNAs.

Molecular Subtypes Based on Cell Differentiation Trajectories in Head and Neck Squamous Cell Carcinoma: Differential Prognosis and Immunotherapeutic Responses

Validation of Two Biologically Distinct Hpv-Associated Head and Neck Cancer Subtypes and Correlation with E5 Expression

Identification of Immune Subtypes for Predicting the Prognosis of Patients in Head and Neck Squamous Cell Carcinoma

Systematic Analysis of Molecular Subtypes and Immune Prediction Based on CD8 T Cell Pattern Genes Based on Head and Neck Cancer

Subtyping of Human Papillomavirus-Positive Cervical Cancers Based on the Expression Profiles of 50 Genes

Comprehensive analysis reveals novel gene signature in head and neck squamous cell carcinoma: predicting is associated with poor prognosis in patients

Subclass-Specific Prognosis and Treatment Efficacy Inference in Head and Neck Squamous Carcinoma.

Integrative Analysis of Head and Neck Cancer Identifies Two Biologically Distinct HPV and Three Non-Hpv Subtypes.

Distinguishable Prognostic miRNA Signatures of Head and Neck Squamous Cell Cancer With or Without HPV Infection

Identification of molecular classification and gene signature for predicting prognosis and immunotherapy response in HNSCC using cell differentiation trajectories

Integrative Single-Cell and Bulk Transcriptomes Analyses Identify Intrinsic HNSCC Subtypes with Distinct Prognoses and Therapeutic Vulnerabilities

Identification of Potential Biomarkers and Survival Analysis for Head and Neck Squamous Cell Carcinoma Using Bioinformatics Strategy: A Study Based on TCGA and GEO Datasets

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

Identification of a Seven-Gene Signature and Establishment of a Nomogram Predicting Overall Survival in Head and Neck Squamous Cell Carcinoma

Identification of an 8 HPV-related RNA signature as a novel prognostic biomarker for squamous cell carcinoma of the head and neck