A General Definition and Nomenclature for Alternative Splicing Events
[['Alternative splicing ( AS ) is a fundamental molecular process regulating eukaryotic gene expression and involved in numerous human diseases [1]–[3] .', 'It is usually postulated as the main mechanism to augment protein diversity from a somehow limited set of protein coding genes [4] .', 'Consequently , over the recent years various large scale studies have been undertaken aiming at the exhaustive identification and analysis of AS events ( for recent reviews , see [5]–[7] ) .', 'Current estimations claim around 60–75% of human multi-exonic genes to undergo AS [4] , [8] , [9] .', 'Surprisingly , to some extent , the rigorous formalization of the concept of AS event and its categorization has received relatively little attention .', 'Traditionally , terms for only five kinds of AS events have been coined: exon skipping ( ES ) , mutually exclusive exons ( ME ) , intron retention ( IR ) , alternative donor ( AD ) and acceptor ( AA ) sites [10] .', 'However , currently available transcript evidence shows a plethora of variations in splicing patterns that involve multiple instances of these classical events in various combinations [11] .', 'Figure 1 and Figure S1 give some examples of AS patterns observed in the manually curated RefSeq annotation [12] .', 'Despite the ever growing availability of gene annotations the lack of a universal reference definition of AS and hence of the corresponding categories of AS events prevent AS databases ( e . g . , AEdb [13] , ASD [14] , ATD [15] , Hollywood [16] , PASDB [17] , SpliceNest [18] , PALS db [19] , SpliceDB [20] , AsMamDB [21] , HASDB [22] , ProSplicer [23] , EuSplice [24] , ASAPII [25] etc . … ) , from the automatic identification and update of the AS landscape that characterizes the transcriptome from a particular cell type or condition .', "Such a specific landscape may be revealing the underlying biological mechanisms responsible for the cell's phenotype .", 'Towards that end the challenges to be addressed are', '( i ) to define and identify single instances of AS events in complex exon–intron variations ,', '( ii ) to find an intuitive vocabulary to adequately characterize different AS events , and', '( iii ) to develop methods to efficiently identify and classify AS events from sets of annotated transcripts .', 'Concerning challenge', '( i ) Malko and co-workers proposed to combine the classical terms for each exon observed in a given annotation [26] .', 'While variations of each exon across the compared transcripts can be sufficiently described by this procedure , it does not permit an easy extension for splicing variations across the adjacent introns .', 'However , some splicing evidence ( e . g . , the structures depicted in Figure 1A and 1C ) suggests a relation between splice sites across the intron as by means of intron definition mechanisms [27] , [28] .', 'In another study , Nagasaki et al . propose to delineate AS events at exonic regions common to the compared transcripts [29] .', 'Whereas in this approach intron-defined splicing variations are coherently modelled , events that could be connected by exon definition [30] are by definition assumed to be independent and are treated separately ( as for instance in Figure 1B ) .', 'Furthermore , the separation at overlapping exonic positions does not allow for investigation of potential relations between transcription and splicing , i . e . , relative position of the initiation and polyadenylation site ( Figure 1D ) , for which increasing evidence is reported in literature [31] , [32] .', 'The EnsEmbl databases on splicing , transcript , and exon diversity ( ASD , ATD , and AEdb ) have recently extended their definition of AS events ( e . g . , “complex intron retention” ) in order to allow for additional modifications upstream and/or downstream of a central event .', 'However , this notation system still remains limited and fails to capture structures depicted in Figure 1E or larger .', 'Addressing problem', '( ii ) , only few attempts have been undertaken to univocally denote AS events .', 'Malko et al . [26] proposed strings composed of 5 letters identifying each classical event to redundantly describe the variability separately for each exon observed in a certain annotation ( e . g . , “—AD” for combined variable acceptor and donor sites , Figure 1B ) .', 'These 5-component strings naturally bear a high degree of redundancy as one is required for each different form of exon .', 'Furthermore , the picture of the exon–intron structure can not be inferred solely from these strings , as can be seen by the structures in Figure 1A and 1C producing identical results .', 'Nagasaki and co-workers proposed the so-called “bit matrices” , binary matrices to describe AS events where each row represents a transcript variant and each column represents a genomic position .', 'Each position of the matrix is filled by “1” and “0” according to whether the respective transcript variant exhibits an exon or not at the corresponding position .', 'Neighboring identical columns then are collapsed , such that variations in the exon–intron structure are represented non-redundantly as flip-flop changes .', 'This representation draws a pictorial “bitmap” of the exon–intron structure from compared transcripts .', 'Disadvantages are that the number of “bits” that have to be shown ( i . e . , the matrix area ) is relatively large even for simple events ( e . g . , 14 for the event in Figure 1B ) .', 'Therefore , a condensed encoding of the bits in 2-dimensional integer vectors has been proposed , which however looses transparency of the exon–intron structure .', 'Alternatively , the nomenclature of ASD/ATD/AEdb focuses up to a certain degree on the location of variations around a centric intron/exon up to a certain degree , but does not describe the relative connection between these variations .', 'For instance , a name as “II-5p3p” ( i . e . , “intron isoform with modification at the 3′- and 5′-end” ) cannot distinguish the cases depicted in Figure 1A and 1C .', 'Correspondingly , the term “EI-5p3p” is ambiguous considering the structures in Figure 1B and 1D .', 'The number of such ambiguities grows with the number of concatenated terms: four different splicing structures for instance match the term “CIR-EB-5p3p” .', 'Also , the identification of a “central event” becomes problematic in large splicing variations ( Figure 1E ) .', 'With respect to issue', '( iii ) , splicing graphs as a non-redundant data structure have gained popularity in AS over the recent years , but definitions vary across literature .', 'Capturing the 5′→3′ directionality of transcription , they naturally all form directed acyclic graphs ( DAGs ) .', 'Going back to [33] , matching ( parts of ) ESTs [22] , [34] , [35] have been used as nodes connected by edges representing the EST evidence , in order to cluster them and/or to allow the analysis of AS .', 'Heber and co-workers [35] subsequently collapse ( remove ) vertices with indegree ( i . e . , the number of inedges ) =\u200aoutdegree ( the number of outedges ) =\u200a1 .', 'Later on , two works from the same year proposed a graph structure where every vertex corresponds to a splice site and the connecting edges represent the intermediate exon/intron [36] , [37] , labelled according to the mRNA or EST evidence .', 'Another kind of graph uses exons as nodes instead of splice sites [38] .', 'Whereas intuitive for visualization , the graph structure may redundantly contain common exon flanks .', 'Other graph-based approaches on exon–intron structures described in literature use similar techniques [25] , [39]–[41] .', 'However , all these analyses focus exclusively on the four types of traditional AS events , and thus capture only a limited fraction of the splicing variation encompassed in the transcriptome—probably a main consequence of problem', '( i ) .', 'Indeed , without a universal definition of AS event , the retrieval of a single type of splicing variation requires to define its corresponding sub-graph pattern and to localize all occurrences of this pattern in the whole splicing graph .', 'Consequently , a comprehensive characterization of AS needs an exhaustive set of such ad hoc patterns , which explains why usually only 4–6 types of events are considered .', 'In this work , we propose a general definition of “AS event” and we present a novel notation based on the relative position of alternative exon boundaries to flexibly describe such events .', 'Unlike traditional nomenclatures , this generic notation system allows the assignment of a univocal “AS code” to identify any possible variation of the exon–intron structure between two or more transcripts , and thus provides a platform for the automatic and exhaustive extraction of such variations from a dataset of annotated genes .', 'Here , we also describe in detail the method implemented in AStalavista ( Alternative Splicing transcriptional landscape visualization tool ) for the dynamic characterization of AS events in splicing graphs .', 'AStalavista is accessible as a web server at ( http://genome . imim . es/astalavista ) [42] .', 'We have used AStalavista to characterize and compare the “landscape” of AS in different human reference annotations as well as in annotations of other metazoan species , i . e . , chimp ( Pan troglodytes ) , mouse ( Mus musculus ) , rat ( Rattus norvegicus ) , dog ( Canis familiaris ) , cow ( Bos taurus ) , chicken ( Gallus gallus ) , frog ( Xenopus tropicalis ) , zebrafish ( Danio rerio ) , honeybee ( Apis mellifera ) , fruitfly ( Drosophila melanogaster ) , and worm ( Caenorhabditis elegans ) .', 'In contrast to previous large-scale studies , our approach focuses on splicing structure variations rather than on ( sequence ) attributes of alternative exons/introns [43] , [44] .', 'Results indicate that while most AS events can be assigned to a few categories , the categorization of AS events in different structures is quite complex , with a plethora of minor AS configurations .', 'Relative frequencies of particular patterns change with respect to the corresponding annotation protocol , species-specific attributes and coding constraints of the respective locus , and we present computational studies that investigate the reasons behind these fluctuations .'], ['The concurrent and regulated molecular mechanisms of exon and intron definition are generally responsible for the splicing structure in a certain transcript variant .', 'Although case studies for the mechanics of intron and exon recognition are given in literature [27] , [28] , [30] , no general rule could ( yet ) be deduced .', 'Therefore , neither of the mechanisms can be excluded from occurring during the splicing process and both are to be considered in a generally robust definition of AS event that is applicable to any organism without being a priori restricted to exon or intron definition .', 'In order to allow for possible interactions of parts of the splicing machinery across all exons and introns when delimiting AS events in exon–intron variations , our definition of AS events is based on sites: given an annotation , i . e . , transcript sequences aligned to the genome , we use the terminus “site” to describe genomic locations of aligned exon boundaries ( Definition 1 ) .', 'We propose a novel notation system to allow a complete classification of AS events .', 'The general idea is to assign to any AS event a string-based “AS code” that describes the structure of the splicing variation in a concise and univocal manner .', 'AS events of the same type ( e . g . , exon skipping ) are given an identical code and thus can be classified in the same structural group .', 'The codes are built dynamically with respect to each observed splicing variation without the requirement of an a priori defined catalogue of putative AS events .', 'Our notation system is based on the relative position of the variable sites that are involved in the AS event and proceeds as follows: first , all the variable sites of an AS event ( see Definition 4 ) are considered in the order of their genomic position from 5′ to 3′ .', 'The indices i ∈ N+ defined by this relative order are assigned to the corresponding variable sites .', 'In addition , a symbol is attributed to each site depending on its type .', 'We use the alphabet Σ\u200a=\u200a{[ , ˆ , - , ]} , where “[” denotes a TSS , “ ˆ” a splice donor , “-” an acceptor , and “]” a PAS .', 'Therefore , each site is represented by a number ( the relative position i ) and a symbol ( identifying the type ) .', 'To describe one of the splicing structures resulting from an AS event , the number and the symbol of all of the sites that are used by the corresponding mRNA within the event are concatenated into a string .', 'The digit “0” is used if the transcript does not use any variable site ( for instance by skipping an exon ) .', 'The AS code of the event corresponds to the concatenation of these strings , separating the descriptions of the variants by a comma .', 'We order the strings according to the relative position of their first site .', 'Examples are presented in Figures 1 and 2 .', 'Using this notation , AS events with identical codes are structurally equivalent , e . g . , all exon skipping or all alternative donor events .', 'Moreover , a specific AS code can always be defined for any splicing variation , which guarantees the exhaustiveness of the notation system .', 'For instance , the nonconventional events observed in Figure 1 are assigned the codes ( A ) 1ˆ3- , 2ˆ4- , ( B ) 1–3ˆ , 2–4ˆ , and ( C ) 1ˆ4–5ˆ6–7ˆ8- , 2ˆ3- .', 'Globally , the distribution of AS events into distinct structural classes forms the landscape of alternative splicing encompassed in a given annotation .', 'AStalavista is a JAVA-based tool designed to extract and visualize the structural landscape of AS events as reflected by a given annotation .', 'The input is provided in GTF format , containing the genomic coordinates of exons in the transcripts ( and , optionally , the coordinates of the coding regions ) .', 'AStalavista can be applied to any species for delineating the AS landscape from a whole genome annotation , or to a subset of genes composed according to custom criteria .', 'The output depicts the AS landscape by giving a summary of all pairwise AS events grouped into structurally equal classes which are ranked according to their observed abundances .', 'The web server [42] ( http://genome . imim . es/astalavista ) has been upgraded and depicts the spectrum of AS structures as described in this manuscript , including variable TSSs/PASs as pointed out by Definition 3 and Definition 4 .', 'This means that it is now possible to investigate for instance potential correlations between AS and alternative transcription initiation .', 'Also , the number of species and reference annotations that are supported has been increased .', 'To assess the agreement of AS events predicted according to our definition with data available from public sources , we compared the output of AStalavista for 5 well studied genes with the events classified for these in recently published or updated databases ( Table 1 ) .', 'Since AStalavista is a method rather than a fixed database , the number of AS events that are predicted crucially depends on the transcript annotation ( s ) under consideration .', 'Therefore , we conducted a first comparison of events extracted by AStalavista from mRNA annotations in Genbank [46] with the EuSplice database that is based on gene annotations .', 'In another run , we enriched the input data by ESTs from dbEST [47] and compared the corresponding results to the EST-based databases ASD , ATD and Hollywood .', 'In order to make the number of events in AStalavista quantitatively comparable with the number of events from public databases , we disregarded in either case AS events predicted in correlation with alternative transcription initiation or polyadenylation .', 'Table 1 shows that AStalavista clearly finds more bona fide events in either dataset than is available from public databases .', 'We additionally set off to investigate the overlap of the events in a case study ( Figure S2 ) and found that in the FOXP2 gene AStalavista ( Figure S2A ) finds 5 out of 6 events reported by Hollywood ( Figure S2B ) and 2 out of 3 events in EuSplice ( Figure S2C ) : in one instance Hollywood marked an alternative splice donor with a very untypical sequence that is supported exclusively by 2 ESTs ( Figure_S2B ) , and in the other case EuSplice predicted a cryptic exon based on the alignment of 2 nt in an intronic stretch which subsequently is tagged with the warning “short exon” and excluded from the analysis on splice site sequences ( Figure S2C ) .', 'For those AStalavista events that are not retrieved from both reference databases ( 8 out of 10 for EuSplice and 19 out of 24 for Hollywood ) , we found in total 4 cases that—although the evidence is present in the reference database—have not been reported , probably due to a limitation of the applied classification scheme .', 'These cases are: 0 , 1–2ˆ3–4ˆ ( i . e . , the skipping of two consecutive exons in events 14 and 15 ) , 1–2ˆ , 3–4ˆ ( the mutually exclusive exons in event 23 ) and 1–2ˆ3- , 4- ( the skipping of an exon when an alternative downstream acceptor is used , event 24 ) .', 'We ran AStalavista on three human popular annotation datasets , namely RefSeq [12] , EnsEmbl [48] and Gencode [49] .', 'With our clustering method ( see Materials and Methods ) , the 25 , 170 RefSeq transcripts clustered into 18 , 334 loci , the 43 , 102 EnsEmbl transcripts into 22 , 303 loci , and the 1 , 352 coding transcripts of Gencode into 381 loci ( Table 2 ) .', 'The differences in the average number of coding transcripts per locus between these annotations ( 1 . 4 for RefSeq , 1 . 9 for EnsEmbl , and 3 . 6 for Gencode ) reflect the differences in exhaustiveness among them .', 'We extracted all variations of the exon–intron structures according to Definition 4 .', 'To compensate for artefacts that may occur in automatic annotation pipelines , we omitted AS events that involved introns with no canonical splice site dinucleotides ( i . e . , not GT/AG ) .', 'Note that this filtering step consumes a considerable part of the observed running time ( Table 2 ) , since for each intron the splice site nucleotides are extracted from the genomic sequence .', 'As expected , the observed running times reflect the number and distribution of transcripts in each input annotation and the longest run ( for EnsEmbl ) took a bit more than a minute ( Table 2 ) .', 'Next , we analyzed the transcript diversity by characterizing the AS landscapes produced by AStalavista from the different annotations ( Figure 3 ) .', 'To compare the results with other studies , we focused on the traditional AS events that present a “simple” splicing pattern—involving at most two alternative splice sites and not correlated with variable TSS/PAS .', 'Agreeing with previously reported observations [29] , these simple events are equally ranked from the most abundant to the less in all annotations data sets in the order: exon skipping ( ES ) , alternate donor ( AD ) , alternate acceptors ( AA ) and intron retention ( IR ) .', 'All other AS events are pooled together ( Figure 3 , grey sectors in the pie diagrams ) .', 'These “complex” events form as a whole a substantial part of the AS landscape ( from 23 . 18% in RefSeq up to 35 . 4% in EnsEmbl ) , and each of them can be unambiguously described by the notation proposed herein .', 'The composition of these events varies ( Table S1 ) : the 1 , 070 AS events detected in RefSeq correspond to 85 structural distinct classes ( A ) , whereas the 4 , 321 events in EnsEmbl show 388 classes ( C ) .', 'The fairly most abundant of these complex event ( from 25 . 6% of them in EnsEmbl to 32 . 6% in Gencode ) is the skipping of two exons in a row ( 0 , 1–2ˆ3–4ˆ ) .', 'Mutually exclusive exons ( 1–2ˆ , 3–4ˆ ) are less frequent ( from 12% to 14 . 5% ) , probably due to a more complex molecular mechanism that regulates them .', 'As expected , the higher the complexity of an event—as measured by the number of splice sites involved—the lower its relative abundance .', 'For instance , the “triple exon skipping” ( 0 , 1–2ˆ3–4ˆ5–6ˆ ) forms ∼7–9% of the complex events .', 'The fact that this event still represents 93 reported cases in the RefSeq annotation ( Table S1 ) illustrates the need for an exhaustive AS notation system and for the corresponding retrieval method .', 'Obviously , there are differences in the AS landscape between the different reference annotations .', 'This probably reflects the differences in biological data and in the annotation process: manually reviewed full-length cDNA sequences in RefSeq , automatically annotated proteins/cDNAs in EnsEmbl and manually annotated transcripts including ESTs evidence augmented by experimentally verified computational predictions in Gencode .', 'Nevertheless , the different proportions of events agrees with previous results ( e . g . , [29] , [37] ) and their ranking is consistent across the sets , which illustrates the general consistence in the AS taxonomies reflected by these annotation systems .', 'Particularly relevant is , in our opinion , the consistency in the AS landscape between the RefSeq and the much richer Gencode annotation .', 'Even though Gencode contains 2 . 5-fold the number of alternative transcripts per locus , it includes only a marginally larger proportion of the “other” complex AS events than the conservative RefSeq , indicating that while only a fraction of the protein coding transcripts in the human genome may be currently known , the broad AS landscape characterizing the RefSeq annotation is also likely to characterize the entire human transcript complement .', 'We have investigated the differences in the type of AS events occurring in the CDS ( coding sequence ) from those occurring only in the 5′ UTR ( 5′ untranslated region ) .', 'Figure 4 shows the distribution of the simple AS events in 5′ UTRs and in CDSs from the RefSeq annotation .', 'The distribution in 3′ UTRs ( 3′ untranslated regions ) is not shown because of the low frequency of ( alternative ) splicing in these regions .', 'The analysis here focuses on events completely included in a certain region ( see Methods ) —i . e . , in the 5′ UTRs or in the CDSs—but the same trends can be observed for events overlapping the 5′ UTR and the CDS ( Figure S3 ) .', 'The distributions differ even in the ranking of the four most abundant events .', 'In agreement with [29] , the proportion of ES is significantly higher in CDS ( 50 . 9% of the landscape ) than in the 5′ UTR ( 37 . 9% , p-value<10−4 , χ2 test ) .', 'A straightforward explanation is the fact that ES requires at least two introns , which are present in a minority of 5′ UTRs .', 'Coherently with this explanation , we observe the following low proportions of complex AS events in 5′ UTRs vs . CDS: 26 . 6% vs . 33 . 1% for 0 , 1–2ˆ3–4ˆ events ( skipping of two exons ) , 10 . 6% vs . 17 . 3% for 1–2ˆ , 3–4ˆ events ( mutually exclusive exons ) and—more drastically—1 . 3% vs . 10 . 7% for 0 , 1–2ˆ3–4ˆ5–6ˆ ( the joint skipping of 3 neighboring exons ) events .', 'Expectedly , since retained introns in CDSs are likely to introduce in-frame stop codons , the relative proportion of IR is much higher in 5′ UTR ( 8 . 4% vs . 2 . 1% , p-value<10−4 , χ2 test ) .', 'Strikingly , the relative frequency of AA and AD events shows a “reciprocal asymmetry” between the CDS and the 5′ UTRs .', 'In the CDS , the proportion of AAs is nearly twice as high as the proportion of ADs ( 14 . 8% vs . 8% ) , while in the 5′ UTR regions the ratio is the other way around ( 13 . 7% vs . 22 . 5% ) .', 'Considering findings on the possibly differing molecular mechanism for short range variations at the donor and acceptor site [50]–[52] , we repeated the analysis disregarding variations between AD or AA shorter than 5 bp and found a comparable asymmetry ( data not shown ) .', 'The bias against AAs in 5′ UTRs can be explained by the shorter sequence span where alternate acceptor sites can appear without disrupting the downstream protein sequence .', 'Indeed , if we consider the 5′ UTRs that contain exactly one intron ( 75% of the spliced 5′ UTRs ) , the length of the potential target for alternative upstream donor site creation , that is the first exon , is significantly larger than the length of the potential target for alternative downstream acceptor sites creation in 5′ UTR , that is from the acceptor site to the ATG codon ( 260 vs . 47 nucleotides on average ) .', 'In order to confirm that the bias against AAs in the 5′ UTR is mainly due to constraints of the start codon , we considered in multi-intronic 5′ UTRs the AS events that do not affect the last intron .', 'Then , the AD/AA ratio drops from factor >1 . 64 to factor 1 . 2 ( 30 AD events compared to 25 AA events in RefSeq ) .', 'In our opinion , the remaining polarity stems from the fact that the first exon is significantly longer than the second ( median 149 vs . 137 , p-value ∼3e-6 , Kolmogorov-Smirnov-Test ) , probably resulting from differences in the mechanism for exon definition [27] .', 'On the other hand , the observed asymmetry against ADs in the CDS can be explained by the propensity towards the creation of stop codons when considering alternative downstream donor sites , due to the peculiar composition of the donor site consensus sequence .', 'As already reported in the past [53] , splicing consensus sequences harbor a high content of intrinsic stop codons ( shaded grey in Figure 5A and 5B ) .', 'To test this hypothesis , we have artificially extended constitutively chosen exon boundaries into the intronic flanks and measured the frequency of in-frame stop codon occurrence separately for the 5′ and the 3′ end .', 'As summarized in Figure 5 , the inclusion of one additional codon from the intronic sequence already interrupts the CDS at the donor site ∼50% more often than at the acceptor site .', 'Interestingly , another—though lower—peak of potential stops at the acceptor site is observed after ∼9 codons of extension and coincides with the common location of the branch point consensus ( Figure 5C ) .', 'This difference of potentially included stop codons biases against ADs up to 22 codons of extension ( Figure S4 ) and therefore gives strong evidence for the more frequent use of AAs at flanks of coding exons—albeit more complex mechanisms are also expected to play an additional role .', 'Additional evidence of the strong effects of the protein coding constraints in shaping the AS landscape comes from the comparison of AS in protein coding and noncoding transcripts .', 'For this comparison , the Gencode annotation is particularly appropriate: it contains many non protein-coding transcripts ( 2 , 247 vs . 1 , 332 coding transcripts ) , most of them actually occurring also in protein coding loci .', 'In other words , protein coding loci seem to be able to encode both , protein coding and noncoding transcripts .', 'Figure 6 shows the distribution of the AS events in protein coding regions ( i . e . , in the CDSs ) and in noncoding transcripts .', 'The differences are substantial , interestingly also in comparison to the AS events in 5′ UTRs ( Figure 4 ) , not biased by the difference in size between the datasets ( Figure S5 ) .', 'Almost one third ( 31 . 5% ) of the AS events observed in noncoding transcripts correspond to complex splice events , compared to only about one fourth ( 24 . 3% ) in CDSs .', 'Also , the composition of the complex fraction in noncoding transcripts is richer ( 57 structural different classes vs . 22 in CDSs ) .', 'Consequently , simple events that are frequently reported in the CDSs of Gencode transcripts are relatively less abundant in noncoding transcripts ( e . g . , from 48 . 5% to 34 . 4% for exon skipping ) .', 'Naturally , we observe a relaxation of selective constraints against retained introns that make up ∼12% of the landscape in transcripts without an annotated reading frame .', 'The AA/AD ratio is more balanced in noncoding transcripts ( 1 . 6 vs . 2 . 6 in CDSs ) .', 'The remaining polarity stems from asymmetries in the first compared to the last intron: whereas an alternative TSS in the first exon is often associated with an alternative first donor site ( 87 instances ) , an alternative acceptor site in the last exon is less frequently observed with a different PAS ( 56 cases ) .', 'When taking into account such events , the numbers for variable 5′ and 3′ flanks of exons are about equal ( 150 ADs and 159 AAs ) .', 'This indeed underlines the very different selective constrains acting on coding and noncoding transcripts—even though they may be extensively sharing the same genomic space .', 'To investigate the evolution of the AS landscape , we have applied AStalavista to the annotation of 12 different metazoan genomes: human ( Homo sapiens ) , chimp ( Pan troglodytes ) , mouse ( Mus musculus ) , rat ( Rattus norvegicus ) , dog ( Canis familiaris ) , cow ( Bos taurus ) , chicken ( Gallus gallus ) , frog ( Xenopus tropicalis ) , zebrafish ( Danio rerio ) , honeybee ( Apis mellifera ) , fruitfly ( Drosophila melanogaster ) , and worm ( Caenorhabditis elegans ) .', 'While many of the fluctuations observed are likely due to the species-specific differences in amount and quality of the transcriptional data from which the annotations have been derived , our study reveals some interesting trends , suggesting overall that AS patterns did not change gradually but rather abruptly during metazoan evolution ( Figure 7 ) .', 'More specifically , IR events are clearly more abundant in invertebrates than in vertebrates .', 'This is consistent with the fact that invertebrates have much shorter introns .', 'Indeed , one could think that IR events involving short introns are less likely to be negatively selected , since the probability for the protein sequence to get disrupted by the introduction of a stop codon is lower than with long introns ( Table S2 ) .', 'On the other hand , vertebrates—and especially mammals—exhibit a higher proportion of ES events , while , in contrast , relying relatively less on the usage of alternative donors and acceptors .', 'This may reflect a higher level of regulation of AS in vertebrates , possibly correlated with a higher frequency of exon shuffling and protein domains rearrangements [54] .', 'Finally , we observe an accumulation of complex events in vertebrate genomes compared to the invertebrates ( Figure 7 ) .', 'This could be due to the larger number of exons per gene on average in vertebrate genomes ( Table S3 ) , which allows to increase the combinatory level , but it also suggest a higher level of sophistication in the control of AS in vertebrate genomes when compared to invertebrates .'], ['Alternative Splicing increases enormously the encoding capacity of the genome of the higher eukaryotic organisms .', 'Its differential regulation is likely to play a substantial role in defining the phenotype of a given cell type , or cell state .', 'We have developed a method to automatically catalogue the patterns of AS events occurring in a given gene/transcript annotation .', 'The method ( and the resulting ) taxonomy relies on a precise definition of AS event .', 'We have implemented the method in a publicly available software system , named AStalavista .', 'As a proof of concept , the application of AStalavista to a number of popular annotations of the human genomes has revealed the existence of a plethora of AS types that are usually ignored in published analyses .', 'Indeed , about one quarter of all AS events in these collections belong to this category .', 'Some of these complex AS events , like double exon skipping or mutually exclusive exons , are likely to be under specific regulation .', 'In addition , we report notable differences in the AS landscape between coding and noncoding regions and transcripts , with the landscape in coding regions being largely modelled by protein coding constraints and the landscape in noncoding transcripts suggesting a relaxation of selective constraints .', 'Our comparison of the AS landscape across 12 metazoan genomes reveals strong differences between vertebrate and non-vertebrate genomes .', 'We observe a higher fraction of intron retention events in invertebrates , while in contrast exon skipping and complex splicing events are more prevalent in vertebrates .', 'While the latter could simply reflect the richer transcript data available for vert
['Introduction', 'Results', 'Discussion', 'Materials and Methods']
['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 .']
['computational biology/alternative splicing', 'computational biology/genomics', 'computational biology/transcriptional regulation']