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.     

Analysis and recognition of 5' UTR intron splice sites in human pre-mRNA

Prediction of splice sites in non-coding regions of genes is one of the most challenging aspects of gene structure recognition. We perform a rigorous analysis of such splice sites embedded in human 5' untranslated regions (UTRs), and investigate correlations between this class of splice sites and other features found in the adjacent exons and introns. By restricting the training of neural network algorithms to 'pure' UTRs (not extending partially into protein coding regions), we for the first time investigate the predictive power of the splicing signal proper, in contrast to conventional splice site prediction, which typically relies on the change in sequence at the transition from protein coding to non-coding. By doing so, the algorithms were able to pick up subtler splicing signals that were otherwise masked by 'coding' noise, thus enhancing significantly the prediction of 5' UTR splice sites. For example, the non-coding splice site predicting networks pick up compositional and positional bias in the 3' ends of non-coding exons and 5' non-coding intron ends, where cytosine and guanine are over-represented. This compositional bias at the true UTR donor sites is also visible in the synaptic weights of the neural networks trained to identify UTR donor sites. Conventional splice site prediction methods perform poorly in UTRs because the reading frame pattern is absent. The NetUTR method presented here performs 2-3-fold better compared with NetGene2 and GenScan in 5' UTRs. We also tested the 5' UTR trained method on protein coding regions, and discovered, surprisingly, that it works quite well (although it cannot compete with NetGene2). This indicates that the local splicing pattern in UTRs and coding regions is largely the same. The NetUTR method is made publicly available at www.cbs.dtu.dk/services/NetUTR.     

Using estimative reaction free energy to predict splice sites and their flanking competitors

A crucial part in the gene structure prediction is to identify the accurate splice sites, not only constitutive but also alternative ones. Here, we use the maximum information principle (MIP) to analyze the conservative segments around splice sites. According to the MIP, a reaction free energy (RFE) expression is deduced, which can be employed to estimate the free energy change during splicing reaction involving a donor or acceptor site. The expression contains not only the background probability factors, but also all kinds of dependencies among both adjacent and non-adjacent bases. We apply the RFE expression to recognize splice sites and their flanking competitors in human genes, the results show high sensitivity and specificity, so the RFE expression accords well with the splicing reaction process. Moreover, the RFE expression is better than previous methods for predicting competitors of splice sites, and it outperforms the reaction free energy subtraction (RFES), that implies RFE competition between a given splice site and its flanking competitor may not be an only primary factor for alternative splice site selection. The work is helpful to not only the understanding of splicing reaction from its relation to MIP, but also the research on computational recognition of splicing sites and alternative splice events.     

基于RNA-Seq数据识别果蝇剪接位点和可变剪接事件

完整基因结构的预测是当前生命科学研究的一个重要基础课题,其中一个关键环节是剪接位点和各种可变剪接事件的精确识别.基于转录组测序（RNA-seq）数据,识别剪接位点和可变剪接事件是近几年随着新一代测序技术发展起来的新技术策略和方法.本工作基于黑腹果蝇睾丸RNA-seq数据,使用TopHat软件成功识别出39718个果蝇剪接位点,其中有10584个新剪接位点.同时,基于剪接位点的不同组合,针对各类型可变剪接特征开发出计算识别算法,成功识别了8477个可变剪接事件（其中新识别的可变剪接事件3922个）,包括可变供体位点、可变受体位点、内含子保留和外显子缺失4种类型.RT-PCR实验验证了2个果蝇基因上新识别的可变剪接事件,发现了全新的剪接异构体.进一步表明,RNA-seq数据可有效应用于识别剪接位点和可变剪接事件,为深入揭示剪接机制及可变剪接生物学功能提供新思路和新手段.     

Pathways for selection of 5' splice sites by U1 snRNPs and SF2/ASF.

We have used protection against ribonuclease H to investigate the mechanisms by which U1 small nuclear ribonucleoprotein particles (snRNPs) determine the use of two alternative 5' splice sites. The initial binding of U1 snRNPs to alternative consensus splice sites was indiscriminate, and on a high proportion of pre-mRNA molecules both sites were occupied simultaneously. When the sites were close, this inhibited splicing. We propose that double occupancy leads to the use of the downstream site for splicing and that this is the cause of the proximity effect seen with strong alternative splice sites. This model predicts that splicing to an upstream site of any strength requires a low affinity of U1 snRNPs for the downstream site. This prediction was tested both by cleaving the 5' end of U1 snRNA and by altering the sequence of the downstream site of an adenovirus E1A gene. The enhancement of downstream 5' splice site use by splicing factor SF2/ASF appears to be mediated by an increase in the strength of U1 snRNP binding to all sites indiscriminately.     

Intronic sequences and 3' splice sites control Rous sarcoma virus RNA splicing.

cis-acting sequences of Rous sarcoma virus (RSV) RNA involved in control of the incomplete splicing that is part of the retroviral life cycle have been studied. The 5' and two alternative 3' splice sites, as well as negative regulator of splicing element in the intron, have been introduced into chimeric constructs, and their responsive roles in splicing inhibition have been evaluated by transient transfection experiments. Although the RSV 5' splice site was used efficiently in these assays, substrates containing either the RSV env or the RSV src 3' splice site were not spliced completely, resulting in 40 to 50% unspliced RNA. Addition of the negative regulator of splicing element to substrates containing RSV 3' splice sites resulted in greater inhibition of splicing (70 to 80% unspliced RNA), suggesting that the two elements function independently and additively. Deletion of sequences more than 70 nucleotides upstream of the src 3' splice site resulted in efficient splicing at this site, suggesting that inefficient usage is not inherent in this splice site but is instead due to to sequences upstream of it. Insertion of these upstream sequences into the intron of a heterologous pre-mRNA resulted in partial inhibition of its splicing. In addition, secondary structure interactions were predicted to occur between the src 3' splice site and the inhibitory sequences upstream of it. Thus, RSV splicing control involves both intronic sequences and 3' splice sites, with different mechanisms involved in the underutilization of the env and src splice acceptor sites.     

Efficient internal exon recognition depends on near equal contributions from the 3' and 5' splice sites

Pre-mRNA splicing is carried out by the spliceosome, which identifies exons and removes intervening introns. In vertebrates, most splice sites are initially recognized by the spliceosome across the exon, because most exons are small and surrounded by large introns. This gene architecture predicts that efficient exon recognition depends largely on the strength of the flanking 3' and 5' splice sites. However, it is unknown if the 3' or the 5' splice site dominates the exon recognition process. Here, we test the 3' and 5' splice site contributions towards efficient exon recognition by systematically replacing the splice sites of an internal exon with sequences of different splice site strengths. We show that the presence of an optimal splice site does not guarantee exon inclusion and that the best predictor for exon recognition is the sum of both splice site scores. Using a genome-wide approach, we demonstrate that the combined 3' and 5' splice site strengths of internal exons provide a much more significant separator between constitutive and alternative exons than either the 3' or the 5' splice site strength alone.     

Influences of separation and adjacent sequences on the use of alternative 5' splice sites.

Single nucleotide changes to the sequence between two alternative 5' splice sites, separated by 25 nucleotides in a beta-globin gene derivative, caused substantial shifts in pre-mRNA splicing preferences, both in vivo and in vitro. An activating sequence for splicing was located. Models for the recognition by U1 small nuclear ribonucleoproteins (snRNPs) of competing 5' splice sites were tested by altering the distance separating the two sites. Use of the upstream splice site declined sharply when it was separated from the downstream (natural) site by distances of 40 nucleotides or more. This effect was reversed in vivo, but not in vitro, by altering the upstream sequence to that of a consensus 5' splice site sequence. Dilution of an extract used for splicing in vitro shifted preferences when the sites were close towards the downstream site. We conclude that the mechanism of selection depends on the distance apart of the potential splice sites and that with close sites steric interference between factors bound to both sites may impede splicing and affect splicing preferences.     

Hierarchy for 5' splice site preference determined in vivo

The relationship between preferences among alternative 5' splice sites and their sequences has been investigated for 37 sequences by assessing their use in splicing relative to the 5' splice site of IVS-2 of rabbit beta-globin. There are strong correlations between the intrinsic strength of a 5' splice site and both the extent to which it resembles the consensus sequence and the calculated stability of its interactions with U1 small nuclear RNA. However, present methods of calculating either of the latter values do not allow predictions to be made of the relative preferences among a small number of sequences.     

A novel protein factor is required for use of distal alternative 5' splice sites in vitro.

Adenovirus E1A pre-mRNA was used as a model to examine alternative 5' splice site selection during in vitro splicing reactions. Strong preference for the downstream 13S 5' splice site over the upstream 12S or 9S 5' splice sites was observed. However, the 12S 5' splice site was used efficiently when a mutant pre-mRNA lacking the 13S 5' splice site was processed, and 12S splicing from this substrate was not reduced by 13S splicing from a separate pre-mRNA, demonstrating that 13S splicing reduced 12S 5' splice site selection through a bona fide cis-competition. DEAE-cellulose chromatography of nuclear extract yielded two fractions with different splicing activities. The bound fraction contained all components required for efficient splicing of simple substrates but was unable to utilize alternative 5' splice sites. In contrast, the flow-through fraction, which by itself was inactive, contained an activity required for alternative splicing and was shown to stimulate 12S and 9S splicing, while reducing 13S splicing, when added to reactions carried out by the bound fraction. Furthermore, the activity, which we have called distal splicing factor (DSF), enhanced utilization of an upstream 5' splice site on a simian virus 40 early pre-mRNA, suggesting that the factor acts in a position-dependent, substrate-independent fashion. Several lines of evidence are presented suggesting that DSF is a non-small nuclear ribonucleoprotein protein. Finally, we describe a functional interaction between DSF and ASF, a protein that enhances use of downstream 5' splice sites.     

Sequestering of the 3' splice site in a theophylline-responsive riboswitch allows ligand-dependent control of alternative splicing

Despite the important role of alternative splicing in various aspects of biological processes, our ability to regulate this process at will remains a challenge. In this report, we asked whether a theophylline-responsive riboswitch could be adapted to manipulate alternative splicing. We constructed a pre-mRNA containing a single upstream 5' splice site and two 3' splice sites, of which the proximal 3' splice site is embedded in theophylline-responsive riboswitch. We show that this pre-mRNA spliced with preferential utilization of proximal 3' splice site in vitro. However, addition of theophylline to the splicing reaction promoted splicing at distal 3' splice site thereby changing the ratio of distal-to-proximal 3' splice site usage by more than twofold. Our data suggest that theophylline influenced 3' splice site choice without affecting the kinetics of the splicing reaction. We conclude that an in vitro selected riboswitch can be adapted to control alternative splicing, which may find many applications in basic, biotechnological, and biomedical research.     

Optimization of a weak 3' splice site counteracts the function of a bovine papillomavirus type 1 exonic splicing suppressor in vitro and in vivo

Alternative splicing is a critical component of the early to late switch in papillomavirus gene expression. In bovine papillomavirus type 1 (BPV-1), a switch in 3' splice site utilization from an early 3' splice site at nucleotide (nt) 3225 to a late-specific 3' splice site at nt 3605 is essential for expression of the major capsid (L1) mRNA. Three viral splicing elements have recently been identified between the two alternative 3' splice sites and have been shown to play an important role in this regulation. A bipartite element lies approximately 30 nt downstream of the nt 3225 3' splice site and consists of an exonic splicing enhancer (ESE), SE1, followed immediately by a pyrimidine-rich exonic splicing suppressor (ESS). A second ESE (SE2) is located approximately 125 nt downstream of the ESS. We have previously demonstrated that the ESS inhibits use of the suboptimal nt 3225 3' splice site in vitro through binding of cellular splicing factors. However, these in vitro studies did not address the role of the ESS in the regulation of alternative splicing. In the present study, we have analyzed the role of the ESS in the alternative splicing of a BPV-1 late pre-mRNA in vivo. Mutation or deletion of just the ESS did not significantly change the normal splicing pattern where the nt 3225 3' splice site is already used predominantly. However, a pre-mRNA containing mutations in SE2 is spliced predominantly using the nt 3605 3' splice site. In this context, mutation of the ESS restored preferential use of the nt 3225 3' splice site, indicating that the ESS also functions as a splicing suppressor in vivo. Moreover, optimization of the suboptimal nt 3225 3' splice site counteracted the in vivo function of the ESS and led to preferential selection of the nt 3225 3' splice site even in pre-mRNAs with SE2 mutations. In vitro splicing assays also showed that the ESS is unable to suppress splicing of a pre-mRNA with an optimized nt 3225 3' splice site. These data confirm that the function of the ESS requires a suboptimal upstream 3' splice site. A surprising finding of our study is the observation that SE1 can stimulate both the first and the second steps of splicing.     

One parameter to describe the mechanism of splice sites competition

The choice of a splice site is not only related to its own intrinsic strength, but also is influenced by its flanking competitors. Splice site competition is an important mechanism for splice site prediction, especially, it is a new insight for alternative splice site prediction. In this paper, the position weight matrix scoring function is used to represent splice site strength, and the mechanism of splice site competition is described by only one parameter: scoring function subtraction. While applying on the alternative splice site prediction, based on the only one parameter, 68.22% of donor sites and 70.86% of acceptor sites are correctly classified into alternative and constitutive. The prediction abilities are approximately equal to the recent method which is based on the mechanism of splice site competition. The results reveal that the scoring function subtraction is the best parameter to describe the mechanism of splice sites competition.     

TIA-1 and TIAR activate splicing of alternative exons with weak 5' splice sites followed by a U-rich stretch on their own pre-mRNAs

TIA-1 has recently been shown to activate splicing of specific pre-mRNAs transcribed from transiently transfected minigenes, and of some 5' splice sites in vitro, but has not been shown to activate splicing of any endogenous pre-mRNA. We show here that overexpression of TIA-1 or the related protein TIAR has little effect on splicing of several endogenous pre-mRNAs containing alternative exons, but markedly activates splicing of some normally rarely used alternative exons on the TIA-1 and TIAR pre-mRNAs. These exons have weak 5' splice sites followed by U-rich stretches. When the U-rich stretch following the 5' splice site of a TIA-1 alternative exon was deleted, TIAR overexpression induced use of a cryptic 5' splice site also followed by a U-rich stretch in place of the original splice site. Using in vitro splicing assays, we have shown that TIA-1 is directly involved in activating the 5' splice sites of the TIAR alternative exons. Activation requires a downstream U-rich stretch of at least 10 residues. Our results confirm that TIA-1 activates 5' splice sites followed by U-rich sequences and show that TIAR exerts a similar activity. They suggest that both proteins may autoregulate their expression at the level of splicing.     

The essential pre-mRNA splicing factor SF2 influences 5' splice site selection by activating proximal sites

SF2 is a 33 kd protein factor required for 5' splice site cleavage and lariat formation during pre-mRNA splicing in HeLa cell extracts. In addition to its essential role in constitutive splicing, SF2 can strongly influence 5' splice site selection. When pre-mRNAs containing multiple cis-competing 5' splice sites are spliced in vitro, high concentrations of purified SF2 promote the use of the 5' splice site closest to the 3' splice site. However, SF2 discriminates properly between authentic and cryptic splice sites. These effects of SF2 on splice site selection may reflect the cellular mechanisms that prevent exon skipping and ensure the accuracy of splicing. In addition, alterations in the concentration or activity of SF2, and of other general splicing factors, may serve to regulate alternative splicing in vivo.     

Imprecise excision of the Caenorhabditis elegans transposon Tc1 creates functional 5' splice sites.

Imprecise excision of the Caenorhabditis elegans transposon Tc1 from a specific site of insertion within the unc-54 myosin heavy chain gene generates either wild-type or partial phenotypic revertants. Wild-type revertants and one class of partial revertants contain insertions of four nucleotides in the unc-54 third exon (Tc1 "footprints"). Such revertants express large amounts of functional unc-54 myosin despite having what would appear to be frameshifting insertions in the unc-54 third exon. We demonstrate that these Tc1 footprints act as efficient 5' splice sites for removal of the unc-54 third intron. Splicing of these new 5' splice sites to the normal third intron splice acceptor removes the Tc1 footprint from the mature mRNA and restores the normal translational reading frame. Partial revertant unc-54(r661), which contains a single nucleotide substitution relative to the wild-type gene, is spliced similarly, except that the use of its new 5' splice site creates a frameshift in the mature mRNA rather than removing one. In all of these revertants, two alternative 5' splice sites are available to remove intron 3. We determined the relative efficiency with which each alternative 5' splice site is used by stabilizing frameshifted mRNAs with smg(-) genetic backgrounds. In all cases, the upstream member of the two alternative sites is used preferentially (> 75% utilization). This may reflect an inherent preference of the splicing machinery for the upstream member of two closely spaced 5' splice sites. Creation of new 5' splice sites may be a general characteristic of Tc1 insertion and excision events.     

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.     

Natural base-pairing interactions between 5' splice site and branch site sequences affect mammalian 5' splice site selection.

In the murine gene encoding the neuronal cell adhesion molecule (NCAM), the integrity of the 5' splice site of exon 18 (E18) is essential for regulation of alternative splicing. To further study the contribution of 5' splice site sequences, we used a simple NCAM pre-mRNA containing a portion of E18 fused to E19 and separated by a shortened intron. This RNA is spliced in vitro to produce five sets of lariat intermediates and products, the most abundant set displaying aberrant migration in acrylamide/urea gels. Base pairing interactions between positions +5 and +8 of the intron and positions -3 and -6 from the branch point were responsible for the faster migration of this set of lariat molecules. To test whether the duplex structure forms earlier and contributes to 5' splice site selection, we used NCAM substrates carrying the 5' splice sites of E17 and E18 in competition for the 3' splice site of E19. Mutations upstream of the major branch site improve E18/E19 splicing in NIH3T3 extracts, whereas compensatory mutations at positions +7 and +8 neutralize the effect of branch site mutations and curtail E18/E19 splicing. Our data suggest that duplex formation occurs early and interferes with the assembly of complexes initiated on the 5' splice site of NCAM E18. This novel type of intron interaction may exist in the introns of other mammalian pre-mRNAs.     

Alternative splicing regulation at tandem 3' splice sites

Alternative splicing (AS) constitutes a major mechanism creating protein diversity in humans. Previous bioinformatics studies based on expressed sequence tag and mRNA data have identified many AS events that are conserved between humans and mice. Of these events, ~25% are related to alternative choices of 3′ and 5′ splice sites. Surprisingly, half of all these events involve 3′ splice sites that are exactly 3 nt apart. These tandem 3′ splice sites result from the presence of the NAGNAG motif at the acceptor splice site, recently reported to be widely spread in the human genome. Although the NAGNAG motif is common in human genes, only a small subset of sites with this motif is confirmed to be involved in AS. We examined the NAGNAG motifs and observed specific features such as high sequence conservation of the motif, high conservation of ~30 bp at the intronic regions flanking the 3′ splice site and overabundance of cis-regulatory elements, which are characteristic of alternatively spliced tandem acceptor sites and can distinguish them from the constitutive sites in which the proximal NAG splice site is selected. Our findings imply that AS at tandem splice sites and constitutive splicing of the distal NAG are highly regulated.