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.     

An LKB1 AT-AC intron mutation causes Peutz-Jeghers syndrome via splicing at noncanonical cryptic splice sites

Peutz-Jeghers syndrome (PJS) is an autosomal dominant disorder associated with gastrointestinal polyposis and an increased cancer risk. PJS is caused by germline mutations in the tumor suppressor gene LKB1. One such mutation, IVS2+1A>G, alters the second intron 5' splice site, which has sequence features of a U12-type AT-AC intron. We report that in patients, LKB1 RNA splicing occurs from the mutated 5' splice site to several cryptic, noncanonical 3' splice sites immediately adjacent to the normal 3' splice site. In vitro splicing analysis demonstrates that this aberrant splicing is mediated by the U12-dependent spliceosome. The results indicate that the minor spliceosome can use a variety of 3' splice site sequences to pair to a given 5' splice site, albeit with tight constraints for maintaining the 3' splice site position. The unusual splicing defect associated with this PJS-causing mutation uncovers differences in splice-site recognition between the major and minor pre-mRNA splicing pathways.     

Identification,characterization and molecular phylogeny of U12-dependent introns in the Arabidopsis thaliana genome

U12-dependent introns are spliced by the minor U12-type spliceosome and occur in a variety of eukaryotic organisms, including Arabidopsis. In this study, a set of putative U12-dependent introns was compiled from a large collection of cDNA/EST- confirmed introns in the Arabidopsis thaliana genome by means of high-throughput bioinformatic analysis combined with manual scrutiny. A total of 165 U12-type introns were identified based upon stringent criteria. This number of sequences well exceeds the total number of U12-type introns previously reported for plants and allows a more thorough statistical analysis of U12-type signals. Of particular note is the discovery that the distance between the branch site adenosine and the acceptor site ranges from 10 to 39 nt, significantly longer than the previously postulated limit of 21 bp. Further analysis indicates that, in addition to the spacing constraint, the sequence context of the potential acceptor site may have an important role in 3′ splice site selection. Several alternative splicing events involving U12-type introns were also captured in this study, providing evidence that U12-dependent acceptor sites can also be recognized by the U2-type spliceosome. Furthermore, phylogenetic analysis suggests that both U12-type AT-AC and U12-type GT-AG introns occurred in Na⁺/H⁺ antiporters in a progenitor of animals and plants.     

Splicing of a divergent subclass of AT-AC introns requires the major spliceosomal snRNAs.

AT-AC introns constitute a minor class of eukaryotic pre-mRNA introns, characterized by 5''-AT and AC-3'' boundaries, in contrast to the 5''-GT and AG-3'' boundaries of the much more prevalent conventional introns. In addition to the AT-AC borders, most known AT-AC introns have highly conserved 5'' splice site and branch site sequence elements of 7-8 nt. Intron 6 of the nucleolar P120 gene and intron 2 of the SCN4A voltage-gated skeletal muscle sodium channel are AT-AC introns that have been shown recently to be processed via a unique splicing pathway involving several minor U snRNAs. Interestingly, intron 21 of the same SCN4A gene and the corresponding intron 25 of the SCN5A cardiac muscle sodium channel gene also have 5''-AT and AC-3'' boundaries, but they have divergent 5'' splice site and presumptive branch site sequences. Here, we report the accurate in vitro processing of these two divergent AT-AC introns and show that they belong to a functionally distinct subclass of AT-AC introns. Splicing of these introns does not require U12, U4atac, and U6atac snRNAs, but instead requires the major spliceosomal snRNAs U1, U2, U4, U5, and U6. Previous studies showed that G --> A mutation at the first position and G --> C mutation at the last position of a conventional yeast or mammalian GT-AG intron suppress each other in vivo, suggesting that the first and last bases participate in an essential non-Watson-Crick interaction. Our results show that such introns, hereafter termed AT-AC II introns, occur naturally and are spliced by a mechanism distinct from that responsible for processing of the apparently more common AT-AC I introns.     

The U11-48K protein contacts the 5' splice site of U12-type introns and the U11-59K protein

Little is currently known about proteins that make contact with the pre-mRNA in the U12-dependent spliceosome and thereby contribute to intron recognition. Using site-specific cross-linking, we detected an interaction between the U11-48K protein and U12-type 5' splice sites (5'ss). This interaction did not require branch point recognition and was sensitive to 5'ss mutations, suggesting that 48K interacts with the 5'ss during the first steps of prespliceosome assembly in a sequence-dependent manner. RNA interference-induced knockdown of 48K in HeLa cells led to reduced cell growth and the inhibition of U12-type splicing, as well as the activation of cryptic, U2-type splice sites, suggesting that 48K plays a critical role in U12-type intron recognition. 48K knockdown also led to reduced levels of U11/U12 di-snRNP, indicating that 48K contributes to the stability and/or formation of this complex. In addition to making contact with the 5'ss, 48K interacts with the U11-59K protein, a protein at the interface of the U11/U12 di-snRNP. These studies provide important insights into the protein-mediated recognition of the U12-type 5'ss, as well as functionally important interactions within the U11/U12 di-snRNP.     

Purine-rich enhancers function in the AT-AC pre-mRNA splicing pathway and do so independently of intact U1 snRNP.

A rare class of introns in higher eukaryotes is processed by the recently discovered AT-AC spliceosome. AT-AC introns are processed inefficiently in vitro, but the reaction is stimulated by exon-definition interactions involving binding of U1 snRNP to the 5'' splice site of the downstream conventional intron. We report that purine-rich exonic splicing enhancers also strongly stimulate sodium channel AT-AC splicing. Intact U2, U4, or U6 snRNAs are not required for enhancer function or for exon definition. Enhancer function is independent of U1 snRNP, showing that splicing stimulation by a downstream 5'' splice site and by an exonic enhancer differ mechanistically.     

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.     

Using Profiles Based on Nucleotide Hydrophobicity to Define Essential Regions for Splicing

The splice-site sequences of U2-type introns are highly degenerate, so many different sequences can function as U2-type splice sites. Using our new profiles based on hydrophobicity properties we pointed out specific properties for regions surrounding splice sites. We built a set T of flanking regions of genes with 1-3 introns from 21st and 22nd chromosomes extracted from GenBank to define positions having conserved properties, namely hydrophobicity, that are potentially essential for recognition by spliceosome.     

Alternative splicing and bioinformatic analysis of human U12-type introns

U12-type introns exist, albeit rarely, in a variety of multicellular organisms. Splicing of U12 intron-containing precursor mRNAs takes place in the U12-type spliceosome that is distinct from the major U2-type spliceosome. Due to incompatibility of these two spliceosomes, alternative splicing involving a U12-type intron may give rise to a relatively complicated impact on gene expression. We studied alternative U12-type intron splicing in an attempt to gain more mechanistic insights. First, we characterized mutually exclusive exon selection of the human JNK2 gene, which involves an unusual intron possessing the U12-type 5′ splice site and the U2-type 3′ splice site. We demonstrated that the long and evolutionary conserved polypyrimidine tract of this hybrid intron provides important signals for inclusion of its downstream alternative exon. In addition, we examined the effects of single nucleotide polymorphisms in the human WDFY1 U12-type intron on pre-mRNA splicing. These results provide mechanistic implications on splice-site selection of U12-type intron splicing. We finally discuss the potential effects of splicing of a U12-type intron with genetic defects or within a set of genes encoding RNA processing factors on global gene expression.     

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.     

Comparative analysis detects dependencies among the 5' splice-site positions

Human-mouse comparative genomics is an informative tool to assess sequence functionality as inferred from its conservation level. We used this approach to examine dependency among different positions of the 5' splice site. We compiled a data set of 50,493 homologous human-mouse internal exons and analyzed the frequency of changes among different positions of homologous human-mouse 5' splice-site pairs. We found mutual relationships between positions +4 and +5, +5 and +6, -2 and +5, and -1 and +5. We also demonstrated the association between the exonic and the intronic positions of the 5' splice site, in which a stronger interaction of U1 snRNA and the intronic portion of the 5' splice site compensates for weak interaction of U1 snRNA and the exonic portion of the 5' splice site, and vice versa. By using an ex vivo system that mimics the effect of mutation in the 5' splice site leading to familial dysautonomia, we demonstrated that U1 snRNA base-pairing with positions +6 and -1 is the only functional requirement for mRNA splicing of this 5' splice site. Our findings indicate the importance of U1 snRNA base-pairing to the exonic portion of the 5' splice site.     

RS domain-splicing signal interactions in splicing of U12-type and U2-type introns

Serine-arginine (SR) proteins are general metazoan splicing factors that contain an essential arginine/serine-rich (RS) domain. On typical U2-type introns, RS domains contact the branchpoint and 5' splice site to promote base-pairing with U small nuclear RNAs (snRNAs). Here we analyze the role of SR proteins in splicing of U12-type introns and in the second step of U2-type intron splicing. We show that RS domains contact the branchpoint and 5' splice site of a U12-type intron. On a U2-type intron, we find that the RS domain contacts the site of the U6 snRNA-5' splice site interaction during the first step of splicing and shifts to contact the site of the U5 snRNA-exon 1 interaction during the second step. Our results reveal alternative interactions between the RS domain and 5' splice site region that coincide with remodeling of the spliceosome between the two catalytic steps.     

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.     

An exonic splicing silencer downstream of the 3' splice site A2 is required for efficient human immunodeficiency virus type 1 replication

Alternative splicing of the human immunodeficiency virus type 1 (HIV-1) genomic mRNA produces more than 40 unique viral mRNA species, of which more than half remain incompletely spliced within an HIV-1-infected cell. Regulation of splicing at HIV-1 3' splice sites (3'ss) requires suboptimal polypyrimidine tracts, and positive or negative regulation of splicing occurs through binding of cellular factors to cis-acting splicing regulatory elements. We have previously shown that splicing at HIV-1 3'ss A2, which produces vpr mRNA and promotes inclusion of HIV-1 exon 3, is repressed by the hnRNP A/B-dependent exonic splicing silencer ESSV. Here we show that ESSV activity downstream of 3'ss A2 is localized to a 16-nucleotide element within HIV-1 exon 3. HIV-1 replication was reduced by 95% when ESSV was inactivated by mutagenesis. Reduced replication was concomitant with increased inclusion of exon 3 within spliced viral mRNA and decreased accumulation of unspliced viral mRNA, resulting in decreased cell-associated p55 Gag. Prolonged culture of ESSV mutant viruses resulted in two independent second-site reversions disrupting the splice sites that define exon 3, 3'ss A2 and 5' splice site D3. Either of these changes restored both HIV-1 replication and regulated viral splicing. Therefore, inhibition of HIV-1 3'ss A2 splicing is necessary for HIV-1 replication.     

Role of the 3' splice site in U12-dependent intron splicing

U12-dependent introns containing alterations of the 3' splice site AC dinucleotide or alterations in the spacing between the branch site and the 3' splice site were examined for their effects on splice site selection in vivo and in vitro. Using an intron with a 5' splice site AU dinucleotide, any nucleotide could serve as the 3'-terminal nucleotide, although a C residue was most active, while a U residue was least active. The penultimate A residue, by contrast, was essential for 3' splice site function. A branch site-to-3' splice site spacing of less than 10 or more than 20 nucleotides strongly activated alternative 3' splice sites. A strong preference for a spacing of about 12 nucleotides was observed. The combined in vivo and in vitro results suggest that the branch site is recognized in the absence of an active 3' splice site but that formation of the prespliceosomal complex A requires an active 3' splice site. Furthermore, the U12-type spliceosome appears to be unable to scan for a distal 3' splice site.     

The human 18S U11/U12 snRNP contains a set of novel proteins not found in the U2-dependent spliceosome

U11 and U12 snRNPs bind U12-type pre-mRNAs as a preformed di-snRNP complex, simultaneously recognizing the 5' splice site and branchpoint sequence. Thus, within the U12-type prespliceosome, U11/U12 components form a molecular bridge connecting both ends of the intron. We have affinity purified human 18S U11/U12 and 12S U11 snRNPs, and identified their protein components by using mass spectrometry. U11/U12 snRNPs lack all known U1 snRNP proteins but contain seven novel proteins (i.e., 65K, 59K, 48K, 35K, 31K, 25K, 20K) not found in the major spliceosome, four of which (59K, 48K, 35K, and 25K) are U11-associated. Thus, protein-protein and protein-RNA interactions contributing to 5' splice site recognition and/or intron bridging appear to differ significantly in the minor versus major prespliceosome. The majority of U11/U12 proteins are highly conserved in organisms known to contain U12-type introns. However, homologs of those associated with U11 were not detected in Drosophila melanogaster, consistent with the presence of a divergent U11 snRNP in flies. RNAi experiments revealed that several U11/U12 proteins are essential for cell viability, suggesting they play key roles in U12-type splicing. The presence of unique U11/U12 snRNP proteins in the U12-type spliceosome provides insight into potential evolutionary relationships between the major and minor spliceosome.     

A computational scan for U12-dependent introns in the human genome sequence

U12-dependent introns are found in small numbers in most eukaryotic genomes, but their scarcity makes accurate characterisation of their properties challenging. A computational search for U12-dependent introns was performed using the draft version of the human genome sequence. Human expressed sequences confirmed 404 U12-dependent introns within the human genome, a 6-fold increase over the total number of non-redundant U12-dependent introns previously identified in all genomes. Although most of these introns had AT-AC or GT-AG terminal dinucleotides, small numbers of introns with a surprising diversity of termini were found, suggesting that many of the non-canonical introns found in the human genome may be variants of U12-dependent introns and, thus, spliced by the minor spliceosome. Comparisons with U2-dependent introns revealed that the U12-dependent intron set lacks the ‘short intron’ peak characteristic of U2-dependent introns. Analysis of this U12-dependent intron set confirmed reports of a biased distribution of U12-dependent introns in the genome and allowed the identification of several alternative splicing events as well as a surprising number of apparent splicing errors. This new larger reference set of U12-dependent introns will serve as a resource for future studies of both the properties and evolution of the U12 spliceosome.