Exonic splicing enhancers contribute to the use of both 3' and 5' splice site usage of rat beta-tropomyosin pre-mRNA

The rat beta-tropomyosin gene encodes two tissue-specific isoforms that contain the internal, mutually exclusive exons 6 (nonmuscle/smooth muscle) and 7 (skeletal muscle). We previously demonstrated that the 3' splice site of exon 6 can be activated by introducing a 9-nt polyuridine tract at its 3' splice site, or by strengthening the 5' splice site to a U1 consensus binding site, or by joining exon 6 to the downstream common exon 8. Examination of sequences within exons 6 and 8 revealed the presence of two purine-rich motifs in exon 6 and three purine-rich motifs in exon 8 that could potentially represent exonic splicing enhancers (ESEs). In this report we carried out substitution mutagenesis of these elements and show that some of them play a critical role in the splice site usage of exon 6 in vitro and in vivo. Using UV crosslinking, we have identified SF2/ASF as one of the cellular factors that binds to these motifs. Furthermore, we show that substrates that have mutated ESEs are blocked prior to A-complex formation, supporting a role for SF2/ASF binding to the ESEs during the commitment step in splicing. Using pre-mRNA substrates containing exons 5 through 8, we show that the ESEs within exon 6 also play a role in cooperation between the 3' and 5' splice sites flanking this exon. The splicing of exon 6 to 8 (i.e., 5' splice site usage of exon 6) was enhanced with pre-mRNAs containing either the polyuridine tract in the 3' splice site or consensus sequence in the 5' splice site around exon 6. We show that the ESEs in exon 6 are required for this effect. However, the ESEs are not required when both the polyuridine and consensus splice site sequences around exon 6 were present in the same pre-mRNA. These results support and extend the exon-definition hypothesis and demonstrate that sequences at the 3' splice site can facilitate use of a downstream 5' splice site. In addition, the data support the hypothesis that ESEs can compensate for weak splice sites, such as those found in alternatively spliced exons, thereby providing a target for regulation.     

Distribution of exonic splicing enhancer elements in human genes

Ab initio prediction of functional exon splicing enhancer (ESE) elements based on RNA sequences present a challenge in the evaluation of the functional impacts of human genetic polymorphisms on splicing. To better understand the behavior of ESEs, we studied their distribution in human exons and introns for four known SR protein-binding motifs: SF2/SAF, SC35, SRp40, and SRp55. ESEs are enriched in regions in exons that are close to the splice sites, especially in the region 80 to 120 bases away from the ends of splice acceptor sites. Significant enrichment of ESEs is associated with weak splice acceptor sites but not weak donor sites. ESE density decreases at the 3 ends of long exons. ESEs are also enriched in introns with weak donor or acceptor sites. These characteristics of ESEs may help to predict functional ESE sites in RNA sequences.     

Inference of splicing regulatory activities by sequence neighborhood analysis

Sequence-specific recognition of nucleic-acid motifs is critical to many cellular processes. We have developed a new and general method called Neighborhood Inference (NI) that predicts sequences with activity in regulating a biochemical process based on the local density of known sites in sequence space. Applied to the problem of RNA splicing regulation, NI was used to predict hundreds of new exonic splicing enhancer (ESE) and silencer (ESS) hexanucleotides from known human ESEs and ESSs. These predictions were supported by cross-validation analysis, by analysis of published splicing regulatory activity data, by sequence-conservation analysis, and by measurement of the splicing regulatory activity of 24 novel predicted ESEs, ESSs, and neutral sequences using an in vivo splicing reporter assay. These results demonstrate the ability of NI to accurately predict splicing regulatory activity and show that the scope of exonic splicing regulatory elements is substantially larger than previously anticipated. Analysis of orthologous exons in four mammals showed that the NI score of ESEs, a measure of function, is much more highly conserved above background than ESE primary sequence. This observation indicates a high degree of selection for ESE activity in mammalian exons, with surprisingly frequent interchangeability between ESE sequences.     

Regulation of alternative RNA splicing by exon definition and exon sequences in viral and mammalian gene expression

Intron removal from a pre-mRNA by RNA splicing was once thought to be controlled mainly by intron splicing signals. However, viral and other eukaryotic RNA exon sequences have recently been found to regulate RNA splicing, polyadenylation, export, and nonsense-mediated RNA decay in addition to their coding function. Regulation of alternative RNA splicing by exon sequences is largely attributable to the presence of two majorcis-acting elements in the regulated exons, the exonic splicing enhancer (ESE) and the suppressor or silencer (ESS). Two types of ESEs have been verified from more than 50 genes or exons: purine-rich ESEs, which are the more common, and non-purine-rich ESEs. In contrast, the sequences of ESSs identified in approximately 20 genes or exons are highly diverse and show little similarity to each other. Through interactions with cellular splicing factors, an ESE or ESS determines whether or not a regulated splice site, usually an upstream 3 splice site, will be used for RNA splicing. However, how these elements function precisely in selecting a regulated splice site is only partially understood. The balance between positive and negative regulation of splice site selection likely depends on thecis-element's identity and changes in cellular splicing factors under physiological or pathological conditions.     

Exon inclusion is dependent on predictable exonic splicing enhancers

We have previously formulated a list of approximately 2,000 RNA octamers as putative exonic splicing enhancers (PESEs) based on a statistical comparison of human exonic and nonexonic sequences (X. H. Zhang and L. A. Chasin, Genes Dev. 18:1241-1250, 2004). When inserted into a poorly spliced test exon, all eight tested octamers stimulated splicing, a result consistent with their identification as exonic splicing enhancers (ESEs). Here we present a much more stringent test of the validity of this list of PESEs. Twenty-two naturally occurring examples of nonoverlapping PESEs or PESE clusters were identified in six mammalian exons; five of the six exons tested are constitutively spliced. Each of the 22 individual PESEs or PESE clusters was disrupted by site-directed mutagenesis, usually by a single-base substitution. Eighteen of the 22 disruptions (82%) resulted in decreased splicing efficiency. In contrast, 24 control mutations had little or no effect on splicing. This high rate of success suggests that most PESEs function as ESEs in their natural context. Like most exons, these exons contain several PESEs. Since knocking out any one of several could produce a severalfold decrease in splicing efficiency, we conclude that there is little redundancy among ESEs in an exon and that they must work in concert to optimize splicing.     

ESEfinder: A web resource to identify exonic splicing enhancers

Point mutations frequently cause genetic diseases by disrupting the correct pattern of pre-mRNA splicing. The effect of a point mutation within a coding sequence is traditionally attributed to the deduced change in the corresponding amino acid. However, some point mutations can have much more severe effects on the structure of the encoded protein, for example when they inactivate an exonic splicing enhancer (ESE), thereby resulting in exon skipping. ESEs also appear to be especially important in exons that normally undergo alternative splicing. Different classes of ESE consensus motifs have been described, but they are not always easily identified. ESEfinder (http://exon.cshl.edu/ESE/) is a web-based resource that facilitates rapid analysis of exon sequences to identify putative ESEs responsive to the human SR proteins SF2/ASF, SC35, SRp40 and SRp55, and to predict whether exonic mutations disrupt such elements.     

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.     

A general role for splicing enhancers in exon definition

Exonic splicing enhancers (ESEs) facilitate exon definition by assisting in the recruitment of splicing factors to the adjacent intron. Here we demonstrate that suboptimal 5' and 3' splice sites are activated independently by ESEs when they are located on different exons. However, when they are situated within a single exon, the same weak 5' and 3' splice sites are activated simultaneously by a single ESE. These findings demonstrate that a single ESE promotes the recognition of both exon/intron junctions within the same step during exon definition. Our results suggest that ESEs recruit a multicomponent complex that minimally contains components of the splicing machinery required for 5' and 3' splice site selection.     

Evidence for purifying selection against synonymous mutations in mammalian exonic splicing enhancers

Silent sites in mammals have classically been assumed to be free from selective pressures. Consequently, the synonymous substitution rate (Ks) is often used as a proxy for the mutation rate. Although accumulating evidence demonstrates that the assumption is not valid, the mechanism by which selection acts remain unclear. Recent work has revealed that the presence of exonic splicing enhancers (ESEs) in coding sequence might influence synonymous evolution. ESEs are predominantly located near intron-exon junctions, which may explain the reduced single-nucleotide polymorphism (SNP) density in these regions. Here we show that synonymous sites in putative ESEs evolve more slowly than the remaining exonic sequence. Differential mutabilities of ESEs do not appear to explain this difference. We observe that substitution frequency at fourfold synonymous sites decreases as one approaches the ends of exons, consistent with the existing SNP data. This gradient is at least in part explained by ESEs being more abundant near junctions. Between-gene variation in Ks is hence partly explained by the proportion of the gene that acts as an ESE. Given the relative abundance of ESEs and the reduced rates of synonymous divergence within them, we estimate that constraints on synonymous evolution within ESEs causes the true mutation rate to be underestimated by not more than approximately 8%. We also find that Ks outside of ESEs is much lower in alternatively spliced exons than in constitutive exons, implying that other causes of selection on synonymous mutations exist. Additionally, selection on ESEs appears to affect nonsynonymous sites and may explain why amino acid usage near intron-exon junctions is nonrandom.     

Specific binding of an exonic splicing enhancer by the pre-mRNA splicing factor SRp55.

Purine-rich exonic splicing enhancers (ESEs) have been identified in many alternatively spliced exons. Alternative splicing of several ESE-containing exons has been shown to depend on subsets of the SR protein family of pre-mRNA splicing factors. In this report, we show that purified SR protein family member SRp55 by itself binds a 30-nt ESE-containing exon, the alternatively spliced exon 5 of avian cardiac troponin T. We show that purified SRp55 binds specifically to this RNA sequence with an apparent Kd of 60 nM as assayed by gel mobility retardation experiments. Mutations in the exon 5 sequence that increase or decrease exon 5 inclusion in vivo and in vitro have correspondingly different affinities for SRp55 in our assays. The exon 5 sequence contains two purine-rich motifs, common to many ESEs, and both are required for SRp55 binding. Hill plot analysis of binding titration reactions indicates that there is a cooperative binding of at least two SRp55 proteins to the exon sequence. Chemical modification interference studies using kethoxal show that SRp55 binding to exon 5 requires the N1 and/or the N2 of almost every G residue in the exon. Dimethylsulfate modification interference studies indicate that none of the N1 positions of A residues in the exon are important for binding. We postulate that SRp55 may recognize both primary sequence and RNA secondary structural elements within pre-mRNA.     

A strong exonic splicing enhancer in dystrophin exon 19 achieve proper splicing without an upstream polypyrimidine tract

Proper splicing is known to proceed under the control of conserved cis-elements located at exon-intron boundaries. Recently, it was shown that additional elements, such as exonic splicing enhancers (ESEs), are essential for the proper splicing of certain exons, in addition to the splice donor and acceptor site sequences; however, the relationship between these cis-elements is still unclear. In this report, we utilize dystrophin exon 19 to analyse the relationship between the ESE and its upstream acceptor site sequences. Dystrophin exon 19, which maintains adequate splicing donor and acceptor consensus sequences, encodes exonic splicing enhancer (dys-ESE19) sequences. Splice pattern analysis, using a minigene reporter expressed in HeLa cells, showed that either a strong polypyrimidine tract (PPT) or a fully active dys-ESE19 is sufficient for proper splicing. Each of these two cis-elements has enough activity for proper exon 19 splicing suggesting that the PPT, which is believed to be an essential cis-element for splicing, is dispensable when the downstream exon contains a strong ESE. This compensation was only seen in living cells but not in 'in vitro splicing'. This suggests the possibility that the previous splicing experiments using an in vitro splicing system could underestimate the activity of ESEs.     

Molecular design of a splicing switch responsive to the RNA binding protein Tra2β

Tra2β regulates a number of splicing switches including activation of the human testis-specific exon HIPK3-T in the Homeodomain Interacting Protein Kinase 3 gene. By testing HIPK3-T exons of different intrinsic strengths, we found Tra2β most efficiently activated splicing inclusion of intrinsically weak exons, although these were spliced at a lower overall level. Both the RRM and N-terminal RS-rich region of Tra2β were required for splicing activation. Bioinformatic searches for splicing enhancers and repressors mapped four physically distinct exonic splicing enhancers (ESEs) within HIPK3-T, each containing the known Tra2β AGAA-rich binding site. Surprisingly disruption of each single ESE prevented Tra2β-mediated activation, although single mutated exons could still bind Tra2β protein by gel shifts and functional splicing analyses. Titration experiments indicate an additive model of HIPK3-T splicing activation, requiring availability of an array of four distinct ESEs to enable splicing activation. To enable this efficient Tra2β-mediated splicing switch to operate, a closely adjacent downstream and potentially competitive stronger 5'-splice site is actively repressed. Our data indicate that a novel arrangement of multiple mono-specific AGAA-rich ESEs coupled to a weak 5'-splice site functions as a responsive gauge. This gauge monitors changes in the specific nuclear concentration of the RNA binding protein Tra2β, and co-ordinately regulates HIPK3-T exon splicing inclusion.     

Identification and experimental validation of splicing regulatory elements in Drosophila melanogaster reveals functionally conserved splicing enhancers in metazoans

RNA sequence elements involved in the regulation of pre-mRNA splicing have previously been identified in vertebrate genomes by computational methods. Here, we apply such approaches to predict splicing regulatory elements in Drosophila melanogaster and compare them with elements previously found in the human, mouse, and pufferfish genomes. We identified 99 putative exonic splicing enhancers (ESEs) and 231 putative intronic splicing enhancers (ISEs) enriched near weak 5' and 3' splice sites of constitutively spliced introns, distinguishing between those found near short and long introns. We found that a significant proportion (58%) of fly enhancer sequences were previously reported in at least one of the vertebrates. Furthermore, 20% of putative fly ESEs were previously identified as ESEs in human, mouse, and pufferfish; while only two fly ISEs, CTCTCT and TTATAA, were identified as ISEs in all three vertebrate species. Several putative enhancer sequences are similar to characterized binding-site motifs for Drosophila and mammalian splicing regulators. To provide additional evidence for the function of putative ISEs, we separately identified 298 intronic hexamers significantly enriched within sequences phylogenetically conserved among 15 insect species. We found that 73 putative ISEs were among those enriched in conserved regions of the D. melanogaster genome. The functions of nine enhancer sequences were verified in a heterologous splicing reporter, demonstrating that these sequences are sufficient to enhance splicing in vivo. Taken together, these data identify a set of predicted positive-acting splicing regulatory motifs in the Drosophila genome and reveal regulatory sequences that are present in distant metazoan genomes.     

Exonic splicing regulatory elements skew synonymous codon usage near intron-exon boundaries in mammals

In mammals there is a bias in amino acid usage near splice sites that is explained, in large part, by the high density of exonic splicing enhancers (ESEs) in these regions. Is there a similar bias for the relative use of synonymous codons, and can any such bias be predicted by their abundance in ESEs? Prior reports suggested that such trends may exist. From analysis of human exons, we find that 47 of the 59 codons with at least one synonym show differential usage in the proximity of exon ends, of which 42 remain significant after correction for multiple testing. Within sets of synonymous codons those more preferred near splice sites are generally those that are relatively more abundant within the ESEs. However, the examples given previously appear exceptionally good fits and there exist many exceptions, the usage of lysine's codons being a case in point. Similar results are observed in mouse exons. We conclude that splice regulation impacts on the choice of synonymous codons in mammals, but the magnitude of this effect is less than might at first have been supposed.     

Missense mutations in cancer suppressor gene TP53 are colocalized with exonic splicing enhancers (ESEs)

Mutation databases can be viewed as footprints of functional organization of a gene and thus can be used to infer its functional organization. We studied the association of exonic splicing enhancers (ESEs) with missense mutations in the tumor suppressor gene TP53 using the International Agency for Research on Cancer (IARC) mutation database. The goals of the study were: (i) to verify the hypothesis that deleterious missense mutations are colocalized with ESEs; (ii) to identify potentially functional ESE sites in the open reading frame (ORF) of the TP53. If some sequence functions as a splicing enhancer, then nucleotide substitutions in the site will disturb splicing, abrogate p53 function, and cause an increased susceptibility to cancer. Therefore, among cancers showing p53 mutations, more missense mutations are expected within functional ESE sites as compared to non-functional ESE motifs. Using several statistical tests, we found that missense mutations in TP53 are strongly colocalized with ESEs, and that only a small fraction of ESE sites contributes to the association. There are usually one or two ESEs per exon showing a statistically significant association with missense mutations--so-called significant ESE sites. In many respects significant ESE sites are different from those that do not show association with missense mutations. We found that positions of significant ESE sites are codon-dependent--significant ESEs preferentially start from the first position of a codon, whereas non-significant ESEs show no position dependence. Significant ESEs showed a more limited set of sequences compared to non-significant ESEs. These findings suggest that there is a limited number of missense mutations that influence ESE sites and our analysis provides further insight into the types of sites that harbor exonic enhancer elements.     

Pre-mRNA secondary structure and the regulation of splicing

Nuclear pre-mRNAs must be precisely processed to give rise to mature cytoplasmic mRNAs. This maturation process, known as splicing, involves excision of intron sequences and ligation of the exon sequences. One of the major problems in understanding this process is how splice sites, the sequences which form the boundaries between introns and exons, can be accurately selected. A number of studies have defined conserved sequences within introns which were later shown to interact with small nuclear ribonucleoproteins (snRNPs). However, due to the simplicity of these conserved sequences it has become clear that other elements must be involved and a number of studies have indicated the importance of secondary structures within pre-mRNAs. Using various examples, we shall show that such structures can help to specify splice sites by modifying physical distances within introns or by being involved in the definition of exons and, lastly, that they can be part of the regulation of alternative splicing.