An SF1 affinity model to identify branch point sequences in human introns

Splicing factor 1 (SF1) binds to the branch point sequence (BPS) of mammalian introns and is believed to be important for the splicing of some, but not all, introns. To help identify BPSs, particularly those that depend on SF1, we generated a BPS profile model in which SF1 binding affinity data, validated by branch point mapping, were iteratively incorporated into computational models. We searched a data set of 117,499 human introns for best matches to the SF1 Affinity Model above a threshold, and counted the number of matches at each intronic position. After subtracting a background value, we found that 87.9% of remaining high-scoring matches identified were located in a region upstream of 3'-splice sites where BPSs are typically found. Since U2AF65 recognizes the polypyrimidine tract (PPT) and forms a cooperative RNA complex with SF1, we combined the SF1 model with a PPT model computed from high affinity binding sequences for U2AF65. The combined model, together with binding site location constraints, accurately identified introns bound by SF1 that are candidates for SF1-dependent splicing.     

Branch point identification and sequence requirements for intron splicing in Plasmodium falciparum

Splicing of mRNA is an ancient and evolutionarily conserved process in eukaryotic organisms, but intron-exon structures vary. Plasmodium falciparum has an extreme AT nucleotide bias (>80%), providing a unique opportunity to investigate how evolutionary forces have acted on intron structures. In this study, we developed an in vivo luciferase reporter splicing assay and employed it in combination with lariat isolation and sequencing to characterize 5' and 3' splicing requirements and experimentally determine the intron branch point in P. falciparum. This analysis indicates that P. falciparum mRNAs have canonical 5' and 3' splice sites. However, the 5' consensus motif is weakly conserved and tolerates nucleotide substitution, including the fifth nucleotide in the intron, which is more typically a G nucleotide in most eukaryotes. In comparison, the 3' splice site has a strong eukaryotic consensus sequence and adjacent polypyrimidine tract. In four different P. falciparum pre-mRNAs, multiple branch points per intron were detected, with some at U instead of the typical A residue. A weak branch point consensus was detected among 18 identified branch points. This analysis indicates that P. falciparum retains many consensus eukaryotic splice site features, despite having an extreme codon bias, and possesses flexibility in branch point nucleophilic attack.     

Multiple alignment using simulated annealing: branch point definition in human mRNA splicing.

A method for the simultaneous alignment of a very large number of sequences using simulated annealing is presented. The total running time of the algorithm does not depend explicitly on the number of sequences treated. The method has been used for the simultaneous alignment of 1462 human intron sequences upstream of the intron-exon boundary. The consensus sequence of the aligned set together with a calculation of the Shannon information clearly shows that several sequence motives are conserved: (i) a previously undetected guanosine rich region, (ii) the branch point and (iii) the polypyrimidine tract. The nucleotide frequencies at each position of the branch point consensus sequence qualitatively reproduce the frequencies of the experimentally determined branch points.     

Polypyrimidine tract sequences direct selection of alternative branch sites and influence protein binding.

IVS1, an intron derived from the rat fibronectin gene, is spliced inefficiently in vitro, involving the use of three alternative branch sites. Mutation of one branch point site, BP3, so as to increase complementarity to U2 snRNA resulted in exclusive use of that site and improved splicing efficiency, indicating that the wild type BP3 site is one determinant of poor IVS1 splicing. Deletions within the polypyrimidine tract had a variable effect on splicing efficiency and altered the pattern of branch site usage. Selection of each branch site was influenced negatively by purine substitutions ca. 20 nucleotides downstream. It is proposed that all three IVS1 branch sites are pyrimidine tract-dependent. Pyrimidine tract deletions also influenced the crosslinking of PTB (the polypyrimidine tract-binding protein), hnRNP C, and splicing factor U2AF65. All three proteins bound preferentially to distinct regions within the polypyrimidine tract and thus are candidates for mediating pyrimidine tract-dependent branch site selection. The findings indicate the complexity of the IVS1 polypyrimidine tract and suggest a crucial role for this region in modulating branch site selection and IVS1 splicing.     

The polypyrimidine tract binding protein (PTB) represses splicing of exon 6B from the beta-tropomyosin pre-mRNA by directly interfering with the binding of the U2AF65 subunit

Splicing of exon 6B from the beta-tropomyosin pre-mRNA is repressed in nonmuscle cells and myoblasts by a complex array of intronic elements surrounding the exon. In this study, we analyzed the proteins that mediate splicing repression of exon 6B through binding to the upstream element. We identified the polypyrimidine tract binding protein (PTB) as a component of complexes isolated from myoblasts that assemble onto the branch point region and the pyrimidine tract. In vitro splicing assays and PTB knockdown experiments by RNA interference demonstrated that PTB acts as a repressor of splicing of exon 6B. Using psoralen experiments, we showed that PTB acts at an early stage of spliceosome assembly by preventing the binding of U2 snRNA on the branch point. Using UV cross-linking and immunoprecipitation experiments with site-specific labeled RNA in PTB-depleted nuclear extracts, we found that the decrease in PTB was correlated with an increase in U2AF65. In addition, competition experiments showed that PTB is able to displace the binding of U2AF65 on the polypyrimidine tract. Our results strongly support a model whereby PTB competes with U2AF65 for binding to the polypyrimidine tract.     

mRNA processing in Dictyostelium: sequence requirements for termination and splicing

Criteria for the identification of termination regions in Dictyostelium discoideum genes have been established and the sequence requirements for termination in 33 genes have been analyzed. A canonical hexamer signal AATAAA was present 15-30 nucleotides upstream of the cleavage site, usually a TA, and was embedded in a particularly A-rich environment. T- or GT-rich downstream elements characteristic of animal cells could not be identified. In a sample of 102 introns we have established the consensus AG/GTAAGT and ATAG/ for the 5' and 3' splice sites, respectively. Most introns are 75-150 nucleotides long and the A+T content is high (90%). A putative branch point was identified in half of the introns 20-60 nucleotides upstream of the 3' splice site and the consensus TACTAAY was derived. A polypyrimidine tract required for branching in vertebrates was not identified, but weak preference for pyrimidine was found 10-45 nucleotides upstream of the 3' splice site.     

A class of human exons with predicted distant branch points revealed by analysis of AG dinucleotide exclusion zones

Background

The three consensus elements at the 3' end of human introns - the branch point sequence, the polypyrimidine tract, and the 3' splice site AG dinucleotide - are usually closely spaced within the final 40 nucleotides of the intron. However, the branch point sequence and polypyrimidine tract of a few known alternatively spliced exons lie up to 400 nucleotides upstream of the 3' splice site. The extended regions between the distant branch points (dBPs) and their 3' splice site are marked by the absence of other AG dinucleotides. In many cases alternative splicing regulatory elements are located within this region.

Results

We have applied a simple algorithm, based on AG dinucleotide exclusion zones (AGEZ), to a large data set of verified human exons. We found a substantial number of exons with large AGEZs, which represent candidate dBP exons. We verified the importance of the predicted dBPs for splicing of some of these exons. This group of exons exhibits a higher than average prevalence of observed alternative splicing, and many of the exons are in genes with some human disease association.

Conclusion

The group of identified probable dBP exons are interesting first because they are likely to be alternatively spliced. Second, they are expected to be vulnerable to mutations within the entire extended AGEZ. Disruption of splicing of such exons, for example by mutations that lead to insertion of a new AG dinucleotide between the dBP and 3' splice site, could be readily understood even though the causative mutation might be remote from the conventional locations of splice site sequences.     

A Potential Role for U2AF-SAP 155 Interactions in Recruiting U2 snRNP to the Branch Site

Base pairing between U2 snRNA and the branchpoint sequence (BPS) is essential for pre-mRNA splicing. Because the metazoan BPS is short and highly degenerate, this interaction alone is insufficient for specific binding of U2 snRNP. The splicing factor U2AF binds to the pyrimidine tract at the 3′ splice site in the earliest spliceosomal complex, E, and is essential for U2 snRNP binding in the spliceosomal complex A. We show that the U2 snRNP protein SAP 155 UV cross-links to pre-mRNA on both sides of the BPS in the A complex. SAP 155’s downstream cross-linking site is immediately adjacent to the U2AF binding site, and the two proteins interact directly in protein-protein interaction assays. Using UV cross-linking, together with functional analyses of pre-mRNAs containing duplicated BPSs, we show a direct correlation between BPS selection and UV cross-linking of SAP 155 on both sides of the BPS. Together, our data are consistent with a model in which U2AF binds to the pyrimidine tract in the E complex and then interacts with SAP 155 to recruit U2 snRNP to the BPS.     

Atypical 5' splice sites cause CFTR exon 9 to be vulnerable to skipping

The molecular basis of the skipping of constitutive exons in many messenger RNAs is not fully understood. A well-studied example is exon 9 of the human cystic fibrosis transmembrane conductance regulator gene (CFTR), in which an abbreviated polypyrimidine tract between the branch point A and the 3' splice site is associated with increased exon skipping and disease. However, many exons, both in CFTR and in other genes and have short polypyrimidine tracts in their 3' splice sites, yet they are not skipped. Inspection of the 5' splice sites immediately up- and downstream of exon 9 revealed deviations from consensus sequence, so we hypothesized that this exon may be inherently vulnerable to skipping. To test this idea, we constructed a CFTR minigene and replicated exon 9 skipping associated with the length of the polypyrimidine tract upstream of exon 9. We then mutated the flanking 5' splice sites and determined the effect on exon skipping. Conversion of the upstream 5' splice site to consensus by replacing a pyrimidine at position +3 with a purine resulted in increased exon skipping. In contrast, conversion of the downstream 5' splice site to consensus by insertion of an adenine at position +4 resulted in a substantial reduction in exon 9 skipping, regardless of whether the upstream 5' splice site was consensus or not. These results suggested that the native downstream 5' splice site plays an important role in CFTR exon 9 skipping, a hypothesis that was supported by data from sheep and mouse genomes. Although CFTR exon 9 in sheep is preceded by a long polypyrimidine tract (Y(14)), it skips exon 9 in vivo and has a nonconsensus downstream 5' splice site identical to that in humans. On the other hand, CFTR exon 9 in mice is preceded by a short polypyrimidine tract (Y(5)) but is not skipped in vivo. Its downstream 5' splice site differs from that in humans by a 2-nt insertion, which, when introduced into the human CFTR minigene, abolished exon 9 skipping. Taken together, these observations place renewed emphasis on deviations at 5' splice sites in nucleotides other than the invariant GT, particularly when such changes are found in conjunction with other altered splicing sequences, such as a shortened polypyrimidine tract. Thus, careful inspection of entire 5' splice sites may identify constitutive exons that are vulnerable to skipping.     

Both the polypyrimidine tract and the 3' splice site function prior to the first step of splicing in fission yeast.

While it is known that several trans -acting splicing factors are highly conserved between Schizosaccharomyces pombe and mammals, the roles of cis -acting signals have received comparatively little attention. In Saccharomyces cerevisiae, sequences downstream from the branch point are not required prior to the first transesterification reaction, whereas in mammals the polypyrimidine tract and, in some introns, the 3' AG dinucleotide are critical for initial recognition of an intron. We have investigated the contribution of these two sequence elements to splicing in S.pombe. To determine the stage at which the polypyrimidine tract functions, we analyzed the second intron of the cdc2 gene (cdc 2-Int2), in which pyrimidines span the entire interval between the branch point and 3' splice site. Our data indicate that substitution of a polypurine tract results in accumulation of linear pre-mRNA, while expanding the polypyrimidine tract enhances splicing efficiency, as in mammals. To examine the role of the AG dinucleotide in cdc 2-Int2 splicing, we mutated the 3' splice junction in both the wild-type and pyrimidine tract variant RNAs. These changes block the first transesterification reaction, as in a subset of mammalian introns. However, in contrast to the situation in mammals, we were unable to rescue the first step of splicing in a 3' splice site mutant by expanding the polypyrimidine tract. Mutating the terminal G in the third intron of the nda 3 gene (nda 3-Int3) also blocks the first transesterification reaction, suggesting that early recognition of the 3' splice site is a general property of fission yeast introns. Counter to earlier work with an artificial intron, it is not possible to restore the first step of splicing in cdc 2-Int2 and nda 3-Int3 3' splice site mutants by introducing compensatory changes in U1 snRNA. These results highlight the diversity and probable redundancy of mechanisms for identifying the 3' ends of introns.     

Mutational analysis of the U12-dependent branch site consensus sequence

Highly conserved sequences at the 5′ splice site and branch site of U12-dependent introns are important determinants for splicing by U12-dependent spliceosomes. This study investigates the in vivo splicing phenotypes of mutations in the branch site consensus sequence of the U12-dependent intron F from a human NOL1 (P120) minigene. Intron F contains a fully consensus branch site sequence (UUCCUUAAC). Mutations at each position were analyzed for their effects on U12-dependent splicing in vivo. Mutations at most positions resulted in a significant reduction of correct U12-dependent splicing. Defects observed included increased unspliced RNA levels, the activation of cryptic U2-dependent 5′ and 3′ splice sites, and the activation of cryptic U12-dependent branch/3′ splice sites. A strong correlation was observed between the predicted thermodynamic stability of the branch site: U12 snRNA interaction and correct U12-dependent splicing. The lack of a polypyrimidine tract between the branch site and 3′ splice site of U12-dependent introns and the observed reliance on base-pairing interactions for correct U12-dependent splicing emphasize the importance of RNA/RNA interactions during U12-dependent intron recognition and proper splice site selection.     

Alternative modes of binding by U2AF65 at the polypyrimidine tract

During initial recognition of an intron in pre-mRNA, the 3' end of the intron is bound by essential splicing factors. Notably, the consensus RNA sequences bound by these proteins are highly degenerate in humans. This raises the question of 3' splicing factor function in introns lacking canonical binding sites. Investigating the introns of the model organism Neurospora crassa revealed a different organization at the 3' end of the intron compared to most eukaryotic organisms. The predicted branch point sequences of Neurospora introns are much closer to the 3' splice site compared to those in human introns. In addition, Neurospora introns lack the canonical polypyrimidine tract found at the end of introns in most eukaryotic organisms. The large subunit of the U2 snRNP associated factor (U2AF65), which is essential for splicing of human introns and specifically recognizes the polypyrimidine tract, is also present in Neurospora. We show that Neurospora U2AF65 binds RNA with low affinity and specificity, apparently evolving with its disappearing binding site. The arginine/serine rich domain at the N-terminus of Neurospora U2AF65 regulates its RNA binding. We find that this regulated binding can be recapitulated in human U2AF65 which has been mutated to decrease both affinity and overall charge. Finally, we show that the addition of the small U2AF subunit (U2AF35) to U2AF65 with weakened RNA binding affinity significantly enhances the affinity of the resulting U2AF heterodimer.     

A conditional role of U2AF in splicing of introns with unconventional polypyrimidine tracts

Recognition of polypyrimidine (Py) tracts typically present between the branch point and the 3' splice site by the large subunit of the essential splicing factor U2AF is a key early step in pre-mRNA splicing. Diverse intronic sequence arrangements exist, however, including 3' splice sites lacking recognizable Py tracts, which raises the question of how general the requirement for U2AF is for various intron architectures. Our analysis of fission yeast introns in vivo has unexpectedly revealed that whereas introns lacking Py tracts altogether remain dependent on both subunits of U2AF, introns with long Py tracts, unconventionally positioned upstream of branch points, are unaffected by U2AF inactivation. Nevertheless, mutation of these Py tracts causes strong dependence on the large subunit U2AF59. We also find that Py tract diversity influences the requirement for the conserved C-terminal domain of U2AF59 (RNA recognition motif 3), which has been implicated in protein-protein interactions with other splicing factors. Together, these results suggest that in addition to Py tract binding by U2AF, supplementary mechanisms of U2AF recruitment and 3' splice site identification exist to accommodate diverse intron architectures, which have gone unappreciated in biochemical studies of model pre-mRNAs.     

Polypyrimidine tract binding protein inhibits IgM pre-mRNA splicing by diverting U2 snRNA base-pairing away from the branch point

The mouse immunoglobulin (IgM) pre-mRNA contains a splicing inhibitor that bears multiple binding sites for the splicing repressor polypyrimidine tract binding protein (PTB). Here we show that the inhibitor directs assembly of an ATP-dependent complex that contains PTB and U1 and U2 small nuclear RNAs (snRNAs). Unexpectedly, although U2 snRNA is present in the inhibitor complex, it is not base-paired to the branch point. We present evidence that inhibitor-bound PTB contacts U2 snRNA to promote base-pairing to an adjacent branch point–like sequence within the inhibitor, thereby preventing the U2 snRNA–branch point interaction and resulting in splicing repression. Our studies reveal a novel mechanism by which PTB represses splicing.     

Functional analysis of the polypyrimidine tract in pre-mRNA splicing.

The polypyrimidine tract is one of the important cis-acting sequence elements directing intron removal in pre-mRNA splicing. Progressive deletions of the polypyrimidine tract have been found to abolish correct lariat formation, spliceosome assembly and splicing. In addition, the polypyrimidine tract can alter 3'-splice site selection by promoting alternative branch site selection. However, there appears to be great flexibility in the specific sequence of a given tract. Not only the optimal composition of the polypyrimidine tract, but also the role of the tract in introns with no apparent polypyrimidine tracts or where changes in the tract are apparently harmless are uncertain. Accordingly, we have designed a series of cis-competition splicing constructs to test the functional competitive efficiency of a variety of systematically mutated polypyrimidine tracts. An RT/PCR assay was used to detect spliced product formation as a result of differential branch point selection dependent on direct competition between two opposing polypyrimidine tracts. We found that pyrimidine tracts containing 11 continuous uridines are the strongest pyrimidine tracts. In such cases, the position of the uridine stretch between the branch point and 3'-splice site AG is unimportant. In contrast, decreasing the continuous uridine stretch to five or six residues requires that the tract be located immediately adjacent to the AG for optimal competitive efficiency. The block to splicing with decreasing polypyrimidine tract strength is primarily prior to the first step of splicing. While lengthy continuous uridine tracts are the most competitive, tracts with decreased numbers of consecutive uridines and even tracts with alternating purine/pyrimidine residues can still function to promote branch point selection, but are far less effective competitors in 3'-splice site selection assays.     

Sequential recognition of the pre-mRNA branch point by U2AF65 and a novel spliceosome-associated 28-kDa protein.

Splicing of pre-mRNAs occurs via a lariat intermediate in which an intronic adenosine, embedded within a branch point sequence, forms a 2'',5''-phosphodiester bond (RNA branch) with the 5'' end of the intron. How the branch point is recognized and activated remains largely unknown. Using site-specific photochemical cross-linking, we have identified two proteins that specifically interact with the branch point during the splicing reaction. U2AF65, an essential splicing factor that binds to the adjacent polypyrimidine tract, crosslinks to the branch point at the earliest stage of spliceosome formation in an ATP-independent manner. A novel 28-kDa protein, which is a constituent of the mature spliceosome, contacts the branch point after the first catalytic step. Our results indicate that the branch point is sequentially recognized by distinct splicing factors in the course of the splicing reaction.     

Dissection of splicing regulation at an endogenous locus by zinc-finger nuclease-mediated gene editing

Sequences governing RNA splicing are difficult to study in situ due to the great difficulty of traditional targeted mutagenesis. Zinc-finger nuclease (ZFN) technology allows for the rapid and efficient introduction of site-specific mutations into mammalian chromosomes. Using a ZFN pair along with a donor plasmid to manipulate the outcomes of DNA repair, we introduced several discrete, targeted mutations into the fourth intron of the endogenous BAX gene in Chinese hamster ovary cells. Putative lariat branch points, the polypyrimidine tract, and the splice acceptor site were targeted. We recovered numerous otherwise isogenic clones carrying the intended mutations and analyzed the effect of each on BAX pre-mRNA splicing. Mutation of one of three possible branch points, the polypyrimidine tract, and the splice acceptor site all caused exclusion of exon five from BAX mRNA. Interestingly, these exon-skipping mutations allowed usage of cryptic splice acceptor sites within intron four. These data demonstrate that ZFN-mediated gene editing is a highly effective tool for dissection of pre-mRNA splicing regulatory sequences in their endogenous context.     

Interaction of hnRNP A1 with snRNPs and pre-mRNAs: evidence for a possible role of A1 RNA annealing activity in the first steps of spliceosome assembly.

The in vitro interaction of recombinant hnRNP A1 with purified snRNPs and with pre-mRNAs was investigated. We show that protein A1 can stably bind U2 and U4 snRNP but not U1. Oligo-RNAse H cleavage of U2 nucleotides involved in base pairing with the branch site, totally eliminates the A1-U2 interaction. RNase T1 protection and immunoprecipitation experiments demonstrate that recombinant protein A1 specifically binds the 3'-end regions of both beta-globin and Ad-2 introns. However, while on the beta-globin intron only binding to the polypyrimidine tract was observed, on the Ad-2 intron a 32 nt fragment encompassing the branch point and the AG splice-site dinucleotide was bound and protected. Such protection was drastically reduced in the presence of U2 snRNP. Altogether these results indicate that protein A1 can establish a different pattern of association with different pre-mRNAs and support the hypothesis that this protein could play a role in the annealing of U2 to the branch site and hence in the early events of pre-splicing complex assembly.     

A catalogue of splice junction and putative branch point sequences from plant introns.

Splice junction and possible branch point sequences have been collected from 177 plant introns. Consensus sequences for the 5' and 3' splice junctions and for possible branch points have been derived. The splice junction consensus sequences were virtually identical to those of animal introns except that the polypyrimidine stretch at the 3' splice junction was less pronounced in the plant introns. A search for possible branch points with sequences related to the yeast, vertebrate and fungal consensus sequences revealed a similar sequence in plant introns.