The architecture of pre-mRNAs affects mechanisms of splice-site pairing

Multimodal screen identifies noise-regulatory proteins

Gene-expression noise can influence cell-fate choices across pathology and physiology. However, a crucial question persists: do regulatory proteins or pathways exist that control noise independently of mean expression levels? Our integrative approach, combining single-cell RNA sequencing with proteomics and regulator enrichment analysis, identifies 32 putative noise regulators. SON, a nuclear speckle-associated protein, alters transcriptional noise without changing mean expression levels. Furthermore, SON's noise control can propagate to the protein level. Long-read and total RNA sequencing shows that SON's noise control does not significantly change isoform usage or splicing efficiency. Moreover, SON depletion reduces state switching in pluripotent mouse embryonic stem cells and impacts their fate choice during differentiation. Collectively, we demonstrate a class of proteins that control noise orthogonally to mean expression levels. This work serves as a proof of concept that can identify other functional noise regulators throughout development and disease progression.

Timing is everything: advances in quantifying splicing kinetics

Splicing is a highly regulated process critical for proper pre-mRNA maturation and the maintenance of a healthy cellular environment. Splicing events are impacted by ongoing transcription, neighboring splicing events, and cis and trans regulatory factors on the respective pre-mRNA transcript. Within this complex regulatory environment, splicing kinetics have the potential to influence splicing outcomes but have historically been challenging to study in vivo. In this review, we highlight recent technological advancements that have enabled measurements of global splicing kinetics and of the variability of splicing kinetics at single introns. We demonstrate how identifying features that are correlated with splicing kinetics has increased our ability to form potential models for how splicing kinetics may be regulated in vivo.

Emerging and re-emerging themes in co-transcriptional pre-mRNA splicing

Proper gene expression requires the collaborative effort of multiple macromolecular machines to produce functional messenger RNA. As RNA polymerase II (RNA Pol II) transcribes DNA, the nascent pre-messenger RNA is heavily modified by other complexes such as 5' capping enzymes, the spliceosome, the cleavage, and polyadenylation machinery as well as RNA-modifying/editing enzymes. Recent evidence has demonstrated that pre-mRNA splicing and 3' end cleavage can occur on similar timescales as transcription and significantly cross-regulate. In this review, we discuss recent advances in co-transcriptional processing and how it contributes to gene regulation. We highlight how emerging areas-including coordinated splicing events, physical interactions between the RNA synthesis and modifying machinery, rapid and delayed splicing, and nuclear organization-impact mRNA isoforms. Coordination among RNA-processing choices yields radically different mRNA and protein products, foreshadowing the likely regulatory importance of co-transcriptional RNA folding and co-transcriptional modifications that have yet to be characterized in detail.

Expression of Truncated Products at the 5'-Terminal Region of RIPK2 and Evolutive Aspects that Support Their Biological Importance

Alternative splicing is the process of generating different mRNAs from the same primary transcript, which contributes to increase the transcriptome and proteome diversity. Abnormal splicing has been associated with the development of several diseases including cancer. Given that mutations and abnormal levels of the RIPK2 transcript and RIP-2 protein are frequent in tumors, and that RIP-2 modulates immune and inflammatory responses, we investigated alternative splicing events that result in partial deletions of the kinase domain at the N-terminus of RIP-2. We also investigated the structure and expression of the RIPK2 truncated variants and isoforms in different environments. In addition, we searched data throughout Supraprimates evolution that could support the biological importance of RIPK2 alternatively spliced products. We observed that human variants and isoforms were differentially regulated following temperature stress, and that the truncated transcript was more expressed than the long transcript in tumor samples. The inverse was found for the longer protein isoform. The truncated variant was also detected in chimpanzee, gorilla, hare, pika, mouse, rat, and tree shrew. The fact that the same variant has been preserved in mammals with divergence times up to 70 million years raises the hypothesis that it may have a functional significance.

Convergent and/or parallel evolution of RNA-binding proteins in angiosperms after polyploidization

Increasing studies suggest that the biased retention of stress-related transcription factors (TFs) after whole-genome duplications (WGDs) could rewire gene transcriptional networks, facilitating plant adaptation to challenging environments. However, the role of posttranscriptional factors (e.g. RNA-binding proteins, RBPs) following WGDs has been largely ignored. Uncovering thousands of RBPs in 21 representative angiosperm species, we integrate genomic, transcriptomic, regulatomic, and paleotemperature datasets to unravel their evolutionary trajectories and roles in adapting to challenging environments. We reveal functional enrichments of RBP genes in stress responses and identify their convergent retention across diverse angiosperms from independent WGDs, coinciding with global cooling periods. Numerous RBP duplicates derived from WGDs are then identified as cold-induced. A significant overlap of 29 orthogroups between WGD-derived and cold-induced RBP genes across diverse angiosperms highlights a correlation between WGD and cold stress. Notably, we unveil an orthogroup (Glycine-rich RNA-binding Proteins 7/8, GRP7/8) and relevant TF duplicates (CCA1/LHY, RVE4/8, CBF2/4, etc.), co-retained in different angiosperms post-WGDs. Finally, we illustrate their roles in rewiring circadian and cold-regulatory networks at both transcriptional and posttranscriptional levels during global cooling. Altogether, we underline the adaptive evolution of RBPs in angiosperms after WGDs during global cooling, improving our understanding of plants surviving periods of environmental turmoil.

Splicing defects in rare diseases: transcriptomics and machine learning strategies towards genetic diagnosis

Genomic variants affecting pre-messenger RNA splicing and its regulation are known to underlie many rare genetic diseases. However, common workflows for genetic diagnosis and clinical variant interpretation frequently overlook splice-altering variants. To better serve patient populations and advance biomedical knowledge, it has become increasingly important to develop and refine approaches for detecting and interpreting pathogenic splicing variants. In this review, we will summarize a few recent developments and challenges in using RNA sequencing technologies for rare disease investigation. Moreover, we will discuss how recent computational splicing prediction tools have emerged as complementary approaches for revealing disease-causing variants underlying splicing defects. We speculate that continuous improvements to sequencing technologies and predictive modeling will not only expand our understanding of splicing regulation but also bring us closer to filling the diagnostic gap for rare disease patients.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green plants.

Pre-mRNA splicing and its cotranscriptional connections

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

FUBP1 is a general splicing factor facilitating 3' splice site recognition and splicing of long introns

Splicing of pre-mRNAs critically contributes to gene regulation and proteome expansion in eukaryotes, but our understanding of the recognition and pairing of splice sites during spliceosome assembly lacks detail. Here, we identify the multidomain RNA-binding protein FUBP1 as a key splicing factor that binds to a hitherto unknown cis-regulatory motif. By collecting NMR, structural, and in vivo interaction data, we demonstrate that FUBP1 stabilizes U2AF2 and SF1, key components at the 3' splice site, through multivalent binding interfaces located within its disordered regions. Transcriptional profiling and kinetic modeling reveal that FUBP1 is required for efficient splicing of long introns, which is impaired in cancer patients harboring FUBP1 mutations. Notably, FUBP1 interacts with numerous U1 snRNP-associated proteins, suggesting a unique role for FUBP1 in splice site bridging for long introns. We propose a compelling model for 3' splice site recognition of long introns, which represent 80% of all human introns.

Introns: the "dark matter" of the eukaryotic genome

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the "dark matter" of the eukaryotic genome and discuss some of the outstanding questions in the field.

Integrative multi-omics analysis of three early diverged rosid species reveals an ancient hierarchical cold-responsive regulatory network

Elucidating regulators, including transcription factors (TFs) and RNA-binding proteins (RBPs), underlying gene transcriptional and post-transcriptional co-regulatory network is key to understand plant cold responses. Previous studies were mainly conducted on single species, and whether the regulators are conserved across different species remains elusive. Here, we selected three species that diverged at the early evolution of rosids (~99-113 million years ago), performed cold-responsive phylotranscriptome experiments, and integrated chromatin immunoprecipitation- and DNA affinity purification-sequencing (ChIP/DAP-seq) analysis to explore cold-responsive regulators and their regulatory networks. First, we detected over 10,000 cold-induced differentially expressed genes (DEGs) and alternative splicing genes (DASGs) in each species. Among the DEGs, a set of TFs and RBPs were conserved in rosid cold response. Compared to TFs, RBPs displayed a delayed cold-responsive pattern, implying a hierarchical regulation of DEGs and DASGs. By integrating DEGs and DASGs, we identified 259 overlapping DE-DASG orthogroups (closely-related homologs) that were cold-regulated at both transcriptional and post-transcriptional levels in all three studied species. Notably, pathway analysis on each of the DEGs, DASGs, and DE-DASGs in the three species showed a common enrichment connected to the circadian rhythm. Evidently, 226 cold-responsive genes were directly targeted by at least two circadian rhythm components (CCA1, LHY, RV4, RVE7, and RVE8). Finally, we revealed an ancient hierarchy of cold-responsive regulatory networks at transcriptional and post-transcriptional levels launched by circadian components in rosids. Altogether, this study sheds light on conserved regulators underlying cold-responsive regulatory networks across rosid species, despite a long evolutionary history after their divergence.

Profiling lariat intermediates reveals genetic determinants of early and late co-transcriptional splicing

Long introns with short exons in vertebrate genes are thought to require spliceosome assembly across exons (exon definition), rather than introns, thereby requiring transcription of an exon to splice an upstream intron. Here, we developed CoLa-seq (co-transcriptional lariat sequencing) to investigate the timing and determinants of co-transcriptional splicing genome wide. Unexpectedly, 90% of all introns, including long introns, can splice before transcription of a downstream exon, indicating that exon definition is not obligatory for most human introns. Still, splicing timing varies dramatically across introns, and various genetic elements determine this variation. Strong U2AF2 binding to the polypyrimidine tract predicts early splicing, explaining exon definition-independent splicing. Together, our findings question the essentiality of exon definition and reveal features beyond intron and exon length that are determinative for splicing timing.

The nucleolus is the site for inflammatory RNA decay during infection

Inflammatory cytokines are key signaling molecules that can promote an immune response, thus their RNA turnover must be tightly controlled during infection. Most studies investigate the RNA decay pathways in the cytosol or nucleoplasm but never focused on the nucleolus. Although this organelle has well-studied roles in ribosome biogenesis and cellular stress sensing, the mechanism of RNA decay within the nucleolus is not completely understood. Here, we report that the nucleolus is an essential site of inflammatory pre-mRNA instability during infection. RNA-sequencing analysis reveals that not only do inflammatory genes have higher intronic read densities compared with non-inflammatory genes, but their pre-mRNAs are highly enriched in nucleoli during infection. Notably, nucleolin (NCL) acts as a guide factor for recruiting cytosine or uracil (C/U)-rich sequence-containing inflammatory pre-mRNAs and the Rrp6-exosome complex to the nucleolus through a physical interaction, thereby enabling targeted RNA delivery to Rrp6-exosomes and subsequent degradation. Consequently, Ncl depletion causes aberrant hyperinflammation, resulting in a severe lethality in response to LPS. Importantly, the dynamics of NCL post-translational modifications determine its functional activity in phases of LPS. This process represents a nucleolus-dependent pathway for maintaining inflammatory gene expression integrity and immunological homeostasis during infection.

The nucleolus is the traditional site for ribosomal RNA biogenesis. Here, the authors find that the nucleolus is a site of inflammatory pre-mRNA turnover and elucidated how immune homeostasis can be maintained by controlling inflammatory gene expression.

Gene architecture directs splicing outcome in separate nuclear spatial regions

How the splicing machinery defines exons or introns as the spliced unit has remained a puzzle for 30 years. Here, we demonstrate that peripheral and central regions of the nucleus harbor genes with two distinct exon-intron GC content architectures that differ in the splicing outcome. Genes with low GC content exons, flanked by long introns with lower GC content, are localized in the periphery, and the exons are defined as the spliced unit. Alternative splicing of these genes results in exon skipping. In contrast, the nuclear center contains genes with a high GC content in the exons and short flanking introns. Most splicing of these genes occurs via intron definition, and aberrant splicing leads to intron retention. We demonstrate that the nuclear periphery and center generate different environments for the regulation of alternative splicing and that two sets of splicing factors form discrete regulatory subnetworks for the two gene architectures. Our study connects 3D genome organization and splicing, thus demonstrating that exon and intron definition modes of splicing occur in different nuclear regions.

Transcription and splicing dynamics during early Drosophila development

Widespread cotranscriptional splicing has been demonstrated from yeast to human. However, most studies to date addressing the kinetics of splicing relative to transcription used either Saccharomyces cerevisiae or metazoan cultured cell lines. Here, we adapted native elongating transcript sequencing technology (NET-seq) to measure cotranscriptional splicing dynamics during the early developmental stages of Drosophila melanogaster embryos. Our results reveal the position of RNA polymerase II (Pol II) when both canonical and recursive splicing occur. We found heterogeneity in splicing dynamics, with some RNAs spliced immediately after intron transcription, whereas for other transcripts no splicing was observed over the first 100 nt of the downstream exon. Introns that show splicing completion before Pol II has reached the end of the downstream exon are necessarily intron-defined. We studied the splicing dynamics of both nascent pre-mRNAs transcribed in the early embryo, which have few and short introns, as well as pre-mRNAs transcribed later in embryonic development, which contain multiple long introns. As expected, we found a relationship between the proportion of spliced reads and intron size. However, intron definition was observed at all intron sizes. We further observed that genes transcribed in the early embryo tend to be isolated in the genome whereas genes transcribed later are often overlapped by a neighboring convergent gene. In isolated genes, transcription termination occurred soon after the polyadenylation site, while in overlapped genes, Pol II persisted associated with the DNA template after cleavage and polyadenylation of the nascent transcript. Taken together, our data unravel novel dynamic features of Pol II transcription and splicing in the developing Drosophila embryo.

Exon-independent recruitment of SRSF1 is mediated by U1 snRNP stem-loop 3

SRSF1 protein and U1 snRNPs are closely connected splicing factors. They both stimulate exon inclusion, SRSF1 by binding to exonic splicing enhancer sequences (ESEs) and U1 snRNPs by binding to the downstream 5' splice site (SS), and both factors affect 5' SS selection. The binding of U1 snRNPs initiates spliceosome assembly, but SR proteins such as SRSF1 can in some cases substitute for it. The mechanistic basis of this relationship is poorly understood. We show here by single-molecule methods that a single molecule of SRSF1 can be recruited by a U1 snRNP. This reaction is independent of exon sequences and separate from the U1-independent process of binding to an ESE. Structural analysis and cross-linking data show that SRSF1 contacts U1 snRNA stem-loop 3, which is required for splicing. We suggest that the recruitment of SRSF1 to a U1 snRNP at a 5'SS is the basis for exon definition by U1 snRNP and might be one of the principal functions of U1 snRNPs in the core reactions of splicing in mammals.

Splice site proximity influences alternative exon definition

Alternative splicing enables higher eukaryotes to expand mRNA diversity from a finite number of genes through highly combinatorial splice site selection mechanisms that are influenced by the sequence of competing splice sites, cis-regulatory elements binding trans-acting factors, the length of exons and introns harbouring alternative splice sites and RNA secondary structures at putative splice junctions. To test the hypothesis that the intron definition or exon definition modes of splice site recognition direct the selection of alternative splice patterns, we created a database of alternative splice site usage (ALTssDB). When alternative splice sites are embedded within short introns (intron definition), the 5' and 3' splice sites closest to each other across the intron preferentially pair, consistent with previous observations. However, when alternative splice sites are embedded within large flanking introns (exon definition), the 5' and 3' splice sites closest to each other across the exon are preferentially selected. Thus, alternative splicing decisions are influenced by the intron and exon definition modes of splice site recognition. The results demonstrate that the spliceosome pairs splice sites that are closest in proximity within the unit of initial splice site selection.

Alternative splicing level related to intron size and organism complexity

Background

Alternative splicing is the process of selecting different combinations of splice sites to produce variably spliced mRNAs. However, the relationships between alternative splicing prevalence and level (ASP/L) and variations of intron size and organism complexity (OC) remain vague. Here, we developed a robust protocol to analyze the relationships between ASP/L and variations of intron size and OC. Approximately 8 Tb raw RNA-Seq data from 37 eumetazoan species were divided into three sets of species based on variations in intron size and OC.

Results

We found a strong positive correlation between ASP/L and OC, but no correlation between ASP/L and intron size across species. Surprisingly, ASP/L displayed a positive correlation with mean intron size of genes within individual genomes. Moreover, our results revealed that four ASP/L-related pathways contributed to the differences in ASP/L that were associated with OC. In particular, the spliceosome pathway displayed distinct genomic features, such as the highest gene expression level, conservation level, and fraction of disordered regions. Interestingly, lower or no obvious correlations were observed among these genomic features.

Conclusions

The positive correlation between ASP/L and OC ubiquitously exists in eukaryotes, and this correlation is not affected by the mean intron size of these species. ASP/L-related splicing factors may play an important role in the evolution of OC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08172-2.

Influenza A virus segments five and six can harbor artificial introns allowing expanded coding capacity

Influenza A viruses encode their genomes across eight, negative sense RNA segments. The six largest segments produce mRNA transcripts that do not generally splice; however, the two smallest segments are actively spliced to produce the essential viral proteins NEP and M2. Thus, viral utilization of RNA splicing effectively expands the viral coding capacity without increasing the number of genomic segments. As a first step towards understanding why splicing is not more broadly utilized across genomic segments, we designed and inserted an artificial intron into the normally nonsplicing NA segment. This insertion was tolerated and, although viral mRNAs were incompletely spliced, we observed only minor effects on viral fitness. To take advantage of the unspliced viral RNAs, we encoded a reporter luciferase gene in frame with the viral ORF such that when the intron was not removed the reporter protein would be produced. This approach, which we also show can be applied to the NP encoding segment and in different viral genetic backgrounds, led to high levels of reporter protein expression with minimal effects on the kinetics of viral replication or the ability to cause disease in experimentally infected animals. These data together show that the influenza viral genome is more tolerant of splicing than previously appreciated and this knowledge can be leveraged to develop viral genetic platforms with utility for biotechnology applications.

Unlike most host mRNAs, some viral mRNAs encode multiple discrete, functional proteins. One method influenza A viruses use to increase the protein products from two of their eight RNA genome segments is splicing. Splicing requires host machinery to remove part of the viral mRNA, the intron, to generate a different mRNA product. Although only certain influenza viral segments naturally splice, we were interested in whether additional segments could splice to produce multiple proteins. We inserted artificial introns harboring reporter genes into otherwise nonsplicing genomic segments of an H1N1 influenza A virus and found that this modification was well tolerated by the virus. We further demonstrated that an unrelated H3N2 influenza A virus could similarly support splicing and express a reporter protein from an artificial intron. These findings have implications for our understanding of how viruses expand their coding capacity with a limited genome. Additionally, encoding reporter proteins in spliced intronic sequences also represents a new method of generating reporter viruses requiring limited manipulation of the viral RNA.

High Concentration of an ISS-N1-Targeting Antisense Oligonucleotide Causes Massive Perturbation of the Transcriptome

Intronic splicing silencer N1 (ISS-N1) located within Survival Motor Neuron 2 (SMN2) intron 7 is the target of a therapeutic antisense oligonucleotide (ASO), nusinersen (Spinraza), which is currently being used for the treatment of spinal muscular atrophy (SMA), a leading genetic disease associated with infant mortality. The discovery of ISS-N1 as a promising therapeutic target was enabled in part by Anti-N1, a 20-mer ASO that restored SMN2 exon 7 inclusion by annealing to ISS-N1. Here, we analyzed the transcriptome of SMA patient cells treated with 100 nM of Anti-N1 for 30 h. Such concentrations are routinely used to demonstrate the efficacy of an ASO. While 100 nM of Anti-N1 substantially stimulated SMN2 exon 7 inclusion, it also caused massive perturbations in the transcriptome and triggered widespread aberrant splicing, affecting expression of essential genes associated with multiple cellular processes such as transcription, splicing, translation, cell signaling, cell cycle, macromolecular trafficking, cytoskeletal dynamics, and innate immunity. We validated our findings with quantitative and semiquantitative PCR of 39 candidate genes associated with diverse pathways. We also showed a substantial reduction in off-target effects with shorter ISS-N1-targeting ASOs. Our findings are significant for implementing better ASO design and dosing regimens of ASO-based drugs.

Spliceosomal component PRP-40 is a central regulator of microexon splicing

Microexons (≤27 nt) play critical roles in nervous system development and function but create unique challenges for the splicing machinery. The mechanisms of microexon regulation are therefore of great interest. We performed a genetic screen for alternative splicing regulators in the C. elegans nervous system and identify PRP-40, a core component of the U1 snRNP. RNA-seq reveals that PRP-40 is required for inclusion of alternatively spliced, but not constitutively spliced, exons. PRP-40 is particularly required for inclusion of neuronal microexons, and our data indicate that PRP-40 is a central regulator of microexon splicing. Microexons can be relieved from PRP-40 dependence by artificially increasing exon size or reducing flanking intron size, indicating that PRP-40 is specifically required for microexons surrounded by conventionally sized introns. Knockdown of the orthologous PRPF40A in mouse neuroblastoma cells causes widespread dysregulation of microexons but not conventionally sized exons. PRP-40 regulation of neuronal microexons is therefore a widely conserved phenomenon.

Induction of cryptic pre-mRNA splice-switching by antisense oligonucleotides

Antisense oligomers (AOs) are increasingly being used to modulate RNA splicing in live cells, both for research and for the development of therapeutics. While the most common intended effect of these AOs is to induce skipping of whole exons, rare examples are emerging of AOs that induce skipping of only part of an exon, through activation of an internal cryptic splice site. In this report, we examined seven AO-induced cryptic splice sites in six genes. Five of these cryptic splice sites were discovered through our own experiments, and two originated from other published reports. We modelled the predicted effects of AO binding on the secondary structure of each of the RNA targets, and how these alterations would in turn affect the accessibility of the RNA to splice factors. We observed that a common predicted effect of AO binding was disruption of the exon definition signal within the exon's excluded segment.

Nonproductive Splicing Prevents Expression of MYH7b Protein in the Mammalian Heart

Background

Although the roles of alpha‐myosin heavy chain (α‐MyHC) and beta‐myosin heavy chain (β‐MyHC) proteins in cardiac contractility have long been appreciated, the biological contribution of another closely related sarcomeric myosin family member, MYH7b (myosin heavy chain 7b), has become a matter of debate. In mammals, MYH7b mRNA is transcribed but undergoes non‐productive alternative splicing that prevents protein expression in a tissue‐specific manner, including in the heart. However, several studies have recently linked MYH7b variants to different cardiomyopathies or have reported MYH7b protein expression in mammalian hearts.

Methods and Results

By analyzing mammalian cardiac transcriptome and proteome data, we show that the vast majority of MYH7b RNA is subject to exon skipping and cannot be translated into a functional myosin molecule. Notably, we discovered a lag in the removal of introns flanking the alternatively spliced exon, which could retain the non‐coding RNA in the nucleus. This process could play a significant role in controlling MYH7b expression as well as the activity of other cardiac genes. Consistent with the negligible level of full‐length protein coding mRNA, no MYH7b protein expression was detected in adult mouse, rat, and human hearts by Western blot analysis. Furthermore, proteome surveys including quantitative mass spectrometry analyses revealed only traces of cardiac MYH7b protein and even then, only in a subset of individual samples.

Conclusions

The comprehensive analysis presented here suggests that previous studies showing cardiac MYH7b protein expression were likely attributable to antibody cross‐reactivity. More importantly, our data predict that the MYH7b disease‐associated variants may operate through the alternately spliced RNA itself.

Base editing repairs an SGCA mutation in human primary muscle stem cells

Skeletal muscle can regenerate from muscle stem cells and their myogenic precursor cell progeny, myoblasts. However, precise gene editing in human muscle stem cells for autologous cell replacement therapies of untreatable genetic muscle diseases has not yet been reported. Loss-of-function mutations in SGCA, encoding α-sarcoglycan, cause limb-girdle muscular dystrophy 2D/R3, an early-onset, severe, and rapidly progressive form of muscular dystrophy affecting both male and female patients. Patients suffer from muscle degeneration and atrophy affecting the limbs, respiratory muscles, and heart. We isolated human muscle stem cells from 2 donors, with the common SGCA c.157G>A mutation affecting the last coding nucleotide of exon 2. We found that c.157G>A is an exonic splicing mutation that induces skipping of 2 coregulated exons. Using adenine base editing, we corrected the mutation in the cells from both donors with > 90% efficiency, thereby rescuing the splicing defect and α-sarcoglycan expression. Base-edited patient cells regenerated muscle and contributed to the Pax7⁺ satellite cell compartment in vivo in mouse xenografts. Here, we provide the first evidence to our knowledge that autologous gene–repaired human muscle stem cells can be harnessed for cell replacement therapies of muscular dystrophies.

Disruption of exon-bridging interactions between the minor and major spliceosomes results in alternative splicing around minor introns

Vertebrate genomes contain major (>99.5%) and minor (<0.5%) introns that are spliced by the major and minor spliceosomes, respectively. Major intron splicing follows the exon-definition model, whereby major spliceosome components first assemble across exons. However, since most genes with minor introns predominately consist of major introns, formation of exon-definition complexes in these genes would require interaction between the major and minor spliceosomes. Here, we report that minor spliceosome protein U11-59K binds to the major spliceosome U2AF complex, thereby supporting a model in which the minor spliceosome interacts with the major spliceosome across an exon to regulate the splicing of minor introns. Inhibition of minor spliceosome snRNAs and U11-59K disrupted exon-bridging interactions, leading to exon skipping by the major spliceosome. The resulting aberrant isoforms contained a premature stop codon, yet were not subjected to nonsense-mediated decay, but rather bound to polysomes. Importantly, we detected elevated levels of these alternatively spliced transcripts in individuals with minor spliceosome-related diseases such as Roifman syndrome, Lowry–Wood syndrome and early-onset cerebellar ataxia. In all, we report that the minor spliceosome informs splicing by the major spliceosome through exon-definition interactions and show that minor spliceosome inhibition results in aberrant alternative splicing in disease.

Transcription initiation of distant core promoters in a large-sized genome of an insect

Background

Core promoters have a substantial influence on various steps of transcription, including initiation, elongation, termination, polyadenylation, and finally, translation. The characterization of core promoters is crucial for exploring the regulatory code of transcription initiation. However, the current understanding of insect core promoters is focused on those of Diptera (especially Drosophila) species with small genome sizes.

Results

Here, we present an analysis of the transcription start sites (TSSs) in the migratory locust, Locusta migratoria, which has a genome size of 6.5 Gb. The genomic differences, including lower precision of transcription initiation and fewer constraints on the distance from transcription factor binding sites or regulatory elements to TSSs, were revealed in locusts compared with Drosophila insects. Furthermore, we found a distinct bimodal log distribution of the distances from the start codons to the core promoters of locust genes. We found stricter constraints on the exon length of mRNA leaders and widespread expression activity of the distant core promoters in locusts compared with fruit flies. We further compared core promoters in seven arthropod species across a broad range of genome sizes to reinforce our results on the emergence of distant core promoters in large-sized genomes.

Conclusions

In summary, our results provide novel insights into the effects of genome size expansion on distant transcription initiation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01004-5.

Role of promoters in regulating alternative splicing

Alternative splicing (AS) plays a critical role in enhancing proteome complexity in higher eukaryotes. Almost all the multi intron-containing genes undergo AS in humans. Splicing mainly occurs co-transcriptionally, where RNA polymerase II (RNA pol II) plays a crucial role in coordinating transcription and pre-mRNA splicing. Aberrant AS leads to non-functional proteins causative in various pathophysiological conditions such as cancers, neurodegenerative diseases, and muscular dystrophies. Transcription and pre-mRNA splicing are deeply interconnected and can influence each other's functions. Several studies evinced that specific promoters employed by RNA pol II dictate the RNA processing decisions. Promoter-specific recruitment of certain transcriptional factors or transcriptional coactivators influences splicing, and the extent to which these factors affect splicing has not been discussed in detail. Here, in this review, various DNA-binding proteins and their influence on promoter-specific AS are extensively discussed. Besides, this review highlights how the promoter-specific epigenetic changes might regulate AS.

NineTeen Complex-subunit Salsa is required for efficient splicing of a subset of introns and dorsal-ventral patterning

The NineTeen Complex (NTC), also known as pre-mRNA-processing factor 19 (Prp19) complex, regulates distinct spliceosome conformational changes necessary for splicing. During Drosophila midblastula transition, splicing is particularly sensitive to mutations in NTC-subunit Fandango, which suggests differential requirements of NTC during development. We show that NTC-subunit Salsa, the Drosophila ortholog of human RNA helicase Aquarius, is rate-limiting for splicing of a subset of small first introns during oogenesis, including the first intron of gurken Germline depletion of Salsa and splice site mutations within gurken first intron impair both adult female fertility and oocyte dorsal-ventral patterning, due to an abnormal expression of Gurken. Supporting causality, the fertility and dorsal-ventral patterning defects observed after Salsa depletion could be suppressed by the expression of a gurken construct without its first intron. Altogether, our results suggest that one of the key rate-limiting functions of Salsa during oogenesis is to ensure the correct expression and efficient splicing of the first intron of gurken mRNA. Retention of gurken first intron compromises the function of this gene most likely because it undermines the correct structure and function of the transcript 5'UTR.

STRAP regulates alternative splicing fidelity during lineage commitment of mouse embryonic stem cells

Alternative splicing (AS) is involved in cell fate decisions and embryonic development. However, regulation of these processes is poorly understood. Here, we have identified the serine threonine kinase receptor-associated protein (STRAP) as a putative spliceosome-associated factor. Upon Strap deletion, there are numerous AS events observed in mouse embryoid bodies (EBs) undergoing a neuroectoderm-like state. Global mapping of STRAP-RNA binding in mouse embryos by enhanced-CLIP sequencing (eCLIP-seq) reveals that STRAP preferably targets transcripts for nervous system development and regulates AS through preferred binding positions, as demonstrated for two neuronal-specific genes, Nnat and Mark3. We have found that STRAP involves in the assembly of 17S U2 snRNP proteins. Moreover, in Xenopus, loss of Strap leads to impeded lineage differentiation in embryos, delayed neural tube closure, and altered exon skipping. Collectively, our findings reveal a previously unknown function of STRAP in mediating the splicing networks of lineage commitment, alteration of which may be involved in early embryonic lethality in mice.

Global regulatory features of alternative splicing across tissues and within the nervous system of C. elegans

Alternative splicing plays a major role in shaping tissue-specific transcriptomes. Among the broad tissue types present in metazoans, the central nervous system contains some of the highest levels of alternative splicing. Although many documented examples of splicing differences between broad tissue types exist, there remains much to be understood about the splicing factors and the cis sequence elements controlling tissue and neuron subtype-specific splicing patterns. By using translating ribosome affinity purification coupled with deep-sequencing (TRAP-seq) in Caenorhabditis elegans, we have obtained high coverage profiles of ribosome-associated mRNA for three broad tissue classes (nervous system, muscle, and intestine) and two neuronal subtypes (dopaminergic and serotonergic neurons). We have identified hundreds of splice junctions that exhibit distinct splicing patterns between tissue types or within the nervous system. Alternative splicing events differentially regulated between tissues are more often frame-preserving, are more highly conserved across Caenorhabditis species, and are enriched in specific cis regulatory motifs, when compared with other types of exons. By using this information, we have identified a likely mechanism of splicing repression by the RNA-binding protein UNC-75/CELF via interactions with cis elements that overlap a 5′ splice site. Alternatively spliced exons also overlap more frequently with intrinsically disordered peptide regions than constitutive exons. Moreover, regulated exons are often shorter than constitutive exons but are flanked by longer intron sequences. Among these tissue-regulated exons are several highly conserved microexons <27 nt in length. Collectively, our results indicate a rich layer of tissue-specific gene regulation at the level of alternative splicing in C. elegans that parallels the evolutionary forces and constraints observed across metazoa.

Large-scale comparative evaluation of user-friendly tools for predicting variant-induced alterations of splicing regulatory elements

Discriminating which nucleotide variants cause disease or contribute to phenotypic traits remains a major challenge in human genetics. In theory, any intragenic variant can potentially affect RNA splicing by altering splicing regulatory elements (SREs). However, these alterations are often ignored mainly because pioneer SRE predictors have proved inefficient. Here, we report the first large-scale comparative evaluation of four user-friendly SRE-dedicated algorithms (QUEPASA, HEXplorer, SPANR, and HAL) tested both as standalone tools and in multiple combined ways based on two independent benchmark datasets adding up to >1,300 exonic variants studied at the messenger RNA level and mapping to 89 different disease-causing genes. These methods display good predictive power, based on decision thresholds derived from the receiver operating characteristics curve analyses, with QUEPASA and HAL having the best accuracies either as standalone or in combination. Still, overall there was a tight race between the four predictors, suggesting that all methods may be of use. Additionally, QUEPASA and HEXplorer may be beneficial as well for predicting variant-induced creation of pseudoexons deep within introns. Our study highlights the potential of SRE predictors as filtering tools for identifying disease-causing candidates among the plethora of variants detected by high-throughput DNA sequencing and provides guidance for their use in genomic medicine settings.

Comprehensive database and evolutionary dynamics of U12-type introns

During nuclear maturation of most eukaryotic pre-messenger RNAs and long non-coding RNAs, introns are removed through the process of RNA splicing. Different classes of introns are excised by the U2-type or the U12-type spliceosomes, large complexes of small nuclear ribonucleoprotein particles and associated proteins. We created intronIC, a program for assigning intron class to all introns in a given genome, and used it on 24 eukaryotic genomes to create the Intron Annotation and Orthology Database (IAOD). We then used the data in the IAOD to revisit several hypotheses concerning the evolution of the two classes of spliceosomal introns, finding support for the class conversion model explaining the low abundance of U12-type introns in modern genomes.

Hypoxia-induced alternative splicing: the 11th Hallmark of Cancer

Hypoxia-induced alternative splicing is a potent driving force in tumour pathogenesis and progression. In this review, we update currents concepts of hypoxia-induced alternative splicing and how it influences tumour biology. Following brief descriptions of tumour-associated hypoxia and the pre-mRNA splicing process, we review the many ways hypoxia regulates alternative splicing and how hypoxia-induced alternative splicing impacts each individual hallmark of cancer. Hypoxia-induced alternative splicing integrates chemical and cellular tumour microenvironments, underpins continuous adaptation of the tumour cellular microenvironment responsible for metastatic progression and plays clear roles in oncogene activation and autonomous tumour growth, tumor suppressor inactivation, tumour cell immortalization, angiogenesis, tumour cell evasion of programmed cell death and the anti-tumour immune response, a tumour-promoting inflammatory response, adaptive metabolic re-programming, epithelial to mesenchymal transition, invasion and genetic instability, all of which combine to promote metastatic disease. The impressive number of hypoxia-induced alternative spliced protein isoforms that characterize tumour progression, classifies hypoxia-induced alternative splicing as the 11th hallmark of cancer, and offers a fertile source of potential diagnostic/prognostic markers and therapeutic targets.

The florigen interactor BdES43 represses flowering in the model temperate grass Brachypodium distachyon

FLOWERING LOCUS T (FT) protein, physiologically florigen, has been identified as a system integrator of numerous flowering time pathways in many studies, and its homologs are found throughout the plant lineage. It is important to uncover how precisely florigenic homologs contribute to flowering initiation and how these factors interact genetically. Here we dissected the function of Brachypodium FT orthologs BdFTL1 and BdFTL2 using overexpression and gene-editing experiments. Transgenic assays showed that both BdFTL1 and BdFTL2 could promote flowering, whereas BdFTL2 was essential for flowering initiation. Notably, BdFTL1 is subject to alternative splicing (AS), and its transcriptional level and AS are significantly affected by BdFTL2. Additionally, BdFTL2 could bind with the PHD-containing protein BdES43, an H3K4me3 reader. Furthermore, BdES43 was antagonistic to BdFTL2 in flowering initiation in a transcription-dependent manner and significantly affected BdFTL1 expression. BdFTL2, BdES43 and H3K4me3 also had highly similar distribution patterns within the BdFTL1 locus, indicating their interplay in regulating target genes. Taken together, florigen BdFTL2 functions as a potential epigenetic effector of BdFTL1 by interacting with a BdES43-H3K4me3 complex. This finding provides an additional insight for the regulatory mechanism underlying the multifaceted roles of florigen.

Splicing Kinetics and Coordination Revealed by Direct Nascent RNA Sequencing through Nanopores

Understanding how splicing events are coordinated across numerous introns in metazoan RNA transcripts requires quantitative analyses of transient RNA processing events in living cells. We developed nanopore analysis of co-transcriptional processing (nano-COP), in which nascent RNAs are directly sequenced through nanopores, exposing the dynamics and patterns of RNA splicing without biases introduced by amplification. Long nano-COP reads reveal that, in human and Drosophila cells, splicing occurs after RNA polymerase II transcribes several kilobases of pre-mRNA, suggesting that metazoan splicing transpires distally from the transcription machinery. Inhibition of the branch-site recognition complex SF3B rapidly diminished global co-transcriptional splicing. We found that splicing order does not strictly follow the order of transcription and is associated with cis-acting elements, alternative splicing, and RNA-binding factors. Further, neighboring introns in human cells tend to be spliced concurrently, implying that splicing of these introns occurs cooperatively. Thus, nano-COP unveils the organizational complexity of RNA processing.

Exon Definition Facilitates Reliable Control of Alternative Splicing in the RON Proto-Oncogene

Alternative splicing is a key step in eukaryotic gene expression that allows for the production of multiple transcript and protein isoforms from the same gene. Even though splicing is perturbed in many diseases, we currently lack insights into regulatory mechanisms promoting its precision and efficiency. We analyze high-throughput mutagenesis data obtained for an alternatively spliced exon in the proto-oncogene RON and determine the functional units that control this splicing event. Using mathematical modeling of distinct splicing mechanisms, we show that alternative splicing is based in RON on a so-called "exon definition" mechanism. Here, the recognition of the adjacent exons by the spliceosome is required for removal of an intron. We use our model to analyze the differences between the exon and intron definition scenarios and find that exon definition prevents the accumulation of deleterious, partially spliced retention products during alternative splicing regulation. Furthermore, it modularizes splicing control, as multiple regulatory inputs are integrated into a common net input, irrespective of the location and nature of the corresponding cis-regulatory elements in the pre-messenger RNA. Our analysis suggests that exon definition promotes robust and reliable splicing outcomes in RON splicing.

The Function of Pre-mRNA Alternative Splicing in Mammal Spermatogenesis

Alternative pre-mRNA splicing plays important roles in co-transcriptional and post-transcriptional regulation of gene expression functioned during many developmental processes, such as spermatogenesis. The studies focusing on alternative splicing on spermatogenesis supported the notion that the development of testis is regulated by a higher level of alternative splicing than other tissues. Here, we aim to review the mechanisms underlying alternative splicing, particularly the splicing variants functioned in the process of spermatogenesis and the male infertility. There are five points regarding the alternative splicing including ⅰ) a brief introduction of alternative pre-mRNA splicing; ⅱ) the alternative splicing events in spermatogenesis-associated genes enriched in different stages of spermatogenesis; ⅲ) the mechanisms of alternative splicing regulation, such as splicing factors and m⁶A demethylation; ⅳ) the splice site recognition and alternative splicing, including the production and degradation of abnormal transcripts caused by gene variations and nonsense-mediated mRNA decay, respectively; ⅴ) abnormal alternative splicing correlated with male infertility. Taking together, this review highlights the impacts of alternative splicing and splicing variants in mammal spermatogenesis and provides new insights of the potential application of the alternative splicing into the therapy of male infertility.

Comparative Genomics Analysis of Ciliates Provides Insights on the Evolutionary History Within "Nassophorea-Synhymenia-Phyllopharyngea" Assemblage

Ciliated protists (ciliates) are widely used for investigating evolution, mostly due to their successful radiation after their early evolutionary branching. In this study, we employed high-throughput sequencing technology to reveal the phylogenetic position of Synhymenia, as well as two classes Nassophorea and Phyllopharyngea, which have been a long-standing puzzle in the field of ciliate systematics and evolution. We obtained genomic and transcriptomic data from single cells of one synhymenian (Chilodontopsis depressa) and six other species of phyllopharyngeans (Chilodochona sp., Dysteria derouxi, Hartmannula sinica, Trithigmostoma cucullulus, Trochilia petrani, and Trochilia sp.). Phylogenomic analysis based on 157 orthologous genes comprising 173,835 amino acid residues revealed the affiliation of C. depressa within the class Phyllopharyngea, and the monophyly of Nassophorea, which strongly support the assignment of Synhymenia as a subclass within the class Phyllopharyngea. Comparative genomic analyses further revealed that C. depressa shares more orthologous genes with the class Nassophorea than with Phyllopharyngea, and the stop codon usage in C. depressa resembles that of Phyllopharyngea. Functional enrichment analysis demonstrated that biological pathways in C. depressa are more similar to Phyllopharyngea than Nassophorea. These results suggest that genomic and transcriptomic data can be used to provide insights into the evolutionary relationships within the "Nassophorea-Synhymenia-Phyllopharyngea" assemblage.

Exon size and sequence conservation improves identification of splice-altering nucleotides

Pre-mRNA splicing is regulated through multiple trans-acting splicing factors. These regulators interact with the pre-mRNA at intronic and exonic positions. Given that most exons are protein coding, the evolution of exons must be modulated by a combination of selective coding and splicing pressures. It has previously been demonstrated that selective splicing pressures are more easily deconvoluted when phylogenetic comparisons are made for exons of identical size, suggesting that exon size-filtered sequence alignments may improve identification of nucleotides evolved to mediate efficient exon ligation. To test this hypothesis, an exon size database was created, filtering 76 vertebrate sequence alignments based on exon size conservation. In addition to other genomic parameters, such as splice-site strength, gene position, or flanking intron length, this database permits the identification of exons that are size- and/or sequence-conserved. Highly size-conserved exons are always sequence-conserved. However, sequence conservation does not necessitate exon size conservation. Our analysis identified evolutionarily young exons and demonstrated that length conservation is a strong predictor of alternative splicing. A published data set of approximately 5000 exonic SNPs associated with disease was analyzed to test the hypothesis that exon size-filtered sequence comparisons increase detection of splice-altering nucleotides. Improved splice predictions could be achieved when mutations occur at the third codon position, especially when a mutation decreases exon inclusion efficiency. The results demonstrate that coding pressures dominate nucleotide composition at invariable codon positions and that exon size-filtered sequence alignments permit identification of splice-altering nucleotides at wobble positions.

Genome-wide identification and characterization of the Hsp70 gene family in allopolyploid rapeseed (Brassica napus L.) compared with its diploid progenitors

Heat shock protein 70 (Hsp70) plays an essential role in plant growth and development, as well as stress response. Rapeseed (Brassica napus L.) originated from recently interspecific hybridization between Brassica rapa and Brassica oleracea. In this study, a total of 47 Hsp70 genes were identified in B. napus (A_nA_nC_nC_n genome), including 22 genes from A_n subgenome and 25 genes from C_n subgenome. Meanwhile, 29 and 20 Hsp70 genes were explored in B. rapa (A_rA_r genome) and B. oleracea (C_oC_o genome), respectively. Based on phylogenetic analysis, 114 Hsp70 proteins derived from B. napus, B. rapa, B. oleracea and Arabidopsis thaliana, were divided into 6 subfamilies containing 16 A_r-A_n and 11 C_o-C_n reliable orthologous pairs. The homology and synteny analysis indicated whole genome triplication and segmental duplication may be the major contributor for the expansion of Hsp70 gene family. Intron gain of BnHsp70 genes and domain loss of BnHsp70 proteins also were found in B. napus, associating with intron evolution and module evolution of proteins after allopolyploidization. In addition, transcriptional profiles analyses indicated that expression patterns of most BnHsp70 genes were tissue-specific. Moreover, Hsp70 orthologs exhibited different expression patterns in the same tissue and C_n subgenome biased expression was observed in leaf. These findings contribute to exploration of the evolutionary adaptation of polyploidy and will facilitate further application of BnHsp70 gene functions.

Deep Splicing Code: Classifying Alternative Splicing Events Using Deep Learning

Alternative splicing (AS) is the process of combining different parts of the pre-mRNA to produce diverse transcripts and eventually different protein products from a single gene. In computational biology field, researchers try to understand AS behavior and regulation using computational models known as “Splicing Codes”. The final goal of these algorithms is to make an in-silico prediction of AS outcome from genomic sequence. Here, we develop a deep learning approach, called Deep Splicing Code (DSC), for categorizing the well-studied classes of AS namely alternatively skipped exons, alternative 5’ss, alternative 3’ss, and constitutively spliced exons based only on the sequence of the exon junctions. The proposed approach significantly improves the prediction and the obtained results reveal that constitutive exons have distinguishable local characteristics from alternatively spliced exons. Using the motif visualization technique, we show that the trained models learned to search for competitive alternative splice sites as well as motifs of important splicing factors with high precision. Thus, the proposed approach greatly expands the opportunities to improve alternative splicing modeling. In addition, a web-server for AS events prediction has been developed based on the proposed method.

Mathematical modeling identifies potential gene structure determinants of co-transcriptional control of alternative pre-mRNA splicing

The spliceosome catalyzes the removal of introns from pre-messenger RNA (mRNA) and subsequent pairing of exons with remarkable fidelity. Some exons are known to be skipped or included in the mature mRNA in a cell type- or context-dependent manner (cassette exons), thereby contributing to the diversification of the human proteome. Interestingly, splicing is initiated (and sometimes completed) co-transcriptionally. Here, we develop a kinetic mathematical modeling framework to investigate alternative co-transcriptional splicing (CTS) and, specifically, the control of cassette exons' inclusion. We show that when splicing is co-transcriptional, default splice patterns of exon inclusion are more likely than when splicing is post-transcriptional, and that certain exons are more likely to be regulatable (i.e. cassette exons) than others, based on the exon-intron structure context. For such regulatable exons, transcriptional elongation rates may affect splicing outcomes. Within the CTS paradigm, we examine previously described hypotheses of co-operativity between splice sites of short introns (i.e. 'intron definition') or across short exons (i.e. 'exon definition'), and find that models encoding these faithfully recapitulate observations in the fly and human genomes, respectively.

Combinatorial regulation of alternative splicing

The generation of protein coding mRNAs from pre-mRNA is a fundamental biological process that is required for gene expression. Alternative pre-mRNA splicing is responsible for much of the transcriptomic and proteomic diversity observed in higher order eukaryotes. Aberrations that disrupt regular alternative splicing patterns are known to cause human diseases, including various cancers. Alternative splicing is a combinatorial process, meaning many factors affect which two splice sites are ligated together. The features that dictate exon inclusion are comprised of splice site strength, intron-exon architecture, RNA secondary structure, splicing regulatory elements, promoter use and transcription speed by RNA polymerase and the presence of post-transcriptional nucleotide modifications. A comprehensive view of all of the factors that influence alternative splicing decisions is necessary to predict splicing outcomes and to understand the molecular basis of disease. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Context matters: Regulation of splice donor usage

Elaborate research on splicing, starting in the late seventies, evolved from the discovery that 5' splice sites are recognized by their complementarity to U1 snRNA towards the realization that RNA duplex formation cannot be the sole basis for 5'ss selection. Rather, their recognition is highly influenced by a number of context factors including transcript architecture as well as splicing regulatory elements (SREs) in the splice site neighborhood. In particular, proximal binding of splicing regulatory proteins highly influences splicing outcome. The importance of SRE integrity especially becomes evident in the light of human pathogenic mutations where single nucleotide changes in SREs can severely affect the resulting transcripts. Bioinformatics tools nowadays greatly assist in the computational evaluation of 5'ss, their neighborhood and the impact of pathogenic mutations. Although predictions are already quite robust, computational evaluation of the splicing regulatory landscape still faces challenges to increase future reliability. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Identification and characterization of a new long noncoding RNA iab-1 in the Hox cluster of silkworm, Bombyx mori identification of iab-1

Long noncoding RNAs (lncRNAs) areinvolvedin a variety of biological processes. In silkworm, numerous lncRNAs have been predicted through deep transcriptome sequencing, but no functional role has been experimentally validated yet. Here, we characterized a new lncRNA iab-1 that was mainly encoded by the intergenic region between Bmabd-A and Bmabd-B in the Homeobox (Hox) cluster of the silkworm, Bombyx mori. More than seven alternative splicing isoforms of lncRNA iab-1 were cloned, which were subgrouped into types 1 and 2 based on the location of the 3'-ends. The iab-1 was expressed at a low level, but the expression of iab-1 peaked at several specific development stages, including 3 to 4 days during the embryonic stage, stages before fourth molting, and the sixth hour after the fourth molting, and early stages during metamorphosis. It was highly expressed in the nervus and epidermis, especially the epidermis of the posterior abdomen at the fourth instar premolting stage. The relationship between iab-1 and nearby Hox genes was analyzed at different developmental stages. Iab-1 expression was highly associated with Bmabd-A as well as Bmabd-B in the embryonic and larval stages, while this association was decreased at the metamorphic stage; iab-1 expression was highly associated with BmUbx only in the embryonic stage. Downregulation of iab-1 expression by small interfering RNA led to the death of most of the treated individuals at the larval stage, suggesting that iab-1 transcript expression might be involved in certain relevant physiological processes. The expression of Bmabd-A and Bmabd-B did not change in iab-1 downregulated individuals, indicating that the relevance between the two genes and iab-1 was not induced by iab-1 transcript. Collectively, the results showed that the newly identified iab-1 may be involved in some physiological processes, and the interaction between iab-1 and Hox genes was also preliminarily analyzed.

Genome-wide identification of cancer-specific alternative splicing in circRNA

Circular RNA (circRNA) is a group of RNA families generated by RNA circularization, which was discovered ubiquitously across different cancers. However, the internal structure of circRNA is difficult to determine due to alternative splicing that occurs in its exons and introns. Furthermore, cancer-specific alternative splicing of circRNA is less likely to be identified. Here, we proposed a de novo algorithm, CircSplice, that could identify internal alternative splicing in circRNA and compare differential circRNA splicing events between different conditions (http://gb.whu.edu.cn/CircSplice or https://github.com/GeneFeng/CircSplice). By applying CircSplice in clear cell renal cell carcinoma and bladder cancer, we detected 4498 and 2977 circRNA alternative splicing (circ-AS) events in the two datasets respectively and confirmed the expression of circ-AS events by RT-PCR. We further inspected the distributions and patterns of circ-AS in cancer and adjacent normal tissues. To further understand the potential functions of cancer-specific circ-AS, we classified those events into tumor suppressors and oncogenes and performed pathway enrichment analysis. This study is the first comprehensive view of cancer-specific circRNA alternative splicing, which could contribute significantly to regulation and functional research of circRNAs in cancers.

Spliceosome Profiling Visualizes Operations of a Dynamic RNP at Nucleotide Resolution

Tools to understand how the spliceosome functions in vivo have lagged behind advances in the structural biology of the spliceosome. Here, methods are described to globally profile spliceosome-bound pre-mRNA, intermediates, and spliced mRNA at nucleotide resolution. These tools are applied to three yeast species that span 600 million years of evolution. The sensitivity of the approach enables the detection of canonical and non-canonical events, including interrupted, recursive, and nested splicing. This application of statistical modeling uncovers independent roles for the size and position of the intron and the number of introns per transcript in substrate progression through the two catalytic stages. These include species-specific inputs suggestive of spliceosome-transcriptome coevolution. Further investigations reveal the ATP-dependent discard of numerous endogenous substrates after spliceosome assembly in vivo and connect this discard to intron retention, a form of splicing regulation. Spliceosome profiling is a quantitative, generalizable global technology used to investigate an RNP central to eukaryotic gene expression.

Predominant patterns of splicing evolution on human, chimpanzee and macaque evolutionary lineages

Although splicing is widespread and evolves rapidly among species, the mechanisms driving this evolution, as well as its functional implications, are not yet fully understood. We analyzed the evolution of splicing patterns based on transcriptome data from five tissues of humans, chimpanzees, rhesus macaques and mice. In total, 1526 exons and exon sets from 1236 genes showed significant splicing differences among primates. More than 60% of these differences represent constitutive-to-alternative exon transitions while an additional 25% represent changes in exon inclusion frequency. These two dominant evolutionary patterns have contrasting conservation, regulation and functional features. The sum of these features indicates that, despite their prevalence, constitutive-to-alternative exon transitions do not substantially contribute to long-term functional transcriptome changes. Conversely, changes in exon inclusion frequency appear to be functionally relevant, especially for changes taking place in the brain on the human evolutionary lineage.

Splicing of long non-coding RNAs primarily depends on polypyrimidine tract and 5' splice-site sequences due to weak interactions with SR proteins

Many nascent long non-coding RNAs (lncRNAs) undergo the same maturation steps as pre-mRNAs of protein-coding genes (PCGs), but they are often poorly spliced. To identify the underlying mechanisms for this phenomenon, we searched for putative splicing inhibitory sequences using the ncRNA-a2 as a model. Genome-wide analyses of intergenic lncRNAs (lincRNAs) revealed that lincRNA splicing efficiency positively correlates with 5'ss strength while no such correlation was identified for PCGs. In addition, efficiently spliced lincRNAs have higher thymidine content in the polypyrimidine tract (PPT) compared to efficiently spliced PCGs. Using model lincRNAs, we provide experimental evidence that strengthening the 5'ss and increasing the T content in PPT significantly enhances lincRNA splicing. We further showed that lincRNA exons contain less putative binding sites for SR proteins. To map binding of SR proteins to lincRNAs, we performed iCLIP with SRSF2, SRSF5 and SRSF6 and analyzed eCLIP data for SRSF1, SRSF7 and SRSF9. All examined SR proteins bind lincRNA exons to a much lower extent than expression-matched PCGs. We propose that lincRNAs lack the cooperative interaction network that enhances splicing, which renders their splicing outcome more dependent on the optimality of splice sites.

RNA Splicing and Disease: Animal Models to Therapies

Alternative splicing of pre-mRNA increases genetic diversity, and recent studies estimate that most human multiexon genes are alternatively spliced. If this process is not highly regulated and accurate, it leads to mis-splicing events, which may result in proteins with altered function. A growing body of work has implicated mis-splicing events in a range of diseases, including cancer, neurodegenerative diseases, and muscular dystrophies. Understanding the mechanisms that cause aberrant splicing events and how this leads to disease is vital for designing effective therapeutic strategies. In this review, we focus on advances in therapies targeting splicing, and highlight the animal models developed to recapitulate disease phenotypes as a model for testing these therapies.