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

Multifactorial interplay controls the splicing profile of Alu-derived exons

Exonization of Alu elements creates primate-specific genomic diversity. Here we combine bioinformatic and experimental methodologies to reconstruct the molecular changes leading to exon selection. Our analyses revealed an intricate network involved in Alu exonization. A typical Alu element contains multiple sites with the potential to serve as 5' splice sites (5'ss). First, we demonstrated the role of 5'ss strength in controlling exonization events. Second, we found that a cryptic 5'ss enhances the selection of a more upstream site and demonstrate that this is mediated by binding of U1 snRNA to the cryptic splice site, challenging the traditional role attributed to U1 snRNA of binding the 5'ss only. Third, we used a simple algorithm to identify specific sequences that determine splice site selection within specific Alu exons. Finally, by inserting identical exons within different sequences, we demonstrated the importance of flanking genomic sequences in determining whether an Alu exon will undergo exonization. Overall, our results demonstrate the complex interplay between at least four interacting layers that affect Alu exonization. These results shed light on the mechanism through which Alu elements enrich the primate transcriptome and allow a better understanding of the exonization process in general.

Requirements for gene silencing mediated by U1 snRNA binding to a target sequence

U1 interference (U1i) is a novel method to block gene expression. U1i requires expression of a 5'-end-mutated U1 snRNA designed to base pair to the 3'-terminal exon of the target gene's pre-mRNA that leads to inhibition of polyadenylation. Here, we show U1i is robust (> or =95%) and a 10-nt target length is sufficient for good silencing. Surprisingly, longer U1 snRNAs, which could increase annealing to the target, fail to improve silencing. Extensive mutagenesis of the 10-bp U1 snRNA:target duplex shows that any single mismatch different from GU at positions 3-8, destroys silencing. However, mismatches within the other positions give partial silencing, suggesting that off-target inhibition could occur. The specificity of U1i may be enhanced, however, by the fact that silencing is impaired by RNA secondary structure or by splicing factors binding nearby, the latter mediated by Arginine-Serine (RS) domains. U1i inhibition can be reconstituted in vivo by tethering of RS domains of U1-70K and U2AF65. These results help to: (i) define good target sites for U1i; (ii) identify and understand natural cellular examples of U1i; (iii) clarify the contribution of hydrogen bonding to U1i and to U1 snRNP binding to 5' splice sites and (iv) understand the mechanism of U1i.

Alternative splicing of Alu exons--two arms are better than one

Alus, primate-specific retroelements, are the most abundant repetitive elements in the human genome. They are composed of two related but distinct monomers, left and right arms. Intronic Alu elements may acquire mutations that generate functional splice sites, a process called exonization. Most exonizations occur in right arms of antisense Alu elements, and are alternatively spliced. Here we show that without the left arm, exonization of the right arm shifts from alternative to constitutive splicing. This eliminates the evolutionary conserved isoform and may thus be selected against. We further show that insertion of the left arm downstream of a constitutively spliced non-Alu exon shifts splicing from constitutive to alternative. Although the two arms are highly similar, the left arm is characterized by weaker splicing signals and lower exonic splicing regulatory (ESR) densities. Mutations that improve these potential splice signals activate exonization and shift splicing from the right to the left arm. Collaboration between two or more putative splice signals renders the intronic left arm with a pseudo-exon function. Thus, the dimeric form of the Alu element fortuitously provides it with an evolutionary advantage, allowing enrichment of the primate transcriptome without compromising its original repertoire.

Positive selection acting on splicing motifs reflects compensatory evolution

We have used comparative genomics to characterize the evolutionary behavior of predicted splicing regulatory motifs. Using base substitution rates in intronic regions as a calibrator for neutral change, we found a strong avoidance of synonymous substitutions that disrupt predicted exonic splicing enhancers or create predicted exonic splicing silencers. These results attest to the functionality of the hexameric motif set used and suggest that they are subject to purifying selection. We also found that synonymous substitutions in constitutive exons tend to create exonic splicing enhancers and to disrupt exonic splicing silencers, implying positive selection for these splicing promoting events. We present evidence that this positive selection is the result of splicing-positive events compensating for splicing-negative events as well as for mutations that weaken splice-site sequences. Such compensatory events include nonsynonymous mutations, synonymous mutations, and mutations at splice sites. Compensation was also seen from the fact that orthologous exons tend to maintain the same number of predicted splicing motifs. Our data fit a splicing compensation model of exon evolution, in which selection for splicing-positive mutations takes place to counter the effect of an ongoing splicing-negative mutational process, with the exon as a whole being conserved as a unit of splicing. In the course of this analysis, we observed that synonymous positions in general are conserved relative to intronic sequences, suggesting that messenger RNA molecules are rich in sequence information for functions beyond protein coding and splicing.

Plant SR proteins and their functions

SR proteins are a family of splicing factors important for splice site recognition and spliceosome assembly. Their ability to bind to RNA and to interact with proteins as well identifies them as important players in splice site choice and alternative splicing. Plants possess twice as many SR proteins as animals, and some of the subfamilies are plant specific. Arabidopsis SR proteins are involved in different aspects of plant growth and development as well as in responses to environmental cues. The plant-specific subfamilies have been shown to be regulated by alternative splicing events, which are highly conserved in evolution. The tight regulation of splicing factors by alternative splicing might allow coordinated responses of their target genes.

Alternative approach to a heavy weight problem

Obesity is reaching epidemic proportions in developed countries and represents a significant risk factor for hypertension, heart disease, diabetes, and dyslipidemia. Splicing mutations constitute at least 14% of disease-causing mutations, thus implicating polymorphisms that affect splicing as likely candidates for disease susceptibility. A recent study suggested that genes associated with obesity were significantly enriched for rare nucleotide variants. Here, we examined these variants and revealed that they are located near splice junctions and tend to affect exonic splicing regulatory sequences. We also show that the majority of the exons that harbor these SNPs are constitutively spliced, yet they exhibit weak splice sites, typical to alternatively spliced exons, and are hence suboptimal for recognition by the splicing machinery and prone to become alternatively spliced. Using ex vivo assays, we tested a few representative variants and show that they indeed affect splicing by causing a shift from a constitutive to an alternative pattern, suggesting a possible link between extreme body mass index and abnormal splicing patterns.

Multisite and bidirectional exonic splicing enhancer in CD44 alternative exon v3

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

Biased exonization of transposed elements in duplicated genes: A lesson from the TIF-IA gene

BACKGROUND

Gene duplication and exonization of intronic transposed elements are two mechanisms that enhance genomic diversity. We examined whether there is less selection against exonization of transposed elements in duplicated genes than in single-copy genes.

RESULTS

Genome-wide analysis of exonization of transposed elements revealed a higher rate of exonization within duplicated genes relative to single-copy genes. The gene for TIF-IA, an RNA polymerase I transcription initiation factor, underwent a humanoid-specific triplication, all three copies of the gene are active transcriptionally, although only one copy retains the ability to generate the TIF-IA protein. Prior to TIF-IA triplication, an Alu element was inserted into the first intron. In one of the non-protein coding copies, this Alu is exonized. We identified a single point mutation leading to exonization in one of the gene duplicates. When this mutation was introduced into the TIF-IA coding copy, exonization was activated and the level of the protein-coding mRNA was reduced substantially. A very low level of exonization was detected in normal human cells. However, this exonization was abundant in most leukemia cell lines evaluated, although the genomic sequence is unchanged in these cancerous cells compared to normal cells.

CONCLUSION

The definition of the Alu element within the TIF-IA gene as an exon is restricted to certain types of cancers; the element is not exonized in normal human cells. These results further our understanding of the delicate interplay between gene duplication and alternative splicing and of the molecular evolutionary mechanisms leading to genetic innovations. This implies the existence of purifying selection against exonization in single copy genes, with duplicate genes free from such constrains.

Genome-wide analysis of alternative pre-mRNA splicing

Alternative splicing of mRNA precursors allows the synthesis of multiple mRNAs from a single primary transcript, significantly expanding the information content and regulatory possibilities of higher eukaryotic genomes. High-throughput enabling technologies, particularly large-scale sequencing and splicing-sensitive microarrays, are providing unprecedented opportunities to address key questions in this field. The picture emerging from these pioneering studies is that alternative splicing affects most human genes and a significant fraction of the genes in other multicellular organisms, with the potential to greatly influence the evolution of complex genomes. A combinatorial code of regulatory signals and factors can deploy physiologically coherent programs of alternative splicing that are distinct from those regulated at other steps of gene expression. Pre-mRNA splicing and its regulation play important roles in human pathologies, and genome-wide analyses in this area are paving the way for improved diagnostic tools and for the identification of novel and more specific pharmaceutical targets.

U2 snRNP: not just for poly(A) mRNAs

Polyadenylated mRNAs and replication-dependent histone mRNAs, which are not polyadenylated, share factors involved in 3' end formation. Friend et al. (2007) now report that the U2 snRNP, required for pre-mRNA splicing, is also required for histone mRNA 3' end formation.

Features generated for computational splice-site prediction correspond to functional elements

Background

Accurate selection of splice sites during the splicing of precursors to messenger RNA requires both relatively well-characterized signals at the splice sites and auxiliary signals in the adjacent exons and introns. We previously described a feature generation algorithm (FGA) that is capable of achieving high classification accuracy on human 3' splice sites. In this paper, we extend the splice-site prediction to 5' splice sites and explore the generated features for biologically meaningful splicing signals.

Results

We present examples from the observed features that correspond to known signals, both core signals (including the branch site and pyrimidine tract) and auxiliary signals (including GGG triplets and exon splicing enhancers). We present evidence that features identified by FGA include splicing signals not found by other methods.

Conclusion

Our generated features capture known biological signals in the expected sequence interval flanking splice sites. The method can be easily applied to other species and to similar classification problems, such as tissue-specific regulatory elements, polyadenylation sites, promoters, etc.

ADAM15 gene structure and differential alternative exon use in human tissues

Background

ADAM15 is a metalloprotease-disintegrin implicated in ectodomain shedding and cell adhesion. Aberrant ADAM15 expression has been associated with human cancer and other disorders. We have previously shown that the alternative splicing of ADAM15 transcripts is mis-regulated in cancer cells. To gain a better understanding of ADAM15 regulation, its genomic organization and regulatory elements as well as the alternative exon use in human tissues were characterized.

Results

Human ADAM15, flanked by the FLJ32785/DCST1 and ephrin-A4 genes, spans 11.4 kb from the translation initiation codon to the polyadenylation signal, being the shortest multiple-exon ADAM gene. The gene contains 23 exons varying from 63 to 316 bp and 22 introns from 79 to 1283 bp. The gene appeared to have several transcription start sites and their location suggested the promoter location within a CpG island proximal to the translation start. Reporter expression experiments confirmed the location of functional GC-rich, TATAless and CAATless promoter, with the most critical transcription-supporting elements located -266 to -23 bp relative to the translation start. Normal human tissues showed different complex patterns of at least 13 different ADAM15 splice variants arising from the alternative use of the cytosolic-encoding exons 19, 20a/b, and 21a/b. The deduced ADAM15 protein isoforms have different combinations of cytosolic regulatory protein interaction motifs.

Conclusion

Characterization of human ADAM15 gene and identification of elements involved in the regulation of transcription and alternative splicing provide important clues for elucidation of physiological and pathological roles of ADAM15. The present results also show that the alternative exon use is a physiological post-transcriptional mechanism regulating ADAM15 expression in human tissues.

Global control of aberrant splice-site activation by auxiliary splicing sequences: evidence for a gradient in exon and intron definition

Auxiliary splicing signals play a major role in the regulation of constitutive and alternative pre-mRNA splicing, but their relative importance in selection of mutation-induced cryptic or de novo splice sites is poorly understood. Here, we show that exonic sequences between authentic and aberrant splice sites that were activated by splice-site mutations in human disease genes have lower frequencies of splicing enhancers and higher frequencies of splicing silencers than average exons. Conversely, sequences between authentic and intronic aberrant splice sites have more enhancers and less silencers than average introns. Exons that were skipped as a result of splice-site mutations were smaller, had lower SF2/ASF motif scores, a decreased availability of decoy splice sites and a higher density of silencers than exons in which splice-site mutation activated cryptic splice sites. These four variables were the strongest predictors of the two aberrant splicing events in a logistic regression model. Elimination or weakening of predicted silencers in two reporters consistently promoted use of intron-proximal splice sites if these elements were maintained at their original positions, with their modular combinations producing expected modification of splicing. Together, these results show the existence of a gradient in exon and intron definition at the level of pre-mRNA splicing and provide a basis for the development of computational tools that predict aberrant splicing outcomes.

Regulation of retention of FosB intron 4 by PTB

One effect of stressors such as chronic drug administration is that sequence within the terminal exon of the transcription factor FosB is recognized as intronic and removed by alternative splicing. This results in an open-reading-frame shift that produces a translation stop codon and ultimately a truncated protein, termed DeltaFosB. In vitro splicing assays with control and mutated transcripts generated from a fosB mini-gene construct indicated a CU-rich sequence at the 3' end of intron 4 (I4) plays an important role in regulating fosB pre-mRNA splicing due to its binding of polypyrimidine tract binding protein (PTB). PTB binding to this sequence is dependent upon phosphorylation by protein kinase A and is blocked if the CU-rich sequence is mutated to a U-rich region. When this mutated fosB minigene is expressed in HeLa cells, the splicing efficiency of its product is increased compared to wild type. Moreover, transient transfection of PTB-1 in HeLa cells decreased the splicing efficiency of a wild type fosB minigene transcript. Depletion of PTB from nuclear extracts facilitated U2AF65 binding to wild type sequence in vitro, suggesting these proteins function in a dynamic equilibrium to modulate fosB pre-mRNA alternative splicing. These results demonstrate for the first time that phosphorylated PTB promotes intron retention and thereby silences the splicing of fosB I4.

Splicing in disease: disruption of the splicing code and the decoding machinery

Human genes contain a dense array of diverse cis-acting elements that make up a code required for the expression of correctly spliced mRNAs. Alternative splicing generates a highly dynamic human proteome through networks of coordinated splicing events. Cis- and trans-acting mutations that disrupt the splicing code or the machinery required for splicing and its regulation have roles in various diseases, and recent studies have provided new insights into the mechanisms by which these effects occur. An unexpectedly large fraction of exonic mutations exhibit a primary pathogenic effect on splicing. Furthermore, normal genetic variation significantly contributes to disease severity and susceptibility by affecting splicing efficiency.

Towards realistic codon models: among site variability and dependency of synonymous and non-synonymous rates

Codon evolutionary models are widely used to infer the selection forces acting on a protein. The non-synonymous to synonymous rate ratio (denoted by Ka/Ks) is used to infer specific positions that are under purifying or positive selection. Current evolutionary models usually assume that only the non-synonymous rates vary among sites while the synonymous substitution rates are constant. This assumption ignores the possibility of selection forces acting at the DNA or mRNA levels. Towards a more realistic description of sequence evolution, we present a model that accounts for among-site-variation of both synonymous and non-synonymous substitution rates. Furthermore, we alleviate the widespread assumption that positions evolve independently of each other. Thus, possible sources of bias caused by random fluctuations in either the synonymous or non-synonymous rate estimations at a single site is removed. Our model is based on two hidden Markov models that operate on the spatial dimension: one describes the dependency between adjacent non-synonymous rates while the other describes the dependency between adjacent synonymous rates. The presented model is applied to study the selection pressure across the HIV-1 genome. The new model better describes the evolution of all HIV-1 genes, as compared to current codon models. Using both simulations and real data analyses, we illustrate that accounting for synonymous rate variability and dependency greatly increases the accuracy of Ka/Ks estimation and in particular of positively selected sites. Finally, we discuss the applicability of the developed model to infer the selection forces in regulatory and overlapping regions of the HIV-1 genome.

Role of purine-rich exonic splicing enhancers in nuclear retention of pre-mRNAs

Intron-containing pre-mRNAs are normally retained in the nucleus until they are spliced to produce mature mRNAs that are exported to the cytoplasm. Although the detailed mechanism is not well understood, the formation of splicing-related complexes on pre-mRNAs is thought to be responsible for the nuclear retention. Therefore, pre-mRNAs containing suboptimal splice sites should tend to leak out to the cytoplasm. Such pre-mRNAs often contain purine-rich exonic splicing enhancers (ESEs) that stimulate splicing of the adjacent intron. Here, we show that ESEs per se possess an activity to retain RNAs in the nucleus through a saturable nuclear retention factor. Cross-competition experiments revealed that intron-containing pre-mRNAs (without ESEs) used the same saturable nuclear retention factor as ESEs. Interestingly, although intronless mRNAs containing ESEs were also poorly exported, spliced mRNAs produced from ESE-containing pre-mRNAs were efficiently exported to the cytoplasm. Thus, the splicing reaction can reset the nuclear retention state caused by ESEs, allowing nuclear export of mature mRNAs. Our results reveal a novel aspect of ESE activity that should contribute to gene expression and RNA quality control.

A comprehensive computational characterization of conserved mammalian intronic sequences reveals conserved motifs associated with constitutive and alternative splicing

Orthologous mammalian introns contain many highly conserved sequences. Of these sequences, many are likely to represent protein binding sites that are under strong positive selection. In order to identify conserved protein binding sites that are important for splicing, we analyzed the composition of intronic sequences that are conserved between human and six eutherian mammals. We focused on all completely conserved sequences of seven or more nucleotides located in the regions adjacent to splice-junctions. We found that these conserved intronic sequences are enriched in specific motifs, and that many of these motifs are statistically associated with either alternative or constitutive splicing. In validation of our methods, we identified several motifs that are known to play important roles in alternative splicing. In addition, we identified several novel motifs containing GCT that are abundant and are associated with alternative splicing. Furthermore, we demonstrate that, for some of these motifs, conservation is a strong indicator of potential functionality since conserved instances are associated with alternative splicing while nonconserved instances are not. A surprising outcome of this analysis was the identification of a large number of AT-rich motifs that are strongly associated with constitutive splicing. Many of these appear to be novel and may represent conserved intronic splicing enhancers (ISEs). Together these data show that conservation provides important insights into the identification and possible roles of cis-acting intronic sequences important for alternative and constitutive splicing.

The emergence of alternative 3' and 5' splice site exons from constitutive exons

Alternative 3' and 5' splice site (ss) events constitute a significant part of all alternative splicing events. These events were also found to be related to several aberrant splicing diseases. However, only few of the characteristics that distinguish these events from alternative cassette exons are known currently. In this study, we compared the characteristics of constitutive exons, alternative cassette exons, and alternative 3'ss and 5'ss exons. The results revealed that alternative 3'ss and 5'ss exons are an intermediate state between constitutive and alternative cassette exons, where the constitutive side resembles constitutive exons, and the alternative side resembles alternative cassette exons. The results also show that alternative 3'ss and 5'ss exons exhibit low levels of symmetry (frame-preserving), similar to constitutive exons, whereas the sequence between the two alternative splice sites shows high symmetry levels, similar to alternative cassette exons. In addition, flanking intronic conservation analysis revealed that exons whose alternative splice sites are at least nine nucleotides apart show a high conservation level, indicating intronic participation in the regulation of their splicing, whereas exons whose alternative splice sites are fewer than nine nucleotides apart show a low conservation level. Further examination of these exons, spanning seven vertebrate species, suggests an evolutionary model in which the alternative state is a derivative of an ancestral constitutive exon, where a mutation inside the exon or along the flanking intron resulted in the creation of a new splice site that competes with the original one, leading to alternative splice site selection. This model was validated experimentally on four exons, showing that they indeed originated from constitutive exons that acquired a new competing splice site during evolution.

SplicePort--an interactive splice-site analysis tool

SplicePort is a web-based tool for splice-site analysis that allows the user to make splice-site predictions for submitted sequences. In addition, the user can also browse the rich catalog of features that underlies these predictions, and which we have found capable of providing high classification accuracy on human splice sites. Feature selection is optimized for human splice sites, but the selected features are likely to be predictive for other mammals as well. With our interactive feature browsing and visualization tool, the user can view and explore subsets of features used in splice-site prediction (either the features that account for the classification of a specific input sequence or the complete collection of features). Selected feature sets can be searched, ranked or displayed easily. The user can group features into clusters and frequency plot WebLogos can be generated for each cluster. The user can browse the identified clusters and their contributing elements, looking for new interesting signals, or can validate previously observed signals. The SplicePort web server can be accessed at http://www.cs.umd.edu/projects/SplicePort and http://www.spliceport.org.

Identification of common genetic variation that modulates alternative splicing

Alternative splicing of genes is an efficient means of generating variation in protein function. Several disease states have been associated with rare genetic variants that affect splicing patterns. Conversely, splicing efficiency of some genes is known to vary between individuals without apparent ill effects. What is not clear is whether commonly observed phenotypic variation in splicing patterns, and hence potential variation in protein function, is to a significant extent determined by naturally occurring DNA sequence variation and in particular by single nucleotide polymorphisms (SNPs). In this study, we surveyed the splicing patterns of 250 exons in 22 individuals who had been previously genotyped by the International HapMap Project. We identified 70 simple cassette exon alternative splicing events in our experimental system; for six of these, we detected consistent differences in splicing pattern between individuals, with a highly significant association between splice phenotype and neighbouring SNPs. Remarkably, for five out of six of these events, the strongest correlation was found with the SNP closest to the intron-exon boundary, although the distance between these SNPs and the intron-exon boundary ranged from 2 bp to greater than 1,000 bp. Two of these SNPs were further investigated using a minigene splicing system, and in each case the SNPs were found to exert cis-acting effects on exon splicing efficiency in vitro. The functional consequences of these SNPs could not be predicted using bioinformatic algorithms. Our findings suggest that phenotypic variation in splicing patterns is determined by the presence of SNPs within flanking introns or exons. Effects on splicing may represent an important mechanism by which SNPs influence gene function.

Factors associated with a purine-rich exonic splicing enhancer sequence in Xenopus oocyte nucleus

Purine-rich exonic splicing enhancers (ESEs) stimulate splicing of the adjacent introns with suboptimal splice sites. To elucidate the mechanism regarding ESEs, factors specifically associated with ESEs in HeLa cell nuclear extracts were previously investigated, and shown to include SR (serine/arginine-rich) proteins. However, factors associated with ESEs in vivo have not yet been explored. Here we show that a GAA repeat RNA sequence, a typical ESE, is associated in Xenopus oocyte nuclei with at least one SR protein, SF2/ASF, as was expected. Moreover, components of SF3a/b complexes, U2 snRNA, and U2AF(65) were also found to be associated with the ESE in the nucleus. Since SF3a/b complexes are the constituents of the 17S U2 snRNP, these results suggest that the 17S U2 snRNP is associated with the ESE in the nucleus, probably through bridging interactions of U2AF and SR proteins. The identified factors may represent a functional splicing enhancer complex in vivo.

A computational survey of candidate exonic splicing enhancer motifs in the model plant Arabidopsis thaliana

BACKGROUND

Algorithmic approaches to splice site prediction have relied mainly on the consensus patterns found at the boundaries between protein coding and non-coding regions. However exonic splicing enhancers have been shown to enhance the utilization of nearby splice sites.

RESULTS

We have developed a new computational technique to identify significantly conserved motifs involved in splice site regulation. First, 84 putative exonic splicing enhancer hexamers are identified in Arabidopsis thaliana. Then a Gibbs sampling program called ELPH was used to locate conserved motifs represented by these hexamers in exonic regions near splice sites in confirmed genes. Oligomers containing 35 of these motifs have been shown experimentally to induce significant inclusion of A. thaliana exons. Second, integration of our regulatory motifs into two different splice site recognition programs significantly improved the ability of the software to correctly predict splice sites in a large database of confirmed genes. We have released GeneSplicerESE, the improved splice site recognition code, as open source software.

CONCLUSION

Our results show that the use of the ESE motifs consistently improves splice site prediction accuracy.

Discovery and analysis of evolutionarily conserved intronic splicing regulatory elements

Knowledge of the functional cis-regulatory elements that regulate constitutive and alternative pre-mRNA splicing is fundamental for biology and medicine. Here we undertook a genome-wide comparative genomics approach using available mammalian genomes to identify conserved intronic splicing regulatory elements (ISREs). Our approach yielded 314 ISREs, and insertions of ~70 ISREs between competing splice sites demonstrated that 84% of ISREs altered 5′ and 94% altered 3′ splice site choice in human cells. Consistent with our experiments, comparisons of ISREs to known splicing regulatory elements revealed that 40%–45% of ISREs might have dual roles as exonic splicing silencers. Supporting a role for ISREs in alternative splicing, we found that 30%–50% of ISREs were enriched near alternatively spliced (AS) exons, and included almost all known binding sites of tissue-specific alternative splicing factors. Further, we observed that genes harboring ISRE-proximal exons have biases for tissue expression and molecular functions that are ISRE-specific. Finally, we discovered that for Nova1, neuronal PTB, hnRNP C, and FOX1, the most frequently occurring ISRE proximal to an alternative conserved exon in the splicing factor strongly resembled its own known RNA binding site, suggesting a novel application of ISRE density and the propensity for splicing factors to auto-regulate to associate RNA binding sites to splicing factors. Our results demonstrate that ISREs are crucial building blocks in understanding general and tissue-specific AS regulation and the biological pathways and functions regulated by these AS events.

During RNA splicing, sequences (introns) in a pre-mRNA are excised and discarded, and the remaining sequences (exons) are joined to form the mature RNA. Splicing is regulated not only by the binding of the basic splicing machinery to splice sites located at the exon–intron boundaries, but also by the combined effects of various other splicing factors that bind to a multitude of sequence elements located both in the exons as well as the flanking introns. Instances of alternative splicing, where usage of splice site(s) is incomplete or different between tissues, cell types, or lineages, can be created by the interaction of sequence elements and tissue, cell type, and stage-specific splicing factors. To better understand constitutive and alternative pre-mRNA splicing, the authors describe a comparative genomics approach, using available mammalian genomes, to systematically identify splicing regulatory elements located in the introns proximal to exons. A quarter of the elements were tested experimentally, and most of them altered splicing in human cells. The authors also showed that that the intronic elements are close to tissue-specific alternative exons and are more likely to be located in specific positions in the introns, suggestive of potential regulatory function. These elements are also frequently found in tissue-specific genes, suggesting a coupling between expression and alternative splicing of these genes. Finally, the authors propose a strategy using the elements to identify the binding sites of several splicing factors.

Depolarization and CaM kinase IV modulate NMDA receptor splicing through two essential RNA elements

Alternative splicing controls the activity of many proteins important for neuronal excitation, but the signal-transduction pathways that affect spliced isoform expression are not well understood. One particularly interesting system of alternative splicing is exon 21 (E21) of the NMDA receptor 1 (NMDAR1 E21), which controls the trafficking of NMDA receptors to the plasma membrane and is repressed by Ca(++)/calmodulin-dependent protein kinase (CaMK) IV signaling. Here, we characterize the splicing of NMDAR1 E21. We find that E21 splicing is reversibly repressed by neuronal depolarization, and we identify two RNA elements within the exon that function together to mediate the inducible repression. One of these exonic elements is similar to an intronic CaMK IV-responsive RNA element (CaRRE) originally identified in the 3' splice site of the BK channel STREX exon, but not previously observed within an exon. The other element is a new RNA motif. Introduction of either of these two motifs, called CaRRE type 1 and CaRRE type 2, into a heterologous constitutive exon can confer CaMK IV-dependent repression on the new exon. Thus, either exonic CaRRE can be sufficient for CaMK IV-induced repression. Single nucleotide scanning mutagenesis defined consensus sequences for these two CaRRE motifs. A genome-wide motif search and subsequent RT-PCR validation identified a group of depolarization-regulated alternative exons carrying CaRRE consensus sequences. Many of these exons are likely to alter neuronal function. Thus, these two RNA elements define a group of co-regulated splicing events that respond to a common stimulus in neurons to alter their activity.

Improving the Caenorhabditis elegans genome annotation using machine learning

For modern biology, precise genome annotations are of prime importance, as they allow the accurate definition of genic regions. We employ state-of-the-art machine learning methods to assay and improve the accuracy of the genome annotation of the nematode Caenorhabditis elegans. The proposed machine learning system is trained to recognize exons and introns on the unspliced mRNA, utilizing recent advances in support vector machines and label sequence learning. In 87% (coding and untranslated regions) and 95% (coding regions only) of all genes tested in several out-of-sample evaluations, our method correctly identified all exons and introns. Notably, only 37% and 50%, respectively, of the presently unconfirmed genes in the C. elegans genome annotation agree with our predictions, thus we hypothesize that a sizable fraction of those genes are not correctly annotated. A retrospective evaluation of the Wormbase WS120 annotation [] of C. elegans reveals that splice form predictions on unconfirmed genes in WS120 are inaccurate in about 18% of the considered cases, while our predictions deviate from the truth only in 10%-13%. We experimentally analyzed 20 controversial genes on which our system and the annotation disagree, confirming the superiority of our predictions. While our method correctly predicted 75% of those cases, the standard annotation was never completely correct. The accuracy of our system is further corroborated by a comparison with two other recently proposed systems that can be used for splice form prediction: SNAP and ExonHunter. We conclude that the genome annotation of C. elegans and other organisms can be greatly enhanced using modern machine learning technology.

Widespread mRNA polyadenylation events in introns indicate dynamic interplay between polyadenylation and splicing

mRNA polyadenylation and pre-mRNA splicing are two essential steps for the maturation of most human mRNAs. Studies have shown that some genes generate mRNA variants involving both alternative polyadenylation and alternative splicing. Polyadenylation in introns can lead to conversion of an internal exon to a 3' terminal exon, which is termed composite terminal exon, or usage of a 3' terminal exon that is otherwise skipped, which is termed skipped terminal exon. Using cDNA/EST and genome sequences, we identified polyadenylation sites in introns for all currently known human genes. We found that approximately 20% human genes have at least one intronic polyadenylation event that can potentially lead to mRNA variants, most of which encode different protein products. The conservation of human intronic poly(A) sites in mouse and rat genomes is lower than that of poly(A) sites in 3'-most exons. Quantitative analysis of a number of mRNA variants generated by intronic poly(A) sites suggests that the intronic polyadenylation activity can vary under different cellular conditions for most genes. Furthermore, we found that weak 5' splice site and large intron size are the determining factors controlling the usage of composite terminal exon poly(A) sites, whereas skipped terminal exon poly(A) sites tend to be associated with strong polyadenylation signals. Thus, our data indicate that dynamic interplay between polyadenylation and splicing leads to widespread polyadenylation in introns and contributes to the complexity of transcriptome in the cell.

Functional and mechanistic insights from genome-wide studies of splicing regulation in the brain

We review here results arising from the systematic functional analysis of Nova, a neuron-specific RNA binding protein targeted in an autoimmune neurological disorder associated with cancer. We have developed a combination of biochemical, genetic and bioinformatic methods to generate a global understanding of Nova's role as a splicing regulator. Genome-wide identification and validation of Nova target RNAs has yielded unexpected insights into the protein's mechanism of action and into the functionally coherent role of Nova in the biology of the neuronal synapse. These studies provide us with a platform for understanding the role of RNA binding proteins in tissue-specific splicing regulation and in disease.

RNA-editing-mediated exon evolution

BACKGROUND

Alu retroelements are specific to primates and abundant in the human genome. Through mutations that create functional splice sites within intronic Alus, these elements can become new exons in a process denoted exonization. It was recently shown that Alu elements are also heavily changed by RNA editing in the human genome.

RESULTS

Here we show that the human nuclear prelamin A recognition factor contains a primate-specific Alu-exon that exclusively depends on RNA editing for its exonization. We demonstrate that RNA editing regulates the exonization in a tissue-dependent manner, through both the creation of a functional AG 3' splice site, and alteration of functional exonic splicing enhancers within the exon. Furthermore, a premature stop codon within the Alu-exon is eliminated by an exceptionally efficient RNA editing event. The sequence surrounding this editing site is important not only for editing of that site but also for editing in other neighboring sites as well.

CONCLUSION

Our results show that the abundant RNA editing of Alu sequences can be recruited as a mechanism supporting the birth of new exons in the human genome.

Alternative splicing of pre-messenger RNAs in plants in the genomic era

Primary transcripts (precursor-mRNAs) with introns can undergo alternative splicing to produce multiple transcripts from a single gene by differential use of splice sites, thereby increasing the transcriptome and proteome complexity within and between cells and tissues. Alternative splicing in plants is largely an unexplored area of gene expression, as this phenomenon used to be considered rare. However, recent genome-wide computational analyses have revealed that alternative splicing in flowering plants is far more prevalent than previously thought. Interestingly, pre-mRNAs of many spliceosomal proteins, especially serine/arginine-rich (SR) proteins, are extensively alternatively spliced. Furthermore, stresses have a dramatic effect on alternative splicing of pre-mRNAs including those that encode many spliceosomal proteins. Although the mechanisms that regulate alternative splicing in plants are largely unknown, several reports strongly suggest a key role for SR proteins in spliceosome assembly and regulated splicing. Recent studies suggest that alternative splicing in plants is an important posttranscriptional regulatory mechanism in modulating gene expression and eventually plant form and function.

SR proteins and related factors in alternative splicing

SR proteins are a family of RNA binding proteins that contain a signature RS domain enriched with serine/arginine repeats. The RS domain is also found in many other proteins, which are collectively referred to as SR-related proteins. Several prototypical SR proteins are essential splicing factors, but the majority of RS domain-containing factors are characterized by their ability to alter splice site selection in vitro or in transfected cells. SR proteins and SR-related proteins are generally believed to modulate splice site selection via RNA recognition motif-mediated binding to exonic splicing enhancers and RS domain-mediated protein-protein and protein-RNA interactions during spliceosome assembly. However, the biological function of individual RS domain-containing splicing regulators is complex because of redundant as well as competitive functions, context-dependent effects and regulation by cotranscriptional and post-translational events. This chapter will focus on our current mechanistic understanding of alternative splicing regulation by SR proteins and SR-related proteins and will discuss some of the questions that remain to be addressed in future research.

Reduced splicing efficiency induced by synonymous substitutions may generate a substrate for natural selection of new splicing isoforms: the case of CFTR exon 12

Alternative splicing has been associated with increased evolutionary changes and with recent exon creation or loss. The addition of a new exon can be explained by its inclusion in only a fraction of the transcripts leaving the original form intact and giving to the new form the possibility to evolve independently but the exon loss phenomenon is less clear. To explore the mechanism that could be involved in CFTR exon 12 lower splicing efficiency in primates, we have analyzed the effect of multiple synonymous variations. Random patterns of synonymous variations were created in CFTR exon12 and the majority of them induced exon inclusion, suggesting a suboptimal splicing efficiency of the human gene. In addition, the effect of each single synonymous substitution on splicing is strongly dependent on the exonic context and does not correlate with available in silico exon splicing prediction programs. We propose that casual synonymous substitutions may lead to a reduced splicing efficiency that can result in a variable proportion of exon loss. If this phenomenon happens in in-frame exons and to an extent tolerated by the cells it can have an important evolutionary effect since it may generate a substrate for natural selection of new splicing isoforms.

An RNA map predicting Nova-dependent splicing regulation

Nova proteins are a neuron-specific alternative splicing factors. We have combined bioinformatics, biochemistry and genetics to derive an RNA map describing the rules by which Nova proteins regulate alternative splicing. This map revealed that the position of Nova binding sites (YCAY clusters) in a pre-messenger RNA determines the outcome of splicing. The map correctly predicted Nova's effect to inhibit or enhance exon inclusion, which led us to examine the relationship between the map and Nova's mechanism of action. Nova binding to an exonic YCAY cluster changed the protein complexes assembled on pre-mRNA, blocking U1 snRNP (small nuclear ribonucleoprotein) binding and exon inclusion, whereas Nova binding to an intronic YCAY cluster enhanced spliceosome assembly and exon inclusion. Assays of splicing intermediates of Nova-regulated transcripts in mouse brain revealed that Nova preferentially regulates removal of introns harbouring (or closest to) YCAY clusters. These results define a genome-wide map relating the position of a cis-acting element to its regulation by an RNA binding protein, namely that Nova binding to YCAY clusters results in a local and asymmetric action to regulate spliceosome assembly and alternative splicing in neurons.

Splicing regulation in neurologic disease

The importance of alternative splicing in the regulation of diverse biological processes is reflected in the growing list of human diseases associated with known or suspected splicing defects. It is becoming evident that alternative splicing plays a particularly important role in neurologic disease, which is perhaps not surprising given the important role splicing plays in generating complexity and function in the brain. This review considers the evidence that defects in regulation of splicing may underlie many types of human neurologic diseases.

Identification of hnRNPs K, L and A2/B1 as candidate proteins involved in the nutritional regulation of mRNA splicing

Nutrient regulation of glucose-6-phosphate dehydrogenase (G6PD) expression occurs through changes in the rate of splicing of G6PD pre-mRNA. This posttranscriptional mechanism accounts for the 12- to 15-fold increase in G6PD expression in livers of mice that were starved and then refed a high-carbohydrate diet. Regulation of G6PD pre-mRNA splicing requires a cis-acting element in exon 12 of the pre-mRNA. Using RNA probes to exon 12 and nuclear extracts from livers of mice that were starved or refed, proteins of 60 kDa and 37 kDa were detected bound to nucleotides 65-79 of exon 12 and this binding was decreased by 50% with nuclear extracts from refed mice. The proteins were identified as hnRNPs K, L, and A2/B1 by LC-MS/MS. The decrease in binding of these proteins to exon 12 during refeeding was not accompanied by a decrease in the total amount of these proteins in total nuclear extract. HnRNPs K, L and A2/B1 have known roles in the regulation of mRNA splicing. The decrease in binding of these proteins during treatments that increase G6PD expression is consistent with a role for these proteins in the inhibition of G6PD mRNA splicing.

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.

Gene expression involves a series of steps in which specific short DNA or RNA segments are recognized by nucleic acid–binding proteins. One step that is particularly prominent and complex in humans and other vertebrates is the removal of introns and the ligation of exons in the process of pre-mRNA splicing. To better understand the sequences in exons that regulate this process, the authors have developed a method termed Neighborhood Inference that predicts the splicing regulatory activity of RNA segments based on the known splicing enhancer or silencer activity of other segments that have closely neighboring sequences. This method is applied to predict hundreds of new exonic splicing regulatory elements, as well as splicing-neutral sequences. A number of these predictions were validated experimentally, indicating that the number of exonic splicing regulatory sequences is larger than previously suspected. Neighborhood Inference scoring is also used to show that selection on exonic splicing enhancers (ESEs) frequently allows conversion of one ESE sequence to another over evolutionary time periods, suggesting that ESEs are, to at least some degree, interchangeable in constitutively spliced exons. The methods described may also find application in the study of other biomolecular processes that involve sequence-specific nucleic acid–binding proteins.

In silico and in vivo splicing analysis of MLH1 and MSH2 missense mutations shows exon- and tissue-specific effects

BACKGROUND

Abnormalities of pre-mRNA splicing are increasingly recognized as an important mechanism through which gene mutations cause disease. However, apart from the mutations in the donor and acceptor sites, the effects on splicing of other sequence variations are difficult to predict. Loosely defined exonic and intronic sequences have been shown to affect splicing efficiency by means of silencing and enhancement mechanisms. Thus, nucleotide substitutions in these sequences can induce aberrant splicing. Web-based resources have recently been developed to facilitate the identification of nucleotide changes that could alter splicing. However, computer predictions do not always correlate with in vivo splicing defects. The issue of unclassified variants in cancer predisposing genes is very important both for the correct ascertainment of cancer risk and for the understanding of the basic mechanisms of cancer gene function and regulation. Therefore we aimed to verify how predictions that can be drawn from in silico analysis correlate with results obtained in an in vivo splicing assay.

RESULTS

We analysed 99 hMLH1 and hMSH2 missense mutations with six different algorithms. Transfection of three different cell lines with 20 missense mutations, showed that a minority of them lead to defective splicing. Moreover, we observed that some exons and some mutations show cell-specific differences in the frequency of exon inclusion.

CONCLUSION

Our results suggest that the available algorithms, while potentially helpful in identifying splicing modulators especially when they are located in weakly defined exons, do not always correspond to an obvious modification of the splicing pattern. Thus caution must be used in assessing the pathogenicity of a missense or silent mutation with prediction programs. The variations observed in the splicing proficiency in three different cell lines suggest that nucleotide changes may dictate alternative splice site selection in a tissue-specific manner contributing to the widely observed phenotypic variability in inherited cancers.

Unweaving the Meanings of Messenger RNA Sequences

In addition to protein-coding information, mRNAs harbor regulatory sequences necessary for appropriate processing of their precursors. Goren et al. (2006) and Wang et al. (2006) explore the diversity of these signals and the rules by which they function.