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.

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.

Hydroxyurea-stalled replication forks become progressively inactivated and require two different RAD51-mediated pathways for restart and repair.

Faithful DNA replication is essential to all life. Hydroxyurea (HU) depletes the cells of dNTPs, which initially results in stalled replication forks that, after prolonged treatment, collapse into DSBs. Here, we report that stalled replication forks are efficiently restarted in a RAD51-dependent process that does not trigger homologous recombination (HR). The XRCC3 protein, which is required for RAD51 foci formation, is also required for replication restart of HU-stalled forks, suggesting that RAD51-mediated strand invasion supports fork restart. In contrast, replication forks collapsed by prolonged replication blocks do not restart, and global replication is rescued by new origin firing. We find that RAD51-dependent HR is triggered for repair of collapsed replication forks, without apparent restart. In conclusion, our data suggest that restart of stalled replication forks and HR repair of collapsed replication forks require two distinct RAD51-mediated pathways.

TRAF3/CYLD mutations identify a distinct subset of human papillomavirus-associated head and neck squamous cell carcinoma.

BACKGROUND: The incidence of human papillomavirus (HPV)-associated (HPV-positive) head and neck squamous cell carcinoma (HNSCC) of the oropharynx has dramatically increased over the last decade and continues to rise. Newly diagnosed HPV-positive HNSCCs in the United States currently outnumber any other HPV-associated cancers, including cervical cancer. Despite introduction of the HPV vaccine, the epidemic of HPV-positive HNSCC is expected to continue for approximately 60 years. Compared with patients who have tobacco-associated HNSCC, those who have HPV-positive HNSCC have better overall survival and response to treatment. Current treatment, including chemotherapy and radiation therapy, is associated with lifelong morbidity, and there are limited treatments and no curative options for patients who develop recurrent metastatic disease. Therapeutic de-escalation (decreased radiation dose) is being tested through clinical trials; however, those studies select patients based solely on tumor and patient smoking characteristics. Mechanisms of HPV-driven carcinogenesis in HNSCC are not well understood, which limits new therapeutic strategies and hinders the appropriate selection of patients for de-escalation therapy. METHODS: The authors analyzed HNSCC data from The Cancer Genome Atlas to identify molecular characteristics that correlate with outcomes and integration status of the HPV genome. RESULTS: The current investigations identified a subset of HPV-positive HNSCCs with mutations in the genes TRAF3 (tumor necrosis factor receptor-associated factor 3) and CYLD (cylindromatosis lysine 63 deubiquitinase). Defects in TRAF3 and CYLD correlated with the activation of transcriptional factor nuclear factor κB, episomal HPV status of tumors, and improved patient survival. CONCLUSIONS: Defects in TRAF3/CYLD were accompanied with the activation of nuclear factor κB signaling and maintenance of episomal HPV in tumors, suggesting that these mutations may support an alternative mechanism of HPV tumorigenesis in head and neck tumors. Cancer 2017;123:1778-1790. © 2017 The Authors. Cancer published by Wiley Periodicals, Inc. on behalf of American Cancer Society. This is an open access article under the terms of the Creative Commons Attribution NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.

Characterization of HPV and host genome interactions in primary head and neck cancers.

Previous studies have established that a subset of head and neck tumors contains human papillomavirus (HPV) sequences and that HPV-driven head and neck cancers display distinct biological and clinical features. HPV is known to drive cancer by the actions of the E6 and E7 oncoproteins, but the molecular architecture of HPV infection and its interaction with the host genome in head and neck cancers have not been comprehensively described. We profiled a cohort of 279 head and neck cancers with next generation RNA and DNA sequencing and show that 35 (12.5%) tumors displayed evidence of high-risk HPV types 16, 33, or 35. Twenty-five cases had integration of the viral genome into one or more locations in the human genome with statistical enrichment for genic regions. Integrations had a marked impact on the human genome and were associated with alterations in DNA copy number, mRNA transcript abundance and splicing, and both inter- and intrachromosomal rearrangements. Many of these events involved genes with documented roles in cancer. Cancers with integrated vs. nonintegrated HPV displayed different patterns of DNA methylation and both human and viral gene expressions. Together, these data provide insight into the mechanisms by which HPV interacts with the human genome beyond expression of viral oncoproteins and suggest that specific integration events are an integral component of viral oncogenesis.

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.

PARP is activated at stalled forks to mediate Mre11-dependent replication restart and recombination.

If replication forks are perturbed, a multifaceted response including several DNA repair and cell cycle checkpoint pathways is activated to ensure faithful DNA replication. Here, we show that poly(ADP-ribose) polymerase 1 (PARP1) binds to and is activated by stalled replication forks that contain small gaps. PARP1 collaborates with Mre11 to promote replication fork restart after release from replication blocks, most likely by recruiting Mre11 to the replication fork to promote resection of DNA. Both PARP1 and PARP2 are required for hydroxyurea-induced homologous recombination to promote cell survival after replication blocks. Together, our data suggest that PARP1 and PARP2 detect disrupted replication forks and attract Mre11 for end processing that is required for subsequent recombination repair and restart of replication forks.

Tumour suppressor Fus1 provides a molecular link between inflammatory response and mitochondrial homeostasis.

Fus1, encoded by a 3p21.3 tumour suppressor gene, is down-regulated, mutated or lost in the majority of inflammatory thoracic malignancies. The mitochondrial localization of Fus1 stimulated us to investigate how Fus1 modulates inflammatory response and mitochondrial function in a mouse model of asbestos-induced peritoneal inflammation. Asbestos treatment resulted in a decreased Fus1 expression in wild-type (WT) peritoneal immune cells, suggesting that asbestos exposure may compromise the Fus1-mediated inflammatory response. Untreated Fus1(-/-) mice had an ~eight-fold higher proportion of peritoneal granulocytes than Fus1(+/+) mice, pointing at ongoing chronic inflammation. Fus1(-/-) mice exhibited a perturbed inflammatory response to asbestos, reflected in decreased immune organ weight and peritoneal fluid protein concentration, along with an increased proportion of peritoneal macrophages. Fus1(-/-) immune cells showed augmented asbestos-induced activation of key inflammatory, anti-oxidant and genotoxic stress response proteins ERK1/2, NFκB, SOD2, γH2AX, etc. Moreover, Fus1(-/-) mice demonstrated altered dynamics of pro- and anti-inflammatory cytokine expression, such as IFNγ, TNFα, IL-1A, IL-1B and IL-10. 'Late' response cytokine Ccl5 was persistently under-expressed in Fus1(-/-) immune cells at both basal and asbestos-activated states. We observed an asbestos-related difference in the size of CD3(+) CD4(-) CD8(-) DN T cell subset that was expanded four-fold in Fus1(-/-) mice. Finally, we demonstrated Fus1-dependent basal and asbestos-induced changes in major mitochondrial parameters (ROS production, mitochondrial potential and UCP2 expression) in Fus1(-/-) immune cells and in Fus1-depleted cancer cells, thus supporting our hypothesis that Fus1 establishes its immune- and tumour-suppressive activities via regulation of mitochondrial homeostasis.

Oncogene-induced senescence is part of the tumorigenesis barrier imposed by DNA damage checkpoints.

Recent studies have indicated the existence of tumorigenesis barriers that slow or inhibit the progression of preneoplastic lesions to neoplasia. One such barrier involves DNA replication stress, which leads to activation of the DNA damage checkpoint and thereby to apoptosis or cell cycle arrest, whereas a second barrier is mediated by oncogene-induced senescence. The relationship between these two barriers, if any, has not been elucidated. Here we show that oncogene-induced senescence is associated with signs of DNA replication stress, including prematurely terminated DNA replication forks and DNA double-strand breaks. Inhibiting the DNA double-strand break response kinase ataxia telangiectasia mutated (ATM) suppressed the induction of senescence and in a mouse model led to increased tumour size and invasiveness. Analysis of human precancerous lesions further indicated that DNA damage and senescence markers cosegregate closely. Thus, senescence in human preneoplastic lesions is a manifestation of oncogene-induced DNA replication stress and, together with apoptosis, provides a barrier to malignant progression.

Tumor suppressor Lzap regulates cell cycle progression, doming, and zebrafish epiboly.

Initial stages of embryonic development rely on rapid, synchronized cell divisions of the fertilized egg followed by a set of morphogenetic movements collectively called epiboly and gastrulation. Lzap is a putative tumor suppressor whose expression is lost in 30% of head and neck squamous cell carcinomas. Lzap activities include regulation of cell cycle progression and response to therapeutic agents. Here, we explore developmental roles of the lzap gene during zebrafish morphogenesis. Lzap is highly conserved among vertebrates and is maternally deposited. Expression is initially ubiquitous during gastrulation, and later becomes more prominent in the pharyngeal arches, digestive tract, and brain. Antisense morpholino-mediated depletion of Lzap resulted in delayed cell divisions and apoptosis during blastomere formation, resulting in fewer, larger cells. Cell cycle analysis suggested that Lzap loss in early embryonic cells resulted in a G2/M arrest. Furthermore, the Lzap-deficient embryos failed to initiate epiboly--the earliest morphogenetic movement in animal development--which has been shown to be dependent on cell adhesion and migration of epithelial sheets. Our results strongly implicate Lzap in regulation of cell cycle progression, adhesion and migratory activity of epithelial cell sheets during early development. These functions provide further insight into Lzap activity that may contribute not only to development, but also to tumor formation.

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.

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.

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.

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.].

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.

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.

Small molecule RITA binds to p53, blocks p53-HDM-2 interaction and activates p53 function in tumors.

In tumors that retain wild-type p53, its tumor-suppressor function is often impaired as a result of the deregulation of HDM-2, which binds to p53 and targets it for proteasomal degradation. We have screened a chemical library and identified a small molecule named RITA (reactivation of p53 and induction of tumor cell apoptosis), which bound to p53 and induced its accumulation in tumor cells. RITA prevented p53-HDM-2 interaction in vitro and in vivo and affected p53 interaction with several negative regulators. RITA induced expression of p53 target genes and massive apoptosis in various tumor cells lines expressing wild-type p53. RITA suppressed the growth of human fibroblasts and lymphoblasts only upon oncogene expression and showed substantial p53-dependent antitumor effect in vivo. RITA may serve as a lead compound for the development of an anticancer drug that targets tumors with wild-type p53.

Restoration of the tumor suppressor function to mutant p53 by a low-molecular-weight compound.

The tumor suppressor p53 inhibits tumor growth primarily through its ability to induce apoptosis. Mutations in p53 occur in at least 50% of human tumors. We hypothesized that reactivation of mutant p53 in such tumors should trigger massive apoptosis and eliminate the tumor cells. To test this, we screened a library of low-molecular-weight compounds in order to identify compounds that can restore wild-type function to mutant p53. We found one compound capable of inducing apoptosis in human tumor cells through restoration of the transcriptional transactivation function to mutant p53. This molecule, named PRIMA-1, restored sequence-specific DNA binding and the active conformation to mutant p53 proteins in vitro and in living cells. PRIMA-1 rescued both DNA contact and structural p53 mutants. In vivo studies in mice revealed an antitumor effect with no apparent toxicity. This molecule may serve as a lead compound for the development of anticancer drugs targeting mutant p53.

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.