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

A high proportion of DNA variants of BRCA1 and BRCA2 is associated with aberrant splicing in breast/ovarian cancer patients

PURPOSE

Most BRCA1/2 mutations are of unknown clinical relevance. An increasing amount of evidence indicates that there can be deleterious effects through the disruption of the splicing process. We have investigated the effect of aberrant splicing of BRCA1/2 on hereditary breast/ovarian cancer (HBOC).

EXPERIMENTAL DESIGN

DNA variants were analyzed with splicing prediction programs to select putative splicing mutations. Splicing assays of 57 genetic variants were done by lymphocyte reverse transcription-PCR and/or hybrid minigenes in HeLa and nontumor breast epithelial cells.

RESULTS

Twenty-four BRCA1/2 variants of Spanish HBOC patients were bioinformatically preselected. Functional assays showed that 12 variants induced anomalous splicing patterns, 6 of which accounted for 58.5% of BRCA1 families. To further evaluate the defective splicing of BRCA1/2, we analyzed 31 Breast Cancer Information Core Database (BIC) and two artificial variants that were generated by mutagenesis. Sixteen variants induced different degrees of aberrant splicing. Altogether, anomalous splicing was caused by 28 BRCA1/2 variants of all types, indicating that any DNA change can disrupt pre-mRNA processing. We show that a wide range of regulatory elements can be involved, including the canonical and cryptic splice sites, the polypyrimidine tract, and splicing enhancers/silencers. Twenty mutations were predicted to truncate the BRCA proteins and/or to delete essential domains, thus supporting a role in HBOC.

CONCLUSIONS

An important fraction of DNA variants of BRCA1/2 presents splicing aberrations that may represent a relevant disease-causing mechanism in HBOC. The identification of splicing disruptions by functional assays is a valuable tool to discriminate between benign polymorphisms and pathogenic mutations.

Expansion of the eukaryotic proteome by alternative splicing

The collection of components required to carry out the intricate processes involved in generating and maintaining a living, breathing and, sometimes, thinking organism is staggeringly complex. Where do all of the parts come from? Early estimates stated that about 100,000 genes would be required to make up a mammal; however, the actual number is less than one-quarter of that, barely four times the number of genes in budding yeast. It is now clear that the 'missing' information is in large part provided by alternative splicing, the process by which multiple different functional messenger RNAs, and therefore proteins, can be synthesized from a single gene.

Genomic features defining exonic variants that modulate splicing

A comparative analysis of SNPs and their exonic and intronic environments identifies the features predictive of splice affecting variants.

Background

Single point mutations at both synonymous and non-synonymous positions within exons can have severe effects on gene function through disruption of splicing. Predicting these mutations in silico purely from the genomic sequence is difficult due to an incomplete understanding of the multiple factors that may be responsible. In addition, little is known about which computational prediction approaches, such as those involving exonic splicing enhancers and exonic splicing silencers, are most informative.

Results

We assessed the features of single-nucleotide genomic variants verified to cause exon skipping and compared them to a large set of coding SNPs common in the human population, which are likely to have no effect on splicing. Our findings implicate a number of features important for their ability to discriminate splice-affecting variants, including the naturally occurring density of exonic splicing enhancers and exonic splicing silencers of the exon and intronic environment, extensive changes in the number of predicted exonic splicing enhancers and exonic splicing silencers, proximity to the splice junctions and evolutionary constraint of the region surrounding the variant. By extending this approach to additional datasets, we also identified relevant features of variants that cause increased exon inclusion and ectopic splice site activation.

Conclusions

We identified a number of features that have statistically significant representation among exonic variants that modulate splicing. These analyses highlight putative mechanisms responsible for splicing outcome and emphasize the role of features important for exon definition. We developed a web-tool, Skippy, to score coding variants for these relevant splice-modulating features.

Competing upstream 5' splice sites enhance the rate of proximal splicing

Alternative 5' splice site selection is one of the major pathways resulting in mRNA diversification. Regulation of this type of alternative splicing depends on the presence of regulatory elements that activate or repress the use of competing splice sites, usually leading to the preferential use of the proximal splice site. However, the mechanisms involved in proximal splice site selection and the thermodynamic advantage realized by proximal splice sites are not well understood. Here, we have carried out a systematic analysis of alternative 5' splice site usage using in vitro splicing assays. We show that observed rates of splicing correlate well with their U1 snRNA base pairing potential. Weak U1 snRNA interactions with the 5' splice site were significantly rescued by the proximity of the downstream exon, demonstrating that the intron definition mode of splice site recognition is highly efficient. In the context of competing splice sites, the proximity to the downstream 3' splice site was more influential in dictating splice site selection than the actual 5' splice site/U1 snRNA base pairing potential. Surprisingly, the kinetic analysis also demonstrated that an upstream competing 5' splice site enhances the rate of proximal splicing. These results reveal the discovery of a new splicing regulatory element, an upstream 5' splice site functioning as a splicing enhancer.

Overlapping splicing regulatory motifs--combinatorial effects on splicing

Regulation of splicing in eukaryotes occurs through the coordinated action of multiple splicing factors. Exons and introns contain numerous putative binding sites for splicing regulatory proteins. Regulation of splicing is presumably achieved by the combinatorial output of the binding of splicing factors to the corresponding binding sites. Although putative regulatory sites often overlap, no extensive study has examined whether overlapping regulatory sequences provide yet another dimension to splicing regulation. Here we analyzed experimentally-identified splicing regulatory sequences using a computational method based on the natural distribution of nucleotides and splicing regulatory sequences. We uncovered positive and negative interplay between overlapping regulatory sequences. Examination of these overlapping motifs revealed a unique spatial distribution, especially near splice donor sites of exons with weak splice donor sites. The positively selected overlapping splicing regulatory motifs were highly conserved among different species, implying functionality. Overall, these results suggest that overlap of two splicing regulatory binding sites is an evolutionary conserved widespread mechanism of splicing regulation. Finally, over-abundant motif overlaps were experimentally tested in a reporting minigene revealing that overlaps may facilitate a mode of splicing that did not occur in the presence of only one of the two regulatory sequences that comprise it.

Comparative in silico analyses and experimental validation of novel splice site and missense mutations in the genes MLH1 and MSH2

Hereditary non-polyposis colorectal cancer, an autosomal dominant predisposition to colorectal cancer and other malignancies, is caused by inactivating mutations of DNA mismatch repair genes, mainly MLH1 and MSH2. Missense mutations affect protein structure or function, but may also cause aberrant splicing, if located within splice sites (ss) or cis-acting sequences of splicing regulatory proteins, i.e., exonic splicing enhancers or exonic splicing silencers. Despite significant progress of ss scoring algorithms, the prediction for the impact of mutations on splicing is still unsatisfactory. For this study, we assessed ten ss and nine missense mutations outside ss in MLH1 and MSH2, including eleven newly identified mutations, and experimentally analyzed their effect at the RNA level. We additionally tested and compared the reliability of several web-based programs for the prediction of splicing outcome for these mutations.

COMIT: identification of noncoding motifs under selection in coding sequences

COMIT is presented; an algorithm for detecting functional non-coding motifs in coding regions, separating nucleotide and amino acid effects.

Coding nucleotide sequences contain myriad functions independent of their encoded protein sequences. We present the COMIT algorithm to detect functional noncoding motifs in coding regions using sequence conservation, explicitly separating nucleotide from amino acid effects. COMIT concurs with diverse experimental datasets, including splicing enhancers, silencers, replication motifs, and microRNA targets, and predicts many novel functional motifs. Intriguingly, COMIT scores are well-correlated to scores uncalibrated for amino acids, suggesting that nucleotide motifs often override peptide-level constraints.

The pivotal roles of TIA proteins in 5' splice-site selection of alu exons and across evolution

More than 5% of alternatively spliced internal exons in the human genome are derived from Alu elements in a process termed exonization. Alus are comprised of two homologous arms separated by an internal polypyrimidine tract (PPT). In most exonizations, splice sites are selected from within the same arm. We hypothesized that the internal PPT may prevent selection of a splice site further downstream. Here, we demonstrate that this PPT enhanced the selection of an upstream 5′ splice site (5′ss), even in the presence of a stronger 5′ss downstream. Deletion of this PPT shifted selection to the stronger downstream 5′ss. This enhancing effect depended on the strength of the downstream 5′ss, on the efficiency of base-pairing to U1 snRNA, and on the length of the PPT. This effect of the PPT was mediated by the binding of TIA proteins and was dependent on the distance between the PPT and the upstream 5′ss. A wide-scale evolutionary analysis of introns across 22 eukaryotes revealed an enrichment in PPTs within ∼20 nt downstream of the 5′ss. For most metazoans, the strength of the 5′ss inversely correlated with the presence of a downstream PPT, indicative of the functional role of the PPT. Finally, we found that the proteins that mediate this effect, TIA and U1C, and in particular their functional domains, are highly conserved across evolution. Overall, these findings expand our understanding of the role of TIA1/TIAR proteins in enhancing recognition of exons, in general, and Alu exons, in particular.

Human genes are composed of functional regions, termed exons, separated by non-functional regions, termed introns. Intronic sequences may gradually accumulate mutations and subsequently become recognized by the splicing machinery as exons, a process termed exonization. Alu elements are prone to undergo exonization: more than 5% of alternatively spliced internal exons in the human genome originate from Alu elements. A typical Alu element is ∼300 nucleotides long, consisting of two arms separated by a polypyrimdine tract (PPT). Interestingly, in most cases, exonization occurs almost exclusively within either the right arm or the left, not both. Here we found that the PPT between the two arms serves as a binding site for TIA proteins and prevents the exon selection process from expanding into downstream regions. To obtain a wider overview of TIA function, we performed a cross-evolutionary analysis within 22 eukaryotes of this protein and of U1C, a protein known to interact with it, and found that functional regions of both these proteins were highly conserved. These findings highlight the pivotal role of TIA proteins in 5′ splice-site selection of Alu exons and exon recognition in general.

Computational prediction of splicing regulatory elements shared by Tetrapoda organisms

Background

Auxiliary splicing sequences play an important role in ensuring accurate and efficient splicing by promoting or repressing recognition of authentic splice sites. These cis-acting motifs have been termed splicing enhancers and silencers and are located both in introns and exons. They co-evolved into an intricate splicing code together with additional functional constraints, such as tissue-specific and alternative splicing patterns. We used orthologous exons extracted from the University of California Santa Cruz multiple genome alignments of human and 22 Tetrapoda organisms to predict candidate enhancers and silencers that have reproducible and statistically significant bias towards annotated exonic boundaries.

Results

A total of 2,546 Tetrapoda enhancers and silencers were clustered into 15 putative core motifs based on their Markov properties. Most of these elements have been identified previously, but 118 putative silencers and 260 enhancers (~15%) were novel. Examination of previously published experimental data for the presence of predicted elements showed that their mutations in 21/23 (91.3%) cases altered the splicing pattern as expected. Predicted intronic motifs flanking 3' and 5' splice sites had higher evolutionary conservation than other sequences within intronic flanks and the intronic enhancers were markedly differed between 3' and 5' intronic flanks.

Conclusion

Difference in intronic enhancers supporting 5' and 3' splice sites suggests an independent splicing commitment for neighboring exons. Increased evolutionary conservation for ISEs/ISSs within intronic flanks and effect of modulation of predicted elements on splicing suggest functional significance of found elements in splicing regulation. Most of the elements identified were shown to have direct implications in human splicing and therefore could be useful for building computational splicing models in biomedical research.

Divergence of exonic splicing elements after gene duplication and the impact on gene structures

An analysis of human exonic splicing elements in duplicated genes reveals their important role in the generation of new gene structures.

Background

The origin of new genes and their contribution to functional novelty has been the subject of considerable interest. There has been much progress in understanding the mechanisms by which new genes originate. Here we examine a novel way that new gene structures could originate, namely through the evolution of new alternative splicing isoforms after gene duplication.

Results

We studied the divergence of exonic splicing enhancers and silencers after gene duplication and the contributions of such divergence to the generation of new splicing isoforms. We found that exonic splicing enhancers and exonic splicing silencers diverge especially fast shortly after gene duplication. About 10% and 5% of paralogous exons undergo significantly asymmetric evolution of exonic splicing enhancers and silencers, respectively. When compared to pre-duplication ancestors, we found that there is a significant overall loss of exonic splicing enhancers and the magnitude increases with duplication age. Detailed examination reveals net gains and losses of exonic splicing enhancers and silencers in different copies and paralog clusters after gene duplication. Furthermore, we found that exonic splicing enhancer and silencer changes are mainly caused by synonymous mutations, though nonsynonymous changes also contribute. Finally, we found that exonic splicing enhancer and silencer divergence results in exon splicing state transitions (from constitutive to alternative or vice versa), and that the proportion of paralogous exon pairs with different splicing states also increases over time, consistent with previous predictions.

Conclusions

Our results suggest that exonic splicing enhancer and silencer changes after gene duplication have important roles in alternative splicing divergence and that these changes contribute to the generation of new gene structures.

Characterization of conserved tandem donor sites and intronic motifs required for alternative splicing in corticosteroid receptor genes

Alternative splicing events from tandem donor sites result in mRNA variants coding for additional amino acids in the DNA binding domain of both the glucocorticoid (GR) and mineralocorticoid (MR) receptors. We now show that expression of both splice variants is extensively conserved in mammalian species, providing strong evidence for their functional significance. An exception to the conservation of the MR tandem splice site (an A at position +5 of the MR+12 donor site in the mouse) was predicted to decrease U1 small nuclear RNA binding. In accord with this prediction, we were unable to detect the MR+12 variant in this species. The one exception to the conservation of the GR tandem splice site, an A at position +3 of the platypus GRgamma donor site that was predicted to enhance binding of U1 snRNA, was unexpectedly associated with decreased expression of the variant from the endogenous gene as well as a minigene. An intronic pyrimidine motif present in both GR and MR genes was found to be critical for usage of the downstream donor site, and overexpression of TIA1/TIAL1 RNA binding proteins, which are known to bind such motifs, led to a marked increase in the proportion of GRgamma and MR+12. These results provide striking evidence for conservation of a complex splicing mechanism that involves processes other than stochastic spliceosome binding and identify a mechanism that would allow regulation of variant expression.

Alternative splicing attenuates transgenic expression directed by the apolipoprotein E promoter-enhancer based expression vector pLIV11

The plasmid vector pLIV11 is used commonly to achieve liver-specific expression of genes of interest in transgenic mice and rabbits. Expression is driven by the human apolipoprotein (apo)E 5' proximal promoter, which includes 5 kb of upstream sequence, exon 1, intron 1, and 5 bp of exon 2. A 3.8 kb 3' hepatic control region, derived from a region approximately 18 kb downstream of the apoE gene, enhances liver-specific expression. Here, we report that cDNA sequences inserted into the multiple cloning site (MCS) of pLIV11, which is positioned just downstream of truncated exon 2, can cause exon 2 skipping. Hence, splicing is displaced to downstream cryptic 3' splice acceptor sites causing deletion of cloned 5' untranslated mRNA sequences and, in some cases, deletion of the 5' end of an open reading frame. To prevent use of cryptic splice sites, the pLIV11 vector was modified with an engineered 3' splice acceptor site inserted immediately downstream of truncated apoE exon 2. Presence of this sequence fully shifted splicing of exon 1 from the native intron 1-exon 2 splice acceptor site to the engineered site. This finding confirmed that sequences inserted into the MCS of the vector pLIV11 can affect exon 2 recognition and provides a strategy to protect cloned sequences from alternative splicing and possible attenuation of transgenic expression.

Chromatin organization marks exon-intron structure

An increasing body of evidence indicates that transcription and splicing are coupled, and it is accepted that chromatin organization regulates transcription. Little is known about the cross-talk between chromatin structure and exon-intron architecture. By analysis of genome-wide nucleosome-positioning data sets from humans, flies and worms, we found that exons show increased nucleosome-occupancy levels with respect to introns, a finding that we link to differential GC content and nucleosome-disfavoring elements between exons and introns. Analysis of genome-wide chromatin immunoprecipitation data in humans and mice revealed four specific post-translational histone modifications enriched in exons. Our findings indicate that previously described enrichment of H3K36me3 modifications in exons reflects a more fundamental phenomenon, namely increased nucleosome occupancy along exons. Our results suggest an RNA polymerase II-mediated cross-talk between chromatin structure and exon-intron architecture, implying that exon selection may be modulated by chromatin structure.

Disease-causing mutations improving the branch site and polypyrimidine tract: pseudoexon activation of LINE-2 and antisense Alu lacking the poly(T)-tail

Cryptic exons or pseudoexons are typically activated by point mutations that create GT or AG dinucleotides of new 5' or 3' splice sites in introns, often in repetitive elements. Here we describe two cases of tetrahydrobiopterin deficiency caused by mutations improving the branch point sequence and polypyrimidine tracts of repeat-containing pseudoexons in the PTS gene. In the first case, we demonstrate a novel pathway of antisense Alu exonization, resulting from an intronic deletion that removed the poly(T)-tail of antisense AluSq. The deletion brought a favorable branch point sequence within proximity of the pseudoexon 3' splice site and removed an upstream AG dinucleotide required for the 3' splice site repression on normal alleles. New Alu exons can thus arise in the absence of poly(T)-tails that facilitated inclusion of most transposed elements in mRNAs by serving as polypyrimidine tracts, highlighting extraordinary flexibility of Alu repeats in shaping intron-exon structure. In the other case, a PTS pseudoexon was activated by an A>T substitution 9 nt upstream of its 3' splice site in a LINE-2 sequence, providing the first example of a disease-causing exonization of the most ancient interspersed repeat. These observations expand the spectrum of mutational mechanisms that introduce repetitive sequences in mature transcripts and illustrate the importance of intronic mutations in alternative splicing and phenotypic variability of hereditary disorders.

Genome-wide prediction of cis-acting RNA elements regulating tissue-specific pre-mRNA alternative splicing

Background

Human genes undergo various patterns of pre-mRNA splicing across different tissues. Such variation is primarily regulated by trans-acting factors that bind on exonic and intronic cis-acting RNA elements (CAEs). Here we report a computational method to mechanistically identify cis-acting RNA elements that contribute to the tissue-specific alternative splicing pattern. This method is an extension of our previous model, SplicingModeler, which predicts the significant CAEs that contribute to the splicing differences between two tissues. In this study, we introduce tissue-specific functional levels estimation step, which allows evaluating regulatory functions of predicted CAEs that are involved in more than two tissues.

Results

Using a publicly available Affymetrix Genechip^® Human Exon Array dataset, our method identifies 652 cis-acting RNA elements (CAEs) across 11 human tissues. About one third of predicted CAEs can be mapped to the known RBP (RNA binding protein) binding sites or match with other predicted exonic splicing regulator databases. Interestingly, the vast majority of predicted CAEs are in intronic regulatory regions. A noticeable exception is that many exonic elements are found to regulate the alternative splicing between cerebellum and testes. Most identified elements are found to contribute to the alternative splicing between two tissues, while some are important in multiple tissues. This suggests that genome-wide alternative splicing patterns are regulated by a combination of tissue-specific cis-acting elements and "general elements" whose functional activities are important but differ across multiple tissues.

Conclusion

In this study, we present a model-based computational approach to identify potential cis-acting RNA elements by considering the exon splicing variation as the combinatorial effects of multiple cis-acting regulators. This methodology provides a novel evaluation on the functional levels of cis-acting RNA elements by estimating their tissue-specific functions on various tissues.

Evolution of alternative splicing regulation: changes in predicted exonic splicing regulators are not associated with changes in alternative splicing levels in primates

Alternative splicing is tightly regulated in a spatio-temporal and quantitative manner. This regulation is achieved by a complex interplay between spliceosomal (trans) factors that bind to different sequence (cis) elements. cis-elements reside in both introns and exons and may either enhance or silence splicing. Differential combinations of cis-elements allows for a huge diversity of overall splicing signals, together comprising a complex ‘splicing code’. Many cis-elements have been identified, and their effects on exon inclusion levels demonstrated in reporter systems. However, the impact of interspecific differences in these elements on the evolution of alternative splicing levels has not yet been investigated at genomic level. Here we study the effect of interspecific differences in predicted exonic splicing regulators (ESRs) on exon inclusion levels in human and chimpanzee. For this purpose, we compiled and studied comprehensive datasets of predicted ESRs, identified by several computational and experimental approaches, as well as microarray data for changes in alternative splicing levels between human and chimpanzee. Surprisingly, we found no association between changes in predicted ESRs and changes in alternative splicing levels. This observation holds across different ESR exon positions, exon lengths, and 5′ splice site strengths. We suggest that this lack of association is mainly due to the great importance of context for ESR functionality: many ESR-like motifs in primates may have little or no effect on splicing, and thus interspecific changes at short-time scales may primarily occur in these effectively neutral ESRs. These results underscore the difficulties of using current computational ESR prediction algorithms to identify truly functionally important motifs, and provide a cautionary tale for studies of the effect of SNPs on splicing in human disease.

SROOGLE: webserver for integrative, user-friendly visualization of splicing signals

Exons are typically only 140 nt in length and are surrounded by intronic oceans that are thousands of nucleotides long. Four core splicing signals, aided by splicing-regulatory sequences (SRSs), direct the splicing machinery to the exon/intron junctions. Many different algorithms have been developed to identify and score the four splicing signals and thousands of putative SRSs have been identified, both computationally and experimentally. Here we describe SROOGLE, a webserver that makes splicing signal sequence and scoring data available to the biologist in an integrated, visual, easily interpretable, and user-friendly format. SROOGLE's input consists of the sequence of an exon and flanking introns. The graphic browser output displays the four core splicing signals with scores based on nine different algorithms and highlights sequences belonging to 13 different groups of SRSs. The interface also offers the ability to examine the effect of point mutations at any given position, as well a range of additional metrics and statistical measures regarding each potential signal. SROOGLE is available at http://sroogle.tau.ac.il, and may also be downloaded as a desktop version.

In vivo selection of kinase-responsive RNA elements controlling alternative splicing

Alternative pre-mRNA splicing is often controlled by cell signals, for example, those activating the cAMP-dependent protein kinase (PKA) or the Ca2+/calmodulin-dependent protein kinase IV (CaMKIV). We have shown that CaMKIV regulates alternative splicing through short CA repeats and hnRNP L. Here we use a splicing reporter that shows PKA/CaMKIV promotion of exon inclusion to select from exons containing random 13-nt sequences for RNA elements responsive to the kinases in cultured cells. This selection not only identified both PKA- and CaMKIV-responsive elements that are similar to the CaMKIV-responsive RNA element 1 (CaRRE1) or CA repeats, but also A-rich elements not previously known to respond to these kinases. Consistently, hnRNP L is identified as a factor binding the CA-rich elements. Analyses of the motifs in the highly responsive elements indicate that they are indeed critical for the kinase effect and are enriched in alternative exons. Interestingly, a CAAAAAA motif is sufficient for the PKA/CaMKIV-regulated splicing of the exon 16 of the CaMK kinase beta1 (CaMKK2) transcripts, implying a role of this motif in signaling cross-talk or feedback regulation between these kinases through alternative splicing. Therefore, these experiments identified a group of RNA elements responsive to PKA and CaMKIV from in vivo selection. This also provides an approach for selecting RNA elements similarly responsive to other cell signals controlling alternative splicing.

Genome-wide analysis of alternative pre-mRNA splicing and RNA-binding specificities of the Drosophila hnRNP A/B family members

Heterogeneous nuclear ribonucleoproteins (hnRNPs) have been traditionally seen as proteins packaging RNA nonspecifically into ribonucleoprotein particles (RNPs), but evidence suggests specific cellular functions on discrete target pre-mRNAs. Here we report genome-wide analysis of alternative splicing patterns regulated by four Drosophila homologs of the mammalian hnRNP A/B family (hrp36, hrp38, hrp40, and hrp48). Analysis of the global RNA-binding distributions of each protein revealed both small and extensively bound regions on target transcripts. A significant subset of RNAs were bound and regulated by more than one hnRNP protein, revealing a combinatorial network of interactions. In vitro RNA-binding site selection experiments (SELEX) identified distinct binding motif specificities for each protein, which were overrepresented in their respective regulated and bound transcripts. These results indicate that individual heterogeneous ribonucleoproteins have specific affinities for overlapping, but distinct, populations of target pre-mRNAs controlling their patterns of RNA processing.

Conserved alternative splicing and expression patterns of arthropod N-cadherin

Metazoan development requires complex mechanisms to generate cells with diverse function. Alternative splicing of pre-mRNA not only expands proteomic diversity but also provides a means to regulate tissue-specific molecular expression. The N-Cadherin gene in Drosophila contains three pairs of mutually-exclusive alternatively-spliced exons (MEs). However, no significant differences among the resulting protein isoforms have been successfully demonstrated in vivo. Furthermore, while the N-Cadherin gene products exhibit a complex spatiotemporal expression pattern within embryos, its underlying mechanisms and significance remain unknown. Here, we present results that suggest a critical role for alternative splicing in producing a crucial and reproducible complexity in the expression pattern of arthropod N-Cadherin. We demonstrate that the arthropod N-Cadherin gene has maintained the three sets of MEs for over 400 million years using in silico and in vivo approaches. Expression of isoforms derived from these MEs receives precise spatiotemporal control critical during development. Both Drosophila and Tribolium use ME-13a and ME-13b in “neural” and “mesodermal” splice variants, respectively. As proteins, either ME-13a- or ME-13b-containing isoform can cell-autonomously rescue the embryonic lethality caused by genetic loss of N-Cadherin. Ectopic muscle expression of either isoform beyond the time it normally ceases leads to paralysis and lethality. Together, our results offer an example of well-conserved alternative splicing increasing cellular diversity in metazoans.

Animal development requires complex mechanisms to generate many different types of cells. Alternative splicing is a process by which a single gene could produce several protein variants under particular circumstances. It is a useful means to generate a diversified set of proteins in different cell types. In this report, we showed that the alternative splicing of the arthropod N-Cadherin gene has been maintained for over 400 million years. The switch of expression patterns of two distinct variants is also well conserved in arthropods. As proteins, these two N-Cadherin splice variants have similar ability to rescue the embryonic lethality caused by genetic loss of N-Cadherin. However, when the expression of either isoforms was prolonged in muscles where the endogenous expression ceased beyond certain stages, it leads to larval lethality, suggesting the importance of precise spatiotemporal regulation of N-Cadherin splice-variant expression. This finding is particularly important because it offers an example of well-conserved alternative splicing increasing cellular diversity in animals.

Predicting potentially functional SNPs in drug-response genes

SNPs are known to contribute to variations in drug response and there are more than 14 million polymorphisms spanning the human genome. However, not all of these SNPs are functional. It would be impractical and costly to evaluate every individual SNP for functionality experimentally. Consequently, one of the major challenges for researchers has been to seek out functional SNPs from all the SNPs in the human genome. In silico or bioinformatic methods are economical, less labor intensive, yet powerful approaches to filter out potentially functional SNPs in drug-response genes for further study. This allows researchers to prioritize which SNPs to subsequently evaluate experimentally for drug-response studies, as well as potentially providing insights into possible mechanisms underlying how SNPs may affect drug-response genes.

Human Splicing Finder: an online bioinformatics tool to predict splicing signals

Thousands of mutations are identified yearly. Although many directly affect protein expression, an increasing proportion of mutations is now believed to influence mRNA splicing. They mostly affect existing splice sites, but synonymous, non-synonymous or nonsense mutations can also create or disrupt splice sites or auxiliary cis-splicing sequences. To facilitate the analysis of the different mutations, we designed Human Splicing Finder (HSF), a tool to predict the effects of mutations on splicing signals or to identify splicing motifs in any human sequence. It contains all available matrices for auxiliary sequence prediction as well as new ones for binding sites of the 9G8 and Tra2-β Serine-Arginine proteins and the hnRNP A1 ribonucleoprotein. We also developed new Position Weight Matrices to assess the strength of 5′ and 3′ splice sites and branch points. We evaluated HSF efficiency using a set of 83 intronic and 35 exonic mutations known to result in splicing defects. We showed that the mutation effect was correctly predicted in almost all cases. HSF could thus represent a valuable resource for research, diagnostic and therapeutic (e.g. therapeutic exon skipping) purposes as well as for global studies, such as the GEN2PHEN European Project or the Human Variome Project.

A computational approach for genome-wide mapping of splicing factor binding sites

A computational method is presented for genome-wide mapping of splicing factor binding sites that considers both the genomic environment and evolutionary conservation.

Alternative splicing is regulated by splicing factors that serve as positive or negative effectors, interacting with regulatory elements along exons and introns. Here we present a novel computational method for genome-wide mapping of splicing factor binding sites that considers both the genomic environment and the evolutionary conservation of the regulatory elements. The method was applied to study the regulation of different alternative splicing events, uncovering an interesting network of interactions among splicing factors.

Rational design of antisense oligomers to induce dystrophin exon skipping

Duchenne muscular dystrophy (DMD), one of the most severe neuromuscular disorders of childhood, is caused by the absence of a functional dystrophin. Antisense oligomer (AO) induced exon skipping is being investigated to restore functional dystrophin expression in models of muscular dystrophy and DMD patients. One of the major challenges will be in the development of clinically relevant oligomers and exon skipping strategies to address many different mutations. Various models, including cell-free extracts, cells transfected with artificial constructs, or mice with a human transgene, have been proposed as tools to facilitate oligomer design. Despite strong sequence homology between the human and mouse dystrophin genes, directing an oligomer to the same motifs in both species does not always induce comparable exon skipping. We report substantially different levels of exon skipping induced in normal and dystrophic human myogenic cell lines and propose that animal models or artificial assay systems useful in initial studies may be of limited relevance in designing the most efficient compounds to induce targeted skipping of human dystrophin exons for therapeutic outcomes.

A missense mutation in the APC tumor suppressor gene disrupts an ASF/SF2 splicing enhancer motif and causes pathogenic skipping of exon 14

A missense mutation at codon 640 in the APC gene was identified in a familial adenomatous polyposis (FAP) patient, however, its pathological consequence remained unclear. Here we found that this missense mutation interferes at the nucleotide level with an exonic splicing regulatory element and leads to aberrant splicing of the mutant APC transcript rather than exerting its effect through the observed amino acid change. Analysis of the patient RNA revealed complete skipping of exon 14 in transcripts from the mutant APC allele, leading to a frameshift and a premature stop codon. When cloned into a splicing reporter minigene and transfected into colorectal cell lines, the exon 14 point mutation c.1918C>G (pR640G) was found sufficient to cause the observed exon skipping. Bioinformatic analysis predicted the mutation to change SRp55, hnRNP A1 or ASF/SF2 splicing factor binding sites. Using RNA interference methodology these predictions were experimentally validated and revealed that only ASF/SF2 was required for exon 14 inclusion. These research data identify APC mutation c.1918C>G (pR640G) as pathogenic and indicate a mechanism involving disruption of an ASF/SF2 exonic splicing enhancer element. The results allow genetic diagnosis of a hereditary tumour predisposition but also illustrate the need to complement in silico prediction by splicing reporter assays.

Alu exonization events reveal features required for precise recognition of exons by the splicing machinery

Despite decades of research, the question of how the mRNA splicing machinery precisely identifies short exonic islands within the vast intronic oceans remains to a large extent obscure. In this study, we analyzed Alu exonization events, aiming to understand the requirements for correct selection of exons. Comparison of exonizing Alus to their non-exonizing counterparts is informative because Alus in these two groups have retained high sequence similarity but are perceived differently by the splicing machinery. We identified and characterized numerous features used by the splicing machinery to discriminate between Alu exons and their non-exonizing counterparts. Of these, the most novel is secondary structure: Alu exons in general and their 5′ splice sites (5′ss) in particular are characterized by decreased stability of local secondary structures with respect to their non-exonizing counterparts. We detected numerous further differences between Alu exons and their non-exonizing counterparts, among others in terms of exon–intron architecture and strength of splicing signals, enhancers, and silencers. Support vector machine analysis revealed that these features allow a high level of discrimination (AUC = 0.91) between exonizing and non-exonizing Alus. Moreover, the computationally derived probabilities of exonization significantly correlated with the biological inclusion level of the Alu exons, and the model could also be extended to general datasets of constitutive and alternative exons. This indicates that the features detected and explored in this study provide the basis not only for precise exon selection but also for the fine-tuned regulation thereof, manifested in cases of alternative splicing.

A typical human gene consists of 9 exons around 150 nucleotides in length, separated by introns that are ∼3,000 nucleotides long. The challenge of the splicing machinery is to precisely identify and ligate the exons, while removing the introns. We aimed to understand how the splicing machinery meets this momentous challenge, based on Alu exonization events. Alus are transposable elements, of which approximately one million copies exist in the human genome, a large portion of which within introns. Throughout evolution, some intronic Alus accumulated mutations and became recognized by the splicing machinery as exons, a process termed exonization. Such Alus remain highly similar to their non-exonizing counterparts but are perceived as different by the splicing machinery. By comparing exonizing Alus to their non-exonizing counterparts, we were able to identify numerous features in which they differ and which presumably lead to the recognition only of the former by the splicing machinery. Our findings reveal insights regarding the role of local RNA secondary structures, exon–intron architecture constraints, and splicing regulatory signals. We integrated these features in a computational model, which was able to successfully mimic the function of the splicing machinery and discriminate between true Alu exons and their intronic counterparts, highlighting the functional importance of these features.

Splicing of designer exons reveals unexpected complexity in pre-mRNA splicing

Pre-messengerRNA (mRNA) splicing requires the accurate recognition of splice sites by the cellular RNA processing machinery. In addition to sequences that comprise the branchpoint and the 3' and 5' splice sites, the cellular splicing machinery relies on additional information in the form of exonic and intronic splicing enhancer and silencer sequences. The high abundance of these motifs makes it difficult to investigate their effects using standard genetic perturbations, since their disruption often leads to the formation of yet new elements. To lessen this problem, we have designed synthetic exons comprised of multiple copies of a single prototypical exonic enhancer and a single prototypical exonic silencer sequence separated by neutral spacer sequences. The spacer sequences buffer the exon against the formation of new elements as the number and order of the original elements are varied. Over 100 such designer exons were constructed by random ligation of enhancer, silencer, and neutral elements. Each exon was positioned as the central exon in a 3-exon minigene and tested for exon inclusion after transient transfection. The level of inclusion of the test exons was seen to be dependent on the provision of enhancers and could be decreased by the provision of silencers. In general, there was a good quantitative correlation between the proportion of enhancers and splicing. However, widely varying inclusion levels could be produced by different permutations of the enhancer and silencer elements, indicating that even in this simplified system splicing decisions rest on complex interplays of yet to be determined parameters.

Exon creation and establishment in human genes

BACKGROUND

A large proportion of species-specific exons are alternatively spliced. In primates, Alu elements play a crucial role in the process of exon creation but many new exons have appeared through other mechanisms. Despite many recent studies, it is still unclear which are the splicing regulatory requirements for de novo exonization and how splicing regulation changes throughout an exon's lifespan.

RESULTS

Using comparative genomics, we have defined sets of exons with different evolutionary ages. Younger exons have weaker splice-sites and lower absolute values for the relative abundance of putative splicing regulators between exonic and adjacent intronic regions, indicating a less consolidated splicing regulation. This relative abundance is shown to increase with exon age, leading to higher exon inclusion. We show that this local difference in the density of regulators might be of biological significance, as it outperforms other measures in real exon versus pseudo-exon classification. We apply this new measure to the specific case of the exonization of anti-sense Alu elements and show that they are characterized by a general lack of exonic splicing silencers.

CONCLUSIONS

Our results suggest that specific sequence environments are required for exonization and that these can change with time. We propose a model of exon creation and establishment in human genes, in which splicing decisions depend on the relative local abundance of regulatory motifs. Using this model, we provide further explanation as to why Alu elements serve as a major substrate for exon creation in primates. Finally, we discuss the benefits of integrating such information in gene prediction.

Ab initio prediction of mutation-induced cryptic splice-site activation and exon skipping

Mutations that affect splicing of precursor messenger RNAs play a major role in the development of hereditary diseases. Most splicing mutations have been found to eliminate GT or AG dinucleotides that define the 5' and 3' ends of introns, leading to exon skipping or cryptic splice-site activation. Although accurate description of the mis-spliced transcripts is critical for predicting phenotypic consequences of these alterations, their exact nature in affected individuals cannot often be determined experimentally. Using a comprehensive collection of exons that sustained cryptic splice-site activation or were skipped as a result of splice-site mutations, we have developed a multivariate logistic discrimination procedure that distinguishes the two aberrant splicing outcomes from DNA sequences. The new algorithm was validated using an independent sample of exons and implemented as a free online utility termed CRYP-SKIP (http://www.dbass.org.uk/cryp-skip/). The web application takes up one or more mutated alleles, each consisting of one exon and flanking intronic sequences, and provides a list of important predictor variables and their values, the overall probability of activating cryptic splice vs exon skipping, and the location and intrinsic strength of predicted cryptic splice sites in the input sequence. These results will facilitate phenotypic prediction of splicing mutations and provide further insights into splicing enhancer and silencer elements and their relative importance for splice-site selection in vivo.

An RNA code for the FOX2 splicing regulator revealed by mapping RNA-protein interactions in stem cells

The elucidation of a code for regulated splicing has been a long-standing goal in understanding the control of post-transcriptional gene expression events that are crucial for cell survival, differentiation and development. We decoded functional RNA elements in vivo by constructing an RNA map for the cell type-specific splicing regulator FOX2 (also known as RBM9) via cross-linking immunoprecipitation coupled with high-throughput sequencing (CLIP-seq) in human embryonic stem cells. The map identified a large cohort of specific FOX2 targets, many of which are themselves splicing regulators, and comparison between the FOX2 binding profile and validated splicing events revealed a general rule for FOX2-regulated exon inclusion or skipping in a position-dependent manner. These findings suggest that FOX2 functions as a critical regulator of a splicing network, and we further show that FOX2 is important for the survival of human embryonic stem cells.

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 alpha-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.

A regulatory module embedded in the coding region of Hoxa2 controls expression in rhombomere 2

Here, we define a gene regulatory network for Hoxa2, responsible for temporal and spatial expression in hindbrain development. Hoxa2 plays an important role in regulating the regional identity of rhombomere 2 (r2) and is the only Hox gene expressed in this segment. In this study, we found that a Hoxa2 cis-regulatory module consists of five elements that direct expression in r2 of the developing hindbrain. Surprisingly, the module is imbedded in the second coding exon of Hoxa2 and therefore may be constrained by both protein coding and gene regulatory requirements. This highly conserved enhancer consists of two consensus Sox binding sites and several additional elements that act in concert to direct strong r2 specific expression. Our findings provide important insight into the regulation of segmental identity in the anterior hindbrain. Furthermore, they have broader implications in designing arrays and interpreting data from global analyses of gene regulation because regulatory input from coding regions needs to be considered.

Prediction of alternatively skipped exons and splicing enhancers from exon junction arrays

BACKGROUND

Alternative splicing of exons in a pre-mRNA transcript is an important mechanism which contributes to protein diversity in human. Arrays for detecting alternative splicing are available using several different probe designs, including those based on exon-junctions. In this work, we introduce a new method for predicting alternatively skipped exons from exon-junction arrays. Predictions based on our method are compared against controls and their sequences are analyzed to identify motifs important for regulating alternative splicing.

RESULTS

Our comparison of several alternative methods shows that an exon-skipping score based on neighboring junctions best discriminates between positive and negative controls. Sequence analysis of our predicted exons confirms the presence of known splicing regulatory sequences. In addition, we also derive a set of development-related alternatively spliced genes based on fetal versus adult tissue comparisons and find that our predictions are consistent with their functional annotations. Ab initio motif finding algorithms are applied to identify several motifs that may be relevant for splicing during development.

CONCLUSION

This work describes a new method for analyzing exon-junction arrays, identifies sequence motifs that are specific for alternative and constitutive splicing and suggests a role for several known splicing factors and their motifs in developmental regulation.

Improved identification of conserved cassette exons using Bayesian networks

Background

Alternative splicing is a major contributor to the diversity of eukaryotic transcriptomes and proteomes. Currently, large scale detection of alternative splicing using expressed sequence tags (ESTs) or microarrays does not capture all alternative splicing events. Moreover, for many species genomic data is being produced at a far greater rate than corresponding transcript data, hence in silico methods of predicting alternative splicing have to be improved.

Results

Here, we show that the use of Bayesian networks (BNs) allows accurate prediction of evolutionary conserved exon skipping events. At a stringent false positive rate of 0.5%, our BN achieves an improved true positive rate of 61%, compared to a previously reported 50% on the same dataset using support vector machines (SVMs). Incorporating several novel discriminative features such as intronic splicing regulatory elements leads to the improvement. Features related to mRNA secondary structure increase the prediction performance, corroborating previous findings that secondary structures are important for exon recognition. Random labelling tests rule out overfitting. Cross-validation on another dataset confirms the increased performance. When using the same dataset and the same set of features, the BN matches the performance of an SVM in earlier literature. Remarkably, we could show that about half of the exons which are labelled constitutive but receive a high probability of being alternative by the BN, are in fact alternative exons according to the latest EST data. Finally, we predict exon skipping without using conservation-based features, and achieve a true positive rate of 29% at a false positive rate of 0.5%.

Conclusion

BNs can be used to achieve accurate identification of alternative exons and provide clues about possible dependencies between relevant features. The near-identical performance of the BN and SVM when using the same features shows that good classification depends more on features than on the choice of classifier. Conservation based features continue to be the most informative, and hence distinguishing alternative exons from constitutive ones without using conservation based features remains a challenging problem.

Expression of 24,426 human alternative splicing events and predicted cis regulation in 48 tissues and cell lines

Alternative pre–messenger RNA splicing impacts development, physiology, and disease, but its regulation in humans is not well understood, partially due to the limited scale to which the expression of specific splicing events has been measured. We generated the first genome-scale expression compendium of human alternative splicing events using custom whole-transcript microarrays monitoring expression of 24,426 alternative splicing events in 48 diverse human samples. Over 11,700 genes and 9,500 splicing events were differentially expressed, providing a rich resource for studying splicing regulation. An unbiased, systematic screen of 21,760 4-mer to 7-mer words for cis-regulatory motifs identified 143 RNA 'words' enriched near regulated cassette exons, including six clusters of motifs represented by UCUCU, UGCAUG, UGCU, UGUGU, UUUU, and AGGG, which map to trans-acting regulators PTB, Fox, Muscleblind, CELF/CUG-BP, TIA-1, and hnRNP F/H, respectively. Each cluster showed a distinct pattern of genomic location and tissue specificity. For example, UCUCU occurs 110 to 35 nucleotides preceding cassette exons upregulated in brain and striated muscle but depleted in other tissues. UCUCU and UGCAUG appear to have similar function but independent action, occurring 5' and 3', respectively, of 33% of the cassette exons upregulated in skeletal muscle but co-occurring for only 2%.

Diverse splicing patterns of exonized Alu elements in human tissues

Exonization of Alu elements is a major mechanism for birth of new exons in primate genomes. Prior analyses of expressed sequence tags show that almost all Alu-derived exons are alternatively spliced, and the vast majority of these exons have low transcript inclusion levels. In this work, we provide genomic and experimental evidence for diverse splicing patterns of exonized Alu elements in human tissues. Using Exon array data of 330 Alu-derived exons in 11 human tissues and detailed RT-PCR analyses of 38 exons, we show that some Alu-derived exons are constitutively spliced in a broad range of human tissues, and some display strong tissue-specific switch in their transcript inclusion levels. Most of such exons are derived from ancient Alu elements in the genome. In SEPN1, mutations of which are linked to a form of congenital muscular dystrophy, the muscle-specific inclusion of an Alu-derived exon may be important for regulating SEPN1 activity in muscle. Realtime qPCR analysis of this SEPN1 exon in macaque and chimpanzee tissues indicates human-specific increase in its transcript inclusion level and muscle specificity after the divergence of humans and chimpanzees. Our results imply that some Alu exonization events may have acquired adaptive benefits during the evolution of primate transcriptomes.

New exons have been created and added to existing functional genes during eukaryotic genome evolution. Alu elements, a class of primate-specific retrotransposons, are a major source of new exons in primates. However, recent analyses of expressed sequence tags suggest that the vast majority of Alu-derived exons are low-abundance splice forms and represent non-functional evolutionary intermediates. In order to elucidate the evolutionary impact of Alu-derived exons, we investigated the splicing of 330 Alu-derived exons in 11 human tissues using data from high-density exon arrays with multiple oligonucleotide probes for every exon in the human genome. Our exon array analysis and further RT-PCR experiments reveal surprisingly diverse splicing patterns of these exons. Some Alu-derived exons are constitutively spliced, and some are strongly tissue-specific. In SEPN1, a gene implicated in a form of congenital muscular dystrophy, our data suggest that the muscle-specific inclusion of an Alu-derived exon results from a human-specific splicing change after the divergence of humans and chimpanzees. Our study provides novel insight into the evolutionary significance of Alu exonization events. A subset of Alu-derived exons, especially those derived from more ancient Alu elements in the genome, may have contributed to functional novelties during primate evolution.

Intronic Alus influence alternative splicing

Examination of the human transcriptome reveals higher levels of RNA editing than in any other organism tested to date. This is indicative of extensive double-stranded RNA (dsRNA) formation within the human transcriptome. Most of the editing sites are located in the primate-specific retrotransposed element called Alu. A large fraction of Alus are found in intronic sequences, implying extensive Alu-Alu dsRNA formation in mRNA precursors. Yet, the effect of these intronic Alus on splicing of the flanking exons is largely unknown. Here, we show that more Alus flank alternatively spliced exons than constitutively spliced ones; this is especially notable for those exons that have changed their mode of splicing from constitutive to alternative during human evolution. This implies that Alu insertions may change the mode of splicing of the flanking exons. Indeed, we demonstrate experimentally that two Alu elements that were inserted into an intron in opposite orientation undergo base-pairing, as evident by RNA editing, and affect the splicing patterns of a downstream exon, shifting it from constitutive to alternative. Our results indicate the importance of intronic Alus in influencing the splicing of flanking exons, further emphasizing the role of Alus in shaping of the human transcriptome.

Phosphorylation switches the general splicing repressor SRp38 to a sequence-specific activator

SRp38 is an atypical SR protein that functions as a general splicing repressor when dephosphorylated. We now show that phosphorylated SRp38 functions as a sequence-specific splicing activator. Unlike characterized splicing activators, SRp38 functions in the absence of other SR proteins but requires a cofactor for activity. SRp38 was able to induce formation of splicing complex A in the absence of the cofactor, but this factor was necessary for progression to complexes B and C. Mechanistically, SRp38 strengthens the ability of the U1 and U2 small nuclear ribonucleoproteins to stably recognize the pre-mRNA. Extending these findings, analysis of alternative splicing of pre-mRNA encoding the glutamate receptor B revealed that SRp38 alters its splicing pattern in a sequence-specific manner. Together, our data demonstrate that SRp38, in addition to its role as a splicing repressor, can function as an unusual sequence-specific splicing activator.

Genome-wide analyses of alternative splicing in plants: opportunities and challenges

Alternative splicing (AS) creates multiple mRNA transcripts from a single gene. While AS is known to contribute to gene regulation and proteome diversity in animals, the study of its importance in plants is in its early stages. However, recently available plant genome and transcript sequence data sets are enabling a global analysis of AS in many plant species. Results of genome analysis have revealed differences between animals and plants in the frequency of alternative splicing. The proportion of plant genes that have one or more alternative transcript isoforms is approximately 20%, indicating that AS in plants is not rare, although this rate is approximately one-third of that observed in human. The majority of plant AS events have not been functionally characterized, but evidence suggests that AS participates in important plant functions, including stress response, and may impact domestication and trait selection. The increasing availability of plant genome sequence data will enable larger comparative analyses that will identify functionally important plant AS events based on their evolutionary conservation, determine the influence of genome duplication on the evolution of AS, and discover plant-specific cis-elements that regulate AS. This review summarizes recent analyses of AS in plants, discusses the importance of further analysis, and suggests directions for future efforts.

Quantitative analysis of clinically relevant mutations occurring in lymphoid cells harboring gamma-retrovirus-encoded hsvtk suicide genes

The in vivo regulation of T lymphocyte activity by the activation of a suicide mechanism is an essential paradigm for the safety of adoptive cell therapies. In light of reports showing that gamma-retroviral vector-encoded herpes simplex virus thymidine kinase (hsvtk) undergoes recombination, we undertook a thorough investigation of the genomic stability of SFG-based vectors using two variants of the wild-type hsvtk gene. In a large panel of independent clones, we demonstrate that both hsvtk genes undergo recombination with molecular signatures indicative of template switching in GC-rich regions displaying homology at the deletion junctions or RNA splicing. In the absence of ganciclovir selection, the frequency of recombination is 3% per retroviral replication cycle. Our results underscore the importance of the five nucleotide difference between the two hsvtk genes that account for the presence of recombinogenic hot spots in one variant and not the other, indicating that the probability of RNA splicing is influenced by minute nucleotide changes in sequences adjacent to the splice donor and acceptor sites. Furthermore, our mutational analysis in an unbiased panel of human lymphoid cells (that is, without immune or ganciclovir-mediated selective pressure) provides a robust in vitro assay to predict and quantify clinically relevant mutations in hsvtk suicide genes, which can be applied to studying and improving the stability of any transgene expressed in gamma-retroviral or lentiviral vectors.

Pre-mRNA secondary structures influence exon recognition

The secondary structure of a pre-mRNA influences a number of processing steps including alternative splicing. Since most splicing regulatory proteins bind to single-stranded RNA, the sequestration of RNA into double strands could prevent their binding. Here, we analyzed the secondary structure context of experimentally determined splicing enhancer and silencer motifs in their natural pre-mRNA context. We found that these splicing motifs are significantly more single-stranded than controls. These findings were validated by transfection experiments, where the effect of enhancer or silencer motifs on exon skipping was much more pronounced in single-stranded conformation. We also found that the structural context of predicted splicing motifs is under selection, suggesting a general importance of secondary structures on splicing and adding another level of evolutionary constraints on pre-mRNAs. Our results explain the action of mutations that affect splicing and indicate that the structural context of splicing motifs is part of the mRNA splicing code.

Almost all human protein-coding genes contain several exons and introns. Prior to translation, introns have to be removed and exons have to be joined, which happens in a processing step called splicing that generates the mature mRNA. For most genes, certain exons can be either included or excluded from the mature mRNA. It is currently not fully understood which signals are needed to accurately recognize the boundaries of exons in the intron-containing primary transcript. As in transcriptional regulation, enhancer and silencer sequence motifs are crucial for the correct recognition of exons. Splicing regulatory proteins identify these motifs in a sequence-specific manner. In general, these proteins bind to single-stranded RNA. Here, we analyzed local secondary structures of primary transcripts and found that known splicing motifs are preferentially located in a single-stranded context. Experimental tests demonstrated that motifs in single-stranded contexts have a stronger effect on splice site selection than those located in double-stranded regions. These results help to understand the action of human mutations that change the splicing pattern and indicate that local pre-mRNA secondary structures influence exon recognition.

The "alternative" choice of constitutive exons throughout evolution

Alternative cassette exons are known to originate from two processes—exonization of intronic sequences and exon shuffling. Herein, we suggest an additional mechanism by which constitutively spliced exons become alternative cassette exons during evolution. We compiled a dataset of orthologous exons from human and mouse that are constitutively spliced in one species but alternatively spliced in the other. Examination of these exons suggests that the common ancestors were constitutively spliced. We show that relaxation of the 5′ splice site during evolution is one of the molecular mechanisms by which exons shift from constitutive to alternative splicing. This shift is associated with the fixation of exonic splicing regulatory sequences (ESRs) that are essential for exon definition and control the inclusion level only after the transition to alternative splicing. The effect of each ESR on splicing and the combinatorial effects between two ESRs are conserved from fish to human. Our results uncover an evolutionary pathway that increases transcriptome diversity by shifting exons from constitutive to alternative splicing.

Alternative splicing is believed to play a major role in the creation of transcriptomic diversification leading to higher order of organismal complexity, especially in mammals. As much as 80% of human genes generate more than one type of mRNA by alternative splicing. Thus, alternative splicing can bridge the low number of protein coding genes (∼24,500) and the total number of proteins generated in the human proteome (∼90,000). The correlation between the higher order of phenotypic diversity and alternative splicing was recently demonstrated and thus the origin of alternative splicing is of great interest. There are currently two models regarding the origin of alternatively spliced exons—exonization of intronic sequences and exon shuffling. According to these two mechanisms, a protein-coding gene was first established and only then a new alternative exon appeared within it or was added to the gene. Our current study provides evidences for a new mechanism indicating that during evolution constitutively spliced exons became alternatively spliced. Large-scale bioinformatic analyses reveal the magnitude of this process and experimental validation systems provide insights into its mechanisms.

Splicing regulation: from a parts list of regulatory elements to an integrated splicing code

Alternative splicing of pre-mRNAs is a major contributor to both proteomic diversity and control of gene expression levels. Splicing is tightly regulated in different tissues and developmental stages, and its disruption can lead to a wide range of human diseases. An important long-term goal in the splicing field is to determine a set of rules or "code" for splicing that will enable prediction of the splicing pattern of any primary transcript from its sequence. Outside of the core splice site motifs, the bulk of the information required for splicing is thought to be contained in exonic and intronic cis-regulatory elements that function by recruitment of sequence-specific RNA-binding protein factors that either activate or repress the use of adjacent splice sites. Here, we summarize the current state of knowledge of splicing cis-regulatory elements and their context-dependent effects on splicing, emphasizing recent global/genome-wide studies and open questions.

Functional importance of different patterns of correlation between adjacent cassette exons in human and mouse

Background

Alternative splicing expands transcriptome diversity and plays an important role in regulation of gene expression. Previous studies focus on the regulation of a single cassette exon, but recent experiments indicate that multiple cassette exons within a gene may interact with each other. This interaction can increase the potential to generate various transcripts and adds an extra layer of complexity to gene regulation. Several cases of exon interaction have been discovered. However, the extent to which the cassette exons coordinate with each other remains unknown.

Results

Based on EST data, we employed a metric of correlation coefficients to describe the interaction between two adjacent cassette exons and then categorized these exon pairs into three different groups by their interaction (correlation) patterns. Sequence analysis demonstrates that strongly-correlated groups are more conserved and contain a higher proportion of pairs with reading frame preservation in a combinatorial manner. Multiple genome comparison further indicates that different groups of correlated pairs have different evolutionary courses: (1) The vast majority of positively-correlated pairs are old, (2) most of the weakly-correlated pairs are relatively young, and (3) negatively-correlated pairs are a mixture of old and young events.

Conclusion

We performed a large-scale analysis of interactions between adjacent cassette exons. Compared with weakly-correlated pairs, the strongly-correlated pairs, including both the positively and negatively correlated ones, show more evidence that they are under delicate splicing control and tend to be functionally important. Additionally, the positively-correlated pairs bear strong resemblance to constitutive exons, which suggests that they may evolve from ancient constitutive exons, while negatively and weakly correlated pairs are more likely to contain newly emerging exons.

Searching for splicing motifs

Intron removal during pre-mRNA splicing in higher eukaryotes requires the accurate identification of the two splice sites at the ends of the exons, or exon definition. The sequences constituting the splice sites provide insufficient information to distinguish true splice sites from the greater number of false splice sites that populate transcripts. Additional information used for exon recognition resides in a large number of positively or negatively acting elements that lie both within exons and in the adjacent introns. The identification of such sequence motifs has progressed rapidly in recent years, such that extensive lists are now available for exonic splicing enhancers and exonic splicing silencers. These motifs have been identified both by empirical experiments and by computational predictions, the validity of the latter being confirmed by experimental verification. Molecular searches have been carried out either by the selection of sequences that bind to splicing factors, or enhance or silence splicing in vitro or in vivo. Computational methods have focused on sequences of 6 or 8 nucleotides that are over- or under-represented in exons, compared to introns or transcripts that do not undergo splicing. These various methods have sought to provide global definitions of motifs, yet the motifs are distinctive to the method used for identification and display little overlap. Astonishingly, at least three-quarters of a typical mRNA would be comprised of these motifs. A present challenge lies in understanding how the cell integrates this surfeit of information to generate what is usually a binary splicing decision.

A G-tract element in apoptotic agents-induced alternative splicing

Alternative splicing of a single pre-mRNA transcript can produce protein isoforms that promote either cell growth or death. Here we show that Ro-31-8220 (Ro), an apoptotic agent that inhibits protein kinase C and activates the c-Jun N terminal kinase, decreased the proportion of the cell growth-promoting Bcl-xL splice variant. Targeted mutagenesis analyses narrowed down a critical sequence to a 16-nt G-tract element (Gt16). Transferring this element to a heterologous gene conferred Ro response on an otherwise constitutive exon. The Ro effect was reduced by okadaic acid, an inhibitor of protein phosphatases PP1 and PP2A, in a concentration-dependent manner. Search in the human genome followed by RT–PCR identified a group of genes that contain similar exonic G-tract elements and are responsive to Ro. Moreover, the Gt16 element also mediates the regulation of alternative splicing by other cell apoptosis-inducers particularly retinoic acid. Therefore, the G-tract element likely plays a role in the apoptotic agents-induced alternative splicing of a group of genes. The functions of these genes imply that this regulation will have impact on cell growth/death.

RNA landscape of evolution for optimal exon and intron discrimination

Accurate pre-mRNA splicing requires primary splicing signals, including the splice sites, a polypyrimidine tract, and a branch site, other splicing-regulatory elements (SREs). The SREs include exonic splicing enhancers (ESEs), exonic splicing silencers (ESSs), intronic splicing enhancers (ISEs), and intronic splicing silencers (ISSs), which are typically located near the splice sites. However, it is unclear to what extent splicing-driven selective pressure constrains exonic and intronic sequences, especially those distant from the splice sites. Here, we studied the distribution of SREs in human genes in terms of DNA strand-asymmetry patterns. Under a neutral evolution model, each mononucleotide or oligonucleotide should have a symmetric (Chargaff's second parity rule), or weakly asymmetric yet uniform, distribution throughout a pre-mRNA transcript. However, we found that large sets of unbiased, experimentally determined SREs show a distinct strand-asymmetry pattern that is inconsistent with the neutral evolution model, and reflects their functional roles in splicing. ESEs are selected in exons and depleted in introns and vice versa for ESSs. Surprisingly, this trend extends into deep intronic sequences, accounting for one third of the genome. Selection is detectable even at the mononucleotide level, so that the asymmetric base compositions of exons and introns are predictive of ESEs and ESSs. We developed a method that effectively predicts SREs based on strand asymmetry, expanding the current catalog of SREs. Our results suggest that human genes have been optimized for exon and intron discrimination through an RNA landscape shaped during evolution.

Alternative splicing: current perspectives

Alternative splicing is a well-characterized mechanism by which multiple transcripts are generated from a single mRNA precursor. By allowing production of several protein isoforms from one pre-mRNA, alternative splicing contributes to proteomic diversity. But what do we know about the origin of this mechanism? Do the same evolutionary forces apply to alternatively and constitutively splice exons? Do similar forces act on all types of alternative splicing? Are the products generated by alternative splicing functional? Why is "improper" recognition of exons and introns allowed by the splicing machinery? In this review, we summarize the current knowledge regarding these issues from an evolutionary perspective.