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.

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.

Purpose: Treating refractory or relapsed neuroblastoma remains challenging. Monitoring body fluids for tumor-derived molecular information indicating minimal residual disease supports more frequent diagnostic surveillance and may have the power to detect resistant subclones before they give rise to relapses. If actionable targets are identified from liquid biopsies, targeted treatment options can be considered earlier. Experimental Design: Droplet digital PCR assays assessing MYCN and ALK copy numbers and allelic frequencies of ALK p.F1174L and ALK p.R1275Q mutations were applied to longitudinally collected liquid biopsies and matched tumor tissue samples from 31 patients with high-risk neuroblastoma. Total cell-free DNA (cfDNA) levels and marker detection were compared with data from routine clinical diagnostics. Results: Total cfDNA concentrations in blood plasma from patients with high-risk neuroblastoma were higher than in healthy controls and consistently correlated with neuron-specific enolase levels and lactate dehydrogenase activity but not with 123I-meta-iodobenzylguanidine scores at relapse diagnosis. Targeted cfDNA diagnostics proved superior for early relapse detection to all current diagnostics in 2 patients. Marker analysis in cfDNA indicated intratumor heterogeneity for cell clones harboring MYCN amplifications and druggable ALK alterations that were not detectable in matched tumor tissue samples in 17 patients from our cohort. Proof of concept is provided for molecular target detection in cerebrospinal fluid from patients with isolated central nervous system relapses. Conclusions: Tumor-specific alterations can be identified and monitored during disease course in liquid biopsies from pediatric patients with high-risk neuroblastoma. This approach to cfDNA surveillance warrants further prospective validation and exploitation for diagnostic purposes and to guide therapeutic decisions.

Targeting tachykinin receptors in neuroblastoma Targeting tachykinin receptors in neuroblastoma

Circular extrachromosomal DNA (ecDNA) is a form of oncogene amplification found across cancer types and associated with poor outcome in patients. ecDNA can be structurally complex and can contain rearranged DNA sequences derived from multiple chromosome locations. As the structure of ecDNA can impact oncogene regulation and may indicate mechanisms of its formation, disentangling it at high resolution from sequencing data is essential. Even though methods have been developed to identify and reconstruct ecDNA in cancer genome sequencing, it remains challenging to resolve complex ecDNA structures, in particular amplicons with shared genomic footprints. We here introduce Decoil, a computational method that combines a breakpoint-graph approach withLASSOregression to reconstruct complex ecDNA and deconvolve co-occurring ecDNA elements with overlapping genomic footprints from long-read nanopore sequencing. Decoil outperforms de novo assembly and alignment-based methods in simulated long-read sequencing data for both simple and complex ecDNAs. Applying Decoil on whole-genome sequencing data uncovered different ecDNA topologies and explored ecDNA structure heterogeneity in neuroblastoma tumors and cell lines, indicating that this method may improve ecDNA structural analyses in cancer.

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.