TERT gene rearrangement with transcriptional superenhancers leads to TERT overexpression and neuroblastoma. No targeted therapy is available for clinical trials in patients with TERT-rearranged neuroblastoma.

Glucocorticoids play a critical role in the treatment of lymphoid malignancies. While glucocorticoid efficacy can be largely attributed to lymphocyte-specific apoptosis, its molecular basis remains elusive. Here, we studied genome-wide lymphocyte-specific open chromatin domains (LSOs), and integrated LSOs with glucocorticoid-induced RNA transcription and chromatin modulation using an in vivo patient-derived xenograft model of acute lymphoblastic leukemia (ALL). This led to the identification of LSOs critical for glucocorticoid-induced apoptosis. Glucocorticoid receptor cooperated with CTCF at these LSOs to mediate DNA looping, which was inhibited by increased DNA methylation in glucocorticoid-resistant ALL and non-lymphoid cell types. Our study demonstrates that lymphocyte-specific epigenetic modifications pre-determine glucocorticoid resistance in ALL and may account for the lack of glucocorticoid sensitivity in other cell types.

Acute lymphoblastic leukaemia (ALL) is the most common paediatric malignancy. Glucocorticoids form a critical component of chemotherapy regimens and resistance to glucocorticoid therapy is predictive of poor outcome. We have previously shown that glucocorticoid resistance is associated with upregulation of the oncogene C-MYC and failure to induce the proapoptotic gene BIM.

Diffuse intrinsic pontine glioma (DIPG) is an incurable malignant childhood brain tumor, with no active systemic therapies and a 5-year survival of less than 1%. Polyamines are small organic polycations that are essential for DNA replication, translation and cell proliferation. Ornithine decarboxylase 1 (ODC1), the rate-limiting enzyme in polyamine synthesis, is irreversibly inhibited by difluoromethylornithine (DFMO). Herein we show that polyamine synthesis is upregulated in DIPG, leading to sensitivity to DFMO. DIPG cells compensate for ODC1 inhibition by upregulation of the polyamine transporter SLC3A2. Treatment with the polyamine transporter inhibitor AMXT 1501 reduces uptake of polyamines in DIPG cells, and co-administration of AMXT 1501 and DFMO leads to potent in vitro activity, and significant extension of survival in three aggressive DIPG orthotopic animal models. Collectively, these results demonstrate the potential of dual targeting of polyamine synthesis and uptake as a therapeutic strategy for incurable DIPG.

The ability of the N-MYC transcription factor to drive cancer progression is well demonstrated in neuroblastoma, the most common extracranial pediatric solid tumor, where MYCN amplification heralds a poor prognosis, with only 11% of high-risk patients surviving past 5 years. However, decades of attempts of direct inhibition of N-MYC or its paralogues has led to the conclusion that this protein is "undruggable." Therefore, targeting pathways upregulated by N-MYC signaling presents an alternative therapeutic approach. Here, we show that MYCN-amplified neuroblastomas are characterized by elevated rates of protein synthesis and that high expression of ABCE1, a translation factor directly upregulated by N-MYC, is itself a strong predictor of poor clinical outcome. Despite the potent ability of N-MYC in heightening protein synthesis and malignant characteristics in cancer cells, suppression of ABCE1 alone selectively negated this effect, returning the rate of translation to baseline levels and significantly reducing the growth, motility, and invasiveness of MYCN-amplified neuroblastoma cells and patient-derived xenograft tumors in vivo. The growth of nonmalignant cells or MYCN-nonamplified neuroblastoma cells remained unaffected by reduced ABCE1, supporting a therapeutic window associated with targeting ABCE1. Neuroblastoma cells with c-MYC overexpression also required ABCE1 to maintain cell proliferation and translation. Taken together, ABCE1-mediated translation constitutes a critical process in the progression of N-MYC-driven and c-MYC-driven cancers that warrants investigations into methods of its therapeutic inhibition. SIGNIFICANCE: These findings demonstrate that N-MYC-driven cancers are reliant on elevated rates of protein synthesis driven by heightened expression of ABCE1, a vulnerability that can be exploited through suppression of ABCE1.

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.

Diffuse midline gliomas (DMG) harbouring H3K27M mutation are paediatric tumours with a dismal outcome. Recently, a new subtype of midline gliomas has been described with similar features to DMG, including loss of H3K27 trimethylation, but lacking the canonical H3K27M mutation (H3-WT). Here, we report a cohort of five H3-WT tumours profiled by whole-genome sequencing, RNA sequencing and DNA methylation profiling and combine their analysis with previously published cases. We show that these tumours have recurrent and mutually exclusive mutations in either ACVR1 or EGFR and are characterised by high expression of EZHIP associated to its promoter hypomethylation. Affected patients share a similar poor prognosis as patients with H3K27M DMG. Global molecular analysis of H3-WT and H3K27M DMG reveal distinct transcriptome and methylome profiles including differential methylation of homeobox genes involved in development and cellular differentiation. Patients have distinct clinical features, with a trend demonstrating ACVR1 mutations occurring in H3-WT tumours at an older age. This in-depth exploration of H3-WT tumours further characterises this novel DMG, H3K27-altered sub-group, characterised by a specific immunohistochemistry profile with H3K27me3 loss, wild-type H3K27M and positive EZHIP. It also gives new insights into the possible mechanism and pathway regulation in these tumours, potentially opening new therapeutic avenues for these tumours which have no known effective treatment. This study has been retrospectively registered on clinicaltrial.gov on 8 November 2017 under the registration number NCT03336931 ( https://clinicaltrials.gov/ct2/show/NCT03336931 ).

The causes of pediatric cancers' distinctiveness compared to adult-onset tumors of the same type are not completely clear and not fully explained by their genomes. In this study, we used an optimized multilevel RNA clustering approach to derive molecular definitions for most childhood cancers. Applying this method to 13,313 transcriptomes, we constructed a pediatric cancer atlas to explore age-associated changes. Tumor entities were sometimes unexpectedly grouped due to common lineages, drivers or stemness profiles. Some established entities were divided into subgroups that predicted outcome better than current diagnostic approaches. These definitions account for inter-tumoral and intra-tumoral heterogeneity and have the potential of enabling reproducible, quantifiable diagnostics. As a whole, childhood tumors had more transcriptional diversity than adult tumors, maintaining greater expression flexibility. To apply these insights, we designed an ensemble convolutional neural network classifier. We show that this tool was able to match or clarify the diagnosis for 85% of childhood tumors in a prospective cohort. If further validated, this framework could be extended to derive molecular definitions for all cancer types.

Glutamine is a conditionally essential nutrient for many cancer cells, but it remains unclear how consuming glutamine in excess of growth requirements confers greater fitness to glutamine-addicted cancers. By contrasting two breast cancer subtypes with distinct glutamine dependencies, we show that glutamine-indispensable triple-negative breast cancer (TNBC) cells rely on a non-canonical glutamine-to-glutamate overflow, with glutamine carbon routed once through the TCA cycle. Importantly, this single-pass glutaminolysis increases TCA cycle fluxes and replenishes TCA cycle intermediates in TNBC cells, a process that achieves net oxidation of glucose but not glutamine. The coupling of glucose and glutamine catabolism appears hard-wired via a distinct TNBC gene expression profile biased to strip and then sequester glutamine nitrogen, but hampers the ability of TNBC cells to oxidise glucose when glutamine is limiting. Our results provide a new understanding of how metabolically rigid TNBC cells are sensitive to glutamine deprivation and a way to select vulnerable TNBC subtypes that may be responsive to metabolic-targeted therapies.

RNA-sequencing (RNA-seq) efforts in acute lymphoblastic leukaemia (ALL) have identified numerous prognostically significant genomic alterations which can guide diagnostic risk stratification and treatment choices when detected early. However, integrating RNA-seq in a clinical setting requires rapid detection and accurate reporting of clinically relevant alterations. Here we present RaScALL, an implementation of the k-mer based variant detection tool km, capable of identifying more than 100 prognostically significant lesions observed in ALL, including gene fusions, single nucleotide variants and focal gene deletions. We compared genomic alterations detected by RaScALL and those reported by alignment-based de novo variant detection tools in a study cohort of 180 Australian patient samples. Results were validated using 100 patient samples from a published North American cohort. RaScALL demonstrated a high degree of accuracy for reporting subtype defining genomic alterations. Gene fusions, including difficult to detect fusions involving EPOR and DUX4, were accurately identified in 98% of reported cases in the study cohort (n = 164) and 95% of samples (n = 63) in the validation cohort. Pathogenic sequence variants were correctly identified in 75% of tested samples, including all cases involving subtype defining variants PAX5 p.P80R (n = 12) and IKZF1 p.N159Y (n = 4). Intragenic IKZF1 deletions resulting in aberrant transcript isoforms were also detectable with 98% accuracy. Importantly, the median analysis time for detection of all targeted alterations averaged 22 minutes per sample, significantly shorter than standard alignment-based approaches. The application of RaScALL enables rapid identification and reporting of previously identified genomic alterations of known clinical relevance.

For one-third of patients with pediatric cancer enrolled in precision medicine programs, molecular profiling does not result in a therapeutic recommendation. To identify potential strategies for treating these high-risk pediatric patients, we performed in vitro screening of 125 patient-derived samples against a library of 126 anticancer drugs. Tumor cell expansion did not influence drug responses, and 82% of the screens on expanded tumor cells were completed while the patients were still under clinical care. High-throughput drug screening (HTS) confirmed known associations between activating genomic alterations in NTRK, BRAF, and ALK and responses to matching targeted drugs. The in vitro results were further validated in patient-derived xenograft models in vivo and were consistent with clinical responses in treated patients. In addition, effective combinations could be predicted by correlating sensitivity profiles between drugs. Furthermore, molecular integration with HTS identified biomarkers of sensitivity to WEE1 and MEK inhibition. Incorporating HTS into precision medicine programs is a powerful tool to accelerate the improved identification of effective biomarker-driven therapeutic strategies for treating high-risk pediatric cancers.

The prognosis of recurrent malignant peripheral nerve sheath tumors (MPNST) is dismal, with surgical resection being the only definitive salvage therapy. Treatment with chemoradiation approaches has not significantly improved patient outcomes. Similarly, trials of therapies targeting MPNST genomic drivers have thus far been unsuccessful. Improved understanding of the molecular pathogenesis of MPNST indicates frequent activation of the mitogen-activated protein kinase (MAPK) cell signaling pathway. MEK inhibitors have shown activity in preclinical studies; however, their clinical efficacy has not been reported to date. We describe here a case of sustained complete response to MEK inhibition in an adolescent patient with a recurrent metastatic MPNST with multiple alterations in the MAPK pathway, guided by a precision oncology approach.

Minimal residual disease (MRD) measurement is a cornerstone of contemporary acute lymphoblastic leukaemia (ALL) treatment. The presence of immunoglobulin (Ig) and T cell receptor (TCR) gene recombinations in leukaemic clones allows widespread use of patient-specific, DNA-based MRD assays. In contrast, paediatric solid tumour MRD remains experimental and has focussed on generic assays targeting tumour-specific messenger RNA, methylated DNA or microRNA.

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.

Approximately 50%-55% of high-grade serous ovarian carcinoma (HGSOC) patients have MYC oncogenic pathway activation. Because MYC is not directly targetable, we have analyzed molecular pathways enriched in MYC-high HGSOC tumors to identify potential therapeutic targets. Here, we report that MYC-high HGSOC tumors show enrichment in genes controlled by NRF2, an antioxidant signaling pathway, along with increased thioredoxin redox activity. Treatment of MYC-high HGSOC tumors cells with US Food and Drug Administration (FDA)-approved thioredoxin reductase 1 (TrxR1) inhibitor auranofin resulted in significant growth suppression and apoptosis in MYC-high HGSOC cells in vitro and also significantly reduced tumor growth in an MYC-high HGSOC patient-derived tumor xenograft. We found that auranofin treatment inhibited glycolysis in MYC-high cells via oxidation-induced GAPDH inhibition. Interestingly, in response to auranofin-induced glycolysis inhibition, MYC-high HGSOC cells switched to glutamine metabolism for survival. Depletion of glutamine with either glutamine starvation or glutaminase (GLS1) inhibitor CB-839 exerted synergistic anti-tumor activity with auranofin in HGSOC cells and OVCAR-8 cell line xenograft. These findings suggest that applying a combined therapy of GLS1 inhibitor and TrxR1 inhibitor could effectively treat MYC-high HGSOC patients.

Symptomatic venous thromboembolism (VTE) occurs in five percent of children treated for acute lymphoblastic leukemia (ALL), but whether a genetic predisposition exists across different ALL treatment regimens has not been well studied.

Chromosome 17q21-ter is commonly gained in neuroblastoma, but it is unclear which gene in the region is important for tumorigenesis. The JMJD6 gene at 17q21-ter activates gene transcription. Here we show that JMJD6 forms protein complexes with N-Myc and BRD4, and is important for E2F2, N-Myc and c-Myc transcription. Knocking down JMJD6 reduces neuroblastoma cell proliferation and survival in vitro and tumor progression in mice, and high levels of JMJD6 expression in human neuroblastoma tissues independently predict poor patient prognosis. In addition, JMJD6 gene is associated with transcriptional super-enhancers. Combination therapy with the CDK7/super-enhancer inhibitor THZ1 and the histone deacetylase inhibitor panobinostat synergistically reduces JMJD6, E2F2, N-Myc, c-Myc expression, induces apoptosis in vitro and leads to neuroblastoma tumor regression in mice, which are significantly reversed by forced JMJD6 over-expression. Our findings therefore identify JMJD6 as a neuroblastoma tumorigenesis factor, and the combination therapy as a treatment strategy.

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. Hypothetical model for SNRPD3 as a co-factor for MYCN oncogenesis. SNRPD3 and MYCN participate in a regulatory loop to balance splicing fidelity in neuroblastoma cells. First MYCN transactivates SNRPD3 to lead to high-level expression. Second, SNRPD3 and MYCN form a protein complex involving PRMT5. Third, this leads to balanced alterative splicing (AS) activitiy that is favorable to neuroblastoma. Together this forms as a therapeutic vulnerability where SNRPD3 perturbation or PRMT5 inhibitors are selectively toxic to neuroblastoma by conditionally disturbing splicing activity.

Ependymomas (EPN) are the third most common malignant brain cancer in children. Treatment strategies for pediatric EPN have remained unchanged over recent decades, with 10-year survival rates stagnating at just 67% for children aged 0-14 years. Moreover, a proportion of patients who survive treatment often suffer long-term neurological side effects as a result of therapy. It is evident that there is a need for safer, more effective treatments for pediatric EPN patients. There are ten distinct subgroups of EPN, each with their own molecular and prognostic features. To identify and facilitate the testing of new treatments for EPN, in vivo laboratory models representative of the diverse molecular subtypes are required. Here, we describe the establishment of a patient-derived orthotopic xenograft (PDOX) model of posterior fossa A (PFA) EPN, derived from a metastatic cranial lesion.