SOX11 regulates SWI/SNF complex components as member of the adrenergic neuroblastoma core regulatory circuitry.

The pediatric extra-cranial tumor neuroblastoma displays a low mutational burden while recurrent copy number alterations are present in most high-risk cases. Here, we identify SOX11 as a dependency transcription factor in adrenergic neuroblastoma based on recurrent chromosome 2p focal gains and amplifications, specific expression in the normal sympatho-adrenal lineage and adrenergic neuroblastoma, regulation by multiple adrenergic specific (super-)enhancers and strong dependency on high SOX11 expression in adrenergic neuroblastomas. SOX11 regulated direct targets include genes implicated in epigenetic control, cytoskeleton and neurodevelopment. Most notably, SOX11 controls chromatin regulatory complexes, including 10 SWI/SNF core components among which SMARCC1, SMARCA4/BRG1 and ARID1A. Additionally, the histone deacetylase HDAC2, PRC1 complex component CBX2, chromatin-modifying enzyme KDM1A/LSD1 and pioneer factor c-MYB are regulated by SOX11. Finally, SOX11 is identified as a core transcription factor of the core regulatory circuitry (CRC) in adrenergic high-risk neuroblastoma with a potential role as epigenetic master regulator upstream of the CRC.

Publisher Correction: In silico discovery of a FOXM1 driven embryonal signaling pathway in therapy resistant neuroblastoma tumors.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

In silico discovery of a FOXM1 driven embryonal signaling pathway in therapy resistant neuroblastoma tumors.

Chemotherapy resistance is responsible for high mortality rates in neuroblastoma. MYCN, an oncogenic driver in neuroblastoma, controls pluripotency genes including LIN28B. We hypothesized that enhanced embryonic stem cell (ESC) gene regulatory programs could mark tumors with high pluripotency capacity and subsequently increased risk for therapy failure. An ESC miRNA signature was established based on publicly available data. In addition, an ESC mRNA signature was generated including the 500 protein coding genes with the highest positive expression correlation with the ESC miRNA signature score in 200 neuroblastomas. High ESC m(i)RNA expression signature scores were significantly correlated with poor neuroblastoma patient outcome specifically in the subgroup of MYCN amplified tumors and stage 4 nonamplified tumors. Further data-mining identified FOXM1, as the major predicted driver of this ESC signature, controlling a large set of genes implicated in cell cycle control and DNA damage response. Of further interest, re-analysis of published data showed that MYCN transcriptionally activates FOXM1 in neuroblastoma cells. In conclusion, a novel ESC m(i)RNA signature stratifies neuroblastomas with poor prognosis, enabling the identification of therapy-resistant tumors. The finding that this signature is strongly FOXM1 driven, warrants for drug design targeted at FOXM1 or key components controlling this pathway.

Network Modeling of microRNA-mRNA Interactions in Neuroblastoma Tumorigenesis Identifies miR-204 as a Direct Inhibitor of MYCN.

Neuroblastoma is a pediatric cancer of the sympathetic nervous system where MYCN amplification is a key indicator of poor prognosis. However, mechanisms by which MYCN promotes neuroblastoma tumorigenesis are not fully understood. In this study, we analyzed global miRNA and mRNA expression profiles of tissues at different stages of tumorigenesis from TH-MYCN transgenic mice, a model of MYCN-driven neuroblastoma. On the basis of a Bayesian learning network model in which we compared pretumor ganglia from TH-MYCN+/+ mice to age-matched wild-type controls, we devised a predicted miRNA-mRNA interaction network. Among the miRNA-mRNA interactions operating during human neuroblastoma tumorigenesis, we identified miR-204 as a tumor suppressor miRNA that inhibited a subnetwork of oncogenes strongly associated with MYCN-amplified neuroblastoma and poor patient outcome. MYCN bound to the miR-204 promoter and repressed miR-204 transcription. Conversely, miR-204 directly bound MYCN mRNA and repressed MYCN expression. miR-204 overexpression significantly inhibited neuroblastoma cell proliferation in vitro and tumorigenesis in vivo Together, these findings identify novel tumorigenic miRNA gene networks and miR-204 as a tumor suppressor that regulates MYCN expression in neuroblastoma tumorigenesis.Significance: Network modeling of miRNA-mRNA regulatory interactions in a mouse model of neuroblastoma identifies miR-204 as a tumor suppressor and negative regulator of MYCN. Cancer Res; 78(12); 3122-34. ©2018 AACR.

Genome wide expression profiling of p53 regulated miRNAs in neuroblastoma.

Restoration of the antitumor activity of p53 could offer a promising approach for the treatment of neuroblastoma. MicroRNAs (miRNAs) are important mediators of p53 activity, but their role in the p53 response has not yet been comprehensively addressed in neuroblastoma. Therefore, we set out to characterize alterations in miRNA expression that are induced by p53 activation in neuroblastoma cells. Genome-wide miRNA expression analysis showed that miR-34a-5p, miR-182-5p, miR-203a, miR-222-3p, and miR-432-5p are upregulated following nutlin-3 treatment in a p53 dependent manner. The function of miR-182-5p, miR-203a, miR-222-3p, and miR-432-5p was analyzed by ectopic overexpression of miRNA mimics. We observed that these p53-regulated miRNAs inhibit the proliferation of neuroblastoma cells to varying degrees, with the most profound growth inhibition recorded for miR-182-5p. Overexpression of miR-182-5p promoted apoptosis in some neuroblastoma cell lines and induced neuronal differentiation of NGP cells. Using Chromatin Immunoprecipitation-qPCR (ChIP-qPCR), we did not observe direct binding of p53 to MIR182, MIR203, MIR222, and MIR432 in neuroblastoma cells. Taken together, our findings yield new insights in the network of p53-regulated miRNAs in neuroblastoma.

Upregulation of MAPK Negative Feedback Regulators and RET in Mutant ALK Neuroblastoma: Implications for Targeted Treatment.

PURPOSE: Activating ALK mutations are present in almost 10% of primary neuroblastomas and mark patients for treatment with small-molecule ALK inhibitors in clinical trials. However, recent studies have shown that multiple mechanisms drive resistance to these molecular therapies. We anticipated that detailed mapping of the oncogenic ALK-driven signaling in neuroblastoma can aid to identify potential fragile nodes as additional targets for combination therapies. EXPERIMENTAL DESIGN: To achieve this goal, transcriptome profiling was performed in neuroblastoma cell lines with the ALK(F1174L) or ALK(R1275Q) hotspot mutations, ALK amplification, or wild-type ALK following pharmacologic inhibition of ALK using four different compounds. Next, we performed cross-species genomic analyses to identify commonly transcriptionally perturbed genes in MYCN/ALK(F1174L) double transgenic versus MYCN transgenic mouse tumors as compared with the mutant ALK-driven transcriptome in human neuroblastomas. RESULTS: A 77-gene ALK signature was established and successfully validated in primary neuroblastoma samples, in a neuroblastoma cell line with ALK(F1174L) and ALK(R1275Q) regulable overexpression constructs and in other ALKomas. In addition to the previously established PI3K/AKT/mTOR, MAPK/ERK, and MYC/MYCN signaling branches, we identified that mutant ALK drives a strong upregulation of MAPK negative feedback regulators and upregulates RET and RET-driven sympathetic neuronal markers of the cholinergic lineage. CONCLUSIONS: We provide important novel insights into the transcriptional consequences and the complexity of mutant ALK signaling in this aggressive pediatric tumor. The negative feedback loop of MAPK pathway inhibitors may affect novel ALK inhibition therapies, whereas mutant ALK induced RET signaling can offer novel opportunities for testing ALK-RET oriented molecular combination therapies.

Effective Alu repeat based RT-Qpcr normalization in cancer cell perturbation experiments.

BACKGROUND: Measuring messenger RNA (mRNA) levels using the reverse transcription quantitative polymerase chain reaction (RT-qPCR) is common practice in many laboratories. A specific set of mRNAs as internal control reference genes is considered as the preferred strategy to normalize RT-qPCR data. Proper selection of reference genes is a critical issue, especially in cancer cells that are subjected to different in vitro manipulations. These manipulations may result in dramatic alterations in gene expression levels, even of assumed reference genes. In this study, we evaluated the expression levels of 11 commonly used reference genes as internal controls for normalization of 19 experiments that include neuroblastoma, T-ALL, melanoma, breast cancer, non small cell lung cancer (NSCL), acute myeloid leukemia (AML), prostate cancer, colorectal cancer, and cervical cancer cell lines subjected to various perturbations. RESULTS: The geNorm algorithm in the software package qbase+ was used to rank the candidate reference genes according to their expression stability. We observed that the stability of most of the candidate reference genes varies greatly in perturbation experiments. Expressed Alu repeats show relatively stable expression regardless of experimental condition. These Alu repeats are ranked among the best reference assays in all perturbation experiments and display acceptable average expression stability values (M<0.5). CONCLUSIONS: We propose the use of Alu repeats as a reference assay when performing cancer cell perturbation experiments.

Focal DNA copy number changes in neuroblastoma target MYCN regulated genes.

Neuroblastoma is an embryonic tumor arising from immature sympathetic nervous system cells. Recurrent genomic alterations include MYCN and ALK amplification as well as recurrent patterns of gains and losses of whole or large partial chromosome segments. A recent whole genome sequencing effort yielded no frequently recurring mutations in genes other than those affecting ALK. However, the study further stresses the importance of DNA copy number alterations in this disease, in particular for genes implicated in neuritogenesis. Here we provide additional evidence for the importance of focal DNA copy number gains and losses, which are predominantly observed in MYCN amplified tumors. A focal 5 kb gain encompassing the MYCN regulated miR-17~92 cluster as sole gene was detected in a neuroblastoma cell line and further analyses of the array CGH data set demonstrated enrichment for other MYCN target genes in focal gains and amplifications. Next we applied an integrated genomics analysis to prioritize MYCN down regulated genes mediated by MYCN driven miRNAs within regions of focal heterozygous or homozygous deletion. We identified RGS5, a negative regulator of G-protein signaling implicated in vascular normalization, invasion and metastasis, targeted by a focal homozygous deletion, as a new MYCN target gene, down regulated through MYCN activated miRNAs. In addition, we expand the miR-17~92 regulatory network controlling TGFß signaling in neuroblastoma with the ring finger protein 11 encoding gene RNF11, which was previously shown to be targeted by the miR-17~92 member miR-19b. Taken together, our data indicate that focal DNA copy number imbalances in neuroblastoma (1) target genes that are implicated in MYCN signaling, possibly selected to reinforce MYCN oncogene addiction and (2) serve as a resource for identifying new molecular targets for treatment.

N-cadherin in neuroblastoma disease: expression and clinical significance.

One of the first and most important steps in the metastatic cascade is the loss of cell-cell and cell-matrix interactions. N-cadherin, a crucial mediator of homotypic and heterotypic cell-cell interactions, might play a central role in the metastasis of neuroblastoma (NB), a solid tumor of neuroectodermal origin. Using Reverse Transcription Quantitative PCR (RT-qPCR), Western blot, immunocytochemistry and Tissue MicroArrays (TMA) we demonstrate the expression of N-cadherin in neuroblastoma tumors and cell lines. All neuroblastic tumors (n = 356) and cell lines (n = 10) expressed various levels of the adhesion protein. The N-cadherin mRNA expression was significantly lower in tumor samples from patients suffering metastatic disease. Treatment of NB cell lines with the N-cadherin blocking peptide ADH-1 (Exherin, Adherex Technologies Inc.), strongly inhibited tumor cell proliferation in vitro by inducing apoptosis. Our results suggest that N-cadherin signaling may play a role in neuroblastoma disease, marking involvement of metastasis and determining neuroblastoma cell viability.

Dickkopf-3 is regulated by the MYCN-induced miR-17-92 cluster in neuroblastoma.

Neuroblastoma (NB) is a paediatric tumour with a remarkable diverse clinical behaviour. Approximately half of the high stage aggressive tumours are characterized by MYCN gene amplification but our understanding of the role of MYCN in NB oncogenesis is incomplete. Previous studies have shown that MYCN expression is inversely correlated with expression of Dickkopf-3 (DKK3), a gene encoding an extracellular protein with presumed tumour suppressor activity, but direct MYCN regulation of DKK3 was excluded leaving the mechanism of regulation unexplained. Given the recently established role of MYCN-regulated miRNAs in downregulation of protein-coding genes and predicted seeds for miR-17-92 cluster members within the DKK3 3'UTR, we hypothesized that this mechanism would act in MYCN regulation of DKK3. To investigate this, we used a validated miR-17-92-inducible cellular system and could demonstrate robust downregulation of DKK3 mRNA and protein levels upon miR-17-92 overexpression. Next, two of the three predicted miRNAs, miR-19b and miR-92a, were shown to lower DKK3 protein levels, in addition to measurable DKK3 mRNA knock-down by miR-92a. Direct interaction between miR-19b or miR-92a and the 3'UTR of DKK3 was validated using luciferase reporter assays. In conclusion, this study demonstrates that the MYCN-induced downregulation of DKK3 results from direct upregulation of miR-17-92 components effecting both DKK3 mRNA stability and translation which further contributes to the pleiotropic oncogenic effect of elevated MYCN levels. The strict MYCN-mediated regulation of DKK3 is suggestive for an important downstream function of the MYCN protein and thus warrants further investigations to unravel the role of DKK3 in NB.

Meta-analysis of neuroblastomas reveals a skewed ALK mutation spectrum in tumors with MYCN amplification.

PURPOSE: Activating mutations of the anaplastic lymphoma kinase (ALK) were recently described in neuroblastoma. We carried out a meta-analysis of 709 neuroblastoma tumors to determine their frequency and mutation spectrum in relation to genomic and clinical parameters, and studied the prognostic significance of ALK copy number and expression. EXPERIMENTAL DESIGN: The frequency and type of ALK mutations, copy number gain, and expression were analyzed in a new series of 254 neuroblastoma tumors. Data from 455 published cases were used for further in-depth analysis. RESULTS: ALK mutations were present in 6.9% of 709 investigated tumors, and mutations were found in similar frequencies in favorable [International Neuroblastoma Staging System (INSS) 1, 2, and 4S; 5.7%] and unfavorable (INSS 3 and 4; 7.5%) neuroblastomas (P = 0.087). Two hotspot mutations, at positions R1275 and F1174, were observed (49% and 34.7% of the mutated cases, respectively). Interestingly, the F1174 mutations occurred in a high proportion of MYCN-amplified cases (P = 0.001), and this combined occurrence was associated with a particular poor outcome, suggesting a positive cooperative effect between both aberrations. Furthermore, the F1174L mutant was characterized by a higher degree of autophosphorylation and a more potent transforming capacity as compared with the R1275Q mutant. Chromosome 2p gains, including the ALK locus (91.8%), were associated with a significantly increased ALK expression, which was also correlated with poor survival. CONCLUSIONS: ALK mutations occur in equal frequencies across all genomic subtypes, but F1174L mutants are observed in a higher frequency of MYCN-amplified tumors and show increased transforming capacity as compared with the R1275Q mutants.

Meta-mining of neuroblastoma and neuroblast gene expression profiles reveals candidate therapeutic compounds.

PURPOSE: Neuroblastoma is a heterogeneous childhood tumor with poor survival outcome for the aggressive type despite intensive multimodal therapies. In this study, we aimed to identify new treatment options for neuroblastoma based on integrative genomic analysis. EXPERIMENTAL DESIGN: The Connectivity Map is a database comprising expression profiles in response to known therapeutic compounds. This renders it a useful tool in the search for potential therapeutic compounds based on comparison of gene expression profiles of diseased cells and a database of profiles in response to known therapeutic compounds. We have used this strategy in the search for new therapeutic molecules for neuroblastoma based on data of an integrative meta-analysis of gene copy number and expression profiles from 146 primary neuroblastoma tumors and normal fetal neuroblasts. RESULTS: In a first step, a 132-gene classifier was established that discriminates three major genomic neuroblastoma subgroups, reflecting inherent differences in gene expression between these subgroups. Subsequently, we screened the Connectivity Map database using gene lists generated by comparing expression profiles of fetal adrenal neuroblasts and the genomic subgroups of neuroblastomas. A putative therapeutic effect was predicted for several compounds of which six were empirically tested. A significant reduction in cell viability was shown for five of these molecules: 17-allylamino-geldanamycin, monorden, fluphenazine, trichostatin, and rapamycin. CONCLUSIONS: This proof-of-principle study indicates that an integrative genomic meta-analysis approach with inclusion of neuroblast data enables the identification of promising compounds for treatment of children with neuroblastoma. Further studies are warranted to explore in detail the therapeutic potential of these compounds.

Low-cost dedicated mini-arrays for high-throughput analysis of DNA copy-number alterations in neuroblastoma.

ArrayCGH is commonly used for high-resolution detection of copy-number alterations in tumours, allowing identification of chromosomal aberrations with prognostic or diagnostic relevance. Currently available arrayCGH platforms are still very expensive for analysis of large sets of samples. For this purpose, we have constructed a dedicated mini-array that is enriched for BAC/PAC clones in the prognostic important regions for neuroblastoma and that only covers a small area on the slide, allowing down-scaling of the labelling and hybridisation reagents and hence reducing the price. The mini-arrays were validated on neuroblastoma samples and comparison with high-resolution whole-genome arrayCGH data yielded complete concordant results.