Computational comparative analyses of alternative splicing regulation using full-length cDNA of various eukaryotes

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.

POINT technology illuminates the processing of polymerase-associated intact nascent transcripts

Mammalian chromatin is the site of both RNA polymerase II (Pol II) transcription and coupled RNA processing. However, molecular details of such co-transcriptional mechanisms remain obscure, partly because of technical limitations in purifying authentic nascent transcripts. We present a new approach to characterize nascent RNA, called polymerase intact nascent transcript (POINT) technology. This three-pronged methodology maps nascent RNA 5′ ends (POINT-5), establishes the kinetics of co-transcriptional splicing patterns (POINT-nano), and profiles whole transcription units (POINT-seq). In particular, we show by depletion of the nuclear exonuclease Xrn2 that this activity acts selectively on cleaved 5′ P-RNA at polyadenylation sites. Furthermore, POINT-nano reveals that co-transcriptional splicing either occurs immediately after splice site transcription or is delayed until Pol II transcribes downstream sequences. Finally, we connect RNA cleavage and splicing with either premature or full-length transcript termination. We anticipate that POINT technology will afford full dissection of the complexity of co-transcriptional RNA processing.

•

POINT methodology dissects intact nascent RNA processing

•

Specificity of Xrn2 exonuclease in co-transcriptional RNA degradation

•

Splicing suppresses Xrn2-dependent premature termination

•

Different kinetic classes of co-transcriptional splicing in human genes

Dynamic changes in intron retention are tightly associated with regulation of splicing factors and proliferative activity during B-cell development

Intron retention (IR) has been proposed to modulate the delay between transcription and translation. Here, we provide an exhaustive characterization of IR in differentiated white blood cells from both the myeloid and lymphoid lineage where we observed highest levels of IR in monocytes and B-cells, in addition to previously reported granulocytes. During B-cell differentiation, we found an increase in IR from the bone marrow precursors to cells residing in secondary lymphoid organs. B-cells that undergo affinity maturation to become antibody producing plasma cells steadily decrease retention. In general, we found an inverse relationship between global IR levels and both the proliferative state of cells, and the global levels of expression of splicing factors. IR dynamics during B-cell differentiation appear to be conserved between human and mouse, suggesting that IR plays an important biological role, evolutionary conserved, during blood cell differentiation. By correlating the expression of non-core splicing factors with global IR levels, and analyzing RNA binding protein knockdown and eCLIP data, we identify a few splicing factors likely playing an evolutionary conserved role in IR regulation. Our work provides new insights into the role of IR during hematopoiesis, and on the main factors involved in regulating IR.

RegSNPs-intron: a computational framework for predicting pathogenic impact of intronic single nucleotide variants

Single nucleotide variants (SNVs) in intronic regions have yet to be systematically investigated for their disease-causing potential. Using known pathogenic and neutral intronic SNVs (iSNVs) as training data, we develop the RegSNPs-intron algorithm based on a random forest classifier that integrates RNA splicing, protein structure, and evolutionary conservation features. RegSNPs-intron showed excellent performance in evaluating the pathogenic impacts of iSNVs. Using a high-throughput functional reporter assay called ASSET-seq (ASsay for Splicing using ExonTrap and sequencing), we evaluate the impact of RegSNPs-intron predictions on splicing outcome. Together, RegSNPs-intron and ASSET-seq enable effective prioritization of iSNVs for disease pathogenesis.

Activation of cryptic splice sites in three patients with chronic granulomatous disease

Background

Chronic granulomatous disease (CGD) is a primary immune deficiency caused by mutations in the genes encoding the structural components of the phagocyte NADPH oxidase. As a result, the patients cannot generate sufficient amounts of reactive oxygen species required for killing pathogenic microorganisms.

Methods

We analyzed NADPH oxidase activity and component expression in neutrophils, performed genomic DNA and cDNA analysis, and used mRNA splicing prediction tools to evaluate the impact of mutations.

Results

In two patients with CGD, we had previously found mutations that cause aberrant pre‐mRNA splicing. In one patient an exonic mutation in a cryptic donor splice site caused the deletion of the 3' part of exon 6 from the mRNA of CYBB. This patient suffers from X‐linked CGD. The second patient, with autosomal CGD, has a mutation in the donor splice site of intron 1 of CYBA that activates a cryptic donor splice site downstream in intron 1, causing the insertion of intronic sequences in the mRNA. The third patient, recently analyzed, also with autosomal CGD, has a mutation in intron 4 of CYBA, 15 bp from the acceptor splice site. This mutation weakens a branch site and activates a cryptic acceptor splice site, causing the insertion of 14 intronic nucleotides into the mRNA.

Conclusion

We found three different mutations, one exonic, one in a donor splice site and one intronic, that all caused missplicing of pre‐mRNA. We analyzed these mutations with four different splice prediction programs and found that predictions of splice site strength, splice enhancer and splice silencer protein binding and branch site strength are all essential for correct prediction of pre‐mRNA splicing.

Mutation in an exonic splicing enhancer site causing chronic granulomatous disease

In a male patient suffering from X-linked chronic granulomatous disease (CGD) we found a c.389G>T mutation in exon 5 of the CYBB gene. We have analyzed why 95% of the transcripts of this gene lacked exon 5, leading to a frameshift and premature termination codon. The mutation was located in a region comprising three putative exonic splicing enhancer binding sites, for SRSF1, SRFS2 and SRFS6, according to the ESEfinder Tool (http://rulai.cshl.edu/cgi-bin/tools/ESE3/esefinder.cgi). With the Analyser Splice Tool we calculated the probability of skipping of exon 5 in CYBB mRNA, and by means of Sroogle the number of putative binding motifs for splicing enhancer and splicing silencer proteins (http://astlab.tau.ac.il/index.php). These analyses clarify why this exon was skipped in the majority of the mRNA. The normally spliced transcript contains an amino acid change p.Arg130Leu. This poorly expressed transcript gives rise to a protein with low expression but presumably normal activity, leading to a respiratory burst activity in the patient's neutrophils of about 15% of normal.

regSNPs-splicing: a tool for prioritizing synonymous single-nucleotide substitution

While synonymous single-nucleotide variants (sSNVs) have largely been unstudied, since they do not alter protein sequence, mounting evidence suggests that they may affect RNA conformation, splicing, and the stability of nascent-mRNAs to promote various diseases. Accurately prioritizing deleterious sSNVs from a pool of neutral ones can significantly improve our ability of selecting functional genetic variants identified from various genome-sequencing projects, and, therefore, advance our understanding of disease etiology. In this study, we develop a computational algorithm to prioritize sSNVs based on their impact on mRNA splicing and protein function. In addition to genomic features that potentially affect splicing regulation, our proposed algorithm also includes dozens structural features that characterize the functions of alternatively spliced exons on protein function. Our systematical evaluation on thousands of sSNVs suggests that several structural features, including intrinsic disorder protein scores, solvent accessible surface areas, protein secondary structures, and known and predicted protein family domains, show significant differences between disease-causing and neutral sSNVs. Our result suggests that the protein structure features offer an added dimension of information while distinguishing disease-causing and neutral synonymous variants. The inclusion of structural features increases the predictive accuracy for functional sSNV prioritization.

Regulation of the chloroplastic copper chaperone (CCS) and cuprozinc superoxide dismutase (CSD2) by alternative splicing and copper excess in Glycine max

Metal homeostasis is an important aspect of plant physiology, and the copper transport into the chloroplast and its fate after delivery is of special relevance for plants. In this work, the regulation of the chloroplastic copper chaperone for the cuprozinc superoxide dismutase (GmCCS) and its target, the cuprozinc superoxide dismutase (GmCSD2), was investigated in photosynthetic cell suspensions and entire plants from Glycine max (L.) Merr. Both genes were expressed in cell suspensions and in all plant tissues analysed, and their RNAs matured by alternative splicing with intron retention (IntronR). This mechanism generated a spliced and three non-spliced mRNAs in the case of GmCCS but only a spliced and a non-spliced mRNAs in GmCSD2. Copper excess strongly upregulated the expression of both fully spliced mRNAs but mostly unaffected the non-spliced forms. In entire plants, some tissue specificity was also observed depending on copper content status. At the protein level, the GmCCS was mostly unaffected but the GmCSD2 was strongly induced under copper excess in all subcellular fractions analysed, suggesting a post-transcriptional regulation for the former. This different protein regulation of the chaperone and its target may indicate some additional function for the CSD2 protein. In addition to its well-known superoxide dismutase (SOD) activity, it may also function as a metal sink in copper excess availability to avoid metal cell damage. Furthermore, the GmCCS seems to be present in the stroma only but the GmCSD2 was present in both stroma and thylakoids despite the general idea that the SOD enzymes are typically soluble stroma proteins. The presence of the SOD enzyme on the surface of the thylakoid membranes is reasonable considering that the superoxide radical (O2-) is preferentially formed at the acceptor side of the PSI.

Genome-wide analysis reveals SR protein cooperation and competition in regulated splicing

SR proteins are well-characterized RNA binding proteins that promote exon inclusion by binding to exonic splicing enhancers (ESEs). However, it has been unclear whether regulatory rules deduced on model genes apply generally to activities of SR proteins in the cell. Here, we report global analyses of two prototypical SR proteins, SRSF1 (SF2/ASF) and SRSF2 (SC35), using splicing-sensitive arrays and CLIP-seq on mouse embryo fibroblasts (MEFs). Unexpectedly, we find that these SR proteins promote both inclusion and skipping of exons in vivo, but their binding patterns do not explain such opposite responses. Further analyses reveal that loss of one SR protein is accompanied by coordinated loss or compensatory gain in the interaction of other SR proteins at the affected exons. Therefore, specific effects on regulated splicing by one SR protein actually depend on a complex set of relationships with multiple other SR proteins in mammalian genomes.

U2AF65 adapts to diverse pre-mRNA splice sites through conformational selection of specific and promiscuous RNA recognition motifs

Degenerate splice site sequences mark the intron boundaries of pre-mRNA transcripts in multicellular eukaryotes. The essential pre-mRNA splicing factor U2AF⁶⁵ is faced with the paradoxical tasks of accurately targeting polypyrimidine (Py) tracts preceding 3′ splice sites while adapting to both cytidine and uridine nucleotides with nearly equivalent frequencies. To understand how U2AF⁶⁵ recognizes degenerate Py tracts, we determined six crystal structures of human U2AF⁶⁵ bound to cytidine-containing Py tracts. As deoxy-ribose backbones were required for co-crystallization with these Py tracts, we also determined two baseline structures of U2AF⁶⁵ bound to the deoxy-uridine counterparts and compared the original, RNA-bound structure. Local structural changes suggest that the N-terminal RNA recognition motif 1 (RRM1) is more promiscuous for cytosine-containing Py tracts than the C-terminal RRM2. These structural differences between the RRMs were reinforced by the specificities of wild-type and site-directed mutant U2AF⁶⁵ for region-dependent cytosine- and uracil-containing RNA sites. Small-angle X-ray scattering analyses further demonstrated that Py tract variations select distinct inter-RRM spacings from a pre-existing ensemble of U2AF⁶⁵ conformations. Our results highlight both local and global conformational selection as a means for universal 3′ splice site recognition by U2AF⁶⁵.

Widespread recognition of 5' splice sites by noncanonical base-pairing to U1 snRNA involving bulged nucleotides

An established paradigm in pre-mRNA splicing is the recognition of the 5' splice site (5'ss) by canonical base-pairing to the 5' end of U1 small nuclear RNA (snRNA). We recently reported that a small subset of 5'ss base-pair to U1 in an alternate register that is shifted by 1 nucleotide. Using genetic suppression experiments in human cells, we now demonstrate that many other 5'ss are recognized via noncanonical base-pairing registers involving bulged nucleotides on either the 5'ss or U1 RNA strand, which we term "bulge registers." By combining experimental evidence with transcriptome-wide free-energy calculations of 5'ss/U1 base-pairing, we estimate that 10,248 5'ss (∼5% of human 5'ss) in 6577 genes use bulge registers. Several of these 5'ss occur in genes with mutations causing genetic diseases and are often associated with alternative splicing. These results call for a redefinition of an essential element for gene expression that incorporates these registers, with important implications for the molecular classification of splicing mutations and for alternative splicing.

Regulation of alternative splicing by the core spliceosomal machinery

Alternative splicing (AS) plays a major role in the generation of proteomic diversity and in gene regulation. However, the role of the basal splicing machinery in regulating AS remains poorly understood. Here we show that the core snRNP (small nuclear ribonucleoprotein) protein SmB/B' self-regulates its expression by promoting the inclusion of a highly conserved alternative exon in its own pre-mRNA that targets the spliced transcript for nonsense-mediated mRNA decay (NMD). Depletion of SmB/B' in human cells results in reduced levels of snRNPs and a striking reduction in the inclusion levels of hundreds of additional alternative exons, with comparatively few effects on constitutive exon splicing levels. The affected alternative exons are enriched in genes encoding RNA processing and other RNA-binding factors, and a subset of these exons also regulate gene expression by activating NMD. Our results thus demonstrate a role for the core spliceosomal machinery in controlling an exon network that appears to modulate the levels of many RNA processing factors.

COMIT: identification of noncoding motifs under selection in coding sequences

COMIT is presented; an algorithm for detecting functional non-coding motifs in coding regions, separating nucleotide and amino acid effects.

Coding nucleotide sequences contain myriad functions independent of their encoded protein sequences. We present the COMIT algorithm to detect functional noncoding motifs in coding regions using sequence conservation, explicitly separating nucleotide from amino acid effects. COMIT concurs with diverse experimental datasets, including splicing enhancers, silencers, replication motifs, and microRNA targets, and predicts many novel functional motifs. Intriguingly, COMIT scores are well-correlated to scores uncalibrated for amino acids, suggesting that nucleotide motifs often override peptide-level constraints.

Large introns in relation to alternative splicing and gene evolution: a case study of Drosophila bruno-3

Background

Alternative splicing (AS) of maturing mRNA can generate structurally and functionally distinct transcripts from the same gene. Recent bioinformatic analyses of available genome databases inferred a positive correlation between intron length and AS. To study the interplay between intron length and AS empirically and in more detail, we analyzed the diversity of alternatively spliced transcripts (ASTs) in the Drosophila RNA-binding Bruno-3 (Bru-3) gene. This gene was known to encode thirteen exons separated by introns of diverse sizes, ranging from 71 to 41,973 nucleotides in D. melanogaster. Although Bru-3's structure is expected to be conducive to AS, only two ASTs of this gene were previously described.

Results

Cloning of RT-PCR products of the entire ORF from four species representing three diverged Drosophila lineages provided an evolutionary perspective, high sensitivity, and long-range contiguity of splice choices currently unattainable by high-throughput methods. Consequently, we identified three new exons, a new exon fragment and thirty-three previously unknown ASTs of Bru-3. All exon-skipping events in the gene were mapped to the exons surrounded by introns of at least 800 nucleotides, whereas exons split by introns of less than 250 nucleotides were always spliced contiguously in mRNA. Cases of exon loss and creation during Bru-3 evolution in Drosophila were also localized within large introns. Notably, we identified a true de novo exon gain: exon 8 was created along the lineage of the obscura group from intronic sequence between cryptic splice sites conserved among all Drosophila species surveyed. Exon 8 was included in mature mRNA by the species representing all the major branches of the obscura group. To our knowledge, the origin of exon 8 is the first documented case of exonization of intronic sequence outside vertebrates.

Conclusion

We found that large introns can promote AS via exon-skipping and exon turnover during evolution likely due to frequent errors in their removal from maturing mRNA. Large introns could be a reservoir of genetic diversity, because they have a greater number of mutable sites than short introns. Taken together, gene structure can constrain and/or promote gene evolution.

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.

Exon creation and establishment in human genes

BACKGROUND

A large proportion of species-specific exons are alternatively spliced. In primates, Alu elements play a crucial role in the process of exon creation but many new exons have appeared through other mechanisms. Despite many recent studies, it is still unclear which are the splicing regulatory requirements for de novo exonization and how splicing regulation changes throughout an exon's lifespan.

RESULTS

Using comparative genomics, we have defined sets of exons with different evolutionary ages. Younger exons have weaker splice-sites and lower absolute values for the relative abundance of putative splicing regulators between exonic and adjacent intronic regions, indicating a less consolidated splicing regulation. This relative abundance is shown to increase with exon age, leading to higher exon inclusion. We show that this local difference in the density of regulators might be of biological significance, as it outperforms other measures in real exon versus pseudo-exon classification. We apply this new measure to the specific case of the exonization of anti-sense Alu elements and show that they are characterized by a general lack of exonic splicing silencers.

CONCLUSIONS

Our results suggest that specific sequence environments are required for exonization and that these can change with time. We propose a model of exon creation and establishment in human genes, in which splicing decisions depend on the relative local abundance of regulatory motifs. Using this model, we provide further explanation as to why Alu elements serve as a major substrate for exon creation in primates. Finally, we discuss the benefits of integrating such information in gene prediction.

Similar selective factors affect both between-gene and between-exon divergence in Drosophila

As a consequence of alternative splicing, a gene's exons will have different frequencies of inclusion into mature mRNA and different patterns of expression. These differences affect their patterns of evolutionary divergence. Using the recently reannotated genome of Drosophila melanogaster and the genome sequences of four closely related species of the melanogaster subgroup, we investigated the effect of alternative splicing, inclusion level (defined as the number of transcripts an exon is found in), and expression pattern on exon evolution across divergence times ranging from 1 to 12.5 Ma. Genes undergoing alternative splicing have a broader pattern of expression associated with a lower divergence rate in comparison with genes with a single annotated protein isoform. Within genes undergoing alternative splicing, we report a significant effect of inclusion level on exon evolution, as alternatively spliced exons are less conserved than constitutively spliced exons. More generally, there are significant negative correlations between inclusion level and exon evolutionary rates that can be associated with relaxation of selection. A significant effect of expression pattern on evolution rates is also observed. Overall, we found that similar selective factors such as the expression level and the pattern of expression are affecting both gene and exon evolution.

In silico analysis of phosphoproteome data suggests a rich-get-richer process of phosphosite accumulation over evolution

Recent phosphoproteome analyses using mass spectrometry-based technologies have provided new insights into the extensive presence of protein phosphorylation in various species and have raised the interesting question of how this protein modification was gained evolutionarily on such a large scale. We investigated this issue by using human and mouse phosphoproteome data. We initially found that phosphoproteins followed a power-law distribution with regard to their number of phosphosites: most of the proteins included only a few phosphosites, but some included dozens of phosphosites. The power-law distribution, unlike more commonly observed distributions such as normal and log-normal distributions, is considered by the field of complex systems science to be produced by a specific rich-get-richer process called preferential attachment growth. Therefore, we explored the factors that may have promoted the rich-get-richer process during phosphosite evolution. We conducted a bioinformatics analysis to evaluate the relationship of amino acid sequences of phosphoproteins with the positions of phosphosites and found an overconcentration of phosphosites in specific regions of protein surfaces and implications that in many phosphoproteins these clusters of phosphosites are activated simultaneously. Multiple phosphosites concentrated in limited spaces on phosphoprotein surfaces may therefore function biologically as cooperative modules that are resistant to selective pressures during phosphoprotein evolution. We therefore proposed a hypothetical model by which the modularization of multiple phosphosites has been resistant to natural selection and has driven the rich-get-richer process of the evolutionary growth of phosphosite numbers.

Fast splice site detection using information content and feature reduction

BACKGROUND

Accurate identification of splice sites in DNA sequences plays a key role in the prediction of gene structure in eukaryotes. Already many computational methods have been proposed for the detection of splice sites and some of them showed high prediction accuracy. However, most of these methods are limited in terms of their long computation time when applied to whole genome sequence data.

RESULTS

In this paper we propose a hybrid algorithm which combines several effective and informative input features with the state of the art support vector machine (SVM). To obtain the input features we employ information content method based on Shannon's information theory, Shapiro's score scheme, and Markovian probabilities. We also use a feature elimination scheme to reduce the less informative features from the input data.

CONCLUSION

In this study we propose a new feature based splice site detection method that shows improved acceptor and donor splice site detection in DNA sequences when the performance is compared with various state of the art and well known methods.

Prediction of alternatively skipped exons and splicing enhancers from exon junction arrays

BACKGROUND

Alternative splicing of exons in a pre-mRNA transcript is an important mechanism which contributes to protein diversity in human. Arrays for detecting alternative splicing are available using several different probe designs, including those based on exon-junctions. In this work, we introduce a new method for predicting alternatively skipped exons from exon-junction arrays. Predictions based on our method are compared against controls and their sequences are analyzed to identify motifs important for regulating alternative splicing.

RESULTS

Our comparison of several alternative methods shows that an exon-skipping score based on neighboring junctions best discriminates between positive and negative controls. Sequence analysis of our predicted exons confirms the presence of known splicing regulatory sequences. In addition, we also derive a set of development-related alternatively spliced genes based on fetal versus adult tissue comparisons and find that our predictions are consistent with their functional annotations. Ab initio motif finding algorithms are applied to identify several motifs that may be relevant for splicing during development.

CONCLUSION

This work describes a new method for analyzing exon-junction arrays, identifies sequence motifs that are specific for alternative and constitutive splicing and suggests a role for several known splicing factors and their motifs in developmental regulation.

The effect of intron length on exon creation ratios during the evolution of mammalian genomes

Recent studies report that alternatively spliced exons tend to occur in longer introns, which is attributed to the length constraints for splice site pairing for the two major splicing mechanisms, intron definition versus exon definition. Using genome-wide studies of EST and microarray data from human and mouse, we have analyzed the distribution of various subsets of alternatively spliced exons, based on their inclusion level and evolutionary history, versus increasing intron length. Alternative exons may be included in either a major or minor fraction of all transcripts (known as major-form and minor-form exons, respectively). We find that major-form exons are seven- to eightfold more likely to be contained in short introns (<400 nt) than minor-form exons, which occur preferentially in longer introns. Since minor-form exons are more likely to be novel (approximately 75%), this implied that novel exons arise more frequently in longer introns. To test this hypothesis, we used whole genome alignments to classify exons according to their phylogenetic age. We find that older exons, i.e., exons that are conserved in all mammals, predominate at shorter intron lengths, for both major- and minor-form exons. In contrast, exons that arose recently during primate evolution are more prevalent at longer intron lengths (>1000 nt). This suggests that the observed correlation of longer intron lengths with alternatively spliced exons may be at least partly due to biases in the probability of exon creation, which is higher in long introns.

Solution conformation and thermodynamic characteristics of RNA binding by the splicing factor U2AF65

The U2 auxiliary factor large subunit (U2AF65) is an essential pre-mRNA splicing factor for the initial stages of spliceosome assembly. Tandem RNA recognition motifs (RRM)s of U2AF65 recognize polypyrimidine tract signals adjacent to 3' splice sites. Despite the central importance of U2AF65 for splice site recognition, the relative arrangement of the U2AF65 RRMs and the energetic forces driving polypyrimidine tract recognition remain unknown. Here, the solution conformation of the U2AF65 RNA binding domain determined using small angle x-ray scattering reveals a bilobal shape without apparent interdomain contacts. The proximity of the N and C termini within the inter-RRM configuration is sufficient to explain the action of U2AF65 on spliceosome components located both 5' and 3' to its binding site. Isothermal titration calorimetry further demonstrates that an unusually large enthalpy-entropy compensation underlies U2AF65 recognition of an optimal polyuridine tract. Qualitative similarities were observed between the pairwise distance distribution functions of the U2AF65 RNA binding domain and those either previously observed for N-terminal RRMs of Py tract-binding protein that lack interdomain contacts or calculated from the high resolution coordinates of a U2AF65 deletion variant bound to RNA. To further test this model, the shapes and RNA interactions of the wild-type U2AF65 RNA binding domain were compared with those of U2AF65 variants containing either Py tract-binding protein linker sequences or a deletion within the inter-RRM linker. Results of these studies suggest inter-RRM conformational plasticity as a possible means for U2AF65 to universally identify diverse pre-mRNA splice sites.

A DNA damage signal activates and derepresses exon inclusion in Drosophila TAF1 alternative splicing

Signal-dependent alternative splicing is important for regulating gene expression in eukaryotes, yet our understanding of how signals impact splicing mechanisms is limited. A model to address this issue is alternative splicing of Drosophila TAF1 pre-mRNA in response to camptothecin (CPT)-induced DNA damage signals. CPT treatment of Drosophila S2 cells causes increased inclusion of TAF1 alternative cassette exons 12a and 13a through an ATR signaling pathway. To evaluate the role of TAF1 pre-mRNA sequences in the alternative splicing mechanism, we developed a TAF1 minigene (miniTAF1) and an S2 cell splicing assay that recapitulated key aspects of CPT-induced alternative splicing of endogenous TAF1. Analysis of miniTAF1 indicated that splice site strength underlies independent and distinct mechanisms that control exon 12a and 13a inclusion. Mutation of the exon 13a weak 5' splice site or weak 3' splice site to a consensus sequence was sufficient for constitutive exon 13a inclusion. In contrast, mutation of the exon 12a strong 5' splice site or moderate 3' splice site to a consensus sequence was only sufficient for constitutive exon 12a inclusion in the presence of CPT-induced signals. Analogous studies of the exon 13 3' splice site suggest that exon 12a inclusion involves signal-dependent pairing between constitutive and alternative splice sites. Finally, intronic elements identified by evolutionary conservation were necessary for full repression of exon 12a inclusion or full activation of exon 13a inclusion and may be targets of CPT-induced signals. In summary, this work defines the role of sequence elements in the regulation of TAF1 alternative splicing in response to a DNA damage signal.

Searching for splicing motifs

Intron removal during pre-mRNA splicing in higher eukaryotes requires the accurate identification of the two splice sites at the ends of the exons, or exon definition. The sequences constituting the splice sites provide insufficient information to distinguish true splice sites from the greater number of false splice sites that populate transcripts. Additional information used for exon recognition resides in a large number of positively or negatively acting elements that lie both within exons and in the adjacent introns. The identification of such sequence motifs has progressed rapidly in recent years, such that extensive lists are now available for exonic splicing enhancers and exonic splicing silencers. These motifs have been identified both by empirical experiments and by computational predictions, the validity of the latter being confirmed by experimental verification. Molecular searches have been carried out either by the selection of sequences that bind to splicing factors, or enhance or silence splicing in vitro or in vivo. Computational methods have focused on sequences of 6 or 8 nucleotides that are over- or under-represented in exons, compared to introns or transcripts that do not undergo splicing. These various methods have sought to provide global definitions of motifs, yet the motifs are distinctive to the method used for identification and display little overlap. Astonishingly, at least three-quarters of a typical mRNA would be comprised of these motifs. A present challenge lies in understanding how the cell integrates this surfeit of information to generate what is usually a binary splicing decision.

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.

Alternative approach to a heavy weight problem

Obesity is reaching epidemic proportions in developed countries and represents a significant risk factor for hypertension, heart disease, diabetes, and dyslipidemia. Splicing mutations constitute at least 14% of disease-causing mutations, thus implicating polymorphisms that affect splicing as likely candidates for disease susceptibility. A recent study suggested that genes associated with obesity were significantly enriched for rare nucleotide variants. Here, we examined these variants and revealed that they are located near splice junctions and tend to affect exonic splicing regulatory sequences. We also show that the majority of the exons that harbor these SNPs are constitutively spliced, yet they exhibit weak splice sites, typical to alternatively spliced exons, and are hence suboptimal for recognition by the splicing machinery and prone to become alternatively spliced. Using ex vivo assays, we tested a few representative variants and show that they indeed affect splicing by causing a shift from a constitutive to an alternative pattern, suggesting a possible link between extreme body mass index and abnormal splicing patterns.

Case study for identification of potentially indel-caused alternative expression isoforms in the rice subspecies japonica and indica by integrative genome analysis

Alternative splicing (AS) is one of the most significant components of the functional complexity of the eukaryote genome, increasing protein diversity, creating isoforms, and affecting mRNA stability. Recently, whole genome sequences and large microarray data sets have become available, making data integration feasible and allowing the study of the possible regulatory mechanism of AS in rice (Oryza sativa) by erecting and testing hypotheses before doing bench studies. We have developed a new strategy and have identified 215 rice genes with alternative expression isoforms related to insertion and deletion (indel) between subspecies indica and subspecies japonica. We did a case study for alternative expression isoforms of the rice peroxidase gene LOC_Os06g48030 to investigate possible mechanisms by which indels caused alternative splicing between the indica and the japonica varieties by mining of array data together with validation by RT-PCR and genome sequencing analysis. Multiple poly(A) signals were detected in the specific indel region for LOC_Os06g48030. We present a new methodology to promote more discoveries of potentially indel-caused AS genes in rice, which may serve as the foundation for research into the regulatory mechanism of alternative expression isoforms between subspecies.

Alternative splicing of pre-mRNA in plants

Alternative splicing is an important cellular mechanism that increases the diversity of gene products. The study of alternatively spliced genes reported so far in plants is far less documented than that in mammals, but considerable results have been reported, showing the role of these genes in regulating mechanisms, influencing factors, and specificities and function of alternative splicing in plants. This review summarizes briefly the major progress made on alternative splicing in plants.

Relating alternative splicing to proteome complexity and genome evolution

Prior to genomics, studies of alternative splicing primarily focused on the function and mechanism of alternative splicing in individual genes and exons. This has changed dramatically since the late 1990s. High-throughput genomics technologies, such as EST sequencing and microarrays designed to detect changes in splicing, led to genome-wide discoveries and quantification of alternative splicing in a wide range of species from human to Arabidopsis. Consensus estimates of AS frequency in the human genome grew from less than 5% in mid-1990s to as high as 60-74% now. The rapid growth in sequence and microarray data for alternative splicing has made it possible to look into the global impact of alternative splicing on protein function and evolution of genomes. In this chapter, we review recent research on alternative splicing's impact on proteomic complexity and its role in genome evolution.

Evolutionary divergence of exon flanks: a dissection of mutability and selection

The intronic sequences flanking exon-intron junctions (i.e., exon flanks) are important for splice site recognition and pre-mRNA splicing. Recent studies show a higher degree of sequence conservation at flanks of alternative exons, compared to flanks of constitutive exons. In this article we performed a detailed analysis on the evolutionary divergence of exon flanks between human and chimpanzee, aiming to dissect the impact of mutability and selection on their evolution. Inside exon flanks, sites that might reside in ancestral CpG dinucleotides evolved significantly faster than sites outside of ancestral CpG dinucleotides. This result reflects a systematic variation of mutation rates (mutability) at exon flanks, depending on the local CpG contexts. Remarkably, we observed a significant reduction of the nucleotide substitution rate in flanks of alternatively spliced exons, independent of the site-by-site variation in mutability due to different CpG contexts. Our data provide concrete evidence for increased purifying selection at exon flanks associated with regulation of alternative splicing.

Differentiated evolutionary rates in alternative exons and the implications for splicing regulation

BACKGROUND

Alternatively spliced exons play an important role in the diversification of gene function in most metazoans and are highly regulated by conserved motifs in exons and introns. Two contradicting properties have been associated to evolutionary conserved alternative exons: higher sequence conservation and higher rate of non-synonymous substitutions, relative to constitutive exons. In order to clarify this issue, we have performed an analysis of the evolution of alternative and constitutive exons, using a large set of protein coding exons conserved between human and mouse and taking into account the conservation of the transcript exonic structure. Further, we have also defined a measure of the variation of the arrangement of exonic splicing enhancers (ESE-conservation score) to study the evolution of splicing regulatory sequences. We have used this measure to correlate the changes in the arrangement of ESEs with the divergence of exon and intron sequences.

RESULTS

We find evidence for a relation between the lack of conservation of the exonic structure and the weakening of the sequence evolutionary constraints in alternative and constitutive exons. Exons in transcripts with non-conserved exonic structures have higher synonymous (dS) and non-synonymous (dN) substitution rates than exons in conserved structures. Moreover, alternative exons in transcripts with non-conserved exonic structure are the least constrained in sequence evolution, and at high EST-inclusion levels they are found to be very similar to constitutive exons, whereas alternative exons in transcripts with conserved exonic structure have a dS significantly lower than average at all EST-inclusion levels. We also find higher conservation in the arrangement of ESEs in constitutive exons compared to alternative ones. Additionally, the sequence conservation at flanking introns remains constant for constitutive exons at all ESE-conservation values, but increases for alternative exons at high ESE-conservation values.

CONCLUSION

We conclude that most of the differences in dN observed between alternative and constitutive exons can be explained by the conservation of the transcript exonic structure. Low dS values are more characteristic of alternative exons with conserved exonic structure, but not of those with non-conserved exonic structure. Additionally, constitutive exons are characterized by a higher conservation in the arrangement of ESEs, and alternative exons with an ESE-conservation similar to that of constitutive exons are characterized by a conservation of the flanking intron sequences higher than average, indicating the presence of more intronic regulatory signals.

Alternative splicing and RNA selection pressure--evolutionary consequences for eukaryotic genomes

Genome-wide analyses of alternative splicing have established its nearly ubiquitous role in gene regulation in many organisms. Genome sequencing and comparative genomics have made it possible to look in detail at the evolutionary history of specific alternative exons or splice sites, resulting in a flurry of publications in recent years. Here, we consider how alternative splicing has contributed to the evolution of modern genomes, and discuss constraints on evolution associated with alternative splicing that might have important medical implications.

Protein modularity of alternatively spliced exons is associated with tissue-specific regulation of alternative splicing

Recent comparative genomic analysis of alternative splicing has shown that protein modularity is an important criterion for functional alternative splicing events. Exons that are alternatively spliced in multiple organisms are much more likely to be an exact multiple of 3 nt in length, representing a class of “modular” exons that can be inserted or removed from the transcripts without affecting the rest of the protein. To understand the precise roles of these modular exons, in this paper we have analyzed microarray data for 3,126 alternatively spliced exons across ten mouse tissues generated by Pan and coworkers. We show that modular exons are strongly associated with tissue-specific regulation of alternative splicing. Exons that are alternatively spliced at uniformly high transcript inclusion levels or uniformly low levels show no preference for protein modularity. In contrast, alternatively spliced exons with dramatic changes of inclusion levels across mouse tissues (referred to as “tissue-switched” exons) are both strikingly biased to be modular and are strongly conserved between human and mouse. The analysis of different subsets of tissue-switched exons shows that the increased protein modularity cannot be explained by the overall exon inclusion level, but is specifically associated with tissue-switched alternative splicing.

Alternative splicing is a biological process that generates multiple mRNA and protein variants through alternative combinations of protein-coding exons. It is a widespread mechanism of gene regulation in higher eukaryotes. In recent years, scientists have found that when an exon is observed to be alternatively spliced in multiple species, its length is much more likely to be an exact multiple of three nucleotides. Since each amino acid is encoded by three nucleotides, these exons can be inserted or removed from the transcript as a “modular” protein-coding unit, without affecting the downstream protein translation. However, the precise roles of these modular exons in gene regulation and genome evolution remain unclear.

Xing and Lee have now investigated these modular exons using high-throughput genomics data. They analyzed the mouse splicing microarray data from the research group of Dr. Benjamin Blencowe at University of Toronto. Exons whose alternative splicing levels vary dramatically across multiple tissues are much more likely to be modular exons and are highly conserved during human and mouse evolution. This study establishes a strong link between protein modularity of alternatively spliced exons and tissue-specific regulation of alternative splicing. It provides new insights into the function and regulation of alternative splicing and how it evolves.

Using mRNAs lengths to accurately predict the alternatively spliced gene products in Caenorhabditis elegans

MOTIVATION

Computational gene prediction methods are an important component of whole genome analyses. While ab initio gene finders have demonstrated major improvements in accuracy, the most reliable methods are evidence-based gene predictors. These algorithms can rely on several different sources of evidence including predictions from multiple ab initio gene finders, matches to known proteins, sequence conservation and partial cDNAs to predict the final product. Despite the success of these algorithms, prediction of complete gene structures, especially for alternatively spliced products, remains a difficult task.

RESULTS

LOCUS (Length Optimized Characterization of Unknown Spliceforms) is a new evidence-based gene finding algorithm which integrates a length-constraint into a dynamic programming-based framework for prediction of gene products. On a Caenorhabditis elegans test set of alternatively spliced internal exons, its performance exceeds that of current ab initio gene finders and in most cases can accurately predict the correct form of all the alternative products. As the length information used by the algorithm can be obtained in a high-throughput fashion, we propose that integration of such information into a gene-prediction pipeline is feasible and doing so may improve our ability to fully characterize the complete set of mRNAs for a genome.

AVAILABILITY

LOCUS is available from http://ural.wustl.edu/software.html

Can RNA selection pressure distort the measurement of Ka/Ks?

Recently, an interesting question has emerged in the evolutionary interpretation of sequence substitution data as evidence of amino acid selection pressure. Specifically, the Ka/Ks metric was designed to measure selection pressure on amino acid substitutions, assuming that the synonymous substitution rate Ks reflects the neutral nucleotide substitution rate. However, there is increasing evidence for selection pressure at silent sites due to constraints of RNA splicing. Is Ka/Ks an appropriate metric for selection pressure on amino acid substitutions, in the presence of other selection pressures acting only at the RNA level (such as selection for exonic splicing enhancers)? Or can the resulting decreases in Ks from such selection pressures introduce bias into the Ka/Ks metric, so that it no longer gives an accurate measure of amino acid level selection pressure? In this review, we present both mathematical models and empirical evidence for these divergent points of view.

Alternatively and constitutively spliced exons are subject to different evolutionary forces

There has been a controversy on whether alternatively spliced exons (ASEs) evolve faster than constitutively spliced exons (CSEs). Although it has been noted that ASEs are subject to weaker selective constraints than CSEs, so they evolve faster, there have also been studies that indicated slower evolution in ASEs than in CSEs. In this study, we retrieve more than 5,000 human-mouse orthologous exons and calculate the synonymous (KS) and nonsynonymous (KA) substitution rates in these exons. Our results show that ASEs have higher KA values and higher KA/KS ratios than CSEs, indicating faster amino acid-level evolution in ASEs. The faster evolution may be in part due to weaker selective constraints. It is also possible that the faster rate is in part due to faster functional evolution in ASEs. On the other hand, the majority of ASEs have lower KS values than CSEs. With reference to the substitution rate in introns, we show that the KS values in ASEs are close to the neutral substitution rate, whereas the synonymous substitution rate in CSEs has likely been accelerated. The elevated synonymous rate in CSEs is not related to CpG dinucleotides or low-complexity regions of protein but may be weakly related to codon usage bias. The overall trends of higher KA and lower KS in ASEs than in CSEs are also observed in human-rat and mouse-rat comparisons. Therefore, our observations hold for mammals of different molecular clocks.

Characteristics and regulatory elements defining constitutive splicing and different modes of alternative splicing in human and mouse

Alternative splicing is a major contributor to genomic complexity, disease, and development. Previous studies have captured some of the characteristics that distinguish alternative splicing from constitutive splicing. However, most published work only focuses on skipped exons and/or a single species. Here we take advantage of the highly curated data in the MAASE database (see related paper in this issue) to analyze features that characterize different modes of splicing. Our analysis confirms previous observations about alternative splicing, including weaker splicing signals at alternative splice sites, higher sequence conservation surrounding orthologous alternative exons, shorter exon length, and more frequent reading frame maintenance in skipped exons. In addition, our study reveals potentially novel regulatory principles underlying distinct modes of alternative splicing and a role of a specific class of repeat elements (transposons) in the origin/evolution of alternative exons. These features suggest diverse regulatory mechanisms and evolutionary paths for different modes of alternative splicing.

Evidence of functional selection pressure for alternative splicing events that accelerate evolution of protein subsequences

Recently, it was proposed that alternative splicing may act as a mechanism for opening accelerated paths of evolution, by reducing negative selection pressure, but there has been little evidence so far that this mechanism could produce adaptive benefit. Here, we use metrics of very different types of selection pressures [e.g., against amino acid mutations (Ka/Ks), against mutations at synonymous sites (Ks), and for protein reading-frame preservation] to address this question by genomewide analyses of human, chimpanzee, mouse, and rat. These data show that alternative splicing relaxes Ka/Ks selection pressure up to 7-fold, but intriguingly this effect is accompanied by a strong increase in selection pressure against synonymous mutations, which propagates into the adjacent intron, and correlates strongly with the alternative splicing level observed for each exon. These effects are highly local to the alternatively spliced exon. Comparisons of these four genomes consistently show an increase in the density of amino acid mutations (Ka) in alternatively spliced exons and a decrease in the density of synonymous mutations (Ks). This selection pressure against synonymous mutations in alternatively spliced exons was accompanied in all four genomes by a striking increase in selection pressure for protein reading-frame preservation, and both increased markedly with increasing evolutionary age. Restricting our analysis to a subset of exons with strong evidence for biologically functional alternative splicing produced identical results. Thus alternative splicing apparently can create evolutionary "hotspots" within a protein sequence, and these events have evidently been selected for during mammalian evolution.

Distribution of SR protein exonic splicing enhancer motifs in human protein-coding genes

Exonic splicing enhancers (ESEs) are pre-mRNA cis-acting elements required for splice-site recognition. We previously developed a web-based program called ESEfinder that scores any sequence for the presence of ESE motifs recognized by the human SR proteins SF2/ASF, SRp40, SRp55 and SC35 (). Using ESEfinder, we have undertaken a large-scale analysis of ESE motif distribution in human protein-coding genes. Significantly higher frequencies of ESE motifs were observed in constitutive internal protein-coding exons, compared with both their flanking intronic regions and with pseudo exons. Statistical analysis of ESE motif frequency distributions revealed a complex relationship between splice-site strength and increased or decreased frequencies of particular SR protein motifs. Comparison of constitutively and alternatively spliced exons demonstrated slightly weaker splice-site scores, as well as significantly fewer ESE motifs, in the alternatively spliced group. Our results underline the importance of ESE-mediated SR protein function in the process of exon definition, in the context of both constitutive splicing and regulated alternative splicing.

Sequence conservation, relative isoform frequencies, and nonsense-mediated decay in evolutionarily conserved alternative splicing

Studies of expressed sequence tag data sets have revealed large numbers of splicing variants for human genes, but it remains challenging to distinguish functionally important variants from aberrant splicing, clarify the nature of the alternative functions, and understand the signals that regulate splicing choices. To help address these issues, we have constructed and analyzed a large data set of 1,478 exon-skipping alternative splicing (AS) variants evolutionarily conserved in human and mouse. In about one-fifth of cases, one isoform appears subject to nonsense-mediated mRNA decay (NMD), supporting the idea that a major role of AS is to regulate gene expression; one-quarter of these NMD-inducing cases involve a conserved exon whose apparent sole purpose is to mediate destruction of the message when included. We explore sequence conservation likely related to splicing regulation, using in part a measure of the overall amount of conserved information in a sequence, and find that the increased conservation that has been observed within AS exons primarily affects synonymous sites, suggesting that regulatory signals significantly constrain synonymous substitution rates. We show that a lower frequency of the inclusion isoform relative to the exclusion isoform tends to be associated with weaker splice site signals, smaller exon size, and higher intronic sequence conservation, and provide evidence that all of these factors are under selection to control relative isoform frequencies. Some conserved instances of AS appear to represent aberrant splicing events that by chance have occurred in both species, and we develop a nonparametric likelihood approach to identify these.

A novel intronic cis element, ISE/ISS-3, regulates rat fibroblast growth factor receptor 2 splicing through activation of an upstream exon and repression of a downstream exon containing a noncanonical branch point sequence

Mutually exclusive splicing of fibroblast growth factor receptor 2 (FGFR2) exons IIIb and IIIc yields two receptor isoforms, FGFR2-IIIb and -IIIc, with distinctly different ligand binding properties. Several RNA cis elements in the intron (intron 8) separating these exons have been described that are required for splicing regulation. Using a heterologous splicing reporter, we have identified a new regulatory element in this intron that confers cell-type-specific inclusion of an unrelated exon that mirrors its ability to promote cell-type-specific inclusion of exon IIIb. This element promoted inclusion of exon IIIb while at the same time silencing exon IIIc inclusion in cells expressing FGFR2-IIIb; hence, we have termed this element ISE/ISS-3 (for "intronic splicing enhancer-intronic splicing silencer 3"). Silencing of exon IIIc splicing by ISE/ISS-3 was shown to require a branch point sequence (BPS) using G as the primary branch nucleotide. Replacing a consensus BPS with A as the primary branch nucleotide resulted in constitutive splicing of exon IIIc. Our results suggest that the branch point sequence constitutes an important component that can contribute to the efficiency of exon definition of alternatively spliced cassette exons. Noncanonical branch points may thus facilitate cell-type-specific silencing of regulated exons by flanking cis elements.

DEDB: a database of Drosophila melanogaster exons in splicing graph form

Background

A wealth of quality genomic and mRNA/EST sequences in recent years has provided the data required for large-scale genome-wide analysis of alternative splicing. We have capitalized on this by constructing a database that contains alternative splicing information organized as splicing graphs, where all transcripts arising from a single gene are collected, organized and classified. The splicing graph then serves as the basis for the classification of the various types of alternative splicing events.

Description

DEDB is a database of Drosophila melanogaster exons obtained from FlyBase arranged in a splicing graph form that permits the creation of simple rules allowing for the classification of alternative splicing events. Pfam domains were also mapped onto the protein sequences allowing users to access the impact of alternative splicing events on domain organization.

Conclusions

DEDB's catalogue of splicing graphs facilitates genome-wide classification of alternative splicing events for genome analysis. The splicing graph viewer brings together genome, transcript, protein and domain information to facilitate biologists in understanding the implications of alternative splicing.