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.

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.

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.

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.

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.

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.

Tax 1-independent induction of vascular endothelial growth factor in adult T-cell leukemia caused by human T-cell leukemia virus type 1.

Adult T-cell leukemia (ATL) is caused by human T-cell leukemia virus type 1 (HTLV-1). Elevated expression of vascular endothelial growth factor (VEGF) in ATL patients is associated with leukemic cell invasion and infiltration in different organs. The regulatory protein Tax 1 encoded by HTLV-1 plays a pivotal role in T-cell transformation by deregulating the function and expression of several cellular factors. In the present study, we examined the effect of Tax 1 on VEGF expression at transcriptional and posttranscriptional levels in order to elucidate the regulatory mechanisms involved. Using functional assays, we demonstrate that Tax 1 downregulates the VEGF promoter through a cluster of Sp1 sites located close to the transcriptional start site. Using gel mobility shift assays, we show that Tax 1 reduced Sp1:DNA complex formation. We demonstrate that the level of secreted VEGF was significantly lower in Tax 1-transfected 293T cells compared to nontransfected cells, which is consistent with the observed downregulatory effect of Tax 1 at the transcription level. We showed that VEGF was secreted by HTLV-1-transformed and nontransformed cells, irrespective of Tax 1 expression. Overall our data indicate that, contrary to a previous report, Tax 1 downregulates VEGF expression and suggest there are Tax 1-independent mechanisms of VEGF activation in ATL.

The Natural Product ß-Escin Targets Cancer and Stromal Cells of the Tumor Microenvironment to Inhibit Ovarian Cancer Metastasis.

The high mortality of OvCa is caused by the wide dissemination of cancer within the abdominal cavity. OvCa cells metastasize to the peritoneum, which is covered by mesothelial cells, and invade into the underlying stroma, composed of extracellular matrices (ECM) and stromal cells. In a study using a three-dimensional quantitative high-throughput screening platform (3D-qHTS), we found that ß-escin, a component of horse chestnut seed extract, inhibited OvCa adhesion/invasion. Here, we determine whether ß-escin and structurally similar compounds have a therapeutic potential against OvCa metastasis. Different sources of ß-escin and horse chestnut seed extract inhibited OvCa cell adhesion/invasion, both in vitro and in vivo. From a collection of 160 structurally similar compounds to ß-escin, we found that cardiac glycosides inhibited OvCa cell adhesion/invasion and proliferation in vitro, and inhibited adhesion/invasion and metastasis in vivo. Mechanistically, ß-escin and the cardiac glycosides inhibited ECM production in mesothelial cells and fibroblasts. The oral administration of ß-escin inhibited metastasis in both OvCa prevention and intervention mouse models. Specifically, ß-escin inhibited ECM production in the omental tumors. Additionally, the production of HIF1α-targeted proteins, lactate dehydrogenase A, and hexokinase 2 in omental tumors was blocked by ß-escin. This study reveals that the natural compound ß-escin has a therapeutic potential because of its ability to prevent OvCa dissemination by targeting both cancer and stromal cells in the OvCa tumor microenvironment.

Mesothelial Cell HIF1α Expression Is Metabolically Downregulated by Metformin to Prevent Oncogenic Tumor-Stromal Crosstalk.

The tumor microenvironment (TME) plays a pivotal role in cancer progression, and, in ovarian cancer (OvCa), the primary TME is the omentum. Here, we show that the diabetes drug metformin alters mesothelial cells in the omental microenvironment. Metformin interrupts bidirectional signaling between tumor and mesothelial cells by blocking OvCa cell TGF-ß signaling and mesothelial cell production of CCL2 and IL-8. Inhibition of tumor-stromal crosstalk by metformin is caused by the reduced expression of the tricarboxylic acid (TCA) enzyme succinyl CoA ligase (SUCLG2). Through repressing this TCA enzyme and its metabolite, succinate, metformin activated prolyl hydroxylases (PHDs), resulting in the degradation of hypoxia-inducible factor 1α (HIF1α) in mesothelial cells. Disruption of HIF1α-driven IL-8 signaling in mesothelial cells by metformin results in reduced OvCa invasion in an organotypic 3D model. These findings indicate that tumor-promoting signaling between mesothelial and OvCa cells in the TME can be targeted using metformin.

Organotypic 3D Models of the Ovarian Cancer Tumor Microenvironment.

Ovarian cancer progression involves multifaceted and variable tumor microenvironments (TMEs), from the in situ carcinoma in the fallopian tube or ovary to dissemination into the peritoneal cavity as single cells or spheroids and attachment to the mesothelial-lined surfaces of the omentum, bowel, and abdominal wall. The TME comprises the tumor vasculature and lymphatics (including endothelial cells and pericytes), in addition to mesothelial cells, fibroblasts, immune cells, adipocytes and extracellular matrix (ECM) proteins. When generating 3D models of the ovarian cancer TME, researchers must incorporate the most relevant stromal components depending on the TME in question (e.g., early or late disease). Such complexity cannot be captured by monolayer 2D culture systems. Moreover, immortalized stromal cell lines, such as mesothelial or fibroblast cell lines, do not always behave the same as primary cells whose response in functional assays may vary from donor to donor; 3D models with primary stromal cells may have more physiological relevance than those using stromal cell lines. In the current review, we discuss the latest developments in organotypic 3D models of the ovarian cancer early metastatic microenvironment. Organotypic culture models comprise two or more interacting cell types from a particular tissue. We focus on organotypic 3D models that include at least one type of primary stromal cell type in an ECM background, such as collagen or fibronectin, plus ovarian cancer cells. We provide an overview of the two most comprehensive current models-a 3D model of the omental mesothelium and a microfluidic model. We describe the cellular and non-cellular components of the models, the incorporation of mechanical forces, and how the models have been adapted and utilized in functional assays. Finally, we review a number of 3D models that do not incorporate primary stromal cells and summarize how integration of current models may be the next essential step in tackling the complexity of the different ovarian cancer TMEs.

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.

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.

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.

Cytokine and growth factor expression by HTLV-1 Lck-tax transgenic cells in SCID mice.

The Tax protein encoded by HTLV-1 plays a key role in the development of ATL in infected individuals. Our previous studies showed that tax transgenic mice develop disease that is almost identical to human ATL, with widespread organ invasion by lymphomatous cells and the development of leukemia. The same pathology develops rapidly in SCID mice engrafted with cells from the transgenic animals. In the present study, we used this SCID model to analyze the expression levels of several cytokines, growth factors, and adhesion molecules to determine their possible involvement in the development of disease. We showed that Tax expression was undetectable at the protein level in the tax-transformed cells used to inoculate the SCID mice and that these cells displayed constitutive NF-kappaB and Akt activity. We demonstrated significant differences in the levels of circulating PDGF-BB, TNF-alpha, sICAM-1, and sVCAM-1 in inoculated animals. Cell-surface staining of the tax transgenic cells showed that they do not express receptors for any of the upregulated growth factors. Significant differences were not found in the secreted levels of bFGF, MMP9, VEGF, or E-selectin, whereas IL-2, IL-15, IL-6, IL-1beta, and IFN-gamma expression was undetectable. Even though the number of factors analyzed is limited, our study identified TNF-alpha, PDGF-BB, and the adhesion molecules sICAM-1 and sVCAM-1 as factors that may contribute to the high levels of organ infiltration by leukemic cells in this tax transgenic SCID model.

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.