Genomic heterogeneity and copy number variant burden are associated with poor recurrence-free survival and 11q loss in human papillomavirus-positive squamous cell carcinoma of the oropharynx

Metabolic pathway-based subtypes associate glycan biosynthesis and treatment response in head and neck cancer

Head and Neck Squamous Cell Carcinoma (HNSCC) is a heterogeneous malignancy that remains a significant challenge in clinical management due to frequent treatment failures and pronounced therapy resistance. While metabolic dysregulation appears to be a critical factor in this scenario, comprehensive analyses of the metabolic HNSCC landscape and its impact on clinical outcomes are lacking. This study utilized transcriptomic data from four independent clinical cohorts to investigate metabolic heterogeneity in HNSCC and define metabolic pathway-based subtypes (MPS). In HPV-negative HNSCCs, MPS1 and MPS2 were identified, while MPS3 was enriched in HPV-positive cases. MPS classification was associated with clinical outcome post adjuvant radio(chemo)therapy, with MPS1 consistently exhibiting the highest risk of therapeutic failure. MPS1 was uniquely characterized by upregulation of glycan (particularly chondroitin/dermatan sulfate) metabolism genes. Immunohistochemistry and pilot mass spectrometry imaging analyses confirmed this at metabolite level. The histological context and single-cell RNA sequencing data identified the malignant cells as key contributors. Globally, MPS1 was distinguished by a unique transcriptomic landscape associated with increased disease aggressiveness, featuring motifs related to epithelial-mesenchymal transition, immune signaling, cancer stemness, tumor microenvironment assembly, and oncogenic signaling. This translated into a distinct histological appearance marked by extensive extracellular matrix remodeling, abundant spindle-shaped cancer-associated fibroblasts, and intimately intertwined populations of malignant and stromal cells. Proof-of-concept data from orthotopic xenotransplants replicated the MPS phenotypes on the histological and transcriptome levels. In summary, this study introduces a metabolic pathway-based classification of HNSCC, pinpointing glycan metabolism-enriched MPS1 as the most challenging subgroup that necessitates alternative therapeutic strategies.

Unravelling the heterogeneity of oral squamous cell carcinoma by integrative analysis of single-cell and bulk transcriptome data

Oral squamous cell carcinoma (OSCC) is a prevalent malignancy of the head and neck with rising global incidence. Despite advances in treatment modalities, OSCC prognosis remains diverse due to the complex molecular and cellular heterogeneity within tumours, as well as the heterogeneity in tumour microenvironment (TME). In this study, we utilized single-cell RNA sequencing (scRNA-seq) analysis to explore distinct subpopulations of tumour cells in OSCC tissues and their interaction with components in TME. We identified four major tumour cell subpopulations (C0, C1, C2 and C3) with unique molecular characteristics and functional features. Pathway enrichment analysis revealed that C0 primarily expressed genes involved in extracellular matrix interactions and C1 showed higher proliferation levels, suggesting that the two cell subpopulations exhibited tumour aggressiveness. Conversely, C2 and C3 displayed features associated with keratinization and cornified envelope formation. Accordingly, C0 and C1 subpopulations were associated with shorter overall and disease-free survival times, while C2 and C3 were weakly correlated with longer survival. Genomic analysis showed that C1 demonstrated a positive correlation with tumour mutation burden. Furthermore, C0 exhibited resistant to cisplatin treatment, while C1 showed more sensitive to cisplatin treatment, indicating that C0 might exhibit more aggressive compared to C1. Additionally, C0 had a higher level of communication with fibroblasts and endothelial cells in TME via integrin-MAPK signalling, suggesting that the function of C0 was maintained by that pathway. In summary, this study provided critical insights into the molecular and cellular heterogeneity of OSCC, with potential implications for prognosis prediction and personalized therapeutic approaches.

Peripheral surrogates of tumor burden to guide chemotherapeutic and immunotherapeutic strategies for HPV-associated malignancies

With the rapid adoption of immunotherapy into clinical practice for HPV-associated malignancies, assessing tumor burden using "liquid biopsies" would further our understanding of clinical outcomes mediated by immunotherapy and allow for tailoring of treatment based on real-time tumor dynamics. In this review, we examine translational studies on peripheral surrogates of tumor burden derived from peripheral blood in HPV-associated malignancies, including levels and methylation of circulating tumor DNA (ctDNA), miRNA derived from extracellular vesicles, circulating tumor cells (CTCs), and HPV-specific antibodies and T cell responses. We review their utility as prognostic and predictive biomarkers of response to chemotherapy and radiation, with a focus on how they may inform and guide immunotherapies to treat locally advanced and metastatic HPV-associated malignancies. We also highlight unanswered questions that must be addressed to translate and integrate these peripheral tumor biomarkers into the clinic.

Circulating Tumor DNA in Head and Neck Squamous Cell Carcinoma

Tumors shed cell-free DNA (cfDNA) into the plasma. “Liquid biopsies” are a diagnostic test to analyze cfDNA in order to detect minimal residual cancer, profile the genomic tumor landscape, and monitor cancers non-invasively over time. This technique may be useful in patients with head and neck squamous cell carcinoma (HNSCC) due to genetic tumor heterogeneity and limitations in imaging sensitivity. However, there are technical challenges that need to be overcome for the widespread use of liquid biopsy in the clinical management of these patients. In this review, we discuss our current understanding of HNSCC genetics and the role of cfDNA genomic analyses as an emerging precision diagnostic tool.

Comprehensive Viral Genotyping Reveals Prognostic Viral Phylogenetic Groups in HPV16-Associated Squamous Cell Carcinoma of the Oropharynx

Human papillomavirus-positive (HPV+) squamous cell carcinoma of the oropharynx (OPSCC) is the most prevalent HPV-associated malignancy in the United States and is primarily caused by HPV subtype 16 (HPV16). Favorable treatment outcomes have led to increasing interest in treatment deescalation to reduce treatment-related morbidity. Prognostic biomarkers are needed to identify appropriately low-risk patients for reduced treatment intensity. Targeted DNA sequencing including all HPV16 open reading frames was performed on tumors from 104 patients with HPV16+ OPSCC treated at a single center. Genotypes closely related to the HPV16-A1 reference were associated with increased numbers of somatic copy-number variants in the human genome and poor recurrence-free survival (RFS). Genotypes divergent from HPV16-A1 were associated with favorable RFS. These findings were independent of tobacco smoke exposure. Total RNA sequencing was performed on a second independent cohort of 89 HPV16+ OPSCC cases. HPV16 genotypes divergent from HPV16-A1 were again validated in this independent cohort, to be prognostic of improved RFS in patients with moderate (less than 30 pack-years) or low (no more than 10 pack-years) of tobacco smoke exposure. In summary, we show in two independent cohorts that viral sequence divergence from the HPV16-A1 reference is correlated with improved RFS in patients with moderate or low tobacco smoke exposure.

IMPLICATIONS

HPV16 genotype is a potential biomarker that could be easily adopted to guide therapeutic decision-making related to deescalation therapy.

Direct Comparison of HPV16 Viral Genomic Integration, Copy Loss, and Structural Variants in Oropharyngeal and Uterine Cervical Cancers Reveal Distinct Relationships to E2 Disruption and Somatic Alteration

Squamous cell carcinoma of the oropharynx caused by HPV type 16 (HPV16+ OPSCC) is the most common HPV-associated malignancy in the USA and has many molecular differences from uterine cervical squamous cell carcinoma (UCSCC). Our understanding of HPV oncogenesis relied on studies of UCSCC revealing a consensus model reliant on HPV integration with a loss of E2. Here, we compare patterns of HPV integration in UCSCC and OPSCC by analysis of affinity capture sequencing of the HPV16 genome in 104 OPSCC and 44 UCSCC tumors. These cohorts were contemporaneously sequenced using an identical strategy. Integration was identified using discordant read pair clustering and assembly-based approaches. Viral integration sites, structural variants, and copy losses were examined. While large-scale deep losses of HPV16 genes were common in UCSCC and were associated with E2 loss, deep copy losses of the HPV16 genome were infrequent in HPV16+ OPSCC. Similarly, structural variants within HPV16 favored E2 loss in UCSCC but not OPSCC. HPV16 integration sites were non-random, with recurrent integration hot-spots identified. OPSCC tumors had many more integration sites per tumor when compared to UCSCC and had more integration sites in genomic regions with high gene density. These data show that viral integration and E2 disruption are distinct in UCSCC and OPSCC. Our findings also add to growing literature suggesting that HPV tumorigenesis in OPSCC does not follow the model developed based on UCSCC.

Human Papillomavirus Integration Strictly Correlates with Global Genome Instability in Head and Neck Cancer

Human papillomavirus (HPV)-positive head and neck cancers, predominantly oropharyngeal squamous cell carcinoma (OPSCC), exhibit epidemiologic, clinical, and molecular characteristics distinct from those OPSCCs lacking HPV. We applied a combination of whole-genome sequencing and optical genome mapping to interrogate the genome structure of HPV-positive OPSCCs. We found that the virus had integrated in the host genome in two thirds of the tumors examined but resided solely extrachromosomally in the other third. Integration of the virus occurred at essentially random sites within the genome. Focal amplification of the virus and the genomic sequences surrounding it often occurred subsequent to integration, with the number of tandem repeats in the chromosome accounting for the increased copy number of the genome sequences flanking the site of integration. In all cases, viral integration correlated with pervasive genome-wide somatic alterations at sites distinct from that of viral integration and comprised multiple insertions, deletions, translocations, inversions, and point mutations. Few or no somatic mutations were present in tumors with only episomal HPV. Our data could be interpreted by positing that episomal HPV is captured in the host genome following an episode of global genome instability during tumor development. Viral integration correlated with higher grade tumors, which may be explained by the associated extensive mutation of the genome and suggests that HPV integration status may inform prognosis.

IMPLICATIONS

Our results indicate that HPV integration in head and neck cancer correlates with extensive pangenomic structural variation, which may have prognostic implications.

Prognostic Value of Copy Number Alteration Burden in Early-Stage Breast Cancer and the Construction of an 11-Gene Copy Number Alteration Model

The increasing burden of breast cancer has prompted a wide range of researchers to search for new prognostic markers. Considering that tumor mutation burden (TMB) is low and copy number alteration burden (CNAB) is high in breast cancer, we built a CNAB-based model using a public database and validated it with a Chinese population. We collected formalin-fixed, paraffin-embedded (FFPE) tissue samples from 31 breast cancer patients who were treated between 2010 and 2014 at the National Cancer Center (CICAMS). METABRIC and TCGA data were downloaded via cBioPortal. In total, 2295 patients with early-stage breast cancer were enrolled in the study, including 1427 in the METABRIC cohort, 837 in the TCGA cohort, and 31 in the CICAMS cohort. Based on the ROC curve, we consider 2.2 CNA/MBp as the threshold for the CNAB-high and CNAB-low groupings. In both the TCGA cohort and the CICAMS cohort, CNAB-high had a worse prognosis than CNAB-low. We further simplified this model by establishing a prognostic nomogram for early breast cancer patients by 11 core genes, and this nomogram was highly effective in both the TCGA cohort and the CICAMS cohort. We hope that this model will subsequently help clinicians with prognostic assessments.

Oropharyngeal cancer outcomes correlate with p16 status, multinucleation and immune infiltration

Oropharyngeal squamous cell carcinoma (OPSCC), largely fueled by the human papillomavirus (HPV), has a complex biological and immunologic phenotype. Although HPV/p16 status can be used to stratify OPSCC patients as a function of survival, it remains unclear what drives an improved treatment response in HPV-associated OPSCC and whether targetable biomarkers exist that can inform a precision oncology approach. We analyzed OPSCC patients treated between 2000 and 2016 and correlated locoregional control (LRC), disease-free survival (DFS) and overall survival (OS) with conventional clinical parameters, risk parameters generated using deep-learning algorithms trained to quantify tumor-infiltrating lymphocytes (TILs) (OP-TIL) and multinucleated tumor cells (MuNI) and targeted transcriptomics. P16 was a dominant determinant of LRC, DFS and OS, but tobacco exposure, OP-TIL and MuNI risk features correlated with clinical outcomes independent of p16 status and the combination of p16, OP-TIL and MuNI generated a better stratification of OPSCC risk compared to individual parameters. Differential gene expression (DEG) analysis demonstrated overlap between MuNI and OP-TIL and identified genes involved in DNA repair, oxidative stress response and tumor immunity as the most prominent correlates with survival. Alteration of inflammatory/immune pathways correlated strongly with all risk features and oncologic outcomes. This suggests that development of OPSCC consists of an intersection between multiple required and permissive oncogenic and immunologic events which may be mechanistically linked. The strong relationship between tumor immunity and oncologic outcomes in OPSCC regardless of HPV status may provide opportunities for further biomarker development and precision oncology approaches incorporating immune checkpoint inhibitors for maximal anti-tumor efficacy.