Positive selection in alternatively spliced exons of human genes

Alternative splicing is a well-recognized mechanism of accelerated genome evolution. We have studied single-nucleotide polymorphisms and human-chimpanzee divergence in the exons of 6672 alternatively spliced human genes, with the aim of understanding the forces driving the evolution of alternatively spliced sequences. Here, we show that alternatively spliced exons and exon fragments (alternative exons) from minor isoforms experience lower selective pressure at the amino acid level, accompanied by selection against synonymous sequence variation. The results of the McDonald-Kreitman test suggest that alternatively spliced exons, unlike exons constitutively included in the mRNA, are also subject to positive selection, with up to 27% of amino acids fixed by positive selection.     

Different alternative splicing patterns are subject to opposite selection pressure for protein reading frame preservation

Background  

Alternative splicing (AS) has been regarded capable of altering selection pressure on protein subsequences. Particularly, the frequency of reading frame preservation (FRFP), as a measure of selection pressure, has been reported to be higher in alternatively spliced exons (ASEs) than in constitutively spliced exons (CSEs). However, recently it has been reported that different ASE types – simple and complex ASEs – may be subject to opposite selection forces. Therefore, it is necessary to re-evaluate the evolutionary effects of such splicing patterns on frame preservation.     

Combinatorial splicing of exon pairs by two-site binding of U1 small nuclear ribonucleoprotein particle.

A two-site model for the binding of U1 small nuclear ribonucleoprotein particle (U1 snRNP) was tested in order to understand how exon partners are selected in complex pre-mRNAs containing alternative exons. In this model, it is proposed that two U1 snRNPs define a functional unit of splicing by base pairing to the 3' boundary of the downstream exon as well as the 5' boundary of the intron to be spliced. Three-exon substrates contained the alternatively spliced exon 4 (E4) region of the preprotachykinin gene. Combined 5' splice site mutations at neighboring exons demonstrate that weakened binding of U1 snRNP at the downstream site and improved U1 snRNP binding at the upstream site result in the failure to rescue splicing of the intron between the mutations. These results indicate the stringency of the requirement for binding a second U1 snRNP to the downstream 5' splice site for these substrates as opposed to an alternative model in which a certain threshold level of U1 snRNP can be provided at either site. Further support for the two-site model is provided by single-site mutations in the 5' splice site of the third exon, E5, that weaken base complementarity to U1 RNA. These mutations block E5 branchpoint formation and, surprisingly, generate novel branchpoints that are specified chiefly by their proximity to a cryptic 5' splice site located at the 3' terminus of the pre-mRNA. The experiments shown here demonstrate a true stimulation of 3' splice site activity by the downstream binding of U1 snRNP and suggest a possible mechanism by which combinatorial patterns of exon selection are achieved for alternatively spliced pre-mRNAs.     

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.     

A nonsense exon in the Tpm1 gene is silenced by hnRNP H and F

As well as generating protein isoform diversity, in some cases alternative splicing generates RNAs that harbor premature termination codons and that are subject to nonsense-mediated decay (NMD). We previously identified an apparent pseudo-exon in the rat α-tropomyosin (Tpm1) gene as a probable genuine alternatively spliced exon that causes NMD when spliced into Tpm1 RNA. Here, we report the analysis of cis-acting splicing regulatory elements within this “nonsense exon.” Guided by the data set of predicted splicing enhancer and silencer elements compiled by Zhang and Chasin, we made a series of mutations through the nonsense exon and found that like authentic exons it is densely packed with enhancer and silencer elements. Strikingly, 11 of 13 tested mutations behaved as predicted computationally. In particular, we found that a G-rich silencer at the 5′ end, which is crucial for skipping of the nonsense exon, functions by binding hnRNP-H and F.     

Cis requirements for alternative splicing of the cardiac troponin T pre-mRNA.

The cardiac troponin T (cTNT) pre-mRNA splices 17 exons contiguously but alternatively splices (includes or excludes) the fifth exon. Because both alternative splice products are processed from the same pre-mRNA species, the cTNT pre-mRNA must contain cis-acting sequences which specify exon 5 as an alternative exon. A cTNT minigene (SM-1) transfected into cultured cells produces mRNAs both including and excluding exon 5. The junctions of exons 4-5-6 and 4-6 in the cTNT minigene mRNAs are identical to those of endogenous cTNT mRNAs and no other exons are alternatively spliced. Thus, the SM-1 pre-mRNA is correctly alternatively spliced in transfected cells. To circumscribe the pre-mRNA regions which are required for the alternative nature of exon 5, we have constructed a systematic series of deletion mutants of SM-1. Transfection of this series demonstrates that a 1200 nt pre-mRNA region containing exons 4, 5, and 6 is sufficient to direct alternative splicing of exon 5. Within this region are two relatively large inverted repeats which potentially sequester the alternative exon via intramolecular base-pairing. Such sequestration of an alternative exon is consistent with models which propose pre-mRNA conformation as being determinative for alternative splicing of some pre-mRNAs. However, deletion mutants which remove the majority of each of the inverted repeats retain the ability to alternatively splice exon 5 demonstrating that neither is required for cTNT alternative splice site selection. Taken together, deletion analysis has limited cis elements required for alternative splicing to three small regions of the pre-mRNA containing exons 4, 5, and 6. In addition, the cTNT minigene pre-mRNA expresses both alternative splice products in a wide variety of cultured non-muscle cells as well as in cultured striated muscle cells, although expression of the cTNT pre-mRNA is normally restricted to striated muscle. This indicates that cis elements involved in defining the cTNT exon 5 as an alternative exon do not require muscle-specific factors in trans to function.     

Alternative splicing and genomic stability

Alternative splicing allows an organism to make different proteins in different cells at different times, all from the same gene. In a cell that uses alternative splicing, the total length of all the exons is much shorter than in a cell that encodes the same set of proteins without alternative splicing. This economical use of exons makes genes more stable during reproduction and development because a genome with a shorter exon length is more resistant to harmful mutations. Genomic stability may be the reason why higher vertebrates splice alternatively. For a broad class of alternatively spliced genes, a formula is given for the increase in their stability.     

A test of translational selection at 'silent' sites in the human genome: base composition comparisons in alternatively spliced genes

Natural selection appears to discriminate among synonymous codons to enhance translational efficiency in a wide range of prokaryotes and eukaryotes. Codon bias is strongly related to gene expression levels in these species. In addition, between-gene variation in silent DNA divergence is inversely correlated with codon bias. However, in mammals, between-gene comparisons are complicated by distinctive nucleotide-content bias (isochores) throughout the genome. In this study, we attempted to identify translational selection by analyzing the DNA sequences of alternatively spliced genes in humans and in Drosophila melanogaster. Among codons in an alternatively spliced gene, those in constitutively expressed exons are translated more often than those in alternatively spliced exons. Thus, translational selection should act more strongly to bias codon usage and reduce silent divergence in constitutive than in alternative exons. By controlling for regional forces affecting base-composition evolution, this within-gene comparison makes it possible to detect codon selection at synonymous sites in mammals. We found that GC-ending codons are more abundant in constitutive than alternatively spliced exons in both Drosophila and humans. Contrary to our expectation, however, silent DNA divergence between mammalian species is higher in constitutive than in alternative exons.     

Intronic alternative splicing regulators identified by comparative genomics in nematodes

Many alternative splicing events are regulated by pentameric and hexameric intronic sequences that serve as binding sites for splicing regulatory factors. We hypothesized that intronic elements that regulate alternative splicing are under selective pressure for evolutionary conservation. Using a Wobble Aware Bulk Aligner genomic alignment of Caenorhabditis elegans and Caenorhabditis briggsae, we identified 147 alternatively spliced cassette exons that exhibit short regions of high nucleotide conservation in the introns flanking the alternative exon. In vivo experiments on the alternatively spliced let-2 gene confirm that these conserved regions can be important for alternative splicing regulation. Conserved intronic element sequences were collected into a dataset and the occurrence of each pentamer and hexamer motif was counted. We compared the frequency of pentamers and hexamers in the conserved intronic elements to a dataset of all C. elegans intron sequences in order to identify short intronic motifs that are more likely to be associated with alternative splicing. High-scoring motifs were examined for upstream or downstream preferences in introns surrounding alternative exons. Many of the high- scoring nematode pentamer and hexamer motifs correspond to known mammalian splicing regulatory sequences, such as (T)GCATG, indicating that the mechanism of alternative splicing regulation is well conserved in metazoans. A comparison of the analysis of the conserved intronic elements, and analysis of the entire introns flanking these same exons, reveals that focusing on intronic conservation can increase the sensitivity of detecting putative splicing regulatory motifs. This approach also identified novel sequences whose role in splicing is under investigation and has allowed us to take a step forward in defining a catalog of splicing regulatory elements for an organism. In vivo experiments confirm that one novel high-scoring sequence from our analysis, (T)CTATC, is important for alternative splicing regulation of the unc-52 gene.     

Finding signals that regulate alternative splicing in the post-genomic era

Alternative splicing of pre-mRNAs is central to the generation of diversity from the relatively small number of genes in metazoan genomes. Auxiliary cis elements and trans-acting factors are required for the recognition of constitutive and alternatively spliced exons and their inclusion in pre-mRNA. Here, we discuss the regulatory elements that direct alternative splicing and how genome-wide analyses can aid in their identification.     

Evolutionary divergence of exon flanks: a dissection of mutability and selection

The intronic sequences flanking exon-intron junctions (i.e., exon flanks) are important for splice site recognition and pre-mRNA splicing. Recent studies show a higher degree of sequence conservation at flanks of alternative exons, compared to flanks of constitutive exons. In this article we performed a detailed analysis on the evolutionary divergence of exon flanks between human and chimpanzee, aiming to dissect the impact of mutability and selection on their evolution. Inside exon flanks, sites that might reside in ancestral CpG dinucleotides evolved significantly faster than sites outside of ancestral CpG dinucleotides. This result reflects a systematic variation of mutation rates (mutability) at exon flanks, depending on the local CpG contexts. Remarkably, we observed a significant reduction of the nucleotide substitution rate in flanks of alternatively spliced exons, independent of the site-by-site variation in mutability due to different CpG contexts. Our data provide concrete evidence for increased purifying selection at exon flanks associated with regulation of alternative splicing.