Single nucleotide polymorphism-based validation of exonic splicing enhancers

Because deleterious alleles arising from mutation are filtered by natural selection, mutations that create such alleles will be underrepresented in the set of common genetic variation existing in a population at any given time. Here, we describe an approach based on this idea called VERIFY (variant elimination reinforces functionality), which can be used to assess the extent of natural selection acting on an oligonucleotide motif or set of motifs predicted to have biological activity. As an application of this approach, we analyzed a set of 238 hexanucleotides previously predicted to have exonic splicing enhancer (ESE) activity in human exons using the relative enhancer and silencer classification by unanimous enrichment (RESCUE)-ESE method. Aligning the single nucleotide polymorphisms (SNPs) from the public human SNP database to the chimpanzee genome allowed inference of the direction of the mutations that created present-day SNPs. Analyzing the set of SNPs that overlap RESCUE-ESE hexamers, we conclude that nearly one-fifth of the mutations that disrupt predicted ESEs have been eliminated by natural selection (odds ratio = 0.82 +/- 0.05). This selection is strongest for the predicted ESEs that are located near splice sites. Our results demonstrate a novel approach for quantifying the extent of natural selection acting on candidate functional motifs and also suggest certain features of mutations/SNPs, such as proximity to the splice site and disruption or alteration of predicted ESEs, that should be useful in identifying variants that might cause a biological phenotype.     

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.     

The paradox of viable <Emphasis Type="Italic">sup45</Emphasis> STOP mutations: a necessary equilibrium between translational readthrough,activity and stability of the protein

The mechanisms leading to non-lethality of nonsense mutations in essential genes are poorly understood. Here, we focus on the factors influencing viability of yeast cells bearing premature termination codons (PTCs) in the essential gene SUP45 encoding translation termination factor eRF1. Using a dual reporter system we compared readthrough efficiency of the natural termination codon of SUP45 gene, spontaneous sup45-n (nonsense) mutations, nonsense mutations obtained by site-directed mutagenesis (76Q → TAA, 242R → TGA, 317L → TAG). The nonsense mutations in SUP45 gene were shown to be situated in moderate contexts for readthrough efficiency. We showed that readthrough efficiency of some of the mutations present in the sup45 mutants is not correlated with full-length Sup45 protein amount. This resulted from modification of both sup45 mRNA stability which varies 3-fold among sup45-n mutants and degradation rate of mutant Sup45 proteins. Our results demonstrate that some substitutions in the place of PTCs decrease Sup45 stability. The viability of sup45 nonsense mutants is therefore supported by diverse mechanisms that control the final amount of functional Sup45 in cells.     

Missense mutations in hMLH1 and hMSH2 are associated with exonic splicing enhancers

There is a critical need to understand why missense mutations are deleterious. The deleterious effects of missense mutations are commonly attributed to their impact on primary amino acid sequence and protein structure. However, several recent studies have shown that some missense mutations are deleterious because they disturb cis-acting splicing elements-so-called "exonic splicing enhancers" (ESEs). It is not clear whether the ESE-related deleterious effects of missense mutations are common. We have evaluated colocalization of pathogenic missense mutations (found in affected individuals) with high-score ESE motifs in the human mismatch-repair genes hMSH2 and hMLH1. We found that pathogenic missense mutations in the hMSH2 and hMLH1 genes are located in ESE sites significantly more frequently than expected. Pathogenic missense mutations also tended to decrease ESE scores, thus leading to a higher propensity for splicing defects. In contrast, nonpathogenic missense mutations (polymorphisms found in unaffected individuals) and nonsense mutations are distributed randomly in relation to ESE sites. Comparison of the observed and expected frequencies of missense mutations in ESE sites shows that pathogenic effects of >/=20% of mutations in hMSH2 result from disruption of ESE sites and disturbed splicing. Similarly, pathogenic effects of >/=16% of missense mutations in the hMLH1 gene are ESE related. The colocalization of pathogenic missense mutations with ESE sites strongly suggests that their pathogenic effects are splicing related.     

General and specific functions of exonic splicing silencers in splicing control

Correct splice site recognition is critical in pre-mRNA splicing. We find that almost all of a diverse panel of exonic splicing silencer (ESS) elements alter splice site choice when placed between competing sites, consistently inhibiting use of intron-proximal 5' and 3' splice sites. Supporting a general role for ESSs in splice site definition, we found that ESSs are both abundant and highly conserved between alternative splice site pairs and that mutation of ESSs located between natural alternative splice site pairs consistently shifted splicing toward the intron-proximal site. Some exonic splicing enhancers (ESEs) promoted use of intron-proximal 5' splice sites, and tethering of hnRNP A1 and SF2/ASF proteins between competing splice sites mimicked the effects of ESS and ESE elements, respectively. Further, we observed that specific subsets of ESSs had distinct effects on a multifunctional intron retention reporter and that one of these subsets is likely preferred for regulation of endogenous intron retention events. Together, our findings provide a comprehensive picture of the functions of ESSs in the control of diverse types of splicing decisions.     

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.     

Frame-disrupting mutations elicit pre-mRNA accumulation independently of frame disruption

The T-cell receptor (TCR) and immunoglobulin (Ig) genes are unique among vertebrate genes in that they undergo programmed rearrangement, a process that allows them to generate an enormous array of receptors with different antigen specificities. While crucial for immune function, this rearrangement mechanism is highly error prone, often generating frameshift or nonsense mutations that render the rearranged TCR and Ig genes defective. Such frame-disrupting mutations have been reported to increase the level of TCRβ and Igµ pre-mRNA, suggesting the hypothesis that RNA processing is blocked when frame disruption is sensed. Using a chimeric gene that contains TCRβ sequences conferring this upregulatory response, we provide evidence that pre-mRNA upregulation is neither frame- nor translation-dependent; instead, several lines of evidence suggested that it is the result of disrupted cis elements necessary for efficient RNA splicing. In particular, we identify the rearranging VDJ_β exon as being uniquely densely packed with exonic-splicing enhancers (ESEs), rendering this exon hypersensitive to mutational disruption. As the chimeric gene that we developed for these studies generates unusually stable nuclear pre-mRNAs that accumulate when challenged with ESE mutations, we suggest it can be used as a sensitive in vivo system to identify and characterize ESEs.     

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.     

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.     

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.     

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.     

Factors associated with a purine-rich exonic splicing enhancer sequence in Xenopus oocyte nucleus

Purine-rich exonic splicing enhancers (ESEs) stimulate splicing of the adjacent introns with suboptimal splice sites. To elucidate the mechanism regarding ESEs, factors specifically associated with ESEs in HeLa cell nuclear extracts were previously investigated, and shown to include SR (serine/arginine-rich) proteins. However, factors associated with ESEs in vivo have not yet been explored. Here we show that a GAA repeat RNA sequence, a typical ESE, is associated in Xenopus oocyte nuclei with at least one SR protein, SF2/ASF, as was expected. Moreover, components of SF3a/b complexes, U2 snRNA, and U2AF(65) were also found to be associated with the ESE in the nucleus. Since SF3a/b complexes are the constituents of the 17S U2 snRNP, these results suggest that the 17S U2 snRNP is associated with the ESE in the nucleus, probably through bridging interactions of U2AF and SR proteins. The identified factors may represent a functional splicing enhancer complex in vivo.     

In vitro and in silico analysis reveals an efficient algorithm to predict the splicing consequences of mutations at the 5' splice sites

We have found that two previously reported exonic mutations in the PINK1 and PARK7 genes affect pre-mRNA splicing. To develop an algorithm to predict underestimated splicing consequences of exonic mutations at the 5′ splice site, we constructed and analyzed 31 minigenes carrying exonic splicing mutations and their derivatives. We also examined 189 249 U2-dependent 5′ splice sites of the entire human genome and found that a new variable, the SD-Score, which represents a common logarithm of the frequency of a specific 5′ splice site, efficiently predicts the splicing consequences of these minigenes. We also employed the information contents (R_i) to improve the prediction accuracy. We validated our algorithm by analyzing 32 additional minigenes as well as 179 previously reported splicing mutations. The SD-Score algorithm predicted aberrant splicings in 198 of 204 sites (sensitivity = 97.1%) and normal splicings in 36 of 38 sites (specificity = 94.7%). Simulation of all possible exonic mutations at positions −3, −2 and −1 of the 189 249 sites predicts that 37.8, 88.8 and 96.8% of these mutations would affect pre-mRNA splicing, respectively. We propose that the SD-Score algorithm is a practical tool to predict splicing consequences of mutations affecting the 5′ splice site.     

RNA splicing promotes translation and RNA surveillance

Aberrant mRNAs harboring premature termination codons (PTCs or nonsense codons) are degraded by the nonsense-mediated mRNA decay (NMD) pathway. mRNAs transcribed from genes that naturally acquire PTCs during lymphocyte development are strongly downregulated by PTCs. Here we show that a signal essential for this robust mRNA downregulatory response is efficient RNA splicing. Strong mRNA downregulation can be conferred on a poor NMD substrate by either strengthening its splicing signals or removing its weak introns. Efficient splicing also strongly promotes translation, providing a molecular explanation for enhanced NMD and suggesting that efficient splicing may have evolved to enhance both protein production and RNA surveillance. Our results suggest simple approaches for increasing protein expression from expression vectors and treating human genetic diseases caused by nonsense and frameshift mutations.