Evidence for a trade-off between translational efficiency and splicing regulation in determining synonymous codon usage in Drosophila melanogaster

In Drosophila melanogaster, synonymous codons corresponding to the most abundant cognate tRNAs are used more frequently, especially in highly expressed genes. Increased use of such "optimal" codons is considered an adaptation for translational efficiency. Need it always be the case that selection should favor the use of a translationally optimal codon? Here, we investigate one possible confounding factor, namely, the need to specify information in exons necessary to enable correct splicing. As expected from such a model, in Drosophila many codons show different usage near intron-exon boundaries versus exon core regions. However, this finding is in principle also consistent with Hill-Robertson effects modulating usage of translationally optimal codons. However, several results support the splice model over the translational selection model: 1) the trends in codon usage are strikingly similar to those in mammals in which codon usage near boundaries correlates with abundance in exonic splice enhancers (ESEs), 2) codons preferred near boundaries tend to be enriched for A and avoid C (conversely those avoided near boundaries prefer C rather than A), as expected were ESEs involved, and 3) codons preferred near boundaries are typically not translationally optimal. We conclude that usage of translationally optimal codons usage is compromised in the vicinity of splice junctions in intron-containing genes, to the effect that we observe higher levels of usage of translationally optimal codons at the center of exons. On the gene level, however, controlling for known correlates of codon bias, the impact on codon usage patterns is quantitatively small. These results have implications for inferring aspects of the mechanism of splicing given nothing more than a well-annotated genome.     

Inference of splicing regulatory activities by sequence neighborhood analysis

Sequence-specific recognition of nucleic-acid motifs is critical to many cellular processes. We have developed a new and general method called Neighborhood Inference (NI) that predicts sequences with activity in regulating a biochemical process based on the local density of known sites in sequence space. Applied to the problem of RNA splicing regulation, NI was used to predict hundreds of new exonic splicing enhancer (ESE) and silencer (ESS) hexanucleotides from known human ESEs and ESSs. These predictions were supported by cross-validation analysis, by analysis of published splicing regulatory activity data, by sequence-conservation analysis, and by measurement of the splicing regulatory activity of 24 novel predicted ESEs, ESSs, and neutral sequences using an in vivo splicing reporter assay. These results demonstrate the ability of NI to accurately predict splicing regulatory activity and show that the scope of exonic splicing regulatory elements is substantially larger than previously anticipated. Analysis of orthologous exons in four mammals showed that the NI score of ESEs, a measure of function, is much more highly conserved above background than ESE primary sequence. This observation indicates a high degree of selection for ESE activity in mammalian exons, with surprisingly frequent interchangeability between ESE sequences.     

Mitochondrial-encoded endonucleases drive recombination of protein-coding genes in yeast

Mitochondrial recombination in yeast is well recognized, yet the underlying genetic mechanisms are not well understood. Recent progress has suggested that mobile introns in mitochondrial genomes (mitogenomes) can facilitate the recombination of their corresponding intron-containing genes through a mechanism known as intron homing. As many mitochondrial genes lack introns, there is a critical need to determine the extent of recombination and underlying mechanism of intron-lacking genes. This study leverages yeast mitogenomes to address these questions. In Saccharomyces cerevisiae, the 3′-end sequences of at least three intron-lacking mitochondrial genes exhibit elevated nucleotide diversity and recombination hotspots. Each of these 3′-end sequences is immediately adjacent to or even fused as overlapping genes with a stand-alone endonuclease. Our findings suggest that SAEs are responsible for recombination and elevated diversity of adjacent intron-lacking genes. SAEs were also evident to drive recombination of intron-lacking genes in Lachancea kluyveri, a yeast species that diverged from S. cerevisiae more than 100 million years ago. These results suggest SAEs as a common driver in recombination of intron-lacking genes during mitogenome evolution. We postulate that the linkage between intron-lacking gene and its adjacent endonuclease gene is the result of co-evolution.     

A strong exonic splicing enhancer in dystrophin exon 19 achieve proper splicing without an upstream polypyrimidine tract

Proper splicing is known to proceed under the control of conserved cis-elements located at exon-intron boundaries. Recently, it was shown that additional elements, such as exonic splicing enhancers (ESEs), are essential for the proper splicing of certain exons, in addition to the splice donor and acceptor site sequences; however, the relationship between these cis-elements is still unclear. In this report, we utilize dystrophin exon 19 to analyse the relationship between the ESE and its upstream acceptor site sequences. Dystrophin exon 19, which maintains adequate splicing donor and acceptor consensus sequences, encodes exonic splicing enhancer (dys-ESE19) sequences. Splice pattern analysis, using a minigene reporter expressed in HeLa cells, showed that either a strong polypyrimidine tract (PPT) or a fully active dys-ESE19 is sufficient for proper splicing. Each of these two cis-elements has enough activity for proper exon 19 splicing suggesting that the PPT, which is believed to be an essential cis-element for splicing, is dispensable when the downstream exon contains a strong ESE. This compensation was only seen in living cells but not in 'in vitro splicing'. This suggests the possibility that the previous splicing experiments using an in vitro splicing system could underestimate the activity of ESEs.     

Regulation of alternative RNA splicing by exon definition and exon sequences in viral and mammalian gene expression

Intron removal from a pre-mRNA by RNA splicing was once thought to be controlled mainly by intron splicing signals. However, viral and other eukaryotic RNA exon sequences have recently been found to regulate RNA splicing, polyadenylation, export, and nonsense-mediated RNA decay in addition to their coding function. Regulation of alternative RNA splicing by exon sequences is largely attributable to the presence of two majorcis-acting elements in the regulated exons, the exonic splicing enhancer (ESE) and the suppressor or silencer (ESS). Two types of ESEs have been verified from more than 50 genes or exons: purine-rich ESEs, which are the more common, and non-purine-rich ESEs. In contrast, the sequences of ESSs identified in approximately 20 genes or exons are highly diverse and show little similarity to each other. Through interactions with cellular splicing factors, an ESE or ESS determines whether or not a regulated splice site, usually an upstream 3 splice site, will be used for RNA splicing. However, how these elements function precisely in selecting a regulated splice site is only partially understood. The balance between positive and negative regulation of splice site selection likely depends on thecis-element's identity and changes in cellular splicing factors under physiological or pathological conditions.     

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.     

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.     

Exonic splicing enhancers: mechanism of action, diversity and role in human genetic diseases

Exonic splicing enhancers (ESEs) are discrete sequences within exons that promote both constitutive and regulated splicing. The precise mechanism by which ESEs facilitate the assembly of splicing complexes has been controversial. However, recent studies have provided insights into this question and have led to a new model for ESE function. Other recent work has suggested that ESEs are comprised of diverse sequences and occur frequently within exons. Ominously, these latter studies predict that many human genetic diseases linked to mutations within exons might be caused by the inactivation of ESEs.     

Human genomic sequences that inhibit splicing

Mammalian genes are characterized by relatively small exons surrounded by variable lengths of intronic sequence. Sequences similar to the splice signals that define the 5' and 3' boundaries of these exons are also present in abundance throughout the surrounding introns. What causes the real sites to be distinguished from the multitude of pseudosites in pre-mRNA is unclear. Much progress has been made in defining additional sequence elements that enhance the use of particular sites. Less work has been done on sequences that repress the use of particular splice sites. To find additional examples of sequences that inhibit splicing, we searched human genomic DNA libraries for sequences that would inhibit the inclusion of a constitutively spliced exon. Genetic selection experiments suggested that such sequences were common, and we subsequently tested randomly chosen restriction fragments of about 100 bp. When inserted into the central exon of a three-exon minigene, about one in three inhibited inclusion, revealing a high frequency of inhibitory elements in human DNA. In contrast, only 1 in 27 Escherichia coli DNA fragments was inhibitory. Several previously identified silencing elements derived from alternatively spliced exons functioned weakly in this constitutively spliced exon. In contrast, a high-affinity site for U2AF65 strongly inhibited exon inclusion. Together, our results suggest that splicing occurs in a background of repression and, since many of our inhibitors contain splice like signals, we suggest that repression of some pseudosites may occur through an inhibitory arrangement of these sites.     

Insights into the selective activation of alternatively used splice acceptors by the human immunodeficiency virus type-1 bidirectional splicing enhancer

The guanosine-adenosine-rich exonic splicing enhancer (GAR ESE) identified in exon 5 of the human immunodeficiency virus type-1 (HIV-1) pre-mRNA activates either an enhancer-dependent 5′ splice site (ss) or 3′ ss in 1-intron reporter constructs in the presence of the SR proteins SF2/ASF2 and SRp40. Characterizing the mode of action of the GAR ESE inside the internal HIV-1 exon 5 we found that this enhancer fulfils a dual splicing regulatory function (i) by synergistically mediating exon recognition through its individual SR protein-binding sites and (ii) by conferring 3′ ss selectivity within the 3′ ss cluster preceding exon 5. Both functions depend upon the GAR ESE, U1 snRNP binding at the downstream 5′ ss D4 and the E42 sequence located between these elements. Therefore, a network of cross-exon interactions appears to regulate splicing of the alternative exons 4a and 5. As the GAR ESE-mediated activation of the upstream 3′ ss cluster also is essential for the processing of intron-containing vpu/env-mRNAs during intermediate viral gene expression, the GAR enhancer substantially contributes to the regulation of viral replication.     

NF1 mRNA biogenesis: effect of the genomic milieu in splicing regulation of the NF1 exon 37 region

We have studied the splicing regulation of NF1 exons 36 and 37. We show that they not only require an intact exonic Splicing Enhancer (ESE) within exon 37, but also need the genomic region stretching from exons 31 to 38. Any nucleotide change in two exon 37 third codon positions disrupts the ESE. The extent of exons 36 and 37 skipping due to a mutated ESE depends on the genomic context. This is a unique example of what may be a more general phenomena involved in the tuning of pre-mRNA processing and gene expression modulation in the chromosomal setting.     

Linking C5 deficiency to an exonic splicing enhancer mutation

As an important component of the innate immune system, complement provides the initial response to prevent infections by pathogenic microorganisms. Patients with dysfunction of C5 display a propensity for severe recurrent infections. In this study, we present a patient with C5 deficiency demonstrated by immunochemical and functional analyses. Direct sequencing of all C5 exons displayed no mutation of obvious functional significance, except for an A to G transition in exon 10 predicting an exchange from lysine to arginine. This sequence alteration was present in only one allele of family members with a reduced serum C5 concentration and in both alleles of the patient with almost complete C5 deficiency, suggesting that this alteration may be producing the phenotype. Recent findings indicate that distinct nucleotide sequences, termed exonic splicing enhancers (ESEs), influence the splicing process. cDNA from all family members harboring the mutated allele showed skipping of exon 10, which resulted in a premature STOP codon, explaining the lack of C5 in the propositus. Sequence analysis of the mutated region revealed the substitution to be located within an ESE, as predicted by the RESCUE-ESE program. The altered ESE sequence is located close to the 5' splicing site and also lowers the predicted strength of the splice site itself. This apparently inconsequential sequence alteration represents a noncanonical splicing mutation altering an ESE. Our finding sheds a new light on the role of putative silent/conservative mutations in disease-associated genes.     

Elements of the rat tropoelastin gene associated with alternative splicing

Multiple isoforms of tropoelastin, the soluble precursor of elastin, are the products of translation of splice-variant mRNAs derived from the single-copy tropoelastin gene. Previous data had demonstrated DNA sequence heterogeneity in three domains of rat tropoelastin mRNA, indicating alternative splicing of several exons of the rat tropoelastin gene. Rat tropoelastin genomic clones encompassing the sites of alternative splicing were isolated and sequenced. Two sites of alternative splicing identified in rat tropoelastin mRNA sequences corresponded to exons 13-15 and exon 33 of the rat tropoelastin gene. Furthermore, the variable inclusion of an alanine codon in exon 16 resulted from two functional acceptor sites separated by three nucleotides. DNA sequences flanking exons subject to alternative splicing were analyzed. These exons contained splicing signals that differed from consensus sequences and from splicing signals of constitutively spliced exons. Introns immediately 5' of exons 14 and 33, for example, lacked typical polypyrimidine tracts and had weak, overlapping branch point sequences. Further, a region of secondary structure encompassing the acceptor site of exon 13 may influence alternative splicing of this exon. These results demonstrate that multiple cis-acting sequence elements may contribute to alternative splicing of rat tropoelastin pre-mRNA.     

Silencer elements as possible inhibitors of pseudoexon splicing

Human pre-mRNAs contain a definite number of exons and several pseudoexons which are located within intronic regions. We applied a computational approach to address the question of how pseudoexons are neglected in favor of exons and to possibly identify sequence elements preventing pseudoexon splicing. A search for possible splicing silencers was carried out on a pseudoexon selection that resembled exons in terms of splice site strength and exon splicing enhancer (ESE) representation; three motifs were retrieved through hexamer composition comparisons. One of these functions as a powerful silencer in transfection-based splicing assays and matches a previously identified silencer sequence with hnRNP H binding ability. The other two motifs are novel and failed to induce skipping of a constitutive exon, indicating that they might act as weak repressors or in synergy with other unidentified elements. All three motifs are enriched in pseudoexons compared with intronic regions and display higher frequencies in intronless gene-coding sequences compared with exons. We consider that a subpopulation of pseudoexons might rely on negative regulators for splicing repression; this hypothesis, if experimentally verified, might improve our understanding of exonic splicing regulatory sequences and provide the identification of a novel mutation target for human genetic diseases.     

ESEfinder: A web resource to identify exonic splicing enhancers

Point mutations frequently cause genetic diseases by disrupting the correct pattern of pre-mRNA splicing. The effect of a point mutation within a coding sequence is traditionally attributed to the deduced change in the corresponding amino acid. However, some point mutations can have much more severe effects on the structure of the encoded protein, for example when they inactivate an exonic splicing enhancer (ESE), thereby resulting in exon skipping. ESEs also appear to be especially important in exons that normally undergo alternative splicing. Different classes of ESE consensus motifs have been described, but they are not always easily identified. ESEfinder (http://exon.cshl.edu/ESE/) is a web-based resource that facilitates rapid analysis of exon sequences to identify putative ESEs responsive to the human SR proteins SF2/ASF, SC35, SRp40 and SRp55, and to predict whether exonic mutations disrupt such elements.     

Testing for Natural Selection in Human Exonic Splicing Regulators Associated with Evolutionary Rate Shifts

Despite evidence that at the interspecific scale, exonic splicing silencers (ESSs) are under negative selection in constitutive exons, little is known about the effects of slightly deleterious polymorphisms on these splicing regulators. Through the application of a modified version of the McDonald–Kreitman test, we compared the normalized proportions of human polymorphisms and human/rhesus substitutions affecting exonic splicing regulators (ESRs) on sequences of constitutive and alternative exons. Our results show a depletion of substitutions and an enrichment of SNPs associated with ESS gain in constitutive exons. Moreover, we show that this evolutionary pattern is also present in a set of ESRs previously involved in the transition from constitutive to skipped exons in the mammalian lineage. The similarity between these two sets of ESRs suggests that the transition from constitutive to skipped exons in mammals is more frequently associated with the inhibition than with the promotion of splicing signals. This is in accordance with the hypothesis of a constitutive origin of exon skipping and corroborates previous findings about the antagonistic role of certain exonic splicing enhancers.     

Alternative splicing of RNA triplets is often regulated and accelerates proteome evolution

Thousands of human genes contain introns ending in NAGNAG (N any nucleotide), where both NAGs can function as 3' splice sites, yielding isoforms that differ by inclusion/exclusion of three bases. However, few models exist for how such splicing might be regulated, and some studies have concluded that NAGNAG splicing is purely stochastic and nonfunctional. Here, we used deep RNA-Seq data from 16 human and eight mouse tissues to analyze the regulation and evolution of NAGNAG splicing. Using both biological and technical replicates to estimate false discovery rates, we estimate that at least 25% of alternatively spliced NAGNAGs undergo tissue-specific regulation in mammals, and alternative splicing of strongly tissue-specific NAGNAGs was 10 times as likely to be conserved between species as was splicing of non-tissue-specific events, implying selective maintenance. Preferential use of the distal NAG was associated with distinct sequence features, including a more distal location of the branch point and presence of a pyrimidine immediately before the first NAG, and alteration of these features in a splicing reporter shifted splicing away from the distal site. Strikingly, alignments of orthologous exons revealed a ～15-fold increase in the frequency of three base pair gaps at 3' splice sites relative to nearby exon positions in both mammals and in Drosophila. Alternative splicing of NAGNAGs in human was associated with dramatically increased frequency of exon length changes at orthologous exon boundaries in rodents, and a model involving point mutations that create, destroy, or alter NAGNAGs can explain both the increased frequency and biased codon composition of gained/lost sequence observed at the beginnings of exons. This study shows that NAGNAG alternative splicing generates widespread differences between the proteomes of mammalian tissues, and suggests that the evolutionary trajectories of mammalian proteins are strongly biased by the locations and phases of the introns that interrupt coding sequences.