To achieve the very high oncoprotein levels required to drive the malignant state cancer cells utilise the ubiquitin proteasome system to upregulate transcription factor levels. Here our analyses identify ALYREF, expressed from the most common genetic copy number variation in neuroblastoma, chromosome 17q21-ter gain as a key regulator of MYCN protein turnover. We show strong co-operativity between ALYREF and MYCN from transgenic models of neuroblastoma in vitro and in vivo. The two proteins form a nuclear coactivator complex which stimulates transcription of the ubiquitin specific peptidase 3, USP3. We show that increased USP3 levels reduce K-48- and K-63-linked ubiquitination of MYCN, thus driving up MYCN protein stability. In the MYCN-ALYREF-USP3 signal, ALYREF is required for MYCN effects on the malignant phenotype and that of USP3 on MYCN stability. This data defines a MYCN oncoprotein dependency state which provides a rationale for future pharmacological studies.

Peripheral neuroblastic tumors (PNTs) represent a spectrum of neural-crest-derived tumors, including neuroblastoma, ganglioneuroblastoma, and ganglioneuroma. Malignant cells in PNTs are theorized to interconvert between adrenergic/noradrenergic and mesenchymal/neural crest cell states. Here, single-cell RNA-sequencing analysis of 10 PNTs demonstrates extensive transcriptomic heterogeneity. Trajectory modeling suggests that malignant neuroblasts move between adrenergic and mesenchymal cell states via an intermediate state that we term "transitional." Transitional cells express programs linked to a sympathoadrenal development and aggressive tumor phenotypes such as rapid proliferation and tumor dissemination. Among primary bulk tumor patient cohorts, high expression of the transitional gene signature is predictive of poor prognosis compared with adrenergic and mesenchymal expression patterns. High transitional gene expression in neuroblastoma cell lines identifies a similar transitional H3K27-acetylation super-enhancer landscape. Collectively, our study supports the concept that PNTs have phenotypic plasticity and uncovers potential biomarkers and therapeutic targets.

Histone deacetylase (HDAC) inhibitors are effective in MYCN-driven cancers, because of a unique need for HDAC recruitment by the MYCN oncogenic signal. However, HDAC inhibitors are much more effective in combination with other anti-cancer agents. To identify novel compounds which act synergistically with HDAC inhibitor, such as suberanoyl hydroxamic acid (SAHA), we performed a cell-based, high-throughput drug screen of 10,560 small molecule compounds from a drug-like diversity library and identified a small molecule compound (SE486-11) which synergistically enhanced the cytotoxic effects of SAHA. Effects of drug combinations on cell viability, proliferation, apoptosis and colony forming were assessed in a panel of neuroblastoma cell lines. Treatment with SAHA and SE486-11 increased MYCN ubiquitination and degradation, and markedly inhibited tumorigenesis in neuroblastoma xenografts, and, MYCN transgenic zebrafish and mice. The combination reduced ubiquitin-specific protease 5 (USP5) levels and increased unanchored polyubiquitin chains. Overexpression of USP5 rescued neuroblastoma cells from the cytopathic effects of the combination and reduced unanchored polyubiquitin, suggesting USP5 is a therapeutic target of the combination. SAHA and SE486-11 directly bound to USP5 and the drug combination exhibited a 100-fold higher binding to USP5 than individual drugs alone in microscale thermophoresis assays. MYCN bound to the USP5 promoter and induced USP5 gene expression suggesting that USP5 and MYCN expression created a forward positive feedback loop in neuroblastoma cells. Thus, USP5 acts as an oncogenic cofactor with MYCN in neuroblastoma and the novel combination of HDAC inhibitor with SE486-11 represents a novel therapeutic approach for the treatment of MYCN-driven neuroblastoma.

Many of the pro-tumorigenic functions of the oncogene MYCN are attributed to its regulation of global gene expression programs. Alternative splicing is another important regulator of gene expression and has been implicated in neuroblastoma development, however, the molecular mechanisms remain unknown. We found that MYCN up-regulated the expression of the core spliceosomal protein, SNRPD3, in models of neuroblastoma initiation and progression. High mRNA expression of SNRPD3 in human neuroblastoma tissues was a strong, independent prognostic factor for poor patient outcome. Repression of SNRPD3 expression correlated with loss of colony formation in vitro and reduced tumorigenicity in vivo. The effect of SNRPD3 on cell viability was in part dependent on MYCN as an oncogenic co-factor. RNA-sequencing revealed a global increase in the number of genes being differentially spliced when MYCN was overexpressed. Surprisingly, depletion of SNRPD3 in the presence of overexpressed MYCN further increased differential splicing, particularly of cell cycle regulators, such as BIRC5 and CDK10. MYCN directly bound SNRPD3, and the protein arginine methyltransferase, PRMT5, consequently increasing SNRPD3 methylation. Indeed, the PRMT5 inhibitor, JNJ-64619178, reduced cell viability and SNRPD3 methylation in neuroblastoma cells with high SNRPD3 and MYCN expression. Our findings demonstrate a functional relationship between MYCN and SNRPD3, which maintains the fidelity of MYCN-driven alternative splicing in the narrow range required for neuroblastoma cell growth. SNRPD3 methylation and its protein-protein interface with MYCN represent novel therapeutic targets.

MYCN oncogene amplification is frequently observed in aggressive childhood neuroblastoma. Using an unbiased large-scale mutagenesis screen in neuroblastoma-prone transgenic mice, we identify a single germline point mutation in the transcriptional corepressor Runx1t1, which abolishes MYCN-driven tumorigenesis. This loss-of-function mutation disrupts a highly conserved zinc finger domain within Runx1t1. Deletion of one Runx1t1 allele in an independent Runx1t1 knockout mouse model is also sufficient to prevent MYCN-driven neuroblastoma development, and reverse ganglia hyperplasia, a known pre-requisite for tumorigenesis. Silencing RUNX1T1 in human neuroblastoma cells decreases colony formation in vitro, and inhibits tumor growth in vivo. Moreover, RUNX1T1 knockdown inhibits the viability of PAX3-FOXO1 fusion-driven rhabdomyosarcoma and MYC-driven small cell lung cancer cells. Despite the role of Runx1t1 in MYCN-driven tumorigenesis neither gene directly regulates the other. We show RUNX1T1 forms part of a transcriptional LSD1-CoREST3-HDAC repressive complex recruited by HAND2 to enhancer regions to regulate chromatin accessibility and cell-fate pathway genes.

Peripheral neuroblastic tumors (PNTs) are the most common extracranial solid tumors in early childhood. They represent a spectrum of neural crest derived tumors including neuroblastoma, ganglioneuroblastoma and ganglioneuroma. PNTs exhibit heterogeneity due to interconverting malignant cell states described as adrenergic/nor-adrenergic or mesenchymal/neural crest cell in origin. The factors determining individual patient levels of tumor heterogeneity, their impact on the malignant phenotype, and the presence of other cell states are unknown. Here, single-cell RNA-sequencing analysis of 4267 cells from 7 PNTs demonstrated extensive transcriptomic heterogeneity. Trajectory modelling showed that malignant neuroblasts move between adrenergic and mesenchymal cell states via a novel state that we termed a “transitional” phenotype. Transitional cells are characterized by gene expression programs linked to a sympathoadrenal development, and aggressive tumor phenotypes such as rapid proliferation and tumor dissemination. Among primary bulk tumor patient cohorts, high expression of the transitional gene signature was highly predictive of poor prognosis when compared to adrenergic and mesenchymal expression patterns. High transitional gene expression in neuroblastoma cell lines identified a similar transitional H3K27-acetylation super-enhancer landscape, supporting the concept that PNTs have phenotypic plasticity and transdifferentiation capacity. Additionally, examination of PNT microenvironments, found that neuroblastomas contained low immune cell infiltration, high levels of non-inflammatory macrophages, and low cytotoxic T lymphocyte levels compared with more benign PNT subtypes. Modeling of cell-cell signaling in the tumor microenvironment predicted specific paracrine effects toward the various subtypes of malignant cells, suggesting further cell-extrinsic influences on malignant cell phenotype. Collectively, our study reveals the presence of a previously unrecognized transitional cell state with high malignant potential and an immune cell architecture which serve both as potential biomarkers and therapeutic targets.