Mutational topography reflects clinical neuroblastoma heterogeneity.

Neuroblastoma is a pediatric solid tumor characterized by strong clinical heterogeneity. Although clinical risk-defining genomic alterations exist in neuroblastomas, the mutational processes involved in their generation remain largely unclear. By examining the topography and mutational signatures derived from all variant classes, we identified co-occurring mutational footprints, which we termed mutational scenarios. We demonstrate that clinical neuroblastoma heterogeneity is associated with differences in the mutational processes driving these scenarios, linking risk-defining pathognomonic variants to distinct molecular processes. Whereas high-risk MYCN-amplified neuroblastomas were characterized by signs of replication slippage and stress, homologous recombination-associated signatures defined high-risk non-MYCN-amplified patients. Non-high-risk neuroblastomas were marked by footprints of chromosome mis-segregation and TOP1 mutational activity. Furthermore, analysis of subclonal mutations uncovered differential activity of these processes through neuroblastoma evolution. Thus, clinical heterogeneity of neuroblastoma patients can be linked to differences in the mutational processes that are active in their tumors.

Defining the landscape of circular RNAs in neuroblastoma unveils a global suppressive function of MYCN.

Circular RNAs (circRNAs) are a regulatory RNA class. While cancer-driving functions have been identified for single circRNAs, how they modulate gene expression in cancer is not well understood. We investigate circRNA expression in the pediatric malignancy, neuroblastoma, through deep whole-transcriptome sequencing in 104 primary neuroblastomas covering all risk groups. We demonstrate that MYCN amplification, which defines a subset of high-risk cases, causes globally suppressed circRNA biogenesis directly dependent on the DHX9 RNA helicase. We detect similar mechanisms in shaping circRNA expression in the pediatric cancer medulloblastoma implying a general MYCN effect. Comparisons to other cancers identify 25 circRNAs that are specifically upregulated in neuroblastoma, including circARID1A. Transcribed from the ARID1A tumor suppressor gene, circARID1A promotes cell growth and survival, mediated by direct interaction with the KHSRP RNA-binding protein. Our study highlights the importance of MYCN regulating circRNAs in cancer and identifies molecular mechanisms, which explain their contribution to neuroblastoma pathogenesis.

Mutations in ALK signaling pathways conferring resistance to ALK inhibitor treatment lead to collateral vulnerabilities in neuroblastoma cells.

BACKGROUND: Development of resistance to targeted therapies has tempered initial optimism that precision oncology would improve poor outcomes for cancer patients. Resistance mechanisms, however, can also confer new resistance-specific vulnerabilities, termed collateral sensitivities. Here we investigated anaplastic lymphoma kinase (ALK) inhibitor resistance in neuroblastoma, a childhood cancer frequently affected by activating ALK alterations. METHODS: Genome-wide forward genetic CRISPR-Cas9 based screens were performed to identify genes associated with ALK inhibitor resistance in neuroblastoma cell lines. Furthermore, the neuroblastoma cell line NBLW-R was rendered resistant by continuous exposure to ALK inhibitors. Genes identified to be associated with ALK inhibitor resistance were further investigated by generating suitable cell line models. In addition, tumor and liquid biopsy samples of four patients with ALK-mutated neuroblastomas before ALK inhibitor treatment and during tumor progression under treatment were genomically profiled. RESULTS: Both genome-wide CRISPR-Cas9-based screens and preclinical spontaneous ALKi resistance models identified NF1 loss and activating NRASQ61K mutations to confer resistance to chemically diverse ALKi. Moreover, human neuroblastomas recurrently developed de novo loss of NF1 and activating RAS mutations after ALKi treatment, leading to therapy resistance. Pathway-specific perturbations confirmed that NF1 loss and activating RAS mutations lead to RAS-MAPK signaling even in the presence of ALKi. Intriguingly, NF1 loss rendered neuroblastoma cells hypersensitive to MEK inhibition. CONCLUSIONS: Our results provide a clinically relevant mechanistic model of ALKi resistance in neuroblastoma and highlight new clinically actionable collateral sensitivities in resistant cells.

Neuroblastoma signalling models unveil combination therapies targeting feedback-mediated resistance.

Very high risk neuroblastoma is characterised by increased MAPK signalling, and targeting MAPK signalling is a promising therapeutic strategy. We used a deeply characterised panel of neuroblastoma cell lines and found that the sensitivity to MEK inhibitors varied drastically between these cell lines. By generating quantitative perturbation data and mathematical modelling, we determined potential resistance mechanisms. We found that negative feedbacks within MAPK signalling and via the IGF receptor mediate re-activation of MAPK signalling upon treatment in resistant cell lines. By using cell-line specific models, we predict that combinations of MEK inhibitors with RAF or IGFR inhibitors can overcome resistance, and tested these predictions experimentally. In addition, phospho-proteomic profiling confirmed the cell-specific feedback effects and synergy of MEK and IGFR targeted treatment. Our study shows that a quantitative understanding of signalling and feedback mechanisms facilitated by models can help to develop and optimise therapeutic strategies. Our findings should be considered for the planning of future clinical trials introducing MEKi in the treatment of neuroblastoma.

Chromothripsis followed by circular recombination drives oncogene amplification in human cancer.

The mechanisms behind the evolution of complex genomic amplifications in cancer have remained largely unclear. Using whole-genome sequencing data of the pediatric tumor neuroblastoma, we here identified a type of amplification, termed 'seismic amplification', that is characterized by multiple rearrangements and discontinuous copy number levels. Overall, seismic amplifications occurred in 9.9% (274 of 2,756) of cases across 38 cancer types, and were associated with massively increased copy numbers and elevated oncogene expression. Reconstruction of the development of seismic amplification showed a stepwise evolution, starting with a chromothripsis event, followed by formation of circular extrachromosomal DNA that subsequently underwent repetitive rounds of circular recombination. The resulting amplicons persisted as extrachromosomal DNA circles or had reintegrated into the genome in overt tumors. Together, our data indicate that the sequential occurrence of chromothripsis and circular recombination drives oncogene amplification and overexpression in a substantial fraction of human malignancies.

Spatial and temporal intratumour heterogeneity has potential consequences for single biopsy-based neuroblastoma treatment decisions.

Intratumour heterogeneity is a major cause of treatment failure in cancer. We present in-depth analyses combining transcriptomic and genomic profiling with ultra-deep targeted sequencing of multiregional biopsies in 10 patients with neuroblastoma, a devastating childhood tumour. We observe high spatial and temporal heterogeneity in somatic mutations and somatic copy-number alterations which are reflected on the transcriptomic level. Mutations in some druggable target genes including ALK and FGFR1 are heterogeneous at diagnosis and/or relapse, raising the issue whether current target prioritization and molecular risk stratification procedures in single biopsies are sufficiently reliable for therapy decisions. The genetic heterogeneity in gene mutations and chromosome aberrations observed in deep analyses from patient courses suggest clonal evolution before treatment and under treatment pressure, and support early emergence of metastatic clones and ongoing chromosomal instability during disease evolution. We report continuous clonal evolution on mutational and copy number levels in neuroblastoma, and detail its implications for therapy selection, risk stratification and therapy resistance.

Neuroblastoma Risk Assessment and Treatment Stratification with Hybrid Capture-Based Panel Sequencing.

For many years, the risk-based therapy stratification of children with neuroblastoma has relied on clinical and molecular covariates. In recent years, genome analysis has revealed further alterations defining risk, tumor biology, and therapeutic targets. The implementation of a robust and scalable method for analyzing traditional and new molecular markers in routine diagnostics is an urgent clinical need. Here, we investigated targeted panel sequencing as a diagnostic approach to analyze all relevant genomic neuroblastoma risk markers in one assay. Our "neuroblastoma hybrid capture sequencing panel" (NB-HCSP) assay employs a technology for the high-coverage sequencing (>1000×) of 55 selected genes and neuroblastoma-relevant genomic regions, which allows for the detection of single nucleotide changes, structural rearrangements, and copy number alterations. We validated our assay by analyzing 15 neuroblastoma cell lines and a cohort of 20 neuroblastomas, for which reference routine diagnostic data and genome sequencing data were available. We observed a high concordance for risk markers identified by the NB-HSCP assay, clinical routine diagnostics, and genome sequencing. Subsequently, we demonstrated clinical applicability of the NB-HCSP assay by analyzing routine clinical samples. We conclude that the NB-HCSP assay may be implemented into routine diagnostics as a single assay that covers all essential covariates for initial neuroblastoma classification, extended risk stratification, and targeted therapy selection.

Enhancer hijacking determines extrachromosomal circular MYCN amplicon architecture in neuroblastoma.

MYCN amplification drives one in six cases of neuroblastoma. The supernumerary gene copies are commonly found on highly rearranged, extrachromosomal circular DNA (ecDNA). The exact amplicon structure has not been described thus far and the functional relevance of its rearrangements is unknown. Here, we analyze the MYCN amplicon structure using short-read and Nanopore sequencing and its chromatin landscape using ChIP-seq, ATAC-seq and Hi-C. This reveals two distinct classes of amplicons which explain the regulatory requirements for MYCN overexpression. The first class always co-amplifies a proximal enhancer driven by the noradrenergic core regulatory circuit (CRC). The second class of MYCN amplicons is characterized by high structural complexity, lacks key local enhancers, and instead contains distal chromosomal fragments harboring CRC-driven enhancers. Thus, ectopic enhancer hijacking can compensate for the loss of local gene regulatory elements and explains a large component of the structural diversity observed in MYCN amplification.

Publisher Correction: Extrachromosomal circular DNA drives oncogenic genome remodeling in neuroblastoma.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Extrachromosomal circular DNA drives oncogenic genome remodeling in neuroblastoma.

Extrachromosomal circularization of DNA is an important genomic feature in cancer. However, the structure, composition and genome-wide frequency of extrachromosomal circular DNA have not yet been profiled extensively. Here, we combine genomic and transcriptomic approaches to describe the landscape of extrachromosomal circular DNA in neuroblastoma, a tumor arising in childhood from primitive cells of the sympathetic nervous system. Our analysis identifies and characterizes a wide catalog of somatically acquired and undescribed extrachromosomal circular DNAs. Moreover, we find that extrachromosomal circular DNAs are an unanticipated major source of somatic rearrangements, contributing to oncogenic remodeling through chimeric circularization and reintegration of circular DNA into the linear genome. Cancer-causing lesions can emerge out of circle-derived rearrangements and are associated with adverse clinical outcome. It is highly probable that circle-derived rearrangements represent an ongoing mutagenic process. Thus, extrachromosomal circular DNAs represent a multihit mutagenic process, with important functional and clinical implications for the origins of genomic remodeling in cancer.

UBQLN4 Represses Homologous Recombination and Is Overexpressed in Aggressive Tumors.

Genomic instability can be a hallmark of both human genetic disease and cancer. We identify a deleterious UBQLN4 mutation in families with an autosomal recessive syndrome reminiscent of genome instability disorders. UBQLN4 deficiency leads to increased sensitivity to genotoxic stress and delayed DNA double-strand break (DSB) repair. The proteasomal shuttle factor UBQLN4 is phosphorylated by ATM and interacts with ubiquitylated MRE11 to mediate early steps of homologous recombination-mediated DSB repair (HRR). Loss of UBQLN4 leads to chromatin retention of MRE11, promoting non-physiological HRR activity in vitro and in vivo. Conversely, UBQLN4 overexpression represses HRR and favors non-homologous end joining. Moreover, we find UBQLN4 overexpressed in aggressive tumors. In line with an HRR defect in these tumors, UBQLN4 overexpression is associated with PARP1 inhibitor sensitivity. UBQLN4 therefore curtails HRR activity through removal of MRE11 from damaged chromatin and thus offers a therapeutic window for PARP1 inhibitor treatment in UBQLN4-overexpressing tumors.

Telomerase Is a Prognostic Marker of Poor Outcome and a Therapeutic Target in Neuroblastoma.

PURPOSE: Telomere maintenance is a hallmark of high-risk neuroblastoma; however, the contribution of telomerase and alternative lengthening of telomeres (ALT) to clinical phenotypes has remained unclear. We aimed to determine the clinical relevance of telomerase activation versus ALT as biomarkers in pretreatment neuroblastoma and to assess the potential value of telomerase as a therapeutic target. MATERIALS AND METHODS: The genomic status of TERT and MYCN was assessed in 457 pretreatment neuroblastomas by fluorescence in situ hybridization. ALT was examined in 273 of 457 tumors by detection of ALT-associated promyelocytic leukemia nuclear bodies, and TERT expression was determined by microarrays in 223 of these. Cytotoxic effects of telomerase-interacting compounds were analyzed in neuroblastoma cell lines in vitro and in vivo. RESULTS: We detected TERT rearrangements in 46 of 457 cases (10.1%), MYCN amplification in 93 of 457 cases (20.4%), and elevated TERT expression in tumors lacking TERT or MYCN alterations in 10 of 223 cases (4.5%). ALT activation was found in 49 of 273 cases (17.9%). All these alterations occurred almost mutually exclusively and were associated with unfavorable prognostic variables and adverse outcome. The presence of activated telomerase (ie, TERT rearrangements, MYCN amplification, or high TERT expression without these alterations) was associated with poorest overall survival and was an independent prognostic marker in multivariable analyses. We also found that the telomerase-interacting compound 6-thio-2'-deoxyguanosine effectively inhibited viability and proliferation of neuroblastoma cells bearing activated telomerase. Similarly, tumor growth was strongly impaired upon 6-thio-2'-deoxyguanosine treatment in telomerase-positive neuroblastoma xenografts in mice. CONCLUSION: Our data suggest telomerase activation and ALT define distinct neuroblastoma subgroups with adverse outcome and that telomerase may represent a promising therapeutic target in many high-risk neuroblastomas.

A mechanistic classification of clinical phenotypes in neuroblastoma.

Neuroblastoma is a pediatric tumor of the sympathetic nervous system. Its clinical course ranges from spontaneous tumor regression to fatal progression. To investigate the molecular features of the divergent tumor subtypes, we performed genome sequencing on 416 pretreatment neuroblastomas and assessed telomere maintenance mechanisms in 208 of these tumors. We found that patients whose tumors lacked telomere maintenance mechanisms had an excellent prognosis, whereas the prognosis of patients whose tumors harbored telomere maintenance mechanisms was substantially worse. Survival rates were lowest for neuroblastoma patients whose tumors harbored telomere maintenance mechanisms in combination with RAS and/or p53 pathway mutations. Spontaneous tumor regression occurred both in the presence and absence of these mutations in patients with telomere maintenance-negative tumors. On the basis of these data, we propose a mechanistic classification of neuroblastoma that may benefit the clinical management of patients.

The mutational landscape of MYCN, Lin28b and ALKF1174L driven murine neuroblastoma mimics human disease.

Genetically engineered mouse models have proven to be essential tools for unraveling fundamental aspects of cancer biology and for testing novel therapeutic strategies. To optimally serve these goals, it is essential that the mouse model faithfully recapitulates the human disease. Recently, novel mouse models for neuroblastoma have been developed. Here, we report on the further genomic characterization through exome sequencing and DNA copy number analysis of four of the currently available murine neuroblastoma model systems (ALK, Th-MYCN, Dbh-MYCN and Lin28b). The murine tumors revealed a low number of genomic alterations - in keeping with human neuroblastoma - and a positive correlation of the number of genetic lesions with the time to onset of tumor formation was observed. Gene copy number alterations are the hallmark of both murine and human disease and frequently affect syntenic genomic regions. Despite low mutational load, the genes mutated in murine disease were found to be enriched for genes mutated in human disease. Taken together, our study further supports the validity of the tested mouse models for mechanistic and preclinical studies of human neuroblastoma.

Sense-Antisense lncRNA Pair Encoded by Locus 6p22.3 Determines Neuroblastoma Susceptibility via the USP36-CHD7-SOX9 Regulatory Axis.

Trait-associated loci often map to genomic regions encoding long noncoding RNAs (lncRNAs), but the role of these lncRNAs in disease etiology is largely unexplored. We show that a pair of sense/antisense lncRNA (6p22lncRNAs) encoded by CASC15 and NBAT1 located at the neuroblastoma (NB) risk-associated 6p22.3 locus are tumor suppressors and show reduced expression in high-risk NBs. Loss of functional synergy between 6p22lncRNAs results in an undifferentiated state that is maintained by a gene-regulatory network, including SOX9 located on 17q, a region frequently gained in NB. 6p22lncRNAs regulate SOX9 expression by controlling CHD7 stability via modulating the cellular localization of USP36, encoded by another 17q gene. This regulatory nexus between 6p22.3 and 17q regions may lead to potential NB treatment strategies.

The putative tumor suppressor gene EphA7 is a novel BMI-1 target.

Bmi1 was originally identified as a gene that contributes to the development of mouse lymphoma by inhibiting MYC-induced apoptosis through repression of Ink4a and Arf. It codes for the Polycomb group protein BMI-1 and acts primarily as a transcriptional repressor via chromatin modifications. Although it binds to a large number of genomic regions, the direct BMI-1 target genes described so far do not explain the full spectrum of BMI-1-mediated effects. Here we identify the putative tumor suppressor gene EphA7 as a novel direct BMI-1 target in neural cells and lymphocytes. EphA7 silencing has been reported in several different human tumor types including lymphomas, and our data suggest BMI1 overexpression as a novel mechanism leading to EphA7 inactivation via H3K27 trimethylation and DNA methylation.

Regulation of Nuclear Hormone Receptors by MYCN-Driven miRNAs Impacts Neural Differentiation and Survival in Neuroblastoma Patients.

MYCN amplification and MYC signaling are associated with high-risk neuroblastoma with poor prognosis. Treating these tumors remains challenging, although therapeutic approaches stimulating differentiation have generated considerable interest. We have previously shown that the MYCN-regulated miR-17∼92 cluster inhibits neuroblastoma differentiation by repressing estrogen receptor alpha. Here, we demonstrate that this microRNA (miRNA) cluster selectively targets several members of the nuclear hormone receptor (NHR) superfamily, and we present a unique NHR signature associated with the survival of neuroblastoma patients. We found that suppressing glucocorticoid receptor (GR) expression in MYCN-driven patient and mouse tumors was associated with an undifferentiated phenotype and decreased survival. Importantly, MYCN inhibition and subsequent reactivation of GR signaling promotes neural differentiation and reduces tumor burden. Our findings reveal a key role for the miR-17∼92-regulated NHRs in neuroblastoma biology, thereby providing a potential differentiation approach for treating neuroblastoma patients.

Transcription factor activating protein 2 beta (TFAP2B) mediates noradrenergic neuronal differentiation in neuroblastoma.

Neuroblastoma is an embryonal pediatric tumor that originates from the developing sympathetic nervous system and shows a broad range of clinical behavior, ranging from fatal progression to differentiation into benign ganglioneuroma. In experimental neuroblastoma systems, retinoic acid (RA) effectively induces neuronal differentiation, and RA treatment has been therefore integrated in current therapies. However, the molecular mechanisms underlying differentiation are still poorly understood. We here investigated the role of transcription factor activating protein 2 beta (TFAP2B), a key factor in sympathetic nervous system development, in neuroblastoma pathogenesis and differentiation. Microarray analyses of primary neuroblastomas (n = 649) demonstrated that low TFAP2B expression was significantly associated with unfavorable prognostic markers as well as adverse patient outcome. We also found that low TFAP2B expression was strongly associated with CpG methylation of the TFAP2B locus in primary neuroblastomas (n = 105) and demethylation with 5-aza-2'-deoxycytidine resulted in induction of TFAP2B expression in vitro, suggesting that TFAP2B is silenced by genomic methylation. Tetracycline inducible re-expression of TFAP2B in IMR-32 and SH-EP neuroblastoma cells significantly impaired proliferation and cell cycle progression. In IMR-32 cells, TFAP2B induced neuronal differentiation, which was accompanied by up-regulation of the catecholamine biosynthesizing enzyme genes DBH and TH, and down-regulation of MYCN and REST, a master repressor of neuronal genes. By contrast, knockdown of TFAP2B by lentiviral transduction of shRNAs abrogated RA-induced neuronal differentiation of SH-SY5Y and SK-N-BE(2)c neuroblastoma cells almost completely. Taken together, our results suggest that TFAP2B is playing a vital role in retaining RA responsiveness and mediating noradrenergic neuronal differentiation in neuroblastoma.

Boolean modeling identifies Greatwall/MASTL as an important regulator in the AURKA network of neuroblastoma.

Aurora Kinase A (AURKA) is often overexpressed in neuroblastoma (NB) with poor outcome. The causes of AURKA overexpression in NB are unknown. Here, we describe a gene regulatory network consisting of core regulators of AURKA protein expression and activation during mitosis to identify potential causes. This network was transformed to a dynamic Boolean model. Simulated activation of the serine/threonine protein kinase Greatwall (GWL, encoded by MASTL) that attenuates the pivotal AURKA inhibitor PP2A, predicted stabilization of AURKA. Consistent with this notion, gene set enrichment analysis showed enrichment of mitotic spindle assembly genes and MYCN target genes in NB with high GWL/MASTL expression. In line with the prediction of GWL/MASTL enhancing AURKA, elevated expression of GWL/MASTL was associated with NB risk factors and poor survival of patients. These results establish Boolean network modeling of oncogenic pathways in NB as a useful means for guided discovery in this enigmatic cancer.