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.

BackgroundMedulloblastoma is the most common malignant brain tumor in children. Despite recent improvements in cure rates, prediction of disease outcome remains a major challenge and survivors suffer from serious therapy-related side-effects. Recent data showed that patients with WNT-activated tumors have a favorable prognosis, suggesting that these patients could be treated less intensively, thereby reducing the side-effects. This illustrates the potential benefits of a robust classification of medulloblastoma patients and a detailed knowledge of associated biological mechanisms.Methods and FindingsTo get a better insight into the molecular biology of medulloblastoma we established mRNA expression profiles of 62 medulloblastomas and analyzed 52 of them also by comparative genomic hybridization (CGH) arrays. Five molecular subtypes were identified, characterized by WNT signaling (A; 9 cases), SHH signaling (B; 15 cases), expression of neuronal differentiation genes (C and D; 16 and 11 cases, respectively) or photoreceptor genes (D and E; both 11 cases). Mutations in β-catenin were identified in all 9 type A tumors, but not in any other tumor. PTCH1 mutations were exclusively identified in type B tumors. CGH analysis identified several fully or partly subtype-specific chromosomal aberrations. Monosomy of chromosome 6 occurred only in type A tumors, loss of 9q mostly occurred in type B tumors, whereas chromosome 17 aberrations, most common in medulloblastoma, were strongly associated with type C or D tumors. Loss of the inactivated X-chromosome was highly specific for female cases of type C, D and E tumors. Gene expression levels faithfully reflected the chromosomal copy number changes. Clinicopathological features significantly different between the 5 subtypes included metastatic disease and age at diagnosis and histology. Metastatic disease at diagnosis was significantly associated with subtypes C and D and most strongly with subtype E. Patients below 3 yrs of age had type B, D, or E tumors. Type B included most desmoplastic cases. We validated and confirmed the molecular subtypes and their associated clinicopathological features with expression data from a second independent series of 46 medulloblastomas.ConclusionsThe new medulloblastoma classification presented in this study will greatly enhance the understanding of this heterogeneous disease. It will enable a better selection and evaluation of patients in clinical trials, and it will support the development of new molecular targeted therapies. Ultimately, our results may lead to more individualized therapies with improved cure rates and a better quality of life.

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.

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.