Search

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., (© 2023. The Author(s).)

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 NRAS Q61K 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., (© 2022. The Author(s).)

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., (© 2021. The Author(s).)

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., Competing Interests: The authors have declared that no competing interests exist.

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 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.

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., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

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., (Copyright © 2016 The Author(s). Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2015 Elsevier Ireland Ltd. All rights reserved.)

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., (Copyright © 2015 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.)

Neuroblastoma is a malignant paediatric tumour of the sympathetic nervous system. Roughly half of these tumours regress spontaneously or are cured by limited therapy. By contrast, high-risk neuroblastomas have an unfavourable clinical course despite intensive multimodal treatment, and their molecular basis has remained largely elusive. Here we have performed whole-genome sequencing of 56 neuroblastomas (high-risk, n = 39; low-risk, n = 17) and discovered recurrent genomic rearrangements affecting a chromosomal region at 5p15.33 proximal of the telomerase reverse transcriptase gene (TERT). These rearrangements occurred only in high-risk neuroblastomas (12/39, 31%) in a mutually exclusive fashion with MYCN amplifications and ATRX mutations, which are known genetic events in this tumour type. In an extended case series (n = 217), TERT rearrangements defined a subgroup of high-risk tumours with particularly poor outcome. Despite a large structural diversity of these rearrangements, they all induced massive transcriptional upregulation of TERT. In the remaining high-risk tumours, TERT expression was also elevated in MYCN-amplified tumours, whereas alternative lengthening of telomeres was present in neuroblastomas without TERT or MYCN alterations, suggesting that telomere lengthening represents a central mechanism defining this subtype. The 5p15.33 rearrangements juxtapose the TERT coding sequence to strong enhancer elements, resulting in massive chromatin remodelling and DNA methylation of the affected region. Supporting a functional role of TERT, neuroblastoma cell lines bearing rearrangements or amplified MYCN exhibited both upregulated TERT expression and enzymatic telomerase activity. In summary, our findings show that remodelling of the genomic context abrogates transcriptional silencing of TERT in high-risk neuroblastoma and places telomerase activation in the centre of transformation in a large fraction of these tumours.

Background: Gene expression profiling is being widely applied in cancer research to identify biomarkers for clinical endpoint prediction. Since RNA-seq provides a powerful tool for transcriptome-based applications beyond the limitations of microarrays, we sought to systematically evaluate the performance of RNA-seq-based and microarray-based classifiers in this MAQC-III/SEQC study for clinical endpoint prediction using neuroblastoma as a model., Results: We generate gene expression profiles from 498 primary neuroblastomas using both RNA-seq and 44 k microarrays. Characterization of the neuroblastoma transcriptome by RNA-seq reveals that more than 48,000 genes and 200,000 transcripts are being expressed in this malignancy. We also find that RNA-seq provides much more detailed information on specific transcript expression patterns in clinico-genetic neuroblastoma subgroups than microarrays. To systematically compare the power of RNA-seq and microarray-based models in predicting clinical endpoints, we divide the cohort randomly into training and validation sets and develop 360 predictive models on six clinical endpoints of varying predictability. Evaluation of factors potentially affecting model performances reveals that prediction accuracies are most strongly influenced by the nature of the clinical endpoint, whereas technological platforms (RNA-seq vs. microarrays), RNA-seq data analysis pipelines, and feature levels (gene vs. transcript vs. exon-junction level) do not significantly affect performances of the models., Conclusions: We demonstrate that RNA-seq outperforms microarrays in determining the transcriptomic characteristics of cancer, while RNA-seq and microarray-based models perform similarly in clinical endpoint prediction. Our findings may be valuable to guide future studies on the development of gene expression-based predictive models and their implementation in clinical practice.

Neuroblastoma is an embryonal tumor of the sympathetic nervous system and the most common extracranial tumor of childhood. By sequencing transcriptomes of low- and high-risk neuroblastomas, we detected differentially expressed annotated and nonannotated long noncoding RNAs (lncRNAs). We identified a lncRNA neuroblastoma associated transcript-1 (NBAT-1) as a biomarker significantly predicting clinical outcome of neuroblastoma. CpG methylation and a high-risk neuroblastoma associated SNP on chromosome 6p22 functionally contribute to NBAT-1 differential expression. Loss of NBAT-1 increases cellular proliferation and invasion. It controls these processes via epigenetic silencing of target genes. NBAT-1 loss affects neuronal differentiation through activation of the neuronal-specific transcription factor NRSF/REST. Thus, loss of NBAT-1 contributes to aggressive neuroblastoma by increasing proliferation and impairing differentiation of neuronal precursors., (Copyright © 2014 Elsevier Inc. All rights reserved.)

The prognostic relevance of chromosome 17 gain in neuroblastoma is still discussed. This investigation specifies the frequency, type, size, and transcriptional relevance in a large patient cohort. Primary tumor material of 202 patients was analyzed using high-resolution oligonucleotide array-based comparative genomic hybridization (aCGH) and correlated with clinical and survival data. A subset (n = 145) was correlated for differentially expressed genes (DEG) by microarray analysis. Chromosome 17 aCGH analysis showed numerical gain in 94/202 patients (47%), partial gain in 93/202 patients (46%), and no gain in 15/202 patients (7%). The frequency of partial gain was higher in stage 4 neuroblastoma (stage 1 15%; stage 2 12%; stage 3 16%; stage 4S 7%; and stage 4 50%). Overall survival (OS) was superior in patients with numerical gain compared with patients with partial gain or no gain (5-y-OS: 0.95 ± 0.02 vs. 0.63 ± 0.05 vs. 0.60 ± 0.13; P < 0.001). Gene expression analysis demonstrated 95/130 DEGs between tumors with numerical or partial chromosome/no gain. Only one DEG (CCKBR) was detected comparing tumors with partial gain and those with no gain. In patients with partial gain, the distribution of breakpoints did not correlate with stage and 11q status, but with MYCN amplification and 1p status. The "best" breakpoints in cases with partial 17q gain were at 42.5 Mb for event-free and 26.6 Mb for OS. Numerical gain of chromosome 17 is associated with a better prognosis than partial and no gain. The group of tumors with partial gain was similar to the group without gain with respect to stage distribution, outcome, and gene expression profile., (© 2014 Wiley Periodicals, Inc.)