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

Context-dependent structure formation of hairpin motifs in bacteriophage MS2 genomic RNA

Many functions of ribonucleic acid (RNA) rely on its ability to assume specific sequence-structure motifs. Packaging signals found in certain RNA viruses are one such prominent example of functional RNA motifs. These signals are short hairpin loops that interact with coat proteins and drive viral self-assembly. As they are found in different positions along the much longer genomic RNA, the formation of their correct structure occurs as a part of a larger context. Any changes to this context can consequently lead to changes in the structure of the motifs themselves. In fact, previous studies have shown that structure and function of RNA motifs can be highly context sensitive to the flanking sequence surrounding them. However, in what ways different flanking sequences influence the structure of an RNA motif they surround has yet to be studied in detail. We focus on a hairpin-rich region of the RNA genome of bacteriophage MS2-a well-studied RNA virus with a wide potential for use in biotechnology-and systematically examine context-dependent structural stability of 14 previously identified hairpin motifs, which include putative and confirmed packaging signals. Combining secondary and tertiary RNA structure prediction of the hairpin motifs placed in different contexts, ranging from the native genomic sequence to random RNA sequences and unstructured poly-U sequences, we determine different measures of motif structural stability. In this way, we show that while some motif structures can be stable in any context, others require specific context provided by the genome. Our results demonstrate the importance of context in RNA structure formation and how changes in the flanking sequence of an RNA motif sometimes lead to drastic changes in its structure. Structural stability of a motif in different contexts could provide additional insights into its functionality as well as assist in determining whether it remains functional when intentionally placed in other contexts.

Understanding species-specific and conserved RNA-protein interactions in vivo and in vitro

While evolution is often considered from a DNA- and protein-centric view, RNA-based regulation can also impact gene expression and protein sequences. Here we examine interspecies differences in RNA-protein interactions using the conserved neuronal RNA-binding protein, Unkempt (UNK) as model. We find that roughly half of mRNAs bound in human are also bound in mouse. Unexpectedly, even when transcript-level binding was conserved across species differential motif usage was prevalent. To understand the biochemical basis of UNK-RNA interactions, we reconstitute the human and mouse UNK-RNA interactomes using a high-throughput biochemical assay. We uncover detailed features driving binding, show that in vivo patterns are captured in vitro, find that highly conserved sites are the strongest bound, and associate binding strength with downstream regulation. Furthermore, subtle sequence differences surrounding motifs are key determinants of species-specific binding. We highlight the complex features driving protein-RNA interactions and how these evolve to confer species-specific regulation.

Genetic regulation of nascent RNA maturation revealed by direct RNA nanopore sequencing

Quantitative trait loci analyses have revealed an important role for genetic variants in regulating alternative splicing (AS) and alternative cleavage and polyadenylation (APA) in humans. Yet, these studies are generally performed with mature mRNA, so they report on the outcome rather than the processes of RNA maturation and thus may overlook how variants directly modulate pre-mRNA processing. The order in which the many introns of a human gene are removed can substantially influence AS, while nascent RNA polyadenylation can affect RNA stability and decay. However, how splicing order and poly(A) tail length are regulated by genetic variation has never been explored. Here, we used direct RNA nanopore sequencing to investigate allele-specific pre-mRNA maturation in 12 human lymphoblastoid cell lines. We found frequent splicing order differences between alleles and uncovered significant single nucleotide polymorphism (SNP)-splicing order associations in 17 genes. This included SNPs located in or near splice sites as well as more distal intronic and exonic SNPs. Moreover, several genes showed allele-specific poly(A) tail lengths, many of which also had a skewed allelic abundance ratio. HLA class I transcripts, which encode proteins that play an essential role in antigen presentation, showed the most allele-specific splicing orders, which frequently co-occurred with allele-specific AS, APA or poly(A) tail length differences. Together, our results expose new layers of genetic regulation of pre-mRNA maturation and highlight the power of long-read RNA sequencing for allele-specific analyses.

Big data and deep learning for RNA biology

The exponential growth of big data in RNA biology (RB) has led to the development of deep learning (DL) models that have driven crucial discoveries. As constantly evidenced by DL studies in other fields, the successful implementation of DL in RB depends heavily on the effective utilization of large-scale datasets from public databases. In achieving this goal, data encoding methods, learning algorithms, and techniques that align well with biological domain knowledge have played pivotal roles. In this review, we provide guiding principles for applying these DL concepts to various problems in RB by demonstrating successful examples and associated methodologies. We also discuss the remaining challenges in developing DL models for RB and suggest strategies to overcome these challenges. Overall, this review aims to illuminate the compelling potential of DL for RB and ways to apply this powerful technology to investigate the intriguing biology of RNA more effectively.

Differences in 5'untranslated regions highlight the importance of translational regulation of dosage sensitive genes

BACKGROUND

Untranslated regions (UTRs) are important mediators of post-transcriptional regulation. The length of UTRs and the composition of regulatory elements within them are known to vary substantially across genes, but little is known about the reasons for this variation in humans. Here, we set out to determine whether this variation, specifically in 5'UTRs, correlates with gene dosage sensitivity.

RESULTS

We investigate 5'UTR length, the number of alternative transcription start sites, the potential for alternative splicing, the number and type of upstream open reading frames (uORFs) and the propensity of 5'UTRs to form secondary structures. We explore how these elements vary by gene tolerance to loss-of-function (LoF; using the LOEUF metric), and in genes where changes in dosage are known to cause disease. We show that LOEUF correlates with 5'UTR length and complexity. Genes that are most intolerant to LoF have longer 5'UTRs, greater TSS diversity, and more upstream regulatory elements than their LoF tolerant counterparts. We show that these differences are evident in disease gene-sets, but not in recessive developmental disorder genes where LoF of a single allele is tolerated.

CONCLUSIONS

Our results confirm the importance of post-transcriptional regulation through 5'UTRs in tight regulation of mRNA and protein levels, particularly for genes where changes in dosage are deleterious and lead to disease. Finally, to support gene-based investigation we release a web-based browser tool, VuTR, that supports exploration of the composition of individual 5'UTRs and the impact of genetic variation within them.

Regulation by the RNA-binding protein Unkempt at its effector interface

How RNA-binding proteins (RBPs) convey regulatory instructions to the core effectors of RNA processing is unclear. Here, we document the existence and functions of a multivalent RBP-effector interface. We show that the effector interface of a conserved RBP with an essential role in metazoan development, Unkempt, is mediated by a novel type of 'dual-purpose' peptide motifs that can contact two different surfaces of interacting proteins. Unexpectedly, we find that the multivalent contacts do not merely serve effector recruitment but are required for the accuracy of RNA recognition by Unkempt. Systems analyses reveal that multivalent RBP-effector contacts can repurpose the principal activity of an effector for a different function, as we demonstrate for the reuse of the central eukaryotic mRNA decay factor CCR4-NOT in translational control. Our study establishes the molecular assembly and functional principles of an RBP-effector interface.

Understanding species-specific and conserved RNA-protein interactions in vivo and in vitro

While evolution is often considered from a DNA- and protein-centric view, RNA-based regulation can also impact gene expression and protein sequences. Here we examined interspecies differences in RNA-protein interactions using the conserved neuronal RNA binding protein, Unkempt (UNK) as model. We find that roughly half of mRNAs bound in human are also bound in mouse. Unexpectedly, even when transcript-level binding was conserved across species differential motif usage was prevalent. To understand the biochemical basis of UNK-RNA interactions, we reconstituted the human and mouse UNK-RNA interactomes using a high-throughput biochemical assay. We uncover detailed features driving binding, show that in vivo patterns are captured in vitro, find that highly conserved sites are the strongest bound, and associate binding strength with downstream regulation. Furthermore, subtle sequence differences surrounding motifs are key determinants of species-specific binding. We highlight the complex features driving protein-RNA interactions and how these evolve to confer species-specific regulation.

A paradigm for regulation at the effector interface with RNA-binding proteins

RNA-binding proteins (RBPs) are key regulators of gene expression, but how RBPs convey regulatory instructions to the core effectors of RNA processing is unclear. Here we document the existence and functions of a multivalent RBP-effector interface. We show that the effector interface of a deeply conserved RBP with an essential role in metazoan development, Unkempt, is mediated by a novel type of 'dual-purpose' peptide motifs that can contact two different surfaces of interacting proteins. Unexpectedly, we find that the multivalent contacts do not merely serve effector recruitment but are required for the accuracy of RNA recognition by the recruiting RBP. Systems analyses reveal that multivalent RBP-effector contacts can repurpose the principal activity of an effector for a different function, as we demonstrate for reuse of the central eukaryotic mRNA decay factor CCR4-NOT in translational control. Our study establishes the molecular assembly and functional principles of an RBP-effector interface, with implications for the evolution and function of RBP-operated regulatory networks.

Cooperation and Competition of RNA Secondary Structure and RNA-Protein Interactions in the Regulation of Alternative Splicing

The regulation of alternative splicing in eukaryotic cells is carried out through the coordinated action of a large number of factors, including RNA-binding proteins and RNA structure. The RNA structure influences alternative splicing by blocking cis-regulatory elements, or bringing them closer or farther apart. In combination with RNA-binding proteins, it generates transcript conformations that help to achieve the necessary splicing outcome. However, the binding of regulatory proteins depends on RNA structure and, vice versa, the formation of RNA structure depends on the interaction with regulators. Therefore, RNA structure and RNA-binding proteins are inseparable components of common regulatory mechanisms. This review highlights examples of alternative splicing regulation by RNA-binding proteins, the regulation through local and long-range RNA structures, as well as how these elements work together, cooperate, and compete.

CRISPR Based Programmable RNA Editing in Primary Hippocampal Neurons

Investigating the RNA regulation landscape primarily relies on understanding how RNA-protein interactions are governed in various cell types, including neurons. Analysis of RNA-protein interactions in physiological environments warrants the development of new tools that rely on RNA manipulation. Recently, a CRISPR-based RNA-editing tool (dCas13b-ADAR2DD ) was developed to mitigate disease-associated point mutations in cell lines. Here, we explored the targeted sequence editing potential of the tool (dCas13b-ADAR2DD system) by adapting it to manipulate RNA function to visualize RNA editing in primary hippocampal neurons. This two-component system includes a programmable guide RNA (gRNA) complementary to the target RNA and a catalytically dead version of the Cas13b enzyme fused to ADAR. The RNA editing protocol outlined in this article relies on gRNA-dependent targeting of the dCas13b-ADAR fusion protein to the mutant form of the Dendra2 transcript. Dendra2 is not required for intrinsic cellular functioning. It was ectopically expressed for fluorescent detection as a proof-of-principle demonstration of targeted RNA editing. We first abrogated the fluorescence of Dendra2 by introducing a nonsense mutation that precludes the formation of the functional protein. To visualize the efficacy of the RNA editing in neurons, we used the dCas13b-ADAR2DD system to edit specific nucleotides within the Dendra2 mRNA to restore the amino acid codes critical for Dendra2 fluorescence. This method lays the foundation for future studies on the dynamics of activity-induced RNA-protein interactions in neurons and can be extended to manipulate the endogenous RNome in diverse neuronal subtypes. Furthermore, this methodology will enable investigators to visualize the spatial and temporal resolution of RNA-protein interactions without altering the genomes via conventional methods. © 2023 Wiley Periodicals LLC. Support Protocol: Preparation of mouse primary hippocampal culture Basic Protocol: Targeted editing of RNA.

Pre-mRNA splicing and its cotranscriptional connections

Transcription of eukaryotic genes by RNA polymerase II (Pol II) yields RNA precursors containing introns that must be spliced out and the flanking exons ligated together. Splicing is catalyzed by a dynamic ribonucleoprotein complex called the spliceosome. Recent evidence has shown that a large fraction of splicing occurs cotranscriptionally as the RNA chain is extruded from Pol II at speeds of up to 5 kb/minute. Splicing is more efficient when it is tethered to the transcription elongation complex, and this linkage permits functional coupling of splicing with transcription. We discuss recent progress that has uncovered a network of connections that link splicing to transcript elongation and other cotranscriptional RNA processing events.

A motif in the 5'untranslated region of messenger RNAs regulates protein synthesis in a S6 kinase-dependent manner

The 5' untranslated regions (UTRs) in messenger RNAs (mRNAs) play an important role in the regulation of protein synthesis. We had previously identified a group of mRNAs that includes human semaphorin 7A (SEMA7A) whose translation is upregulated by the Erk/p90S6K pathway in human eosinophils, with a potential negative impact in asthma and airway inflammation. In the current study, we aimed to find a common 5'UTR regulatory cis-element, and determine its impact on protein synthesis. We identified a common and conserved 5'UTR motif GGCTG-[(C/G)T(C/G)]n-GCC that was present in this group of mRNAs. Mutations of the first two GG bases in this motif in SEMA7A 5'UTR led to a complete loss of S6K activity dependence for maximal translation. In conclusion, the newly identified 5'UTR motif present in SEMA7A has a critical role in regulating S6K-dependent protein synthesis.

Pre-mRNA splicing order is predetermined and maintains splicing fidelity across multi-intronic transcripts

Combinatorially, intron excision within a given nascent transcript could proceed down any of thousands of paths, each of which would expose different dynamic landscapes of cis-elements and contribute to alternative splicing. In this study, we found that post-transcriptional multi-intron splicing order in human cells is largely predetermined, with most genes spliced in one or a few predominant orders. Strikingly, these orders were conserved across cell types and stages of motor neuron differentiation. Introns flanking alternatively spliced exons were frequently excised last, after their neighboring introns. Perturbations to the spliceosomal U2 snRNA altered the preferred splicing order of many genes, and these alterations were associated with the retention of other introns in the same transcript. In one gene, early removal of specific introns was sufficient to induce delayed excision of three proximal introns, and this delay was caused by two distinct cis-regulatory mechanisms. Together, our results demonstrate that multi-intron splicing order in human cells is predetermined, is influenced by a component of the spliceosome and ensures splicing fidelity across long pre-mRNAs.

RNA conformational propensities determine cellular activity

Cellular processes are the product of interactions between biomolecules, which associate to form biologically active complexes1. These interactions are mediated by intermolecular contacts, which if disrupted, lead to alterations in cell physiology. Nevertheless, the formation of intermolecular contacts nearly universally requires changes in the conformations of the interacting biomolecules. As a result, binding affinity and cellular activity crucially depend both on the strength of the contacts and on the inherent propensities to form binding-competent conformational states2,3. Thus, conformational penalties are ubiquitous in biology and must be known in order to quantitatively model binding energetics for protein and nucleic acid interactions4,5. However, conceptual and technological limitations have hindered our ability to dissect and quantitatively measure how conformational propensities affect cellular activity. Here we systematically altered and determined the propensities for forming the protein-bound conformation of HIV-1 TAR RNA. These propensities quantitatively predicted the binding affinities of TAR to the RNA-binding region of the Tat protein and predicted the extent of HIV-1 Tat-dependent transactivation in cells. Our results establish the role of ensemble-based conformational propensities in cellular activity and reveal an example of a cellular process driven by an exceptionally rare and short-lived RNA conformational state.

Modeling CRISPR-Cas13d on-target and off-target effects using machine learning approaches

A major challenge in the application of the CRISPR-Cas13d system is to accurately predict its guide-dependent on-target and off-target effect. Here, we perform CRISPR-Cas13d proliferation screens and design a deep learning model, named DeepCas13, to predict the on-target activity from guide sequences and secondary structures. DeepCas13 outperforms existing methods to predict the efficiency of guides targeting both protein-coding and non-coding RNAs. Guides targeting non-essential genes display off-target viability effects, which are closely related to their on-target efficiencies. Choosing proper negative control guides during normalization mitigates the associated false positives in proliferation screens. We apply DeepCas13 to the guides targeting lncRNAs, and identify lncRNAs that affect cell viability and proliferation in multiple cell lines. The higher prediction accuracy of DeepCas13 over existing methods is extensively confirmed via a secondary CRISPR-Cas13d screen and quantitative RT-PCR experiments. DeepCas13 is freely accessible via http://deepcas13.weililab.org .

Artificial intelligence methods enhance the discovery of RNA interactions

Understanding how RNAs interact with proteins, RNAs, or other molecules remains a challenge of main interest in biology, given the importance of these complexes in both normal and pathological cellular processes. Since experimental datasets are starting to be available for hundreds of functional interactions between RNAs and other biomolecules, several machine learning and deep learning algorithms have been proposed for predicting RNA-RNA or RNA-protein interactions. However, most of these approaches were evaluated on a single dataset, making performance comparisons difficult. With this review, we aim to summarize recent computational methods, developed in this broad research area, highlighting feature encoding and machine learning strategies adopted. Given the magnitude of the effect that dataset size and quality have on performance, we explored the characteristics of these datasets. Additionally, we discuss multiple approaches to generate datasets of negative examples for training. Finally, we describe the best-performing methods to predict interactions between proteins and specific classes of RNA molecules, such as circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs), and methods to predict RNA-RNA or RNA-RBP interactions independently of the RNA type.

Meiotic cDNA libraries reveal gene truncations and mitochondrial proteins important for competitive fitness in Saccharomyces cerevisiae

Gametogenesis is an evolutionarily conserved developmental program whereby a diploid progenitor cell undergoes meiosis and cellular remodeling to differentiate into haploid gametes, the precursors for sexual reproduction. Even in the simple eukaryotic organism Saccharomyces cerevisiae, the meiotic transcriptome is very rich and complex, thereby necessitating new tools for functional studies. Here, we report the construction of 5 stage-specific, inducible complementary DNA libraries from meiotic cells that represent over 84% of the genes found in the budding yeast genome. We employed computational strategies to detect endogenous meiotic transcript isoforms as well as library-specific gene truncations. Furthermore, we developed a robust screening pipeline to test the effect of each complementary DNA on competitive fitness. Our multiday proof-of-principle time course revealed 877 complementary DNAs that were detrimental for competitive fitness when overexpressed. The list included mitochondrial proteins that cause dose-dependent disruption of cellular respiration as well as library-specific gene truncations that expose a dominant negative effect on competitive growth. Together, these high-quality complementary DNA libraries provide an important tool for systematically identifying meiotic genes, transcript isoforms, and protein domains that are important for a specific biological function.

Secondary structures in RNA synthesis, splicing and translation

Even though the functional role of mRNA molecules is primarily decided by the nucleotide sequence, several properties are determined by secondary structure conformations. Examples of secondary structures include long range interactions, hairpins, R-loops and G-quadruplexes and they are formed through interactions of non-adjacent nucleotides. Here, we discuss advances in our understanding of how secondary structures can impact RNA synthesis, splicing, translation and mRNA half-life. During RNA synthesis, secondary structures determine RNA polymerase II (RNAPII) speed, thereby influencing splicing. Splicing is also determined by RNA binding proteins and their binding rates are modulated by secondary structures. For the initiation of translation, secondary structures can control the choice of translation start site. Here, we highlight the mechanisms by which secondary structures modulate these processes, discuss advances in technologies to detect and study them systematically, and consider the roles of RNA secondary structures in disease.

Principles and pitfalls of high-throughput analysis of microRNA-binding thermodynamics and kinetics by RNA Bind-n-Seq

RNA Bind-n-Seq (RBNS) is a cost-effective, high-throughput method capable of identifying the sequence preferences of RNA-binding proteins and of qualitatively defining relative dissociation constants. Although RBNS is often described as an unbiased method, several factors may influence the outcome of the analysis. Here, we discuss these biases and present an analytical strategy to estimate absolute binding affinities from RBNS data, extend RBNS to kinetic studies, and develop a framework to compute relative association and dissociation rate constants. As proof of principle, we measured the equilibrium binding properties of mammalian Argonaute2 (AGO2) guided by eight microRNAs (miRNAs) and kinetic parameters for let-7a. The miRNA-binding site repertoires, dissociation constants, and kinetic parameters calculated from RBNS data using our methods correlate well with values measured by traditional ensemble and single-molecule approaches. Our data provide additional quantitative measurements for Argonaute-bound miRNA binding that should facilitate development of quantitative targeting rules for individual miRNAs.

POSTAR3: an updated platform for exploring post-transcriptional regulation coordinated by RNA-binding proteins

RNA-binding proteins (RBPs) play key roles in post-transcriptional regulation. Accurate identification of RBP binding sites in multiple cell lines and tissue types from diverse species is a fundamental endeavor towards understanding the regulatory mechanisms of RBPs under both physiological and pathological conditions. Our POSTAR annotation processes make use of publicly available large-scale CLIP-seq datasets and external functional genomic annotations to generate a comprehensive map of RBP binding sites and their association with other regulatory events as well as functional variants. Here, we present POSTAR3, an updated database with improvements in data collection, annotation infrastructure, and analysis that support the annotation of post-transcriptional regulation in multiple species including: we made a comprehensive update on the CLIP-seq and Ribo-seq datasets which cover more biological conditions, technologies, and species; we added RNA secondary structure profiling for RBP binding sites; we provided miRNA-mediated degradation events validated by degradome-seq; we included RBP binding sites at circRNA junction regions; we expanded the annotation of RBP binding sites, particularly using updated genomic variants and mutations associated with diseases. POSTAR3 is freely available at http://postar.ncrnalab.org.

The disease-associated proteins Drosophila Nab2 and Ataxin-2 interact with shared RNAs and coregulate neuronal morphology

Nab2 encodes the Drosophila melanogaster member of a conserved family of zinc finger polyadenosine RNA-binding proteins (RBPs) linked to multiple steps in post-transcriptional regulation. Mutation of the Nab2 human ortholog ZC3H14 gives rise to an autosomal recessive intellectual disability but understanding of Nab2/ZC3H14 function in metazoan nervous systems is limited, in part because no comprehensive identification of metazoan Nab2/ZC3H14-associated RNA transcripts has yet been conducted. Moreover, many Nab2/ZC3H14 functional protein partnerships remain unidentified. Here, we present evidence that Nab2 genetically interacts with Ataxin-2 (Atx2), which encodes a neuronal translational regulator, and that these factors coordinately regulate neuronal morphology, circadian behavior, and adult viability. We then present the first high-throughput identifications of Nab2- and Atx2-associated RNAs in Drosophila brain neurons using RNA immunoprecipitation-sequencing (RIP-Seq). Critically, the RNA interactomes of each RBP overlap, and Nab2 exhibits high specificity in its RNA associations in neurons in vivo, associating with a small fraction of all polyadenylated RNAs. The identities of shared associated transcripts (e.g., drk, me31B, stai) and of transcripts specific to Nab2 or Atx2 (e.g., Arpc2 and tea) promise insight into neuronal functions of, and genetic interactions between, each RBP. Consistent with prior biochemical studies, Nab2-associated neuronal RNAs are overrepresented for internal A-rich motifs, suggesting these sequences may partially mediate Nab2 target selection. These data support a model where Nab2 functionally opposes Atx2 in neurons, demonstrate Nab2 shares associated neuronal RNAs with Atx2, and reveal Drosophila Nab2 associates with a more specific subset of polyadenylated mRNAs than its polyadenosine affinity alone may suggest.

A Deep Exon Cryptic Splice Site Promotes Aberrant Intron Retention in a Von Willebrand Disease Patient

A translationally silent single nucleotide mutation in exon 44 (E44) of the von Willebrand factor (VWF) gene is associated with inefficient removal of intron 44 in a von Willebrand disease (VWD) patient. This intron retention (IR) event was previously attributed to reordered E44 secondary structure that sequesters the normal splice donor site. We propose an alternative mechanism: the mutation introduces a cryptic splice donor site that interferes with the function of the annotated site to favor IR. We evaluated both models using minigene splicing reporters engineered to vary in secondary structure and/or cryptic splice site content. Analysis of splicing efficiency in transfected K562 cells suggested that the mutation-generated cryptic splice site in E44 was sufficient to induce substantial IR. Mutations predicted to vary secondary structure at the annotated site also had modest effects on IR and shifted the balance of residual splicing between the cryptic site and annotated site, supporting competition among the sites. Further studies demonstrated that introduction of cryptic splice donor motifs at other positions in E44 did not promote IR, indicating that interference with the annotated site is context dependent. We conclude that mutant deep exon splice sites can interfere with proper splicing by inducing IR.

Sequence determinants as key regulators in gene expression of T cells

T cell homeostasis, T cell differentiation, and T cell effector function rely on the constant fine-tuning of gene expression. To alter the T cell state, substantial remodeling of the proteome is required. This remodeling depends on the intricate interplay of regulatory mechanisms, including post-transcriptional gene regulation. In this review, we discuss how the sequence of a transcript influences these post-transcriptional events. In particular, we review how sequence determinants such as sequence conservation, GC content, and chemical modifications define the levels of the mRNA and the protein in a T cell. We describe the effect of different forms of alternative splicing on mRNA expression and protein production, and their effect on subcellular localization. In addition, we discuss the role of sequences and structures as binding hubs for miRNAs and RNA-binding proteins in T cells. The review thus highlights how the intimate interplay of post-transcriptional mechanisms dictate cellular fate decisions in T cells.

RNAPhaSep: a resource of RNAs undergoing phase separation

Liquid-liquid phase separation (LLPS) partitions cellular contents, underlies the formation of membraneless organelles and plays essential biological roles. To date, most of the research on LLPS has focused on proteins, especially RNA-binding proteins. However, accumulating evidence has demonstrated that RNAs can also function as ‘scaffolds’ and play essential roles in seeding or nucleating the formation of granules. To better utilize the knowledge dispersed in published literature, we here introduce RNAPhaSep (http://www.rnaphasep.cn), a manually curated database of RNAs undergoing LLPS. It contains 1113 entries with experimentally validated RNA self-assembly or RNA and protein co-involved phase separation events. RNAPhaSep contains various types of information, including RNA information, protein information, phase separation experiment information and integrated annotation from multiple databases. RNAPhaSep provides a valuable resource for exploring the relationship between RNA properties and phase behaviour, and may further enhance our comprehensive understanding of LLPS in cellular functions and human diseases.

Spatial correlation statistics enable transcriptome-wide characterization of RNA structure binding

Molecular interactions at identical transcriptomic locations or at proximal but non-overlapping sites can mediate RNA modification and regulation, necessitating tools to uncover these spatial relationships. We present nearBynding, a flexible algorithm and software pipeline that models spatial correlation between transcriptome-wide tracks from diverse data types. nearBynding can process and correlate interval as well as continuous data and incorporate experimentally derived or in silico predicted transcriptomic tracks. nearBynding offers visualization functions for its statistics to identify colocalizations and adjacent features. We demonstrate the application of nearBynding to correlate RNA-binding protein (RBP) binding preferences with other RBPs, RNA structure, or RNA modification. By cross-correlating RBP binding and RNA structure data, we demonstrate that nearBynding recapitulates known RBP binding to structural motifs and provides biological insights into RBP binding preference of G-quadruplexes. nearBynding is available as an R/Bioconductor package and can run on a personal computer, making correlation of transcriptomic features broadly accessible.

Effects of sequence motifs in the yeast 3' untranslated region determined from massively parallel assays of random sequences

BACKGROUND

The 3' untranslated region (UTR) plays critical roles in determining the level of gene expression through effects on activities such as mRNA stability and translation. Functional elements within this region have largely been identified through analyses of native genes, which contain multiple co-evolved sequence features.

RESULTS

To explore the effects of 3' UTR sequence elements outside of native sequence contexts, we analyze hundreds of thousands of random 50-mers inserted into the 3' UTR of a reporter gene in the yeast Saccharomyces cerevisiae. We determine relative protein expression levels from the fitness of transformants in a growth selection. We find that the consensus 3' UTR efficiency element significantly boosts expression, independent of sequence context; on the other hand, the consensus positioning element has only a small effect on expression. Some sequence motifs that are binding sites for Puf proteins substantially increase expression in the library, despite these proteins generally being associated with post-transcriptional downregulation of native mRNAs. Our measurements also allow a systematic examination of the effects of point mutations within efficiency element motifs across diverse sequence backgrounds. These mutational scans reveal the relative in vivo importance of individual bases in the efficiency element, which likely reflects their roles in binding the Hrp1 protein involved in cleavage and polyadenylation.

CONCLUSIONS

The regulatory effects of some 3' UTR sequence features, like the efficiency element, are consistent regardless of sequence context. In contrast, the consequences of other 3' UTR features appear to be strongly dependent on their evolved context within native genes.

Matrix-screening reveals a vast potential for direct protein-protein interactions among RNA binding proteins

RNA-binding proteins (RBPs) are crucial factors of post-transcriptional gene regulation and their modes of action are intensely investigated. At the center of attention are RNA motifs that guide where RBPs bind. However, sequence motifs are often poor predictors of RBP-RNA interactions in vivo. It is hence believed that many RBPs recognize RNAs as complexes, to increase specificity and regulatory possibilities. To probe the potential for complex formation among RBPs, we assembled a library of 978 mammalian RBPs and used rec-Y2H matrix screening to detect direct interactions between RBPs, sampling > 600 K interactions. We discovered 1994 new interactions and demonstrate that interacting RBPs bind RNAs adjacently in vivo. We further find that the mRNA binding region and motif preferences of RBPs deviate, depending on their adjacently binding interaction partners. Finally, we reveal novel RBP interaction networks among major RNA processing steps and show that splicing impairing RBP mutations observed in cancer rewire spliceosomal interaction networks. The dataset we provide will be a valuable resource for understanding the combinatorial interactions of RBPs with RNAs and the resulting regulatory outcomes.

Discovery of a pre-mRNA structural scaffold as a contributor to the mammalian splicing code

The specific recognition of splice signals at or near exon-intron junctions is not explained by their weak conservation and instead is postulated to require a multitude of features embedded in the pre-mRNA strand. We explored the possibility of 3D structural scaffold of AdML-a model pre-mRNA substrate-guiding early spliceosomal components to the splice signal sequences. We find that mutations in the non-cognate splice signal sequences impede recruitment of early spliceosomal components due to disruption of the global structure of the pre-mRNA. We further find that the pre-mRNA segments potentially interacting with the early spliceosomal component U1 snRNP are distributed across the intron, that there is a spatial proximity of 5' and 3' splice sites within the pre-mRNA scaffold, and that an interplay exists between the structural scaffold and splicing regulatory elements in recruiting early spliceosomal components. These results suggest that early spliceosomal components can recognize a 3D structural scaffold beyond the short splice signal sequences, and that in our model pre-mRNA, this scaffold is formed across the intron involving the major splice signals. This provides a conceptual basis to analyze the contribution of recognizable 3D structural scaffolds to the splicing code across the mammalian transcriptome.

Affecting RNA biology genome-wide by binding small molecules and chemically induced proximity

The ENCODE and genome-wide association projects have shown that much of the genome is transcribed into RNA and much less is translated into protein. These and other functional studies suggest that the druggable transcriptome is much larger than the druggable proteome. This review highlights approaches to define druggable RNA targets and structure-activity relationships across genomic RNA. Binding compounds can be identified and optimized into structure-specific ligands by using sequence-based design with various modes of action, for example, inhibiting translation or directing pre-mRNA splicing outcomes. In addition, strategies to direct protein activity against an RNA of interest via chemically induced proximity is a burgeoning area that has been validated both in cells and in preclinical animal models, and we describe that it may allow rapid access to new avenues to affect RNA biology. These approaches and the unique modes of action suggest that more RNAs are potentially amenable to targeting than proteins.

Conserved long-range base pairings are associated with pre-mRNA processing of human genes

The ability of nucleic acids to form double-stranded structures is essential for all living systems on Earth. Current knowledge on functional RNA structures is focused on locally-occurring base pairs. However, crosslinking and proximity ligation experiments demonstrated that long-range RNA structures are highly abundant. Here, we present the most complete to-date catalog of conserved complementary regions (PCCRs) in human protein-coding genes. PCCRs tend to occur within introns, suppress intervening exons, and obstruct cryptic and inactive splice sites. Double-stranded structure of PCCRs is supported by decreased icSHAPE nucleotide accessibility, high abundance of RNA editing sites, and frequent occurrence of forked eCLIP peaks. Introns with PCCRs show a distinct splicing pattern in response to RNAPII slowdown suggesting that splicing is widely affected by co-transcriptional RNA folding. The enrichment of 3'-ends within PCCRs raises the intriguing hypothesis that coupling between RNA folding and splicing could mediate co-transcriptional suppression of premature pre-mRNA cleavage and polyadenylation.

Structural basis for RNA recognition by the N-terminal tandem RRM domains of human RBM45

RBM45 is an RNA-binding protein involved in neural development, whose aggregation is associated with neurodegenerative diseases, such as amyotrophic lateral sclerosis (ALS) and frontotemporal lobar dementia (FTLD). However, the mechanisms of RNA-binding and aggregation of RBM45 remain unelucidated. Here, we report the crystal structure of the N-terminal tandem RRM domains of human RBM45 in complex with single-stranded DNA (ssDNA). Our structural and biochemical results revealed that both the RRM1 and RRM2 of RBM45 recognized the GAC sequence of RNA/ssDNA. Two aromatic residues and an arginine residue in each RRM were critical for RNA-binding, and the interdomain linker was also involved in RNA-binding. Two RRMs formed a pair of antiparallel RNA-binding sites, indicating that the N-terminal tandem RRM domains of RBM45 bound separate GAC motifs in one RNA strand or GAC motifs in different RNA strands. Our findings will be helpful in the identification of physiologic targets of RBM45 and provide evidence for understanding the physiologic and pathologic functions of RBM45.

Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

Interaction of OIP5-AS1 with MEF2C mRNA promotes myogenic gene expression

Long noncoding (lnc)RNAs potently regulate gene expression programs in physiology and disease. Here, we describe a key function for lncRNA OIP5-AS1 in myogenesis, the process whereby myoblasts differentiate into myotubes during muscle development and muscle regeneration after injury. In human myoblasts, OIP5-AS1 levels increased robustly early in myogenesis, and its loss attenuated myogenic differentiation and potently reduced the levels of the myogenic transcription factor MEF2C. This effect relied upon the partial complementarity of OIP5-AS1 with MEF2C mRNA and the presence of HuR, an RNA-binding protein (RBP) with affinity for both transcripts. Remarkably, HuR binding to MEF2C mRNA, which stabilized MEF2C mRNA and increased MEF2C abundance, was lost after OIP5-AS1 silencing, suggesting that OIP5-AS1 might serve as a scaffold to enhance HuR binding to MEF2C mRNA, in turn increasing MEF2C production. These results highlight a mechanism whereby a lncRNA promotes myogenesis by enhancing the interaction of an RBP and a myogenic mRNA.

Identification of 3' UTR motifs required for mRNA localization to myelin sheaths in vivo

Myelin is a specialized membrane produced by oligodendrocytes that insulates and supports axons. Oligodendrocytes extend numerous cellular processes, as projections of the plasma membrane, and simultaneously wrap multiple layers of myelin membrane around target axons. Notably, myelin sheaths originating from the same oligodendrocyte are variable in size, suggesting local mechanisms regulate myelin sheath growth. Purified myelin contains ribosomes and hundreds of mRNAs, supporting a model that mRNA localization and local protein synthesis regulate sheath growth and maturation. However, the mechanisms by which mRNAs are selectively enriched in myelin sheaths are unclear. To investigate how mRNAs are targeted to myelin sheaths, we tested the hypothesis that transcripts are selected for myelin enrichment through consensus sequences in the 3' untranslated region (3' UTR). Using methods to visualize mRNA in living zebrafish larvae, we identified candidate 3' UTRs that were sufficient to localize mRNA to sheaths and enriched near growth zones of nascent membrane. We bioinformatically identified motifs common in 3' UTRs from 3 myelin-enriched transcripts and determined that these motifs are required and sufficient in a context-dependent manner for mRNA transport to myelin sheaths. Finally, we show that 1 motif is highly enriched in the myelin transcriptome, suggesting that this sequence is a global regulator of mRNA localization during developmental myelination.

Poison Exon Splicing Regulates a Coordinated Network of SR Protein Expression during Differentiation and Tumorigenesis

The RNA isoform repertoire is regulated by splicing factor (SF) expression, and alterations in SF levels are associated with disease. SFs contain ultraconserved poison exon (PE) sequences that exhibit greater identity across species than nearby coding exons, but their physiological role and molecular regulation is incompletely understood. We show that PEs in serine-arginine-rich (SR) proteins, a family of 14 essential SFs, are differentially spliced during induced pluripotent stem cell (iPSC) differentiation and in tumors versus normal tissues. We uncover an extensive cross-regulatory network of SR proteins controlling their expression via alternative splicing coupled to nonsense-mediated decay. We define sequences that regulate PE inclusion and protein expression of the oncogenic SF TRA2β using an RNA-targeting CRISPR screen. We demonstrate location dependency of RS domain activity on regulation of TRA2β-PE using CRISPR artificial SFs. Finally, we develop splice-switching antisense oligonucleotides to reverse the increased skipping of TRA2β-PE detected in breast tumors, altering breast cancer cell viability, proliferation, and migration.

Principles of mRNA control by human PUM proteins elucidated from multimodal experiments and integrative data analysis

The human PUF-family proteins, PUM1 and PUM2, posttranscriptionally regulate gene expression by binding to a PUM recognition element (PRE) in the 3'-UTR of target mRNAs. Hundreds of PUM1/2 targets have been identified from changes in steady-state RNA levels; however, prior studies could not differentiate between the contributions of changes in transcription and RNA decay rates. We applied metabolic labeling to measure changes in RNA turnover in response to depletion of PUM1/2, showing that human PUM proteins regulate expression almost exclusively by changing RNA stability. We also applied an in vitro selection workflow to precisely identify the binding preferences of PUM1 and PUM2. By integrating our results with prior knowledge, we developed a "rulebook" of key contextual features that differentiate functional versus nonfunctional PREs, allowing us to train machine learning models that accurately predict the functional regulation of RNA targets by the human PUM proteins.

Two-step mechanism for selective incorporation of lncRNA into a chromatin modifier

The MLE DExH helicase and the roX lncRNAs are essential components of the chromatin modifying Dosage Compensation Complex (DCC) in Drosophila. To explore the mechanism of ribonucleoprotein complex assembly, we developed vitRIP, an unbiased, transcriptome-wide in vitro assay that reveals RNA binding specificity. We found that MLE has intrinsic specificity for U-/A-rich sequences and tandem stem-loop structures and binds many RNAs beyond roX in vitro. The selectivity of the helicase for physiological substrates is further enhanced by the core DCC. Unwinding of roX2 by MLE induces a highly selective RNA binding surface in the unstructured C-terminus of the MSL2 subunit and triggers-specific association of MLE and roX2 with the core DCC. The exquisite selectivity of roX2 incorporation into the DCC thus originates from intimate cooperation between the helicase and the core DCC involving two distinct RNA selection principles and their mutual refinement.

RNA structures in alternative splicing and back-splicing

Alternative splicing greatly expands the transcriptomic and proteomic diversities related to physiological and developmental processes in higher eukaryotes. Splicing of long noncoding RNAs, and back- and trans- splicing further expanded the regulatory repertoire of alternative splicing. RNA structures were shown to play an important role in regulating alternative splicing and back-splicing. Application of novel sequencing technologies made it possible to identify genome-wide RNA structures and interaction networks, which might provide new insights into RNA splicing regulation in vitro to in vivo. The emerging transcription-folding-splicing paradigm is changing our understanding of RNA alternative splicing regulation. Here, we review the insights into the roles and mechanisms of RNA structures in alternative splicing and back-splicing, as well as how disruption of these structures affects alternative splicing and then leads to human diseases. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Structural disruption of exonic stem-loops immediately upstream of the intron regulates mammalian splicing

Recognition of highly degenerate mammalian splice sites by the core spliceosomal machinery is regulated by several protein factors that predominantly bind exonic splicing motifs. These are postulated to be single-stranded in order to be functional, yet knowledge of secondary structural features that regulate the exposure of exonic splicing motifs across the transcriptome is not currently available. Using transcriptome-wide RNA structural information we show that retained introns in mouse are commonly flanked by a short (≲70 nucleotide), highly base-paired segment upstream and a predominantly single-stranded exonic segment downstream. Splicing assays with select pre-mRNA substrates demonstrate that loops immediately upstream of the introns contain pre-mRNA-specific splicing enhancers, the substitution or hybridization of which impedes splicing. Additionally, the exonic segments flanking the retained introns appeared to be more enriched in a previously identified set of hexameric exonic splicing enhancer (ESE) sequences compared to their spliced counterparts, suggesting that base-pairing in the exonic segments upstream of retained introns could be a means for occlusion of ESEs. The upstream exonic loops of the test substrate promoted recruitment of splicing factors and consequent pre-mRNA structural remodeling, leading up to assembly of the early spliceosome. These results suggest that disruption of exonic stem-loop structures immediately upstream (but not downstream) of the introns regulate alternative splicing events, likely through modulating accessibility of splicing factors.

Concentration-dependent splicing is enabled by Rbfox motifs of intermediate affinity

The Rbfox family of splicing factors regulate alternative splicing during animal development and in disease, impacting thousands of exons in the maturing brain, heart, and muscle. Rbfox proteins have long been known to bind to the RNA sequence GCAUG with high affinity, but just half of Rbfox binding sites contain a GCAUG motif in vivo. We incubated recombinant RBFOX2 with over 60,000 mouse and human transcriptomic sequences to reveal substantial binding to several moderate-affinity, non-GCAYG sites at a physiologically relevant range of RBFOX concentrations. We find that many of these “secondary motifs” bind Rbfox robustly in cells and that several together can exert regulation comparable to GCAUG in a trichromatic splicing reporter assay. Furthermore, secondary motifs regulate RNA splicing in neuronal development and in neuronal subtypes where cellular Rbfox concentrations are highest, enabling a second wave of splicing changes as Rbfox levels increase.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

Structural basis for mRNA recognition by human RBM38

RNA-binding protein RBM38 was reported to bind the mRNA of several p53-related genes through its RRM domain and to up-regulate or down-regulate protein translation by increasing mRNA stability or recruitment of other effector proteins. The recognition mechanism, however, for RNA-binding of RBM38 remains unclear. Here, we report the crystal structure of the RRM domain of human RBM38 in complex with a single-stranded RNA. Our structural and biological results revealed that RBM38 recognizes G(U/C/A)GUG sequence single-stranded RNA in a sequence-specific and structure-specific manner. Two phenylalanine stacked with bases of RNA were crucial for RNA binding, and a series of hydrogen bonds between the base atoms of RNA and main-chain or side-chain atoms of RBM38 determine the sequence-specific recognition. Our results revealed the RNA-recognition mechanism of human RBM38 and provided structural information for understanding the RNA-binding property of RBM38.

FMRP promotes RNA localization to neuronal projections through interactions between its RGG domain and G-quadruplex RNA sequences

The sorting of RNA molecules to subcellular locations facilitates the activity of spatially restricted processes. We have analyzed subcellular transcriptomes of FMRP-null mouse neuronal cells to identify transcripts that depend on FMRP for efficient transport to neurites. We found that these transcripts contain an enrichment of G-quadruplex sequences in their 3' UTRs, suggesting that FMRP recognizes them to promote RNA localization. We observed similar results in neurons derived from Fragile X Syndrome patients. We identified the RGG domain of FMRP as important for binding G-quadruplexes and the transport of G-quadruplex-containing transcripts. Finally, we found that the translation and localization targets of FMRP were distinct and that an FMRP mutant that is unable to bind ribosomes still promoted localization of G-quadruplex-containing messages. This suggests that these two regulatory modes of FMRP may be functionally separated. These results provide a framework for the elucidation of similar mechanisms governed by other RNA-binding proteins.

High throughput approaches to study RNA-protein interactions in vitro

To understand the regulation of gene expression it is critical to determine how proteins interact with and discriminate between different RNAs. In this review, we discuss experimental techniques that utilize high throughput approaches to characterize the interactions of proteins with large numbers of RNAs in vitro. We describe the underlying principles for the main methods, briefly discuss their scope and limitations, and outline how insight from the techniques contributes to our understanding of specificity for RNA-protein interactions.

RNA-protein coevolution study of Gemin5 uncovers the role of the PXSS motif of RBS1 domain for RNA binding

Regulation of protein synthesis is an essential step of gene expression. This process is under the control of cis-acting RNA elements and trans-acting factors. Gemin5 is a multifunctional RNA-binding protein organized in distinct domains. The protein bears a non-canonical RNA-binding site, designated RBS1, at the C-terminal end. Among other cellular RNAs, the RBS1 region recognizes a sequence located within the coding region of Gemin5 mRNA, termed H12. Expression of RBS1 stimulates translation of RNA reporters carrying the H12 sequence, counteracting the negative effect of Gemin5 on global protein synthesis. A computational analysis of RBS1 protein and H12 RNA variability across the evolutionary scale predicts coevolving pairs of amino acids and nucleotides. RBS1 footprint and gel-shift assays indicated a positive correlation between the identified coevolving pairs and RNA-protein interaction. The coevolving residues of RBS1 contribute to the recognition of stem-loop SL1, an RNA structural element of H12 that contains the coevolving nucleotides. Indeed, RBS1 proteins carrying substitutions on the coevolving residues P₁₂₉₇ or S₁₂₉₉S₁₃₀₀, drastically reduced SL1-binding. Unlike the wild type RBS1 protein, expression of these mutant proteins in cells failed to enhance translation stimulation of mRNA reporters carrying the H12 sequence. Therefore, the PXSS motif within the RBS1 domain of Gemin5 and the RNA structural motif SL1 of its mRNA appears to play a key role in fine-tuning the expression level of this essential protein.

3'UTRs Regulate Mouse Ntrk2 mRNA Distribution in Cortical Neurons

There are two major isoforms of NTRK2 (neurotrophic receptor tyrosine kinase 2, or TrkB), full-length isoform with tyrosine kinase (TK) domain intact (+) and spliced isoform without tyrosine kinase domain (TK(−)). Within each isoform, there exist subtypes with minor modifications of the protein sequences. In human, the NTRK2 mRNA transcripts encoding TK(+) have same 3′UTRs, while the transcripts encoding subtypes of NTRK2 TK(−) have two completely different 3′UTRs. In mouse, the mRNA transcripts encoding same NTRK2 protein sequence for either TK(+) or TK(−) have long or short 3′UTRs, respectively. The physiological functions of these different 3′UTRs are still unknown. Pilocarpine stimulation increased Ntrk2 mRNA levels in soma, while the increase in synaptosome was smaller. FISH results further showed that mouse Ntrk2 transcripts with different 3′UTRs were distributed differently in cultured cortical neurons. The transcripts with long 3′UTR were distributed more in apical dendrites compared with transcripts with short 3′UTR. Our results provide evidence of non-coding 3′UTR function in regulating mRNA distribution in neurons.

RNA structural dynamics regulate early embryogenesis through controlling transcriptome fate and function

BACKGROUND

Vertebrate early embryogenesis is initially directed by a set of maternal RNAs and proteins, yet the mechanisms controlling this program remain largely unknown. Recent transcriptome-wide studies on RNA structure have revealed its pervasive and crucial roles in RNA processing and functions, but whether and how RNA structure regulates the fate of the maternal transcriptome have yet to be determined.

RESULTS

Here we establish the global map of four nucleotide-based mRNA structures by icSHAPE during zebrafish early embryogenesis. Strikingly, we observe that RNA structurally variable regions are enriched in the 3' UTR and contain cis-regulatory elements important for maternal-to-zygotic transition (MZT). We find that the RNA-binding protein Elavl1a stabilizes maternal mRNAs by binding to the cis-elements. Conversely, RNA structure formation suppresses Elavl1a's binding leading to the decay of its maternal targets.

CONCLUSIONS

Our study finds that RNA structurally variable regions are enriched in mRNA 3' UTRs and contain cis-regulatory elements during zebrafish early embryogenesis. We reveal that Elavl1a regulates maternal RNA stability in an RNA structure-dependent fashion. Overall, our findings reveal a broad and fundamental role of RNA structure-based regulation in vertebrate early embryogenesis.

dSreg: a Bayesian model to integrate changes in splicing and RNA-binding protein activity

MOTIVATION

Alternative splicing (AS) is an important mechanism in the generation of transcript diversity across mammals. AS patterns are dynamically regulated during development and in response to environmental changes. Defects or perturbations in its regulation may lead to cancer or neurological disorders, among other pathological conditions. The regulatory mechanisms controlling AS in a given biological context are typically inferred using a two-step framework: differential AS analysis followed by enrichment methods. These strategies require setting rather arbitrary thresholds and are prone to error propagation along the analysis.

RESULTS

To overcome these limitations, we propose dSreg, a Bayesian model that integrates RNA-seq with data from regulatory features, e.g. binding sites of RNA-binding proteins. dSreg identifies the key underlying regulators controlling AS changes and quantifies their activity while simultaneously estimating the changes in exon inclusion rates. dSreg increased both the sensitivity and the specificity of the identified AS changes in simulated data, even at low read coverage. dSreg also showed improved performance when analyzing a collection of knock-down RNA-binding proteins' experiments from ENCODE, as opposed to traditional enrichment methods, such as over-representation analysis and gene set enrichment analysis. dSreg opens the possibility to integrate a large amount of readily available RNA-seq datasets at low coverage for AS analysis and allows more cost-effective RNA-seq experiments.

AVAILABILITY AND IMPLEMENTATION

dSreg was implemented in python using stan and is freely available to the community at https://bitbucket.org/cmartiga/dsreg.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The roles of structural dynamics in the cellular functions of RNAs

RNAs fold into 3D structures that range from simple helical elements to complex tertiary structures and quaternary ribonucleoprotein assemblies. The functions of many regulatory RNAs depend on how their 3D structure changes in response to a diverse array of cellular conditions. In this Review, we examine how the structural characterization of RNA as dynamic ensembles of conformations, which form with different probabilities and at different timescales, is improving our understanding of RNA function in cells. We discuss the mechanisms of gene regulation by microRNAs, riboswitches, ribozymes, post-transcriptional RNA modifications and RNA-binding proteins, and how the cellular environment and processes such as liquid-liquid phase separation may affect RNA folding and activity. The emerging RNA-ensemble-function paradigm is changing our perspective and understanding of RNA regulation, from in vitro to in vivo and from descriptive to predictive.