Characterization of the splice sites in GT-AG and GC-AG introns in higher eukaryotes using full-length cDNAs

For the purpose of analyzing the relation between the splice sites and the order of introns, we conducted the following analysis for the GT-AG and GC-AG splice site groups. First, the pre-mRNAs of H. sapiens, M. musculus, D. melanogaster, A. thaliana and O. sativa were sampled by mapping the full-length cDNA to the genomes. Next, the consensus sequences at different regions of pre-mRNAs were analyzed in the five species. We also investigated the mononucleotide and dinucleotide frequencies in the extensive regions around the 5' splice sites (5'ss) and 3' splice sites (3'ss). As a result, differential frequencies of nucleotides at the first 5'ss in both the GT-AG and GC-AG splice site groups were observed in A. thaliana and O. sativa pre-mRNAs. The trend, which indicates that GC 5'ss possess strong consensus sequences, was observed not only in mammalian pre-mRNAs but also in the pre-mRNAs of D. melanogaster, A. thaliana and O. sativa. Furthermore, we examined the consensus sequences of the constitutive and alternative splice sites. It was suggested that in the case of the alternative GC-AG introns, the tendency to have a weak consensus sequence at 5'ss is different between H. sapiens and M. musculus pre-mRNAs.     

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.     

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.     

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.     

Natural selection affects frequencies of AG and GT dinucleotides at the 5' and 3' ends of exons

GT and AG, located at the 5' and 3' ends of introns, are important for correct splicing. It is anticipated that natural selection decreases frequency of AG and GT near the 5' and 3' ends of exons, preventing appearance of cryptic splicing sites. The data presented in this article support the expectation.     

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.     

Global control of aberrant splice-site activation by auxiliary splicing sequences: evidence for a gradient in exon and intron definition

Auxiliary splicing signals play a major role in the regulation of constitutive and alternative pre-mRNA splicing, but their relative importance in selection of mutation-induced cryptic or de novo splice sites is poorly understood. Here, we show that exonic sequences between authentic and aberrant splice sites that were activated by splice-site mutations in human disease genes have lower frequencies of splicing enhancers and higher frequencies of splicing silencers than average exons. Conversely, sequences between authentic and intronic aberrant splice sites have more enhancers and less silencers than average introns. Exons that were skipped as a result of splice-site mutations were smaller, had lower SF2/ASF motif scores, a decreased availability of decoy splice sites and a higher density of silencers than exons in which splice-site mutation activated cryptic splice sites. These four variables were the strongest predictors of the two aberrant splicing events in a logistic regression model. Elimination or weakening of predicted silencers in two reporters consistently promoted use of intron-proximal splice sites if these elements were maintained at their original positions, with their modular combinations producing expected modification of splicing. Together, these results show the existence of a gradient in exon and intron definition at the level of pre-mRNA splicing and provide a basis for the development of computational tools that predict aberrant splicing outcomes.     

Involvement of SR proteins in mRNA surveillance

Nonsense mutations influence several aspects of gene expression, including mRNA stability and splicing fidelity, but the mechanism by which premature termination codons (PTCs) can apparently affect splice-site selection remains elusive. We used a model human beta-globin gene with duplicated 5' splice sites (5'ss) and found that PTCs inserted between the two 5'ss do not directly influence splicing in this system. Instead, their apparent effect on 5'ss selection in vivo is an indirect result of nonsense-mediated mRNA decay (NMD), as conditions that eliminated NMD also abrogated the effect on splicing. Remarkably, we found an unexpected function of SR proteins in targeting several mRNAs with PTCs to the NMD pathway. Overexpression of various SR proteins strongly enhanced NMD, and this effect required an RS domain. Our data argue against a universal role of PTCs in regulating pre-mRNA splicing and reveal an additional function of SR proteins in eukaryotic gene expression.     

Intronic CA-repeat and CA-rich elements: a new class of regulators of mammalian alternative splicing

We have recently identified an intronic polymorphic CA-repeat region in the human endothelial nitric oxide synthase (eNOS) gene as an important determinant of the splicing efficiency, requiring specific binding of hnRNP L. Here, we analyzed the position requirements of this CA-repeat element, which revealed its potential role in alternative splicing. In addition, we defined the RNA binding specificity of hnRNP L by SELEX: not only regular CA repeats are recognized with high affinity but also certain CA-rich clusters. Therefore, we have systematically searched the human genome databases for CA-repeat and CA-rich elements associated with alternative 5' splice sites (5'ss), followed by minigene transfection assays. Surprisingly, in several specific human genes that we tested, intronic CA RNA elements could function either as splicing enhancers or silencers, depending on their proximity to the alternative 5'ss. HnRNP L was detected specifically bound to these diverse CA elements. These data demonstrated that intronic CA sequences constitute novel and widespread regulatory elements of alternative splicing.     

Determinants of the inherent strength of human 5' splice sites

We previously showed that the authentic 5' splice site (5'ss) of the first exon in the human beta-globin gene is intrinsically stronger than a cryptic 5'ss located 16 nucleotides upstream. Here we examined by mutational analysis the contribution of individual 5'ss nucleotides to discrimination between these two 5'ss. Based on the in vitro splicing efficiencies of a panel of 26 wild-type and mutant substrates in two separate 5'ss competition assays, we established a hierarchy of 5'ss and grouped them into three functional subclasses: strong, intermediate, and weak. Competition between two 5'ss from different subclasses always resulted in selection of the 5'ss that belongs to the stronger subclass. Moreover, each subclass has different characteristic features. Strong and intermediate 5'ss can be distinguished by their predicted free energy of base-pairing to the U1 snRNA 5' terminus (DeltaG). Whereas the extent of splicing via the strong 5'ss correlates well with the DeltaG, this is not the case for competition between intermediate 5'ss. Weak 5'ss were used only when the competing authentic 5'ss was inactivated by mutation. These results indicate that extensive complementarity to U1 snRNA exerts a dominant effect for 5'ss selection, but in the case of competing 5'ss with similarly modest complementarity to U1, the role of other 5'ss features is more prominent. This study reveals the importance of additional submotifs present in certain 5'ss sequences, whose characterization will be critical for understanding 5'ss selection in human genes.     

RNA secondary structure mediates alternative 3'ss selection in Saccharomyces cerevisiae

Alternative splicing is the mechanism by which different combinations of exons in the pre-mRNA give rise to distinct mature mRNAs. This process is mediated by splicing factors that bind the pre-mRNA and affect the recognition of its splicing signals. Saccharomyces species lack many of the regulatory factors present in metazoans. Accordingly, it is generally assumed that the amount of alternative splicing is limited. However, there is recent compelling evidence that yeast have functional alternative splicing, mainly in response to environmental conditions. We have previously shown that sequence and structure properties of the pre-mRNA could explain the selection of 3' splice sites (ss) in Saccharomyces cerevisiae. In this work, we extend our previous observations to build a computational classifier that explains most of the annotated 3'ss in the CDS and 5' UTR of this organism. Moreover, we show that the same rules can explain the selection of alternative 3'ss. Experimental validation of a number of predicted alternative 3'ss shows that their usage is low compared to annotated 3'ss. The majority of these alternative 3'ss introduce premature termination codons (PTCs), suggesting a role in expression regulation. Furthermore, a genome-wide analysis of the effect of temperature, followed by experimental validation, yields only a small number of changes, indicating that this type of regulation is not widespread. Our results are consistent with the presence of alternative 3'ss selection in yeast mediated by the pre-mRNA structure, which can be responsive to external cues, like temperature, and is possibly related to the control of gene expression.     

Functional selection of splicing enhancers that stimulate trans-splicing in vitro.

The role of exonic sequences in naturally occurring trans-splicing has not been explored in detail. Here, we have identified trans-splicing enhancers through the use of an iterative selection scheme. Several classes of enhancer sequences were identified that led to dramatic increases in trans-splicing efficiency. Two sequence families were investigated in detail. These include motifs containing the element (G/C)GAC(G/C) and also 5' splice site-like sequences. Distinct elements were tested for their ability to function as splicing enhancers and in competition experiments. In addition, discrete trans-acting factors were identified. This work demonstrates that splicing enhancers are able to effect a large increase in trans-splicing efficiency and that the process of exon definition is able to positively enhance trans-splicing even though the reaction itself is independent of the need for the 5' end of U1 snRNA. Due to the presence of internal introns in messages that are trans-spliced, the natural arrangement of 5' splice sites downstream of trans-splicing acceptors may lead to a general promotion of this unusual reaction.     

The emergence of alternative 3' and 5' splice site exons from constitutive exons

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

Accumulation of GC donor splice signals in mammals

Abstract

The GT dinucleotide in the first two intron positions is the most conserved element of the U2 donor splice signals. However, in a small fraction of donor sites, GT is replaced by GC. A substantial enrichment of GC in donor sites of alternatively spliced genes has been observed previously in human, nematode and Arabidopsis, suggesting that GC signals are important for regulation of alternative splicing. We used parsimony analysis to reconstruct evolution of donor splice sites and inferred 298 GT > GC conversion events compared to 40 GC > GT conversion events in primate and rodent genomes. Thus, there was substantive accumulation of GC donor splice sites during the evolution of mammals. Accumulation of GC sites might have been driven by selection for alternative splicing.

Reviewers

This article was reviewed by Jerzy Jurka and Anton Nekrutenko. For the full reviews, please go to the Reviewers' Reports section.     

Secondary structure is required for 3' splice site recognition in yeast

Higher order RNA structures can mask splicing signals, loop out exons, or constitute riboswitches all of which contributes to the complexity of splicing regulation. We identified a G to A substitution between branch point (BP) and 3' splice site (3'ss) of Saccharomyces cerevisiae COF1 intron, which dramatically impaired its splicing. RNA structure prediction and in-line probing showed that this mutation disrupted a stem in the BP-3'ss region. Analyses of various COF1 intron modifications revealed that the secondary structure brought about the reduction of BP to 3'ss distance and masked potential 3'ss. We demonstrated the same structural requisite for the splicing of UBC13 intron. Moreover, RNAfold predicted stable structures for almost all distant BP introns in S. cerevisiae and for selected examples in several other Saccharomycotina species. The employment of intramolecular structure to localize 3'ss for the second splicing step suggests the existence of pre-mRNA structure-based mechanism of 3'ss recognition.     

The 5′ UTR intron-mediated enhancement of constitutive splicing of the tobacco microsome ω-3 fatty acid desaturase gene

Several plant genes have their first intron in the 5′ untranslated region (5′ UTR), and such 5′ UTR introns often show several biological functions, including the intron-mediated enhancement of protein expression through an increase of mRNA level (IME), intron-dependent spatial expression, and intron-mediated enhancement of translation. Here, we show another function of the 5′ UTR intron, i.e., the 5′ UTR intron-mediated enhancement of constitutive splicing. The NtFAD3 gene, which encodes a tobacco microsome ω-3 fatty acid desaturase, has a 552 nucleotide-long 5′ UTR intron (intron 1), and the other seven introns are located in the coding sequence. The splicing of the 5′ half region of the NtFAD3 was studied through an in vivo splicing assay using Arabidopsis leaf explants. The low splicing efficiency of intron 2 was much improved when the assay construct harbored intron 1. Deletion of intron 1 and the replacement of intron 1 to the NtFAD3 intron 8 decreased the splicing efficiency of intron 2. The splicing enhancers were redundant and dispersed in the 5′ splice site-proximal, 284-nucleotides region of intron 1. In addition, the interaction among the cis-elements, i.e., the splicing enhancers in the intron 1 and exon 2, were necessary for the efficient splicing of intron 2. The 5′ UTR intron-mediated constitutive splicing was partially inhibited when an SR-like protein, SR45, was deficient. These results indicated a novel function of the 5′ UTR intron, namely an enhancement of the constitutive splicing.