Transcriptional profiling of early differentiation of primary human mesenchymal stem cells into chondrocytes.

Articular cartilage has only very limited regenerative capacities in humans. Tissue engineering techniques for cartilage damage repair are limited in the production of hyaline cartilage. Mesenchymal stem/stromal cells (MSCs) are multipotent stem cells and can be differentiated into mature cartilage cells, chondrocytes, which could be used for repairing damaged cartilage. Chondrogenesis is a highly complex, relatively inefficient process lasting over 3 weeks in vitro. Methods: In order to better understand chondrogenic differentiation, especially the commitment phase, we have performed transcriptional profiling of MSC differentiation into chondrocytes from early timepoints starting 15 minutes after induction to 16 hours and fully differentiated chondrocytes at 21 days in triplicates.

Riboregulation of Enolase 1 activity controls glycolysis and embryonic stem cell differentiation.

Differentiating stem cells must coordinate their metabolism and fate trajectories. Here, we report that the catalytic activity of the glycolytic enzyme Enolase 1 (ENO1) is directly regulated by RNAs leading to metabolic rewiring in mouse embryonic stem cells (mESCs). We identify RNA ligands that specifically inhibit ENO1's enzymatic activity in vitro and diminish glycolysis in cultured human cells and mESCs. Pharmacological inhibition or RNAi-mediated depletion of the protein deacetylase SIRT2 increases ENO1's acetylation and enhances its RNA binding. Similarly, induction of mESC differentiation leads to increased ENO1 acetylation, enhanced RNA binding, and inhibition of glycolysis. Stem cells expressing mutant forms of ENO1 that escape or hyper-activate this regulation display impaired germ layer differentiation. Our findings uncover acetylation-driven riboregulation of ENO1 as a physiological mechanism of glycolytic control and of the regulation of stem cell differentiation. Riboregulation may represent a more widespread principle of biological control.

RNA-binding proteins in human genetic disease.

RNA-binding proteins (RBPs) are critical effectors of gene expression, and as such their malfunction underlies the origin of many diseases. RBPs can recognize hundreds of transcripts and form extensive regulatory networks that help to maintain cell homeostasis. System-wide unbiased identification of RBPs has increased the number of recognized RBPs into the four-digit range and revealed new paradigms: from the prevalence of structurally disordered RNA-binding regions with roles in the formation of membraneless organelles to unsuspected and potentially pervasive connections between intermediary metabolism and RNA regulation. Together with an increasingly detailed understanding of molecular mechanisms of RBP function, these insights are facilitating the development of new therapies to treat malignancies. Here, we provide an overview of RBPs involved in human genetic disorders, both Mendelian and somatic, and discuss emerging aspects in the field with emphasis on molecular mechanisms of disease and therapeutic interventions.

Plasticity of nuclear and cytoplasmic stress responses of RNA-binding proteins.

Cellular stress causes multifaceted reactions to trigger adaptive responses to environmental cues at all levels of the gene expression pathway. RNA-binding proteins (RBP) are key contributors to stress-induced regulation of RNA fate and function. Here, we uncover the plasticity of the RNA interactome in stressed cells, differentiating between responses in the nucleus and in the cytoplasm. We applied enhanced RNA interactome capture (eRIC) analysis preceded by nucleo-cytoplasmic fractionation following arsenite-induced oxidative stress. The data reveal unexpectedly compartmentalized RNA interactomes and their responses to stress, including differential responses of RBPs in the nucleus versus the cytoplasm, which would have been missed by whole cell analyses.

An Integrated Global Analysis of Compartmentalized HRAS Signaling.

Modern omics technologies allow us to obtain global information on different types of biological networks. However, integrating these different types of analyses into a coherent framework for a comprehensive biological interpretation remains challenging. Here, we present a conceptual framework that integrates protein interaction, phosphoproteomics, and transcriptomics data. Applying this method to analyze HRAS signaling from different subcellular compartments shows that spatially defined networks contribute specific functions to HRAS signaling. Changes in HRAS protein interactions at different sites lead to different kinase activation patterns that differentially regulate gene transcription. HRAS-mediated signaling is the strongest from the cell membrane, but it regulates the largest number of genes from the endoplasmic reticulum. The integrated networks provide a topologically and functionally resolved view of HRAS signaling. They reveal distinct HRAS functions including the control of cell migration from the endoplasmic reticulum and TP53-dependent cell survival when signaling from the Golgi apparatus.

Discovery of RNA-binding proteins and characterization of their dynamic responses by enhanced RNA interactome capture.

Following the realization that eukaryotic RNA-binding proteomes are substantially larger than anticipated, we must now understand their detailed composition and dynamics. Methods such as RNA interactome capture (RIC) have begun to address this need. However, limitations of RIC have been reported. Here we describe enhanced RNA interactome capture (eRIC), a method based on the use of an LNA-modified capture probe, which yields numerous advantages including greater specificity and increased signal-to-noise ratios compared to existing methods. In Jurkat cells, eRIC reduces the rRNA and DNA contamination by >10-fold compared to RIC and increases the detection of RNA-binding proteins. Due to its low background, eRIC also empowers comparative analyses of changes of RNA-bound proteomes missed by RIC. For example, in cells treated with dimethyloxalylglycine, which inhibits RNA demethylases, eRIC identifies m6A-responsive RNA-binding proteins that escape RIC. eRIC will facilitate the unbiased characterization of RBP dynamics in response to biological and pharmacological cues.

Retinoic acid and TGF-ß signalling cooperate to overcome MYCN-induced retinoid resistance.

BACKGROUND: Retinoid therapy is widely employed in clinical oncology to differentiate malignant cells into their more benign counterparts. However, certain high-risk cohorts, such as patients with MYCN-amplified neuroblastoma, are innately resistant to retinoid therapy. Therefore, we employed a precision medicine approach to globally profile the retinoid signalling response and to determine how an excess of cellular MYCN antagonises these signalling events to prevent differentiation and confer resistance. METHODS: We applied RNA sequencing (RNA-seq) and interaction proteomics coupled with network-based systems level analysis to identify targetable vulnerabilities of MYCN-mediated retinoid resistance. We altered MYCN expression levels in a MYCN-inducible neuroblastoma cell line to facilitate or block retinoic acid (RA)-mediated neuronal differentiation. The relevance of differentially expressed genes and transcriptional regulators for neuroblastoma outcome were then confirmed using existing patient microarray datasets. RESULTS: We determined the signalling networks through which RA mediates neuroblastoma differentiation and the inhibitory perturbations to these networks upon MYCN overexpression. We revealed opposing regulation of RA and MYCN on a number of differentiation-relevant genes, including LMO4, CYP26A1, ASCL1, RET, FZD7 and DKK1. Furthermore, we revealed a broad network of transcriptional regulators involved in regulating retinoid responsiveness, such as Neurotrophin, PI3K, Wnt and MAPK, and epigenetic signalling. Of these regulators, we functionally confirmed that MYCN-driven inhibition of transforming growth factor beta (TGF-ß) signalling is a vulnerable node of the MYCN network and that multiple levels of cross-talk exist between MYCN and TGF-ß. Co-targeting of the retinoic acid and TGF-ß pathways, through RA and kartogenin (KGN; a TGF-ß signalling activating small molecule) combination treatment, induced the loss of viability of MYCN-amplified retinoid-resistant neuroblastoma cells. CONCLUSIONS: Our approach provides a powerful precision oncology tool for identifying the driving signalling networks for malignancies not primarily driven by somatic mutations, such as paediatric cancers. By applying global omics approaches to the signalling networks regulating neuroblastoma differentiation and stemness, we have determined the pathways involved in the MYCN-mediated retinoid resistance, with TGF-ß signalling being a key regulator. These findings revealed a number of combination treatments likely to improve clinical response to retinoid therapy, including co-treatment with retinoids and KGN, which may prove valuable in the treatment of high-risk MYCN-amplified neuroblastoma.

Next Generation RNA Sequencing Analysis Reveals Expression of a Transient EMT Profile During Early Organization of MCF10A Cells in 3D.

RNA sequencing is a technique widely used to identify and characterize gene expression patterns. We demonstrate that this method can be applied to screen expression profiles in mammary epithelial cells cultured in 3D, supported by a natural laminin-rich extracellular matrix, but requires several specific steps in the preparation of the RNA samples. Here we describe the use of RNA sequencing to analyze mRNA patterns in MCF10A human mammary epithelial cells cultured under 3D conditions in a laminin-rich extracellular matrix. We focus on our methods for total RNA extraction at early time points during the formation and maturation of 3D acinus structures in these cultures and provide examples of our results and downstream analysis.

Wnt signalling is a bi-directional vulnerability of cancer cells.

Wnt signalling is involved in the formation, metastasis and relapse of a wide array of cancers. However, there is ongoing debate as to whether activation or inhibition of the pathway holds the most promise as a therapeutic treatment for cancer, with conflicting evidence from a variety of tumour types. We show that Wnt/ß-catenin signalling is a bi-directional vulnerability of neuroblastoma, malignant melanoma and colorectal cancer, with hyper-activation or repression of the pathway both representing a promising therapeutic strategy, even within the same cancer type. Hyper-activation directs cancer cells to undergo apoptosis, even in cells oncogenically driven by ß-catenin. Wnt inhibition blocks proliferation of cancer cells and promotes neuroblastoma differentiation. Wnt and retinoic acid co-treatments synergise, representing a promising combination treatment for MYCN-amplified neuroblastoma. Additionally, we report novel cross-talks between MYCN and ß-catenin signalling, which repress normal ß-catenin mediated transcriptional regulation. A ß-catenin target gene signature could predict patient outcome, as could the expression level of its DNA binding partners, the TCF/LEFs. This ß-catenin signature provides a tool to identify neuroblastoma patients likely to benefit from Wnt-directed therapy. Taken together, we show that Wnt/ß-catenin signalling is a bi-directional vulnerability of a number of cancer entities, and potentially a more broadly conserved feature of malignant cells.

Integrative omics reveals MYCN as a global suppressor of cellular signalling and enables network-based therapeutic target discovery in neuroblastoma.

Despite intensive study, many mysteries remain about the MYCN oncogene's functions. Here we focus on MYCN's role in neuroblastoma, the most common extracranial childhood cancer. MYCN gene amplification occurs in 20% of cases, but other recurrent somatic mutations are rare. This scarcity of tractable targets has hampered efforts to develop new therapeutic options. We employed a multi-level omics approach to examine MYCN functioning and identify novel therapeutic targets for this largely un-druggable oncogene. We used systems medicine based computational network reconstruction and analysis to integrate a range of omic techniques: sequencing-based transcriptomics, genome-wide chromatin immunoprecipitation, siRNA screening and interaction proteomics, revealing that MYCN controls highly connected networks, with MYCN primarily supressing the activity of network components. MYCN's oncogenic functions are likely independent of its classical heterodimerisation partner, MAX. In particular, MYCN controls its own protein interaction network by transcriptionally regulating its binding partners.Our network-based approach identified vulnerable therapeutically targetable nodes that function as critical regulators or effectors of MYCN in neuroblastoma. These were validated by siRNA knockdown screens, functional studies and patient data. We identified ß-estradiol and MAPK/ERK as having functional cross-talk with MYCN and being novel targetable vulnerabilities of MYCN-amplified neuroblastoma. These results reveal surprising differences between the functioning of endogenous, overexpressed and amplified MYCN, and rationalise how different MYCN dosages can orchestrate cell fate decisions and cancerous outcomes. Importantly, this work describes a systems-level approach to systematically uncovering network based vulnerabilities and therapeutic targets for multifactorial diseases by integrating disparate omic data types.

The RNA-binding proteomes from yeast to man harbour conserved enigmRBPs.

RNA-binding proteins (RBPs) exert a broad range of biological functions. To explore the scope of RBPs across eukaryotic evolution, we determined the in vivo RBP repertoire of the yeast Saccharomyces cerevisiae and identified 678 RBPs from yeast and additionally 729 RBPs from human hepatocytic HuH-7 cells. Combined analyses of these and recently published data sets define the core RBP repertoire conserved from yeast to man. Conserved RBPs harbour defined repetitive motifs within disordered regions, which display striking evolutionary expansion. Only 60% of yeast and 73% of the human RBPs have functions assigned to RNA biology or structural motifs known to convey RNA binding, and many intensively studied proteins surprisingly emerge as RBPs (termed 'enigmRBPs'), including almost all glycolytic enzymes, pointing to emerging connections between gene regulation and metabolism. Analyses of the mitochondrial hydroxysteroid dehydrogenase (HSD17B10) uncover the RNA-binding specificity of an enigmRBP.

Measuring Transcription Rate Changes via Time-Course 4-Thiouridine Pulse-Labelling Improves Transcriptional Target Identification.

Identifying changes in the transcriptional regulation of target genes from high-throughput studies is important for unravelling molecular mechanisms controlled by a given perturbation. When measuring global transcript levels only, the effect of the perturbation [e.g., transcription factor (TF) overexpression or drug treatment] on its target genes is often obscured by delayed feedback and secondary effects until the changes are fully propagated. As a proof of principle, we show that selective measuring of transcripts that are only synthesised after a perturbation [4-thiouridine (4sU) sequencing (4sU-seq)] is a more sensitive method to identify targets and time-dependent transcriptional responses than global transcript profiling. By metabolically labelling RNA in a time-course setup, we could vastly increase the sensitivity of MYCN target gene detection compared to traditional RNA sequencing. The validity of targets identified by 4sU-seq was demonstrated using chromatin immunoprecipitation sequencing and neuroblastoma microarray tumour data. Here, we describe the methodology, both molecular biology and computational aspects, required to successfully apply this 4sU-seq approach.

p19-INK4d inhibits neuroblastoma cell growth, induces differentiation and is hypermethylated and downregulated in MYCN-amplified neuroblastomas.

Uncontrolled cell cycle entry, resulting from deregulated CDK-RB1-E2F pathway activity, is a crucial determinant of neuroblastoma cell malignancy. Here we identify neuroblastoma-suppressive functions of the p19-INK4d CDK inhibitor and uncover mechanisms of its repression in high-risk neuroblastomas. Reduced p19-INK4d expression was associated with poor event-free and overall survival and neuroblastoma risk factors including amplified MYCN in a set of 478 primary neuroblastomas. High MYCN expression repressed p19-INK4d mRNA and protein levels in different neuroblastoma cell models with conditional MYCN expression. MassARRAY and 450K methylation analyses of 105 primary neuroblastomas uncovered a differentially methylated region within p19-INK4d. Hypermethylation of this region was associated with reduced p19-INK4d expression. In accordance, p19-INK4d expression was activated upon treatment with the demethylating agent, 2'-deoxy-5-azacytidine, in neuroblastoma cell lines. Ectopic p19-INK4d expression decreased viability, clonogenicity and the capacity for anchorage-independent growth of neuroblastoma cells, and shifted the cell cycle towards the G1/0 phase. p19-INK4d also induced neurite-like processes and markers of neuronal differentiation. Moreover, neuroblastoma cell differentiation, induced by all-trans retinoic acid or NGF-NTRK1-signaling, activated p19-INK4d expression. Our findings pinpoint p19-INK4d as a neuroblastoma suppressor and provide evidence for MYCN-mediated repression and for epigenetic silencing of p19-INK4d by DNA hypermethylation in high-risk neuroblastomas.

GSK3 inhibitors regulate MYCN mRNA levels and reduce neuroblastoma cell viability through multiple mechanisms, including p53 and Wnt signaling.

Neuroblastoma is an embryonal tumor accounting for approximately 15% of childhood cancer deaths. There exists a clinical need to identify novel therapeutic targets, particularly for treatment-resistant forms of neuroblastoma. Therefore, we investigated the role of the neuronal master regulator GSK3 in controlling neuroblastoma cell fate. We identified novel GSK3-mediated regulation of MYC (c-MYC and MYCN) mRNA levels, which may have implications for numerous MYC-driven cancers. In addition, we showed that certain GSK3 inhibitors induced large-scale cell death in neuroblastoma cells, primarily through activating apoptosis. mRNA-seq of GSK3 inhibitor-treated cells was performed and subsequent pathway analysis revealed that multiple signaling pathways contributed to the loss of neuroblastoma cell viability. The contribution of two of the signaling pathways highlighted by the mRNA-seq analysis was functionally validated. Inhibition of the p53 tumor suppressor partly rescued the cell death phenotype, whereas activation of canonical Wnt signaling contributed to the loss of viability, in a p53-independent manner. Two GSK3 inhibitors (BIO-acetoxime and LiCl) and one small-molecule Wnt agonist (Wnt Agonist 1) demonstrated therapeutic potential for neuroblastoma treatment. These inhibitors reduced the viability of numerous neuroblastoma cell lines, even those derived from high-risk MYCN-amplified metastatic tumors, for which effective therapeutics are currently lacking. Furthermore, although LiCl was lethal to neuroblastoma cells, it did not reduce the viability of differentiated neurons. Taken together our data suggest that these small molecules may hold potential as effective therapeutic agents for the treatment of neuroblastoma and other MYC-driven cancers.

c-Myc regulates RNA splicing of the A-Raf kinase and its activation of the ERK pathway.

A-Raf kinase can inhibit apoptosis by binding to the proapoptotic mammalian sterile 20-like kinase (MST2). This function relies on expression of hnRNP H, which ensures the correct splicing of a-raf mRNA needed to produce full-length A-Raf protein. Here, we showed that expression of hnRNP H and production of full-length A-Raf is positively controlled by c-Myc. Low c-Myc reduces hnRNP H expression and switches a-raf splicing to produce A-Raf(short), a truncated protein. Importantly, A-Raf(short) fails to regulate MST2 but retains the Ras-binding domain such that it functions as a dominant negative mutant suppressing Ras activation and transformation. Human colon and head and neck cancers exhibit high hnRNP H and high c-Myc levels resulting in enhanced A-Raf expression and reduced expression of A-Raf(short). Conversely, in normal cells and tissues in which c-Myc and hnRNP H are low, A-Raf(short) suppresses extracellular signal regulated kinase activation such that it may act as a safeguard against oncogenic transformation. Our findings offered a new paradigm to understand how c-Myc coordinates diverse cell functions by directly affecting alternate splicing of key signaling components.

GeneSigDB--a curated database of gene expression signatures.

The primary objective of most gene expression studies is the identification of one or more gene signatures; lists of genes whose transcriptional levels are uniquely associated with a specific biological phenotype. Whilst thousands of experimentally derived gene signatures are published, their potential value to the community is limited by their computational inaccessibility. Gene signatures are embedded in published article figures, tables or in supplementary materials, and are frequently presented using non-standard gene or probeset nomenclature. We present GeneSigDB (http://compbio.dfci.harvard.edu/genesigdb) a manually curated database of gene expression signatures. GeneSigDB release 1.0 focuses on cancer and stem cells gene signatures and was constructed from more than 850 publications from which we manually transcribed 575 gene signatures. Most gene signatures (n = 560) were successfully mapped to the genome to extract standardized lists of EnsEMBL gene identifiers. GeneSigDB provides the original gene signature, the standardized gene list and a fully traceable gene mapping history for each gene from the original transcribed data table through to the standardized list of genes. The GeneSigDB web portal is easy to search, allows users to compare their own gene list to those in the database, and download gene signatures in most common gene identifier formats.