Circular RNAs (circRNAs) are a relatively new class of RNA molecules, and knowledge about their biogenesis and function is still in its infancy. It was recently shown that alternative splicing underlies the formation of circular RNAs (circRNA) arising from the Titin (TTN) gene. Since the main mechanism by which circRNAs are formed is still unclear, we hypothesized that alternative splicing, and in particular exon skipping, is a major driver of circRNA production. We performed RNA sequencing on human and mouse hearts, mapped alternative splicing events, and overlaid these with expressed circRNAs at exon-level resolution. In addition, we performed RNA sequencing on hearts of Rbm20 KO mice to address how important Rbm20-mediated alternative splicing is in the production of cardiac circRNAs. In human and mouse hearts, we show that cardiac circRNAs are mostly (∼90%) produced from constitutive exons and less (∼10%) from alternatively spliced exons. In Rbm20 KO hearts, we identified 38 differentially expressed circRNAs of which 12 were produced from the Ttn gene. Even though Ttn appeared the most prominent target of Rbm20 for circularization, we also detected Rbm20-dependent circRNAs arising from other genes including Fan1, Stk39, Xdh, Bcl2l13, and Sorbs1 Interestingly, only Ttn circRNAs seemed to arise from Rbm20-mediated skipped exons. In conclusion, cardiac circRNAs are mostly derived from constitutive exons, suggesting that these circRNAs are generated at the expense of their linear counterpart and that circRNA production impacts the accumulation of the linear mRNA.

Alternative 3' and 5' splice site (ss) events constitute a significant part of all alternative splicing events. These events were also found to be related to several aberrant splicing diseases. However, only few of the characteristics that distinguish these events from alternative cassette exons are known currently. In this study, we compared the characteristics of constitutive exons, alternative cassette exons, and alternative 3'ss and 5'ss exons. The results revealed that alternative 3'ss and 5'ss exons are an intermediate state between constitutive and alternative cassette exons, where the constitutive side resembles constitutive exons, and the alternative side resembles alternative cassette exons. The results also show that alternative 3'ss and 5'ss exons exhibit low levels of symmetry (frame-preserving), similar to constitutive exons, whereas the sequence between the two alternative splice sites shows high symmetry levels, similar to alternative cassette exons. In addition, flanking intronic conservation analysis revealed that exons whose alternative splice sites are at least nine nucleotides apart show a high conservation level, indicating intronic participation in the regulation of their splicing, whereas exons whose alternative splice sites are fewer than nine nucleotides apart show a low conservation level. Further examination of these exons, spanning seven vertebrate species, suggests an evolutionary model in which the alternative state is a derivative of an ancestral constitutive exon, where a mutation inside the exon or along the flanking intron resulted in the creation of a new splice site that competes with the original one, leading to alternative splice site selection. This model was validated experimentally on four exons, showing that they indeed originated from constitutive exons that acquired a new competing splice site during evolution.

The fidelity of RNA splicing is regulated by a network of splicing enhancers and repressors, although the rules that govern this process are not yet fully understood. One mechanism that contributes to splicing fidelity is the repression of nonconserved cryptic exons by splicing factors that recognize dinucleotide repeats. We previously identified that TDP-43 and PTBP1/PTBP2 are capable of repressing cryptic exons utilizing UG and CU repeats, respectively. Here we demonstrate that hnRNP L (HNRNPL) also represses cryptic exons by utilizing exonic CA repeats, particularly near the 5'SS. We hypothesize that hnRNP L regulates CA repeat repression for both cryptic exon repression and developmental processes such as T cell differentiation.

The fidelity of RNA splicing is maintained by a network of factors, but the molecular mechanisms that govern this process have yet to be fully elucidated. We previously found that TDP-43, an RNA-binding protein implicated in neurodegenerative disease, utilizes UG microsatellites to repress nonconserved cryptic exons and prevent their incorporation into mRNA. Here, we report that two well-characterized splicing factors, polypyrimidine tract-binding protein 1 (PTBP1) and polypyrimidine tract-binding protein 2 (PTBP2), are also nonconserved cryptic exon repressors. In contrast to TDP-43, PTBP1 and PTBP2 utilize CU microsatellites to repress both conserved tissue-specific exons and nonconserved cryptic exons. Analysis of these conserved splicing events suggests that PTBP1 and PTBP2 repression is titrated to generate the transcriptome diversity required for neuronal differentiation. We establish that PTBP1 and PTBP2 are members of a family of cryptic exon repressors.

The split structure of most mammalian protein-coding genes allows for the potential to produce multiple different mRNA and protein isoforms from a single gene locus through the process of alternative splicing (AS). We propose a computational approach called UNCOVER based on a pair hidden Markov model to discover conserved coding exonic sequences subject to AS that have so far gone undetected. Applying UNCOVER to orthologous introns of known human and mouse genes predicts skipped exons or retained introns present in both species, while discriminating them from conserved noncoding sequences. The accuracy of the model is evaluated on a curated set of genes with known conserved AS events. The prediction of skipped exons in the approximately 1% of the human genome represented by the ENCODE regions leads to more than 50 new exon candidates. Five novel predicted AS exons were validated by RT-PCR and sequencing analysis of 15 introns with strong UNCOVER predictions and lacking EST evidence. These results imply that a considerable number of conserved exonic sequences and associated isoforms are still completely missing from the current annotation of known genes. UNCOVER also identifies a small number of candidates for conserved intron retention.

The spliceosomal factor TRAP150 is essential for pre-mRNA splicing in vivo and, when overexpressed, it enhances splicing efficiency. In this study, we found that TRAP150 interacted with the cleavage and polyadenylation specificity factor (CPSF) and co-fractionated with CPSF and RNA polymerase II. Moreover, TRAP150 preferentially associated with the U1 small ribonucleoprotein (snRNP). However, our data do not support a role for TRAP150 in alternative 5' splice site or exon selection or in alternative polyadenylation. Because U1 snRNP participates in premature cleavage and polyadenylation (PCPA), we tested whether TRAP150 is a cofactor in the control of PCPA. Although TRAP150 depletion had no significant effect on PCPA, overexpression of TRAP150 forced activation of a cryptic 3' splice site, yielding spliced PCPA transcripts. Mechanistic studies showed that TRAP150-activated splicing occurred in composite but not authentic terminal exons, and such an activity was enhanced by debilitation of U1 snRNP or interference with transcription elongation or termination. Together, these results indicate that TRAP150 provides an additional layer of PCPA regulation, through which it may increase the diversity of abortive RNA transcripts under conditions of compromised gene expression.

Alternative splicing accounts for a considerable portion of transcriptomic diversity, as most protein-coding genes are spliced into multiple mRNA isoforms. However, errors in splicing patterns can give rise to mis-splicing with pathological consequences, such as the congenital diseases familial dysautonomia, Duchenne muscular dystrophy, and spinal muscular atrophy. Small nuclear RNA (snRNA) components of the U snRNP family have been proposed as a therapeutic modality for the treatment of mis-splicing. U1 snRNAs offer great promise, with prior studies demonstrating in vivo efficacy, suggesting additional preclinical development is merited. Improvements in enabling technologies, including screening methodologies, gene delivery vectors, and relevant considerations from gene editing approaches justify further advancement of U1 snRNA as a therapeutic and research tool. With the goal of providing a user-friendly protocol, we compile and demonstrate a methodological toolkit for sequence-specific targeted perturbation of alternatively spliced pre-mRNA with engineered U1 snRNAs. We observe robust modulation of endogenous pre-mRNA transcripts targeted in two contrasting splicing contexts, SMN2 exon 7 and FAS exon 6, exhibiting the utility and applicability of engineered U1 snRNA to both inclusion and exclusion of targeted exons. We anticipate that these demonstrations will contribute to the usability of U1 snRNA in investigating splicing modulation in eukaryotic cells, increasing accessibility to the broader research community.

High-throughput splicing assays have demonstrated that many exonic variants can disrupt splicing; however, splice-disrupting variants distribute non-uniformly across genes. We propose the existence of exons that are particularly susceptible to splice-disrupting variants, which we refer to as hotspot exons. Hotspot exons are also more susceptible to splicing perturbation through drug treatment and knock-down of RNA-binding proteins. We develop a classifier for exonic splice-disrupting variants and use it to infer hotspot exons. We estimate that 1400 exons in the human genome are hotspots. Using panels of splicing reporters, we demonstrate how the ability of an exon to tolerate a mutation is inversely proportional to the strength of its neighboring splice sites.

Alternative splicing is an important mechanism to generate transcriptomic and phenotypic diversity. Existing methods have limited power to detect orthologous isoforms.

Alternative cassette exons are known to originate from two processes-exonization of intronic sequences and exon shuffling. Herein, we suggest an additional mechanism by which constitutively spliced exons become alternative cassette exons during evolution. We compiled a dataset of orthologous exons from human and mouse that are constitutively spliced in one species but alternatively spliced in the other. Examination of these exons suggests that the common ancestors were constitutively spliced. We show that relaxation of the 5' splice site during evolution is one of the molecular mechanisms by which exons shift from constitutive to alternative splicing. This shift is associated with the fixation of exonic splicing regulatory sequences (ESRs) that are essential for exon definition and control the inclusion level only after the transition to alternative splicing. The effect of each ESR on splicing and the combinatorial effects between two ESRs are conserved from fish to human. Our results uncover an evolutionary pathway that increases transcriptome diversity by shifting exons from constitutive to alternative splicing.

Core RNA-processing reactions in eukaryotic cells occur cotranscriptionally in a chromatin context, but the relationship between chromatin structure and pre-mRNA processing is poorly understood. We observed strong nucleosome depletion around human polyadenylation sites (PAS) and nucleosome enrichment just downstream of PAS. In genes with multiple alternative PAS, higher downstream nucleosome affinity was associated with higher PAS usage, independently of known PAS motifs that function at the RNA level. Conversely, exons were associated with distinct peaks in nucleosome density. Exons flanked by long introns or weak splice sites exhibited stronger nucleosome enrichment, and incorporation of nucleosome density data improved splicing simulation accuracy. Certain histone modifications, including H3K36me3 and H3K27me2, were specifically enriched on exons, suggesting active marking of exon locations at the chromatin level. Together, these findings provide evidence for extensive functional connections between chromatin structure and RNA processing.

Mammalian mRNA and lncRNA exons are often small compared to introns. The exon definition model predicts that exons splice autonomously, dependent on proximal exon sequence features, explaining their delineation within large introns. This model has not been examined on a genome-wide scale, however, leaving open the question of how often mRNA and lncRNA exons are autonomous. It is also unknown how frequently such exons can arise by chance. Here, we directly assayed large fragments (500-1000 bp) of the human genome by exon trapping, which detects exons spliced into a heterologous transgene, here designed with a large intron context. We define the trapped exons as "autonomous." We obtained ∼1.25 million trapped exons, including most known mRNA and well-annotated lncRNA internal exons, demonstrating that human exons are predominantly autonomous. mRNA exons are trapped with the highest efficiency. Nearly a million of the trapped exons are unannotated, most located in intergenic regions and antisense to mRNA, with depletion from the forward strand of introns. These exons are not conserved, suggesting they are nonfunctional and arose from random mutations. They are nonetheless highly enriched with known splicing promoting sequence features that delineate known exons. Novel autonomous exons are more numerous than annotated lncRNA exons, and computational models also indicate they will occur with similar frequency in any randomly generated sequence. These results show that most human coding exons splice autonomously, and provide an explanation for the existence of many unconserved lncRNAs, as well as a new annotation and inclusion levels of spliceable loci in the human genome.

Transcriptional isoforms are not just random combinations of exons. What has caused exons to be differentially spliced and whether exons with different splicing frequencies are subjected to divergent regulation by potential elements or splicing signals? Beyond the conventional classification for alternatively spliced exons (ASEs) and constitutively spliced exons (CSEs), we have classified exons from alternatively spliced human genes and their mouse orthologs (12,314 and 5,464, respectively) into four types based on their splicing frequencies. Analysis has indicated that different groups of exons presented divergent compositional and regulatory properties. Interestingly, with the decrease of splicing frequency, exons tend to have greater lengths, higher GC content, and contain more splicing elements and repetitive elements, which seem to imply that the splicing frequency is influenced by such factors. Comparison of non-alternatively spliced (NAS) mouse genes with alternatively spliced human orthologs also suggested that exons with lower splicing frequencies may be newly evolved ones which gained functions with splicing frequencies altered through the evolution. Our findings have revealed for the first time that certain factors may have critical influence on the splicing frequency, suggesting that exons with lower splicing frequencies may originate from old repetitive sequences, with splicing sites altered by mutation, gaining novel functions and become more frequently spliced.

Translation of nucleotides into a numeric form has been approached in many ways and has allowed researchers to investigate the properties of protein-coding sequences and noncoding sequences. Typically, more pronounced long-range correlations and increased regularity were found in intron-containing genes and in non-transcribed regulatory DNA sequences, compared to cDNA sequences or intron-less genes. The regularity is assessed by spectral tools defined on numerical translates. In most popular approaches of numerical translation the resulting spectra depend on the assignment of numerical values to nucleotides. Our contribution is to propose and illustrate a spectra which remains invariant to the translation rules used in traditional approaches.

COL4A5 is a causative gene of X-linked Alport syndrome (XLAS). Male patients with XLAS with nonsense variants have the most severe phenotypes of early onset end-stage kidney disease (ESKD); those with splicing variants have middle phenotypes and those with missense variants have the mildest phenotypes. Therefore, genotyping for male patients with XLAS can be used to predict kidney prognosis. Single-base substitutions at the last nucleotide position in each exon are known to affect splicing patterns and could be splicing variants. Nevertheless, in XLAS, these variants are generally considered to be missense variants, without conducting a transcript analysis, which underestimates some patients as having mild phenotypes. This study aimed to investigate whether single-base substitutions at the last nucleotide position of COL4A5 exons cause aberrant splicing.

Genetic analyses and systematic mutagenesis have revealed that synonymous, non-synonymous and intronic mutations frequently alter the inclusion levels of alternatively spliced exons, consistent with the concept that altered splicing might be a common mechanism by which mutations cause disease. However, most exons expressed in any cell are highly-included in mature mRNAs. Here, by performing deep mutagenesis of highly-included exons and by analysing the association between genome sequence variation and exon inclusion across the transcriptome, we report that mutations only very rarely alter the inclusion of highly-included exons. This is true for both exonic and intronic mutations as well as for perturbations in trans. Therefore, mutations that affect splicing are not evenly distributed across primary transcripts but are focussed in and around alternatively spliced exons with intermediate inclusion levels. These results provide a resource for prioritising synonymous and other variants as disease-causing mutations.

Alternative splicing is a major contributor to the diversity of eukaryotic transcriptomes and proteomes. Currently, large scale detection of alternative splicing using expressed sequence tags (ESTs) or microarrays does not capture all alternative splicing events. Moreover, for many species genomic data is being produced at a far greater rate than corresponding transcript data, hence in silico methods of predicting alternative splicing have to be improved.

Spatial restriction of mRNA to distinct subcellular locations enables local regulation and synthesis of proteins. However, the organizing principles of mRNA localization remain poorly understood. Here we analyzed subcellular transcriptomes of neural projections and soma of primary mouse cortical neurons and two neuronal cell lines and found that alternative last exons (ALEs) often confer isoform-specific localization. Surprisingly, gene-distal ALE isoforms were four times more often localized to neurites than gene-proximal isoforms. Localized isoforms were induced during neuronal differentiation and enriched for motifs associated with muscleblind-like (Mbnl) family RNA-binding proteins. Depletion of Mbnl1 and/or Mbnl2 reduced localization of hundreds of transcripts, implicating Mbnls in localization of mRNAs to neurites. We provide evidence supporting a model in which the linkage between genomic position of ALEs and subcellular localization enables coordinated induction of localization-competent mRNA isoforms through a post-transcriptional regulatory program that is induced during differentiation and reversed in cellular reprogramming and cancer.

Alu elements are major contributors to lineage-specific new exons in primate and human genomes. Recent studies indicate that some Alu exons have high transcript inclusion levels or tissue-specific splicing profiles, and may play important regulatory roles in modulating mRNA degradation or translational efficiency. However, the contribution of Alu exons to the human proteome remains unclear and controversial. The prevailing view is that exons derived from young repetitive elements, such as Alu elements, are restricted to regulatory functions and have not had adequate evolutionary time to be incorporated into stable, functional proteins.

Alternative splicing (AS) contributes significantly to protein diversity, by selectively using different combinations of exons of the same gene under certain circumstances. One particular type of AS is the use of alternative first exons (AFEs), which can have consequences far beyond the fine-tuning of protein functions. For example, AFEs may change the N-termini of proteins and thereby direct them to different cellular compartments. When alternative first exons are distant, they are usually associated with alternative promoters, thereby conferring an extra level of gene expression regulation. However, only few studies have examined the patterns of AFEs, and these analyses were mainly focused on mammalian genomes. Recent studies have shown that AFEs exist in the rice genome, and are regulated in a tissue-specific manner. Our current understanding of AFEs in plants is still limited, including important issues such as their regulation, contribution to protein diversity, and evolutionary conservation.