A network of conserved co-occurring motifs for the regulation of alternative splicing

Cis-acting short sequence motifs play important roles in alternative splicing. It is now possible to identify such sequence motifs as conserved sequence patterns in genome sequence alignments. Here, we report the systematic search for motifs in the neighboring introns of alternatively spliced exons by using comparative analysis of mammalian genome alignments. We identified 11 conserved sequence motifs that might be involved in the regulation of alternative splicing. These motifs are not only significantly overrepresented near alternatively spliced exons, but they also co-occur with each other, thus, forming a network of cis-elements, likely to be the basis for context-dependent regulation. Based on this finding, we applied the motif co-occurrence to predict alternatively skipped exons. We verified exon skipping in 29 cases out of 118 predictions (25%) by EST and mRNA sequences in the databases. For the predictions not verified by the database sequences, we confirmed exon skipping in 10 additional cases by using both RT–PCR experiments and the publicly available RNA-Seq data. These results indicate that even more alternative splicing events will be found with the progress of large-scale and high-throughput analyses for various tissue samples and developmental stages.     

An intronic signal for alternative splicing in the human genome

An important level at which the expression of programmed cell death (PCD) genes is regulated is alternative splicing. Our previous work identified an intronic splicing regulatory element in caspase-2 (casp-2) gene. This 100-nucleotide intronic element, In100, consists of an upstream region containing a decoy 3' splice site and a downstream region containing binding sites for splicing repressor PTB. Based on the signal of In100 element in casp-2, we have detected the In100-like sequences as a family of sequence elements associated with alternative splicing in the human genome by using computational and experimental approaches. A survey of human genome reveals the presence of more than four thousand In100-like elements in 2757 genes. These In100-like elements tend to locate more frequent in intronic regions than exonic regions. EST analyses indicate that the presence of In100-like elements correlates with the skipping of their immediate upstream exons, with 526 genes showing exon skipping in such a manner. In addition, In100-like elements are found in several human caspase genes near exons encoding the caspase active domain. RT-PCR experiments show that these caspase genes indeed undergo alternative splicing in a pattern predicted to affect their functional activity. Together, these results suggest that the In100-like elements represent a family of intronic signals for alternative splicing in the human genome.     

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

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

Modulation of exon skipping and inclusion by heterogeneous nuclear ribonucleoprotein A1 and pre-mRNA splicing factor SF2/ASF.

The essential splicing factor SF2/ASF and the heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1) modulate alternative splicing in vitro of pre-mRNAs that contain 5' splice sites of comparable strengths competing for a common 3' splice site. Using natural and model pre-mRNAs, we have examined whether the ratio of SF2/ASF to hnRNP A1 also regulates other modes of alternative splicing in vitro. We found that an excess of SF2/ASF effectively prevents inappropriate exon skipping and also influences the selection of mutually exclusive tissue-specific exons in natural beta-tropomyosin pre-mRNA. In contrast, an excess of hnRNP A1 does not cause inappropriate exon skipping in natural constitutively or alternatively spliced pre-mRNAs. Although hnRNP A1 can promote alternative exon skipping, this effect is not universal and is dependent, e.g., on the size of the internal alternative exon and on the strength of the polypyrimidine tract in the preceding intron. With appropriate alternative exons, an excess of SF2/ASF promotes exon inclusion, whereas an excess of hnRNP A1 causes exon skipping. We propose that in some cases the ratio of SF2/ASF to hnRNP A1 may play a role in regulating alternative splicing by exon inclusion or skipping through the antagonistic effects of these proteins on alternative splice site selection.     

Differential GC content between exons and introns establishes distinct strategies of splice-site recognition

During evolution segments of homeothermic genomes underwent a GC content increase. Our analyses reveal that two exon-intron architectures have evolved from an ancestral state of low GC content exons flanked by short introns with a lower GC content. One group underwent a GC content elevation that abolished the differential exon-intron GC content, with introns remaining short. The other group retained the overall low GC content as well as the differential exon-intron GC content, and is associated with longer introns. We show that differential exon-intron GC content regulates exon inclusion level in this group, in which disease-associated mutations often lead to exon skipping. This group's exons also display higher nucleosome occupancy compared to flanking introns and exons of the other group, thus "marking" them for spliceosomal recognition. Collectively, our results reveal that differential exon-intron GC content is a previously unidentified determinant of exon selection and argue that the two GC content architectures reflect the two mechanisms by which splicing signals are recognized: exon definition and intron definition.     

Evidence for splice site pairing via intron definition in Schizosaccharomyces pombe

Schizosaccharomyces pombe pre-mRNAs are generally multi-intronic and share certain features with pre-mRNAs from Drosophila melanogaster, in which initial splice site pairing can occur via either exon or intron definition. Here, we present three lines of evidence suggesting that, despite these similarities, fission yeast splicing is most likely restricted to intron definition. First, mutating either or both splice sites flanking an internal exon in the S. pombe cdc2 gene produced almost exclusively intron retention, in contrast to the exon skipping observed in vertebrates. Second, we were unable to induce skipping of the internal microexon in fission yeast cgs2, whereas the default splicing pathway excludes extremely small exons in mammals. Because nearly quantitative removal of the downstream intron in cgs2 could be achieved by expanding the microexon, we propose that its retention is due to steric occlusion. Third, several cryptic 5' junctions in the second intron of fission yeast cdc2 are located within the intron, in contrast to their generally exonic locations in metazoa. The effects of expanding and contracting this intron are as predicted by intron definition; in fact, even highly deviant 5' junctions can compete effectively with the standard 5' splice site if they are closer to the 3' splicing signals. Taken together, our data suggest that pairing of splice sites in S. pombe most likely occurs exclusively across introns in a manner that favors excision of the smallest segment possible.     

Protein domains correlate strongly with exons in multiple eukaryotic genomes--evidence of exon shuffling?

We conducted a multi-genome analysis correlating protein domain organization with the exon-intron structure of genes in nine eukaryotic genomes. We observed a significant correlation between the borders of exons and domains on a genomic scale for both invertebrates and vertebrates. In addition, we found that the more complex organisms displayed consistently stronger exon-domain correlation, with substantially more significant correlations detected in vertebrates compared with invertebrates. Our observations concur with the principles of exon shuffling theory, including the prediction of predominantly symmetric phase of introns flanking the borders of correlating exons. These results suggest that extensive exon shuffling events during evolution significantly contributed to the shaping of eukaryotic proteomes.     

Conserved sequence elements associated with exon skipping

One of the major forms of alternative splicing, which generates multiple mRNA isoforms differing in the precise combinations of their exon sequences, is exon skipping. While in constitutive splicing all exons are included, in the skipped pattern(s) one or more exons are skipped. The regulation of this process is still not well understood; so far, cis- regulatory elements (such as exonic splicing enhancers) were identified in individual cases. We therefore set to investigate the possibility that exon skipping is controlled by sequences in the adjacent introns. We employed a computer analysis on 54 sequences documented as undergoing exon skipping, and identified two motifs both in the upstream and downstream introns of the skipped exons. One motif is highly enriched in pyrimidines (mostly C residues), and the other motif is highly enriched in purines (mostly G residues). The two motifs differ from the known cis-elements present at the 5′ and 3′ splice site. Interestingly, the two motifs are complementary, and their relative positional order is conserved in the flanking introns. These suggest that base pairing interactions can underlie a mechanism that involves secondary structure to regulate exon skipping. Remarkably, the two motifs are conserved in mouse orthologous genes that undergo exon skipping.     

Efficient internal exon recognition depends on near equal contributions from the 3' and 5' splice sites

Pre-mRNA splicing is carried out by the spliceosome, which identifies exons and removes intervening introns. In vertebrates, most splice sites are initially recognized by the spliceosome across the exon, because most exons are small and surrounded by large introns. This gene architecture predicts that efficient exon recognition depends largely on the strength of the flanking 3' and 5' splice sites. However, it is unknown if the 3' or the 5' splice site dominates the exon recognition process. Here, we test the 3' and 5' splice site contributions towards efficient exon recognition by systematically replacing the splice sites of an internal exon with sequences of different splice site strengths. We show that the presence of an optimal splice site does not guarantee exon inclusion and that the best predictor for exon recognition is the sum of both splice site scores. Using a genome-wide approach, we demonstrate that the combined 3' and 5' splice site strengths of internal exons provide a much more significant separator between constitutive and alternative exons than either the 3' or the 5' splice site strength alone.     

Organization and chromosomal localization of the human interleukin 5 receptor alpha-chain gene.

The gene for the hIL-5R alpha subunit, which is present in a single copy in the human genome, has been analysed in detail. It is located on chromosome 3 in the region 3p26. The gene organization reflects its relationship to the cytokine/haematopoietin receptor superfamily. Three introns are located in the 5' untranslated region. The subsequent exons determine the functional domains of the hIL-5R alpha protein: the signal peptide, three fibronectin type III-like (FN-like) modules, each built up by two exons, the membrane anchor and two exons forming the cytoplasmic tail, the first of which contains the proline cluster region. In addition, a specific exon generating a soluble isoform is located before the membrane anchor exon. This specific exon contains an in frame TAA stop codon, followed by a polyadenylation signal. Hence, a normal splicing event leads to a soluble IL-5R alpha variant, whereas alternative splicing is required for cell membrane anchoring. A second area of alternative splicing is found in the 5' leader sequence, and possibly relates to the presence of short open reading frames preceding the main ATG. All intron-exon junctions meet the GT-AG rule. The gene structures of all cytokine/haematopoietin receptors documented so far have also been compared with respect to intron phasing. This shows that all introns between the FN-III-like modules are of the +1 type, but in addition, splice sites within the Cys-module and WS-WS-module are invariably of the +2 and 0 type, respectively.     

Evolution of Nova-dependent splicing regulation in the brain

A large number of alternative exons are spliced with tissue-specific patterns, but little is known about how such patterns have evolved. Here, we study the conservation of the neuron-specific splicing factors Nova1 and Nova2 and of the alternatively spliced exons they regulate in mouse brain. Whereas Nova RNA binding domains are 94% identical across vertebrate species, Nova-dependent splicing silencer and enhancer elements (YCAY clusters) show much greater divergence, as less than 50% of mouse YCAY clusters are conserved at orthologous positions in the zebrafish genome. To study the relation between the evolution of tissue-specific splicing and YCAY clusters, we compared the brain-specific splicing of Nova-regulated exons in zebrafish, chicken, and mouse. The presence of YCAY clusters in lower vertebrates invariably predicted conservation of brain-specific splicing across species, whereas their absence in lower vertebrates correlated with a loss of alternative splicing. We hypothesize that evolution of Nova-regulated splicing in higher vertebrates proceeds mainly through changes in cis-acting elements, that tissue-specific splicing might in some cases evolve in a single step corresponding to evolution of a YCAY cluster, and that the conservation level of YCAY clusters relates to the functions encoded by the regulated RNAs.     

Intron definition in splicing of small Drosophila introns.

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