Reconstructing extrachromosomal DNA structural heterogeneity from long-read sequencing data using Decoil.

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 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 which combines a breakpoint-graph approach with regression 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 analyzes in cancer.

Subclonal NT5C2 mutations are associated with poor outcomes after relapse of pediatric acute lymphoblastic leukemia.

Activating mutations in cytosolic 5'-nucleotidase II (NT5C2) are considered to drive relapse formation in acute lymphoblastic leukemia (ALL) by conferring purine analog resistance. To examine the clinical effects of NT5C2 mutations in relapsed ALL, we analyzed NT5C2 in 455 relapsed B-cell precursor ALL patients treated within the ALL-REZ BFM 2002 relapse trial using sequencing and sensitive allele-specific real-time polymerase chain reaction. We detected 110 NT5C2 mutations in 75 (16.5%) of 455 B-cell precursor ALL relapses. Two-thirds of relapses harbored subclonal mutations and only one-third harbored clonal mutations. Event-free survival after relapse was inferior in patients with relapses with clonal and subclonal NT5C2 mutations compared with those without (19% and 25% vs 53%, P < .001). However, subclonal, but not clonal, NT5C2 mutations were associated with reduced event-free survival in multivariable analysis (hazard ratio, 1.89; 95% confidence interval, 1.28-2.69; P = .001) and with an increased rate of nonresponse to relapse treatment (subclonal 32%, clonal 12%, wild type 9%, P < .001). Nevertheless, 27 (82%) of 33 subclonal NT5C2 mutations became undetectable at the time of nonresponse or second relapse, and in 10 (71%) of 14 patients subclonal NT5C2 mutations were undetectable already after relapse induction treatment. These results show that subclonal NT5C2 mutations define relapses associated with high risk of treatment failure in patients and at the same time emphasize that their role in outcome is complex and goes beyond mutant NT5C2 acting as a targetable driver during relapse progression. Sensitive, prospective identification of NT5C2 mutations is warranted to improve the understanding and treatment of this aggressive ALL relapse subtype.

High sensitivity and clonal stability of the genomic fusion as single marker for response monitoring in ETV6-RUNX1-positive acute lymphoblastic leukemia.

BACKGROUND: Assessment of minimal residual disease (MRD) is an integral component for response monitoring and treatment stratification in acute lymphoblastic leukemia (ALL). We aimed to evaluate the genomic ETV6-RUNX1 fusion sites as a single marker for MRD quantification. PROCEDURE: In a representative, uniformly treated cohort of pediatric relapsed ALL patients (n = 52), ETV6-RUNX1 fusion sites were compared to the current gold standard, immunoglobulin/T-cell receptor (Ig/TCR) gene rearrangements. RESULTS: Primer/probe sets designed to ETV6-RUNX1 fusions achieved significantly more frequent a sensitivity and a quantitative range of at least 10-4 compared to the gold standard with 100% and 73% versus 76% and 47%, respectively. The breakpoint sequence was identical at diagnosis and relapse in all tested cases. There was a high degree of concordance between quantitative MRD results assessed using ETV6-RUNX1 and the highest Ig/TCR marker (Spearman's 0.899, P < .01) with differences >½ log-step in only 6% of patients. A high proportion of ETV6-RUNX1-positive ALL relapses (40%) in our cohort showed a poor response to induction treatment at relapse, and therefore had an indication for hematopoietic stem cell transplantation, demonstrating the need of accurate identification of this subgroup. CONCLUSIONS: ETV6-RUNX1 fusion sites are highly sensitive and reliable MRD markers. Our data confirm that they are unaffected by clonal evolution and selection during front-line and second-line chemotherapy in contrast to Ig/TCR rearrangements, which require several markers per patient to compensate for the observed loss of target clones. In future studies, the genomic ETV6-RUNX1 fusion can be used as single MRD marker.

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.

Targeted Analysis of Cell-free Circulating Tumor DNA is Suitable for Early Relapse and Actionable Target Detection in Patients with Neuroblastoma.

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.

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.

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.

Small-Molecule Dual PLK1 and BRD4 Inhibitors are Active Against Preclinical Models of Pediatric Solid Tumors.

Simultaneous inhibition of multiple molecular targets is an established strategy to improve the continuance of clinical response to therapy. Here, we screened 49 molecules with dual nanomolar inhibitory activity against BRD4 and PLK1, best classified as dual kinase-bromodomain inhibitors, in pediatric tumor cell lines for their antitumor activity. We identified two candidate dual kinase-bromodomain inhibitors with strong and tumor-specific activity against neuroblastoma, medulloblastoma, and rhabdomyosarcoma tumor cells. Dual PLK1 and BRD4 inhibitor treatment suppressed proliferation and induced apoptosis in pediatric tumor cell lines at low nanomolar concentrations. This was associated with reduced MYCN-driven gene expression as assessed by RNA sequencing. Treatment of patient-derived xenografts with dual inhibitor UMB103 led to significant tumor regression. We demonstrate that concurrent inhibition of two central regulators of MYC protein family of protooncogenes, BRD4, and PLK1, with single small molecules has strong and specific antitumor effects in preclinical pediatric cancer models.

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.

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.

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.

Targeting tachykinin receptors in neuroblastoma.

Neuroblastoma is the most common extracranial tumor in children. Despite aggressive multimodal treatment, high-risk neuroblastoma remains a clinical challenge with survival rates below 50%. Adding targeted drugs to first-line therapy regimens is a promising approach to improve survival in these patients. TACR1 activation by substance P has been reported to be mitogenic in cancer cell lines. Tachykinin receptor (TACR1) antagonists are approved for clinical use as an antiemetic remedy since 2003. Tachykinin receptor inhibition has recently been shown to effectively reduce growth of several tumor types. Here, we report that neuroblastoma cell lines express TACR1, and that targeting TACR1 activity significantly reduced cell viability and induced apoptosis in neuroblastoma cell lines. Gene expression profiling revealed that TACR1 inhibition repressed E2F2 and induced TP53 signaling. Treating mice harboring established neuroblastoma xenograft tumors with Aprepitant also significantly reduced tumor burden. Thus, we provide evidence that the targeted inhibition of tachykinin receptor signaling shows therapeutic efficacy in preclinical models for high-risk neuroblastoma.