A large-scale binding and functional map of human RNA-binding proteins.

Many proteins regulate the expression of genes by binding to specific regions encoded in the genome1. Here we introduce a new data set of RNA elements in the human genome that are recognized by RNA-binding proteins (RBPs), generated as part of the Encyclopedia of DNA Elements (ENCODE) project phase III. This class of regulatory elements functions only when transcribed into RNA, as they serve as the binding sites for RBPs that control post-transcriptional processes such as splicing, cleavage and polyadenylation, and the editing, localization, stability and translation of mRNAs. We describe the mapping and characterization of RNA elements recognized by a large collection of human RBPs in K562 and HepG2 cells. Integrative analyses using five assays identify RBP binding sites on RNA and chromatin in vivo, the in vitro binding preferences of RBPs, the function of RBP binding sites and the subcellular localization of RBPs, producing 1,223 replicated data sets for 356 RBPs. We describe the spectrum of RBP binding throughout the transcriptome and the connections between these interactions and various aspects of RNA biology, including RNA stability, splicing regulation and RNA localization. These data expand the catalogue of functional elements encoded in the human genome by the addition of a large set of elements that function at the RNA level by interacting with RBPs.

Sequence, Structure, and Context Preferences of Human RNA Binding Proteins.

RNA binding proteins (RBPs) orchestrate the production, processing, and function of mRNAs. Here, we present the affinity landscapes of 78 human RBPs using an unbiased assay that determines the sequence, structure, and context preferences of these proteins in vitro by deep sequencing of bound RNAs. These data enable construction of "RNA maps" of RBP activity without requiring crosslinking-based assays. We found an unexpectedly low diversity of RNA motifs, implying frequent convergence of binding specificity toward a relatively small set of RNA motifs, many with low compositional complexity. Offsetting this trend, however, we observed extensive preferences for contextual features distinct from short linear RNA motifs, including spaced "bipartite" motifs, biased flanking nucleotide composition, and bias away from or toward RNA structure. Our results emphasize the importance of contextual features in RNA recognition, which likely enable targeting of distinct subsets of transcripts by different RBPs that recognize the same linear motif.

RNA Sequence Context Effects Measured In Vitro Predict In Vivo Protein Binding and Regulation.

Many RNA binding proteins (RBPs) bind specific RNA sequence motifs, but only a small fraction (∼15%-40%) of RBP motif occurrences are occupied in vivo. To determine which contextual features discriminate between bound and unbound motifs, we performed an in vitro binding assay using 12,000 mouse RNA sequences with the RBPs MBNL1 and RBFOX2. Surprisingly, the strength of binding to motif occurrences in vitro was significantly correlated with in vivo binding, developmental regulation, and evolutionary age of alternative splicing. Multiple lines of evidence indicate that the primary context effect that affects binding in vitro and in vivo is RNA secondary structure. Large-scale combinatorial mutagenesis of unfavorable sequence contexts revealed a consistent pattern whereby mutations that increased motif accessibility improved protein binding and regulatory activity. Our results indicate widespread inhibition of motif binding by local RNA secondary structure and suggest that mutations that alter sequence context commonly affect RBP binding and regulation.

Antagonistic regulation of mRNA expression and splicing by CELF and MBNL proteins.

RNA binding proteins of the conserved CUGBP1, Elav-like factor (CELF) family contribute to heart and skeletal muscle development and are implicated in myotonic dystrophy (DM). To understand their genome-wide functions, we analyzed the transcriptome dynamics following induction of CELF1 or CELF2 in adult mouse heart and of CELF1 in muscle by RNA-seq, complemented by crosslinking/immunoprecipitation-sequencing (CLIP-seq) analysis of mouse cells and tissues to distinguish direct from indirect regulatory targets. We identified hundreds of mRNAs bound in their 3' UTRs by both CELF1 and the developmentally induced MBNL1 protein, a threefold greater overlap in target messages than expected, including messages involved in development and cell differentiation. The extent of 3' UTR binding by CELF1 and MBNL1 predicted the degree of mRNA repression or stabilization, respectively, following CELF1 induction. However, CELF1's RNA binding specificity in vitro was not detectably altered by coincubation with recombinant MBNL1. These findings support a model in which CELF and MBNL proteins bind independently to mRNAs but functionally compete to specify down-regulation or localization/stabilization, respectively, of hundreds of mRNA targets. Expression of many alternative 3' UTR isoforms was altered following CELF1 induction, with 3' UTR binding associated with down-regulation of isoforms and genes. The splicing of hundreds of alternative exons was oppositely regulated by these proteins, confirming an additional layer of regulatory antagonism previously observed in a handful of cases. The regulatory relationships between CELFs and MBNLs in control of both mRNA abundance and splicing appear to have evolved to enhance developmental transitions in major classes of heart and muscle genes.

Alternative splicing of RNA triplets is often regulated and accelerates proteome evolution.

Thousands of human genes contain introns ending in NAGNAG (N any nucleotide), where both NAGs can function as 3' splice sites, yielding isoforms that differ by inclusion/exclusion of three bases. However, few models exist for how such splicing might be regulated, and some studies have concluded that NAGNAG splicing is purely stochastic and nonfunctional. Here, we used deep RNA-Seq data from 16 human and eight mouse tissues to analyze the regulation and evolution of NAGNAG splicing. Using both biological and technical replicates to estimate false discovery rates, we estimate that at least 25% of alternatively spliced NAGNAGs undergo tissue-specific regulation in mammals, and alternative splicing of strongly tissue-specific NAGNAGs was 10 times as likely to be conserved between species as was splicing of non-tissue-specific events, implying selective maintenance. Preferential use of the distal NAG was associated with distinct sequence features, including a more distal location of the branch point and presence of a pyrimidine immediately before the first NAG, and alteration of these features in a splicing reporter shifted splicing away from the distal site. Strikingly, alignments of orthologous exons revealed a ∼15-fold increase in the frequency of three base pair gaps at 3' splice sites relative to nearby exon positions in both mammals and in Drosophila. Alternative splicing of NAGNAGs in human was associated with dramatically increased frequency of exon length changes at orthologous exon boundaries in rodents, and a model involving point mutations that create, destroy, or alter NAGNAGs can explain both the increased frequency and biased codon composition of gained/lost sequence observed at the beginnings of exons. This study shows that NAGNAG alternative splicing generates widespread differences between the proteomes of mammalian tissues, and suggests that the evolutionary trajectories of mammalian proteins are strongly biased by the locations and phases of the introns that interrupt coding sequences.

The conformation of microRNA seed regions in native microRNPs is prearranged for presentation to mRNA targets.

MicroRNAs control gene expression by post-transcriptional down-regulation of their target mRNAs. Complementarity between the seed region (nucleotides 2-8) of a microRNA and the 3'-UTR of its target mRNA is the key determinant in recognition. However, the structural basis of the ability of the seed region to dominate target recognition in eukaryotic argonaute complexes has not been directly demonstrated. To better understand this problem, we performed chemical probing of microRNAs held in native argonaute-containing complexes isolated from Caenorhabditis elegans. Direct probing of the RNA backbone in isolated native microRNP complexes shows that the conformation of the seed region is uniquely constrained, while the rest of the microRNA structure is conformationally flexible. Probing the Watson-Crick edges of the bases shows that bases 2-4 are largely inaccessible to solvent, while seed region bases 5-8 are readily modified; collectively our probing results suggest a model in which these bases are primed for initiating base pairing with the target mRNA. In addition, an unusual DMS reactivity with U at position 6 is observed. We propose that interaction of miRNAs with argonaute proteins pre-organizes the structure of the seed sequence for specific recognition of target mRNAs.

Global analysis of alternative splicing uncovers developmental regulation of nonsense-mediated decay in C. elegans.

Alternative splicing coupled to nonsense-mediated decay (AS-NMD) is a mechanism for post-transcriptional regulation of gene expression. We analyzed the global effects of mutations in seven genes of the C. elegans NMD pathway on AS isoform ratios. We find that mutations in two NMD factors, smg-6 and smg-7, have weaker global effects relative to mutations in other smg genes. We did an in-depth analysis of 12 pre-mRNA splicing factor genes that are subject to AS-NMD. For four of these, changes in the ratio of alternatively spliced isoforms during development are caused by developmentally regulated inhibition of NMD, and not by changes in alternative splicing. Using sucrose gradient analysis of mRNAs undergoing translation, we find several examples of NMD-dependent enrichment of premature termination codon (PTC) isoforms in the monosome fraction. In contrast, we present evidence of two genes for which the PTC-containing isoforms are found in polysomes and have a translational profile similar to non-PTC-containing transcripts from the same gene. We propose that NMD of certain alternatively spliced isoforms is regulated, and that some stabilized NMD targets may be translated.

The sarcin-ricin loop of 23S rRNA is essential for assembly of the functional core of the 50S ribosomal subunit.

The sarcin-ricin loop (SRL) of 23S rRNA in the large ribosomal subunit is a factor-binding site that is essential for GTP-catalyzed steps in translation, but its precise functional role is thus far unknown. Here, we replaced the 15-nucleotide SRL with a GAAA tetraloop and affinity purified the mutant 50S subunits for functional and structural analysis in vitro. The SRL deletion caused defects in elongation-factor-dependent steps of translation and, unexpectedly, loss of EF-Tu-independent A-site tRNA binding. Detailed chemical probing analysis showed disruption of a network of rRNA tertiary interactions that hold together the 23S rRNA elements of the functional core of the 50S subunit, accompanied by loss of ribosomal protein L16. Our results reveal an influence of the SRL on the higher-order structure of the 50S subunit, with implications for its role in translation.

Early Assessment of Molecular Progression and Response by Whole-genome Circulating Tumor DNA in Advanced Solid Tumors.

Treatment response assessment for patients with advanced solid tumors is complex and existing methods require greater precision. Current guidelines rely on imaging, which has known limitations, including the time required to show a deterministic change in target lesions. Serial changes in whole-genome (WG) circulating tumor DNA (ctDNA) were used to assess response or resistance to treatment early in the treatment course. Ninety-six patients with advanced cancer were prospectively enrolled (91 analyzed and 5 excluded), and blood was collected before and after initiation of a new, systemic treatment. Plasma cell-free DNA libraries were prepared for either WG or WG bisulfite sequencing. Longitudinal changes in the fraction of ctDNA were quantified to retrospectively identify molecular progression (MP) or major molecular response (MMR). Study endpoints were concordance with first follow-up imaging (FFUI) and stratification of progression-free survival (PFS) and overall survival (OS). Patients with MP (n = 13) had significantly shorter PFS (median 62 days vs. 310 days) and OS (255 days vs. not reached). Sensitivity for MP to identify clinical progression was 54% and specificity was 100%. MP calls were from samples taken a median of 28 days into treatment and 39 days before FFUI. Patients with MMR (n = 27) had significantly longer PFS and OS compared with those with neither call (n = 51). These results demonstrated that ctDNA changes early after treatment initiation inform response to treatment and correlate with long-term clinical outcomes. Once validated, molecular response assessment can enable early treatment change minimizing side effects and costs associated with additional cycles of ineffective treatment.

Enhanced CLIP Uncovers IMP Protein-RNA Targets in Human Pluripotent Stem Cells Important for Cell Adhesion and Survival.

Human pluripotent stem cells (hPSCs) require precise control of post-transcriptional RNA networks to maintain proliferation and survival. Using enhanced UV crosslinking and immunoprecipitation (eCLIP), we identify RNA targets of the IMP/IGF2BP family of RNA-binding proteins in hPSCs. At the broad region and binding site levels, IMP1 and IMP2 show reproducible binding to a large and overlapping set of 3' UTR-enriched targets. RNA Bind-N-seq applied to recombinant full-length IMP1 and IMP2 reveals CA-rich motifs that are enriched in eCLIP-defined binding sites. We observe that IMP1 loss in hPSCs recapitulates IMP1 phenotypes, including a reduction in cell adhesion and increase in cell death. For cell adhesion, we find IMP1 maintains levels of integrin mRNA specifically regulating RNA stability of ITGB5 in hPSCs. Additionally, we show that IMP1 can be linked to hPSC survival via direct target BCL2. Thus, transcriptome-wide binding profiles identify hPSC targets modulating well-characterized IMP1 roles.

Distinct and shared functions of ALS-associated proteins TDP-43, FUS and TAF15 revealed by multisystem analyses.

The RNA-binding protein (RBP) TAF15 is implicated in amyotrophic lateral sclerosis (ALS). To compare TAF15 function to that of two ALS-associated RBPs, FUS and TDP-43, we integrate CLIP-seq and RNA Bind-N-Seq technologies, and show that TAF15 binds to ∼4,900 RNAs enriched for GGUA motifs in adult mouse brains. TAF15 and FUS exhibit similar binding patterns in introns, are enriched in 3' untranslated regions and alter genes distinct from TDP-43. However, unlike FUS and TDP-43, TAF15 has a minimal role in alternative splicing. In human neural progenitors, TAF15 and FUS affect turnover of their RNA targets. In human stem cell-derived motor neurons, the RNA profile associated with concomitant loss of both TAF15 and FUS resembles that observed in the presence of the ALS-associated mutation FUS R521G, but contrasts with late-stage sporadic ALS patients. Taken together, our findings reveal convergent and divergent roles for FUS, TAF15 and TDP-43 in RNA metabolism.

RNA Bind-n-Seq: Measuring the Binding Affinity Landscape of RNA-Binding Proteins.

RNA-binding proteins (RBPs) coordinate post-transcriptional control of gene expression, often through sequence-specific recognition of primary transcripts or mature messenger RNAs. Hundreds of RBPs are encoded in the human genome, most with undefined or incompletely defined biological roles. Understanding the function of these factors will require the identification of each RBP's distinct RNA binding specificity. RNA Bind-n-Seq (RBNS) is a high-throughput, cost-effective in vitro method capable of resolving sequence and secondary structure preferences of RBPs. Dissociation constants can also be inferred from RBNS data when provided with additional experimental information. Here, we describe the experimental procedures to perform RBNS and discuss important parameters of the method and ways that the experiment can be tailored to the specific RBP under study. Additionally, we present the conceptual framework and execution of the freely available RBNS computational pipeline and describe the outputs of the pipeline. Different approaches to quantify binding specificity, quality control metrics, and estimation of binding constants are also covered.

Musashi proteins are post-transcriptional regulators of the epithelial-luminal cell state.

The conserved Musashi (Msi) family of RNA binding proteins are expressed in stem/progenitor and cancer cells, but generally absent from differentiated cells, consistent with a role in cell state regulation. We found that Msi genes are rarely mutated but frequently overexpressed in human cancers and are associated with an epithelial-luminal cell state. Using ribosome profiling and RNA-seq analysis, we found that Msi proteins regulate translation of genes implicated in epithelial cell biology and epithelial-to-mesenchymal transition (EMT), and promote an epithelial splicing pattern. Overexpression of Msi proteins inhibited the translation of Jagged1, a factor required for EMT, and repressed EMT in cell culture and in mammary gland in vivo. Knockdown of Msis in epithelial cancer cells promoted loss of epithelial identity. Our results show that mammalian Msi proteins contribute to an epithelial gene expression program in neural and mammary cell types.