Danny A. Zwijnenburg

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.

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.

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.

Evaluation of the similarity of gene expression data estimated with SAGE and Affymetrix GeneChips

BackgroundSerial Analysis of Gene Expression (SAGE) and microarrays have found awidespread application, but much ambiguity exists regarding the evaluation of these technologies. Cross-platform utilization of gene expression data from the SAGE and microarray technology could reduce the need for duplicate experiments and facilitate a more extensive exchange of data within the research community. This requires a measure for the correspondence of the different gene expression platforms. To date, a number of cross-platform evaluations (including a few studies using SAGE and Affymetrix GeneChips) have been conducted showing a variable, but overall low, concordance. This study evaluates these overall measures and introduces the between-ratio difference as a concordance measure pergene.ResultsIn this study, gene expression measurements of Unigene clusters represented by both Affymetrix GeneChips HG-U133A and SAGE were compared using two independent RNA samples. After matching of the data sets the final comparison contains a small data set of 1094 unique Unigene clusters, which is unbiased with respect to expression level. Different overall correlation approaches, like Up/Down classification, contingency tables and correlation coefficients were used to compare both platforms. In addition, we introduce a novel approach to compare two platforms based on the calculation of differences between expression ratios observed in each platform for each individual transcript. This approach results in a concordance measure per gene (with statistical probability value), as opposed to the commonly used overall concordance measures between platforms.ConclusionWe can conclude that intra-platform correlations are generally good, but that overall agreement between the two platforms is modest. This might be due to the binomially distributed sampling variation in SAGE tag counts, SAGE annotation errors and the intensity variation between probe sets of a single gene in Affymetrix GeneChips. We cannot identify or advice which platform performs better since both have their (dis)-advantages. Therefore it is strongly recommended to perform follow-up studies of interesting genes using additional techniques. The newly introduced between-ratio difference is a filtering-independent measure for between-platform concordance. Moreover, the between-ratio difference per gene can be used to detect transcripts with similar regulation on both platforms.

OTX2 directly activates cell cycle genes and inhibits differentiation in medulloblastoma cells

The transcription factor OTX2 has been implicated as an oncogene in medulloblastoma, which is the most common malignant brain tumor in children. It is highly expressed in most medulloblastomas and amplified in a subset of them. To study the role OTX2 has in medulloblastoma we investigated the downstream pathway of OTX2. We generated D425 medulloblastoma cells in which endogenous OTX2 can be silenced by inducible shRNA. Silencing of OTX2 strongly inhibited cell proliferation and resulted in a neuronal‐like differentiation. Expression profiling of time courses after silencing showed a progressive change in gene expression for many cellular processes. Downregulated genes were highly enriched for cell cycle and visual perception genes, while upregulated genes were enriched for genes involved in development and differentiation. This shift is reminiscent of expression changes described during normal cerebellum development where proliferating granule progenitor cells have high OTX2 expression, which diminishes when these cells exit the cell cycle and start to differentiate. ChIP‐on‐chip analyses of OTX2 in D425 cells identified cell cycle and perception genes as direct OTX2 targets, while regulation of most differentiation genes appeared to be indirect. The expression of many directly regulated genes correlated to OTX2 expression in primary tumors, suggesting the in vivo relevance of these genes and their potential as targets for therapeutic intervention. These analyses provide more insight in the molecular network of OTX2, demonstrating that OTX2 is essential in medulloblastoma and directly drives proliferation by regulation of cell cycle genes.

Silencing of SPRY1 Triggers Complete Regression of Rhabdomyosarcoma Tumors Carrying a Mutated RAS Gene

RAS oncogenes are among the most frequently mutated genes in human cancer, but effective strategies for therapeutic inhibition of the RAS pathway have been elusive. Sprouty1 (SPRY1) is an upstream antagonist of RAS that is activated by extracellular signal-related kinase (ERK), providing a negative feedback loop for RAS signaling, and other evidence suggests that SPRY1 may have a tumor suppressor function. Studies of RAS status in the human childhood tumor rhabdomyosarcoma (RMS) indicated mutations in approximately half of the tumors of the embryonal rhabdomyosarcoma subtype (ERMS) but not the alveolar subtype (ARMS). ERMS tumors also showed overexpression of SPRY1, which was indeed upregulated by mutant RAS. However, we found that, in the presence of mutant RAS, the function of SPRY1 was changed from an antagonist to an agonist of RAS signaling. Thus, SPRY1 supported formation of activated ERK and mitogen-activated protein/ERK kinase and was essential for ERMS cell proliferation and survival. Conversely, silencing of SPRY1 in ERMS cells (but not ARMS cells) abolished their tumorigenicity in mice. Moreover, silencing of SPRY1 caused regression of established ERMS tumors (but not ARMS tumors) formed in xenograft settings. Our findings argue that SPRY1 inhibition can offer a therapeutic strategy to treat childhood RMS and possibly other tumors carrying oncogenic RAS mutations.

Regulation of Cell Cycle Genes and Induction of Senescence by Overexpression of OTX2 in Medulloblastoma Cell Lines

The transcription factor orthodenticle homeobox 2 (OTX2) has been implicated in the pathogenesis of medulloblastoma, as it is often highly expressed and sometimes amplified in these tumors. Little is known of the downstream pathways regulated by OTX2. We therefore generated MED8A and DAOY medulloblastoma cell lines with doxycycline-inducible OTX2 expression. In both cell lines, OTX2 inhibited proliferation and induced a senescence-like phenotype with senescence-associated β-galactosidase activity. Expression profiles of time series after OTX2 induction in MED8A showed early upregulation of cell cycle genes related to the G2-M phase, such as AURKA, CDC25C, and CCNG2. Paradoxically, G1-S phase genes such as MYC, CDK4, CDK6, CCND1, and CCND2 were strongly downregulated, in line with the observed G1 arrest. ChIP-on-chip analyses of OTX2 binding to promoter regions in MED8A and DAOY showed a strong enrichment for binding to the G2-M genes, suggesting a direct activation. Their mRNA expression correlated with OTX2 expression in primary tumors, underscoring the in vivo relevance of this regulation. OTX2 induction activated the P53 pathway in MED8A, but not in DAOY, which carries a mutated P53 gene. In DAOY cells, senescence-associated secretory factors, such as interleukin-6 and insulin-like growth factor binding protein 7, were strongly upregulated after OTX2 induction. We hypothesize that the imbalance in cell cycle stimulation by OTX2 leads to cellular senescence either by activating the P53 pathway or through the induction of secretory factors. Our data indicate that OTX2 directly induces a series of cell cycle genes but requires cooperating genes for an oncogenic acceleration of the cell cycle. Mol Cancer Res; 8(10); 1344–57. ©2010 AACR.

OTX2 sustains a bivalent-like state of OTX2-bound promoters in medulloblastoma by maintaining their H3K27me3 levels

Recent studies showed frequent mutations in histone H3 lysine 27 (H3K27) demethylases in medulloblastomas of Group 3 and Group 4, suggesting a role for H3K27 methylation in these tumors. Indeed, trimethylated H3K27 (H3K27me3) levels were shown to be higher in Group 3 and 4 tumors compared to WNT and SHH medulloblastomas, also in tumors without detectable mutations in demethylases. Here, we report that polycomb genes, required for H3K27 methylation, are consistently upregulated in Group 3 and 4 tumors. These tumors show high expression of the homeobox transcription factor OTX2. Silencing of OTX2 in D425 medulloblastoma cells resulted in downregulation of polycomb genes such as EZH2, EED, SUZ12 and RBBP4 and upregulation of H3K27 demethylases KDM6A, KDM6B,JARID2 and KDM7A. This was accompanied by decreased H3K27me3 and increased H3K27me1 levels in promoter regions. Strikingly, the decrease of H3K27me3 was most prominent in promoters that bind OTX2. OTX2-bound promoters showed high levels of the H3K4me3 and H3K9ac activation marks and intermediate levels of the H3K27me3 inactivation mark, reminiscent of a bivalent modification. After silencing of OTX2, H3K27me3 levels strongly dropped, but H3K4me3 and H3K9ac levels remained high. OTX2-bound bivalent genes showed high expression levels in D425, but the expression of most of these genes did not change after OTX2 silencing and loss of the H3K27me3 mark. Maintaining promoters in a bivalent state by sustaining H3K27 trimethylation therefore seems to be an important function of OTX2 in medulloblastoma, while other transcription factors might regulate the actual expression levels of these genes.

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.

Joint binding of OTX2 and MYC in promotor regions is associated with high gene expression in medulloblastoma.

Both OTX2 and MYC are important oncogenes in medulloblastoma, the most common malignant brain tumor in childhood. Much is known about MYC binding to promoter regions, but OTX2 binding is hardly investigated. We used ChIP-on-chip data to analyze the binding patterns of both transcription factors in D425 medulloblastoma cells. When combining the data for all promoter regions in the genome, OTX2 binding showed a remarkable bi-modal distribution pattern with peaks around −250 bp upstream and +650 bp downstream of the transcription start sites (TSSs). Indeed, 40.2% of all OTX2-bound TSSs had more than one significant OTX2-binding peak. This OTX2-binding pattern was very different from the TSS-centered single peak binding pattern observed for MYC and other known transcription factors. However, in individual promoter regions, OTX2 and MYC have a strong tendency to bind in proximity of each other. OTX2-binding sequences are depleted near TSSs in the genome, providing an explanation for the observed bi-modal distribution of OTX2 binding. This contrasts to the enrichment of E-box sequences at TSSs. Both OTX2 and MYC binding independently correlated with higher gene expression. Interestingly, genes of promoter regions with multiple OTX2 binding as well as MYC binding showed the highest expression levels in D425 cells and in primary medulloblastomas. Genes within this class of promoter regions were enriched for medulloblastoma and stem cell specific genes. Our data suggest an important functional interaction between OTX2 and MYC in regulating gene expression in medulloblastoma.