HPV’s hidden hand – RNA sequencing and spatial transcriptomics reveal key details in head and neck cancer microenvironment

Google+ Pinterest LinkedIn Tumblr +


Using RNA sequencing and spatial transcriptomics, the roles of PKM gene and EPHA2 pathway in the development of HNSCC were discovered.

Head and neck squamous cell carcinoma (HNSCC) is a type of cancer, and its prevalence can be associated with oncogenic human papillomavirus (HPV) strains

Head and neck squamous cell carcinoma (HNSCC) is a type of cancer that affects the mucous membranes of the mouth, nose, and throat. HNSCC is typically associated with tobacco exposure, alcohol abuse, and viral infections. The links between human papillomavirus (HPV) infection status and the molecular characteristics of HNSCC are not clearly defined.

In a novel and pioneering research work, a team of scientists led by Yun Hak Kim, Assistant Professor in the Department of Anatomy and Department of Biomedical Informatics at the School of Medicine at Pusan National University, Korea, has recently studied the cellular diversity and molecular attributes within the complex tumor microenvironment of HNSCC with respect to the status of HPV infection. Their research findings were recently published in the Journal of Medical Virology on 18 January 2024.

“The incidence of HNSCC has gone up considerably in the recent years and yet very little is known about the exact molecular mechanisms regarding the role of HPV infection and its link with HNSCC”, shares Dr. Kim, revealing the inspiration behind the present research work.

Using samples from both HPV negative (HPV) and HPV positive (HPV+) patients for single-cell RNA profiling, a genetic sequencing technique, and spatial transcriptomics (ST), an advanced molecular profiling technique, the primary tumors and lymph node metastatic tumors (LNMTs) were analyzed. The in-depth research findings provide remarkable imagery into the complex links between HPV infection and the development of HNSCC.

Single-cell transcriptional profiling of HNSCC patients

(A) Schematic overview of the overall study design. Eight samples were collected from the primary sites and lymph nodes of patients with HNSCC for scRNA-seq. (B) Copy number of HNSCC samples from the scRNA expression data. Epithelial cells were used to identify tumor cells, whereas immune cells were used as reference cells. Red and blue indicate chromosomal amplifications and deletions, respectively. (C) UMAP plot of 29,987 cells from 4 patients; each of the colors represents the 17 cell types and their transcript counts. (D) Cell proportion of 17 cell types for each tissue origin. (E) Cell proportion of 17 cell types according to HPV infection status. (F) Cell distribution in individual tissues of patients with HNSCC. PT1 and PT2 tested negative for HPV, whereas PT3 and PT4 tested positive for HPV. (G) UMAP plot of monocyte/macrophage subclusters. (H) Bar plot of monocyte/macrophage subclusters. The percentage of monocyte/macrophage subclusters at each site separated by HPV infection status is shown in the left panel, and the number of cells is shown in the right panel.

In the case of HPV+ HNSCC, immune cells actively participated in the development of tumors. The molecular characteristics of HNSCC were highly influenced by the status of HPV infection rather than the actual location of tumor. The analyses of RNA sequencing and ST indicated the crucial role of pyruvate kinase muscle (PKM) gene in the development of cancer stem cell-like populations in the LNMT of HPV tissue. In HPV HNSCC cell lines, the expression of PKM gene was consistent and indicated a role in the spheroid structure formation ability, which was studied using genetic knockdown of PKM2 gene. Interestingly, during the ST analysis of the HPV+ patients, the presence of such ectopic lymphoid structure was confirmed.

Dr. Kim and their team went a step further and focused their efforts to study the critical angiogenesis-new blood vessel formation pathway. They used single-cell RNA and ST profiles of HPV patients and found that Ephrin-A (EPHA2) was mainly involved in angiogenesis and cell migration pathways.

“Our novel research findings of PKM and EPHA2 molecular pathways can be utilized as therapeutic targets in suppressing the growth and spread of HNSCC cancer cells. Also, the tumor microenvironment is affected by the status of HPV infection and thus, precision medicine can be developed against HNSCC” concludes Dr. Kim, underlining the potential real-life applications of the research work.

Overall, this study significantly enhances our understanding of HNSCC microenvironment. Thus, it can guide new approaches to HNSCC treatment, especially in inhibiting cancer cell growth and migration. Therefore, the HPV infection could be considered a criterion for developing treatment strategies, thereby enabling precision medicine tailored to the patient’s condition!


Lee H, Park S, Yun JH, Seo C, Ahn JM, Cha HY, Shin YS, Park HR, Lee D, Roh J, Heo HJ, Baek SE, Kim EK, Lee HS, Kim CH, Kim YH, Jang JY. (2024) Deciphering head and neck cancer microenvironment: Single-cell and spatial transcriptomics reveals human papillomavirus-associated differences. J Med Virol 96(1):e29386. [article]
Share.