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.

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.

Enhanced crosslinking and immunoprecipitation (eCLIP) sequencing is a method for transcriptome-wide detection of binding sites of RNA-binding proteins (RBPs). However, identified crosslink sites can deviate from experimentally established functional elements of even well-studied RBPs. Current peak-calling strategies result in low replication and high false positive rates. Here, we present the R/Bioconductor package DEWSeq that makes use of replicate information and size-matched input controls. We benchmarked DEWSeq on 107 RBPs for which both eCLIP data and RNA sequence motifs are available and were able to more than double the number of motif-containing binding regions relative to standard eCLIP processing. The improvement not only relates to the number of binding sites (3.1-fold with known motifs for RBFOX2), but also their subcellular localization (1.9-fold of mitochondrial genes for FASTKD2) and structural targets (2.2-fold increase of stem-loop regions for SLBP. On several orthogonal CLIP-seq datasets, DEWSeq recovers a larger number of motif-containing binding sites (3.3-fold). DEWSeq is a well-documented R/Bioconductor package, scalable to adequate numbers of replicates, and tends to substantially increase the proportion and total number of RBP binding sites containing biologically relevant features.

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.

RNA-protein interactions are central to cardiac function, but how activity of individual RNA-binding protein is regulated through signaling cascades in cardiomyocytes during heart failure development is largely unknown. The mechanistic target of rapamycin kinase is a central signaling hub that controls mRNA translation in cardiomyocytes; however, a direct link between mTOR signaling and RNA-binding proteins in the heart has not been established. Integrative transcriptome and translatome analysis revealed mTOR dependent translational upregulation of the RNA binding protein Ybx1 during early pathological remodeling independent of mRNA levels. Ybx1 is necessary for pathological cardiomyocyte growth by regulating protein synthesis. To identify the molecular mechanisms how Ybx1 regulates cellular growth and protein synthesis, we identified mRNAs bound to Ybx1. We discovered that eucaryotic elongation factor 2 (Eef2) mRNA is bound to Ybx1, and its translation is upregulated during cardiac hypertrophy dependent on Ybx1 expression. Eef2 itself is sufficient to drive pathological growth by increasing global protein translation. Finally, Ybx1 depletion in vivo preserved heart function during pathological cardiac hypertrophy. Thus, activation of mTORC1 links pathological signaling cascades to altered gene expression regulation by activation of Ybx1 which in turn promotes translation through increased expression of Eef2.

Hepatocyte death is an important contributing factor in a number of diseases of the liver. PHD1 confers hypoxic sensitivity upon transcription factors including the hypoxia inducible factor (HIF) and nuclear factor-kappaB (NF-κB). Reduced PHD1 activity is linked to decreased apoptosis. Here, we investigated the underlying mechanism(s) in hepatocytes. Basal NF-κB activity was elevated in PHD1(-/-) hepatocytes compared to wild type controls. ChIP-seq analysis confirmed enhanced binding of NF-κB to chromatin in regions proximal to the promoters of genes involved in the regulation of apoptosis. Inhibition of NF-κB (but not knock-out of HIF-1 or HIF-2) reversed the anti-apoptotic effects of pharmacologic hydroxylase inhibition. We hypothesize that PHD1 inhibition leads to altered expression of NF-κB-dependent genes resulting in reduced apoptosis. This study provides new information relating to the possible mechanism of therapeutic action of hydroxylase inhibitors that has been reported in pre-clinical models of intestinal and hepatic disease.

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.

Small noncoding RNAs fulfill key functions in cellular and organismal biology, typically working in concert with RNA-binding proteins (RBPs). While proteome-wide methodologies have enormously expanded the repertoire of known RBPs, these methods do not distinguish RBPs binding to small noncoding RNAs from the rest. To specifically identify this relevant subclass of RBPs, we developed small noncoding RNA interactome capture (snRIC2C) based on the differential RNA-binding capacity of silica matrices (2C). We define the S. cerevisiae proteome of nearly 300 proteins that specifically binds to RNAs smaller than 200 nt in length (snRBPs), identifying informative distinctions from the total RNA-binding proteome determined in parallel. Strikingly, the snRBPs include most glycolytic enzymes from yeast. With further methodological developments using silica matrices, 12 tRNAs were identified as specific binders of the glycolytic enzyme GAPDH. We show that tRNA engagement of GAPDH is carbon source-dependent and regulated by the RNA polymerase III repressor Maf1, suggesting a regulatory interaction between glycolysis and RNA polymerase III activity. We conclude that snRIC2C and other 2C-derived methods greatly facilitate the study of RBPs, revealing previously unrecognized interactions.

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.

Proteins and RNA functionally and physically intersect in multiple biological processes, however, currently no universal method is available to purify protein-RNA complexes. Here, we introduce XRNAX, a method for the generic purification of protein-crosslinked RNA, and demonstrate its versatility to study the composition and dynamics of protein-RNA interactions by various transcriptomic and proteomic approaches. We show that XRNAX captures all RNA biotypes and use this to characterize the sub-proteomes that interact with coding and non-coding RNAs (ncRNAs) and to identify hundreds of protein-RNA interfaces. Exploiting the quantitative nature of XRNAX, we observe drastic remodeling of the RNA-bound proteome during arsenite-induced stress, distinct from autophagy-related changes in the total proteome. In addition, we combine XRNAX with crosslinking immunoprecipitation sequencing (CLIP-seq) to validate the interaction of ncRNA with lamin B1 and EXOSC2. Thus, XRNAX is a resourceful approach to study structural and compositional aspects of protein-RNA interactions to address fundamental questions in RNA-biology.

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.

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.

Vault RNAs (vtRNA) are small non-coding RNAs transcribed by RNA polymerase III found in many eukaryotes. Although they have been linked to drug resistance, apoptosis, and viral replication, their molecular functions remain unclear. Here, we show that vault RNAs directly bind the autophagy receptor sequestosome-1/p62 in human and murine cells. Overexpression of human vtRNA1-1 inhibits, while its antisense LNA-mediated knockdown enhances p62-dependent autophagy. Starvation of cells reduces the steady-state and p62-bound levels of vault RNA1-1 and induces autophagy. Mechanistically, p62 mutants that fail to bind vtRNAs display increased p62 homo-oligomerization and augmented interaction with autophagic effectors. Thus, vtRNA1-1 directly regulates selective autophagy by binding p62 and interference with oligomerization, a critical step of p62 function. Our data uncover a striking example of the potential of RNA to control protein functions directly, as previously recognized for protein-protein interactions and post-translational modifications.

Sufficient amino acid supplies are critical for protein synthesis and, thus, cell growth and proliferation. Specialized transporters mediate amino acid exchange across membranes and their regulation is critical for amino acid homeostasis. Here, we report that the DNA- and RNA-binding protein YBX3 regulates the expression of amino acid transporters. To investigate the functions of YBX3, we integrated proteomic and transcriptomic data from cells depleted of YBX3 with analyses of YBX3 RNA binding sites to identify RNAs directly regulated by YBX3. The data implicate YBX3 as a RNA-binding protein that regulates distinct sets of mRNAs by discrete mechanisms, including mRNA abundance. Among direct YBX3 targets, two solute carrier (SLC) amino acid transporters (SLC7A5 and SLC3A2) were identified. We show that YBX3 stabilizes these SLC mRNAs and that YBX3 depletion diminishes the expression of SLC7A5/SLC3A2, which specifically reduces SLC7A5/SLC3A2 amino acid substrates. Thus, YBX3 emerges as a key regulator of amino acid levels.

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.

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.

RNA-binding proteins (RBPs) control critical aspects of cardiomyocyte function, but the repertoire of active RBPs in cardiomyocytes during the growth response is largely unknown. We define RBPs in healthy and diseased cardiomyocytes at a system-wide level by RNA interactome capture. This identifies 67 cardiomyocyte-specific RBPs, including several contractile proteins. Furthermore, we identify the cytoplasmic polyadenylation element-binding protein 4 (Cpeb4) as a dynamic RBP, regulating cardiac growth both in vitro and in vivo. We identify mRNAs bound to and regulated by Cpeb4 in cardiomyocytes. Cpeb4 regulates cardiac remodeling by differential expression of transcription factors. Among Cpeb4 target mRNAs, two zinc finger transcription factors (Zeb1 and Zbtb20) are discovered. We show that Cpeb4 regulates the expression of these mRNAs and that Cpeb4 depletion increases their expression. Thus, Cpeb4 emerges as a critical regulator of cardiomyocyte function by differential binding to specific mRNAs in response to pathological growth stimulation.

System-wide approaches have unveiled an unexpected breadth of the RNA-bound proteomes of cultured cells. Corresponding information regarding RNA-binding proteins (RBPs) of mammalian organs is still missing, largely due to technical challenges. Here, we describe ex vivo enhanced RNA interactome capture (eRIC) to characterize the RNA-bound proteomes of three different mouse organs. The resulting organ atlases encompass more than 1300 RBPs active in brain, kidney or liver. Nearly a quarter (291) of these had formerly not been identified in cultured cells, with more than 100 being metabolic enzymes. Remarkably, RBP activity differs between organs independent of RBP abundance, suggesting organ-specific levels of control. Similarly, we identify systematic differences in RNA binding between animal organs and cultured cells. The pervasive RNA binding of enzymes of intermediary metabolism in organs points to tightly knit connections between gene expression and metabolism, and displays a particular enrichment for enzymes that use nucleotide cofactors. We describe a generically applicable refinement of the eRIC technology and provide an instructive resource of RBPs active in intact mammalian organs, including the brain.

Individual-nucleotide resolution crosslinking and immunoprecipitation (iCLIP) allows the determination of crosslinking sites of RNA-binding proteins (RBPs) on RNAs. iCLIP is based on ultraviolet light crosslinking of RBPs to RNA, reverse transcription and high-throughput sequencing of fragments terminating at the site of crosslinking. As a result, start sites of iCLIP fragments are expected to cluster with a narrow distribution, typically representing the site of direct interaction between the RBP and the RNA. Here we show that for several RBPs (eIF4A3, PTB, SRSF3, SRSF4 and hnRNP L), the start sites of iCLIP fragments show a fragment length-dependent broader distribution that can be shifted to positions upstream of the known RNA-binding site. We developed an analysis tool that identifies these shifts and can improve the positioning of RBP binding sites.

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.