Prospective derivation of a living organoid biobank of colorectal cancer patients.

In Rspondin-based 3D cultures, Lgr5 stem cells from multiple organs form ever-expanding epithelial organoids that retain their tissue identity. We report the establishment of tumor organoid cultures from 20 consecutive colorectal carcinoma (CRC) patients. For most, organoids were also generated from adjacent normal tissue. Organoids closely recapitulate several properties of the original tumor. The spectrum of genetic changes within the "living biobank" agrees well with previous large-scale mutational analyses of CRC. Gene expression analysis indicates that the major CRC molecular subtypes are represented. Tumor organoids are amenable to high-throughput drug screens allowing detection of gene-drug associations. As an example, a single organoid culture was exquisitely sensitive to Wnt secretion (porcupine) inhibitors and carried a mutation in the negative Wnt feedback regulator RNF43, rather than in APC. Organoid technology may fill the gap between cancer genetics and patient trials, complement cell-line- and xenograft-based drug studies, and allow personalized therapy design. PAPERCLIP.

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.

Dissecting the genomic complexity underlying medulloblastoma.

Medulloblastoma is an aggressively growing tumour, arising in the cerebellum or medulla/brain stem. It is the most common malignant brain tumour in children, and shows tremendous biological and clinical heterogeneity. Despite recent treatment advances, approximately 40% of children experience tumour recurrence, and 30% will die from their disease. Those who survive often have a significantly reduced quality of life. Four tumour subgroups with distinct clinical, biological and genetic profiles are currently identified. WNT tumours, showing activated wingless pathway signalling, carry a favourable prognosis under current treatment regimens. SHH tumours show hedgehog pathway activation, and have an intermediate prognosis. Group 3 and 4 tumours are molecularly less well characterized, and also present the greatest clinical challenges. The full repertoire of genetic events driving this distinction, however, remains unclear. Here we describe an integrative deep-sequencing analysis of 125 tumour-normal pairs, conducted as part of the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. Tetraploidy was identified as a frequent early event in Group 3 and 4 tumours, and a positive correlation between patient age and mutation rate was observed. Several recurrent mutations were identified, both in known medulloblastoma-related genes (CTNNB1, PTCH1, MLL2, SMARCA4) and in genes not previously linked to this tumour (DDX3X, CTDNEP1, KDM6A, TBR1), often in subgroup-specific patterns. RNA sequencing confirmed these alterations, and revealed the expression of what are, to our knowledge, the first medulloblastoma fusion genes identified. Chromatin modifiers were frequently altered across all subgroups. These findings enhance our understanding of the genomic complexity and heterogeneity underlying medulloblastoma, and provide several potential targets for new therapeutics, especially for Group 3 and 4 patients.

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.

Copy number defects of G1-cell cycle genes in neuroblastoma are frequent and correlate with high expression of E2F target genes and a poor prognosis.

The tightly controlled network of cell cycle genes consists of a core of cyclin dependent kinases (CDKs) that are activated by periodically expressed cyclins. The activity of the cyclin-CDK complexes is regulated by cyclin dependent kinase inhibitors (CDKIs) and multiple signal transduction routes that converge on the cell cycle. Neuroblastoma are pediatric tumors that belong to the group of small round blue cell tumors, characterized by a fast proliferation. Here, we present high throughput analyses of cell cycle regulating genes in neuroblastoma. We analyzed a series of 82 neuroblastomas by comparative genomic hybridization arrays, single nucleotide polymorphism arrays, and Affymetrix expression arrays and analyzed the datasets in parallel with the R2 bioinformatic tool (http://r2.amc.nl). About 30% of the tumors had genomic amplifications, gains, or losses with shortest regions of overlap that suggested implication of a series of G1 cell cycle regulating genes. CCND1 (cyclin D1) and CDK4 were amplified or gained and the chromosomal regions containing the CDKN2 (INK4) group of CDKIs were frequently deleted. Cluster analysis showed that tumors with genomic aberrations in G1 regulating genes over-expressed E2F target genes, which regulate S and G2/M phase progression. These tumors have a poor prognosis. Our findings suggest that pharmacological inhibition of cell cycle genes might bear therapeutic promises for patients with high risk neuroblastoma.

LIN28A immunoreactivity is a potent diagnostic marker of embryonal tumor with multilayered rosettes (ETMR).

Embryonal tumor with multilayered rosettes (ETMR, previously known as ETANTR) is a highly aggressive embryonal CNS tumor, which almost exclusively affects infants and is associated with a dismal prognosis. Accurate diagnosis is of critical clinical importance because of its poor response to current treatment protocols and its distinct biology. Amplification of the miRNA cluster at 19q13.42 has been identified previously as a genetic hallmark for ETMR, but an immunohistochemistry-based assay for clinical routine diagnostics [such as INI-1 for atypical teratoid rhabdoid tumor (AT/RT)] is still lacking. In this study, we screened for an ETMR-specific marker using a gene-expression profiling dataset of more than 1,400 brain tumors and identified LIN28A as a highly specific marker for ETMR. The encoded protein binds small RNA and has been implicated in stem cell pluripotency, metabolism and tumorigenesis. Using an LIN28A specific antibody, we carried out immunohistochemical analysis of LIN28A in more than 800 childhood brain-tumor samples and confirmed its high specificity for ETMR. Strong LIN28A immunoexpression was found in all 37 ETMR samples tested, whereas focal reactivity was only present in a small (6/50) proportion of AT/RT samples. All other pediatric brain tumors were completely LIN28A-negative. In summary, we established LIN28A immunohistochemistry as a highly sensitive and specific, rapid, inexpensive diagnostic tool for routine pathological verification of ETMR.

Identification of BIRC6 as a novel intervention target for neuroblastoma therapy.

BACKGROUND: Neuroblastoma are pediatric tumors of the sympathetic nervous system with a poor prognosis. Apoptosis is often deregulated in cancer cells, but only a few defects in apoptotic routes have been identified in neuroblastoma. METHODS: Here we investigated genomic aberrations affecting genes of the intrinsic apoptotic pathway in neuroblastoma. We analyzed DNA profiling data (CGH and SNP arrays) and mRNA expression data of 31 genes of the intrinsic apoptotic pathway in a dataset of 88 neuroblastoma tumors using the R2 bioinformatic platform ( http://r2.amc.nl). BIRC6 was selected for further analysis as a tumor driving gene. Knockdown experiments were performed using BIRC6 lentiviral shRNA and phenotype responses were analyzed by Western blot and MTT-assays. In addition, DIABLO levels and interactions were investigated with immunofluorescence and co-immunoprecipitation. RESULTS: We observed frequent gain of the BIRC6 gene on chromosome 2, which resulted in increased mRNA expression. BIRC6 is an inhibitor of apoptosis protein (IAP), that can bind and degrade the cytoplasmic fraction of the pro-apoptotic protein DIABLO. DIABLO mRNA expression was exceptionally high in neuroblastoma but the protein was only detected in the mitochondria. Upon silencing of BIRC6 by shRNA, DIABLO protein levels increased and cells went into apoptosis. Co-immunoprecipitation confirmed direct interaction between DIABLO and BIRC6 in neuroblastoma cell lines. CONCLUSION: Our findings indicate that BIRC6 may have a potential oncogenic role in neuroblastoma by inactivating cytoplasmic DIABLO. BIRC6 inhibition may therefore provide a means for therapeutic intervention in neuroblastoma.

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.

Dissection of the oncogenic MYCN transcriptional network reveals a large set of clinically relevant cell cycle genes as drivers of neuroblastoma tumorigenesis.

Amplification of the oncogenic transcription factor MYCN plays a major role in the pathogenesis of several pediatric cancers, including neuroblastoma, medulloblastoma, and rhabodomyosarcoma. For neuroblastoma, MYCN amplification is the most powerful genetic predictor of poor patient survival, yet the mechanism by which MYCN drives tumorigenesis is only partially understood. To gain an insight into the distribution of MYCN binding and to identify clinically relevant MYCN target genes, we performed an integrated analysis of MYCN ChIP-chip and mRNA expression using the MYCN repressible SHEP-21N neuroblastoma cell line. We hypothesized that genes exclusively MYCN bound in SHEP-21N cells over-expressing MYCN would be enriched for direct targets which contribute to the process of disease progression. Integrated analysis revealed that MYCN drives tumorigenesis predominantly as a positive regulator of target gene transcription. A high proportion of genes (24%) that are MYCN bound and up-regulated in the SHEP-21N model are significantly associated with poor overall patient survival (OS) in a set of 88 tumors. In contrast, the proportion of genes down-regulated when bound by MYCN in the SHEP-21N model and which are significantly associated with poor overall patient survival when under-expressed in primary tumors was significantly lower (5%). Gene ontology analysis determined a highly statistically significant enrichment for cell cycle related genes within the over-expressed MYCN target group which were also associated with poor OS. We conclude that the over-expression of MYCN leads to aberrant binding and over-expression of genes associated with cell cycle regulation which are significantly correlated with poor OS and MYCN amplification.

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.

RAS-MAPK Pathway-Driven Tumor Progression Is Associated with Loss of CIC and Other Genomic Aberrations in Neuroblastoma.

Mutations affecting the RAS-MAPK pathway frequently occur in relapsed neuroblastoma tumors, which suggests that activation of this pathway is associated with a more aggressive phenotype. To explore this hypothesis, we generated several model systems to define a neuroblastoma RAS-MAPK pathway signature. Activation of this pathway in primary tumors indeed correlated with poor survival and was associated with known activating mutations in ALK and other RAS-MAPK pathway genes. Integrative analysis showed that mutations in PHOX2B, CIC, and DMD were also associated with an activated RAS-MAPK pathway. Mutation of PHOX2B and deletion of CIC in neuroblastoma cell lines induced activation of the RAS-MAPK pathway. This activation was independent of phosphorylated ERK in CIC knockout systems. Furthermore, deletion of CIC caused a significant increase in tumor growth in vivo These results show that the RAS-MAPK pathway is involved in tumor progression and establish CIC as a powerful tumor suppressor that functions downstream of this pathway in neuroblastoma.Significance: This work identifies CIC as a powerful tumor suppressor affecting the RAS-MAPK pathway in neuroblastoma and reinforces the importance of mutation-driven activation of this pathway in cancer. Cancer Res; 78(21); 6297-307. ©2018 AACR.

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.

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.

OTX1 and OTX2 expression correlates with the clinicopathologic classification of medulloblastomas.

OTX1 and OTX2 are transcription factors with an essential role in the development of the cerebellum. We previously described a high OTX2 expression in medulloblastoma. Here, we analyzed amplification and mRNA expression of OTX1 and OTX2 in a series of human medulloblastomas. In addition, OTX2 protein expression was analyzed on tissue arrays. The OTX2 gene was amplified in the medulloblastoma cell line D425 and mRNA and protein data showed expression in 114 of 152 medulloblastomas (75%), but not in postnatal cerebellum. Northern blot (n = 10) and reverse transcriptase-polymerase chain reaction (n = 45) analyses demonstrated that virtually all medulloblastomas expressed OTX1, OTX2, or both. OTX2 mRNA expression correlated with a classic medulloblastoma histology (29 of 34 cases), whereas expression of OTX1 mRNA only was correlated with a nodular/desmoplastic histology (9 of 11 cases). Immunohistochemical analysis of a series of classic medulloblastomas detected OTX2 protein expression in 83 of 107 (78%) cases. The OTX2-positive tumors of this series were preferentially localized in the vermis of the cerebellum, whereas OTX2-negative tumors more frequently occurred in the hemispheres of the cerebellum. In addition, OTX2-positive tumors were mainly found in children, but OTX2-negative tumors occurred in 2 patient groups: very young patients (<5 years) and adults (>20 years). Nodular/desmoplastic medulloblastomas are thought to arise from the external granular layer (EGL). However, it is unclear whether classic medulloblastomas also originate from the EGL or from the ventricular matrix. Analysis of human fetal brain showed OTX2 protein expression in a small number of presumptive neuronal precursor cells of the EGL, but not in precursor cells of the ventricular matrix. Combined with data from rodents, our results therefore suggest that both nodular/desmoplastic and at least part of the classic medulloblastomas originate from cells of the EGL, albeit from different regions.

Atypical Teratoid/Rhabdoid Tumors Are Comprised of Three Epigenetic Subgroups with Distinct Enhancer Landscapes.

Atypical teratoid/rhabdoid tumor (ATRT) is one of the most common brain tumors in infants. Although the prognosis of ATRT patients is poor, some patients respond favorably to current treatments, suggesting molecular inter-tumor heterogeneity. To investigate this further, we genetically and epigenetically analyzed 192 ATRTs. Three distinct molecular subgroups of ATRTs, associated with differences in demographics, tumor location, and type of SMARCB1 alterations, were identified. Whole-genome DNA and RNA sequencing found no recurrent mutations in addition to SMARCB1 that would explain the differences between subgroups. Whole-genome bisulfite sequencing and H3K27Ac chromatin-immunoprecipitation sequencing of primary tumors, however, revealed clear differences, leading to the identification of subgroup-specific regulatory networks and potential therapeutic targets.

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.

Mutations in the SIX1/2 pathway and the DROSHA/DGCR8 miRNA microprocessor complex underlie high-risk blastemal type Wilms tumors.

Blastemal histology in chemotherapy-treated pediatric Wilms tumors (nephroblastoma) is associated with adverse prognosis. To uncover the underlying tumor biology and find therapeutic leads for this subgroup, we analyzed 58 blastemal type Wilms tumors by exome and transcriptome sequencing and validated our findings in a large replication cohort. Recurrent mutations included a hotspot mutation (Q177R) in the homeo-domain of SIX1 and SIX2 in tumors with high proliferative potential (18.1% of blastemal cases); mutations in the DROSHA/DGCR8 microprocessor genes (18.2% of blastemal cases); mutations in DICER1 and DIS3L2; and alterations in IGF2, MYCN, and TP53, the latter being strongly associated with dismal outcome. DROSHA and DGCR8 mutations strongly altered miRNA expression patterns in tumors, which was functionally validated in cell lines expressing mutant DROSHA.

Relapsed neuroblastomas show frequent RAS-MAPK pathway mutations.

The majority of patients with neuroblastoma have tumors that initially respond to chemotherapy, but a large proportion will experience therapy-resistant relapses. The molecular basis of this aggressive phenotype is unknown. Whole-genome sequencing of 23 paired diagnostic and relapse neuroblastomas showed clonal evolution from the diagnostic tumor, with a median of 29 somatic mutations unique to the relapse sample. Eighteen of the 23 relapse tumors (78%) showed mutations predicted to activate the RAS-MAPK pathway. Seven of these events were detected only in the relapse tumor, whereas the others showed clonal enrichment. In neuroblastoma cell lines, we also detected a high frequency of activating mutations in the RAS-MAPK pathway (11/18; 61%), and these lesions predicted sensitivity to MEK inhibition in vitro and in vivo. Our findings provide a rationale for genetic characterization of relapse neuroblastomas and show that RAS-MAPK pathway mutations may function as a biomarker for new therapeutic approaches to refractory disease.