Intronic motif pairs cooperate across exons to promote pre-mRNA splicing
© Ke and Chasin; licensee BioMed Central Ltd. 2010
Received: 9 May 2010
Accepted: 12 August 2010
Published: 12 August 2010
A very early step in splice site recognition is exon definition, a process that is as yet poorly understood. Communication between the two ends of an exon is thought to be required for this step. We report genome-wide evidence for exons being defined through the combinatorial activity of motifs located in flanking intronic regions.
Strongly co-occurring motifs were found to specifically reside in four intronic regions surrounding a large number of human exons. These paired motifs occur around constitutive and alternative exons but not pseudo exons. Most co-occurring motifs are limited to intronic regions within 100 nucleotides of the exon. They are preferentially associated with weaker exons. Their pairing is conserved in evolution and they exhibit a lower frequency of single nucleotide polymorphism when paired. Paired motifs display specificity with respect to distance from the exon borders and in constitutive versus alternative splicing. Many resemble binding sites for heterogeneous nuclear ribonucleoproteins. Specific pairs are associated with tissue-specific genes, the higher expression of which coincides with that of the pertinent RNA binding proteins. Tested pairs acted synergistically to enhance exon inclusion, and this enhancement was found to be exon-specific.
The exon-flanking sequence pairs identified here by genomic analysis promote exon inclusion and may play a role in the exon definition step in pre-mRNA splicing. We propose a model in which multiple concerted interactions are required between exonic sequences and flanking intronic sequences to effect exon definition.
All pre-mRNA splicing reactions involve the removal of an intron from between two exons and so require the pairing of the splice sites at the two ends of the intron; such pairing can be considered as a mandatory 'intron definition' step in splicing. However, it is likely that the initial recognition of most splice sites also involves 'exon definition,' the identification of two splice sites across an exon. This idea was first put forth to explain the observation that appending a 5' splice site downstream of the second exon in a two-exon pre-mRNA molecule greatly enhances splicing of the upstream intron in vitro . There has since been a wealth of genetic evidence supporting this idea: the common consequence of mutating one splice site in an internal exon is the skipping of the entire exon, leaving the wild-type splice site at the other end of the exon unused . One can imagine exon definition as serving a quality control function, preventing splicing from occurring at an isolated splice site unless it results in the inclusion of a bona fide exon. Despite the wide acceptance of this idea, especially in metazoans where intron size is much greater than exon size, most biochemical investigations of splicing have focused on protein-protein interaction across introns, rather than on complexes that form across exons [3, 4].
It is possible that spliceosomal components themselves mediate this concurrent recognition of splice sites [5, 6]. For instance, a mutation in a 5' splice site that eliminates splicing can be suppressed by a mutation in the upstream 3' splice site that improves its agreement to the consensus . However, given the surfeit of splicing regulatory motifs , it seems likely that exonic and/or intronic enhancers play a role in exon definition as well. Evolutionary changes that weaken a splice site can be compensated by changes in exonic splicing enhancer (ESE) or silencer (ESS) content and vice versa [9, 10], implying that the exon in its entirety represents an evolutionary unit. Downstream intronic splicing enhancers (ISEs) show specificity for different classes of 5' splice site sequences  and could be contributing to exon definition. Specific and widespread combinations of motifs can also act negatively to promote exon skipping .
A simple first step in the end-to-end recognition of an exon could be the binding of proteins at the two ends of the exon that are capable of specifically interacting with each other. If there is a limited repertoire of such proteins, then their existence should be signaled by the occurrence of specific combinations of sequences that serve as binding sites for these putative exon definition factors. Such pair-wise combinations can act to promote the intron definition step in splicing. The binding of the same heterogeneous nuclear ribonucleoprotein (hnRNP) at the two ends of a long intron can promote splicing [13, 14]. A computational search revealed motifs that co-occur at intron ends and such motif pairs were shown to promote intron removal .
Here we have sought evidence for cis-acting elements that act in combination at an earlier step in splicing, interacting from the two ends of an exon to mediate exon definition. Whereas most past computational searches for cis-acting splicing elements have focused on single motifs [16–18], here we have sought pairs of motifs that demonstrate an unusually strong tendency to co-occur across exons. We have limited ourselves to intronic motifs that are paired across exons for two reasons. First, there is increasing evidence that the intronic flanks of exons can play an important role in splice site recognition [15, 19–21]. Second, a search for motif combination within protein-coding exons is complicated by the possibility of correlation due to the non-random association of protein motifs [22, 23].
We found that more than 15% of exons harbor flanking motif pairs that are strongly associated with each other. These pairs are found around constitutive and alternative exons but not pseudo exons and their pairing is evolutionarily conserved. They are also associated most frequently with exons that appear relatively weak by other criteria. Specific pairs are also associated with tissue-specific genes. When tested in a heterologous context, these motif pairs were found to synergistically enhance exon inclusion. This enhancement proved to be context dependent, with specificity that was imparted by exonic sequences. Thus, the communication between exon ends may involve multiple interactions across the exon and its intronic flanks.
Results and discussion
Co-occurring motifs are found in the intronic flanks of exons
There are about one million pentamer combinations to consider when comparing two regions (45 × 45 = 410). If we set the P-value cutoff at 1/410 (referred to hereafter as 10-6), we expect to see around one pentamer pair having a P-value smaller than this cutoff if pentamers in one intronic region are independent of those in the other region. Examining about 80,000 human constitutive exons, we found more than 60,000 pentamer pairs (approximately 6% of 410) that passed this P-value cutoff. The top motif pairs detected all shared similar GC contents, being either GC-rich or AT-rich (shown for the UpDp region pairs in Figure S1a in Additional file 1 and Table S1 in Additional file 1). A GC content correlation between intronic regions flanking exons was expected due to the widespread occurrence of GC isochores in the human genome  and the exaggeration of this dichotomy in and around exons [20, 26]. This GC content correlation is illustrated for the UpDp intron region pairing as an example in Figure 1b; the correlation coefficient (r) is 0.73. This strong correlation of the two intronic flanking regions is not observed for GA or GT content (r = -0.01 and 0.04, respectively; Figure S1b,c in Additional file 1).
To confirm our suspicion that these pairings were not specific, we performed a control experiment: the 50-nucleotide Up intronic region upstream of each exon was randomly exchanged with that of another exon having the same regional GC content; this procedure was then repeated for the downstream Dp region. This shuffling should greatly decrease the correlation if there were specific intronic pentamer pairs in the original pairings. No such decrease occurred and once again almost all of the pairs passing the 10-6 P-value cutoff were either GC-rich or AT-rich. The P-value distribution of this shuffled control was quite close to that of the original constitutive exons, and both were substantially different from the null hypothesis model (Figure 1C).
As an additional control, we asked whether the co-occurrence of pentamer pairs could also be found around other genomic sequences of a similar size. We examined pseudo exons , defined as deep intronic sequences of typical internal exon size (50 to 250 nucleotides) bounded by sequences resembling 3' and 5' splice sites, but which are never spliced. We applied the BBC-COOC algorithm to a large set (approximately 100,000) of nonredundant pseudo exons, using the same combinations of Ud, Up, Dp, Dd regions as for real exons. All four searches for correlations produced only numbers close to that expected for the null hypothesis (Figure 2, row 3). As a further control we examined the flanks of pseudo splice sites located upstream or downstream of real constitutive exons. That is, we searched the upstream intronic region of constitutive exons and found sequences with better 3' splice site scores than those of the exon and confirmed that these pseudo 3' splice sites were not used for splicing based on EST databases. Pseudo 5' splice sites were defined in the same way. We re-defined Ud, Up, Dp and Dd for these extended constitutive exons and checked the motif correlations of the four regional combinations with BBC-COOC. All four cases generated only background numbers of co-occurring motif pairs (Figure 2, rows 4 and 5). These results support the idea that the co-occurring motif pairs discovered in constitutive exon intronic flanks are involved in splicing and are not general features of the nonrandomness of the human genome. The discovery of particular significantly correlated intronic motif pairs located close to splice sites suggests that they may be working cooperatively across exons to promote exon definition and exon splicing. It may also be worth noting that the absence of co-occurring pairs around pseudo exons argues against such combinations being used to silence these false splice sites.
We next analyzed alternatively spliced exons using BBC-COOC and again found significantly co-occurring motifs. For three of the four regional classes, alternative exons gave rise to only about 40% of the number of motif pairs yielded by constitutive exons. This result might be attributable to the lower statistical power afforded by the smaller number of the former (approximately 35,000) compared to the latter (approximately 80,000). Interestingly, in the regional class UpDd, alternative exons yielded more co-occurring pairs than constitutive exons. This excess of alternative splicing motifs associated with a downstream distal region (more than +50 nucleotides) echoes the discovery of intronic elements regulating the alternative splicing of individual exons (for example, in the control of N-src splicing ) as well as with the global mapping of predicted Nova binding sites . For most of the characterization of co-occurring motif pairs described below, we used the constitutive set to focus on exons with equally strong splicing.
Motif pairs occur close to splice sites
Co-occurring motif pairs exhibit regional specificity
We asked whether the lower number of motifs pairs passing the cutoff of P ≤ 10-6 in the other three regional combinations was due to a lower number of motifs and a consequent loss of statistical power. Such was not the case, as the expected number of motifs pairs (based on the number of individual motifs) was comparable in almost all cases; for 98% of the co-occurring pairs, the lowest number of expected pairs (based on the null hypothesis) was within a factor of two of that for the defining region (UpDp in this case). The same was true for the other three regional combinations shown in Figure 4a.
Motif pairs around alternative and constitutive exons differ
In the same way, we asked to what extent motif pairs discovered around constitutive exons overlapped with those found around alternative exons. Here again we saw specificity: most of the pairs from constitutive exons that passed the 10-6 cutoff were not among those that passed the cutoff from alternative exons and vice versa (Figure 4b). Because the cutoff is quite stringent, this result does not necessarily mean that the constitutive motif pairs are not found around alternative exons. But it could be interpreted to mean that alternative exons make greater use of special motif pairs. An interesting possibility is that the motif pairs found around alternative exons are actually acting negatively to promote alternative exon skipping. We explore this idea further below. Alternatively, the distinction may be secondary to tissue specificity, which is likely to be higher among alternatively spliced exons. The idea that the genes that harbor these constitutive exons are confined to just a few functional classes was ruled out by the observation that they comprise a very wide variety of Gene Ontology classes (data not shown).
Motif pairs are conserved in evolution
Co-occurring pairs have a lower SNP density
The co-occurring pair hypothesis predicts that mutations that occur in these motifs should have a higher likelihood of disrupting exon splicing than those that occur in the same motifs when they are alone. Therefore, the former would be more likely to be eliminated by purifying selection. Thus, the motifs of co-occurring pairs should have a lower SNP density. Consistent with this prediction, for all four regional classes the SNP density was significantly lower when the motifs were paired than when they were unpaired for both human constitutive exons and alternative exons (Figure 6c,d). This observation suggests that motifs of co-occurring pairs have been subject to purifying selection as pairs in recent human evolution and reinforces the conclusion from the human-macaque comparison. SNPs that disrupt a co-occurring pair could result in decreased exon inclusion, a lower level of the protein product and a mutant phenotype. In this way they may provide a class of functional markers for the identification of quantitative traits affecting human phenotypes, including disease associations.
Motif pairs are associated with weaker exons
If intronic co-occurring pairs act to promote splicing, then they might be expected to contribute more frequently to exons that are otherwise relatively deficient in splicing signals. We compared all constitutive exons that contain co-occurring motif pairs of a particular class (that is, UpDp, UpDd, and so on) to the constitutive exons of a set that did not contain such pairs. The exons of the second set were exactly matched to the first set in the GC content of the relevant paired intronic regions so as to minimize the influence of base composition on any correlations seen. For instance, regions high in GC content will tend to be associated with splice sites that are high in GC content , which in turn are associated with poorer splice site consensus scores.
Motif pairs associated with alternatively spliced exons might not have shown a correlation with weak exons if many motif pairs were acting to help silence rather than enhance splicing. However, the statistically significant results in the case of alternatively spliced exons also showed an association with weaker exons (Figure 7b), consistent with motif pairs enhancement of splicing for alternative as well as constitutive exons.
Sequence characteristics of motif pairs
Most of the co-occurring pentamers are GC-rich (Figure 5; Table S2 in Additional file 2) and approximately 90% contain at least one CpG dinucleotide. This high CpG content is notable in light of the low general abundance of CpG in introns due to the mutational vulnerability of the oft-methylated C. Somewhat less than half of exons with co-occurring pairs harbor these CpG-containing motifs (41%). We considered the possibility that the high incidence of CpG dinucleotides in co-occurring pairs might be an artifact caused by internal exons that are located close to the 5' ends of transcripts. The transcription of most human genes is driven by CpG islands that lie upstream of the transcription start site, but that often extend several kilobases beyond it. If so, then pseudo exons should be subject to the same bias, as many of them would also be located near the 5' ends of genes, especially since first introns tend to be long  and would therefore be major contributors to the pseudo exon pool. The absence of co-occurring pairs from around pseudo exons (Figure 2) argues strongly against the possibility that these co-occurring pairs arise from CpG island transcription signals rather than from splicing signals. It should be noted that CpG-rich motifs are characteristic of the binding site of RBM4, a multifunctional RNA binding protein .
Despite the high GC content of most of these pentamers and their attendant sequence simplification, we saw no evidence for complementarity among them; perfectly complementary pairs appear at a frequency (7/211) no greater than that seen among random pentamers with the same overall base composition (for example, 10/211). Thus, secondary structure does not seem to be playing a role in the selection of these motif pairs.
Comparison with previously generated intronic motifs
If the intronic motifs discovered here function to promote splicing, they may overlap with previously reported motifs computationally predicted to do the same. We compared the 38 unique downstream motifs from the constitutive and alternative UpDp classes with pentamers located in downstream intronic flanks that were predicted to be ISEs based on their relative abundance and/or evolutionary conservation [19–21, 24]. There was little overlap among the ISEs (Table S3 in Additional file 1), perhaps because the co-occurring motifs are distinctive in their pairing rather than their individual relative abundances or conservation.
Genomic distribution of motif pairs
The co-occurring motif pairs are abundant: overall, 17% of internal constitutive exons have co-occurring motif pairs in their intronic flank regions. The proximal UpDp combination yielded the greatest number of co-occurring pairs, but all combinations were substantially represented: UpDp, 7.6%; UpDd, 5.0%; UdDp, 3.5%; UdDd, 4.6% (these numbers add up to more than 17% because many exons have more than one class of pairs). Because we set a stringent P-value threshold for detecting these co-occurring pairs, the actual proportion of human constitutive exons with functioning co-occurring pairs may be much higher. This abundance would allow co-occurring motif pairs to play a role in the splicing of many human constitutive exons. For constitutive exons, motif pairs that originate from proximal regions tend to be clustered at the proximal end of the 50-nucleotide region (closer to the splice site); on the contrary, motifs from distal regions are spread throughout the distal region (Figure 4c). Interestingly, the clustering close to the 3' splice site is not seen among alternative exon motifs. Although the Up region spans the usual position of branch points, none of the Up motifs resembles that consensus (Figure 5; Table S2 in Additional file 2).
Many co-occurring motifs resemble hnRNP binding sites
Many of the motifs in the co-occurring pairs resemble the binding sites of hnRNPs or other RNA binding proteins, including hnRNPs A1/A2, C, D, F/H,G, I (PTB), K, L, M, and 9G8 (Table S2 in Additional file 2); more than 30% of the individual motifs fall into this category. Almost all of these RNA binding site motifs are more characteristic of introns than of exons. While hnRNPs have been most often associated with splicing silencing, many of those examples involve binding within exons, and there are many other examples in which hnRNPs play a positive role in splicing from positions outside the exon . The position of such binding sites relative to the exon can play a determining role in their mode of action, as exemplified by Nova sites, which are generally inhibitory downstream of exons but stimulatory upstream . Computationally defined  or experimentally selected  exonic silencer sequences are enriched in the intronic flanks surrounding splice sites, where they may aid in accurate splicing by silencing nearby pseudo sites [16, 36]. Chabot and colleagues have shown that two hnRNP A1 molecules can promote intron definition by binding to the two ends of an intron, with the idea that the interacting proteins bring those ends together . It is tempting to speculate that hnRNPs may function in an analogous manner to bring the two ends of an exon together to effect exon definition.
Genes containing exons with motif pairs show tissue-specific bias
Transcripts subject to tissue-specific alternative splicing often contain intronic regulatory sequences corresponding to the binding sites of splicing factors that are preferentially expressed in the corresponding tissues, such as Nova in brain  and Fox in muscle and brain [38, 39]. Tissue-specific differences in the levels of ubiquitously expressed factors such as serine/arginine-rich proteins (SR proteins) and hnRNPs [40, 41] (see Figure S3 in Additional file 1 for examples) can also play a role in tissue-specific alternative splicing [3, 42]. Because many motifs of co-occurring pairs resemble binding sites of hnRNP proteins, we asked whether genes containing exons with those motif pairs would demonstrate tissue-specific biases.
Motif pairs act synergistically to promote splicing
The abundance of co-occurring motif pairs flanking constitutive exons and their absence around pseudo exons, their exclusive residence within 200 nucleotides of the exon borders, their conservation as pairs in evolution, their association with weaker exons and the resemblance of many of them to known RNA binding protein binding sites suggested that the two motifs of a pair may cooperate across an exon to promote exon splicing. To test this idea we used a three-exon minigene with terminal exons and introns derived from the Chinese hamster dhfr gene. As a central exon we used exon 2 of the human beta globin gene (Hb2), including the 3' and 5' splice sites, that is, the 223-nucleotide exon plus 14 nucleotides upstream and 6 nucleotides downstream. Deprived of its natural intronic flanks, the Hb2 exon undergoes very little inclusion (approximately 2%) in this context. Tandem pentamers from five co-occurring pairs (numbered 1 to 5) were then inserted in various combinations on either side of the Hb2 exon. Each of the first four pairs (TCCCT/GGAGG, CCCCT/AGGGA, TTTCT/TGGTG, TTTCT/GGTGG) were chosen because they are found in over 1,000 human constitutive exons. The fifth pair (CGCCG/CGCGC) has the most significant P-value of association (2.9 × 10-12).
To explore the generality of this result, we repeated the experiment using a different central exon, a weakened version of exon 5 of the Wilm's tumor gene (WT1-5: the 51-nucleotide exon body plus 26 nucleotides upstream and 6 nucleotides downstream). Remarkably, the results of co-occurring pair inclusion were now the reciprocal of what was seen for Hb2: co-occurring pairs 1 to 4 were without effect but pair 5 yielded a robust and synergistic 20% exon inclusion (Figure 9b). A third non-co-occurring pair (N3) was again without effect, as were the pairings of the effective 5U with the ineffective (here) 1 D to 4 D motifs. We next asked whether the distinct behavior of the Hb2 and WT1-5 exons was dictated by the exon body sequences or by the different sequences of the splice sites of these two exons. The exon body sequences (from +2 to -4) were swapped between the two exons without changing the 3' and 5' splice sites; these exons are denoted WHW and HWH, where the central letter indicates the exon body sequence and the outside letters the splice sites (Figure 9c). The results are shown in Figure 9d in terms of proportion skipping rather than inclusion. The simple provision of WT-15 splice sites to the Hb2 exon (WHW) greatly promoted splicing, but about 30% exon skipping still took place. Adding co-occurring pair 2 motifs to this exon reduced skipping to only 4%, with the individual motifs being much less effective (Figure 9d). Once again, provision of co-occurring pair 5 motifs was without effect on this chimeric exon. Thus, the effectiveness of the ISE pair was associated with the Hb2 exon body and not its splice sites. Exon HWH failed to splice even when the co-occurring pair 5 was added; apparently the Hb2 3' splice site (splice site sequence consensus value = 72.7 compared to 85.7 for the WT1-5 3' splice site) is too weak to function in this context. The 5' splice site sequence consensus values are strong and comparable at 86.1 for Hb2 and 83.4 for WT1-5.
Since alternative exons by definition often fail to be included, we explored the possibility that co-occurring pairs around alternative exons cooperate to prevent splicing, rather than enhance it. The WHW exon exhibits 70% inclusion in the absence of any added sequences and so provides a sensitive assay for silencing. Three such co-occurring pairs (A1, TGGGG/CTGGG; A2, CAGTG/CTTCT; A3, GGGCG/GCGCG), chosen because they resemble hnRNP binding sites, were inserted in combinations on either side of the WHW exon. All three pairs demonstrated synergy in enhancing exon inclusion, lending no support to the silencing idea (Figure 9d). Interestingly, A1U and A3U silenced exon inclusion as individual motifs, but splicing was re-established when their partners were included (Figure 9d).
These results lead to two conclusions. The first is that the two motifs of a pair are acting not additively but synergistically, suggesting that they act in a concerted manner. The second is that these motifs act in a context-dependent manner, and the specificity of this context can be provided by exon body sequences.
Exon definition can be thought about at several levels. Perhaps the most fundamental level is as an observable phenomenon: the competence of a splice site at one end of an exon is dependent on competence at the other end [1, 2, 7, 30, 46]. At another level it can be thought of as a concept, as a way to impose quality control on splicing; potential exons can be tried out, but if both ends do not pass the test, then the exon is rejected and another is tried. Exon definition can also be considered as a functional complex with biochemically identified components [1, 3, 4]. In the work described here we found evidence for the coordinated use of ISEs that could play a role in exon definition. A variety of computational evidence supported the conclusion that these elements are functional in splicing and experimental tests validated the idea that pairs of motifs were acting synergistically from the two ends of the exon to promote splicing. These experimental validations also revealed considerable specificity not only in motif pairing but also in interaction with exon body sequences. While it is easy to imagine that interactions between elements spanning an exon are effecting exon definition, exactly what is meant by this process mechanistically remains elusive. Moreover, it may be that these motif pairs are acting at a later step in splicing. For instance, a motif at the 5' end of the downstream intron in a lariat structure may direct interaction with a partner motif just upstream of the exon to hold the cut ends of the RNA together to allow the second step in splicing.
Materials and methods
Human mRNA sequences and ESTs were downloaded from the UniGene database  and were aligned to the assembled genomic sequences (hg18) obtained from the human genome sequence  using Sim4. Only ESTs that spanned at least two exon-exon joints were considered. Genes that exhibited no intron-exon joints were excluded. Exons with no evidence of skipping or alternative splice site usage were identified as constitutive exons. An exon that was excluded in one or more transcripts and present in at least one transcript was defined as an alternative cassette exon; an exon that demonstrated either alternative 3' or 5' splice site usage in at least one transcript was designated as an alternative 3' or 5' splice site exon. Only exons flanked by canonical AG and GT dinucleotides were included.
For all exons that we used for BBC-COOC analysis, we first removed any repeats defined by RepeatMasker . If two alternative 3' splice sites were present within 50 nts we kept only one (chosen randomly). The same procedure was applied to alternative 5' splice site exons. We further removed highly similar subsets of exons by running BLASTN on the related intronic regions of each exon against those of all other exons in the set. For example, in the UpDp intronic region analysis, if two exons showed similarity in both Up and Dp regions, both exons were purged from the qualified exon sets. Finally, we removed exons having related intronic regions with extreme GC contents of less than 20% or greater than 80%. After application of these three filters we were left with approximately 80,000 constitutive exons and 35,000 alternative exons (alternative cassette, alternative 3' splice site exons and alternative 5' splice site exons) for UpDp, UpDd, UdDp and UdDd analysis.
Pseudo exons were defined as having lengths between 50 and 250 nucleotides and consensus values of ≥75 for 3' splice sites and ≥78 for 5' splice sites. The consensus values were based on a position-specific weight matrix and were calculated essentially according to Shapiro and Senapathy . In addition, pseudo exons had to be at least 100 nucleotides away from the closest real exon. We defined pseudo 3'/authentic 5' splice site exons by searching for a pseudo 3' splice site located at least 100 nucleotides upstream of the authentic 3' splice site of constitutive exons to ensure no overlap with the Ud intronic region. We required the consensus value of a pseudo 3' splice site to be at least as high as that of the authentic 3' splice site of the exon. Authentic 3'/pseudo 5' splice site exons were defined analogously. After removing repeats by RepeatMasker, filtering similar sequences and those with extreme GC contents as above, approximately 100,000 pseudo exons, 18,000 pseudo 3'/authentic 5' splice site exons and 13,000 authentic 3'/pseudo 5' splice site exons were obtained for BBC-COOC analysis.
and where the expected value of k is .
From this distribution we can calculate the P-values for all observed values of ka, where ka values greater than the mean indicate a positive correlation and vice versa. This method of dealing with the biased GC contents is similar to one used by the Burge laboratory to discover co-occurring pairs across introns  except for our use of the hypergeometric distribution, rather than a Poisson approximation.
Co-conservation of co-occurring motif pairing
Based on the coordinates of the human exons, orthologous macaque exons were extracted from a 17-genome multi-alignment . In this way were able to survey a total of 59,221 constitutive exon pairs. We then identified human exons that had co-occurring pairs in their related intronic regions. We considered a human intronic motif to be conserved if it also existed anywhere within the related intronic region of the macaque orthologous exon (that is, performing no sequence alignment). We defined conservation of upstream pairing as the proportion of those exons with a conserved downstream motif that have also conserved the upstream partner. Conservation of downstream pairing was similarly defined. As a control, for each regional class we used two groups of paired motifs that are non-co-occurring. Non-co-occurring here signifies a motif that has the same GC content as the co-occurring motif, but that generates a P-value of >0.05 with its original partner. For example, for the non-co-occurring motif pair UpDpnon-pairing, Dpnon-pairing had the same GC content as the corresponding Dp but formed a non-co-occurring pair with Up; and for Upnon-pairingDp, Upnon-pairing had the same GC content as the corresponding Up but formed a non-co-occurring pair with Dp.
Human SNP data were downloaded from dbSNP build 130 . SNPs that mapped to multiple genomic regions or known repetitive elements were excluded. The density of SNPs was calculated separately for each class of motifs when co-occurring and when alone.
Tissue expression bias analysis
From BioGPS  we obtained the Affymetrix HG-133A and GNF1B microarray gene expression data from 79 human tissues and cell lines . Mappings for Affymetrix probe identifiers were downloaded from . Genes expressed in a given tissue or cell line at greater than two standard deviations above the median expression across all 79 tissues and cell lines were defined as tissue-specifically expressed. For each co-occurring motif pair, we identified all exons that have the pair in the related intronic regions, and the corresponding genes that contain these exons. For each tissue/cell line, the fraction of genes showing tissue-specific expression, FE,Tissue, was calculated. As a control, we collected the same number of randomly chosen exons that did not contain the co-occurring pair (10,000 iterations). The number of exons in each control gene was matched to the same number of exons in the genes with co-occurring motifs. The fraction of control genes associated with each tissue/cell line was then calculated, FC,Tissue. The P-value for the significance of enrichment (or depletion) for each tissue/cell line was computed for the fraction of times that FE,Tissue was lower than or equal to FC,Tissue (FE,Tissue was higher than or equal to FC,Tissue for depletion). The P-value cutoff for significance of tissue/cell lines bias was 0.001.
Experimental validation of co-occurring motifs
To examine the cooperative activity of co-occurring pairs in splicing, we used a hamster dhfr minigene modified from that previously described . The starting minigene contains a central merged exon 2-3 flanked by exon 1 and a merged exon 4-6. This mini-gene is driven by a cytomegalovirus (CMV) promoter and is terminated by the SV40 late poly A site. Around the middle exon, a NotI site has been inserted in intron 1 and a natural NheI site occurs in the second intron. We tested two exons as the middle exon: human beta globin gene exon 2 (Hb2), and the Wilms tumor gene 1 exon 5 (WT1-5) with two point mutations in its exon body. Hb2 and WT1-5 were amplified without their flanking intronic sequences beyond the splice site consensus sequences (upstream 14 nucleotides and downstream 6 nucleotides for Hb2, upstream 26 nucleotides and downstream 6 nucleotides for WT1-5) using primers tailed with NotI or NheI restriction sites. The amplified fragments were cut by NotI and NheI and inserted as the middle exon of the mini-gene construct. Pentamers U and D of a pair were inserted as tandem duplicates between the NotI site and the 3' splice site region and between 5' splice site region and the NheI site of the tested exons using synthetic double-stranded oligomers. Neutral non-co-occurring pentamer pairs were chosen so as neither to form co-occurring pairs nor to create any known splicing motifs after insertion (that is, also considering the overlaps formed). We used the merged hexamer set (hexESEs and hexESSs) from our previous study  and the ISRE list by Yeo and colleagues  as the known splicing motifs to avoid.
HEK 293 cells cultured in 35 mm dishes were transfected using Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) following the manufacturer's protocol. After cells were incubated for 24 hours, total RNA was extracted using illustra RNAspin Mini Kits (GE Healthcare, Piscataway, NJ, USA). A sample of 400 ng of RNA was reverse transcribed (RT) using Omniscript (Qiagen, Valencia, CA, USA) and a specific primer, AGAGTCTGAGATGGCCTGGCT, which pairs with a region in the third exon. One-tenth of the RT product was used as template in the following PCR amplification: forward primer, GTCAGATCCGCCTCCGCGTA; reverse primer, GTAAACGGAACTGCCTCCAA; initial denaturation, 94°C for 2 minutes; denaturation, 94°C for 45 s; annealing, 60°C, 1 minute; extension, 72°C, 1 minute; 20 cycles; final extension, 72°C, 5 minutes. Splicing products were separated in 1.8% agarose gels stained with ethidium bromide; the intensity of each splicing product was quantified with ImageJ. At least two independent transfections were performed for each construct. Proportion included (I) was defined as Included product/(Skipped product + Included product) in molar quantities.
Definition of synergy index (SI)
By analogy to previous synergy calculations [44, 45], we focused on the effect of our experimental perturbations on exon skipping. Proportion skipped (S) was defined as Skipped product/(Skipped product + Included product). For the insertion of motif U, skipping 'fitness' was defined as WU = SU/Sno-insertion. For motif pair U:D, we evaluated the level of synergy by comparing the skipping fitness WUD of the pair insertion construct with the product of the skipping fitness values WU and WD of the corresponding single motif insertion constructs. The synergy index for promoting splicing inclusion is defined as SI = WUD - WU × WD. SI <0 indicates synergy, SI = 0 indicates no synergy, and SI >0 indicates anti-synergy [44, 45].
base bias corrected co-occurrence
exonic splicing enhancer
exonic splicing silencer
expressed sequence tag
human beta globin gene
heterogeneous nuclear ribonucleoprotein
intronic splicing enhancer
exon 5 of the Wilm's tumor gene 1.
We thank Mauricio Arias for providing the framework minigene into which we inserted different central exons and Juan Valcarcel for suggesting the exon body swap experiment. We are grateful for many useful discussions with Dennis Weiss and with present and former laboratory members Mauricio Arias, Jon Cacciatore, Ron Gejman, Ashira Lubkin, Shulian Shang, and Xiang Zhang. This work was supported by a grant from the NIH (GM072740) to LAC.
- Robberson BL, Cote GJ, Berget SM: Exon definition may facilitate splice site selection in RNAs with multiple exons. Mol Cell Biol. 1990, 10: 84-94.PubMedPubMed CentralView ArticleGoogle Scholar
- O'Neill JP, Rogan PK, Cariello N, Nicklas JA: Mutations that alter RNA splicing of the human HPRT gene: a review of the spectrum. Mutat Res. 1998, 411: 179-214. 10.1016/S1383-5742(98)00013-1.PubMedView ArticleGoogle Scholar
- Sharma S, Kohlstaedt LA, Damianov A, Rio DC, Black DL: Polypyrimidine tract binding protein controls the transition from exon definition to an intron defined spliceosome. Nat Struct Mol Biol. 2008, 15: 183-191. 10.1038/nsmb.1375.PubMedPubMed CentralView ArticleGoogle Scholar
- House AE, Lynch KW: An exonic splicing silencer represses spliceosome assembly after ATP-dependent exon recognition. Nat Struct Mol Biol. 2006, 13: 937-944. 10.1038/nsmb1149.PubMedView ArticleGoogle Scholar
- Chabot B: Directing alternative splicing: cast and scenarios. Trends Genet. 1996, 12: 472-478. 10.1016/0168-9525(96)10037-8.PubMedView ArticleGoogle Scholar
- Reed R: Mechanisms of fidelity in pre-mRNA splicing. Curr Opin Cell Biol. 2000, 12: 340-345. 10.1016/S0955-0674(00)00097-1.PubMedView ArticleGoogle Scholar
- Carothers AM, Urlaub G, Grunberger D, Chasin LA: Splicing mutants and their second-site suppressors at the dihydrofolate reductase locus in Chinese hamster ovary cells. Mol Cell Biol. 1993, 13: 5085-5098.PubMedPubMed CentralView ArticleGoogle Scholar
- Chasin LA: Searching for splicing motifs. Adv Exp Med Biol. 2007, 623: 85-106. full_text.PubMedView ArticleGoogle Scholar
- Xiao X, Wang Z, Jang M, Burge CB: Coevolutionary networks of splicing cis-regulatory elements. Proc Natl Acad Sci USA. 2007, 104: 18583-18588. 10.1073/pnas.0707349104.PubMedPubMed CentralView ArticleGoogle Scholar
- Ke S, Zhang XH, Chasin LA: Positive selection acting on splicing motifs reflects compensatory evolution. Genome Res. 2008, 18: 533-543. 10.1101/gr.070268.107.PubMedPubMed CentralView ArticleGoogle Scholar
- Xiao X, Wang Z, Jang M, Nutiu R, Wang ET, Burge CB: Splice site strength-dependent activity and genetic buffering by poly-G runs. Nat Struct Mol Biol. 2009, 16: 1094-1100. 10.1038/nsmb.1661.PubMedPubMed CentralView ArticleGoogle Scholar
- Han K, Yeo G, An P, Burge CB, Grabowski PJ: A combinatorial code for splicing silencing: UAGG and GGGG motifs. PLoS Biol. 2005, 3: e158-10.1371/journal.pbio.0030158.PubMedPubMed CentralView ArticleGoogle Scholar
- Blanchette M, Chabot B: Modulation of exon skipping by high-affinity hnRNP A1-binding sites and by intron elements that repress splice site utilization. EMBO J. 1999, 18: 1939-1952. 10.1093/emboj/18.7.1939.PubMedPubMed CentralView ArticleGoogle Scholar
- Martinez-Contreras R, Fisette JF, Nasim FU, Madden R, Cordeau M, Chabot B: Intronic binding sites for hnRNP A/B and hnRNP F/H proteins stimulate pre-mRNA splicing. PLoS Biol. 2006, 4: e21-10.1371/journal.pbio.0040021.PubMedPubMed CentralView ArticleGoogle Scholar
- Friedman BA, Stadler MB, Shomron N, Ding Y, Burge CB: Ab initio identification of functionally interacting pairs of cis-regulatory elements. Genome Res. 2008, 18: 1643-1651. 10.1101/gr.080085.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang XH, Chasin LA: Computational definition of sequence motifs governing constitutive exon splicing. Genes Dev. 2004, 18: 1241-1250. 10.1101/gad.1195304.PubMedPubMed CentralView ArticleGoogle Scholar
- Fairbrother WG, Yeh RF, Sharp PA, Burge CB: Predictive identification of exonic splicing enhancers in human genes. Science. 2002, 297: 1007-1013. 10.1126/science.1073774.PubMedView ArticleGoogle Scholar
- Goren A, Ram O, Amit M, Keren H, Lev-Maor G, Vig I, Pupko T, Ast G: Comparative analysis identifies exonic splicing regulatory sequences - the complex definition of enhancers and silencers. Mol Cell. 2006, 22: 769-781. 10.1016/j.molcel.2006.05.008.PubMedView ArticleGoogle Scholar
- Yeo GW, Van Nostrand EL, Liang TY: Discovery and analysis of evolutionarily conserved intronic splicing regulatory elements. PLoS Genet. 2007, 3: e85-10.1371/journal.pgen.0030085.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang XH, Leslie CS, Chasin LA: Dichotomous splicing signals in exon flanks. Genome Res. 2005, 15: 768-779. 10.1101/gr.3217705.PubMedPubMed CentralView ArticleGoogle Scholar
- Voelker RB, Berglund JA: A comprehensive computational characterization of conserved mammalian intronic sequences reveals conserved motifs associated with constitutive and alternative splicing. Genome Res. 2007, 17: 1023-1033. 10.1101/gr.6017807.PubMedPubMed CentralView ArticleGoogle Scholar
- Craik CS, Sprang S, Fletterick R, Rutter WJ: Intron-exon splice junctions map at protein surfaces. Nature. 1982, 299: 180-182. 10.1038/299180a0.PubMedView ArticleGoogle Scholar
- Liu S, Altman RB: Large scale study of protein domain distribution in the context of alternative splicing. Nucleic Acids Res. 2003, 31: 4828-4835. 10.1093/nar/gkg668.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang XH, Heller KA, Hefter I, Leslie CS, Chasin LA: Sequence information for the splicing of human pre-mRNA identified by support vector machine classification. Genome Res. 2003, 13: 2637-2650. 10.1101/gr.1679003.PubMedPubMed CentralView ArticleGoogle Scholar
- Costantini M, Cammarano R, Bernardi G: The evolution of isochore patterns in vertebrate genomes. BMC Genomics. 2009, 10: 146-10.1186/1471-2164-10-146.PubMedPubMed CentralView ArticleGoogle Scholar
- Duret L, Hurst LD: The elevated GC content at exonic third sites is not evidence against neutralist models of isochore evolution. Mol Biol Evol. 2001, 18: 757-762.PubMedView ArticleGoogle Scholar
- Zhang XH, Leslie CS, Chasin LA: Computational searches for splicing signals. Methods. 2005, 37: 292-305. 10.1016/j.ymeth.2005.07.011.PubMedView ArticleGoogle Scholar
- Modafferi EF, Black DL: Combinatorial control of a neuron-specific exon. Rna. 1999, 5: 687-706. 10.1017/S1355838299990155.PubMedPubMed CentralView ArticleGoogle Scholar
- Ule J, Stefani G, Mele A, Ruggiu M, Wang X, Taneri B, Gaasterland T, Blencowe BJ, Darnell RB: An RNA map predicting Nova-dependent splicing regulation. Nature. 2006, 444: 580-586. 10.1038/nature05304.PubMedView ArticleGoogle Scholar
- Web supplement to: Intronic motif pairs cooperate across exons to promote pre-mRNA splicing. [http://www.columbia.edu/cu/biology/faculty/chasin/BBC_COOC]
- Berget SM: Exon recognition in vertebrate splicing. J Biol Chem. 1995, 270: 2411-2414.PubMedView ArticleGoogle Scholar
- Clark F, Thanaraj TA: Categorization and characterization of transcript-confirmed constitutively and alternatively spliced introns and exons from human. Hum Mol Genet. 2002, 11: 451-464. 10.1093/hmg/11.4.451.PubMedView ArticleGoogle Scholar
- Bradnam KR, Korf I: Longer first introns are a general property of eukaryotic gene structure. PLoS One. 2008, 3: e3093-10.1371/journal.pone.0003093.PubMedPubMed CentralView ArticleGoogle Scholar
- Martinez-Contreras R, Cloutier P, Shkreta L, Fisette JF, Revil T, Chabot B: hnRNP proteins and splicing control. Adv Exp Med Biol. 2007, 623: 123-147. full_text.PubMedView ArticleGoogle Scholar
- Wang Z, Rolish ME, Yeo G, Tung V, Mawson M, Burge CB: Systematic identification and analysis of exonic splicing silencers. Cell. 2004, 119: 831-845. 10.1016/j.cell.2004.11.010.PubMedView ArticleGoogle Scholar
- Wang Z, Xiao X, Van Nostrand E, Burge CB: General and specific functions of exonic splicing silencers in splicing control. Mol Cell. 2006, 23: 61-70. 10.1016/j.molcel.2006.05.018.PubMedPubMed CentralView ArticleGoogle Scholar
- Ule J, Ule A, Spencer J, Williams A, Hu JS, Cline M, Wang H, Clark T, Fraser C, Ruggiu M, Zeeberg BR, Kane D, Weinstein JN, Blume J, Darnell RB: Nova regulates brain-specific splicing to shape the synapse. Nat Genet. 2005, 37: 844-852. 10.1038/ng1610.PubMedView ArticleGoogle Scholar
- Zhang C, Zhang Z, Castle J, Sun S, Johnson J, Krainer AR, Zhang MQ: Defining the regulatory network of the tissue-specific splicing factors Fox-1 and Fox-2. Genes Dev. 2008, 22: 2550-2563. 10.1101/gad.1703108.PubMedPubMed CentralView ArticleGoogle Scholar
- Das D, Clark TA, Schweitzer A, Yamamoto M, Marr H, Arribere J, Minovitsky S, Poliakov A, Dubchak I, Blume JE, Conboy JG: A correlation with exon expression approach to identify cis-regulatory elements for tissue-specific alternative splicing. Nucleic Acids Res. 2007, 35: 4845-4857. 10.1093/nar/gkm485.PubMedPubMed CentralView ArticleGoogle Scholar
- Hanamura A, Caceres JF, Mayeda A, Franza BR, Krainer AR: Regulated tissue-specific expression of antagonistic pre-mRNA splicing factors. RNA. 1998, 4: 430-444.PubMedPubMed CentralGoogle Scholar
- Calarco JA, Superina S, O'Hanlon D, Gabut M, Raj B, Pan Q, Skalska U, Clarke L, Gelinas D, van der Kooy D, Zhen M, Ciruna B, Blencowe BJ: Regulation of vertebrate nervous system alternative splicing and development by an SR-related protein. Cell. 2009, 138: 898-910. 10.1016/j.cell.2009.06.012.PubMedView ArticleGoogle Scholar
- Venables JP, Koh CS, Froehlich U, Lapointe E, Couture S, Inkel L, Bramard A, Paquet ER, Watier V, Durand M, Lucier JF, Gervais-Bird J, Tremblay K, Prinos P, Klinck R, Elela SA, Chabot B: Multiple and specific mRNA processing targets for the major human hnRNP proteins. Mol Cell Biol. 2008, 28: 6033-6043. 10.1128/MCB.00726-08.PubMedPubMed CentralView ArticleGoogle Scholar
- Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G, Cooke MP, Walker JR, Hogenesch JB: A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci USA. 2004, 101: 6062-6067. 10.1073/pnas.0400782101.PubMedPubMed CentralView ArticleGoogle Scholar
- Segre D, Deluna A, Church GM, Kishony R: Modular epistasis in yeast metabolism. Nat Genet. 2005, 37: 77-83.PubMedGoogle Scholar
- Elena SF, Lenski RE: Test of synergistic interactions among deleterious mutations in bacteria. Nature. 1997, 390: 395-398. 10.1038/37108.PubMedView ArticleGoogle Scholar
- Krawczak M, Reiss J, Cooper DN: The mutational spectrum of single base-pair substitutions in mRNA splice junctions of human genes: causes and consequences. Hum Genet. 1992, 90: 41-54. 10.1007/BF00210743.PubMedView ArticleGoogle Scholar
- Unigene database. [ftp://ftp.ncbi.nih.gov/repository/UniGene/Homo_sapiens/Hs.seq.all.gz]
- Human genome sequence. [ftp://ftp.ncbi.nih.gov/genomes/H_sapiens/]
- RepeatMasker. [http://repeatmasker.org/]
- Shapiro MB, Senapathy P: RNA splice junctions of different classes of eukaryotes: sequence statistics and functional implications in gene expression. Nucleic Acids Res. 1987, 15: 7155-7174. 10.1093/nar/15.17.7155.PubMedPubMed CentralView ArticleGoogle Scholar
- Multiple genome alignment. [http://hgdownload.cse.ucsc.edu/goldenPath/hg17/multiz17way/]
- Human SNP database, dbSNP. [ftp://ftp.ncbi.nih.gov/snp/]
- BioGPS. [http://biogps.gnf.org/]
- Affymetrix platform GPL96. [http://www.ncbi.nlm.nih.gov/projects/geo/query/acc.cgi?acc=GPL96]
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.