Towards Long-Range RNA Structure Prediction in Eukaryotic Genes

Comprehensive Mapping of Human dsRNAome Reveals Conservation, Neuronal Enrichment, and Intermolecular Interactions

The human transcriptome contains millions of A-to-I editing sites arising from an unclear number of poorly characterized dsRNAs. Editing sites are often used to infer presence of dsRNA, but this method is limited by transcription levels, read depth, ADAR expression and cannot identify unedited dsRNA. To address these limitations, we developed dsRNAscan. Applying dsRNAscan to the human genome predicted 5 million dsRNAs. Genomic distribution of dsRNAs encompassing repetitive elements was widespread, but non-repetitive dsRNAs were sparse and enriched at chromosomal tips. dsRNAscan predicted hundreds of long, highly paired dsRNAs suspected to be immunogenic, but only one was in a 3'UTR, and thus likely to challenge cytoplasmic immune sensors. We observed several thousand editing enriched regions suspected to arise from intermolecular structures, and dozens of neuronally enriched dsRNAs conserved across vertebrates. This study offers the first comprehensive set of dsRNA annotations for the human genome, available as a resource at https://dsrna.chpc.utah.edu/.

BRD2 and BRD3 genes independently evolved RNA structures to control unproductive splicing

The mammalian BRD2 and BRD3 genes encode structurally related proteins from the bromodomain and extraterminal domain protein family. The expression of BRD2 is regulated by unproductive splicing upon inclusion of exon 3b, which is located in the region encoding a bromodomain. Bioinformatic analysis indicated that BRD2 exon 3b inclusion is controlled by a pair of conserved complementary regions (PCCR) located in the flanking introns. Furthermore, we identified a highly conserved element encoding a cryptic poison exon 5b and a previously unknown PCCR in the intron between exons 5 and 6 of BRD3, however, outside of the homologous bromodomain. Minigene mutagenesis and blockage of RNA structure by antisense oligonucleotides demonstrated that RNA structure controls the rate of inclusion of poison exons. The patterns of BRD2 and BRD3 expression and splicing show downregulation upon inclusion of poison exons, which become skipped in response to transcription elongation slowdown, further confirming a role of PCCRs in unproductive splicing regulation. We conclude that BRD2 and BRD3 independently acquired poison exons and RNA structures to dynamically control unproductive splicing. This study describes a convergent evolution of regulatory unproductive splicing mechanisms in these genes, providing implications for selective modulation of their expression in therapeutic applications.

Long-range RNA structures in the human transcriptome beyond evolutionarily conserved regions

RNA structure has been increasingly recognized as a critical player in the biogenesis and turnover of many transcripts classes. In eukaryotes, the prediction of RNA structure by thermodynamic modeling meets fundamental limitations due to the large sizes and complex, discontinuous organization of eukaryotic genes. Signatures of functional RNA structures can be found by detecting compensatory substitutions in homologous sequences, but a comparative approach is applicable only within conserved sequence blocks. Here, we developed a computational pipeline called PHRIC, which is not limited to conserved regions and relies on RNA contacts derived from RNA in situ conformation sequencing (RIC-seq) experiments. It extracts pairs of short RNA fragments surrounded by nested clusters of RNA contacts and predicts long, nearly perfect complementary base pairings formed between these fragments. In application to a panel of RIC-seq experiments in seven human cell lines, PHRIC predicted ~12,000 stable long-range RNA structures with equilibrium free energy below -15 kcal/mol, the vast majority of which fall outside of regions annotated as conserved among vertebrates. These structures, nevertheless, show some level of sequence conservation and remarkable compensatory substitution patterns in other clades. Furthermore, we found that introns have a higher propensity to form stable long-range RNA structures between each other, and moreover that RNA structures tend to concentrate within the same intron rather than connect adjacent introns. These results for the first time extend the application of proximity ligation assays to RNA structure prediction beyond conserved regions.

A Hitchhiker's guide to RNA-RNA structure and interaction prediction tools

RNA biology has risen to prominence after a remarkable discovery of diverse functions of noncoding RNA (ncRNA). Most untranslated transcripts often exert their regulatory functions into RNA-RNA complexes via base pairing with complementary sequences in other RNAs. An interplay between RNAs is essential, as it possesses various functional roles in human cells, including genetic translation, RNA splicing, editing, ribosomal RNA maturation, RNA degradation and the regulation of metabolic pathways/riboswitches. Moreover, the pervasive transcription of the human genome allows for the discovery of novel genomic functions via RNA interactome investigation. The advancement of experimental procedures has resulted in an explosion of documented data, necessitating the development of efficient and precise computational tools and algorithms. This review provides an extensive update on RNA-RNA interaction (RRI) analysis via thermodynamic- and comparative-based RNA secondary structure prediction (RSP) and RNA-RNA interaction prediction (RIP) tools and their general functions. We also highlighted the current knowledge of RRIs and the limitations of RNA interactome mapping via experimental data. Then, the gap between RSP and RIP, the importance of RNA homologues, the relationship between pseudoknots, and RNA folding thermodynamics are discussed. It is hoped that these emerging prediction tools will deepen the understanding of RNA-associated interactions in human diseases and hasten treatment processes.

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.

RNA in situ conformation sequencing reveals novel long-range RNA structures with impact on splicing

Over recent years, long-range RNA structure has emerged as a factor that is fundamental to alternative splicing regulation. An increasing number of human disorders are now being associated with splicing defects; hence it is essential to develop methods that assess long-range RNA structure experimentally. RNA in situ conformation sequencing (RIC-seq) is a method that recapitulates RNA structure within physiological RNA-protein complexes. In this work, we juxtapose pairs of conserved complementary regions (PCCRs) that were predicted in silico with the results of RIC-seq experiments conducted in seven human cell lines. We show statistically that RIC-seq support of PCCRs correlates with their properties, such as equilibrium free energy, presence of compensatory substitutions, and occurrence of A-to-I RNA editing sites and forked eCLIP peaks. Exons enclosed in PCCRs that are supported by RIC-seq tend to have weaker splice sites and lower inclusion rates, which is indicative of post-transcriptional splicing regulation mediated by RNA structure. Based on these findings, we prioritize PCCRs according to their RIC-seq support and show, using antisense nucleotides and minigene mutagenesis, that PCCRs in two disease-associated human genes, PHF20L1 and CASK, and also PCCRs in their murine orthologs, impact alternative splicing. In sum, we demonstrate how RIC-seq experiments can be used to discover functional long-range RNA structures, and particularly those that regulate alternative splicing.

Exploring structural determinants and the role of nucleolin in formation of the long-range interactions between untranslated regions of p53 mRNA

p53 protein is a key regulator of cellular homeostasis by coordinating the framework of antiproliferative pathways as a response to various stress factors. Although the main mechanism of stress-dependent induction of p53 protein relies on post-translational modifications influencing its stability and activity, a growing amount of evidence suggests that complex regulation of p53 expression occurs also at the mRNA level. This study explores structural determinants of long-range RNA-RNA interactions in p53 mRNA, crucial for stress-dependent regulation of p53 protein translation. We demonstrate that the 8-nt bulge motif plays a key structural role in base-pairing of complementary sequences from the 5' and 3' untranslated regions of p53 mRNA. We also show that one of the p53 translation regulators, nucleolin, displays an RNA chaperone activity and facilitates the association of sequences involved in the formation of long-range interactions in p53 mRNA. Nucleolin promotes base-pairing of complementary sequences through the bulge motif, because mutations of this region reduce or inhibit pairing while compensatory mutations restore this interaction. Mutational analysis of nucleolin reveals that all four RNA recognition motifs are indispensable for optimal RNA chaperone activity of nucleolin. These observations help to decipher the unique mechanism of p53 protein translation regulation pointing to bulge motif and nucleolin as the critical factors during intramolecular RNA-RNA recognition in p53 mRNA.

RNAsmc: A integrated tool for comparing RNA secondary structure and evaluating allosteric effects

RNA structure plays a crucial role in gene regulation, in RNA stability and the essential biological processes. RNA secondary structure (RSS) motifs are the basic building blocks for investigating the biological mechanisms of structure. Here, we present a strategy for structural motif-based dynamic alignment, namely, RNA secondary-structural motif-comparing (RNAsmc), to identify structural motifs and quantitatively evaluate their underlying molecular functions. RNAsmc also has strong robustness to sequence length, folding protocol and RNA structural profile by chemical probing. Notably, it is also applicable to quantify structural variation in special RNA editing events (SNVs or SNPs, fragment insertion or deletion, etc.). The findings indicate that RNAsmc can uncover the heterogeneity of RNA secondary structure and score for similarities among components, which provides an impetus to cluster RNA families and evaluate allosteric effects. We find that RNAsmc exhibits remarkable detection efficiency for experimentally-derived RiboSNitches. Finally, the pipeline was assembled into an R software package to serve as an automated toolkit to explore, align, and cluster RSS. It is freely available for download at https://CRAN.R-project.org/package=RNAsmc.

Structural Context of a Critical Exon of Spinal Muscular Atrophy Gene

Humans contain two nearly identical copies of Survival Motor Neuron genes, SMN1 and SMN2. Deletion or mutation of SMN1 causes spinal muscular atrophy (SMA), one of the leading genetic diseases associated with infant mortality. SMN2 is unable to compensate for the loss of SMN1 due to predominant exon 7 skipping, leading to the production of a truncated protein. Antisense oligonucleotide and small molecule-based strategies aimed at the restoration of SMN2 exon 7 inclusion are approved therapies of SMA. Many cis-elements and transacting factors have been implicated in regulation of SMN exon 7 splicing. Also, several structural elements, including those formed by a long-distance interaction, have been implicated in the modulation of SMN exon 7 splicing. Several of these structures have been confirmed by enzymatic and chemical structure-probing methods. Additional structures formed by inter-intronic interactions have been predicted by computational algorithms. SMN genes generate a vast repertoire of circular RNAs through inter-intronic secondary structures formed by inverted Alu repeats present in large number in SMN genes. Here, we review the structural context of the exonic and intronic cis-elements that promote or prevent exon 7 recognition. We discuss how structural rearrangements triggered by single nucleotide substitutions could bring drastic changes in SMN2 exon 7 splicing. We also propose potential mechanisms by which inter-intronic structures might impact the splicing outcomes.

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.

Multiple competing RNA structures dynamically control alternative splicing in the human ATE1 gene

The mammalian Ate1 gene encodes an arginyl transferase enzyme with tumor suppressor function that depends on the inclusion of one of the two mutually exclusive exons (MXE), exons 7a and 7b. We report that the molecular mechanism underlying MXE splicing in Ate1 involves five conserved regulatory intronic elements R1-R5, of which R1 and R4 compete for base pairing with R3, while R2 and R5 form an ultra-long-range RNA structure spanning 30 Kb. In minigenes, single and double mutations that disrupt base pairings in R1R3 and R3R4 lead to the loss of MXE splicing, while compensatory triple mutations that restore RNA structure revert splicing to that of the wild type. In the endogenous Ate1 pre-mRNA, blocking the competing base pairings by LNA/DNA mixmers complementary to R3 leads to the loss of MXE splicing, while the disruption of R2R5 interaction changes the ratio of MXE. That is, Ate1 splicing is controlled by two independent, dynamically interacting, and functionally distinct RNA structure modules. Exon 7a becomes more included in response to RNA Pol II slowdown, however it fails to do so when the ultra-long-range R2R5 interaction is disrupted, indicating that exon 7a/7b ratio depends on co-transcriptional RNA folding. In sum, these results demonstrate that splicing is coordinated both in time and in space over very long distances, and that the interaction of these components is mediated by RNA structure.

Combinatorial metabolic pathway assembly approaches and toolkits for modular assembly

Synthetic Biology is a rapidly growing interdisciplinary field that is primarily built upon foundational advances in molecular biology combined with engineering design principles such as modularity and interoperability. The field considers living systems as programmable at the genetic level and has been defined by the development of new platform technologies and methodological advances. A key concept driving the field is the Design-Build-Test-Learn cycle which provides a systematic framework for building new biological systems. One major application area for synthetic biology is biosynthetic pathway engineering that requires the modular assembly of different genetic regulatory elements and biosynthetic enzymes. In this review we provide an overview of modular DNA assembly and describe and compare the plethora of in vitro and in vivo assembly methods for combinatorial pathway engineering. Considerations for part design and methods for enzyme balancing are also presented, and we briefly discuss alternatives to intracellular pathway assembly including microbial consortia and cell-free systems for biosynthesis. Finally, we describe computational tools and automation for pathway design and assembly and argue that a deeper understanding of the many different variables of genetic design, pathway regulation and cellular metabolism will allow more predictive pathway design and engineering.

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.

The regulation properties of RNA secondary structure in alternative splicing

The RNA secondary structure is important for many functional processes in the cell. The secondary and tertiary structures of cellular RNAs are essential for the activity of these molecules in processes such as transcription, splicing, translation, and localization. New high-throughput analytical methods, including next generation sequencing, have allowed for the in-depth characterization of the 'RNA structurome': a new term describing how the RNA structure controls the activity of RNA by itself and how it regulates the expression of genes. In this review, we present many examples of the influence of structural motifs of RNA, long range interactions and global RNA structure on the alternative splicing processes. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Integrative transcriptomic analysis suggests new autoregulatory splicing events coupled with nonsense-mediated mRNA decay

Nonsense-mediated decay (NMD) is a eukaryotic mRNA surveillance system that selectively degrades transcripts with premature termination codons (PTC). Many RNA-binding proteins (RBP) regulate their expression levels by a negative feedback loop, in which RBP binds its own pre-mRNA and causes alternative splicing to introduce a PTC. We present a bioinformatic analysis integrating three data sources, eCLIP assays for a large RBP panel, shRNA inactivation of NMD pathway, and shRNA-depletion of RBPs followed by RNA-seq, to identify novel such autoregulatory feedback loops. We show that RBPs frequently bind their own pre-mRNAs, their exons respond prominently to NMD pathway disruption, and that the responding exons are enriched with nearby eCLIP peaks. We confirm previously proposed models of autoregulation in SRSF7 and U2AF1 genes and present two novel models, in which (i) SFPQ binds its mRNA and promotes switching to an alternative distal 3'-UTR that is targeted by NMD, and (ii) RPS3 binding activates a poison 5'-splice site in its pre-mRNA that leads to a frame shift and degradation by NMD. We also suggest specific splicing events that could be implicated in autoregulatory feedback loops in RBM39, HNRNPM, and U2AF2 genes. The results are available through a UCSC Genome Browser track hub.

Circular exonic RNAs: When RNA structure meets topology

Although RNA circularization was first documented in the 1990s, the extent to which it occurs was not known until recent advances in high-throughput sequencing enabled the widespread identification of circular RNAs (circRNAs). Despite this, many aspects of circRNA biogenesis, structure, and function yet remain obscure. This review focuses on circular exonic RNAs, a subclass of circRNAs that are generated through backsplicing. Here, I hypothesize that RNA secondary structure can be the common factor that promotes both exon skipping and spliceosomal RNA circularization, and that backsplicing of double-stranded regions could generate topologically linked circRNA molecules. CircRNAs manifest themselves by the presence of tail-to-head exon junctions, which were previously attributed to post-transcriptional exon permutation and repetition. I revisit these observations and argue that backsplicing does not automatically imply RNA circularization because tail-to-head exon junctions give only local information about transcript architecture and, therefore, they are in principle insufficient to determine globally circular topology. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Computational approaches for the discovery of splicing regulatory RNA structures

Global RNA structure and local functional motifs mediate interactions important in determining the rates and patterns of mRNA splicing. In this review, we overview approaches for the computational prediction of RNA secondary structure with a special emphasis on the discovery of motifs important to RNA splicing. The process of identifying and modeling potential splicing regulatory structures is illustrated using a recently-developed approach for RNA structural motif discovery, the ScanFold pipeline, which is applied to the identification of a known splicing regulatory structure in influenza virus.

Special Issue: Computational Analysis of RNA Structure and Function

RNA structure often plays a key role in determining the function of non-coding and coding transcripts [...].

An Evolutionary Mechanism for the Generation of Competing RNA Structures Associated with Mutually Exclusive Exons

Alternative splicing is a commonly-used mechanism of diversifying gene products. Mutually exclusive exons (MXE) represent a particular type of alternative splicing, in which one and only one exon from an array is included in the mature RNA. A number of genes with MXE do so by using a mechanism that depends on RNA structure. Transcripts of these genes contain multiple sites called selector sequences that are all complementary to a regulatory element called the docking site; only one of the competing base pairings can form at a time, which exposes one exon from the cluster to the spliceosome. MXE tend to have similar lengths and sequence content and are believed to originate through tandem genomic duplications. Here, we report that pre-mRNAs of this class of exons have an increased capacity to fold into competing secondary structures. We propose an evolutionary mechanism for the generation of such structures via duplications that affect not only exons, but also their adjacent introns with stem-loop structures. If one of the two arms of a stem-loop is duplicated, it will generate two selector sequences that compete for the same docking site, a pattern that is associated with MXE splicing. A similar partial duplication of two independent stem-loops produces a pattern that is consistent with the so-called bidirectional pairing model. These models explain why tandem exon duplications frequently result in mutually exclusive splicing.