Systematic identification and analysis of exonic splicing silencers

Exonic splicing silencers (ESSs) are cis-regulatory elements that inhibit the use of adjacent splice sites, often contributing to alternative splicing (AS). To systematically identify ESSs, an in vivo splicing reporter system was developed to screen a library of random decanucleotides. The screen yielded 141 ESS decamers, 133 of which were unique. The silencer activity of over a dozen of these sequences was also confirmed in a heterologous exon/intron context and in a second cell type. Of the unique ESS decamers, most could be clustered into groups to yield seven putative ESS motifs, some resembling known motifs bound by hnRNPs H and A1. Potential roles of ESSs in constitutive splicing were explored using an algorithm, ExonScan, which simulates splicing based on known or putative splicing-related motifs. ExonScan and related bioinformatic analyses suggest that these ESS motifs play important roles in suppression of pseudoexons, in splice site definition, and in AS.     

Comparative analysis identifies exonic splicing regulatory sequences--The complex definition of enhancers and silencers

Exonic splicing regulatory sequences (ESRs) are cis-acting factor binding sites that regulate constitutive and alternative splicing. A computational method based on the conservation level of wobble positions and the overabundance of sequence motifs between 46,103 human and mouse orthologous exons was developed, identifying 285 putative ESRs. Alternatively spliced exons that are either short in length or contain weak splice sites show the highest conservation level of those ESRs, especially toward the edges of exons. ESRs that are abundant in those subgroups show a different distribution between constitutively and alternatively spliced exons. Representatives of these ESRs and two SR protein binding sites were shown, experimentally, to display variable regulatory effects on alternative splicing, depending on their relative locations in the exon. This finding signifies the delicate positional effect of ESRs on alternative splicing regulation.     

Simultaneous assessment of the effects of exonic mutations on RNA splicing and protein functions

To simultaneously assess the effects of exonic mutations on RNA splicing and protein functions, we report here an intron-inclusive cDNA (Intinc) expression system. As a test model, twenty-four mutations in exon 9 of the phenylalanine hydroxylase (PAH) gene were examined in an Intinc expression plasmid composed of the PAH cDNA with the exon 9 flanked by its authentic introns. When the PAH enzyme activities from the Intinc plasmid-transfected cells were compared to those of a standard cDNA expression system, five mutations resulted in significant relative differences in PAH activities attributed to altered exon 9-inclusive mRNA levels. Two of the mutations affected exon recognition probably through splice site modifications and the remaining three affected experimentally verified exon splicing enhancer (ESE) motifs. The Intinc expression system allows not only a better link between mutation genotype to disease phenotype but also contributes to further understanding of molecular mechanisms of deleterious effects of mutations.     

The splicing regulator Sam68 binds to a novel exonic splicing silencer and functions in SMN2 alternative splicing in spinal muscular atrophy

Spinal muscular atrophy (SMA) is a neurodegenerative disease caused by loss of motor neurons in patients with null mutations in the SMN1 gene. An almost identical SMN2 gene is unable to compensate for this deficiency because a single C‐to‐T transition at position +6 in exon‐7 causes skipping of the exon by a mechanism not yet fully elucidated. We observed that the C‐to‐T transition in SMN2 creates a putative binding site for the RNA‐binding protein Sam68. RNA pull‐down assays and UV‐crosslink experiments showed that Sam68 binds to this sequence. In vivo splicing assays showed that Sam68 triggers SMN2 exon‐7 skipping. Moreover, mutations in the Sam68‐binding site of SMN2 or in the RNA‐binding domain of Sam68 completely abrogated its effect on exon‐7 skipping. Retroviral infection of dominant‐negative mutants of Sam68 that interfere with its RNA‐binding activity, or with its binding to the splicing repressor hnRNP A1, enhanced exon‐7 inclusion in endogenous SMN2 and rescued SMN protein expression in fibroblasts of SMA patients. Our results thus indicate that Sam68 is a novel crucial regulator of SMN2 splicing.     

Activation and repression functions of an SR splicing regulator depend on exonic versus intronic-binding position

SR proteins and related factors play widespread roles in alternative pre-mRNA splicing and are known to promote splice site recognition through their Arg-Ser-rich effector domains. However, binding of SR regulators to some targets results in repression of splice sites through a distinct mechanism. Here, we investigate how activated and repressed targets of the Drosophila SR regulator Transformer2 elicit its differing effects on splicing. We find that, like activation, repression affects early steps in the recognition of splice sites and spliceosome assembly. Repositioning of regulatory elements reveals that Tra2 complexes that normally repress splicing from intronic positions activate splicing when located in an exon. Protein tethering experiments demonstrate that this position dependence is an intrinsic property of Tra2 and further show that repression and activation are mediated by separate effector domains of this protein. When other Drosophila SR factors (SF2 and Rbp1) that activate splicing from exonic positions were tethered intronically they failed to either activate or repress splicing. Interestingly, both activities of Tra2 favor the exonic identity of the RNA sequences that encompass its binding sites. This suggests a model in which these two opposite functions act in concert to define both the position and extent of alternatively spliced exons.     

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.     

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.     

Recognition of exonic splicing enhancer sequences by the Drosophila splicing repressor RSF1.

The Drosophila repressor splicing factor 1 (RSF1) comprises an N-terminal RNA-binding region and a C-terminal domain rich in glycine, arginine and serine residues, termed the GRS domain. Recently, RSF1 has been shown to antagonize splicing factors of the serine/arginine-rich (SR) family and it is, therefore, expected to play a role in processing of a subset of Drosophila pre-mRNAs through specific interactions with RNA. To investigate the RNA-binding specificity of RSF1, we isolated RSF1-binding RNAs using an in vitro selection approach. We have identified two RNA target motifs recognized by RSF1, designated A (CAACGACGA)- and B (AAACGCGCG)-type sequences. We show here that the A-type cognate sequence behaves as an SR protein-dependent exonic splicing enhancer. Namely, three copies of the A-type ligand bind SR proteins, stimulate the efficiency of splicing of reporter pre-mRNAs several fold and lead to inclusion of a short internal exon both in vitro and in vivo. However, three copies of a B-type ligand were much less active. The finding that RSF1 acts as a potent repressor of pre-mRNA splicing in vitro led us to propose that the equilibrium between a limited number of structurally-related general splicing activators or repressors, competing for common or promiscuous binding sites, may be a major determinant of the underlying mechanisms controlling many alternative pre-mRNA process-ing events.     

Characterization of the intronic splicing silencers flanking FGFR2 exon IIIb

The cell type-specific alternative splicing of FGFR2 pre-mRNA results in the mutually exclusive use of exons IIIb and IIIc, which leads to critically important differences in receptor function. The choice of exon IIIc in mesenchymal cells involves activation of this exon and repression of exon IIIb. This repression is mediated by the function of upstream and downstream intronic splicing silencers (UISS and DISS). Here we present a detailed characterization of the determinants of silencing function within UISS and DISS. We used a systematic mutational analysis, introducing deletions and substitutions to define discrete elements within these two silencers of exon IIIb. We show that UISS requires polypyrimidine tract-binding protein (PTB)-binding sites, which define the UISS1 sub-element, and an eight nucleotide sequence 5'-GCAGCACC-3' (UISS2) that is also required. Even though UISS2 does not bind PTB, the full UISS can be replaced with a synthetic silencer designed to provide optimal PTB binding. DISS is composed of a 5'-conserved sub-element (5'-CE) and two regions that contain multiple PTB sites and are functionally redundant (DISS1 and DISS2). DISS1 and DISS2 are separated by the activator sequence IAS2, and together these opposing elements form the intronic control element. Deletion of DISS in the FGFR2 exon IIIb context resulted in the near full inclusion of exon IIIb, and insertion of this silencer downstream of a heterologous exon with a weak 5' splice site was capable of repressing exon inclusion. Extensive deletion analysis demonstrated that the majority of silencing activity could be mapped to the conserved octamer CUCGGUGC within the 5'CE. Replacement of 5'CE and DISS1 with PTB-binding elements failed to restore repression of exon IIIb. We tested the importance of the relative position of the silencers and of the subelements within each silencer. Whereas UISS1, UISS2, DISS1, and DISS2 appear somewhat malleable, the 5'CE is rigid in terms of relative position and redundancy. Our data defined elements of function within the ISSs flanking exon IIIb and suggested that silencing of this exon is mediated by multiple trans-acting factors.     

Functional properties and evolutionary splicing constraints on a composite exonic regulatory element of splicing in CFTR exon 12

In general, splicing regulatory elements are defined as Enhancers or Silencers depending on their positive or negative effect upon exon inclusion. Often, these sequences are usually present separate from each other in exonic/intronic sequences. The Composite Exonic Splicing Regulatory Elements (CERES) represent an extreme physical overlap of enhancer/silencer activity. As a result, when CERES elements are mutated the consequences on the splicing process are difficult to predict. Here, we show that the functional activity of the CERES2 sequence in CFTR exon 12 is regulated by the binding, in very close proximity to each other, of several SR and hnRNP proteins. Moreover, our results show that practically the entire exon 12 sequence context participate in its definition. The consequences of this situation can be observed at the evolutionary level by comparing changes in conservation of different splicing elements in different species. In conclusion, our study highlights how it is increasingly difficult to define many exonic sequences by simply breaking them down in isolated enhancer/silencer or even neutral elements. The real picture is close to one of continuous competition between positive and negative factors where affinity for the target sequences and other dynamic factors decide the inclusion or exclusion of the exon.     

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.     

Distribution of SR protein exonic splicing enhancer motifs in human protein-coding genes

Exonic splicing enhancers (ESEs) are pre-mRNA cis-acting elements required for splice-site recognition. We previously developed a web-based program called ESEfinder that scores any sequence for the presence of ESE motifs recognized by the human SR proteins SF2/ASF, SRp40, SRp55 and SC35 (http://rulai.cshl.edu/tools/ESE/). Using ESEfinder, we have undertaken a large-scale analysis of ESE motif distribution in human protein-coding genes. Significantly higher frequencies of ESE motifs were observed in constitutive internal protein-coding exons, compared with both their flanking intronic regions and with pseudo exons. Statistical analysis of ESE motif frequency distributions revealed a complex relationship between splice-site strength and increased or decreased frequencies of particular SR protein motifs. Comparison of constitutively and alternatively spliced exons demonstrated slightly weaker splice-site scores, as well as significantly fewer ESE motifs, in the alternatively spliced group. Our results underline the importance of ESE-mediated SR protein function in the process of exon definition, in the context of both constitutive splicing and regulated alternative splicing.     

Alternative splicing of CD200 is regulated by an exonic splicing enhancer and SF2/ASF

CD200, a type I membrane glycoprotein, plays an important role in prevention of inflammatory disorders, graft rejection, autoimmune diseases and spontaneous fetal loss. It also regulates tumor immunity. A truncated CD200 (CD200_tr) resulting from alternative splicing has been identified and characterized as a functional antagonist to full-length CD200. Thus, it is important to explore the mechanism(s) controlling alternative splicing of CD200. In this study, we identified an exonic splicing enhancer (ESE) located in exon 2, which is a putative binding site for a splicing regulatory protein SF2/ASF. Deletion or mutation of the ESE site decreased expression of the full-length CD200. Direct binding of SF2/ASF to the ESE site was confirmed by RNA electrophoretic mobility shift assay (EMSA). Knockdown of expression of SF2/ASF resulted in the same splicing pattern as seen after deletion or mutation of the ESE, whereas overexpression of SF2/ASF increased expression of the full-length CD200. In vivo studies showed that viral infection reversed the alternative splicing pattern of CD200 with increased expression of SF2/ASF and the full-length CD200. Taken together, our data suggest for the first time that SF2/ASF regulates the function of CD200 by controlling CD200 alternative splicing, through direct binding to an ESE located in exon 2 of CD200.     

Comparison of intron-containing and intron-lacking human genes elucidates putative exonic splicing enhancers

Of the rules used by the splicing machinery to precisely determine intron-exon boundaries only a fraction is known. Recent evidence suggests that specific short sequences within exons help in defining these boundaries. Such sequences are known as exonic splicing enhancers (ESE). A possible bioinformatical approach to studying ESE sequences is to compare genes that harbor introns with genes that do not. For this purpose two non-redundant samples of 719 intron-containing and 63 intron-lacking human genes were created. We performed a statistical analysis on these datasets of intron-containing and intron-lacking human coding sequences and found a statistically significant difference (P = 0.01) between these samples in terms of 5-6mer oligonucleotide distributions. The difference is not created by a few strong signals present in the majority of exons, but rather by the accumulation of multiple weak signals through small variations in codon frequencies, codon biases and context-dependent codon biases between the samples. A list of putative novel human splicing regulation sequences has been elucidated by our analysis.     

Selection of the bovine papillomavirus type 1 nucleotide 3225 3' splice site is regulated through an exonic splicing enhancer and its juxtaposed exonic splicing suppressor.

Alternative splicing is an important mechanism for the regulation of bovine papillomavirus type 1 (BPV-1) gene expression during the virus life cycle. However, one 3' splice site, located at nucleotide (nt) 3225, is used for the processing of most BPV-1 pre-mRNAs in BPV-1-transformed C127 cells and at early to intermediate times in productively infected warts. At late stages of the viral life cycle, an alternative 3' splice site at nt 3605 is used for the processing of the late pre-mRNA. In this study, we used in vitro splicing in HeLa cell nuclear extracts to identify cis elements which regulate BPV-1 3' splice site selection. Two purine-rich exonic splicing enhancers were identified downstream of nt 3225. These sequences, designated SE1 (nt 3256 to 3305) and SE2 (nt 3477 to 3526), were shown to strongly stimulate the splicing of a chimeric Drosophila doublesex pre-mRNA, which contains a weak 3' splice site. A BPV-1 late pre-mRNA containing the nt 3225 3' splice site but lacking both SE1 and SE2 was spliced poorly, indicating that this 3' splice site is inherently weak. Analysis of the 3' splice site suggested that this feature is due to both a nonconsensus branch point sequence and a suboptimal polypyrimidine tract. Addition of SE1 to the late pre-mRNA dramatically stimulated splicing, indicating that SE1 also functions as an exonic splicing enhancer in its normal context. However, a late pre-mRNA containing both SE1 and SE2 as well as the sequence in between was spliced inefficiently. Further mapping studies demonstrated that a 48-nt pyrimidine-rich region immediately downstream of SE1 was responsible for this suppression of splicing. Thus, these data suggest that selection of the BPV-1 nt 3225 3' splice site is regulated by both positive and negative exonic sequences.     

Multisite and bidirectional exonic splicing enhancer in CD44 alternative exon v3

The human CD44 gene encodes multiple isoforms of a transmembrane protein that differ in their extracellular domains as a result of alternative splicing of its variable exons. Expression of CD44 is tightly regulated according to the type and physiological status of a cell, with expression of high molecular weight isoforms by inclusion of variable exons and low molecular weight isoforms containing few or no variable exons. Human CD44 variable exon 3 (v3) can follow a specific alternative splicing route different from that affecting other variable exons. Here we map and functionally describe the splicing enhancer element within CD44 exon v3 which regulates its inclusion in the final mRNA. The v3 splicing enhancer is a multisite bipartite element consisting of a tandem nonamer, the XX motif, and an heptamer, the Y motif, located centrally in the exon. Each of the three sites of this multisite enhancer partially retains its splicing enhancing capacity independently from each other in CD44 and shows full enhancing function in gene contexts different from CD44. We further demonstrate that these motifs act cooperatively as at least two motifs are needed to maintain exon inclusion. Their action is differential with respect to the splice-site target abutting v3. The first X motif acts on the 3' splice site, the second X motif acts on both splice sites (as a bidirectional exonic splicing enhancer), and the Y motif acts on the 5' splice site. We also show that the multisite v3 splicing enhancer is functional irrespective of flanking intron length and spatial organization within v3.     

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.