Molecular characterization of the salivary adenoid cystic carcinoma immune landscape by anatomic subsites.

Adenoid cystic carcinoma (AdCC) is a slow-growing salivary gland malignancy that relapses frequently. AdCCs of the submandibular gland exhibit unique differences in prognosis and treatment response to adjuvant radiotherapy compared to other sites, yet the role of tumor anatomic subsite on gene expression and tumor immune microenvironment (TIME) composition remains unclear. We used 87 samples, including 48 samples (27 AdCC and 21 normal salivary gland tissue samples) from 4 publicly available AdCC RNA sequencing datasets, a validation set of 33 minor gland AdCCs, and 39 samples from an in-house cohort (30 AdCC and 9 normal salivary gland samples). RNA sequencing data were used for single sample gene set enrichment analysis and TIME deconvolution. Quantitative PCR and multiplex immunofluorescence were performed on the in-house cohort. Wilcoxon rank-sum, nonparametric equality-of-medians tests and linear regression models were used to evaluate tumor subsite differences. AdCCs of different anatomic subsites including parotid, submandibular, sublingual, and minor salivary glands differed with respect to expression of several key tumorigenic pathways. Among the three major salivary glands, the reactive oxygen species (ROS)/nuclear factor erythroid 2-related factor 2 (NRF2) pathway signature was significantly underexpressed in AdCC of submandibular compared to parotid and sublingual glands while this association was not observed among normal glands. Additionally, the NRF2 pathway, whose expression was associated with favorable overall survival, was overexpressed in AdCCs of parotid gland compared to minor and submandibular glands. The TIME deconvolution identified differences in CD4+ T cell populations between AdCC of major and minor glands and natural killer (NK) cells among AdCC of minor, submandibular, and parotid glands while plasma cells were enriched in normal submandibular glands compared to other normal gland controls. Our data reveal key molecular differences in AdCC of different anatomic subsites. The ROS and NRF2 pathways are underexpressed in submandibular and minor AdCCs compared to parotid gland AdCCs, and NRF2 pathway expression is associated with favorable overall survival. The CD4+ T, NK, and plasma cell populations also vary by tumor subsites, suggesting that the observed submandibular AdCC tumor-intrinsic pathway differences may be responsible for influencing the TIME composition and survival differences.

De-Escalated Therapy and Early Treatment of Recurrences in HPV-Associated Head and Neck Cancer: The Potential for Biomarkers to Revolutionize Personalized Therapy.

Human papillomavirus-associated (HPV+) head and neck squamous cell carcinoma (HNSCC) is the most common HPV-associated cancer in the United States, with a rapid increase in incidence over the last two decades. The burden of HPV+ HNSCC is likely to continue to rise, and given the long latency between infection and the development of HPV+ HNSCC, it is estimated that the effect of the HPV vaccine will not be reflected in HNSCC prevalence until 2060. Efforts have begun to decrease morbidity of standard therapies for this disease, and its improved characterization is being leveraged to identify and target molecular vulnerabilities. Companion biomarkers for new therapies will identify responsive tumors. A more basic understanding of two mechanisms of HPV carcinogenesis in the head and neck has identified subtypes of HPV+ HNSCC that correlate with different carcinogenic programs and that identify tumors with good or poor prognosis. Current development of biomarkers that reliably identify these two subtypes, as well as biomarkers that can detect recurrent disease at an earlier time, will have immediate clinical application.

Noncanonical HPV carcinogenesis drives radiosensitization of head and neck tumors.

We analyzed transcriptional data from 104 HPV+ (Human papillomavirus) HNSCC (head and neck squamous cell carcinoma) tumors together with two publicly available sources to identify highly robust transcriptional programs (modules) which could be detected consistently despite heterogeneous sequencing and quantification methodologies. Among 22 modules identified, we found a single module that naturally subclassifies HPV+ HNSCC tumors based on a bimodal pattern of gene expression, clusters all atypical features of HPV+ HNSCC biology into a single subclass, and predicts patient outcome in four independent cohorts. The subclass-defining gene set was strongly correlated with Nuclear factor kappa B (NF-κB) target expression. Tumors with high expression of this NF-κB module were rarely associated with activating PIK3CA alterations or viral integration, and also expressed higher levels of HPHPV E2 and had decreased APOBEC mutagenesis. Alternatively, they harbored inactivating alterations of key regulators of NF-κB, TNF receptor associated factor 3 (TRAF3), and cylindromatosis (CYLD), as well as retinoblastoma protein (RB1). HPV+ HNSCC cells in culture with experimental depletion of TRAF3 or CYLD displayed increased expression of the subclass-defining genes, as well as robust radio-sensitization, thus recapitulating both the tumor transcriptional state and improved treatment response observed in patient data. Across all gene sets investigated, methylation to expression correlations were the strongest for the subclass-defining, NF-κB-related genes. Increased tumor-infiltrating CD4+ T cells and increased Estrogen receptors alpha (ERα) expression were identified in NF-κB active tumors. Based on the relatively high rates of cure in HPV+ HNSCC, deintensification of therapy to reduce treatment-related morbidity is being studied at many institutions. Tumor subclassification based on oncogenic subtypes may help guide the selection of therapeutic intensity or modality for patients with HPV+ HNSCC.

Erratum: Exploring APOBEC3A and APOBEC3B substrate specificity and their role in HPV positive head and neck cancer.

[This corrects the article DOI: 10.1016/j.isci.2022.105077.].

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.

Exploring ABOBEC3A and APOBEC3B substrate specificity and their role in HPV positive head and neck cancer.

APOBEC3 family members are cytidine deaminases catalyzing conversion of cytidine to uracil. Many studies have established a link between APOBEC3 expression and cancer development and progression, especially APOBEC3A (A3A) and APOBEC3B (A3B). Preclinical studies with human papillomavirus positive (HPV+) head and neck squamous cell carcinoma (HNSCC) and clinical trial specimens revealed induction of A3B, but not A3A expression after demethylation. We examined the kinetic features of the cytidine deaminase activity for full length A3B and found that longer substrates and a purine at -2 position favored by A3B, whereas A3A prefers shorter substrates and an adenine or thymine at -2 position. The importance and biological significance of A3B catalytic activity rather than A3A and a preference for purine at the -2 position was also established in HPV+ HNSCCs. Our study explored factors influencing formation of A3A and A3B-related cancer mutations that are essential for understanding APOBEC3-related carcinogenesis and facilitating drug discovery.

Development of a Novel Molecular Test for Determining HPV Integration Status in HPV-Positive Oropharynx Cancers.

This diagnostic study describes the development of an assay for human papillomavirus­driven cancers of the oropharynx and the role viral integration could play in the process.

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.

Novel Allosteric Mechanism of Dual p53/MDM2 and p53/MDM4 Inhibition by a Small Molecule.

Restoration of the p53 tumor suppressor for personalised cancer therapy is a promising treatment strategy. However, several high-affinity MDM2 inhibitors have shown substantial side effects in clinical trials. Thus, elucidation of the molecular mechanisms of action of p53 reactivating molecules with alternative functional principle is of the utmost importance. Here, we report a discovery of a novel allosteric mechanism of p53 reactivation through targeting the p53 N-terminus which promotes inhibition of both p53/MDM2 (murine double minute 2) and p53/MDM4 interactions. Using biochemical assays and molecular docking, we identified the binding site of two p53 reactivating molecules, RITA (reactivation of p53 and induction of tumor cell apoptosis) and protoporphyrin IX (PpIX). Ion mobility-mass spectrometry revealed that the binding of RITA to serine 33 and serine 37 is responsible for inducing the allosteric shift in p53, which shields the MDM2 binding residues of p53 and prevents its interactions with MDM2 and MDM4. Our results point to an alternative mechanism of blocking p53 interaction with MDM2 and MDM4 and may pave the way for the development of novel allosteric inhibitors of p53/MDM2 and p53/MDM4 interactions.

NF-κB over-activation portends improved outcomes in HPV-associated head and neck cancer.

Evolving understanding of head and neck squamous cell carcinoma (HNSCC) is leading to more specific diagnostic disease classifications. Among HNSCC caused by the human papilloma virus (HPV), tumors harboring defects in TRAF3 or CYLD are associated with improved clinical outcomes and maintenance of episomal HPV. TRAF3 and CYLD are negative regulators of NF-κB and inactivating mutations of either leads to NF-κB overactivity. Here, we developed and validated a gene expression classifier separating HPV+ HNSCCs based on NF-κB activity. As expected, the novel classifier is strongly enriched in NF-κB targets leading us to name it the NF-κB Activity Classifier (NAC). High NF-κB activity correlated with improved survival in two independent cohorts. Using NAC, tumors with high NF-κB activity but lacking defects in TRAF3 or CYLD were identified; thus, while TRAF3 or CYLD gene defects identify the majority of tumors with NF-κB activation, unknown mechanisms leading to NF-kB activity also exist. The NAC correctly classified the functional consequences of two novel CYLD missense mutations. Using a reporter assay, we tested these CYLD mutations revealing that their activity to inhibit NF-kB was equivalent to the wild-type protein. Future applications of the NF-κB Activity Classifier may be to identify HPV+ HNSCC patients with better or worse survival with implications for treatment strategies.

Dual Targeting of Chromatin Stability By The Curaxin CBL0137 and Histone Deacetylase Inhibitor Panobinostat Shows Significant Preclinical Efficacy in Neuroblastoma.

PURPOSE: We investigated whether targeting chromatin stability through a combination of the curaxin CBL0137 with the histone deacetylase (HDAC) inhibitor, panobinostat, constitutes an effective multimodal treatment for high-risk neuroblastoma. EXPERIMENTAL DESIGN: The effects of the drug combination on cancer growth were examined in vitro and in animal models of MYCN-amplified neuroblastoma. The molecular mechanisms of action were analyzed by multiple techniques including whole transcriptome profiling, immune deconvolution analysis, immunofluorescence, flow cytometry, pulsed-field gel electrophoresis, assays to assess cell growth and apoptosis, and a range of cell-based reporter systems to examine histone eviction, heterochromatin transcription, and chromatin compaction. RESULTS: The combination of CBL0137 and panobinostat enhanced nucleosome destabilization, induced an IFN response, inhibited DNA damage repair, and synergistically suppressed cancer cell growth. Similar synergistic effects were observed when combining CBL0137 with other HDAC inhibitors. The CBL0137/panobinostat combination significantly delayed cancer progression in xenograft models of poor outcome high-risk neuroblastoma. Complete tumor regression was achieved in the transgenic Th-MYCN neuroblastoma model which was accompanied by induction of a type I IFN and immune response. Tumor transplantation experiments further confirmed that the presence of a competent adaptive immune system component allowed the exploitation of the full potential of the drug combination. CONCLUSIONS: The combination of CBL0137 and panobinostat is effective and well-tolerated in preclinical models of aggressive high-risk neuroblastoma, warranting further preclinical and clinical investigation in other pediatric cancers. On the basis of its potential to boost IFN and immune responses in cancer models, the drug combination holds promising potential for addition to immunotherapies.

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.

BACKGROUND: Human papillomavirus-positive (HPV+) squamous cell carcinoma of the oropharynx (OPSCC) is the most prevalent HPV-associated malignancy in the United States. Favorable treatment outcomes have led to increased interest in treatment de-escalation to reduce treatment morbidity as well as the development of prognostic markers to identify appropriately low-risk patients. Intratumoral genomic heterogeneity and copy number alteration burden have been demonstrated to be predictive of poor outcomes in many other cancers; therefore, we sought to determine whether intratumor heterogeneity and genomic instability are associated with poor outcomes in HPV+ OPSCC. METHODS: Tumor heterogeneity estimates were made based on targeted exome sequencing of 45 patients with HPV+ OPSCC tumors. Analysis of an additional cohort of HPV+ OPSCC tumors lacking matched normal sequencing allowed copy number analysis of 99 patient tumors. RESULTS: High intratumorally genomic heterogeneity and high numbers of copy number alterations were strongly associated with worse recurrence-free survival. Tumors with higher heterogeneity and frequent copy number alterations were associated with loss of distal 11q, which encodes key genes related to double-strand break repair, including ATM and MRE11A. CONCLUSIONS: Both intratumor genomic heterogeneity and high-burden copy number alterations are strongly associated with poor recurrence-free survival in patients with HPV+ OPSCC. The drivers of genomic instability and heterogeneity in these tumors remains to be elucidated. However, 11q loss and defective DNA double-strand break repair have been associated with genomic instability in other solid tumors. Copy number alteration burden and intratumoral heterogeneity represent promising avenues for risk stratification of patients with HPV+OPSCC.

Global Genome Demethylation Causes Transcription-Associated DNA Double Strand Breaks in HPV-Associated Head and Neck Cancer Cells.

High levels of DNA methylation at CpG loci are associated with transcriptional repression of tumor suppressor genes and dysregulation of DNA repair genes. Human papilloma virus (HPV)-associated head and neck squamous cell carcinomas (HNSCC) have high levels of DNA methylation and methylation has been associated with dampening of an innate immune response in virally infected cells. We have been exploring demethylation as a potential treatment in HPV+ HNSCC and recently reported results of a window clinical trial showing that HNSCCs are particularly sensitive to demethylating agent 5-azacytidine (5-aza). Mechanistically, sensitivity is partially due to downregulation of HPV genes expression and restoration of tumor suppressors p53 and Rb. Here, for the first time, we show that 5-azaC treatment of HPV+ HNSCC induces replication and transcription-associated DNA double strand breaks (DSBs) that occur preferentially at demethylated genomic DNA. Blocking replication or transcription prevented formation of DNA DSBs and reduced sensitivity of HPV-positive head and neck cancer cells to 5-azaC, demonstrating that both replication and active transcription are required for formation of DSBs associated with 5-azaC.

NSD1 mutations by HPV status in head and neck cancer: differences in survival and response to DNA-damaging agents.

BACKGROUND: Compared to HPV-negative head and neck squamous cell carcinomas (HNSCCs), HPV-positive HNSCCs are associated with a favorable prognosis in part due to their improved treatment sensitivity. Inactivating mutations in NSD1 were shown to be a favorable prognostic biomarker in laryngeal cancers. Here, we characterize NSD1 mutations from the expanded The Cancer Genome Atlas (TCGA) HNSCC cohort (n = 522) and examine their prognostic implications based on HPV status of the tumor. We also begin to examine if NSD1 regulates response to platinum-based drugs and other DNA-damaging agents. METHODS: TCGA HNSCC samples were segregated by HPV and NSD1 mutations using cBioPortal and patient survival was determined. Pathogenicity of mutations was predicted using UMD-Predictor. NSD1-depleted cell lines were established by transfection with control or shRNAs against NSD1, followed by puromycin selection, and confirmed by qRT-PCR. Cell sensitivity to DNA damaging agents was assessed using short-term proliferation and long-term clonogenic survival assays. RESULTS: Among 457 HPV(-) tumors, 13% contained alterations in the NSD1 gene. The majority (61.3%) of NSD1 gene alterations in HPV(-) specimens were truncating mutations within or before the enzymatic SET domain. The remaining alterations included homozygous gene deletions (6.7%), missense point mutations (30.7%) and inframe deletions (1.3%). UMD-Predictor categorized 18 of 23 missense point mutations as pathogenic. For HPV(+) HNSCC (n = 65), 6 NSD1 mutations, comprised of two truncating (33%) and 4 missense point (66%) mutations, were identified. Three of the 4 missense point mutations were predicted to be pathogenic or probably pathogenic by UMD-Predictor. Kaplan-Meier survival analysis determined significantly improved survival of HPV(-) HNSCC patients whose tumors harbored NSD1 gene alterations, as compared to patients with wild-type NSD1 tumors. Interestingly, the survival effect of NSD1 mutations observed in HPV-negative HNSCC was reversed in HPV(+) tumors. Proliferation and clonogenic survival of two HPV(-) cell lines stably expressing control or NSD1 shRNAs showed that NSD1-depleted cells were more sensitive to cisplatin and carboplatin, but not to other DNA damaging drugs. CONCLUSIONS: Genetic alterations in NSD1 hold potential as novel prognostic biomarkers in HPV(-) head and neck cancers. NSD1 mutations in HPV(+) cancers may also play a prognostic role, although this effect must be examined in a larger cohort. NSD1 downregulation results in improved sensitivity to cisplatin and carboplatin, but not to other DNA-damaging agents, in epithelial cells. Increased sensitivity to platinum-based chemotherapy agents associated with NSD1 depletion may contribute to improved survival in HPV(-) HNSCCs. Further studies are needed to determine mechanisms through which NSD1 protects HPV(-) HNSCC cells from platinum-based therapy, as well as confirmation of NSD1 effect in HPV(+) HNSCC.

p53 Signaling in Cancers.

This special issue on p53 explores different aspects of the significance of p53 in normal cells and in cancer [...].

Combined Aurora Kinase A (AURKA) and WEE1 Inhibition Demonstrates Synergistic Antitumor Effect in Squamous Cell Carcinoma of the Head and Neck.

PURPOSE: Human papillomavirus (HPV)-negative head and neck squamous cell carcinomas (HNSCC) commonly bear disruptive mutations in TP53, resulting in treatment resistance. In these patients, direct targeting of p53 has not been successful, but synthetic lethal approaches have promise. Although Aurora A kinase (AURKA) is overexpressed and an oncogenic driver, its inhibition has only modest clinical effects in HPV-negative HNSCC. We explored a novel combination of AURKA and WEE1 inhibition to overcome intrinsic resistance to AURKA inhibition.Experimental Design: AURKA protein expression was determined by fluorescence-based automated quantitative analysis of patient specimens and correlated with survival. We evaluated treatment with the AURKA inhibitor alisertib (MLN8237) and the WEE1 inhibitor adavosertib (AZD1775), alone or in combination, using in vitro and in vivo HNSCC models. RESULTS: Elevated nuclear AURKA correlated with worse survival among patients with p16(-) HNSCC. Alisertib caused spindle defects, G2-M arrest and inhibitory CDK1 phosphorylation, and cytostasis in TP53 mutant HNSCC FaDu and UNC7 cells. Addition of adavosertib to alisertib instead triggered mitotic entry and mitotic catastrophe. Moreover, in FaDu and Detroit 562 xenografts, this combination demonstrated synergistic effects on tumor growth and extended overall survival compared with either vehicle or single-agent treatment. CONCLUSIONS: Combinatorial treatment with adavosertib and alisertib leads to synergistic antitumor effects in in vitro and in vivo HNSCC models. These findings suggest a novel rational combination, providing a promising therapeutic avenue for TP53-mutated cancers.

APOBEC-induced mutations and their cancer effect size in head and neck squamous cell carcinoma.

Recent studies have revealed the mutational signatures underlying the somatic evolution of cancer, and the prevalences of associated somatic genetic variants. Here we estimate the intensity of positive selection that drives mutations to high frequency in tumors, yielding higher prevalences than expected on the basis of mutation and neutral drift alone. We apply this approach to a sample of 525 head and neck squamous cell carcinoma exomes, producing a rank-ordered list of gene variants by selection intensity. Our results illustrate the complementarity of calculating the intensity of selection on mutations along with tallying the prevalence of individual substitutions in cancer: while many of the most prevalently-altered genes were heavily selected, their relative importance to the cancer phenotype differs from their prevalence and from their P value, with some infrequent variants exhibiting evidence of strong positive selection. Furthermore, we extend our analysis of effect size by quantifying the degree to which mutational processes (such as APOBEC mutagenesis) contributes mutations that are highly selected, driving head and neck squamous cell carcinoma. We calculate the substitutions caused by APOBEC mutagenesis that make the greatest contribution to cancer phenotype among patients. Lastly, we demonstrate via in vitro biochemical experiments that the APOBEC3B protein can deaminate the cytosine bases at two sites whose mutant states are subject to high net realized selection intensities-PIK3CA E545K and E542K. By quantifying the effects of mutations, we deepen the molecular understanding of carcinogenesis in head and neck squamous cell carcinoma.

Genomic and Transcriptomic Characterization Links Cell Lines with Aggressive Head and Neck Cancers.

Cell lines are important tools for biological and preclinical investigation, and establishing their relationship to genomic alterations in tumors could accelerate functional and therapeutic discoveries. We conducted integrated analyses of genomic and transcriptomic profiles of 15 human papillomavirus (HPV)-negative and 11 HPV-positive head and neck squamous cell carcinoma (HNSCC) lines to compare with 279 tumors from The Cancer Genome Atlas (TCGA). We identified recurrent amplifications on chromosomes 3q22-29, 5p15, 11q13/22, and 8p11 that drive increased expression of more than 100 genes in cell lines and tumors. These alterations, together with loss or mutations of tumor suppressor genes, converge on important signaling pathways, recapitulating the genomic landscape of aggressive HNSCCs. Among these, concurrent 3q26.3 amplification and TP53 mutation in most HPV(-) cell lines reflect tumors with worse survival. Our findings elucidate and validate genomic alterations underpinning numerous discoveries made with HNSCC lines and provide valuable models for future studies.