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.

RNA Bind-n-Seq: quantitative assessment of the sequence and structural binding specificity of RNA binding proteins.

Specific protein-RNA interactions guide posttranscriptional gene regulation. Here, we describe RNA Bind-n-Seq (RBNS), a method that comprehensively characterizes sequence and structural specificity of RNA binding proteins (RBPs), and its application to the developmental alternative splicing factors RBFOX2, CELF1/CUGBP1, and MBNL1. For each factor, we recovered both canonical motifs and additional near-optimal binding motifs. RNA secondary structure inhibits binding of RBFOX2 and CELF1, while MBNL1 favors unpaired Us but tolerates C/G pairing in motifs containing UGC and/or GCU. Dissociation constants calculated from RBNS data using a novel algorithm correlated highly with values measured by surface plasmon resonance. Motifs identified by RBNS were conserved, were bound and active in vivo, and distinguished the subset of motifs enriched by CLIP-Seq that had regulatory activity. Together, our data demonstrate that RBNS complements crosslinking-based methods and show that in vivo binding and activity of these splicing factors is driven largely by intrinsic RNA affinity.

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.

Prevention of oral disease for long-term care and homebound elderly.

Despite the fact that many oral diseases afflicting the long-term care or homebound elderly are preventable or treatable, many older people do not seek available treatment, or their oral health care needs are not being met. The dental profession must, therefore, increase the preventative dental awareness of elders and make preventative and treatment services more accessible to this population. Interdisciplinary training and collaborative efforts among the dental profession, medical profession and caregivers are necessary in preventing oral disease for this geriatric population, which would improve not just oral health, but overall systemic health as well, thereby improving their quality of life.

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.

Evaluation and updating of the Osiris expert system for identification of Escherichia coli beta-lactam resistance phenotypes.

Osiris is a video zone size reader for disk diffusion tests that includes a built-in extended expert system (EES). We evaluated the efficacy of the Osiris EES for the identification of beta-lactam susceptibility phenotypes of Escherichia coli isolates. Fifteen beta-lactam agents and three beta-lactam-beta-lactamase inhibitor combinations were tested by the disk diffusion method against 50 E. coli strains with well-characterized resistance mechanisms. The strains were screened for the production of extended-spectrum beta-lactamase (ESBL) by the double-disk synergy test using a disk of amoxicillin-clavulanic acid with disks of the extended-spectrum cephalosporins and aztreonam. Overall, the EES accurately identified the phenotype for 78% of the strains, indicated an inexact phenotype for 17% of the strains, and could not find a matching phenotype for the remaining 5% of the strains. The percentage of correct identification for each resistance mechanism was 100% for inhibitor-resistant TEM and for TEM plus cephalosporinase, 88.9% for TEM and for ESBL, 70.8% for cephalosporinase overproduction, and 25% for oxacillinase. The main cause of discrepancy was the misidentification of oxacillinase as inhibitor-resistant TEM. The conventional double-disk synergy test failed to detect ESBL production in two strains (one producing VEB-1 and one producing CTX-M-14), but synergy between cefepime and amoxicillin-clavulanic acid was visible after the distance between the disks was reduced to 20 mm. After the interpretative guidelines of the EES were updated according to our results, the percentage of correct phenotype identification increased from 78 to 96%.

Genetic relationship between methicillin-sensitive and methicillin-resistant Staphylococcus aureus strains from France and from international sources: delineation of genomic groups.

Cluster analysis of the SmaI patterns, generated by pulsed-field gel electrophoresis, of 44 methicillin-resistant (MRSA) and 118 methicillin-sensitive (MSSA) Staphylococcus aureus strains isolated in various French hospitals and 61 MRSA and 48 MSSA strains from 20 other countries revealed 20 genomic groups distributed into four distantly related phylogenic branches. Eighty-three of the 105 MRSA strains (79%) were clustered in the six genomic groups of phylogenic branch I; and 154 of the 166 MSSA strains (92.8%) were clustered in the 14 genomic groups of phylogenic branches II, III, and IV. Agreement between genomic group and two other markers, esterase type and phage group, was obtained, emphasizing the clonal structure of the population. The genomic groups were delineated by esterase type. The distribution of the strains within the genomic groups was independent of their geographical origin; French strains were clustered with strains from other countries. The three types of the staphylococcal cassette chromosome mec (SCCmec) complex were distributed according to genomic groups. Most of the time, type I and type II SCCmec complexes were found in the MRSA strains belonging to the same genomic groups. In contrast, the type III SCCmec complex was specific to the MRSA strains belonging to the three genomic groups characterized by a common esterase type.

Prevalence, molecular epidemiology, and clinical significance of heterogeneous glycopeptide-intermediate Staphylococcus aureus in liver transplant recipients.

We investigated the prevalence, molecular epidemiology, and clinical significance of heterogeneous glycopeptide-intermediate Staphylococcus aureus (hGISA) isolates in 48 liver transplant recipients infected or colonized with methicillin-resistant S. aureus over a 5-year period. Strains were screened for hGISA on Mueller-Hinton agar containing 5 mg of teicoplanin per liter. Heterogeneous glycopeptide resistance was confirmed by the E-test method with a dense inoculum and a simplified method of population analysis. hGISA strains were found in 13 (27%) of the 48 patients studied. Eleven of the 13 strains shared a common multiresistant phenotype with homogeneous methicillin resistance and gentamicin resistance, and they were closely related according to the results of pulsed-field gel electrophoresis. Only 2 of the 13 patients infected or colonized with hGISA strains had previously received glycopeptide therapy. Most patients were successfully treated with vancomycin, but one patient who failed to respond to vancomycin subsequently died. These results suggest that the high prevalence of hGISA among our patients was due to the clonal spread of a multiresistant strain.

Evaluation of the Osiris expert system for identification of beta-lactam phenotypes in isolates of Pseudomonas aeruginosa.

Osiris is a video zone size reader for disk diffusion tests featuring a built-in extended expert system (EES). The efficacy of the EES for the identification of the beta-lactam susceptibility phenotypes of Pseudomonas aeruginosa isolates was evaluated. Thirteen beta-lactams were tested in four laboratories by the disk diffusion test with 53 strains with well-characterized resistance mechanisms, including the production of 12 extended-spectrum beta-lactamases (ESBLs). The plates were read with the Osiris system and the results were interpreted with the ESS, and then the phenotype identified by the EES was compared to the resistance mechanism. The strains were also screened for the presence of ESBL production by a double-disk synergy test by placing the strains between an extended-spectrum cephalosporin-containing disk and a clavulanic acid-containing disk at distances of 30, 20, 15, and 10 mm from each other. Overall, the EES accurately identified the phenotypes of 88.2% of the strains and indicated an association with several mechanisms for 3.8% of the strains. No phenotype was identified in four strains with low levels of penicillinase production. Misidentifications were observed for two penicillinase-producing strains: one strain with partially derepressed cephalosporinase production and one strain overexpressing the MexA-MexB-OprM efflux system. The production of only four ESBLs was detected by the standard synergy test with a 30-mm distance between the disks. The production of five further ESBLs was identified by reducing the distance to 20 mm, and the production of the last three ESBLs was detected only at a distance of 15 or 10 mm. Our results indicate that the Osiris EES is an effective tool for the identification of P. aeruginosa beta-lactam phenotypes. A specific double-disk synergy test with reduced disk distances is necessary for the detection of ESBL production by this organism.