Identification of motifs that function in the splicing of non-canonical introns

Background  

While the current model of pre-mRNA splicing is based on the recognition of four canonical intronic motifs (5' splice site, branchpoint sequence, polypyrimidine (PY) tract and 3' splice site), it is becoming increasingly clear that splicing is regulated by both canonical and non-canonical splicing signals located in the RNA sequence of introns and exons that act to recruit the spliceosome and associated splicing factors. The diversity of human intronic sequences suggests the existence of novel recognition pathways for non-canonical introns. This study addresses the recognition and splicing of human introns that lack a canonical PY tract. The PY tract is a uridine-rich region at the 3' end of introns that acts as a binding site for U2AF65, a key factor in splicing machinery recruitment.     

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.     

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.     

Four exons of the serotonin receptor 4 gene are associated with multiple distant branch points

Splicing of vertebrate introns involves recognition of three consensus elements at the 3′ end. The branch point (BP) and polypyrimidine tract (PPT) are usually located within 40 nucleotides (nt) of the 3′ splice site (3′ ss), AG, but can be much more distant. A characteristic of the region between distant BPs (dBPs) and the 3′ ss is the absence of intervening AG dinucleotides, leading to its designation as the “AG exclusion zone” (AGEZ). The human HTR4 gene, which encodes serotonin receptor 4 and has been associated with schizophrenia, bipolar disease, and gastrointestinal disorders, has four exons with extensive AGEZs. We have mapped the BPs for HTR4 exons 3, 4, 5, and g generated by in vitro splicing, and validated them by mutagenesis in exon-trapping vectors. All exons used dBPs up to 273 nt upstream of the exon. Strikingly, exons 4 and 5 used combinations of both distant and conventionally located BPs, suggesting that successful splicing of these exons can occur by distinct pathways. Our results emphasize the importance for single nucleotide polymorphism resequencing projects to take account of potential dBPs, as the extended AGEZs are vulnerable to mutations that could affect splicing itself or regulation of alternative splicing.     

ATP-dependent interaction of yeast U5 snRNA loop 1 with the 5' splice site.

Pre-messenger RNA splicing is a two-step process by which introns are removed and exons joined together. In yeast, the U5 snRNA loop 1 interacts with the 5' exon before the first step of splicing and with the 5' and 3' exons before the second step. In vitro studies revealed that yeast U5 loop 1 is not required for the first step of splicing but is essential for holding the 5' and 3' exons for ligation during the second step. It is critical, therefore, that loop 1 contacts the 5' exon before the first step of splicing to hold this exon following cleavage from the pre-mRNA. At present it is not known how U5 loop 1 is positioned on the 5' exon prior to the first step of splicing. To address this question, we have used site-specific photoactivated crosslinking in yeast spliceosomes to investigate the interaction of U5 loop 1 with the pre-mRNA prior to the first step of splicing. We have found that the highly conserved uridines in loop 1 make ATP-dependent contacts with an approximately 8-nt region at the 5' splice site that includes the invariant GU. These interactions are dependent on functional U2 and U6 snRNAs. Our results support a model where U5 snRNA loop 1 interacts with the 5' exon in two steps during its targeting to the 5' splice site.     

Genomic mRNA profiling reveals compensatory mechanisms for the requirement of the essential splicing factor U2AF

The large subunit of the U2 auxiliary factor (U2AF) recognizes the polypyrimidine tract (Py-tract) located adjacent to the 3' splice site to facilitate U2 snRNP recruitment. While U2AF is considered essential for pre-mRNA splicing, its requirement for splicing on a genome-wide level has not been analyzed. Using Solexa sequencing, we performed mRNA profiling for splicing in the Schizosaccharomyces pombe U2AF(59) (prp2.1) temperature-sensitive mutant. Surprisingly, our analysis revealed that introns show a range of splicing defects in the mutant strain. While U2AF(59) inactivation (nonpermissive) conditions inhibit splicing of some introns, others are spliced apparently normally. Bioinformatics analysis indicated that U2AF(59)-insensitive introns have stronger 5' splice sites and higher A/U content. Most importantly, features that contribute to U2AF(59) insensitivity of an intron unexpectedly reside in its 5'-most 30 nucleotides. These include the 5' splice site, a guanosine at position 7, and the 5' splice site-to-branch point sequence context. A differential requirement (similar to U2AF(59)) for introns may also apply to other general splicing factors (e.g., prp10). Our combined results indicate that U2AF insensitivity is a common phenomenon and that varied intron features support the existence of unrecognized aspects of spliceosome assembly.     

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.     

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.     

Intron definition in splicing of small Drosophila introns.

Approximately half of the introns in Drosophila melanogaster are too small to function in a vertebrate and often lack the pyrimidine tract associated with vertebrate 3' splice sites. Here, we report the splicing and spliceosome assembly properties of two such introns: one with a pyrimidine-poor 3' splice site and one with a pyrimidine-rich 3' splice site. The pyrimidine-poor intron was absolutely dependent on its small size for in vivo and in vitro splicing and assembly. As such, it had properties reminiscent of those of yeast introns. The pyrimidine-rich intron had properties intermediate between those of yeasts and vertebrates. This 3' splice site directed assembly of ATP-dependent complexes when present as either an intron or exon and supported low levels of in vivo splicing of a moderate-length intron. We propose that splice sites can be recognized as pairs across either exons or introns, depending on which distance is shorter, and that a pyrimidine-rich region upstream of the 3' splice site facilitates the exon mode.     

Introns and splicing elements of five diverse fungi

Genomic sequences and expressed sequence tag data for a diverse group of fungi (Saccharomyces cerevisiae, Schizosaccharomyces pombe, Aspergillus nidulans, Neurospora crassa, and Cryptococcus neoformans) provided the opportunity to accurately characterize conserved intronic elements. An examination of large intron data sets revealed that fungal introns in general are short, that 98% or more of them belong to the canonical splice site (ss) class (5'GU...AG3'), and that they have polypyrimidine tracts predominantly in the region between the 5' ss and the branch point. Information content is high in the 5' ss, branch site, and 3' ss regions of the introns but low in the exon regions adjacent to the introns in the fungi examined. The two yeasts have broader intron length ranges and correspondingly higher intron information content than the other fungi. Generally, as intron length increases in the fungi, so does intron information content. Homologs of U2AF spliceosomal proteins were found in all species except for S. cerevisiae, suggesting a nonconventional role for U2AF in the absence of canonical polypyrimidine tracts in the majority of introns. Our observations imply that splicing in fungi may be different from that in vertebrates and may require additional proteins that interact with polypyrimidine tracts upstream of the branch point. Theoretical protein homologs for Nam8p and TIA-1, two proteins that require U-rich regions upstream of the branch point to function, were found. There appear to be sufficient differences between S. cerevisiae and S. pombe introns and the introns of two filamentous members of the Ascomycota and one member of the Basidiomycota to warrant the development of new model organisms for studying the splicing mechanisms of fungi.     

Exon definition may facilitate splice site selection in RNAs with multiple exons.

Interactions at the 3' end of the intron initiate spliceosome assembly and splice site selection in vertebrate pre-mRNAs. Multiple factors, including U1 small nuclear ribonucleoproteins (snRNPs), are involved in initial recognition at the 3' end of the intron. Experiments were designed to test the possibility that U1 snRNP interaction at the 3' end of the intron during early assembly functions to recognize and define the downstream exon and its resident 5' splice site. Splicing precursor RNAs constructed to have elongated second exons lacking 5' splice sites were deficient in spliceosome assembly and splicing activity in vitro. Similar substrates including a 5' splice site at the end of exon 2 assembled and spliced normally as long as the second exon was less than 300 nucleotides long. U2 snRNPs were required for protection of the 5' splice site terminating exon 2, suggesting direct communication during early assembly between factors binding the 3' and 5' splice sites bordering an exon. We suggest that exons are recognized and defined as units during early assembly by binding of factors to the 3' end of the intron, followed by a search for a downstream 5' splice site. In this view, only the presence of both a 3' and a 5' splice site in the correct orientation and within 300 nucleotides of one another will stable exon complexes be formed. Concerted recognition of exons may help explain the 300-nucleotide-length maximum of vertebrate internal exons, the mechanism whereby the splicing machinery ignores cryptic sites within introns, the mechanism whereby exon skipping is normally avoided, and the phenotypes of 5' splice site mutations that inhibit splicing of neighboring introns.     

Kinetic role for mammalian SF1/BBP in spliceosome assembly and function after polypyrimidine tract recognition by U2AF

Two sequences important for pre-mRNA splicing precede the 3' end of introns in higher eukaryotes, the branch point (BP) and the polypyrimidine (Py) tract. Initial recognition of these signals involves cooperative binding of the splicing factor SF1/mammalian branch point binding protein (mBBP) to the BP and of U2AF(65) to the Py tract. Both factors are required for recruitment of the U2 small nuclear ribonucleoprotein particle (U2 snRNP) to the BP in reactions reconstituted from purified components. In contrast, extensive depletion of ST1/BBP in Saccharomyces cerevisiae does not compromise spliceosome assembly or splicing significantly. As BP sequences are less conserved in mammals, these discrepancies could reflect more stringent requirements for SF1/BBP in this system. We report here that extensive depletion of SF1/mBBP from nuclear extracts of HeLa cells results in only modest reduction of their activity in spliceosome assembly and splicing. Some of these effects reflect differences in the kinetics of U2 snRNP binding. Although U2AF(65) binding was reduced in the depleted extracts, the defects caused by SF1/mBBP depletion could not be fully restored by an increase in occupancy of the Py tract by exogenously added U2AF(65), arguing for a role of SF1/mBBP in U2 snRNP recruitment distinct from promoting U2AF(65) binding.     

Branchpoint and polypyrimidine tract mutations mediating the loss and partial recovery of the Moloney murine sarcoma virus MuSVts110 thermosensitive splicing phenotype.

Balanced splicing of retroviral RNAs is mediated by weak signals at the 3' splice site (ss) acting in concert with other cis elements. Moloney murine sarcoma virus MuSVts110 shows a similar balance between unspliced and spliced RNAs, differing only in that the splicing of its RNA is, in addition, growth temperature sensitive. We have generated N-nitroso-N-methylurea (NMU)-treated MuSVts110 revertants in which splicing was virtually complete at all temperatures and have investigated the molecular basis of this reversion on the assumption that the findings would reveal cis-acting elements controlling MuSVts110 splicing thermosensitivity. In a representative revertant (NMU-20), we found that complete splicing was conferred by a G-to-A substitution generating a consensus branchpoint (BP) signal (-CCCUGGC- to -CCCUGAC- [termed G(-25)A]) at -25 relative to the 3' ss. Weakening this BP to -CCCGAC- [G(-25)A,U(-27)C] moderately reduced splicing at the permissive temperature and sharply inhibited splicing at the originally nonpermissive temperature, arguing that MuSVts110 splicing thermosensitivity depends on a suboptimal BP-U2 small nuclear RNA interaction. This conclusion was supported by results indicating that lengthening the short MuSVts110 polypyrimidine tract and altering its uridine content doubled splicing efficiency at permissive temperatures and nearly abrogated splicing thermosensitivity. In vitro splicing experiments showed that MuSVts110 G(-25)A RNA intermediates were far more efficiently ligated than RNAs carrying the wild-type BP, the G(-25)A,U (-27)C BP, or the extended polypyrimidine tract. The efficiency of ligation in vitro roughly paralleled splicing efficiency in vivo [G(-25)A BP > extended polypyrimidine tract > G(-25)A,U(-27)C BP > wild-type BP]. These results suggest that MuSVts110 RNA splicing is balanced by cis elements similar to those operating in other retroviruses and, in addition, that its splicing thermosensitivity is a response to the presence of multiple suboptimal splicing signals.     

Alternative splicing of U12-dependent introns in vivo responds to purine-rich enhancers.

Alternative splicing increases the coding capacity of genes through the production of multiple protein isoforms by the conditional use of splice sites and exons. Many alternative splice sites are regulated by the presence of purine-rich splicing enhancer elements (ESEs) located in the downstream exon. Although the role of ESEs in alternative splicing of the major class U2-dependent introns is well established, no alternatively spliced minor class U12-dependent introns have so far been described. Although in vitro studies have shown that ESEs can stimulate splicing of individual U12-dependent introns, there is no direct evidence that the U12-dependent splicing system can respond to ESEs in vivo. To investigate the ability of U12-dependent introns to use alternative splice sites and to respond to ESEs in an in vivo context, we have constructed two sets of artificial minigenes with alternative splicing pathways and evaluated the effects of ESEs on their alternative splicing patterns. In minigenes with alternative U12-dependent 3' splice sites, a purine-rich ESE promotes splicing to the immediately upstream 3' splice site. As a control, a mutant ESE has no stimulatory effect. In minigene constructs with two adjacent U12-dependent introns, the predominant in vivo splicing pattern results in the skipping of the internal exon. Insertion of a purine-rich ESE into the internal exon promotes the inclusion of the internal exon. These results show that U12-dependent introns can participate in alternative splicing pathways and that U12-dependent splice sites can respond to enhancer elements in vivo.     

Utilization of Splicing Elements and Polyadenylation Signal Elements in the Coupling of Polyadenylation and Last-Intron Removal

Polyadenylation (PA) is the process by which the 3' ends of most mammalian mRNAs are formed. In nature, PA is highly coordinated, or coupled, with splicing. In mammalian systems, the most compelling mechanistic model for coupling arises from data supporting exon definition (2, 34, 37). We have examined the roles of individual functional components of splicing and PA signals in the coupling process by using an in vitro splicing and PA reaction with a synthetic pre-mRNA substrate containing an adenovirus splicing cassette and the simian virus 40 late PA signal. The effects of individually mutating splicing elements and PA elements in this substrate were determined. We found that mutation of the polypyrimidine tract and the 3' splice site significantly reduced PA efficiency and that mutation of the AAUAAA and the downstream elements of the PA signal decreased splicing efficiency, suggesting that these elements are the most significant for the coupling of splicing and PA. Although mutation of the upstream elements (USEs) of the PA signal dramatically decreased PA, splicing was only modestly affected, suggesting that USEs modestly affect coupling. Mutation of the 5' splice site in the presence of a viable polypyrimidine tract and the 3' splice site had no effect on PA, suggesting no effect of this element on coupling. However, our data also suggest that a site for U1 snRNP binding (e.g., a 5' splice site) within the last exon can negatively effect both PA and splicing; hence, a 5' splice site-like sequence in this position appears to be a modulator of coupling. In addition, we show that the RNA-protein complex formed to define an exon may inhibit processing if the definition of an adjacent exon fails. This finding indicates a mechanism for monitoring the appropriate definition of exons and for allowing only pre-mRNAs with successfully defined exons to be processed.     

Ser/Arg-rich protein-mediated communication between U1 and U2 small nuclear ribonucleoprotein particles

Previous work demonstrated that U1 small nuclear ribonucleoprotein particle (snRNP), bound to a downstream 5' splice site, can positively influence utilization of an upstream 3' splice site via exon definition in both trans- and cis-splicing systems. Although exon definition results in the enhancement of splicing of an upstream intron, the nature of the factors involved has remained elusive. We assayed the interaction of U1 snRNP as well as the positive effect of a downstream 5' splice site on trans-splicing in nematode extracts containing either inactive (early in development) or active (later in development) serine/arginine-rich splicing factors (SR proteins). We have determined that U1 snRNP interacts with the 5' splice site in the downstream exon even in the absence of active SR proteins. In addition, we determined that U1 snRNP-directed loading of U2 snRNP onto the branch site as well as efficient trans-splicing in these inactive extracts could be rescued upon the addition of active SR proteins. Identical results were obtained when we examined the interaction of U1 snRNP as well as the requirement for SR proteins in communication across a cis-spliced intron. Weakening of the 3' splice site uncovered distinct differences, however, in the ability of U1 snRNP to promote U2 addition, dependent upon its position relative to the branch site. These results demonstrate that SR proteins are required for communication between U1 and U2 snRNPs whether this interaction is across introns or exons.     

Alternative splicing of RNA triplets is often regulated and accelerates proteome evolution

Thousands of human genes contain introns ending in NAGNAG (N any nucleotide), where both NAGs can function as 3' splice sites, yielding isoforms that differ by inclusion/exclusion of three bases. However, few models exist for how such splicing might be regulated, and some studies have concluded that NAGNAG splicing is purely stochastic and nonfunctional. Here, we used deep RNA-Seq data from 16 human and eight mouse tissues to analyze the regulation and evolution of NAGNAG splicing. Using both biological and technical replicates to estimate false discovery rates, we estimate that at least 25% of alternatively spliced NAGNAGs undergo tissue-specific regulation in mammals, and alternative splicing of strongly tissue-specific NAGNAGs was 10 times as likely to be conserved between species as was splicing of non-tissue-specific events, implying selective maintenance. Preferential use of the distal NAG was associated with distinct sequence features, including a more distal location of the branch point and presence of a pyrimidine immediately before the first NAG, and alteration of these features in a splicing reporter shifted splicing away from the distal site. Strikingly, alignments of orthologous exons revealed a ～15-fold increase in the frequency of three base pair gaps at 3' splice sites relative to nearby exon positions in both mammals and in Drosophila. Alternative splicing of NAGNAGs in human was associated with dramatically increased frequency of exon length changes at orthologous exon boundaries in rodents, and a model involving point mutations that create, destroy, or alter NAGNAGs can explain both the increased frequency and biased codon composition of gained/lost sequence observed at the beginnings of exons. This study shows that NAGNAG alternative splicing generates widespread differences between the proteomes of mammalian tissues, and suggests that the evolutionary trajectories of mammalian proteins are strongly biased by the locations and phases of the introns that interrupt coding sequences.