RNA secondary structure in mutually exclusive splicing

RNA structure in alternative splicing regulation: from mechanism to therapy

Alternative splicing is a highly intricate process that plays a crucial role in post-transcriptional regulation and significantly expands the functional proteome of a limited number of coding genes in eukaryotes. Its regulation is multifactorial, with RNA structure exerting a significant impact. Aberrant RNA conformations lead to dysregulation of splicing patterns, which directly affects the manifestation of disease symptoms. In this review, the molecular mechanisms of RNA secondary structure-mediated splicing regulation are summarized, with a focus on the complex interplay between aberrant RNA conformations and disease phenotypes resulted from splicing defects. This study also explores additional factors that reshape structural conformations, enriching our understanding of the mechanistic network underlying structure-mediated splicing regulation. In addition, an emphasis has been placed on the clinical role of targeting aberrant splicing corrections in human diseases. The principal mechanisms of action behind this phenomenon are described, followed by a discussion of prospective development strategies and pertinent challenges.

Conserved intronic secondary structures with concealed branch sites regulate alternative splicing of poison exons

Alternative splicing (AS) generates multiple RNA isoforms and increases the complexities of transcriptomes and proteomes. However, it remains unclear how RNA structures contribute to AS regulation. Here, we systematically search transcriptomes for secondary structures with concealed branch sites (BSs) in the alternatively spliced introns and predict thousands of them from six organisms, of which many are evolutionarily conserved. Intriguingly, a highly conserved stem-loop structure with concealed BSs is found in animal SF3B3 genes and colocalizes with a downstream poison exon (PE). Destabilization of this structure allows increased usage of the BSs and results in enhanced PE inclusion in human and Drosophila cells, leading to decreased expression of SF3B3. This structure is experimentally validated using an in-cell SHAPE-MaP assay. Through RNA interference screens of 28 RNA-binding proteins, we find that this stem-loop structure is sensitive to U2 factors. Furthermore, we find that SF3B3 also facilitates DNA repair and protects genome stability by enhancing interaction between ERCC6/CSB and arrested RNA polymerase II. Importantly, both Drosophila and human cells with the secondary structure mutated by genome editing exhibit altered DNA repair in vivo. This study provides a novel and common mechanism for AS regulation of PEs and reveals a physiological function of SF3B3 in DNA repair.

RNA structure: implications in viral infections and neurodegenerative diseases

RNA is an intermediary between DNA and protein, a catalyzer of biochemical reactions, and a regulator of genes and transcripts. RNA structures are essential for complicated functions. Recent years have witnessed rapid advancements in RNA secondary structure probing techniques. These technological strides provided comprehensive insights into RNA structures, which significantly contributed to our understanding of diverse cellular regulatory processes, including gene regulation, epigenetic regulation, and post-transactional regulation. Meanwhile, they have facilitated the creation of therapeutic tools for tackling human diseases. Despite their therapeutic applications, RNA structure probing methods also offer a promising avenue for exploring the mechanisms of human diseases, potentially providing the key to overcoming existing research constraints and obtaining the in-depth information necessary for a deeper understanding of disease mechanisms.

Confluence and convergence of Dscam and Pcdh cell-recognition codes

The ability of neurites of the same neuron to avoid each other (self-avoidance) is a conserved feature in both invertebrates and vertebrates. The key to self-avoidance is the generation of a unique subset of cell-surface proteins in individual neurons engaging in isoform-specific homophilic interactions that drive neurite repulsion rather than adhesion. Among these cell-surface proteins are fly Dscam1 and vertebrate clustered protocadherins (cPcdhs), as well as the recently characterized shortened Dscam (sDscam) in the Chelicerata. Herein, we review recent advances in our understanding of how cPcdh, Dscam, and sDscam cell-surface recognition codes are expressed and translated into cellular functions essential for neural wiring.

Cis mutagenesis in vivo reveals extensive noncanonical functions of Dscam1 isoforms in neuronal wiring

Drosophila Down syndrome cell adhesion molecule 1 (Dscam1) encodes tens of thousands of cell recognition molecules via alternative splicing, which are required for neural function. A canonical self-avoidance model seems to provide a central mechanistic basis for Dscam1 functions in neuronal wiring. Here, we reveal extensive noncanonical functions of Dscam1 isoforms in neuronal wiring. We generated a series of allelic cis mutations in Dscam1, encoding a normal number of isoforms, but with an altered isoform composition. Despite normal dendritic self-avoidance and self-/nonself-discrimination in dendritic arborization (da) neurons, which is consistent with the canonical self-avoidance model, these mutants exhibited strikingly distinct spectra of phenotypic defects in the three types of neurons: up to ∼60% defects in mushroom bodies, a significant increase in branching and growth in da neurons, and mild axonal branching defects in mechanosensory neurons. Remarkably, the altered isoform composition resulted in increased dendrite growth yet inhibited axon growth. Moreover, reducing Dscam1 dosage exacerbated axonal defects in mushroom bodies and mechanosensory neurons but reverted dendritic branching and growth defects in da neurons. This splicing-tuned regulation strategy suggests that axon and dendrite growth in diverse neurons cell-autonomously require Dscam1 isoform composition. These findings provide important insights into the functions of Dscam1 isoforms in neuronal wiring.

Thermomorphogenesis: Opportunities and challenges in posttranscriptional regulation

Plants exposed to mildly elevated temperatures display morphological and developmental changes collectively termed thermomorphogenesis. This adaptative process has several undesirable consequences to food production, including yield reduction and increased vulnerability to pathogens. Understanding thermomorphogenesis is, thus, critical for understanding how plants will respond to increasingly warmer temperature conditions, such as those caused by climate change. Recently, we have made major advances in that direction, and it has become apparent that plants resource to a broad range of molecules and molecular mechanisms to perceive and respond to increases in environmental temperature. However, most of our efforts have been focused on regulation of transcription and protein abundance and activity, with an important gap encompassing nearly all processes involving RNA (i.e., posttranscriptional regulation). Here, I summarized our current knowledge of thermomorphogenesis involving transcriptional, posttranscriptional, and posttranslational regulation, focused on opportunities and challenges in understanding posttranscriptional regulation-a fertile field for exciting new discoveries.

Trans-splicing facilitated by RNA pairing greatly expands sDscam isoform diversity but not homophilic binding specificity

The Down syndrome cell adhesion molecule 1 (Dscam1) gene can generate tens of thousands of isoforms via alternative splicing, which is essential for nervous and immune functions. Chelicerates generate approximately 50 to 100 shortened Dscam (sDscam) isoforms by alternative promoters, similar to mammalian protocadherins. Here, we reveal that trans-splicing markedly increases the repository of sDscamβ isoforms in Tetranychus urticae. Unexpectedly, every variable exon cassette engages in trans-splicing with constant exons from another cluster. Moreover, we provide evidence that competing RNA pairing not only governs alternative cis-splicing but also facilitates trans-splicing. Trans-spliced sDscam isoforms mediate cell adhesion ability but exhibit the same homophilic binding specificity as their cis-spliced counterparts. Thus, we reveal a single sDscam locus that generates diverse adhesion molecules through cis- and trans-splicing coupled with alternative promoters. These findings expand understanding of the mechanism underlying molecular diversity and have implications for the molecular control of neuronal and/or immune specificity.

Recent advances in RNA structurome

RNA structures are essential to support RNA functions and regulation in various biological processes. Recently, a range of novel technologies have been developed to decode genome-wide RNA structures and novel modes of functionality across a wide range of species. In this review, we summarize key strategies for probing the RNA structurome and discuss the pros and cons of representative technologies. In particular, these new technologies have been applied to dissect the structural landscape of the SARS-CoV-2 RNA genome. We also summarize the functionalities of RNA structures discovered in different regulatory layers-including RNA processing, transport, localization, and mRNA translation-across viruses, bacteria, animals, and plants. We review many versatile RNA structural elements in the context of different physiological and pathological processes (e.g., cell differentiation, stress response, and viral replication). Finally, we discuss future prospects for RNA structural studies to map the RNA structurome at higher resolution and at the single-molecule and single-cell level, and to decipher novel modes of RNA structures and functions for innovative applications.

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.

Hidden RNA pairings counteract the "first-come, first-served" splicing principle to drive stochastic choice in Dscam1 splice variants

Drosophila melanogaster Dscam1 encodes 38,016 isoforms via mutually exclusive splicing; however, the regulatory mechanism behind this is not fully understood. Here, we found a set of hidden RNA secondary structures that balance the stochastic choice of Dscam1 splice variants (designated balancer RNA secondary structures). In vivo mutational analyses revealed the dual function of these balancer interactions in driving the stochastic choice of splice variants, through enhancement of the inclusion of distal exon 6s by cooperating with docking site–selector pairing to form a stronger multidomain pre-mRNA structure and through simultaneous repression of the inclusion of proximal exon 6s by antagonizing their docking site–selector pairings. Thus, we provide an elegant molecular model based on competition and cooperation between two sets of docking site–selector and balancer pairings, which counteracts the “first-come, first-served” principle. Our findings provide conceptual and mechanistic insight into the dynamics and functions of long-range RNA secondary structures.

May the Odds Be Ever in Your Favor: Non-deterministic Mechanisms Diversifying Cell Surface Molecule Expression

How does the information in the genome program the functions of the wide variety of cells in the body? While the development of biological organisms appears to follow an explicit set of genomic instructions to generate the same outcome each time, many biological mechanisms harness molecular noise to produce variable outcomes. Non-deterministic variation is frequently observed in the diversification of cell surface molecules that give cells their functional properties, and is observed across eukaryotic clades, from single-celled protozoans to mammals. This is particularly evident in immune systems, where random recombination produces millions of antibodies from only a few genes; in nervous systems, where stochastic mechanisms vary the sensory receptors and synaptic matching molecules produced by different neurons; and in microbial antigenic variation. These systems employ overlapping molecular strategies including allelic exclusion, gene silencing by constitutive heterochromatin, targeted double-strand breaks, and competition for limiting enhancers. Here, we describe and compare five stochastic molecular mechanisms that produce variety in pathogen coat proteins and in the cell surface receptors of animal immune and neuronal cells, with an emphasis on the utility of non-deterministic variation.

Dscam in arthropod immune priming: What is known and what remains unknown

A popular view in the current academic circle is that invertebrates have no adaptive immunity. However, the immune memory and specificity of invertebrates pose a serious challenge to this view and constitute immune priming based on innate immunity. The Down syndrome cell adhesion molecule (Dscam) gene of invertebrates, with approximately 10,000 alternatively spliced isoforms, has a unique characteristic: it specifically binds to different types of bacteria and promotes cell phagocytosis; owing to its antibody-like function, Dscam is a key candidate protein for immune priming. However, the high molecular diversity of Dscam and the gaps and inconsistencies in the existing research make the study of regulation of immune priming by Dscam challenging. In recent years, significant research has been conducted on the Dscam-regulated immune functions in insects and crustaceans, providing preliminary results for Dscam-regulated innate immunity and immune priming, but some important questions remain unresolved. In this review, we summarize the existing knowledge about Dscam-regulated immunity and discuss three yet unanswered questions, the study of which may improve the understanding of the role of Dscam-regulated immune priming in invertebrates.

Improved recognition of splice sites in A. thaliana by incorporating secondary structure information into sequence-derived features: a computational study

Identification of splice sites is an important aspect with regard to the prediction of gene structure. In most of the existing splice site prediction studies, machine learning algorithms coupled with sequence-derived features have been successfully employed for splice site recognition. However, the splice site identification by incorporating the secondary structure information is lacking, particularly in plant species. Thus, we made an attempt in this study to evaluate the performance of structural features on the splice site prediction accuracy in Arabidopsis thaliana. Prediction accuracies were evaluated with the sequence-derived features alone as well as by incorporating the structural features into the sequence-derived features, where support vector machine (SVM) was employed as prediction algorithm. Both short (40 base pairs) and long (105 base pairs) sequence datasets were considered for evaluation. After incorporating the secondary structure features, improvements in accuracies were observed only for the longer sequence dataset and the improvement was found to be higher with the sequence-derived features that accounted nucleotide dependencies. On the other hand, either a little or no improvement in accuracies was found for the short sequence dataset. The performance of SVM was further compared with that of LogitBoost, Random Forest (RF), AdaBoost and XGBoost machine learning methods. The prediction accuracies of SVM, AdaBoost and XGBoost were observed to be at par and higher than that of RF and LogitBoost algorithms. While prediction was performed by taking all the sequence-derived features along with the structural features, a little improvement in accuracies was found as compared to the combination of individual sequence-based features and structural features. To the best of our knowledge, this is the first attempt concerning the computational prediction of splice sites using machine learning methods by incorporating the secondary structure information into the sequence-derived features. All the source codes are available at https://github.com/meher861982/SSFeature.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-03036-8.

Reverse complementary matches simultaneously promote both back-splicing and exon-skipping

Background

Circular RNAs (circRNAs) play diverse roles in different biological and physiological environments and are always expressed in a tissue-specific manner. Especially, circRNAs are enriched in the brain tissues of almost all investigated species, including humans, mice, Drosophila, etc. Although circRNAs were found in C. elegans, the neuron-specific circRNA data is not available yet. Exon-skipping is found to be correlated to circRNA formation, but the mechanisms that link them together are not clear.

Results

Here, through large-scale neuron isolation from the first larval (L1) stage of C. elegans followed by RNA sequencing with ribosomal RNA depletion, the neuronal circRNA data in C. elegans were obtained. Hundreds of novel circRNAs were annotated with high accuracy. circRNAs were highly expressed in the neurons of C. elegans and were positively correlated to the levels of their cognate linear mRNAs. Disruption of reverse complementary match (RCM) sequences in circRNA flanking introns effectively abolished circRNA formation. In the zip-2 gene, deletion of either upstream or downstream RCMs almost eliminated the production of both the circular and the skipped transcript. Interestingly, the 13-nt RCM in zip-2 is highly conserved across five nematode ortholog genes, which show conserved exon-skipping patterns. Finally, through in vivo one-by-one mutagenesis of all the splicing sites and branch points required for exon-skipping and back-splicing in the zip-2 gene, I showed that back-splicing still happened without exon-skipping, and vice versa.

Conclusions

Through protocol optimization, total RNA obtained from sorted neurons is increased to hundreds of nanograms. circRNAs highly expressed in the neurons of C. elegans are more likely to be derived from genes also highly expressed in the neurons. RCMs are abundant in circRNA flanking introns, and RCM-deletion is an efficient way to knockout circRNAs. More importantly, these RCMs are not only required for back-splicing but also promote the skipping of exon(s) to be circularized. Finally, RCMs in circRNA flanking introns can directly promote both exon-skipping and back-splicing, providing a new explanation for the correlation between them.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07910-w.

Intron-targeted mutagenesis reveals roles for Dscam1 RNA pairing architecture-driven splicing bias in neuronal wiring

Drosophila melanogaster Down syndrome cell adhesion molecule (Dscam1) can generate 38,016 different isoforms through largely stochastic, yet highly biased, alternative splicing. These isoforms are required for nervous functions. However, the functional significance of splicing bias remains unknown. Here, we provide evidence that Dscam1 splicing bias is required for mushroom body (MB) axonal wiring. We generate mutant flies with normal overall protein levels and an identical number but global changes in exon 4 and 9 isoform bias (DscamΔ4D^-/- and DscamΔ9D^-/-), respectively. In contrast to DscamΔ4D^-/-, DscamΔ9D^-/- exhibits remarkable MB defects, suggesting a variable domain-specific requirement for isoform bias. Importantly, changes in isoform bias cause axonal defects but do not influence the self-avoidance of axonal branches. We conclude that, in contrast to the isoform number that provides the molecular basis for neurite self-avoidance, isoform bias may play a role in MB axonal wiring by influencing non-repulsive signaling.

Unannotated splicing regulatory elements in deep intron space

Deep intron space harbors a diverse array of splicing regulatory elements that cooperate with better-known exon-proximal elements to enforce proper tissue-specific and development-specific pre-mRNA processing. Many deep intron elements have been highly conserved through vertebrate evolution, yet remain poorly annotated in the human genome. Recursive splicing exons (RS-exons) and intraexons promote noncanonical, multistep resplicing pathways in long introns, involving transient intermediate structures that are greatly underrepresented in RNA-seq datasets. Decoy splice sites and decoy exons act at a distance to inhibit splicing catalysis at annotated splice sites, with functional consequences such as exon skipping and intron retention. RNA:RNA bridges can juxtapose distant sequences within or across introns to activate deep intron splicing enhancers and silencers, to loop out exons to be skipped, or to select one member of a mutually exclusive set of exons. Similarly, protein bridges mediated by interactions among transcript-bound RNA binding proteins (RBPs) can modulate splicing outcomes. Experimental disruption of deep intron elements serving any of these functions can abrogate normal splicing, strongly suggesting that natural mutations of deep intron elements can do likewise to cause human disease. Understanding noncanonical splicing pathways and discovering deep intron regulatory signals, many of which map hundreds to many thousands of nucleotides from annotated splice junctions, is of great academic interest for basic scientists studying alternative splicing mechanisms. Hopefully, this knowledge coupled with increased analysis of deep intron sequences will also have important medical applications, as better interpretation of deep intron mutations may reveal new disease mechanisms and suggest new therapies. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing.

Complex RNA Secondary Structures Mediate Mutually Exclusive Splicing of Coleoptera Dscam1

Mutually exclusive splicing is an important mechanism for expanding protein diversity. An extreme example is the Down syndrome cell adhesion molecular (Dscam1) gene of insects, containing four clusters of variable exons (exons 4, 6, 9, and 17), which potentially generates tens of thousands of protein isoforms through mutually exclusive splicing, of which regulatory mechanisms are still elusive. Here, we systematically analyzed the variable exon 4, 6, and 9 clusters of Dscam1 in Coleoptera species. Through comparative genomics and RNA secondary structure prediction, we found apparent evidence that the evolutionarily conserved RNA base pairing mediates mutually exclusive splicing in the Dscam1 exon 4 cluster. In contrast to the fly exon 6, most exon 6 selector sequences in Coleoptera species are partially located in the variable exon region. Besides, bidirectional RNA–RNA interactions are predicted to regulate the mutually exclusive splicing of variable exon 9 of Dscam1. Although the docking sites in exon 4 and 9 clusters are clade specific, the docking sites-selector base pairing is conserved in secondary structure level. In short, our result provided a mechanistic framework for the application of long-range RNA base pairings in regulating the mutually exclusive splicing of Coleoptera Dscam1.

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.

RNA secondary structures in Dscam1 mutually exclusive splicing: unique evolutionary signature from the midge

The Drosophila melanogaster gene Dscam1 potentially generates 38,016 distinct isoforms via mutually exclusive splicing, which are required for both nervous and immune functions. However, the mechanism underlying splicing regulation remains obscure. Here we show apparent evolutionary signatures characteristic of competing RNA secondary structures in exon clusters 6 and 9 of Dscam1 in the two midge species (Belgica antarctica and Clunio marinus). Surprisingly, midge Dscam1 encodes only ∼6000 different isoforms through mutually exclusive splicing. Strikingly, the docking site of the exon 6 cluster is conserved in almost all insects and crustaceans but is specific in the midge; however, the docking site-selector base-pairings are conserved. Moreover, the docking site is complementary to all predicted selector sequences downstream from every variable exon 9 of the midge Dscam1, which is in accordance with the broad spectrum of their isoform expression. This suggests that these cis-elements mainly function through the formation of long-range base-pairings. This study provides a vital insight into the evolution and mechanism of Dscam1 alternative splicing.

Dscam in immunity: A question of diversity in insects and crustaceans

In insects and crustaceans, thousands of Down syndrome cell adhesion molecules (Dscam) can be generated by alternative splicing of variable exons from a single-locus gene, Dscam-hv. This extraordinarily versatile gene (38,016 protein isoforms produced in Drosophila) was first proposed to be involved in exon guidance and subsequently implicated in immunity as a hypervariable immune molecule. Almost 20 y after discovery of Dscam-hv, there have been many studies in insects and crustaceans regarding roles of Dscam in immunity, with many similarities and concurrently, many differences. Here, we review the current status of Dscam-hv, presented as a comparison of similarities and differences in insects and crustaceans and discuss hypotheses of Dscam functions in immunity.

Alternate RNA Structures

RNA molecules fold into complex three-dimensional structures that sample alternate conformations ranging from minor differences in tertiary structure dynamics to major differences in secondary structure. This allows them to form entirely different substructures with each population potentially giving rise to a distinct biological outcome. The substructures can be partitioned along an existing energy landscape given a particular static cellular cue or can be shifted in response to dynamic cues such as ligand binding. We review a few key examples of RNA molecules that sample alternate conformations and how these are capitalized on for control of critical regulatory functions.

Acute thiamethoxam toxicity in honeybees is not enhanced by common fungicide and herbicide and lacks stress-induced changes in mRNA splicing

Securing food supply for a growing population is a major challenge and heavily relies on the use of agrochemicals to maximize crop yield. It is increasingly recognized, that some neonicotinoid insecticides have a negative impact on non-target organisms, including important pollinators such as the European honeybee Apis mellifera. Toxicity of neonicotinoids may be enhanced through simultaneous exposure with additional pesticides, which could help explain, in part, the global decline of honeybee colonies. Here we examined whether exposure effects of the neonicotinoid thiamethoxam on bee viability are enhanced by the commonly used fungicide carbendazim and the herbicide glyphosate. We also analysed alternative splicing changes upon pesticide exposure in the honeybee. In particular, we examined transcripts of three genes: (i) the stress sensor gene X box binding protein-1 (Xbp1), (ii) the Down Syndrome Cell Adhesion Molecule (Dscam) gene and iii) the embryonic lethal/abnormal visual system (elav) gene, which are important for neuronal function. Our results showed that acute thiamethoxam exposure is not enhanced by carbendazim, nor glyphosate. Toxicity of the compounds did not trigger stress-induced, alternative splicing in the analysed mRNAs, thereby leaving dormant a cellular response pathway to these man-made environmental perturbations.

Srrm234, but not canonical SR and hnRNP proteins, drive inclusion of Dscam exon 9 variable exons

Alternative splicing of pre-mRNA is a major mechanism to diversify protein functionality in metazoans from a limited number of genes. The Drosophila melanogaster Down syndrome cell adhesion molecule (Dscam) gene, which is important for neuronal wiring and phagocytosis of bacteria, can generate up to 38,016 isoforms by mutually exclusive alternative splicing in four clusters of variable exons. However, it is not understood how a specific exon is chosen from the many variables and how variable exons are prevented from being spliced together. A main role in the regulation of Dscam alternative splicing has been attributed to RNA binding proteins (RBPs), but how they impact on exon selection is not well understood. Serine-arginine rich (SR) proteins and hnRNP proteins are the two main types of RBPs with major roles in exon definition and splice site selection. Here, we analyzed the role of SR and hnRNP proteins in Dscam exon 9 alternative splicing in mutant Drosophila melanogaster embryos because of their essential function for development. Strikingly, loss or overexpression of canonical SR and hnRNP proteins even when multiple proteins are depleted together, does not affect Dscam alternative exon selection very dramatically. Conversely, noncanonical SR protein Serine-arginine repetitive matrix 2/3/4 (Srrm234) is a main determinant of exon inclusion in the Dscam exon 9 cluster. Since long-range base-pairings are absent in the exon 9 cluster, our data argue for a small complement of regulatory factors as main determinants of exon inclusion in the Dscam exon 9 cluster.

Plasmid-based gap-repair recombineered transgenes reveal a central role for introns in mutually exclusive alternative splicing in Down Syndrome Cell Adhesion Molecule exon 4

Alternative splicing is a key feature of human genes, yet studying its regulation is often complicated by large introns. The Down Syndrome Cell Adhesion Molecule (Dscam) gene from Drosophila is one of the most complex genes generating vast molecular diversity by mutually exclusive alternative splicing. To resolve how alternative splicing in Dscam is regulated, we first developed plasmid-based UAS reporter genes for the Dscam variable exon 4 cluster and show that its alternative splicing is recapitulated by GAL4-mediated expression in neurons. We then developed gap-repair recombineering to very efficiently manipulate these large reporter plasmids in Escherichia coli using restriction enzymes or sgRNA/Cas9 DNA scission to capitalize on the many benefits of plasmids in phiC31 integrase-mediated transgenesis. Using these novel tools, we show that inclusion of Dscam exon 4 variables differs little in development and individual flies, and is robustly determined by sequences harbored in variable exons. We further show that introns drive selection of both proximal and distal variable exons. Since exon 4 cluster introns lack conserved sequences that could mediate robust long-range base-pairing to bring exons into proximity for splicing, our data argue for a central role of introns in mutually exclusive alternative splicing of Dscam exon 4 cluster.

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.

C3a and suPAR drive versican V1 expression in tubular cells of focal segmental glomerulosclerosis

Chronic tubulointerstitial injury impacts the prognosis of focal segmental glomerulosclerosis (FSGS). We found that the level of versican V1 was increased in tubular cells of FSGS patients. Tubular cell-derived versican V1 induced proliferation and collagen synthesis by activating the CD44/Smad3 pathway in fibroblasts. Both urine C3a and suPAR were increased and bound to the tubular cells in FSGS patients. C3a promoted the transcription of versican by activating the AKT/β-catenin pathway. C3aR knockout decreased the expression of versican in Adriamycin-treated (ADR-treated) mice. On the other hand, suPAR bound to integrin β6 and activated Rac1, which bound to SRp40 at the 5' end of exon 7 in versican pre-mRNA. This binding inhibited the 3'-end splicing of intron 6 and the base-pair interactions between intron 6 and intron 8, leading to the formation of versican V1. Cotreatment with ADR and suPAR specifically increased the level of versican V1 in tubulointerstitial tissues and caused more obvious interstitial fibrosis in mice than treatment with only ADR. Altogether, our results show that C3a and suPAR drive versican V1 expression in tubular cells by promoting transcription and splicing, respectively, and the increases in tubular cell-derived versican V1 induce interstitial fibrosis by activating fibroblasts in FSGS.

B52 promotes alternative splicing of Dscam in Chinese mitten crab, Eriocheir sinensis

B52 is a member of the classical serine/arginine (SR)-rich proteins, which are phylogenetically conserved and play significant roles in mRNA maturation, including alternative splicing. In the present study, the docking site, selector sequences and locus control region of the Chinese mitten crab (Eriocheir sinensis) Down syndrome cell adhesion molecule (EsDscam) were identified. Alternative splicing of Dscam is essential to generate different isoforms. We also isolated and characterised the B52 gene from E. sinensis (EsB52). The 876 bp open reading frame of EsB52 encodes a 291 amino acid residue polypeptide, and EsB52 has two RNA recognition motifs (RRMs) at the N-terminus and an arginine/serine-rich domain at the C-terminus. Each RRM contains two degenerate short submotifs, RNP-1 and RNP2. Analysis of tissue distribution revealed that EsB52 mRNA expression was widespread in all tested tissues, and especially high in brain and hemocytes. In hemocytes, EsB52 was upregulated significantly after stimulation with pathogen-associated molecular patterns and bacteria. Furthermore, EsB52 RNAi decreased the number of Ig7 inclusion in mRNA rather than Ig2 or Ig3. Taken together, these findings suggest that EsB52 acts as an alternative splicing activator of EsDscam.

Revisiting Dscam diversity: lessons from clustered protocadherins

The complexity of neuronal wiring relies on the extraordinary recognition diversity of cell surface molecules. Drosophila Dscam1 and vertebrate clustered protocadherins (Pcdhs) are two classic examples of the striking diversity from a complex genomic locus, wherein the former encodes more than 10,000 distinct isoforms via alternative splicing, while the latter employs alternative promoters to attain isoform diversity. These structurally unrelated families show remarkably striking molecular parallels and even similar functions. Recent studies revealed a novel Dscam gene family with tandemly arrayed 5' cassettes in Chelicerata (e.g., the scorpion Mesobuthus martensii and the tick Ixodes scapularis), similar to vertebrate clustered Pcdhs. Likewise, octopus shows a more remarkable expansion of the Pcdh isoform repertoire than human. These discoveries of Dscam and Pcdh diversification reshape the evolutionary landscape of recognition molecule diversity and provide a greater understanding of convergent molecular strategies for isoform diversity. This article reviews new insights into the evolution, regulatory mechanisms, and functions of Dscam and Pcdh isoform diversity. In particular, the convergence of clustered Dscams and Pcdhs is highlighted.

Competing RNA pairings in complex alternative splicing of a 3' variable region

Alternative pre-mRNA splicing remarkably expands protein diversity in eukaryotes. Drosophila PGRP-LC can generate three major 3' splice isoforms that exhibit distinct innate immune recognition and defenses against various microbial infections. However, the regulatory mechanisms underlying the uniquely biased splicing pattern at the 3' variable region remain unclear. Here we show that competing RNA pairings control the unique splicing of the 3' variable region of Drosophila PGRP-LC pre-mRNA. We reveal three roles by which these RNA pairings jointly regulate the 3' variant selection through activating the proximal 3' splice site and concurrently masking the intron-proximal 5' splice site, in combination with physical competition of RNA pairing. We also reveal that competing RNA pairings regulate alternative splicing of the highly complex 3' variable regions of Drosophila CG42235 and Pip Our findings will facilitate a better understanding of the regulatory mechanisms of highly complex alternative splicing as well as highly variable 3' processing.

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.

Towards Long-Range RNA Structure Prediction in Eukaryotic Genes

The ability to form an intramolecular structure plays a fundamental role in eukaryotic RNA biogenesis. Proximate regions in the primary transcripts fold into a local secondary structure, which is then hierarchically assembled into a tertiary structure that is stabilized by RNA-binding proteins and long-range intramolecular base pairings. While the local RNA structure can be predicted reasonably well for short sequences, long-range structure at the scale of eukaryotic genes remains problematic from the computational standpoint. The aim of this review is to list functional examples of long-range RNA structures, to summarize current comparative methods of structure prediction, and to highlight their advances and limitations in the context of long-range RNA structures. Most comparative methods implement the “first-align-then-fold” principle, i.e., they operate on multiple sequence alignments, while functional RNA structures often reside in non-conserved parts of the primary transcripts. The opposite “first-fold-then-align” approach is currently explored to a much lesser extent. Developing novel methods in both directions will improve the performance of comparative RNA structure analysis and help discover novel long-range structures, their higher-order organization, and RNA⁻RNA interactions across the transcriptome.

Pathogen-Specific Binding Soluble Down Syndrome Cell Adhesion Molecule (Dscam) Regulates Phagocytosis via Membrane-Bound Dscam in Crab

The Down syndrome cell adhesion molecule (Dscam) gene is an extraordinary example of diversity that can produce thousands of isoforms and has so far been found only in insects and crustaceans. Cumulative evidence indicates that Dscam may contribute to the mechanistic foundations of specific immune responses in insects. However, the mechanism and functions of Dscam in relation to pathogens and immunity remain largely unknown. In this study, we identified the genome organization and alternative Dscam exons from Chinese mitten crab, Eriocheir sinensis. These variants, designated EsDscam, potentially produce 30,600 isoforms due to three alternatively spliced immunoglobulin (Ig) domains and a transmembrane domain. EsDscam was significantly upregulated after bacterial challenge at both mRNA and protein levels. Moreover, bacterial specific EsDscam isoforms were found to bind specifically with the original bacteria to facilitate efficient clearance. Furthermore, bacteria-specific binding of soluble EsDscam via the complete Ig1–Ig4 domain significantly enhanced elimination of the original bacteria via phagocytosis by hemocytes; this function was abolished by partial Ig1–Ig4 domain truncation. Further studies showed that knockdown of membrane-bound EsDscam inhibited the ability of EsDscam with the same extracellular region to promote bacterial phagocytosis. Immunocytochemistry indicated colocalization of the soluble and membrane-bound forms of EsDscam at the hemocyte surface. Far-Western and coimmunoprecipitation assays demonstrated homotypic interactions between EsDscam isoforms. This study provides insights into a mechanism by which soluble Dscam regulates hemocyte phagocytosis via bacteria-specific binding and specific interactions with membrane-bound Dscam as a phagocytic receptor.

RNA secondary structure profiling in zebrafish reveals unique regulatory features

Background

RNA is known to play diverse roles in gene regulation. The clues for this regulatory function of RNA are embedded in its ability to fold into intricate secondary and tertiary structure.

Results

We report the transcriptome-wide RNA secondary structure in zebrafish at single nucleotide resolution using Parallel Analysis of RNA Structure (PARS). This study provides the secondary structure map of zebrafish coding and non-coding RNAs. The single nucleotide pairing probabilities of 54,083 distinct transcripts in the zebrafish genome were documented. We identified RNA secondary structural features embedded in functional units of zebrafish mRNAs. Translation start and stop sites were demarcated by weak structural signals. The coding regions were characterized by the three-nucleotide periodicity of secondary structure and display a codon base specific structural constrain. The splice sites of transcripts were also delineated by distinct signature signals. Relatively higher structural signals were observed at 3’ Untranslated Regions (UTRs) compared to Coding DNA Sequence (CDS) and 5’ UTRs. The 3′ ends of transcripts were also marked by unique structure signals. Secondary structural signals in long non-coding RNAs were also explored to better understand their molecular function.

Conclusions

Our study presents the first PARS-enabled transcriptome-wide secondary structure map of zebrafish, which documents pairing probability of RNA at single nucleotide precision. Our findings open avenues for exploring structural features in zebrafish RNAs and their influence on gene expression.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4497-0) contains supplementary material, which is available to authorized users.

Mutually exclusive alternative splicing of pre-mRNAs

Pre-mRNA alternative splicing is an important mechanism used to expand protein diversity in higher eukaryotes, and mutually exclusive splicing is a specific type of alternative splicing in which only one of the exons in a cluster is included in functional transcripts. The most extraordinary example of this is the Drosophila melanogaster Down's syndrome cell adhesion molecule gene (Dscam), which potentially encodes 38,016 different isoforms through mutually exclusive splicing. Mutually exclusive splicing is a unique and challenging model that can be used to elucidate the evolution, regulatory mechanism, and function of alternative splicing. The use of new approaches has not only greatly expanded the mutually exclusive exome, but has also enabled the systematic analyses of single-cell alternative splicing during development. Furthermore, the identification of long-range RNA secondary structures provides a mechanistic framework for the regulation of mutually exclusive splicing (i.e., Dscam splicing). This article reviews recent insights into the identification, underlying mechanism, and roles of mutually exclusive splicing. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Intron specificity in pre-mRNA splicing

The occurrence of spliceosomal introns in eukaryotic genomes is highly diverse and ranges from few introns in an organism to multiple introns per gene. Introns vary with respect to their lengths, strengths of splicing signals, and position in resident genes. Higher intronic density and diversity in genetically complex organisms relies on increased efficiency and accuracy of spliceosomes for pre-mRNA splicing. Since intron diversity is critical for functions in RNA stability, regulation of gene expression and alternative splicing, RNA-binding proteins, spliceosomal regulatory factors and post-translational modifications of splicing factors ought to make the splicing process intron-specific. We recently reported function and regulation of a ubiquitin fold harboring splicing regulator, Sde2, which following activation by ubiquitin-specific proteases facilitates excision of selected introns from a subset of multi-intronic genes in Schizosaccharomyces pombe (Thakran et al. EMBO J, https://doi.org/10.15252/embj.201796751 , 2017). By reviewing our findings with understandings of intron functions and regulated splicing processes, we propose possible functions and mechanism of intron-specific pre-mRNA splicing and suggest that this process is crucial to highlight importance of introns in eukaryotic genomes.

Characterize a typically Dscam with alternative splicing in mud crab Scylla paramamosain

As a member of the immunoglobulin superfamily, Down syndrome cell adhesion molecule (Dscam) could function in the innate immunity of invertebrates. Recently, it is shown that arthropod Dscams play similar functions as antibodies in the adaptive immune system. Dscam could produce thousands of isoforms by alternative splicing and specifically bind to various pathogens. In the present study, we cloned the first Dscam from mud crab Scylla paramamosain (SpDscam), with full-length cDNA 7363 bp containing an open reading frame (ORF) of 6069bp and encoding 2022 amino acids, which had typical domain architecture as other arthropods, i.e., 10 immunoglobulin domains (Ig), 6 fibronectin type 3 domains (FN III), transmembrane and cytoplasmic tail. Quantitative real-time PCR revealed that SpDscam was highly expressed in brain, skin, muscle, intestine and hepatopancreas, but weakly expressed in hemolymph, heart and gill. SpDscam had three alternative splicing regions, located at the N-terminal of Ig2 and Ig3 as well as on the whole Ig7. In these regions, 32, 41 and 14 exons were detected, together with the two exon types of transmembrane domain, indicating SpDscam could potentially encode at least 36,736 unique isoforms. SpDscam induced by Vibrio parahaemolyticus challenge had strong binding ability to V. parahaemolyticus. Further, SpDscam induced by V. parahaemolyticus possessed a clearance of V. parahaemolyticus in S. paramamosain. Collectively, the results indicated SpDscam was a hypervariable pattern-recognition receptor (PRR) by alternative splicing in innate immunity system of mud crab S. paramamosain.

The determinants of alternative RNA splicing in human cells

Alternative splicing represents an important level of the regulation of gene function in eukaryotic organisms. It plays a critical role in virtually every biological process within an organism, including regulation of cell division and cell death, differentiation of tissues in the embryo and the adult organism, as well as in cellular response to diverse environmental factors. In turn, studies of the last decade have shown that alternative splicing itself is controlled by different mechanisms. Unfortunately, there is no clear understanding of how these diverse mechanisms, or determinants, regulate and constrain the set of alternative RNA species produced from any particular gene in every cell of the human body. Here, we provide a consolidated overview of alternative splicing determinants including RNA-protein interactions, epigenetic regulation via chromatin remodeling, coupling of transcription-to-alternative splicing, effect of secondary structures in pre-RNA, and function of the RNA quality control systems. We also extensively and critically discuss some mechanistic insights on coordinated inclusion/exclusion of exons during the formation of mature RNA molecules. We conclude that the final structure of RNA is pre-determined by a complex interplay between cis- and trans-acting factors. Altogether, currently available empirical data significantly expand our understanding of the functioning of the alternative splicing machinery of cells in normal and pathological conditions. On the other hand, there are still many blind spots that require further deep investigations.

Role and convergent evolution of competing RNA secondary structures in mutually exclusive splicing

Exon or cassette duplication is an important means of expanding protein and functional diversity through mutually exclusive splicing. However, the mechanistic basis of this process in non-arthropod species remains poorly understood. Here, we demonstrate that MRP1 genes underwent tandem exon duplication in Nematoda, Platyhelminthes, Annelida, Mollusca, Arthropoda, Echinodermata, and early-diverging Chordata but not in late-diverging vertebrates. Interestingly, these events were of independent origin in different phyla, suggesting convergent evolution of alternative splicing. Furthermore, we showed that multiple sets of clade-conserved RNA pairings evolved to guide species-specific mutually exclusive splicing in Arthropoda. Importantly, we also identified a similar structural code in MRP exon clusters of the annelid, Capitella teleta, and chordate, Branchiostoma belcheri, suggesting an evolutionarily conserved competing pairing-guided mechanism in bilaterians. Taken together, these data reveal the molecular determinants and RNA pairing-guided evolution of species-specific mutually exclusive splicing spanning more than 600 million years of bilaterian evolution. These findings have a significant impact on our understanding of the evolution of and mechanism underpinning isoform diversity and complex gene structure.

Alternative splicing of U2AF1 reveals a shared repression mechanism for duplicated exons

The auxiliary factor of U2 small nuclear ribonucleoprotein (U2AF) facilitates branch point (BP) recognition and formation of lariat introns. The gene for the 35-kD subunit of U2AF gives rise to two protein isoforms (termed U2AF35a and U2AF35b) that are encoded by alternatively spliced exons 3 and Ab, respectively. The splicing recognition sequences of exon 3 are less favorable than exon Ab, yet U2AF35a expression is higher than U2AF35b across tissues. We show that U2AF35b repression is facilitated by weak, closely spaced BPs next to a long polypyrimidine tract of exon Ab. Each BP lacked canonical uridines at position -2 relative to the BP adenines, with efficient U2 base-pairing interactions predicted only for shifted registers reminiscent of programmed ribosomal frameshifting. The BP cluster was compensated by interactions involving unpaired cytosines in an upstream, EvoFold-predicted stem loop (termed ESL) that binds FUBP1/2. Exon Ab inclusion correlated with predicted free energies of mutant ESLs, suggesting that the ESL operates as a conserved rheostat between long inverted repeats upstream of each exon. The isoform-specific U2AF35 expression was U2AF65-dependent, required interactions between the U2AF-homology motif (UHM) and the α6 helix of U2AF35, and was fine-tuned by exon Ab/3 variants. Finally, we identify tandem homologous exons regulated by U2AF and show that their preferential responses to U2AF65-related proteins and SRSF3 are associated with unpaired pre-mRNA segments upstream of U2AF-repressed 3′ss. These results provide new insights into tissue-specific subfunctionalization of duplicated exons in vertebrate evolution and expand the repertoire of exon repression mechanisms that control alternative splicing.

RNA structure in splicing: An evolutionary perspective

Pre-mRNA splicing is a key post-transcriptional regulation process in which introns are excised and exons are ligated together. A novel class of structured intron was recently discovered in fish. Simple expansions of complementary AC and GT dimers at opposite boundaries of an intron were found to form a bridging structure, thereby enforcing correct splice site pairing across the intron. In some fish introns, the RNA structures are strong enough to bypass the need of regulatory protein factors for splicing. Here, we discuss the prevalence and potential functions of highly structured introns. In humans, structured introns usually arise through the co-occurrence of C and G-rich repeats at intron boundaries. We explore the potentially instructive example of the HLA receptor genes. In HLA pre-mRNA, structured introns flank the exons that encode the highly polymorphic β sheet cleft, making the processing of the transcript robust to variants that disrupt splicing factor binding. While selective forces that have shaped HLA receptor are fairly atypical, numerous other highly polymorphic genes that encode receptors contain structured introns. Finally, we discuss how the elevated mutation rate associated with the simple repeats that often compose structured intron can make structured introns themselves rapidly evolving elements.

Parallel computation of genome-scale RNA secondary structure to detect structural constraints on human genome

Background

RNA secondary structure around splice sites is known to assist normal splicing by promoting spliceosome recognition. However, analyzing the structural properties of entire intronic regions or pre-mRNA sequences has been difficult hitherto, owing to serious experimental and computational limitations, such as low read coverage and numerical problems.

Results

Our novel software, “ParasoR”, is designed to run on a computer cluster and enables the exact computation of various structural features of long RNA sequences under the constraint of maximal base-pairing distance. ParasoR divides dynamic programming (DP) matrices into smaller pieces, such that each piece can be computed by a separate computer node without losing the connectivity information between the pieces. ParasoR directly computes the ratios of DP variables to avoid the reduction of numerical precision caused by the cancellation of a large number of Boltzmann factors. The structural preferences of mRNAs computed by ParasoR shows a high concordance with those determined by high-throughput sequencing analyses.

Using ParasoR, we investigated the global structural preferences of transcribed regions in the human genome. A genome-wide folding simulation indicated that transcribed regions are significantly more structural than intergenic regions after removing repeat sequences and k-mer frequency bias. In particular, we observed a highly significant preference for base pairing over entire intronic regions as compared to their antisense sequences, as well as to intergenic regions. A comparison between pre-mRNAs and mRNAs showed that coding regions become more accessible after splicing, indicating constraints for translational efficiency. Such changes are correlated with gene expression levels, as well as GC content, and are enriched among genes associated with cytoskeleton and kinase functions.

Conclusions

We have shown that ParasoR is very useful for analyzing the structural properties of long RNA sequences such as mRNAs, pre-mRNAs, and long non-coding RNAs whose lengths can be more than a million bases in the human genome. In our analyses, transcribed regions including introns are indicated to be subject to various types of structural constraints that cannot be explained from simple sequence composition biases. ParasoR is freely available at https://github.com/carushi/ParasoR.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1067-9) contains supplementary material, which is available to authorized users.

Regulatory effects of cotranscriptional RNA structure formation and transitions

RNAs, which play significant roles in many fundamental biological processes of life, fold into sophisticated and precise structures. RNA folding is a dynamic and intricate process, which conformation transition of coding and noncoding RNAs form the primary elements of genetic regulation. The cellular environment contains various intrinsic and extrinsic factors that potentially affect RNA folding in vivo, and experimental and theoretical evidence increasingly indicates that the highly flexible features of the RNA structure are affected by these factors, which include the flanking sequence context, physiochemical conditions, cis RNA-RNA interactions, and RNA interactions with other molecules. Furthermore, distinct RNA structures have been identified that govern almost all steps of biological processes in cells, including transcriptional activation and termination, transcriptional mutagenesis, 5'-capping, splicing, 3'-polyadenylation, mRNA export and localization, and translation. Here, we briefly summarize the dynamic and complex features of RNA folding along with a wide variety of intrinsic and extrinsic factors that affect RNA folding. We then provide several examples to elaborate RNA structure-mediated regulation at the transcriptional and posttranscriptional levels. Finally, we illustrate the regulatory roles of RNA structure and discuss advances pertaining to RNA structure in plants. WIREs RNA 2016, 7:562-574. doi: 10.1002/wrna.1350 For further resources related to this article, please visit the WIREs website.

Structure Prediction: New Insights into Decrypting Long Noncoding RNAs

Long noncoding RNAs (lncRNAs), which form a diverse class of RNAs, remain the least understood type of noncoding RNAs in terms of their nature and identification. Emerging evidence has revealed that a small number of newly discovered lncRNAs perform important and complex biological functions such as dosage compensation, chromatin regulation, genomic imprinting, and nuclear organization. However, understanding the wide range of functions of lncRNAs related to various processes of cellular networks remains a great experimental challenge. Structural versatility is critical for RNAs to perform various functions and provides new insights into probing the functions of lncRNAs. In recent years, the computational method of RNA structure prediction has been developed to analyze the structure of lncRNAs. This novel methodology has provided basic but indispensable information for the rapid, large-scale and in-depth research of lncRNAs. This review focuses on mainstream RNA structure prediction methods at the secondary and tertiary levels to offer an additional approach to investigating the functions of lncRNAs.

RNA structure replaces the need for U2AF2 in splicing

RNA secondary structure plays an integral role in catalytic, ribosomal, small nuclear, micro, and transfer RNAs. Discovering a prevalent role for secondary structure in pre-mRNAs has proven more elusive. By utilizing a variety of computational and biochemical approaches, we present evidence for a class of nuclear introns that relies upon secondary structure for correct splicing. These introns are defined by simple repeat expansions of complementary AC and GT dimers that co-occur at opposite boundaries of an intron to form a bridging structure that enforces correct splice site pairing. Remarkably, this class of introns does not require U2AF2, a core component of the spliceosome, for its processing. Phylogenetic analysis suggests that this mechanism was present in the ancestral vertebrate lineage prior to the divergence of tetrapods from teleosts. While largely lost from land dwelling vertebrates, this class of introns is found in 10% of all zebrafish genes.

Long-range RNA pairings contribute to mutually exclusive splicing

Mutually exclusive splicing is an important means of increasing the protein repertoire, by which the Down's syndrome cell adhesion molecule (Dscam) gene potentially generates 38,016 different isoforms in Drosophila melanogaster. However, the regulatory mechanisms remain obscure due to the complexity of the Dscam exon cluster. Here, we reveal a molecular model for the regulation of the mutually exclusive splicing of the serpent pre-mRNA based on competition between upstream and downstream RNA pairings. Such dual RNA pairings confer fine tuning of the inclusion of alternative exons. Moreover, we demonstrate that the splicing outcome of alternative exons is mediated in relative pairing strength-correlated mode. Combined comparative genomics analysis and experimental evidence revealed similar bidirectional structural architectures in exon clusters 4 and 9 of the Dscam gene. Our findings provide a novel mechanistic framework for the regulation of mutually exclusive splicing and may offer potentially applicable insights into long-range RNA-RNA interactions in gene regulatory networks.

Nucleobase-Modified PNA Suppresses Translation by Forming a Triple Helix with a Hairpin Structure in mRNA In Vitro and in Cells

Compounds that bind specifically to double-stranded regions of RNA have potential as regulators of structure-based RNA function; however, sequence-selective recognition of double-stranded RNA is challenging. The modification of peptide nucleic acid (PNA) with unnatural nucleobases enables the formation of PNA-RNA triplexes. Herein, we demonstrate that a 9-mer PNA forms a sequence-specific PNA-RNA triplex with a dissociation constant of less than 1 nm at physiological pH. The triplex formed within the 5' untranslated region of an mRNA reduces the protein expression levels both in vitro and in cells. A single triplet mismatch destabilizes the complex, and in this case, no translation suppression is observed. The triplex-forming PNAs are unique and potent compounds that hold promise as inhibitors of cellular functions that are controlled by double-stranded RNAs, such as RNA interference, RNA editing, and RNA localization mediated by protein-RNA interactions.

Mechanisms and Regulation of Alternative Pre-mRNA Splicing

Precursor messenger RNA (pre-mRNA) splicing is a critical step in the posttranscriptional regulation of gene expression, providing significant expansion of the functional proteome of eukaryotic organisms with limited gene numbers. Split eukaryotic genes contain intervening sequences or introns disrupting protein-coding exons, and intron removal occurs by repeated assembly of a large and highly dynamic ribonucleoprotein complex termed the spliceosome, which is composed of five small nuclear ribonucleoprotein particles, U1, U2, U4/U6, and U5. Biochemical studies over the past 10 years have allowed the isolation as well as compositional, functional, and structural analysis of splicing complexes at distinct stages along the spliceosome cycle. The average human gene contains eight exons and seven introns, producing an average of three or more alternatively spliced mRNA isoforms. Recent high-throughput sequencing studies indicate that 100% of human genes produce at least two alternative mRNA isoforms. Mechanisms of alternative splicing include RNA-protein interactions of splicing factors with regulatory sites termed silencers or enhancers, RNA-RNA base-pairing interactions, or chromatin-based effects that can change or determine splicing patterns. Disease-causing mutations can often occur in splice sites near intron borders or in exonic or intronic RNA regulatory silencer or enhancer elements, as well as in genes that encode splicing factors. Together, these studies provide mechanistic insights into how spliceosome assembly, dynamics, and catalysis occur; how alternative splicing is regulated and evolves; and how splicing can be disrupted by cis- and trans-acting mutations leading to disease states. These findings make the spliceosome an attractive new target for small-molecule, antisense, and genome-editing therapeutic interventions.