Comprehensive splice-site analysis using comparative genomics

We have collected over half a million splice sites from five species-Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans and Arabidopsis thaliana-and classified them into four subtypes: U2-type GT-AG and GC-AG and U12-type GT-AG and AT-AC. We have also found new examples of rare splice-site categories, such as U12-type introns without canonical borders, and U2-dependent AT-AC introns. The splice-site sequences and several tools to explore them are available on a public website (SpliceRack). For the U12-type introns, we find several features conserved across species, as well as a clustering of these introns on genes. Using the information content of the splice-site motifs, and the phylogenetic distance between them, we identify: (i) a higher degree of conservation in the exonic portion of the U2-type splice sites in more complex organisms; (ii) conservation of exonic nucleotides for U12-type splice sites; (iii) divergent evolution of C.elegans 3' splice sites (3'ss) and (iv) distinct evolutionary histories of 5' and 3'ss. Our study proves that the identification of broad patterns in naturally-occurring splice sites, through the analysis of genomic datasets, provides mechanistic and evolutionary insights into pre-mRNA splicing.     

Human GC-AG alternative intron isoforms with weak donor sites show enhanced consensus at acceptor exon positions

It has been previously observed that the intrinsically weak variant GC donor sites, in order to be recognized by the U2-type spliceosome, possess strong consensus sequences maximized for base pair formation with U1 and U5/U6 snRNAs. However, variability in signal strength is a fundamental mechanism for splice site selection in alternative splicing. Here we report human alternative GC-AG introns (for the first time from any species), and show that while constitutive GC-AG introns do possess strong signals at their donor sites, a large subset of alternative GC-AG introns possess weak consensus sequences at their donor sites. Surprisingly, this subset of alternative isoforms shows strong consensus at acceptor exon positions 1 and 2. The improved consensus at the acceptor exon can facilitate a strong interaction with U5 snRNA, which tethers the two exons for ligation during the second step of splicing. Further, these isoforms nearly always possess alternative acceptor sites and exhibit particularly weak polypyrimidine tracts characteristic of AG-dependent introns. The acceptor exon nucleotides are part of the consensus required for the U2AF³⁵-mediated recognition of AG in such introns. Such improved consensus at acceptor exons is not found in either normal or alternative GT-AG introns having weak donor sites or weak polypyrimidine tracts. The changes probably reflect mechanisms that allow GC-AG alternative intron isoforms to cope with two conflicting requirements, namely an apparent need for differential splice strength to direct the choice of alternative sites and a need for improved donor signals to compensate for the central mismatch base pair (C-A) in the RNA duplex of U1 snRNA and the pre-mRNA. The other important findings include (i) one in every twenty alternative introns is a GC-AG intron, and (ii) three of every five observed GC-AG introns are alternative isoforms.     

The presence of GC-AG introns in Neurospora crassa and other euascomycetes determined from analyses of complete genomes: implications for automated gene prediction

A combination of experimental and computational approaches was employed to identify introns with noncanonical GC-AG splice sites (GC-AG introns) within euascomycete genomes. Evaluation of 2335 cDNA-confirmed introns from Neurospora crassa revealed 27 such introns (1.2%). A similar frequency (1.0%) of GC-AG introns was identified in Fusarium graminearum, in which 3 of 292 cDNA-confirmed introns contained GC-AG splice sites. Computational analyses of the N. crassa genome using a GC-AG intron consensus sequence identified an additional 20 probable GC-AG introns in this fungus. For 8 of the 47 GC-AG introns identified in N. crassa a GC donor site is also present in a homolog from Magnaporthe grisea, F. graminearum, or Aspergillus nidulans. In most cases, however, homologs in these fungi contain a GT-AG intron or no intron at the corresponding position. These findings have important implications for fungal genome annotation, as the automated annotations of euascomycete genomes incorrectly identified intron boundaries for all of the confirmed and probable GC-AG introns reported here.     

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

We previously reported a computational approach to infer alternative splicing patterns from Mus musculus full-length cDNA clones and microarray data. Although we predicted a large number of unreported splice variants, the general mechanisms regulating alternative splicing were yet unknown. In the present study, we compared alternative exons and constitutive exons in terms of splice-site strength and frequency of potential regulatory sequences. These regulatory features were further compared among five different species: Homo sapiens, M. musculus, Arabidopsis thaliana, Oryza sativa, and Drosophila melanogaster. Solid statistical validations of our comparative analyses indicated that alternative exons have (1) weaker splice sites and (2) more potential regulatory sequences than constitutive exons. Based on our observations, we propose a combinatorial model of alternative splicing mechanisms, which suggests that alternative exons contain weak splice sites regulated alternatively by potential regulatory sequences on the exons.     

Analysis of the role of Caenorhabditis elegans GC-AG introns in regulated splicing

GC-AG introns represent 0.7% of total human pre-mRNA introns. To study the function of GC-AG introns in splicing regulation, 196 cDNA-confirmed GC-AG introns were identified in Caenorhabditis elegans. These represent 0.6% of the cDNA- confirmed intron data set for this organism. Eleven of these GC-AG introns are involved in alternative splicing. In a comparison of the genomic sequences of homologous genes between C.elegans and Caenorhabditis briggsae for 26 GC-AG introns, the C at the +2 position is conserved in only five of these introns. A system to experimentally test the function of GC-AG introns in alternative splicing was developed. Results from these experiments indicate that the conserved C at the +2 position of the tenth intron of the let-2 gene is essential for developmentally regulated alternative splicing. This C allows the splice donor to function as a very weak splice site that works in balance with an alternative GT splice donor. A weak GT splice donor can functionally replace the GC splice donor and allow for splicing regulation. These results indicate that while the majority of GC-AG introns appear to be constitutively spliced and have no evolutionary constraints to prevent them from being GT-AG introns, a subset of GC-AG introns is involved in alternative splicing and the C at the +2 position of these introns can have an important role in splicing regulation.     

Conservation of an open-reading frame as an element affecting 5' splice site selection

Splice site selection is a key element of pre-mRNA splicing and involves specific recognition of consensus sequences at the 5(') and 3(') splice sites. Evidently, the compliance of a given sequence with the consensus 5(') splice site sequence is not sufficient to define it as a functional 5(') splice site, because not all sequences that conform with the consensus are used for splicing. We have previously hypothesized that the necessity to avoid the inclusion of premature termination codons within mature mRNAs may serve as a criterion that differentiates normal 5(') splice sites from unused (latent) ones. We further provided experimental support to this idea, by analyzing the splicing of pre-mRNAs in which in-frame stop codons upstream of a latent 5(') splice site were mutated, and showing that splicing using the latent site is indeed activated by such mutations. Here we evaluate this hypothesis by a computerized survey for latent 5(') splice sites in 446 protein-coding human genes. This data set contains 2311 introns, in which we found 10490 latent 5(') splice sites. The utilization of 10045 (95.8%) of these sites for splicing would have led to the inclusion of an in-frame stop codon within the resultant mRNA. The validity of this finding is confirmed here by statistical analyses. This finding, together with our previous experimental results, invokes a nuclear scanning mechanism, as part of the splicing machine, which identifies in-frame stop codons within the pre-mRNA and prevents splicing that could lead to the formation of a prematurely terminated protein.     

SL1 trans-splicing specified by AU-rich synthetic RNA inserted at the 5' end of Caenorhabditis elegans pre-mRNA.

In Caenorhabditis elegans, pre-mRNAs of many genes are trans-spliced to one of two spliced leaders, SL1 or SL2. Some of those that receive exclusively SL1 have been characterized as having at their 5' ends outrons, AU-rich sequences similar to introns followed by conventional 3' splice sites. Comparison of outrons from many different SL1-specific C. elegans genes has not revealed the presence of any consensus sequence that might encode SL1-specificity. In order to determine what parameters influence the splicing of SL1, we performed in vivo experiments with synthetic splice sites. Synthetic AU-rich RNA, 51 nt or longer, placed upstream of a consensus 3' splice site resulted in efficient trans-splicing. With all sequences tested, this trans-splicing was specifically to SL1. Thus, no information beyond the presence of AU-rich RNA at least as long as the minimum-length C. elegans intron, followed by a 3' splice site, is required to specify trans-splicing or for strict SL1 specificity.     

Identification and characterization by antisense oligonucleotides of exon and intron sequences required for splicing.

Certain thalassemic human beta-globin pre-mRNAs carry mutations that generate aberrant splice sites and/or activate cryptic splice sites, providing a convenient and clinically relevant system to study splice site selection. Antisense 2'-O-methyl oligoribonucleotides were used to block a number of sequences in these pre-mRNAs and were tested for their ability to inhibit splicing in vitro or to affect the ratio between aberrantly and correctly spliced products. By this approach, it was found that (i) up to 19 nucleotides upstream from the branch point adenosine are involved in proper recognition and functioning of the branch point sequence; (ii) whereas at least 25 nucleotides of exon sequences at both 3' and 5' ends are required for splicing, this requirement does not extend past the 5' splice site sequence of the intron; and (iii) improving the 5' splice site of the internal exon to match the consensus sequence strongly decreases the accessibility of the upstream 3' splice site to antisense 2'-O-methyl oligoribonucleotides. This result most likely reflects changes in the strength of interactions near the 3' splice site in response to improvement of the 5' splice site and further supports the existence of communication between these sites across the exon.     

SpliceDB: database of canonical and non-canonical mammalian splice sites

A database (SpliceDB) of known mammalian splice site sequences has been developed. We extracted 43 337 splice pairs from mammalian divisions of the gene-centered Infogene database, including sites from incomplete or alternatively spliced genes. Known EST sequences supported 22 815 of them. After discarding sequences with putative errors and ambiguous location of splice junctions the verified dataset includes 22 489 entries. Of these, 98.71% contain canonical GT-AG junctions (22 199 entries) and 0.56% have non-canonical GC-AG splice site pairs. The remainder (0.73%) occurs in a lot of small groups (with a maximum size of 0.05%). We especially studied non-canonical splice sites, which comprise 3.73% of GenBank annotated splice pairs. EST alignments allowed us to verify only the exonic part of splice sites. To check the conservative dinucleotides we compared sequences of human non-canonical splice sites with sequences from the high throughput genome sequencing project (HTG). Out of 171 human non-canonical and EST-supported splice pairs, 156 (91.23%) had a clear match in the human HTG. They can be classified after sequence analysis as: 79 GC-AG pairs (of which one was an error that corrected to GC-AG), 61 errors corrected to GT-AG canonical pairs, six AT-AC pairs (of which two were errors corrected to AT-AC), one case was produced from a non-existent intron, seven cases were found in HTG that were deposited to GenBank and finally there were only two other cases left of supported non-canonical splice pairs. The information about verified splice site sequences for canonical and non-canonical sites is presented in SpliceDB with the supporting evidence. We also built weight matrices for the major splice groups, which can be incorporated into gene prediction programs. SpliceDB is available at the computational genomic Web server of the Sanger Centre: http://genomic.sanger.ac. uk/spldb/SpliceDB.html and at http://www.softberry. com/spldb/SpliceDB.html.     

Nuclear pre-mRNA processing in plants: distinct modes of 3''-splice-site selection in plants and animals.

The report that human growth hormone pre-mRNA is not processed in transgenic plant tissues (A. Barta, K. Sommergruber, D. Thompson, K. Hartmuth, M.A. Matzke, and A.J.M. Matzke, Plant Mol. Biol. 6:347-357, 1986) has suggested that differences in mRNA splicing processes exist between plants and animals. To gain more information about the specificity of plant pre-mRNA processing, we have compared the splicing of the soybean leghemoglobin pre-mRNA with that of the human beta-globin pre-mRNA in transfected plant (Orychophragmus violaceus and Nicotiana tabacum) protoplasts and mammalian (HeLa) cells. Of the three introns of leghemoglobin pre-mRNA, only intron 2 was correctly and efficiently processed in HeLa cells. The 5' splice sites of the remaining two introns were faithfully recognized, but correct processing of the 3' sites took place only rarely (intron 1) or not at all (intron 3); cryptic 3' splice sites were used instead. While the first intron in human beta-globin pre-mRNA was not spliced in transfected plant protoplasts, intron 2 processing occurred at a low level, indicating that some mammalian introns can be recognized by the plant intron-splicing machinery. However, excision of intron 2 proved to be incorrect, involving the authentic 5' splice site and a cryptic 3' splice site. Our results indicate that the mechanism of 3'-splice-site selection during intron excision differs between plants and animals. This conclusion is supported by analysis of the 3'-splice-site consensus sequences in animal and plant introns which revealed that polypyrimidine tracts, characteristic of animal introns, are not present in plant pre-mRNAs. It is proposed that an elevated AU content of plant introns is important for their processing.     

The apoptosis-promoting factor TIA-1 is a regulator of alternative pre-mRNA splicing

We report here that the apoptosis-promoting protein TIA-1 regulates alternative pre-mRNA splicing of the Drosophila melanogaster gene male-specific-lethal 2 and of the human apoptotic gene Fas. TIA-1 associates selectively with pre-mRNAs that contain 5' splice sites followed by U-rich sequences. TIA-1 binding to the U-rich stretches facilitates 5' splice site recognition by U1 snRNP. This activity is critical for activation of the weak 5' splice site of msl-2 and for modulating the choice of splice site partner in Fas. Structural and functional similarities with the Saccharomyces cerevisiae splicing factor Nam8 suggest striking evolutionary conservation of a mechanism of pre-mRNA splicing regulation that controls biological processes as diverse as meiosis in yeast, dosage compensation in fruit flies, or programmed cell death in humans.     

Sensitivity of splice sites to antisense oligonucleotides in vivo

A series of HeLa cell lines which stably express beta-globin pre-mRNAs carrying point mutations at nt 654, 705, or 745 of intron 2 has been developed. The mutations generate aberrant 5' splice sites and activate a common 3' cryptic splice site upstream leading to aberrantly spliced beta-globin mRNA. Antisense oligonucleotides, which in vivo blocked aberrant splice sites and restored correct splicing of the pre-mRNA, revealed major differences in the sensitivity of these sites to antisense probes. Although the targeted pre-mRNAs differed only by single point mutations, the effective concentrations of the oligonucleotides required for correction of splicing varied up to 750-fold. The differences among the aberrant 5' splice sites affected sensitivity of both the 5' and 3' splice sites; in particular, sensitivity of both splice sites was severely reduced by modification of the aberrant 5' splice sites to the consensus sequence. These results suggest large differences in splicing of very similar pre-mRNAs in vivo. They also indicate that antisense oligonucleotides may provide useful tools for studying the interactions of splicing machinery with pre-mRNA.     

Analysis of canonical and non-canonical splice sites in mammalian genomes

A set of 43 337 splice junction pairs was extracted from mammalian GenBank annotated genes. Expressed sequence tag (EST) sequences support 22 489 of them. Of these, 98.71% contain canonical dinucleotides GT and AG for donor and acceptor sites, respectively; 0.56% hold non-canonical GC-AG splice site pairs; and the remaining 0.73% occurs in a lot of small groups (with a maximum size of 0.05%). Studying these groups we observe that many of them contain splicing dinucleotides shifted from the annotated splice junction by one position. After close examination of such cases we present a new classification consisting of only eight observed types of splice site pairs (out of 256 a priori possible combinations). EST alignments allow us to verify the exonic part of the splice sites, but many non-canonical cases may be due to intron sequencing errors. This idea is given substantial support when we compare the sequences of human genes having non-canonical splice sites deposited in GenBank by high throughput genome sequencing projects (HTG). A high proportion (156 out of 171) of the human non-canonical and EST-supported splice site sequences had a clear match in the human HTG. They can be classified after corrections as: 79 GC-AG pairs (of which one was an error that corrected to GC-AG), 61 errors that were corrected to GT-AG canonical pairs, six AT-AC pairs (of which two were errors that corrected to AT-AC), one case was produced from non-existent intron, seven cases were found in HTG that were deposited to GenBank and finally there were only two cases left of supported non-canonical splice sites. If we assume that approximately the same situation is true for the whole set of annotated mammalian non-canonical splice sites, then the 99.24% of splice site pairs should be GT-AG, 0.69% GC-AG, 0.05% AT-AC and finally only 0.02% could consist of other types of non-canonical splice sites. We analyze several characteristics of EST-verified splice sites and build weight matrices for the major groups, which can be incorporated into gene prediction programs. We also present a set of EST-verified canonical splice sites larger by two orders of magnitude than the current one (22 199 entries versus approximately 600) and finally, a set of 290 EST-supported non-canonical splice sites. Both sets should be significant for future investigations of the splicing mechanism.     

Intron splicing: a conserved internal signal in introns of Drosophila pre-mRNAs.

The introns of Drosophila pre-mRNAs have been analysed for conserved internal sequence elements near the 3' intron boundary similar to the T-A-C-T-A-A-C in yeast introns and the C/T-T-A/G-A-C/T in introns of other organisms. Such conserved internal elements are the 3' splice signals recognized in intron splicing. In the lariat splicing mechanism, the G at the 5' end of an intron joins covalently to the last A of a 3' splice signal to form a branch point in a splicing intermediate. Analysis of 39 published sequences of Drosophila introns reveals that potential 3' splice signals with the consensus C/T-T-A/G-A-C/T are present in 18 cases. In 17 of the remaining cases signals are present which vary from this consensus just in the middle or last position. In Drosophila introns the 3' splice signal is usually located in a discrete region between 18 and 35 nucleotides upstream from the 3' splice point. We note that the Drosophila small nuclear U2-RNA has sequences complementary to C-T-G-A-T, one variant of the signal, and to C-A-G, one variant of the 3' terminus of an intron. We also note that the absence of any A-G between -3 and -19 from the 3' splice point may be an essential feature of a strong 3' boundary.