Differing patterns of selection in alternative and constitutive splice sites

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.

Improved modeling of RNA-binding protein motifs in an interpretable neural model of RNA splicing

Sequence-specific RNA-binding proteins (RBPs) play central roles in splicing decisions. Here, we describe a modular splicing architecture that leverages in vitro-derived RNA affinity models for 79 human RBPs and the annotated human genome to produce improved models of RBP binding and activity. Binding and activity are modeled by separate Motif and Aggregator components that can be mixed and matched, enforcing sparsity to improve interpretability. Training a new Adjusted Motif (AM) architecture on the splicing task not only yields better splicing predictions but also improves prediction of RBP-binding sites in vivo and of splicing activity, assessed using independent data.

mRNA splicing is modulated by intronic microRNAs

Splicing of transcripts is catalyzed by the spliceosome, a mega-complex consisting of hundreds of proteins and five snRNAs, which employs direct interactions. When U1 snRNA forms high-affinity binding, namely more than eight base pairs, with the 5'SS, the result is usually a suppressing effect on the splicing activity. This likely occurs due to the inefficient unwinding of U1/5'SS base-pairing or other regulatory obstructions. Here, we show in vitro and in patient-derived cell lines that pre-microRNAs can modulate the splicing reaction by interacting with U1 snRNA. This leads to reduced binding affinity to the 5'SS, and hence promotes the inclusion of exons containing 5'SS, despite sequence-based high affinity to U1. Application of the mechanism resulted in correction of the splicing defect in the disease-causing VCAN gene from an individual with Wagner syndrome. This pre-miRNA/U1 interaction can regulate the expression of alternatively spliced exons, thus extending the scope of mechanisms regulating splicing.

Long non-coding RNAs: definitions, functions, challenges and recommendations

Genes specifying long non-coding RNAs (lncRNAs) occupy a large fraction of the genomes of complex organisms. The term 'lncRNAs' encompasses RNA polymerase I (Pol I), Pol II and Pol III transcribed RNAs, and RNAs from processed introns. The various functions of lncRNAs and their many isoforms and interleaved relationships with other genes make lncRNA classification and annotation difficult. Most lncRNAs evolve more rapidly than protein-coding sequences, are cell type specific and regulate many aspects of cell differentiation and development and other physiological processes. Many lncRNAs associate with chromatin-modifying complexes, are transcribed from enhancers and nucleate phase separation of nuclear condensates and domains, indicating an intimate link between lncRNA expression and the spatial control of gene expression during development. lncRNAs also have important roles in the cytoplasm and beyond, including in the regulation of translation, metabolism and signalling. lncRNAs often have a modular structure and are rich in repeats, which are increasingly being shown to be relevant to their function. In this Consensus Statement, we address the definition and nomenclature of lncRNAs and their conservation, expression, phenotypic visibility, structure and functions. We also discuss research challenges and provide recommendations to advance the understanding of the roles of lncRNAs in development, cell biology and disease.

Synthetic polycistronic sequences in eukaryotes

The need for co-ordinate, high-level, and stable expression of multiple genes is essential for the engineering of biosynthetic circuits and metabolic pathways. This work outlines the functionality and design of IRES- and 2 A-peptide-based constructs by comparing different strategies for co-expression in polycistronic vectors. In particular, 2 A sequences are small peptides, mostly derived from viral polyproteins, that mediate a ribosome-skipping event such that several, different, separate proteins can be generated from a single open reading frame. When applied to metabolic engineering and synthetic gene circuits, 2 A peptides permit to achieve co-regulated and reliable expression of various genes in eukaryotic cells.

Predicting the effect of variants on splicing using Convolutional Neural Networks

Mutations that cause an error in the splicing of a messenger RNA (mRNA) can lead to diseases in humans. Various computational models have been developed to recognize the sequence pattern of the splice sites. In recent studies, Convolutional Neural Network (CNN) architectures were shown to outperform other existing models in predicting the splice sites. However, an insufficient effort has been put into extending the CNN model to predict the effect of the genomic variants on the splicing of mRNAs. This study proposes a framework to elaborate on the utility of CNNs to assess the effect of splice variants on the identification of potential disease-causing variants that disrupt the RNA splicing process. Five models, including three CNN-based and two non-CNN machine learning based, were trained and compared using two existing splice site datasets, Genome Wide Human splice sites (GWH) and a dataset provided at the Deep Learning and Artificial Intelligence winter school 2018 (DLAI). The donor sites were also used to test on the HSplice tool to evaluate the predictive models. To improve the effectiveness of predictive models, two datasets were combined. The CNN model with four convolutional layers showed the best splice site prediction performance with an AUPRC of 93.4% and 88.8% for donor and acceptor sites, respectively. The effects of variants on splicing were estimated by applying the best model on variant data from the ClinVar database. Based on the estimation, the framework could effectively differentiate pathogenic variants from the benign variants (p = 5.9 × 10⁻⁷). These promising results support that the proposed framework could be applied in future genetic studies to identify disease causing loci involving the splicing mechanism. The datasets and Python scripts used in this study are available on the GitHub repository at https://github.com/smiile8888/rna-splice-sites-recognition.

Gallop Racing Shifts Mature mRNA towards Introns: Does Exercise-Induced Stress Enhance Genome Plasticity?

Physical exercise is universally recognized as stressful. Among the "sport species", the horse is probably the most appropriate model for investigating the genomic response to stress due to the homogeneity of its genetic background. The aim of this work is to dissect the whole transcription modulation in Peripheral Blood Mononuclear Cells (PBMCs) after exercise with a time course framework focusing on unexplored regions related to introns and intergenic portions. PBMCs NGS from five 3 year old Sardinian Anglo-Arab racehorses collected at rest and after a 2000 m race was performed. Apart from differential gene expression ascertainment between the two time points the complexity of transcription for alternative transcripts was identified. Interestingly, we noted a transcription shift from the coding to the non-coding regions. We further investigated the possible causes of this phenomenon focusing on genomic repeats, using a differential expression approach and finding a strong general up-regulation of repetitive elements such as LINE. Since their modulation is also associated with the "exonization", the recruitment of repeats that act with regulatory functions, suggesting that there might be an active regulation of this transcriptional shift. Thanks to an innovative bioinformatic approach, our study could represent a model for the transcriptomic investigation of stress.

More than a messenger: Alternative splicing as a therapeutic target

Alternative splicing of pre-mRNA is an essential post- and co-transcriptional mechanism of gene expression regulation that produces multiple mature mRNA transcripts from a single gene. Genetic mutations that affect splicing underlie numerous devastating diseases. The complexity of splicing regulation allows for multiple therapeutic approaches to correct disease-associated mis-splicing events. In this review, we first highlight recent findings from therapeutic strategies that have used splice switching antisense oligonucleotides and small molecules that bind directly to RNA. Second, we summarize different genetic and chemical approaches to target components of the spliceosome to correct splicing defects in pathological conditions. Finally, we present an overview of compounds that target kinases and accessory pathways that intersect with the splicing machinery. Advancements in the understanding of disease-specific defects caused by mis-regulation of alternative splicing will certainly increase the development of therapeutic options for the clinic. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

The matrices and constraints of GT/AG splice sites of more than 1000 species/lineages

To provide a resource for the splice sites (SS) of different species, we calculated the matrices of nucleotide compositions of about 38 million splice sites from >1000 species/lineages. The matrices are enriched of aGGTAAGT (5'SS) or (Y)₆N(C/t)AG(g/a)t (3'SS) overall; however, they are quite diverse among hundreds of species. The diverse matrices remain prominent even under sequence selection pressures, suggesting the existence of diverse constraints as well as U snRNAs and other spliceosomal factors and/or their interactions with the splice sites. Using an algorithm to measure and compare the splice site constraints across all species, we demonstrate their distinct differences quantitatively. As an example of the resource's application to answering specific questions, we confirm that high constraints of particular positions are significantly associated with transcriptome-wide, increased occurrences of alternative splicing when uncommon nucleotides are present. More interestingly, the abundance of alternative splicing in 16 species correlates with the average constraint index of splice sites in a bell curve. This resource will allow users to assess specific sequences/splice sites against the consensus of every Ensembl-annotated species, and to explore the evolutionary changes or relationship to alternative splicing and transcriptome diversity. Web-search or update features are also included.

Changes in the expression of splicing factor transcripts and variations in alternative splicing are associated with lifespan in mice and humans

Dysregulation of splicing factor expression and altered alternative splicing are associated with aging in humans and other species, and also with replicative senescence in cultured cells. Here, we assess whether expression changes of key splicing regulator genes and consequent effects on alternative splicing are also associated with strain longevity in old and young mice, across 6 different mouse strains with varying lifespan (A/J, NOD.B10Sn-H2(b) /J, PWD.Phj, 129S1/SvlmJ, C57BL/6J and WSB/EiJ). Splicing factor expression and changes to alternative splicing were associated with strain lifespan in spleen and to a lesser extent in muscle. These changes mainly involved hnRNP splicing inhibitor transcripts with most changes more marked in spleens of young animals from long-lived strains. Changes in spleen isoform expression were suggestive of reduced cellular senescence and retained cellular proliferative capacity in long-lived strains. Changes in muscle isoform expression were consistent with reduced pro-inflammatory signalling in longer-lived strains. Two splicing regulators, HNRNPA1 and HNRNPA2B1, were also associated with parental longevity in humans, in the InCHIANTI aging study. Splicing factors may represent a driver, mediator or early marker of lifespan in mouse, as expression differences were present in the young animals of long-lived strains. Changes to alternative splicing patterns of key senescence genes in spleen and key remodelling genes in muscle suggest that correct regulation of alternative splicing may enhance lifespan in mice. Expression of some splicing factors in humans was also associated with parental longevity, suggesting that splicing regulation may also influence lifespan in humans.

RNA-Seq Analysis of Differential Splice Junction Usage and Intron Retentions by DEXSeq

Alternative splicing is an important biological process in the generation of multiple functional transcripts from the same genomic sequences. Differential analysis of splice junctions (SJs) and intron retentions (IRs) is helpful in the detection of alternative splicing events. In this study, we conducted differential analysis of SJs and IRs by use of DEXSeq, a Bioconductor package originally designed for differential exon usage analysis in RNA-seq data analysis. We set up an analysis pipeline including mapping of RNA-seq reads, the preparation of count tables of SJs and IRs as the input files, and the differential analysis in DEXSeq. We analyzed the public RNA-seq datasets generated from RNAi experiments on Drosophila melanogaster S2-DRSC cells to deplete RNA-binding proteins (GSE18508). The analysis confirmed previous findings on the alternative splicing of the trol and Ant2 (sesB) genes in the CG8144 (ps)-depletion experiment and identified some new alternative splicing events in other RNAi experiments. We also identified IRs that were confirmed in our SJ analysis. The proposed method used in our study can output the genomic coordinates of differentially used SJs and thus enable sequence motif search. Sequence motif search and gene function annotation analysis helped us infer the underlying mechanism in alternative splicing events. To further evaluate this method, we also applied the method to public RNA-seq data from human breast cancer (GSE45419) and the plant Arabidopsis (SRP008262). In conclusion, our study showed that DEXSeq can be adapted to differential analysis of SJs and IRs, which will facilitate the identification of alternative splicing events and provide insights into the molecular mechanisms of transcription processes and disease development.

The exon junction complex is required for definition and excision of neighboring introns in Drosophila

The exon junction complex (EJC) plays a central role in controlling RNA fate and aids faithful splicing of pre-mRNAs containing large introns via an unknown mechanism. Brennecke and colleagues show that the core EJC plus the accessory factors RnpS1 and Acinus aid in the definition and efficient splicing of neighboring introns. Interestingly, the most highly affected intron belongs to the piwi locus, which explains the reported transposon desilencing in EJC-depleted Drosophila ovaries. Based on transcriptome-wide analysis, the authors propose that the dependency of splicing on the EJC is exploited to control the temporal order of splicing events.

Splicing of pre-mRNAs results in the deposition of the exon junction complex (EJC) upstream of exon–exon boundaries. The EJC plays crucial post-splicing roles in export, translation, localization, and nonsense-mediated decay of mRNAs. It also aids faithful splicing of pre-mRNAs containing large introns, albeit via an unknown mechanism. Here, we show that the core EJC plus the accessory factors RnpS1 and Acinus aid in definition and efficient splicing of neighboring introns. This requires prior deposition of the EJC in close proximity to either an upstream or downstream splicing event. If present in isolation, EJC-dependent introns are splicing-defective also in wild-type cells. Interestingly, the most affected intron belongs to the piwi locus, which explains the reported transposon desilencing in EJC-depleted Drosophila ovaries. Based on a transcriptome-wide analysis, we propose that the dependency of splicing on the EJC is exploited as a means to control the temporal order of splicing events.

NFAT1 and NFAT3 cooperate with HDAC4 during regulation of alternative splicing of PMCA isoforms in PC12 cells

Background

The bulk of human genes undergo alternative splicing (AS) upon response to physiological stimuli. AS is a great source of protein diversity and biological processes and is associated with the development of many diseases. Pheochromocytoma is a neuroendocrine tumor, characterized by an excessive Ca²⁺-dependent secretion of catecholamines. This underlines the importance of balanced control of calcium transport via regulation of gene expression pattern, including different calcium transport systems, such as plasma membrane Ca²⁺-ATPases (PMCAs), abundantly expressed in pheochromocytoma chromaffin cells (PC12 cells). PMCAs are encoded by four genes (Atp2b1, Atp2b2, Atp2b3, Atp2b4), whose transcript products undergo alternative splicing giving almost 30 variants.

Results

In this scientific report, we propose a novel mechanism of regulation of PMCA alternative splicing in PC12 cells through cooperation of the nuclear factor of activated T-cells (NFAT) and histone deacetylases (HDACs). Luciferase assays showed increased activity of NFAT in PC12 cells, which was associated with altered expression of PMCA. RT-PCR experiments suggested that inhibition of the transcriptional activity of NFAT might result in the rearrangement of PMCA splicing variants in PC12 cells. NFAT inhibition led to dominant expression of 2x/c, 3x/a and 4x/a PMCA variants, while in untreated cells the 2w,z/b, 3z,x/b,c,e,f, and 4x/b variants were found as well. Furthermore, chromatin immunoprecipitation experiments showed that NFAT1-HDAC4 or NFAT3-HDAC4 complexes might be involved in regulation of PMCA2x splicing variant generation.

Conclusions

We suggest that the influence of NFAT/HDAC on PMCA isoform composition might be important for altered dopamine secretion by PC12 cells.

Weak negative and positive selection and the drift load at splice sites

Splice sites (SSs) are short sequences that are crucial for proper mRNA splicing in eukaryotic cells, and therefore can be expected to be shaped by strong selection. Nevertheless, in mammals and in other intron-rich organisms, many of the SSs often involve nonconsensus (Nc), rather than consensus (Cn), nucleotides, and beyond the two critical nucleotides, the SSs are not perfectly conserved between species. Here, we compare the SS sequences between primates, and between Drosophila fruit flies, to reveal the pattern of selection acting at SSs. Cn-to-Nc substitutions are less frequent, and Nc-to-Cn substitutions are more frequent, than neutrally expected, indicating, respectively, negative and positive selection. This selection is relatively weak (1 < |4Nes| < 4), and has a similar efficiency in primates and in Drosophila. Within some nucleotide positions, the positive selection in favor of Nc-to-Cn substitutions is weaker than the negative selection maintaining already established Cn nucleotides; this difference is due to site-specific negative selection favoring current Nc nucleotides. In general, however, the strength of negative selection protecting the Cn alleles is similar in magnitude to the strength of positive selection favoring replacement of Nc alleles, as expected under the simple nearly neutral turnover. In summary, although a fraction of the Nc nucleotides within SSs is maintained by selection, the abundance of deleterious nucleotides in this class suggests a substantial genome-wide drift load.

Function of alternative splicing

Almost all polymerase II transcripts undergo alternative pre-mRNA splicing. Here, we review the functions of alternative splicing events that have been experimentally determined. The overall function of alternative splicing is to increase the diversity of mRNAs expressed from the genome. Alternative splicing changes proteins encoded by mRNAs, which has profound functional effects. Experimental analysis of these protein isoforms showed that alternative splicing regulates binding between proteins, between proteins and nucleic acids as well as between proteins and membranes. Alternative splicing regulates the localization of proteins, their enzymatic properties and their interaction with ligands. In most cases, changes caused by individual splicing isoforms are small. However, cells typically coordinate numerous changes in 'splicing programs', which can have strong effects on cell proliferation, cell survival and properties of the nervous system. Due to its widespread usage and molecular versatility, alternative splicing emerges as a central element in gene regulation that interferes with almost every biological function analyzed.

DNA hypermethylation of alternatively spliced and repeat sequences in humans

DNA methylation is presently accepted as a tentative regulatory parameter in splicing. Recently, we reported significant methylation differences among various exonic splicing-enhancing elements and alternative splicing events, based on CpG methylation data from the Human Epigenome Project for chromosomes 6, 20 and 22. Presently, using a different computational approach and the same database, we report: (a) significant increase of hypermethylation in intronic and exonic sequences close to acceptor sites, relative to overall introns and exons, respectively (1,973 CpGs examined); (b) frequent CpGs, mostly hypomethylated, in donors and infrequent CpGs mostly hypermethylated, in acceptors; and (c) hypermethylation in cassette exons which are occasionally spliced and have weaker average splicing potential, relative to constitutive exons (p < 0.0001). CpGs are hypomethylated in non-coding exons (only 16 % hypermethylation). Single-exon genes, similarly to first exons, frequently contain hypomethylated CpGs, while in internal and last exons CpGs are more frequently hypermethylated. Methylation is also more frequent in strange introns and splice sites processed by the minor spliceosome, e.g., ATAC, (p < 0.0001 in all cases), but not in sites of incomplete processing, e.g., retained introns or bleeding exons, (p = 0.706 and p = 0.313, respectively). Most Alus, which are known to contribute to transcript presentation, are heavily methylated, in contrast with other Alus, e.g., AluJo and mammalian interspersed repetitive elements which have been previously associated with alternative expression. These results elucidate the role of intragenic methylation in association with alternative splicing and facilitate the evaluation of genomic variations/polymorphisms and the development of tools for the prediction of alternative splicing events.

Bioinformatic analysis of the human DHRS4 gene cluster and a proposed mechanism for its transcriptional regulation

BACKGROUND

The human DHRS4 gene cluster consists of three genes, DHRS4, DHRS4L2 and DHRS4L1. Among them, DHRS4 encodes NADP(H)-dependent retinol dehydrogenase/reductase. In a previous study, we investigated the alternative splicing of DHRS4 and DHRS4L2. DHRS4L1 was added to the gene cluster recently, but little is known about its structure and expression. To reveal the regulatory mechanism of the DHRS4 gene cluster expression, we studied the structure and transcription of DHRS4L1 in the context of the transcriptional behaviors of the human DHRS4 gene cluster. Based on the results of bioinformatics analysis, we propose a possible mechanism for the transcriptional regulation of the human DHRS4 gene cluster.

RESULTS

The homologous comparison analysis suggests that DHRS4, DHRS4L2 and DHRS4L1 are three homologous genes in human. DHRS4L1 and DHRS4L2 are paralogues of DHRS4, and DHRS4L2 is the most recent member of the DHRS4 gene cluster. In the minus strand of the human DHRS4 gene cluster, a gene transcribed in an antisense direction was found containing a 5' sequence overlapping the region of exon 1 and promoter of DHRS4. By cloning the full length of RNA variants through 5'RACE and 3'RACE, we identified two transcription start sites, within exon a2 and exon 1, of this newly named gene DHRS4L1 using neuroblastoma cell line BE(2)-M17. Analysis of exon composition in the transcripts of DHRS4 gene cluster revealed that exon 1 was absent in all the transcripts initiated from exon a1 of DHRS4L2 and exon a2 of DHRS4L1.

CONCLUSIONS

Alternatively spliced RNA variants are prevalent in the human DHRS4 gene cluster. Based on the analysis of gene transcripts and bioinformatic prediction, we propose here that antisense transcription may be involved in the transcriptional initiation regulation of DHRS4 and in the posttranscriptional splicing of DHRS4L2 and DRHS4L1 for the homologous identity of DHRS4 gene cluster. Beside the alternative transcriptional start sites, the antisense RNA is novel possible factor serving to remove exon 1 from the transcripts initiated from exon a1 and exon a2.

Splice site strength-dependent activity and genetic buffering by poly-G runs

Pre-mRNA splicing is regulated through combinatorial activity of RNA motifs including splice sites and splicing regulatory elements (SREs). Here, we show that the activity of the G-run class of SREs is ∼4-fold higher when adjacent to intermediate strength 5'ss relative to weak 5'ss, and ∼1.3-fold higher relative to strong 5'ss. This dependence on 5'ss strength was observed in splicing reporters and in global microarray and mRNA-Seq analyses of splicing changes following RNAi against heterogeneous nuclear ribonucleoprotein (hnRNP) H, which crosslinked to G-runs adjacent to many regulated exons. An exon’s responsiveness to changes in hnRNP H levels therefore depends in a complex way on G-run abundance and 5'ss strength, and other splicing factors may function similarly. This pattern of activity enables G-runs and hnRNP H to buffer the effects of 5'ss mutations, augmenting the frequency of 5'ss polymorphism and the evolution of new splicing patterns.

Complex selection on 5' splice sites in intron-rich organisms

In contrast to the typically streamlined genomes of prokaryotes, many eukaryotic genomes are riddled with long intergenic regions, spliceosomal introns, and repetitive elements. What explains the persistence of these and other seemingly suboptimal structures? There are three general hypotheses: (1) the structures in question are not actually suboptimal but optimal, being favored by selection, for unknown reasons; (2) the structures are not suboptimal, but of (essentially) equal fitness to "optimal" ones; or (3) the structures are truly suboptimal, but selection is too weak to systematically eliminate them. The 5' splice sites of introns offer a rare opportunity to directly test these hypotheses. Intron-poor species show a clear consensus splice site; most introns begin with the same six nucleotide sequence (typically GTAAGT or GTATGT), indicating efficient selection for this consensus sequence. In contrast, intron-rich species have much less pronounced boundary consensus sequences, and only small minorities of introns in intron-rich species share the same boundary sequence. We studied rates of evolutionary change of 5' splice sites in three groups of closely related intron-rich species--three primates, five Drosophila species, and four Cryptococcus fungi. Surprisingly, the results indicate that changes from consensus-to-variant nucleotides are generally disfavored by selection, but that changes from variant to consensus are neither favored nor disfavored. This evolutionary pattern is consistent with selective differences across introns, for instance, due to compensatory changes at other sites within the gene, which compensate for the otherwise suboptimal consensus-to-variant changes in splice boundaries.

Alternative-splicing-based bicistronic vectors for ratio-controlled protein expression and application to recombinant antibody production

In the last decade polycistronic vectors have become essential tools for both basic science and gene therapy applications. In order to co-express heterologous polypeptides, different systems have been developed from Internal Ribosome Entry Site (IRES) based vectors to the use of the 2A peptide. Unfortunately, these methods are not fully suitable for the efficient and reproducible modulation of the ratio between the proteins of interest. Here we describe a novel bicistronic vector type based on the use of alternative splicing. By modifying the consensus sequence that governs splicing, we demonstrate that the ratio between the synthesized proteins could easily vary from 1 : 10 to 10 : 1. We have established this system with luciferase genes and we extended its application to the production of recombinant monoclonal antibodies. We have shown that these vectors could be used in several typical cell lines with similar efficiencies. We also present an adaptation of these vectors to hybrid alternative splicing/IRES constructs that allow a ratio-controlled expression of proteins of interest in stably transfected cell lines.

Comparative component analysis of exons with different splicing frequencies

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.

Where do introns come from?

The systems for mRNA surveillance, capping, and cleavage/polyadenylation are proposed to play pivotal roles in the physical establishment and distribution of spliceosomal introns along a transcript.

Improved identification of conserved cassette exons using Bayesian networks

Background

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.

Results

Here, we show that the use of Bayesian networks (BNs) allows accurate prediction of evolutionary conserved exon skipping events. At a stringent false positive rate of 0.5%, our BN achieves an improved true positive rate of 61%, compared to a previously reported 50% on the same dataset using support vector machines (SVMs). Incorporating several novel discriminative features such as intronic splicing regulatory elements leads to the improvement. Features related to mRNA secondary structure increase the prediction performance, corroborating previous findings that secondary structures are important for exon recognition. Random labelling tests rule out overfitting. Cross-validation on another dataset confirms the increased performance. When using the same dataset and the same set of features, the BN matches the performance of an SVM in earlier literature. Remarkably, we could show that about half of the exons which are labelled constitutive but receive a high probability of being alternative by the BN, are in fact alternative exons according to the latest EST data. Finally, we predict exon skipping without using conservation-based features, and achieve a true positive rate of 29% at a false positive rate of 0.5%.

Conclusion

BNs can be used to achieve accurate identification of alternative exons and provide clues about possible dependencies between relevant features. The near-identical performance of the BN and SVM when using the same features shows that good classification depends more on features than on the choice of classifier. Conservation based features continue to be the most informative, and hence distinguishing alternative exons from constitutive ones without using conservation based features remains a challenging problem.

Cross-kingdom patterns of alternative splicing and splice recognition

BACKGROUND

Variations in transcript splicing can reveal how eukaryotes recognize intronic splice sites. Retained introns (RIs) commonly appear when the intron definition (ID) mechanism of splice site recognition inconsistently identifies intron-exon boundaries, and cassette exons (CEs) are often caused by variable recognition of splice junctions by the exon definition (ED) mechanism. We have performed a comprehensive survey of alternative splicing across 42 eukaryotes to gain insight into how spliceosomal introns are recognized.

RESULTS

All eukaryotes we studied exhibit RIs, which appear more frequently than previously thought. CEs are also present in all kingdoms and most of the organisms in our analysis. We observe that the ratio of CEs to RIs varies substantially among kingdoms, while the ratio of competing 3' acceptor and competing 5' donor sites remains nearly constant. In addition, we find the ratio of CEs to RIs in each organism correlates with the length of its introns. In all 14 fungi we examined, as well as in most of the 9 protists, RIs far outnumber CEs. This differs from the trend seen in 13 multicellular animals, where CEs occur much more frequently than RIs. The six plants we analyzed exhibit intermediate proportions of CEs and RIs.

CONCLUSION

Our results suggest that most extant eukaryotes are capable of recognizing splice sites via both ID and ED, although ED is most common in multicellular animals and ID predominates in fungi and most protists.

Positive selection acting on splicing motifs reflects compensatory evolution

We have used comparative genomics to characterize the evolutionary behavior of predicted splicing regulatory motifs. Using base substitution rates in intronic regions as a calibrator for neutral change, we found a strong avoidance of synonymous substitutions that disrupt predicted exonic splicing enhancers or create predicted exonic splicing silencers. These results attest to the functionality of the hexameric motif set used and suggest that they are subject to purifying selection. We also found that synonymous substitutions in constitutive exons tend to create exonic splicing enhancers and to disrupt exonic splicing silencers, implying positive selection for these splicing promoting events. We present evidence that this positive selection is the result of splicing-positive events compensating for splicing-negative events as well as for mutations that weaken splice-site sequences. Such compensatory events include nonsynonymous mutations, synonymous mutations, and mutations at splice sites. Compensation was also seen from the fact that orthologous exons tend to maintain the same number of predicted splicing motifs. Our data fit a splicing compensation model of exon evolution, in which selection for splicing-positive mutations takes place to counter the effect of an ongoing splicing-negative mutational process, with the exon as a whole being conserved as a unit of splicing. In the course of this analysis, we observed that synonymous positions in general are conserved relative to intronic sequences, suggesting that messenger RNA molecules are rich in sequence information for functions beyond protein coding and splicing.

The impact of multiple splice sites in human L1 elements

LINE-1 elements represent a significant proportion of mammalian genomes. The impact of their activity on the structure and function of the host genomes has been recognized from the time of their discovery as an endogenous source of insertional mutagenesis. L1 elements contain numerous functional internal polyadenylation signals and splice sites that generate a variety of processed L1 transcripts. These sites are also reported to contribute to the generation of hybrid transcripts between L1 elements and host genes. Using northern blot analysis we demonstrate that L1 splicing, but not L1 polyadenylation, is delayed during the course of L1 expression. L1 splicing can also be negatively regulated by EBV SM protein known to alter this process. These results suggest a potential for L1 mRNA processing to be regulated in a tissue- and/or development-specific manner. The delay in L1 splicing may also serve to protect host genes from the excessive burden of L1 interference with their normal expression via aberrant splicing.

Evidence for a trade-off between translational efficiency and splicing regulation in determining synonymous codon usage in Drosophila melanogaster

In Drosophila melanogaster, synonymous codons corresponding to the most abundant cognate tRNAs are used more frequently, especially in highly expressed genes. Increased use of such "optimal" codons is considered an adaptation for translational efficiency. Need it always be the case that selection should favor the use of a translationally optimal codon? Here, we investigate one possible confounding factor, namely, the need to specify information in exons necessary to enable correct splicing. As expected from such a model, in Drosophila many codons show different usage near intron-exon boundaries versus exon core regions. However, this finding is in principle also consistent with Hill-Robertson effects modulating usage of translationally optimal codons. However, several results support the splice model over the translational selection model: 1) the trends in codon usage are strikingly similar to those in mammals in which codon usage near boundaries correlates with abundance in exonic splice enhancers (ESEs), 2) codons preferred near boundaries tend to be enriched for A and avoid C (conversely those avoided near boundaries prefer C rather than A), as expected were ESEs involved, and 3) codons preferred near boundaries are typically not translationally optimal. We conclude that usage of translationally optimal codons usage is compromised in the vicinity of splice junctions in intron-containing genes, to the effect that we observe higher levels of usage of translationally optimal codons at the center of exons. On the gene level, however, controlling for known correlates of codon bias, the impact on codon usage patterns is quantitatively small. These results have implications for inferring aspects of the mechanism of splicing given nothing more than a well-annotated genome.