The architecture of pre-mRNAs affects mechanisms of splice-site pairing

Patterns of exon-intron architecture variation of genes in eukaryotic genomes

Background

The origin and importance of exon-intron architecture comprises one of the remaining mysteries of gene evolution. Several studies have investigated the variations of intron length, GC content, ordinal position in a gene and divergence. However, there is little study about the structural variation of exons and introns.

Results

We investigated the length, GC content, ordinal position and divergence in both exons and introns of 13 eukaryotic genomes, representing plant and animal. Our analyses revealed that three basic patterns of exon-intron variation were present in nearly all analyzed genomes (P < 0.001 in most cases): an ordinal reduction of length and divergence in both exon and intron, a co-variation between exon and its flanking introns in their length, GC content and divergence, and a decrease of average exon (or intron) length, GC content and divergence as the total exon numbers of a gene increased. In addition, we observed that the shorter introns had either low or high GC content, and the GC content of long introns was intermediate.

Conclusion

Although the factors contributing to these patterns have not been identified, our results provide three important clues: common factor(s) exist and may shape both exons and introns; the ordinal reduction patterns may reflect a time-orderly evolution; and the larger first and last exons may be splicing-required. These clues provide a framework for elucidating mechanisms involved in the organization of eukaryotic genomes and particularly in building exon-intron structures.

Spliceosome assembly pathways for different types of alternative splicing converge during commitment to splice site pairing in the A complex

Differential splice site pairing establishes alternative splicing patterns resulting in the generation of multiple mRNA isoforms. This process is carried out by the spliceosome, which is activated by a series of sequential structural rearrangements of its five core snRNPs. To determine when splice sites become functionally paired, we carried out a series of kinetic trap experiments using pre-mRNAs that undergo alternative 5' splice site selection or alternative exon inclusion. We show that commitment to splice site pairing in both cases occurs in the A complex, which is characterized by the ATP-dependent association of the U2 snRNP with the branch point. Interestingly, the timing of splice site pairing is independent of the intron or exon definition modes of splice site recognition. Using the ATP analog ATPgammaS, we showed that ATP hydrolysis is required for splice site pairing independent from U2 snRNP binding to the pre-mRNA. These results identify the A complex as the spliceosomal assembly step dedicated to splice site pairing and suggest that ATP hydrolysis locks splice sites into a splicing pattern after stable U2 snRNP association to the branch point.

Exon definition as a potential negative force against intron losses in evolution

Background

Previous studies have indicated that the wide variation in intron density (the number of introns per gene) among different eukaryotes largely reflects varying degrees of intron loss during evolution. The most popular model, which suggests that organisms lose introns through a mechanism in which reverse-transcribed cDNA recombines with the genomic DNA, concerns only one mutational force.

Hypothesis

Using exons as the units of splicing-site recognition, exon definition constrains the length of exons. An intron-loss event results in fusion of flanking exons and thus a larger exon. The large size of the newborn exon may cause splicing errors, i.e., exon skipping, if the splicing of pre-mRNAs is initiated by exon definition. By contrast, if the splicing of pre-mRNAs is initiated by intron definition, intron loss does not matter. Exon definition may thus be a selective force against intron loss. An organism with a high frequency of exon definition is expected to experience a low rate of intron loss throughout evolution and have a high density of spliceosomal introns.

Conclusion

The majority of spliceosomal introns in vertebrates may be maintained during evolution not because of potential functions, but because of their splicing mechanism (i.e., exon definition). Further research is required to determine whether exon definition is a negative force in maintaining the high intron density of vertebrates.

Reviewers

This article was reviewed by Dr. Scott W. Roy (nominated by Dr. John Logsdon), Dr. Eugene V. Koonin, and Dr. Igor B. Rogozin (nominated by Dr. Mikhail Gelfand). For the full reviews, please go to the Reviewers' comments section.

The effect of intron length on exon creation ratios during the evolution of mammalian genomes

Recent studies report that alternatively spliced exons tend to occur in longer introns, which is attributed to the length constraints for splice site pairing for the two major splicing mechanisms, intron definition versus exon definition. Using genome-wide studies of EST and microarray data from human and mouse, we have analyzed the distribution of various subsets of alternatively spliced exons, based on their inclusion level and evolutionary history, versus increasing intron length. Alternative exons may be included in either a major or minor fraction of all transcripts (known as major-form and minor-form exons, respectively). We find that major-form exons are seven- to eightfold more likely to be contained in short introns (<400 nt) than minor-form exons, which occur preferentially in longer introns. Since minor-form exons are more likely to be novel (approximately 75%), this implied that novel exons arise more frequently in longer introns. To test this hypothesis, we used whole genome alignments to classify exons according to their phylogenetic age. We find that older exons, i.e., exons that are conserved in all mammals, predominate at shorter intron lengths, for both major- and minor-form exons. In contrast, exons that arose recently during primate evolution are more prevalent at longer intron lengths (>1000 nt). This suggests that the observed correlation of longer intron lengths with alternatively spliced exons may be at least partly due to biases in the probability of exon creation, which is higher in long introns.

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.

A general definition and nomenclature for alternative splicing events

Understanding the molecular mechanisms responsible for the regulation of the transcriptome present in eukaryotic cells is one of the most challenging tasks in the postgenomic era. In this regard, alternative splicing (AS) is a key phenomenon contributing to the production of different mature transcripts from the same primary RNA sequence. As a plethora of different transcript forms is available in databases, a first step to uncover the biology that drives AS is to identify the different types of reflected splicing variation. In this work, we present a general definition of the AS event along with a notation system that involves the relative positions of the splice sites. This nomenclature univocally and dynamically assigns a specific “AS code” to every possible pattern of splicing variation. On the basis of this definition and the corresponding codes, we have developed a computational tool (AStalavista) that automatically characterizes the complete landscape of AS events in a given transcript annotation of a genome, thus providing a platform to investigate the transcriptome diversity across genes, chromosomes, and species. Our analysis reveals that a substantial part—in human more than a quarter—of the observed splicing variations are ignored in common classification pipelines. We have used AStalavista to investigate and to compare the AS landscape of different reference annotation sets in human and in other metazoan species and found that proportions of AS events change substantially depending on the annotation protocol, species-specific attributes, and coding constraints acting on the transcripts. The AStalavista system therefore provides a general framework to conduct specific studies investigating the occurrence, impact, and regulation of AS.

The genome sequence is said to be an organism's blueprint, a set of instructions driving the organism's biology. The unfolding of these instructions—the so-called genes—is initiated by the transcription of DNA into RNA molecules, which subsequently are processed before they can take their functional role. During this processing step, initially identical RNA molecules may result in different products through a process known as alternative splicing (AS). AS therefore allows for widening the diversity from the limited repertoire of genes, and it is often postulated as an explanation for the apparent paradox that complex and simple organisms resemble in their number of genes; it characterizes species, individuals, and developmental and cellular conditions. Comparing the differences of AS products between cells may help to reveal the broad molecular basis underlying phenotypic differences—for instance, between a cancer and a normal cell. An obstacle for such comparisons has been that, so far, no paradigm existed to delineate each single quantum of AS, so-called AS events. Here, we describe a possibility of exhaustively decomposing AS complements into qualitatively different groups of events and a nomenclature to unequivocally denote them. This typological catalogue of AS events along with their observed frequencies represent the AS landscape, and we propose a procedure to automatically identify such landscapes. We use it to describe the human AS landscape and to investigate how it has changed throughout evolution.

"Nought may endure but mutability": spliceosome dynamics and the regulation of splicing

The spliceosome is both compositionally and conformationally dynamic. Each transition along the splicing pathway presents an opportunity for progression, pausing, or discard, allowing splice site choice to be regulated throughout both the assembly and catalytic phases of the reaction.

Conserved RNA secondary structures promote alternative splicing

Pre-mRNA splicing is carried out by the spliceosome, which identifies exons and removes intervening introns. Alternative splicing in higher eukaryotes results in the generation of multiple protein isoforms from gene transcripts. The extensive alternative splicing observed implies a flexibility of the spliceosome to identify exons within a given pre-mRNA. To reach this flexibility, splice-site selection in higher eukaryotes has evolved to depend on multiple parameters such as splice-site strength, splicing regulators, the exon/intron architecture, and the process of pre-mRNA synthesis itself. RNA secondary structures have also been proposed to influence alternative splicing as stable RNA secondary structures that mask splice sites are expected to interfere with splice-site recognition. Using structural and functional conservation, we identified RNA structure elements within the human genome that associate with alternative splice-site selection. Their frequent involvement with alternative splicing demonstrates that RNA structure formation is an important mechanism regulating gene expression and disease.

A DNA damage signal activates and derepresses exon inclusion in Drosophila TAF1 alternative splicing

Signal-dependent alternative splicing is important for regulating gene expression in eukaryotes, yet our understanding of how signals impact splicing mechanisms is limited. A model to address this issue is alternative splicing of Drosophila TAF1 pre-mRNA in response to camptothecin (CPT)-induced DNA damage signals. CPT treatment of Drosophila S2 cells causes increased inclusion of TAF1 alternative cassette exons 12a and 13a through an ATR signaling pathway. To evaluate the role of TAF1 pre-mRNA sequences in the alternative splicing mechanism, we developed a TAF1 minigene (miniTAF1) and an S2 cell splicing assay that recapitulated key aspects of CPT-induced alternative splicing of endogenous TAF1. Analysis of miniTAF1 indicated that splice site strength underlies independent and distinct mechanisms that control exon 12a and 13a inclusion. Mutation of the exon 13a weak 5' splice site or weak 3' splice site to a consensus sequence was sufficient for constitutive exon 13a inclusion. In contrast, mutation of the exon 12a strong 5' splice site or moderate 3' splice site to a consensus sequence was only sufficient for constitutive exon 12a inclusion in the presence of CPT-induced signals. Analogous studies of the exon 13 3' splice site suggest that exon 12a inclusion involves signal-dependent pairing between constitutive and alternative splice sites. Finally, intronic elements identified by evolutionary conservation were necessary for full repression of exon 12a inclusion or full activation of exon 13a inclusion and may be targets of CPT-induced signals. In summary, this work defines the role of sequence elements in the regulation of TAF1 alternative splicing in response to a DNA damage signal.

Identification of motifs that function in the splicing of non-canonical introns

The enrichment of specific intronic splicing enhancers upstream of weak PY tracts suggests a novel mechanism for intron recognition that compensates for a weakened canonical pre-mRNA splicing motif.

Background

While the current model of pre-mRNA splicing is based on the recognition of four canonical intronic motifs (5' splice site, branchpoint sequence, polypyrimidine (PY) tract and 3' splice site), it is becoming increasingly clear that splicing is regulated by both canonical and non-canonical splicing signals located in the RNA sequence of introns and exons that act to recruit the spliceosome and associated splicing factors. The diversity of human intronic sequences suggests the existence of novel recognition pathways for non-canonical introns. This study addresses the recognition and splicing of human introns that lack a canonical PY tract. The PY tract is a uridine-rich region at the 3' end of introns that acts as a binding site for U2AF65, a key factor in splicing machinery recruitment.

Results

Human introns were classified computationally into low- and high-scoring PY tracts by scoring the likely U2AF65 binding site strength. Biochemical studies confirmed that low-scoring PY tracts are weak U2AF65 binding sites while high-scoring PY tracts are strong U2AF65 binding sites. A large population of human introns contains weak PY tracts. Computational analysis revealed many families of motifs, including C-rich and G-rich motifs, that are enriched upstream of weak PY tracts. In vivo splicing studies show that C-rich and G-rich motifs function as intronic splicing enhancers in a combinatorial manner to compensate for weak PY tracts.

Conclusion

The enrichment of specific intronic splicing enhancers upstream of weak PY tracts suggests that a novel mechanism for intron recognition exists, which compensates for a weakened canonical pre-mRNA splicing motif.

Evidence against the energetic cost hypothesis for the short introns in highly expressed genes

BACKGROUND

In animals, the moss Physcomitrella patens and the pollen of Arabidopsis thaliana, highly expressed genes have shorter introns than weakly expressed genes. A popular explanation for this is selection for transcription efficiency, which includes two sub-hypotheses: to minimize the energetic cost or to minimize the time cost.

RESULTS

In an individual human, different organs may differ up to hundreds of times in cell number (for example, a liver versus a hypothalamus). Considered at the individual level, a gene specifically expressed in a large organ is actually transcribed tens or hundreds of times more than a gene with a similar expression level (a measure of mRNA abundance per cell) specifically expressed in a small organ. According to the energetic cost hypothesis, the former should have shorter introns than the latter. However, in humans and mice we have not found significant differences in intron length between large-tissue/organ-specific genes and small-tissue/organ-specific genes with similar expression levels. Qualitative estimation shows that the deleterious effect (that is, the energetic burden) of long introns in highly expressed genes is too negligible to be efficiently selected against in mammals.

CONCLUSION

The short introns in highly expressed genes should not be attributed to energy constraint. We evaluated evidence for the time cost hypothesis and other alternatives.

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.

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.

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.

Cross-kingdom patterns of alternative splicing and splice recognition

BACKGROUND

Variations in transcript splicing can reveal how eukaryotes recognize intronic splice sites. Retained introns (RIs) commonly appear when the intron definition (ID) mechanism of splice site recognition inconsistently identifies intron-exon boundaries, and cassette exons (CEs) are often caused by variable recognition of splice junctions by the exon definition (ED) mechanism. We have performed a comprehensive survey of alternative splicing across 42 eukaryotes to gain insight into how spliceosomal introns are recognized.

RESULTS

All eukaryotes we studied exhibit RIs, which appear more frequently than previously thought. CEs are also present in all kingdoms and most of the organisms in our analysis. We observe that the ratio of CEs to RIs varies substantially among kingdoms, while the ratio of competing 3' acceptor and competing 5' donor sites remains nearly constant. In addition, we find the ratio of CEs to RIs in each organism correlates with the length of its introns. In all 14 fungi we examined, as well as in most of the 9 protists, RIs far outnumber CEs. This differs from the trend seen in 13 multicellular animals, where CEs occur much more frequently than RIs. The six plants we analyzed exhibit intermediate proportions of CEs and RIs.

CONCLUSION

Our results suggest that most extant eukaryotes are capable of recognizing splice sites via both ID and ED, although ED is most common in multicellular animals and ID predominates in fungi and most protists.

Finding exonic islands in a sea of non-coding sequence: splicing related constraints on protein composition and evolution are common in intron-rich genomes

BACKGROUND

In mammals, splice-regulatory domains impose marked trends on the relative abundance of certain amino acids near exon-intron boundaries. Is this a mammalian particularity or symptomatic of exonic splicing regulation across taxa? Are such trends more common in species that a priori have a harder time identifying exon ends, that is, those with pre-mRNA rich in intronic sequence? We address these questions surveying exon composition in a sample of phylogenetically diverse genomes.

RESULTS

Biased amino acid usage near exon-intron boundaries is common throughout the metazoa but not restricted to the metazoa. There is extensive cross-species concordance as to which amino acids are affected, and reduced/elevated abundances are well predicted by knowledge of splice enhancers. Species expected to rely on exon definition for splicing, that is, those with a higher ratio of intronic to coding sequence, more introns per gene and longer introns, exhibit more amino acid skews. Notably, this includes the intron-rich basidiomycete Cryptococcus neoformans, which, unlike intron-poor ascomycetes (Schizosaccharomyces pombe, Saccharomyces cerevisiae), exhibits compositional biases reminiscent of the metazoa. Strikingly, 5 prime ends of nematode exons deviate radically from normality: amino acids strongly preferred near boundaries are strongly avoided in other species, and vice versa. This we suggest is a measure to avoid attracting trans-splicing machinery.

CONCLUSION

Constraints on amino acid composition near exon-intron boundaries are phylogenetically widespread and characteristic of species where exon localization should be problematic. That compositional biases accord with sequence preferences of splice-regulatory proteins and are absent in ascomycetes is consistent with selection on exonic splicing regulation.

Functioning of the Drosophila Wilms'-tumor-1-associated protein homolog, Fl(2)d, in Sex-lethal-dependent alternative splicing

fl(2)d, the Drosophila homolog of Wilms'-tumor-1-associated protein (WTAP), regulates the alternative splicing of Sex-lethal (Sxl), transformer (tra), and Ultrabithorax (Ubx). Although WTAP has been found in functional human spliceosomes, exactly how it contributes to the splicing process remains unknown. Here we attempt to identify factors that interact genetically and physically with fl(2)d. We begin by analyzing the Sxl-Fl(2)d protein-protein interaction in detail and present evidence suggesting that the female-specific fl(2)d(1) allele is antimorphic with respect to the process of sex determination. Next we show that fl(2)d interacts genetically with early acting general splicing regulators and that Fl(2)d is present in immunoprecipitable complexes with Snf, U2AF50, U2AF38, and U1-70K. By contrast, we could not detect Fl(2)d complexes containing the U5 snRNP protein U5-40K or with a protein that associates with the activated B spliceosomal complex SKIP. Significantly, the genetic and molecular interactions observed for Sxl are quite similar to those detected for fl(2)d. Taken together, our findings suggest that Sxl and fl(2)d function to alter splice-site selection at an early step in spliceosome assembly.

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

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

Polypyrimidine tract binding protein controls the transition from exon definition to an intron defined spliceosome

The polypyrimidine tract binding protein (PTB) binds pre-mRNAs to alter splice-site choice. We characterized a series of spliceosomal complexes that assemble on a pre-mRNA under conditions of either PTB-mediated splicing repression or its absence. In the absence of repression, exon definition complexes that were assembled downstream of the regulated exon could progress to pre-spliceosomal A complexes and functional spliceosomes. Under PTB-mediated repression, assembly was arrested at an A-like complex that was unable to transition to spliceosomal complexes. Trans-splicing experiments indicated that, even when the U1 and U2 small nuclear ribonucleoprotein particles (snRNPs) are properly bound to the upstream and downstream exons, the presence of PTB prevents the interaction of the two exon complexes. Proteomic analyses of these complexes provide a new description of exon definition complexes, and indicate that splicing regulators can act on the transition between the exon definition complex and an intron-defined spliceosome.

Coevolutionary networks of splicing cis-regulatory elements

Accurate and efficient splicing of eukaryotic pre-mRNAs requires recognition by trans-acting factors of a complex array of cis-acting RNA elements. Here, we developed a generalized Bayesian network to model the coevolution of splicing cis elements in diverse eukaryotic taxa. Cross-exon but not cross-intron compensatory interactions between the 5' splice site (5'ss) and 3' splice site (3'ss) were observed in human/mouse, indicating that the exon is the primary evolutionary unit in mammals. Studied plants, fungi, and invertebrates exhibited exclusively cross-intron interactions, suggesting that intron definition drives evolution in these organisms. In mammals, 5'ss strength and the strength of several classes of exonic splicing silencers (ESSs) evolved in a correlated way, whereas specific exonic splicing enhancers (ESEs), including motifs associated with hTra2, SRp55, and SRp20, evolved in a compensatory manner relative to the 5'ss and 3'ss. Interactions between specific ESS or ESE motifs were not observed, suggesting that elements bound by different factors are not commonly interchangeable. Thus, the splicing elements defining exons coevolve in a way that preserves overall exon strength, allowing specific elements to substitute for loss or weakening of others.

Functional coupling of last-intron splicing and 3'-end processing to transcription in vitro: the poly(A) signal couples to splicing before committing to cleavage

We have developed an in vitro transcription system, using HeLa nuclear extract, that supports not only efficient splicing of a multiexon transcript but also efficient cleavage and polyadenylation. In this system, both last-intron splicing and cleavage/polyadenylation are functionally coupled to transcription via the tether of nascent RNA that extends from the terminal exon to the transcribing polymerase downstream. Communication between the 3' splice site and the poly(A) site across the terminal exon is established within minutes of their transcription, and multiple steps leading up to 3'-end processing of this exon can be distinguished. First, the 3' splice site establishes connections to enhance 3'-end processing, while the nascent 3'-end processing apparatus makes reciprocal functional connections to enhance splicing. Then, commitment to poly(A) site cleavage itself occurs and the connections of the 3'-end processing apparatus to the transcribing polymerase are strengthened. Finally, the chemical steps in the processing of the terminal exon take place, beginning with poly(A) site cleavage, continuing with polyadenylation of the 3' end, and then finishing with splicing of the last intron.

The effect of intron location on the splicing of BmKK2 in 293T cells

Previously reported results showed that the BmKK2's intron could be recognized and spliced in cultured HEK 293T cells. At the same time, a cryptic splicing site of BmKK2 gene was found in the second exon. Moreover, replacing BmKK2's intron with BmP03's intron (an artificial BmKK2-BmP03 mosaic gene) did not affect the intron's recognition and splicing, but increased the expression level of the toxin-GFP fusion protein (Cao et al., J Biochem Mol Toxicol 2006;20:1-6). In this investigation, the BmKK2's intron with 79 nucleotides length was artificially shifted from the 49th nt (the 17th Gly codon between the first base and the second base) to the 100th nt (the 34th Gly codon between the first base and the second base). Based on the constructed intron-splicing system, the results of RT-PCR and the western blotting analysis showed that the BmKK2's shifted-intron (named BmKK2-s) was not recognized and spliced correctly, but the cryptic splicing site of BmKK2 gene was still spliced in the second exon, which possibly indicated that locations of introns were very important to the recognition and splicing of introns, and splicing of introns was very much associated with the corresponding upstream and downstream exons. This result possibly provides evidence for splice-site recognition across the exons.

Genome organization and gene expression shape the transposable element distribution in the Drosophila melanogaster euchromatin

The distribution of transposable elements (TEs) in a genome reflects a balance between insertion rate and selection against new insertions. Understanding the distribution of TEs therefore provides insights into the forces shaping the organization of genomes. Past research has shown that TEs tend to accumulate in genomic regions with low gene density and low recombination rate. However, little is known about the factors modulating insertion rates across the genome and their evolutionary significance. One candidate factor is gene expression, which has been suggested to increase local insertion rate by rendering DNA more accessible. We test this hypothesis by comparing the TE density around germline- and soma-expressed genes in the euchromatin of Drosophila melanogaster. Because only insertions that occur in the germline are transmitted to the next generation, we predicted a higher density of TEs around germline-expressed genes than soma-expressed genes. We show that the rate of TE insertions is greater near germline- than soma-expressed genes. However, this effect is partly offset by stronger selection for genome compactness (against excess noncoding DNA) on germline-expressed genes. We also demonstrate that the local genome organization in clusters of coexpressed genes plays a fundamental role in the genomic distribution of TEs. Our analysis shows that-in addition to recombination rate-the distribution of TEs is shaped by the interaction of gene expression and genome organization. The important role of selection for compactness sheds a new light on the role of TEs in genome evolution. Instead of making genomes grow passively, TEs are controlled by the forces shaping genome compactness, most likely linked to the efficiency of gene expression or its complexity and possibly their interaction with mechanisms of TE silencing.

Dual-specificity splice sites function alternatively as 5' and 3' splice sites

As a result of large-scale sequencing projects and recent splicing-microarray studies, estimates of mammalian genes expressing multiple transcripts continue to increase. This expansion of transcript information makes it possible to better characterize alternative splicing events and gain insights into splicing mechanisms and regulation. Here, we describe a class of splice sites that we call dual-specificity splice sites, which we identified through genome-wide, high-quality alignment of mRNA/EST and genome sequences and experimentally verified by RT-PCR. These splice sites can be alternatively recognized as either 5' or 3' splice sites, and the dual splicing is conceptually similar to a pair of mutually exclusive exons separated by a zero-length intron. The dual-splice-site sequences are essentially a composite of canonical 5' and 3' splice-site consensus sequences, with a CAG|GURAG core. The relative use of a dual site as a 5' or 3' splice site can be accurately predicted by assuming competition for specific binding between spliceosomal components involved in recognition of 5' and 3' splice sites, respectively. Dual-specificity splice sites exist in human and mouse, and possibly in other vertebrate species, although most sites are not conserved, suggesting that their origin is recent. We discuss the implications of this unusual splicing pattern for the diverse mechanisms of exon recognition and for gene evolution.

Gene expression, intron density, and splice site strength in Drosophila and Caenorhabditis

In this paper we investigate the relationships among intron density (number of introns per kilobase of coding sequence), gene expression level, and strength of splicing signals in two species: Drosophila melanogaster and Caenorhabditis elegans. We report a negative correlation between intron density and gene expression levels, opposite to the effect previously observed in human. An increase in splice site strength has been observed in long introns in D. melanogaster. We show this is also true of C. elegans. We also examine the relationship between intron density and splice site strength. There is an increase in splice site strength as the intron structure becomes less dense. This could suggest that introns are not recognized in isolation but could function in a cooperative manner to ensure proper splicing. This effect remains if we control for the effects of alternative splicing on splice site strength.

Concurrent splicing and transcription are not sufficient to enhance splicing efficiency

The concept of a tight integration of transcription and splicing of mRNA precursors has been supported with increasing evidence in recent years. However, the mechanism and functional consequences of this integration remain largely unknown. We have examined how these processes impact upon one another when they occur together in HeLa nuclear extract. While both processes do in fact occur in parallel reactions in the extracts, we found no evidence that one process affects the other, under a variety of conditions tested. For example, neither the kinetics nor efficiency of splicing is significantly enhanced by de novo RNA polymerase II-mediated transcription, relative to that of presynthesized RNA added exogenously to the extract. Our results indicate that the act of transcription by RNA polymerase II in vitro is not sufficient to enhance splicing of the newly made RNA.

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.

Wobble splicing reveals the role of the branch point sequence-to-NAGNAG region in 3' tandem splice site selection

Alternative splicing involving the 3' tandem splice site NAGNAG sequence may play a role in the structure-function diversity of proteins. However, how 3' tandem splice site utilization is determined is not well understood. We previously demonstrated that 3' NAGNAG-based wobble splicing occurs mostly in a tissue- and developmental stage-independent manner. Bioinformatic analysis reveals that the nucleotide preceding the AG dinucleotide may influence 3' splice site utilization; this is also supported by an in vivo splicing assay. Moreover, we found that the intron sequence plays an important role in 3' splice site selection for NAGNAG wobble splicing. Mutations of the region between the branch site and the NAGNAG 3' splice site, indeed, affected the ratio of the distal/proximal AG selection. Finally, we found that single nucleotide polymorphisms around the NAGNAG motif could affect the splice site choice, which may lead to a change in mRNA patterns and influence protein function. We conclude that the NAGNAG motif and its upstream region to the branch point sequence are required for 3' tandem splice site selection.

A missense mutation in the chloride/proton ClC-5 antiporter gene results in increased expression of an alternative mRNA form that lacks exons 10 and 11. Identification of seven new CLCN5 mutations in patients with Dent’s disease

Mutations in the voltage-gated chloride/proton antiporter ClC-5 gene, CLCN5, are associated with Dent's disease, an X-linked renal tubulopathy. Our interest is to identify and characterize disease-causing CLCN5 mutations, especially those that alter the splicing of the pre-mRNA. We analyzed the CLCN5 gene from nine unrelated Spanish Dent's disease patients and their relatives by DNA sequencing. Pre-mRNA splicing analysis was performed by RT-PCR. Seven new mutations were identified, consisting of three missense mutations (C219R, F273L, and W547G), one splice-site mutation (IVS-2A > G), one deletion (976delG), and two non-sense mutations (Y140X and W314X). We found that missense mutation W547G also led to increased expression of a new alternative isoform lacking exons 10 and 11 that was expressed in several human tissues. In addition, we describe another novel CLCN5 splicing variant lacking exon 11 alone, which was expressed only in human skeletal muscle. We conclude that missense mutation W547G can also alter the expression levels of a CLCN5 mRNA splicing variant. This type of mutation has not been previously described in the CLCN5 gene. Our results support the importance of a routine analysis at the pre-mRNA level of mutations that are commonly assumed to cause single amino acids alterations.

Modulating role of RNA structure in alternative splicing of a critical exon in the spinal muscular atrophy genes

Humans have two nearly identical copies of the survival motor neuron (SMN ) gene, SMN1 and SMN2. Homozygous loss of SMN1 causes spinal muscular atrophy (SMA). SMN2 is unable to prevent the disease due to skipping of exon 7. Using a systematic approach of in vivo selection, we have previously demonstrated that a weak 5' splice site (ss) serves as the major cause of skipping of SMN2 exon 7. Here we show the inhibitory impact of RNA structure on the weak 5' ss of exon 7. We call this structure terminal stem-loop 2 (TSL2). Confirming the inhibitory nature of TSL2, point mutations that destabilize TSL2 promote exon 7 inclusion in SMN2, whereas strengthening of TSL2 promotes exon 7 skipping even in SMN1. We also demonstrate that TSL2 negatively affects the recruitment of U1snRNP at the 5' ss of exon 7. Using enzymatic structure probing, we confirm that the sequence at the junction of exon 7/intron 7 folds into TSL2 and show that mutations in TSL2 cause predicted structural changes in this region. Our findings reveal for the first time the critical role of RNA structure in regulation of alternative splicing of human SMN.

Compensatory relationship between splice sites and exonic splicing signals depending on the length of vertebrate introns

BACKGROUND

The signals that determine the specificity and efficiency of splicing are multiple and complex, and are not fully understood. Among other factors, the relative contributions of different mechanisms appear to depend on intron size inasmuch as long introns might hinder the activity of the spliceosome through interference with the proper positioning of the intron-exon junctions. Indeed, it has been shown that the information content of splice sites positively correlates with intron length in the nematode, Drosophila, and fungi. We explored the connections between the length of vertebrate introns, the strength of splice sites, exonic splicing signals, and evolution of flanking exons.

RESULTS

A compensatory relationship is shown to exist between different types of signals, namely, the splice sites and the exonic splicing enhancers (ESEs). In the range of relatively short introns (approximately, < 1.5 kilobases in length), the enhancement of the splicing signals for longer introns was manifest in the increased concentration of ESEs. In contrast, for longer introns, this effect was not detectable, and instead, an increase in the strength of the donor and acceptor splice sites was observed. Conceivably, accumulation of A-rich ESE motifs beyond a certain limit is incompatible with functional constraints operating at the level of protein sequence evolution, which leads to compensation in the form of evolution of the splice sites themselves toward greater strength. In addition, however, a correlation between sequence conservation in the exon ends and intron length, particularly, in synonymous positions, was observed throughout the entire length range of introns. Thus, splicing signals other than the currently defined ESEs, i.e., potential new classes of ESEs, might exist in exon sequences, particularly, those that flank long introns.

CONCLUSION

Several weak but statistically significant correlations were observed between vertebrate intron length, splice site strength, and potential exonic splicing signals. Taken together, these findings attest to a compensatory relationship between splice sites and exonic splicing signals, depending on intron length.

Different levels of alternative splicing among eukaryotes

Alternative splicing increases transcriptome and proteome diversification. Previous analyses aiming at comparing the rate of alternative splicing between different organisms provided contradicting results. These contradicting results were attributed to the fact that both analyses were dependent on the expressed sequence tag (EST) coverage, which varies greatly between the tested organisms. In this study we compare the level of alternative splicing among eight different organisms. By employing an EST independent approach we reveal that the percentage of genes and exons undergoing alternative splicing is higher in vertebrates compared with invertebrates. We also find that alternative exons of the skipping type are flanked by longer introns compared to constitutive ones, whereas alternative 5′ and 3′ splice sites events are generally not. In addition, although the regulation of alternative splicing and sizes of introns and exons have changed during metazoan evolution, intron retention remained the rarest type of alternative splicing, whereas exon skipping is more prevalent and exhibits a slight increase, from invertebrates to vertebrates. The difference in the level of alternative splicing suggests that alternative splicing may contribute greatly to the mammal higher level of phenotypic complexity, and that accumulation of introns confers an evolutionary advantage as it allows increasing the number of alternative splicing forms.

Structural and functional characterization of the bovine solute carrier family 27 member 1 (SLC27A1) gene

The Solute Carrier Family 27 Member 1 (SLC27A1) is an evolutionarily conserved protein involved in regulating the long chain fatty acid uptake into cells. It has been shown to be expressed in tissues undergoing rapid fatty acid metabolism such as heart, skeletal muscle and adipose tissues, but no expression is detected in liver. Here we report the molecular characterization of the bovine SLC27A1 gene and draw a comparison with orthologous genes of some monogastric species. The bovine SLC27A1 gene is organized in 13 exons and extends over more than 40 kb of genomic DNA. It codes for a protein of 646 amino acids with a predicted molecular weight of 71 kDa which has 92%, 88% and 88% similarity with the human, mouse and rat SLC27A1 proteins respectively. The bovine SLC27A1 RNA expression was high in heart, testis, nervous tissue and muscle and very low in liver. Surprisingly, adipose tissues showed very low RNA expression levels contrary to the results described for both human and mouse genes. On the other hand, discordances observed between the bovine SLC27A1 RNA and protein expression patterns suggest that complex regulation mechanisms may be involved in determining the final SLC27A1 protein levels in each tissue. Finally, we have identified an alternative transcript generated by exon skipping of exon 3 to 7 which could encode a cytosolic SLC27A1 isoform of approximately 37 kDa.

SR proteins: a foot on the exon before the transition from intron to exon definition

Two recent publications illuminate the evolution of alternative splicing, showing that a SR (serine-arginine-rich) protein that regulates alternative splicing in multicellular organisms is also found in a unicellular organism without alternative splicing, in which it can assist in the splicing of weak introns. Moreover, insertion of SR proteins into an organism lacking such proteins can restore the splicing of weak introns. These results imply that SR proteins had already facilitated the splicing of weak introns before the evolution of alternative splicing.

Defective splicing, disease and therapy: searching for master checkpoints in exon definition

The number of aberrant splicing processes causing human disease is growing exponentially and many recent studies have uncovered some aspects of the unexpectedly complex network of interactions involved in these dysfunctions. As a consequence, our knowledge of the various cis- and trans-acting factors playing a role on both normal and aberrant splicing pathways has been enhanced greatly. However, the resulting information explosion has also uncovered the fact that many splicing systems are not easy to model. In fact we are still unable, with certainty, to predict the outcome of a given genomic variation. Nonetheless, in the midst of all this complexity some hard won lessons have been learned and in this survey we will focus on the importance of the wide sequence context when trying to understand why apparently similar mutations can give rise to different effects. The examples discussed in this summary will highlight the fine 'balance of power' that is often present between all the various regulatory elements that define exon boundaries. In the final part, we shall then discuss possible therapeutic targets and strategies to rescue genetic defects of complex splicing systems.

Linking splicing to Pol II transcription stabilizes pre-mRNAs and influences splicing patterns

RNA processing is carried out in close proximity to the site of transcription, suggesting a regulatory link between transcription and pre-mRNA splicing. Using an in vitro transcription/splicing assay, we demonstrate that an association of RNA polymerase II (Pol II) transcription and pre-mRNA splicing is required for efficient gene expression. Pol II-synthesized RNAs containing functional splice sites are protected from nuclear degradation, presumably because the local concentration of the splicing machinery is sufficiently high to ensure its association over interactions with nucleases. Furthermore, the process of transcription influences alternative splicing of newly synthesized pre-mRNAs. Because other RNA polymerases do not provide similar protection from nucleases, and their RNA products display altered splicing patterns, the link between transcription and RNA processing is RNA Pol II-specific. We propose that the connection between transcription by Pol II and pre-mRNA splicing guarantees an extended half-life and proper processing of nascent pre-mRNAs.

A novel in vitro method to study transcription and splicing leads to the proposal that linking transcription by Pol II and pre-mRNA splicing guarantees an extended half-life and proper processing of nascent pre-mRNAs.

Phenotypic consequences of branch point substitutions

The branch point sequence (BPS) is a conserved splicing signal important for spliceosome assembly and lariat intron formation. BPS mutations may result in aberrant pre-mRNA splicing and genetic disorders, but their phenotypic consequences have been difficult to predict, largely due to a highly degenerate nature of the BPS consensus. Here, we have examined the splicing pattern of nine reporter pre-mRNAs that have previously been shown to give rise to human hereditary diseases as a result of single-nucleotide substitutions in the predicted BPS. Increased exon skipping and intron retention observed in vivo were recapitulated for each mutated pre-mRNA, but the reproducibility of cryptic splice site activation was lower. BP mutations in reporter pre-mRNAs frequently induced aberrant 3' splice sites and also activated a cryptic 5' splice site. Systematic mutagenesis of BP adenosines showed that in most pre-mRNAs, the expression of canonical transcripts was lower for BP transitions than BP transversions. Differential splicing outcome for transitions vs. transversions was abrogated or reduced if introns were truncated to 200 nt or less, suggesting that the nature of the BP residue is less critical for interactions across very short introns. Together, these results improve prediction of phenotypic consequences of point mutations upstream of splice acceptor sites and suggest that the overrepresentation of disease-causing adenosine-to-guanosine BP substitutions observed in Mendelian disorders is due to more profound defects of gene expression at the level of pre-mRNA splicing.

Whole-Genome Microarray in Arabidopsis Facilitates Global Analysis of Retained Introns

Alternative splicing (AS) is an important post-transcriptional regulatory mechanism that can increase protein diversity and affect mRNA stability. Different types of AS have been observed; these include exon skipping, alternative donor or acceptor site and intron retention. In humans, exon skipping is the most common type while intron retention is rare. In contrast, in Arabidopsis, intron retention is the most prevalent AS type (approximately 40%). Here we show that direct transcript expression analysis using high-density oligonucleotide-based whole-genome microarrays (WGAs) is particularly amenable for assessing global intron retention in Arabidopsis. By applying a novel algorithm retained introns are detected in 8% of the transcripts examined. A sampling of 14 transcripts showed that 86% can be confirmed by RT-PCR. This rate of detection predicts an overall total AS rate of 20% for Arabidopsis compared with 10-22% based on EST/cDNA-based analysis. These findings will facilitate monitoring constitutive and dynamic whole-genome splicing on the next generation WGA slides.