Sequencing of neuroblastoma identifies chromothripsis and defects in neuritogenesis genes.

Neuroblastoma is a childhood tumour of the peripheral sympathetic nervous system. The pathogenesis has for a long time been quite enigmatic, as only very few gene defects were identified in this often lethal tumour. Frequently detected gene alterations are limited to MYCN amplification (20%) and ALK activations (7%). Here we present a whole-genome sequence analysis of 87 neuroblastoma of all stages. Few recurrent amino-acid-changing mutations were found. In contrast, analysis of structural defects identified a local shredding of chromosomes, known as chromothripsis, in 18% of high-stage neuroblastoma. These tumours are associated with a poor outcome. Structural alterations recurrently affected ODZ3, PTPRD and CSMD1, which are involved in neuronal growth cone stabilization. In addition, ATRX, TIAM1 and a series of regulators of the Rac/Rho pathway were mutated, further implicating defects in neuritogenesis in neuroblastoma. Most tumours with defects in these genes were aggressive high-stage neuroblastomas, but did not carry MYCN amplifications. The genomic landscape of neuroblastoma therefore reveals two novel molecular defects, chromothripsis and neuritogenesis gene alterations, which frequently occur in high-risk tumours.

Functional MYCN signature predicts outcome of neuroblastoma irrespective of MYCN amplification.

Neuroblastoma is a pediatric tumor of the sympathetic nervous system. MYCN (V-myc myelocytomatosis viral-related oncogene, neuroblastoma derived [avian]) is amplified in 20% of neuroblastomas, and these tumors carry a poor prognosis. However, tumors without MYCN amplification also may have a poor outcome. Here, we identified downstream targets of MYCN by shRNA-mediated silencing MYCN in neuroblastoma cells. From these targets, 157 genes showed an expression profile correlating with MYCN mRNA levels in NB88, a series of 88 neuroblastoma tumors, and therefore represent in vivo relevant MYCN pathway genes. This 157-gene signature identified very poor prognosis tumors in NB88 and independent neuroblastoma cohorts and was more powerful than MYCN amplification or MYCN expression alone. Remarkably, this signature also identified poor outcome of a group of tumors without MYCN amplification. Most of these tumors have low MYCN mRNA levels but high nuclear MYCN protein levels, suggesting stabilization of MYCN at the protein level. One tumor has an MYC amplification and high MYC expression. Chip-on-chip analyses showed that most genes in this signature are directly regulated by MYCN. MYCN induces genes functioning in cell cycle and DNA repair while repressing neuronal differentiation genes. The functional MYCN-157 signature recognizes classical neuroblastoma with MYCN amplification, as well as a newly identified group marked by MYCN protein stabilization.

Inactivation of CDK2 is synthetically lethal to MYCN over-expressing cancer cells.

Two genes have a synthetically lethal relationship when the silencing or inhibiting of 1 gene is only lethal in the context of a mutation or activation of the second gene. This situation offers an attractive therapeutic strategy, as inhibition of such a gene will only trigger cell death in tumor cells with an activated second oncogene but spare normal cells without activation of the second oncogene. Here we present evidence that CDK2 is synthetically lethal to neuroblastoma cells with MYCN amplification and over-expression. Neuroblastomas are childhood tumors with an often lethal outcome. Twenty percent of the tumors have MYCN amplification, and these tumors are ultimately refractory to any therapy. Targeted silencing of CDK2 by 3 RNA interference techniques induced apoptosis in MYCN-amplified neuroblastoma cell lines, but not in MYCN single copy cells. Silencing of MYCN abrogated this apoptotic response in MYCN-amplified cells. Inversely, silencing of CDK2 in MYCN single copy cells did not trigger apoptosis, unless a MYCN transgene was activated. The MYCN induced apoptosis after CDK2 silencing was accompanied by nuclear stabilization of P53, and mRNA profiling showed up-regulation of P53 target genes. Silencing of P53 rescued the cells from MYCN-driven apoptosis. The synthetic lethality of CDK2 silencing in MYCN activated neuroblastoma cells can also be triggered by inhibition of CDK2 with a small molecule drug. Treatment of neuroblastoma cells with roscovitine, a CDK inhibitor, at clinically achievable concentrations induced MYCN-dependent apoptosis. The synthetically lethal relationship between CDK2 and MYCN indicates CDK2 inhibitors as potential MYCN-selective cancer therapeutics.

Urinary 3-Methoxytyramine Is a Biomarker for MYC Activity in Patients With Neuroblastoma.

PURPOSE: Elevated urinary 3-methoxytyramine (3MT) level at diagnosis was recently put forward as independent risk factor for poor prognosis in neuroblastoma. Here, we investigated the biologic basis underlying the putative association between elevated 3MT levels and poor prognosis. METHODS: Urinary 3MT levels and prognosis were investigated in both retrospective Italian (N = 90) and prospective Dutch (N = 95) cohorts. From the Dutch Cancer Oncology Group cohort (N = 122), patients with available urinary 3MT and gene expression data (n = 90) were used to generate a 3MT gene signature. The 3MT gene signature score was then used to predict survival outcome in the Children's Oncology Group (N = 247) and German Pediatric Oncology Group (N = 498) cohorts and compared with other known gene signatures. Immunohistochemistry of MYCN and dopamine ß-hydroxylase proteins was performed on primary tumors. RESULTS: Elevated urinary 3MT levels were associated with poor prognosis in a retrospective cohort and a prospective cohort. Moreover, elevated urinary 3MT levels were associated with eight differentially expressed genes, providing a 3MT gene signature that successfully predicted poor clinical outcome. Even among low-risk patients, high 3MT signature score was associated with poor 5-year overall survival (72% v 99% among low-risk patients with a low 3MT signature score), and the 3MT signature score was correlated with MYC activity in the tumor (R = 82%, P < .0001). Finally, a strong MYCN and weak dopamine ß-hydroxylase staining of tumors derived from patients with elevated urinary 3MT levels was observed, linking MYC activity in the tumor to both catecholamine biosynthesis and elevated urinary 3MT levels. CONCLUSION: Elevated urinary 3MT is a promising biomarker for poor prognosis and reflects increased MYC activity in the tumor. Therefore, urinary 3MT levels should be measured at diagnosis and may assist in assessing risk.

Mesenchymal-Type Neuroblastoma Cells Escape ALK Inhibitors.

Cancer therapy frequently fails due to the emergence of resistance. Many tumors include phenotypically immature tumor cells, which have been implicated in therapy resistance. Neuroblastoma cells can adopt a lineage-committed adrenergic (ADRN) or an immature mesenchymal (MES) state. They differ in epigenetic landscape and transcription factors, and MES cells are more resistant to chemotherapy. Here we analyzed the response of MES cells to targeted drugs. Activating anaplastic lymphoma kinase (ALK) mutations are frequently found in neuroblastoma and ALK inhibitors (ALKi) are in clinical trials. ALKi treatment of ADRN neuroblastoma cells with a tumor-driving ALK mutation induced cell death. Conversely, MES cells did not express either mutant or wild-type ALK and were resistant to ALKi, and MES cells formed tumors that progressed under ALKi therapy. In assessing the role of MES cells in relapse development, TRAIL was identified to specifically induce apoptosis in MES cells and to suppress MES tumor growth. Addition of TRAIL to ALKi treatment of neuroblastoma xenografts delayed relapses in a subset of the animals, suggesting a role for MES cells in relapse formation. While ADRN cells resembled normal embryonal neuroblasts, MES cells resembled immature precursor cells, which also lacked ALK expression. Resistance to targeted drugs can therefore be an intrinsic property of immature cancer cells based on their resemblance to developmental precursors. SIGNIFICANCE: In neuroblastoma, mesenchymal tumor cells lack expression of the tumor-driving ALK oncogene and are resistant to ALKi, but dual treatment with ALKi and mesenchymal cell-targeting TRAIL delays tumor relapse.

A NOTCH feed-forward loop drives reprogramming from adrenergic to mesenchymal state in neuroblastoma.

Transition between differentiation states in development occurs swift but the mechanisms leading to epigenetic and transcriptional reprogramming are poorly understood. The pediatric cancer neuroblastoma includes adrenergic (ADRN) and mesenchymal (MES) tumor cell types, which differ in phenotype, super-enhancers (SEs) and core regulatory circuitries. These cell types can spontaneously interconvert, but the mechanism remains largely unknown. Here, we unravel how a NOTCH3 intracellular domain reprogrammed the ADRN transcriptional landscape towards a MES state. A transcriptional feed-forward circuitry of NOTCH-family transcription factors amplifies the NOTCH signaling levels, explaining the swift transition between two semi-stable cellular states. This transition induces genome-wide remodeling of the H3K27ac landscape and a switch from ADRN SEs to MES SEs. Once established, the NOTCH feed-forward loop maintains the induced MES state. In vivo reprogramming of ADRN cells shows that MES and ADRN cells are equally oncogenic. Our results elucidate a swift transdifferentiation between two semi-stable epigenetic cellular states.

Dickkopf-3 expression is a marker for neuroblastic tumor maturation and is down-regulated by MYCN.

Neuroblastoma and ganglioneuroma are neuroblastic tumors originating from the developing sympathetic peripheral nervous system. Ganglioneuromas are usually benign, while neuroblastomas have a variable prognosis and include very aggressive tumors. Examples exist of neuroblastomas regressing to ganglioneuromas and ganglioneuromas progressing to neuroblastomas. Little is known of the molecular differences between the tumor types. Here we report that Dickkopf-3 (DKK3), a putative extra cellular inhibitor of the Wnt/beta-catenin pathway, showed a strongly differential expression between neuroblastoma and ganglioneuroma. Microarray analyses of 109 neuroblastic tumors revealed that DKK3 is strongly expressed in ganglioneuroma but only weakly in neuroblastoma. Low DKK3 expression in neuroblastoma correlated with a poor prognosis. The expression of DKK3 in the tumor series and in neuroblastoma cell lines was inversely correlated with the expression of the MYCN oncogene. Analysis of 2 neuroblastoma cell lines with inducible activity of MYCN showed that DKK3 is down-regulated by MYCN. We subsequently generated cell lines with inducible expression of DKK3, which revealed an inhibitory effect of DKK3 on proliferation. High DKK3 expression in the benign ganglioneuromas and down-regulation of DKK3 by MYCN in neuroblastoma might contribute to the strongly different clinical behavior of both neuroblastic tumor types.

Direct regulation of the minichromosome maintenance complex by MYCN in neuroblastoma.

The c-Myc and MYCN oncogenes strongly induce cell proliferation. Although a limited series of cell cycle genes were found to be induced by the myc transcription factors, it is still unclear how they mediate the proliferative phenotype. We therefore analysed a neuroblastoma cell line with inducible MYCN expression. We found that all members of the minichromosome maintenance complex (MCM2-7) and MCM8 and MCM10 were up-regulated by MYCN. Expression profiling of 110 neuroblastoma tumours revealed that these genes strongly correlated with MYCN expression in vivo. Extensive chromatin immunoprecipitation experiments were performed to investigate whether the MCM genes were primary MYCN targets. MYCN was bound to the proximal promoters of the MCM2 to -8 genes. These data suggest that MYCN stimulates the expression of not only MCM7, which is a well defined MYCN target gene, but also of the complete minichromosome maintenance complex.

Dickkopf-1 is down-regulated by MYCN and inhibits neuroblastoma cell proliferation.

Neuroblastomas are tumors of the developing peripheral sympathetic nervous system, which originates from the neural crest. Twenty percent of neuroblastomas show amplification of the MYCN oncogene, which correlates with poor prognosis. The MYCN transcription factor can activate and repress gene expression. To broaden our insight in the spectrum of genes down-regulated by MYCN, we generated gene expression profiles of the neuroblastoma cell lines SHEP-21N and SKNAS-NmycER, in which MYCN activity can be regulated. In this study, we show that MYCN suppresses the expression of Dickkopf-1 (DKK1) in both cell lines. DKK1 is a potent inhibitor of the wnt/beta-catenin signalling cascade, which is known to function in neural crest cell migration. We generated a DKK1 inducible cell line, IMR32-DKK1, which showed impaired proliferation upon DKK1 expression. Surprisingly, DKK1 expression did not inhibit the canonical wnt/beta-catenin signalling, suggesting a role of DKK1 in an alternative route of the wnt pathway. Gene expression profiling of two IMR32-DKK1 clones showed that only a few genes, amongst which SYNPO2, were up-regulated by DKK1. SYNPO2 encodes an actin-binding protein and was previously found to inhibit proliferation and invasiveness of prostate cancer cells. These results suggest that MYCN might stimulate cell proliferation by inhibiting the expression of DKK1. DKK1 might exert part of its growth suppressive effect by induction of SYNPO2 expression.

Inhibition of a new differentiation pathway in neuroblastoma by copy number defects of N-myc, Cdc42, and nm23 genes.

The best studied oncogenic mechanisms are inactivating defects in both alleles of tumor suppressor genes and activating mutations in oncogenes. Chromosomal gains and losses are frequent in human tumors, but for many regions, like 1p36 and 17q in neuroblastoma, no mutated tumor suppressor genes or oncogenes were identified. Amplification of N-myc in neuroblastoma is strongly correlated with loss of 1p36 and gain of 17q. Here we report that N-myc down-regulates the mRNA expression of many genes with a role in cell architecture. One of them is the 1p36 gene Cdc42. Restoring the Cdc42 expression in neuroblastoma cells strongly induced differentiation. N-myc also inhibited Cdc42 functioning at the protein level. This was mediated by nm23-H1 and nm23-H2, which are located in the amplified 17q region. Nm23-H1 and nm23-H2 are strongly up-regulated downstream targets of N-myc. Nm23-H1 was shown to bind Cdc42 and prevented the induction of differentiation. Overexpression of Nm23 due to gain of 17q and induction by N-myc combined with weak expression of Cdc42 due to loss of 1p36 and down-regulation by N-myc can thus block differentiation. Although this marks Cdc42 as a candidate tumor suppressor gene, no mutations were found. Further silencing of Cdc42 by small interfering RNA induced massive apoptosis, indicating that tumor cell survival requires a minimal Cdc42 activity. Three regions of chromosomal gain and loss thus affect genes functioning in one pathway in neuroblastoma. They converge to bring the pathway out of balance and prevent Cdc42 mediated differentiation.

Neuroblastoma is composed of two super-enhancer-associated differentiation states.

Neuroblastoma and other pediatric tumors show a paucity of gene mutations, which has sparked an interest in their epigenetic regulation. Several tumor types include phenotypically divergent cells, resembling cells from different lineage development stages. It has been proposed that super-enhancer-associated transcription factor (TF) networks underlie lineage identity, but the role of these enhancers in intratumoral heterogeneity is unknown. Here we show that most neuroblastomas include two types of tumor cells with divergent gene expression profiles. Undifferentiated mesenchymal cells and committed adrenergic cells can interconvert and resemble cells from different lineage differentiation stages. ChIP-seq analysis of isogenic pairs of mesenchymal and adrenergic cells identified a distinct super-enhancer landscape and super-enhancer-associated TF network for each cell type. Expression of the mesenchymal TF PRRX1 could reprogram the super-enhancer and mRNA landscapes of adrenergic cells toward a mesenchymal state. Mesenchymal cells were more chemoresistant in vitro and were enriched in post-therapy and relapse tumors. Two super-enhancer-associated TF networks, which probably mediate lineage control in normal development, thus dominate epigenetic control of neuroblastoma and shape intratumoral heterogeneity.

Full transcriptome analysis of rhabdomyosarcoma, normal, and fetal skeletal muscle: statistical comparison of multiple SAGE libraries.

Rhabdomyosarcoma (RMS) is the most frequent soft tissue sarcoma in children. Improved treatment strategies have increased overall survival, but the response of approximately one-third of the patients is still poor. To increase the knowledge of RMS pathogenesis, we performed the first full transcriptome analysis of RMS using serial analysis of gene expression (SAGE). With a G-test for the simultaneous comparison of subsets of SAGE libraries of normal skeletal muscle, embryonal (ERMS) and alveolar (ARMS) RMS, we identified 251 differentially expressed genes. A literature-mining procedure demonstrated that 158 of these genes have not previously been associated with RMS or normal muscle. Gene Ontology (GO) analysis assigned 198 of the 251 genes to muscle-specific classes, including those involved in normal myogenic development, as well as tumor-related classes. Prominent GO classes were those associated with proliferation and actin reorganization, which are processes that play roles during early muscle development, muscle function, and tumor progression. Using custom microarrays, we confirmed the (up- or down-) regulation of 80% of 98 differentially expressed genes. Another SAGE library of 19- to 22-week-old fetal skeletal muscle was compared with the RMS and normal muscle transcriptomes. Cluster analysis showed that the RMS and fetal muscle SAGE libraries formed one cluster distinct from normal muscle samples. Moreover, the expression profile of 86% of the differentially expressed genes between normal muscle and RMS was highly similar in fetal muscle and RMS. In conclusion, the G-test is a robust tool for analyzing groups of SAGE libraries and correctly identifies genes marking the difference between fully differentiated skeletal muscle and RMS. This study not only substantiates the close association between embryonic myogenesis and RMS development but also provides a rich source of candidate genes to further elucidate the etiology of RMS or to identify diagnostic and/or prognostic markers.

Protein arginine methyltransferase 1 is a novel regulator of MYCN in neuroblastoma.

Amplification or overexpression of MYCN is associated with poor prognosis of human neuroblastoma. We have recently defined a MYCN-dependent transcriptional signature, including protein arginine methyltransferase 1 (PRMT1), which identifies a subgroup of patients with high-risk disease. Here we provide several lines of evidence demonstrating PRMT1 as a novel regulator of MYCN and implicating PRMT1 as a potential therapeutic target in neuroblastoma pathogenesis. First, we observed a strong correlation between MYCN and PRMT1 protein levels in primary neuroblastoma tumors. Second, MYCN physically associates with PRMT1 by direct protein-protein interaction. Third, depletion of PRMT1 through siRNA knockdown reduced neuroblastoma cell viability and MYCN expression. Fourth, we showed that PRMT1 regulates MYCN stability and identified MYCN as a novel substrate of PRMT1. Finally, we demonstrated that mutation of putatively methylated arginine R65 to alanine decreased MYCN stability by altering phosphorylation at residues serine 62 and threonine 58. These results provide mechanistic insights into the modulation of MYCN oncoprotein by PRMT1, and suggest that targeting PRMT1 may have a therapeutic impact on MYCN-driven oncogenesis.

TERT rearrangements are frequent in neuroblastoma and identify aggressive tumors.

Whole-genome sequencing detected structural rearrangements of TERT in 17 of 75 high-stage neuroblastomas, with five cases resulting from chromothripsis. Rearrangements were associated with increased TERT expression and targeted regions immediately up- and downstream of TERT, positioning a super-enhancer close to the breakpoints in seven cases. TERT rearrangements (23%), ATRX deletions (11%) and MYCN amplifications (37%) identify three almost non-overlapping groups of high-stage neuroblastoma, each associated with very poor prognosis.

Ataxia-telangiectasia mutated (ATM) silencing promotes neuroblastoma progression through a MYCN independent mechanism.

Neuroblastoma, a childhood cancer with highly heterogeneous biology and clinical behavior, is characterized by genomic aberrations including amplification of MYCN. Hemizygous deletion of chromosome 11q is a well-established, independent marker of poor prognosis. While 11q22-q23 is the most frequently deleted region, the neuroblastoma tumor suppressor in this region remains to be identified. Chromosome bands 11q22-q23 contain ATM, a cell cycle checkpoint kinase and tumor suppressor playing a pivotal role in the DNA damage response. Here, we report that haploinsufficiency of ATM in neuroblastoma correlates with lower ATM expression, event-free survival, and overall survival. ATM loss occurs in high stage neuroblastoma without MYCN amplification. In SK-N-SH, CLB-Ga and GI-ME-N human neuroblastoma cells, stable ATM silencing promotes neuroblastoma progression in soft agar assays, and in subcutaneous xenografts in nude mice. This effect is dependent on the extent of ATM silencing and does not appear to involve MYCN. Our findings identify ATM as a potential haploinsufficient neuroblastoma tumor suppressor, whose inactivation mirrors the increased aggressiveness associated with 11q deletion in neuroblastoma.

Antitumor activity of sustained N-myc reduction in rhabdomyosarcomas and transcriptional block by antigene therapy.

PURPOSE: Rhabdomyosarcomas are a major cause of cancer death in children, described with MYCN amplification and, in the alveolar subtype, transcription driven by the PAX3-FOXO1 fusion protein. Our aim was to determine the prevalence of N-Myc protein expression and the potential therapeutic effects of reducing expression in rhabdomyosarcomas, including use of an antigene strategy that inhibits transcription. EXPERIMENTAL DESIGN: Protein expression was assessed by immunohistochemistry. MYCN expression was reduced in representative cell lines by RNA interference and an antigene peptide nucleic acid (PNA) oligonucleotide conjugated to a nuclear localization signal peptide. Associated gene expression changes, cell viability, and apoptosis were analyzed in vitro. As a paradigm for antigene therapy, the effects of systemic treatment of mice with rhabdomyosarcoma cell line xenografts were determined. RESULTS: High N-Myc levels were significantly associated with genomic amplification, presence of the PAX3/7-FOXO1 fusion genes, and proliferative capacity. Sustained reduction of N-Myc levels in all rhabdomyosarcoma cell lines that express the protein decreased cell proliferation and increased apoptosis. Positive feedback was shown to regulate PAX3-FOXO1 and N-Myc levels in the alveolar subtype that critically decrease PAX3-FOXO1 levels on reducing N-Myc. Pharmacologic systemic administration of the antigene PNA can eliminate alveolar rhabdomyosarcoma xenografts in mice, without relapse or toxicity. CONCLUSION: N-Myc, with its restricted expression in non-fetal tissues, is a therapeutic target to treat rhabdomyosarcomas, and blocking gene transcription using antigene oligonucleotide strategies has therapeutic potential in the treatment of cancer and other diseases that has not been previously realized in vivo.

LIN28B induces neuroblastoma and enhances MYCN levels via let-7 suppression.

LIN28B regulates developmental processes by modulating microRNAs (miRNAs) of the let-7 family. A role for LIN28B in cancer has been proposed but has not been established in vivo. Here, we report that LIN28B showed genomic aberrations and extensive overexpression in high-risk neuroblastoma compared to several other tumor entities and normal tissues. High LIN28B expression was an independent risk factor for adverse outcome in neuroblastoma. LIN28B signaled through repression of the let-7 miRNAs and consequently resulted in elevated MYCN protein expression in neuroblastoma cells. LIN28B-let-7-MYCN signaling blocked differentiation of normal neuroblasts and neuroblastoma cells. These findings were fully recapitulated in a mouse model in which LIN28B expression in the sympathetic adrenergic lineage induced development of neuroblastomas marked by low let-7 miRNA levels and high MYCN protein expression. Interference with this pathway might offer therapeutic perspectives.

Read-through transcript from NM23-H1 into the neighboring NM23-H2 gene encodes a novel protein, NM23-LV.

NM23-H1 and NM23-H2 are neighboring genes on chromosome 17q. They encode nucleoside diphosphate kinases that have additional roles in signal transduction, transcription, and apoptosis. NM23-H1 expression is a strong marker for prognosis and metastatic behavior in many tumor types. A new bioinformatic tool, TranscriptView, identified read-through transcripts that start in the NM23-H1 gene and continue in the neighboring NM23-H2 gene. Splicing results in a transcript containing exons 1 to 4 of NM23-H1 and exons 2 to 5 of NM23-H2. The resulting mRNA encodes a novel and long variant of the NM23 protein family, NM23-LV, which contains part of NM23-H1 and the complete NM23-H2 protein. The transcript was amplified and sequenced from two neuroblastoma cell lines, confirming the presence of the predicted NM23-LV mRNA in vivo. Tissue analysis showed that NM23-LV is ubiquitously expressed, with the exception of the kidney. Neuroblastoma tumors show high-level expression of NM23-H1 and-H2 as well as NM23-LV mRNA. In neuroblastoma cells, the NM23-LV protein has mainly a cytoplasmic localization, but some nuclear staining was observed as well.