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.

Fewer tumour-specific PD-1+CD8+ TILs in high-risk “Infiltrating” HPV− HNSCC

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

Identification of lncRNA Signature of Tumor-Infiltrating T Lymphocytes With Potential Implications for Prognosis and Chemotherapy of Head and Neck Squamous Cell Carcinoma

Bioinformatics Analysis for the Biomarkers of the Tumor-Infiltrating Lymphocytes in Head and Neck Squamous Cell Carcinoma.

Noncanonical HPV carcinogenesis drives radiosensitization of head and neck tumors

The Presence of Tumour-Infiltrating Lymphocytes (tils) and the Ratios Between Different Subsets Serve As Prognostic Factors in Advanced Hypopharyngeal Squamous Cell Carcinoma.

PD-L1 in Combination with CD8+TIL and HIF-1α are Promising Prognosis Predictors of Head and Neck Squamous Cell Carcinoma

Analysis of Immunological Characteristics and Genomic Alterations in HPV-Positive Oropharyngeal Squamous Cell Carcinoma Based on PD-L1 Expression

Correlation of immune makers with HPV 16 infections and the prognosis in oropharyngeal squamous cell carcinoma

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

Functional HPV-specific PD-1+ stem-like CD8 T cells in head and neck cancer

HPV16 Status Might Correlate to Increasing Tumor-Infiltrating Lymphocytes in Hypopharyngeal Cancer.

PD-L1 expression and immune stromal features in HPV-independent cervical adenocarcinoma

The role of tumor infiltrating lymphocytes (TILs) and the ratios between different subsets serve as prognostic factors in hypopharyngeal squamous cell carcinoma

Blocking of EGFR Signaling Is a Latent Strategy for the Improvement of Prognosis of HPV-Induced Cancer

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

Tertiary Lymphoid Structure Raises Survival and Immunotherapy in HPV− HNSCC

Massive PD-L1 and CD8 double positive TILs characterize an immunosuppressive microenvironment with high mutational burden in lung cancer

Prognostic Value of Tumor PD-L1 Expression Combined with CD8+ Tumor Infiltrating Lymphocytes in High Grade Serous Ovarian Cancer

Deciphering the interplay of HPV infection, MHC-II expression, and CXCL13+ CD4+ T cell activation in oropharyngeal cancer: implications for immunotherapy

754 TIGIT-PVR is a Key Immune Checkpoint and Therapeutic Target in HPV-positive Head and Neck Squamous Cell Carcinomas