Karen M. Watters

MicroRNAs 10a and 10b are potent inducers of neuroblastoma cell differentiation through targeting of nuclear receptor corepressor 2

MicroRNAs function as negative regulators of posttranscriptional gene expression, having major roles in cellular differentiation. Several neuroblastoma cell lines can be induced to undergo differentiation by all-trans-retinoic acid (ATRA) and are used for modeling signaling pathways involved in this process. To identify miRNAs contributing to differentiation, we profiled 364 loci following ATRA treatment of neuroblastoma cell lines and found miR-10a and miR-10b to be highly overexpressed in SK-N-BE, LAN5 and SHSY-5Y. Ectopic overexpression of these miRNAs led to a major reprogramming of the transcriptome and a differentiated phenotype that was similar to that induced by ATRA in each of these cell lines. One of the predicted downregulated miR-10a/b targets was nuclear receptor corepressor 2 (NCOR2), a corepressor of gene transcription, which is known to suppress neurite outgrowth. NCOR2 was experimentally validated as a direct target of miR-10a/b, and siRNA-mediated inhibition of this mRNA alone resulted in neural cell differentiation. Moreover, induction of differentiation could be blocked by ectopic upregulation of NCOR2 using an expression construct lacking the miR-10a/b 3′ untranslated region target site. We conclude that miR-10a/b has major roles in the process of neural cell differentiation through direct targeting of NCOR2, which in turn induces a cascade of primary and secondary transcriptional alterations, including the downregulation of MYCN.

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.

MicroRNA Mediates DNA De-methylation 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.

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.

Global MYCN Transcription Factor Binding Analysis in Neuroblastoma Reveals Association with Distinct E-Box Motifs and Regions of DNA Hypermethylation

Background Neuroblastoma, a cancer derived from precursor cells of the sympathetic nervous system, is a major cause of childhood cancer related deaths. The single most important prognostic indicator of poor clinical outcome in this disease is genomic amplification of MYCN, a member of a family of oncogenic transcription factors. Methodology We applied MYCN chromatin immunoprecipitation to microarrays (ChIP-chip) using MYCN amplified/non-amplified cell lines as well as a conditional knockdown cell line to determine the distribution of MYCN binding sites within all annotated promoter regions. Conclusion Assessment of E-box usage within consistently positive MYCN binding sites revealed a predominance for the CATGTG motif (p<0.0016), with significant enrichment of additional motifs CATTTG, CATCTG, CAACTG in the MYCN amplified state. For cell lines over-expressing MYCN, gene ontology analysis revealed enrichment for the binding of MYCN at promoter regions of numerous molecular functional groups including DNA helicases and mRNA transcriptional regulation. In order to evaluate MYCN binding with respect to other genomic features, we determined the methylation status of all annotated CpG islands and promoter sequences using methylated DNA immunoprecipitation (MeDIP). The integration of MYCN ChIP-chip and MeDIP data revealed a highly significant positive correlation between MYCN binding and DNA hypermethylation. This association was also detected in regions of hemizygous loss, indicating that the observed association occurs on the same homologue. In summary, these findings suggest that MYCN binding occurs more commonly at CATGTG as opposed to the classic CACGTG E-box motif, and that disease associated over expression of MYCN leads to aberrant binding to additional weaker affinity E-box motifs in neuroblastoma. The co-localization of MYCN binding and DNA hypermethylation further supports the dual role of MYCN, namely that of a classical transcription factor affecting the activity of individual genes, and that of a mediator of global chromatin structure.

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.

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.

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.

Orchestrating osteogenic differentiation of mesenchymal stem cells—identification of placental growth factor as a mechanosensitive gene with a pro‐osteogenic role

Skeletogenesis is initiated during fetal development and persists through adult life as either a remodeling process in response to homeostatic regulation or as a regenerative process in response to physical injury. Mesenchymal stem cells (MSCs) play a crucial role providing progenitor cells from which osteoblasts, bone matrix forming cells are differentiated. The mechanical environment plays an important role in regulating stem cell differentiation into osteoblasts, however, the mechanisms by which MSCs respond to mechanical stimuli are yet to be fully elucidated. To increase understanding of MSC mechanotransuction and osteogenic differentiation, this study aimed to identify novel, mechanically augmented genes and pathways with pro‐osteogenic functionality. Using collagen glycoaminoglycan scaffolds as mimics of native extracellular matrix, to create a 3D environment more representative of that found in bone, MSC‐seeded constructs were mechanically stimulated in a flow‐perfusion bioreactor. Global gene expression profiling techniques were used to identify potential candidates warranting further investigation. Of these, placental growth factor (PGF) was selected and expression levels were shown to strongly correlate to both the magnitude and duration of mechanical stimulation. We demonstrated that PGF gene expression was modulated through an actin polymerization‐mediated mechanism. The functional role of PGF in modulating MSC osteogenic differentiation was interrogated, and we showed a concentration‐dependent response whereby low concentrations exhibited the strongest pro‐osteogenic effect. Furthermore, pre‐osteoclast migration and differentiation, as well as endothelial cell tubule formation also maintained concentration‐dependent responses to PGF, suggesting a potential role for PGF in bone resorption and angiogenesis, processes key to bone remodeling and fracture repair. Stem Cells 2013;31:2420–2431

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

BackgroundUltra-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.MethodsWe 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.ResultsFollowing 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.ConclusionsOur 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.