Splice-site signals of mRNA precursors as revealed by computer search. Site-specific mutagenesis and thalassemia

Many eukaryotic genes are interrupted by introns, which are removed from mRNA precursors by the RNA splicing mechanism. Although nucleotide sequences around splice sites have considerable homology among various genes, exact splice-site signals are unknown. Previously, applying the computer searching method to primary nucleotide sequences of various genes, we studied what kinds of patterns are the necessary and sufficient sequences for recognition of the 5'-slice site (donor site). We proposed that four common patterns, AG/GTA, /GTAAGT, RG/GTGAG, and AG/GTXXGT, where R = A or G and X = A, T, G, or C, are often used as such signals. In the present paper, we examined a number of experimental results on site-specific mutagenesis around donor sites and on alpha- and beta-thalassemia in globin genes, and found that the above four common patterns could explain almost all of those results.     

Wobble splicing reveals the role of the branch point sequence-to-NAGNAG region in 3' tandem splice site selection

Alternative splicing involving the 3' tandem splice site NAGNAG sequence may play a role in the structure-function diversity of proteins. However, how 3' tandem splice site utilization is determined is not well understood. We previously demonstrated that 3' NAGNAG-based wobble splicing occurs mostly in a tissue- and developmental stage-independent manner. Bioinformatic analysis reveals that the nucleotide preceding the AG dinucleotide may influence 3' splice site utilization; this is also supported by an in vivo splicing assay. Moreover, we found that the intron sequence plays an important role in 3' splice site selection for NAGNAG wobble splicing. Mutations of the region between the branch site and the NAGNAG 3' splice site, indeed, affected the ratio of the distal/proximal AG selection. Finally, we found that single nucleotide polymorphisms around the NAGNAG motif could affect the splice site choice, which may lead to a change in mRNA patterns and influence protein function. We conclude that the NAGNAG motif and its upstream region to the branch point sequence are required for 3' tandem splice site selection.     

Mutation of putative branchpoint consensus sequences in plant introns reduces splicing efficiency

Intron lariat formation between the 5' end of an intron and a branchpoint adenosine is a fundamental aspect of the first step in animal and yeast nuclear pre-mRNA splicing. Despite similarities in intron sequence requirements and the components of splicing, differences exist between the splicing of plant and vertebrate introns. The identification of AU-rich sequences as major functional elements in plant introns and the demonstration that a branchpoint consensus sequence was not required for splicing have led to the suggestion that the transition from AU-rich intron to GC-rich exon is a major potential signal by which plant pre-mRNA splice sites are recognized. The role of putative branchpoint sequences as an internal signal in plant intron recognition/definition has been re-examined. Single nucleotide mutations in putative branchpoint adenosines contained within CUNAN sequences in four different plant introns all significantly reduced splicing efficiency. These results provide the most direct evidence to date for preferred branchpoint sequences being required for the efficient splicing of at least some plant introns in addition to the important role played by AU sequences in dicot intron recognition. The observed patterns of 3' splice site selection in the introns studied are consistent with the scanning model described for animal intron 3' splice site selection. It is suggested that, despite the clear importance of AU sequences for plant intron splicing, the fundamental processes of splice site selection and splicing in plants are similar to those in animals.     

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.     

U2AF binding selects for the high conservation of the C. elegans 3' splice site

Caenorhabditis elegans is unusual among animals in having a highly conserved octamer sequence at the 3' splice site: UUUU CAG/R. This sequence can bind to the essential heterodimeric splicing factor U2AF, with U2AF65 contacting the U tract and U2AF35 contacting the splice site itself (AG/R). Here we demonstrate a strong correspondence between binding to U2AF of RNA oligonucleotides with variant octamer sequences and the frequency with which such variations occur in splice sites. C. elegans U2AF has a strong preference for the octamer sequence and exerts much of the pressure for 3' splice sites to have the precise UUUUCAG/R sequence. At two positions the splice site has a very strong preference for U even though alternative bases can also bind tightly to U2AF, suggesting that evolution can select against sequences that may have a relatively modest reduction in binding. Although pyrimidines are frequently present at the first base in the exon, U2AF has a very strong bias against them, arguing there is a mechanism to compensate for weakened U2AF binding at this position. Finally, the C in the consensus sequence must remain adjacent to the AG/R rather than to the stretch of U's, suggesting this C is recognized by U2AF35.     

Identification of proteins that interact with exon sequences, splice sites, and the branchpoint sequence during each stage of spliceosome assembly.

We have carried out a systematic analysis of the proteins that interact with specific intron and exon sequences during each stage of mammalian spliceosome assembly. This was achieved by site-specifically labeling individual nucleotides within the 5' and 3' splice sites, the branchpoint sequence (BPS), or the exons with 32P and identifying UV-cross-linked proteins in the E, A, B, or C spliceosomal complex. Significantly, two members of the SR family of splicing factors, which are known to promote E-complex assembly, cross-link within exon sequences to a region approximately 25 nucleotides upstream from the 5' splice site. At the 5' splice site, cross-linking of the U5 small nuclear ribonucleoprotein particle protein, U5(200), was detected in both the B and C complexes. As observed in yeast cells, U5(200), also cross-links to intron/exon sequences at the 3' splice site in the C complex and may play a role in aligning the 5' and 3' exons for ligation. With label at the branch site, we detected three distinct proteins, designated BPS72,BpS70, and BPS56, which replace one another in the E, A, and C complexes. Another dynamic exchange was detected with pre-mRNA labeled at the AG dinucleotide of the 3' splice site. In this case, a protein, AG100,cross-links in the A complex and is replaced by another protein, AG75, in the C complex. The observation that these proteins are specifically associated with critical pre-mRNA sequence elements in functional complexes at different stages of spliceosome assembly implicates roles for these factors in key recognition events during the splicing pathway.     

Splice site consensus sequences are preferentially accessible to nucleases in isolated adenovirus RNA.

The conformation of RNA sequences spanning five 3' splice sites and two 5' splice sites in adenovirus mRNA was probed by partial digestion with single-strand specific nucleases. Although cleavage of nucleotides near both 3' and 5' splice sites was observed, most striking was the preferential digestion of sequences near the 3' splice site. At each 3' splice site a region of very strong cleavage is observed at low concentrations of enzyme near the splice site consensus sequence or the upstream branch point consensus sequence. Additional sites of moderately strong cutting near the branch point consensus sequence were observed in those sequences where the splice site was the preferred target. Since recognition of the 3' splice site and branch site appear to be early events in mRNA splicing these observations may indicate that the local conformation of the splice site sequences may play a direct or indirect role in enhancing the accessibility of sequences important for splicing.     

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.     

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.     

Transduction of c-src coding and intron sequences by a transformation-defective deletion mutant of Rous sarcoma virus.

The mechanism of cellular src (c-src) transduction by a transformation-defective deletion mutant, td109, of Rous sarcoma virus was studied by sequence analysis of the recombinational junctions in three td109-derived recovered sarcoma viruses (rASVs). Our results show that two rASVs have been generated by recombination between td109 and c-src at the region between exons 1 and 2 defined previously. Significant homology between td109 and c-src sequences was present at the sites of recombination. The viral and c-src sequence junction of the third rASV was formed by splicing a cryptic donor site at the 5' region of env of td109 to exon 1 of c-src. Various lengths of c-src internal intron 1 sequences were incorporated into all three rASV genomes, which resulted from activation of potential splice donor and acceptor sites. The incorporated intron 1 sequences were absent in the c-src mRNA, excluding its being the precursor for recombination with td109 and implying that initial recombinations most likely took place at the DNA level. A potential splice acceptor site within the incorporated intron 1 sequences in two rASVs was activated and was used for the src mRNA synthesis in infected cells. The normal env mRNA splice acceptor site was used for src mRNA synthesis for the third rASV.