Identifying Patients with Polycythemia Vera at Risk of Thrombosis after Hydroxyurea Initiation: The Polycythemia Vera-Advanced Integrated Models (PV-AIM) Project.

Patients with polycythemia vera (PV) are at significant risk of thromboembolic events (TE). The PV-AIM study used the Optum® de-identified Electronic Health Record dataset and machine learning to identify markers of TE in a real-world population. Data for 82,960 patients with PV were extracted: 3852 patients were treated with hydroxyurea (HU) only, while 130 patients were treated with HU and then changed to ruxolitinib (HU-ruxolitinib). For HU-alone patients, the annualized incidence rates (IR; per 100 patients) decreased from 8.7 (before HU) to 5.6 (during HU) but increased markedly to 10.5 (continuing HU). Whereas for HU-ruxolitinib patients, the IR decreased from 10.8 (before HU) to 8.4 (during HU) and was maintained at 8.3 (after switching to ruxolitinib). To better understand markers associated with TE risk, we built a machine-learning model for HU-alone patients and validated it using an independent dataset. The model identified lymphocyte percentage (LYP), neutrophil percentage (NEP), and red cell distribution width (RDW) as key markers of TE risk, and optimal thresholds for these markers were established, from which a decision tree was derived. Using these widely used laboratory markers, the decision tree could be used to identify patients at high risk for TE, facilitate treatment decisions, and optimize patient management.

Extensive rewiring of the EGFR network in colorectal cancer cells expressing transforming levels of KRASG13D.

Protein-protein-interaction networks (PPINs) organize fundamental biological processes, but how oncogenic mutations impact these interactions and their functions at a network-level scale is poorly understood. Here, we analyze how a common oncogenic KRAS mutation (KRASG13D) affects PPIN structure and function of the Epidermal Growth Factor Receptor (EGFR) network in colorectal cancer (CRC) cells. Mapping >6000 PPIs shows that this network is extensively rewired in cells expressing transforming levels of KRASG13D (mtKRAS). The factors driving PPIN rewiring are multifactorial including changes in protein expression and phosphorylation. Mathematical modelling also suggests that the binding dynamics of low and high affinity KRAS interactors contribute to rewiring. PPIN rewiring substantially alters the composition of protein complexes, signal flow, transcriptional regulation, and cellular phenotype. These changes are validated by targeted and global experimental analysis. Importantly, genetic alterations in the most extensively rewired PPIN nodes occur frequently in CRC and are prognostic of poor patient outcomes.

Contextual Hub Analysis Tool (CHAT): A Cytoscape app for identifying contextually relevant hubs in biological networks.

Highly connected nodes (hubs) in biological networks are topologically important to the structure of the network and have also been shown to be preferentially associated with a range of phenotypes of interest. The relative importance of a hub node, however, can change depending on the biological context. Here, we report a Cytoscape app, the Contextual Hub Analysis Tool (CHAT), which enables users to easily construct and visualize a network of interactions from a gene list of interest, integrate contextual information, such as gene expression data, and identify hub nodes that are more highly connected to contextual nodes (e.g. genes that are differentially expressed) than expected by chance. In a case study, we use CHAT to construct a network of genes that are differentially expressed in Dengue fever, a viral infection. CHAT was used to identify and compare contextual and degree-based hubs in this network. The top 20 degree-based hubs were enriched in pathways related to the cell cycle and cancer, which is likely due to the fact that proteins involved in these processes tend to be highly connected in general. In comparison, the top 20 contextual hubs were enriched in pathways commonly observed in a viral infection including pathways related to the immune response to viral infection. This analysis shows that such contextual hubs are considerably more biologically relevant than degree-based hubs and that analyses which rely on the identification of hubs solely based on their connectivity may be biased towards nodes that are highly connected in general rather than in the specific context of interest. AVAILABILITY: CHAT is available for Cytoscape 3.0+ and can be installed via the Cytoscape App Store (http://apps.cytoscape.org/apps/chat).

MicroRNA profiling of the bovine alveolar macrophage response to Mycobacterium bovis infection suggests pathogen survival is enhanced by microRNA regulation of endocytosis and lysosome trafficking.

Mycobacterium bovis, the causative agent of bovine tuberculosis, a major problem for global agriculture, spreads via an airborne route and is taken up by alveolar macrophages (AM) in the lung. Here, we describe the first next-generation sequencing (RNA-seq) approach to temporally profile miRNA expression in primary bovine AMs post-infection with M. bovis. One, six, and forty miRNAs were identified as significantly differentially expressed at 2, 24 and 48 h post-infection, respectively. The differential expression of three miRNAs (bta-miR-142-5p, bta-miR-146a, and bta-miR-423-3p) was confirmed by RT-qPCR. Pathway analysis of the predicted mRNA targets of differentially expressed miRNAs suggests that these miRNAs preferentially target several pathways that are functionally relevant for mycobacterial pathogenesis, including endocytosis and lysosome trafficking, IL-1 signalling and the TGF-ß pathway. Over-expression studies using a bovine macrophage cell-line (Bomac) reveal the targeting of two key genes in the innate immune response to M. bovis, IL-1 receptor-associated kinase 1 (IRAK1) and TGF-ß receptor 2 (TGFBR2), by miR-146. Taken together, our study suggests that miRNAs play a key role in tuning the complex interplay between M. bovis survival strategies and the host immune response.

Discovery and visualization of miRNA-mRNA functional modules within integrated data using bicluster analysis.

MicroRNAs (miRNAs) are small non-coding RNA molecules that regulate gene expression at a post-transcriptional level. An miRNA may target many messenger RNA (mRNA) transcripts, and each transcript may be targeted by multiple miRNAs. Our understanding of miRNA regulation is evolving to consider modules of miRNAs that regulate groups of functionally related mRNAs. Here we expand the model of miRNA functional modules and use it to guide the integration of miRNA and mRNA expression and target prediction data. We present evidence of cooperativity between miRNA classes within this integrated miRNA-mRNA association matrix. We then apply bicluster analysis to uncover miRNA functional modules within this integrated data set and develop a novel application to visualize and query these results. We show that this wholly unsupervised approach can discover a network of miRNA-mRNA modules that are enriched for both biological processes and miRNA classes. We apply this method to investigate the interplay of miRNAs and mRNAs in integrated data sets derived from neuroblastoma and human immune cells. This study is the first to apply the technique of biclustering to model functional modules within an integrated miRNA-mRNA association matrix. Results provide evidence of an extensive modular miRNA functional network and enable characterization of miRNA function and dysregulation in disease.

Post-transcriptional dysregulation by miRNAs is implicated in the pathogenesis of gastrointestinal stromal tumor [GIST].

In contrast to adult mutant gastrointestinal stromal tumors [GISTs], pediatric/wild-type GISTs remain poorly understood overall, given their lack of oncogenic activating tyrosine kinase mutations. These GISTs, with a predilection for gastric origin in female patients, show limited response to therapy with tyrosine kinase inhibitors and generally pursue a more indolent course, but still may prove fatal. Defective cellular respiration appears to underpin tumor development in these wild-type cases, which as a group lack expression of succinate dehydrogenase [SDH] B, a surrogate marker for respiratory chain metabolism. Yet, only a small subset of the wild-type tumors show mutations in the genes coding for the SDH subunits [SDHx]. To explore additional pathogenetic mechanisms in these wild-type GISTs, we elected to investigate post-transcriptional regulation of these tumors by conducting microRNA (miRNA) profiling of a mixed cohort of 73 cases including 18 gastric pediatric wild-type, 25 (20 gastric, 4 small bowel and 1 retroperitoneal) adult wild-type GISTs and 30 gastric adult mutant GISTs. By this approach we have identified distinct signatures for GIST subtypes which correlate tightly with clinico-pathological parameters. A cluster of miRNAs on 14q32 show strikingly different expression patterns amongst GISTs, a finding which appears to be explained at least in part by differential allelic methylation of this imprinted region. Small bowel and retroperitoneal wild-type GISTs segregate with adult mutant GISTs and express SDHB, while adult wild-type gastric GISTs are dispersed amongst adult mutant and pediatric wild-type cases, clustering in this situation on the basis of SDHB expression. Interestingly, global methylation analysis has recently similarly demonstrated that these wild-type, SDHB-immunonegative tumors show a distinct pattern compared with KIT and PDGFRA mutant tumors, which as a rule do express SDHB. All cases with Carney triad within our cohort cluster together tightly.

Expressional alterations in functional ultra-conserved non-coding RNAs in response to all-trans retinoic acid--induced differentiation in neuroblastoma cells.

BACKGROUND: Ultra-conserved regions (UCRs) are segments of the genome (≥ 200 bp) that exhibit 100% DNA sequence conservation between human, mouse and rat. Transcribed UCRs (T-UCRs) have been shown to be differentially expressed in cancers versus normal tissue, indicating a possible role in carcinogenesis. All-trans-retinoic acid (ATRA) causes some neuroblastoma (NB) cell lines to undergo differentiation and leads to a significant decrease in the oncogenic transcription factor MYCN. Here, we examine the impact of ATRA treatment on T-UCR expression and investigate the biological significance of these changes. METHODS: We designed a custom tiling microarray to profile the expression of 481 T-UCRs in sense and anti-sense orientation (962 potential transcripts) in untreated and ATRA-treated neuroblastoma cell lines (SH-SY5Y, SK-N-BE, LAN-5). Following identification of significantly differentially expressed T-UCRs, we carried out siRNA knockdown and gene expression microarray analysis to investigate putative functional roles for selected T-UCRs. RESULTS: Following ATRA-induced differentiation, 32 T-UCRs were differentially expressed (16 up-regulated, 16 down-regulated) across all three cell lines. Further insight into the possible role of T-UC.300A, an independent transcript whose expression is down-regulated following ATRA was achieved by siRNA knockdown, resulting in the decreased viability and invasiveness of ATRA-responsive cell lines. Gene expression microarray analysis following knockdown of T-UC.300A revealed a number of genes whose expression was altered by changing T-UC.300A levels and that might play a role in the increased proliferation and invasion of NB cells prior to ATRA-treatment. CONCLUSIONS: Our results indicate that significant numbers of T-UCRs have altered expression levels in response to ATRA. While the precise roles that T-UCRs might play in cancer or in normal development are largely unknown and an important area for future study, our findings strongly indicate that the function of non-coding RNA T-UC.300A is connected with proliferation, invasion and the inhibition of differentiation of neuroblastoma cell lines prior to ATRA treatment.

MiR-137 functions as a tumor suppressor in neuroblastoma by downregulating KDM1A.

Neuroblastoma is the most common extracranial solid tumor of childhood, and accounts for ∼15% of all childhood cancer deaths. The histone demethylase, lysine-specific demethylase 1 (KDM1A, previously known as LSD1), is strongly expressed in neuroblastomas, and overexpression correlates with poor patient prognosis. Inducing differentiation in neuroblastoma cells has previously been shown to down regulate KDM1A, and siRNA-mediated KDM1A knockdown inhibited neuroblastoma cell viability. The microRNA, miR-137, has been reported to be downregulated in several human cancers, and KDM1A mRNA was reported as a putative target of miR-137 in colon cancer. We hypothesized that miR-137 might have a tumor-suppressive role in neuroblastoma mediated via downregulation of KDM1A. Indeed, low levels of miR-137 expression in primary neuroblastomas correlated with poor patient prognosis. Re-expressing miR-137 in neuroblastoma cell lines increased apoptosis and decreased cell viability and proliferation. KDM1A mRNA was repressed by miR-137 in neuroblastoma cells, and was validated as a direct target of miR-137 using reporter assays in SHEP and HEK293 cells. Furthermore, siRNA-mediated KDM1A knockdown phenocopied the miR-137 re-expression phenotype in neuroblastoma cells. We conclude that miR-137 directly targets KDM1A mRNA in neuroblastoma cells, and activates cell properties consistent with tumor suppression. Therapeutic strategies to re-express miR-137 in neuroblastomas could be useful to reduce tumor aggressiveness.

miRNA profiles as a predictor of chemoresponsiveness in Wilms' tumor blastema.

The current SIOP treatment protocol for Wilms' tumor involves pre-operative chemotherapy followed by nephrectomy. Not all patients benefit equally from such chemotherapy. The aim of this study was to generate a miRNA profile of chemo resistant blastemal cells in high risk Wilms' tumors which might serve as predictive markers of therapeutic response at the pre-treatment biopsy stage. We have shown here that unsupervised hierarchical clustering of genome-wide miRNA expression profiles can clearly separate intermediate risk tumors from high risk tumors. A total of 29 miRNAs were significantly differentially expressed between post-treatment intermediate risk and high risk groups, including miRNAs that have been previously linked to chemo resistance in other cancer types. Furthermore, 7 of these 29 miRNAs were already at the pre-treatment biopsy stage differentially expressed between cases ultimately deemed intermediate risk compared to high risk. These miRNA alterations include down-regulation in high risk cases of miR-193a.5p, miR-27a and the up-regulation of miR-483.5p, miR-628.5p, miR-590.5p, miR-302a and miR-367. The demonstration of such miRNA markers at the pre-treatment biopsy stage could permit stratification of patients to more tailored treatment regimens.

MiRNA-335 suppresses neuroblastoma cell invasiveness by direct targeting of multiple genes from the non-canonical TGF-ß signalling pathway.

Transforming growth factor-ß (TGF-ß) signaling regulates many diverse cellular activities through both canonical (SMAD-dependent) and non-canonical branches, which includes the mitogen-activated protein kinase (MAPK), Rho-like guanosine triphosphatase and phosphatidylinositol-3-kinase/AKT pathways. Here, we demonstrate that miR-335 directly targets and downregulates genes in the TGF-ß non-canonical pathways, including the Rho-associated coiled-coil containing protein (ROCK1) and MAPK1, resulting in reduced phosphorylation of downstream pathway members. Specifically, inhibition of ROCK1 and MAPK1 reduces phosphorylation levels of the motor protein myosin light chain (MLC) leading to a significant inhibition of the invasive and migratory potential of neuroblastoma cells. Additionally, miR-335 targets the leucine-rich alpha-2-glycoprotein 1 (LRG1) messenger RNA, which similarly results in a significant reduction in the phosphorylation status of MLC and a decrease in neuroblastoma cell migration and invasion. Thus, we link LRG1 to the migratory machinery of the cell, altering its activity presumably by exerting its effect within the non-canonical TGF-ß pathway. Moreover, we demonstrate that the MYCN transcription factor, whose coding sequence is highly amplified in a particularly clinically aggressive neuroblastoma tumor subtype, directly binds to a region immediately upstream of the miR-335 transcriptional start site, resulting in transcriptional repression. We conclude that MYCN contributes to neuroblastoma cell migration and invasion, by directly downregulating miR-335, resulting in the upregulation of the TGF-ß signaling pathway members ROCK1, MAPK1 and putative member LRG1, which positively promote this process. Our results provide novel insight into the direct regulation of TGF-ß non-canonical signaling by miR-335, which in turn is downregulated by MYCN.

miRNA expression profiling of the murine TH-MYCN neuroblastoma model reveals similarities with human tumors and identifies novel candidate miRNAs.

BACKGROUND: MicroRNAs are small molecules which regulate gene expression post-transcriptionally and aberrant expression of several miRNAs is associated with neuroblastoma, a childhood cancer arising from precursor cells of the sympathetic nervous system. Amplification of the MYCN transcription factor characterizes the most clinically aggressive subtype of this disease, and although alteration of p53 signaling is not commonly found in primary tumors, deregulation of proteins involved in this pathway frequently arise in recurrent disease after pharmacological treatment. TH-MYCN is a well-characterized transgenic model of MYCN-driven neuroblastoma which recapitulates many clinicopathologic features of the human disease. Here, we evaluate the dysregulation of miRNAs in tumors from TH-MYCN mice that are either wild-type (TH-MYCN) or deficient (TH-MYCN/p53ER(TAM)) for the p53 tumor suppressor gene. PRINCIPAL FINDINGS: We analyzed the expression of 591 miRNAs in control (adrenal) and neuroblastoma tumor tissues derived from either TH-MYCN or TH-MYCN/p53ER(TAM) mice, respectively wild-type or deficient in p53. Comparing miRNA expression in tumor and control samples, we identified 159 differentially expressed miRNAs. Using data previously obtained from human neuroblastoma samples, we performed a comparison of miRNA expression between murine and human tumors to assess the concordance between murine and human expression data. Notably, the miR-17-5p-92 oncogenic polycistronic cluster, which is over-expressed in human MYCN amplified tumors, was over-expressed in mouse tumors. Moreover, analyzing miRNAs expression in a mouse model (TH-MYCN/p53ER(TAM)) possessing a transgenic p53 allele that drives the expression of an inactive protein, we identified miR-125b-3p and miR-676 as directly or indirectly regulated by the level of functional p53. SIGNIFICANCE: Our study represents the first miRNA profiling of an important mouse model of neuroblastoma. Similarities and differences in miRNAs expression between human and murine neuroblastoma were identified, providing important insight into the efficacy of this mouse model for assessing miRNA involvement in neuroblastoma and their potential effectiveness as therapeutic targets.

miRNA expression profiling enables risk stratification in archived and fresh neuroblastoma tumor samples.

PURPOSE: More accurate assessment of prognosis is important to further improve the choice of risk-related therapy in neuroblastoma (NB) patients. In this study, we aimed to establish and validate a prognostic miRNA signature for children with NB and tested it in both fresh frozen and archived formalin-fixed paraffin-embedded (FFPE) samples. EXPERIMENTAL DESIGN: Four hundred-thirty human mature miRNAs were profiled in two patient subgroups with maximally divergent clinical courses. Univariate logistic regression analysis was used to select miRNAs correlating with NB patient survival. A 25-miRNA gene signature was built using 51 training samples, tested on 179 test samples, and validated on an independent set of 304 fresh frozen tumor samples and 75 archived FFPE samples. RESULTS: The 25-miRNA signature significantly discriminates the test patients with respect to progression-free and overall survival (P < 0.0001), both in the overall population and in the cohort of high-risk patients. Multivariate analysis indicates that the miRNA signature is an independent predictor of patient survival after controlling for current risk factors. The results were confirmed in an external validation set. In contrast to a previously published mRNA classifier, the 25-miRNA signature was found to be predictive for patient survival in a set of 75 FFPE neuroblastoma samples. CONCLUSIONS: In this study, we present the largest NB miRNA expression study so far, including more than 500 NB patients. We established and validated a robust miRNA classifier, able to identify a cohort of high-risk NB patients at greater risk for adverse outcome using both fresh frozen and archived material.

Co-localization of the oncogenic transcription factor MYCN and the DNA methyl binding protein MeCP2 at genomic sites in neuroblastoma.

BACKGROUND: MYCN is a transcription factor that is expressed during the development of the neural crest and its dysregulation plays a major role in the pathogenesis of pediatric cancers such as neuroblastoma, medulloblastoma and rhabdomyosarcoma. MeCP2 is a CpG methyl binding protein which has been associated with a number of cancers and developmental disorders, particularly Rett syndrome. METHODS AND FINDINGS: Using an integrative global genomics approach involving chromatin immunoprecipitation applied to microarrays, we have determined that MYCN and MeCP2 co-localize to gene promoter regions, as well as inter/intragenic sites, within the neuroblastoma genome (MYCN amplified Kelly cells) at high frequency (70.2% of MYCN sites were also positive for MeCP2). Intriguingly, the frequency of co-localization was significantly less at promoter regions exhibiting substantial hypermethylation (8.7%), as determined by methylated DNA immunoprecipitation (MeDIP) applied to the same microarrays. Co-immunoprecipitation of MYCN using an anti-MeCP2 antibody indicated that a MYCN/MeCP2 interaction occurs at protein level. mRNA expression profiling revealed that the median expression of genes with promoters bound by MYCN was significantly higher than for genes bound by MeCP2, and that genes bound by both proteins had intermediate expression. Pathway analysis was carried out for genes bound by MYCN, MeCP2 or MYCN/MeCP2, revealing higher order functions. CONCLUSIONS: Our results indicate that MYCN and MeCP2 protein interact and co-localize to similar genomic sites at very high frequency, and that the patterns of binding of these proteins can be associated with significant differences in transcriptional activity. Although it is not yet known if this interaction contributes to neuroblastoma disease pathogenesis, it is intriguing that the interaction occurs at the promoter regions of several genes important for the development of neuroblastoma, including ALK, AURKA and BDNF.

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.

MicroRNA-542-5p as a novel tumor suppressor in neuroblastoma.

Several studies have implicated the dysregulation of microRNAs in neuroblastoma pathogenesis, an often fatal paediatric cancer arising from precursor cells of the sympathetic nervous system. Our group and others have demonstrated that lower expression of miR-542-5p is highly associated with poor patient survival, indicating a potential tumor suppressive function. Here, we demonstrate that ectopic over-expression of this miRNA decreases the invasive potential of neuroblastoma cell lines in vitro, along with primary tumor growth and metastases in an orthotopic mouse xenograft model, providing the first functional evidence for the involvement of miR-542-5p as a tumor suppressor in any type of cancer.

Genome-wide DNA methylation analysis of neuroblastic tumors reveals clinically relevant epigenetic events and large-scale epigenomic alterations localized to telomeric regions.

The downregulation of specific genes through DNA hypermethylation is a major hallmark of cancer, although the extent and genomic distribution of hypermethylation occurring within cancer genomes is poorly understood. We report on the first genome-wide analysis of DNA methylation alterations in different neuroblastic tumor subtypes and cell lines, revealing higher order organization and clinically relevant alterations of the epigenome. The methylation status of 33,485 discrete loci representing all annotated CpG islands and RefSeq gene promoters was assessed in primary neuroblastic tumors and cell lines. A comparison of genes that were hypermethylated exclusively in the clinically favorable ganglioneuroma/ganglioneuroblastoma tumors revealed that nine genes were associated with poor clinical outcome when overexpressed in the unfavorable neuroblastoma (NB) tumors. Moreover, an integrated DNA methylation and copy number analysis identified 80 genes that were recurrently concomitantly deleted and hypermethylated in NB, with 37 reactivated by 5-aza-deoxycytidine. Lower expression of four of these genes was correlated with poor clinical outcome, further implicating their inactivation in aggressive disease pathogenesis. Analysis of genome-wide hypermethylation patterns revealed 70 recurrent large-scale blocks of contiguously hypermethylated promoters/CpG islands, up to 590 kb in length, with a distribution bias toward telomeric regions. Genome-wide hypermethylation events in neuroblastic tumors are extensive and frequently occur in large-scale blocks with a significant bias toward telomeric regions, indicating that some methylation alterations have occurred in a coordinated manner. Our results indicate that methylation contributes toward the clinicopathological features of neuroblastic tumors, revealing numerous genes associated with poor patient survival in NB.

MicroRNA mediates DNA demethylation events triggered by retinoic acid during neuroblastoma cell differentiation.

Neuroblastoma is an often fatal pediatric cancer arising from precursor cells of the sympathetic nervous system. 13-Cis retinoic acid is included in the treatment regimen for patients with high-risk disease, and a similar derivative, all-trans-retinoic acid (ATRA), causes neuroblastoma cell lines to undergo differentiation. The molecular signaling pathways involved with ATRA-induced differentiation are complex, and the role that DNA methylation changes might play are unknown. The purpose of this study was to evaluate the genome-wide effects of ATRA on DNA methylation using methylated DNA immunoprecipitation applied to microarrays representing all known promoter and CpG islands. Four hundred and two gene promoters became demethylated, whereas 88 were hypermethylated post-ATRA. mRNA expression microarrays revealed that 82 of the demethylated genes were overexpressed by >2-fold, whereas 13 of the hypermethylated genes were underexpressed. Gene ontology analysis indicated that demethylated and re-expressed genes were enriched for signal transduction pathways, including NOS1, which is required for neural cell differentiation. As a potential mechanism for the DNA methylation changes, we show the downregulation of methyltransferases, DNMT1 and DNMT3B, along with the upregulation of endogenous microRNAs targeting them. Ectopic overexpression of miR-152, targeting DNMT1, also negatively affected cell invasiveness and anchorage-independent growth, contributing in part to the differentiated phenotype. We conclude that functionally important, miRNA-mediated DNA demethylation changes contribute to the process of ATRA-induced differentiation resulting in the activation of NOS1, a critical determinant of neural cell differentiation. Our findings illustrate the plasticity and dynamic nature of the epigenome during cancer cell differentiation.

Therapeutic targeting of miRNAs in neuroblastoma.

IMPORTANCE OF THE FIELD: Neuroblastomas arise from precursor cells of the sympathetic nervous system and are noted for highly heterogeneous clinical behavior. These tumors currently account for approximately 15% of all childhood cancer related deaths in spite of intensive multimodal chemotherapy and are a major problem in pediatric oncology. The identification of novel therapeutic targets is urgently required to reduce patient morbidity. AREAS COVERED IN THIS REVIEW: The purpose of this article is to review and synthesize all of the rapidly expanding evidence for the contribution of microRNAs (miRNAs) in neuroblastoma aggressive disease pathogenesis, along with the prospect of using small RNAs as therapeutics. WHAT THE READER WILL GAIN: The reader will obtain insight on the miRNAs that are dysregulated in neuroblastoma along with potential therapeutic strategies and the most promising targets. TAKE HOME MESSAGE: A number of miRNAs which are associated with aggressive disease pathogenesis in neuroblastoma patients have been demonstrated to contribute in major ways to cell proliferation rates, apoptosis, differentiation, invasiveness and tumor growth in vitro and in vivo. Directly or indirectly interfering with the function of these miRNAs may prove to be an important and novel form of therapy.

MicroRNA-184 inhibits neuroblastoma cell survival through targeting the serine/threonine kinase AKT2.

BACKGROUND: Neuroblastoma is a paediatric cancer of the sympathetic nervous system. The single most important genetic indicator of poor clinical outcome is amplification of the MYCN transcription factor. One of many down-stream MYCN targets is miR-184, which is either directly or indirectly repressed by this transcription factor, possibly due to its pro-apoptotic effects when ectopically over-expressed in neuroblastoma cells. The purpose of this study was to elucidate the molecular mechanism by which miR-184 conveys pro-apoptotic effects. RESULTS: We demonstrate that the knock-down of endogenous miR-184 has the opposite effect of ectopic up-regulation, leading to enhanced neuroblastoma cell numbers. As a mechanism of how miR-184 causes apoptosis when over-expressed, and increased cell numbers when inhibited, we demonstrate direct targeting and degradation of AKT2, a major downstream effector of the phosphatidylinositol 3-kinase (PI3K) pathway, one of the most potent pro-survival pathways in cancer. The pro-apoptotic effects of miR-184 ectopic over-expression in neuroblastoma cell lines is reproduced by siRNA inhibition of AKT2, while a positive effect on cell numbers similar to that obtained by the knock-down of endogenous miR-184 can be achieved by ectopic up-regulation of AKT2. Moreover, co-transfection of miR-184 with an AKT2 expression vector lacking the miR-184 target site in the 3'UTR rescues cells from the pro-apoptotic effects of miR-184. CONCLUSIONS: MYCN contributes to tumorigenesis, in part, by repressing miR-184, leading to increased levels of AKT2, a direct target of miR-184. Thus, two important genes with positive effects on cell growth and survival, MYCN and AKT2, can be linked into a common genetic pathway through the actions of miR-184. As an inhibitor of AKT2, miR-184 could be of potential benefit in miRNA mediated therapeutics of MYCN amplified neuroblastoma and other forms of cancer.

Chromosomal and microRNA expression patterns reveal biologically distinct subgroups of 11q- neuroblastoma.

PURPOSE: The purpose of this study was to further define the biology of the 11q- neuroblastoma tumor subgroup by the integration of array-based comparative genomic hybridization with microRNA (miRNA) expression profiling data to determine if improved patient stratification is possible. EXPERIMENTAL DESIGN: A set of primary neuroblastoma (n = 160), which was broadly representative of all genetic subtypes, was analyzed by array-based comparative genomic hybridization and for the expression of 430 miRNAs. A 15-miRNA expression signature previously shown to be predictive of clinical outcome was used to analyze an independent cohort of 11q- tumors (n = 37). RESULTS: Loss of 4p and gain of 7q occurred at a significantly higher frequency in the 11q- tumors, further defining the genetic characteristics of this subtype. The 11q- tumors could be split into two subgroups using a miRNA expression survival signature that differed significantly in clinical outcome and the overall frequency of large-scale genomic imbalances, with the poor survival subgroup having significantly more imbalances. miRNAs from the expression signature, which were upregulated in unfavorable tumors, were predicted to target downregulated genes from a published mRNA expression classifier of clinical outcome at a higher-than-expected frequency, indicating the miRNAs might contribute to the regulation of genes within the signature. CONCLUSION: We show that two distinct biological subtypes of neuroblastoma with loss of 11q occur, which differ in their miRNA expression profiles, frequency of segmental imbalances, and clinical outcome. A miRNA expression signature, combined with an analysis of segmental imbalances, provides greater prediction of event-free survival and overall survival outcomes than 11q status by itself, improving patient stratification.