a substantial fraction of the genomic dna sequence does not directly encode the primary structure of any cellular protein , or any other cellular product ( 1 ) . this is largely due , in eukaryotes , to the division of genes into introns ( noncoding parts of dna ) and exons ( coding parts of dna ) , each of which have pronounced size distributions ( 2 ) . the actual mechanism ( or mechanisms ) of intron insertion is the subject of intense discussion and a currently popular working hypothesis suggests that introns may be largely produced by insertion of transposons dna elements that can move around to different positions within the genome ( 3 ) . very frequently , this point of view implicitly assumes that there is a higher probability of splitting longer exons since they are larger targets for transposons . here , we show that the distributions of exon sizes for different organisms have a general property : the presence of two distinguishable classes of exons with different size distributions . we present an idealized scheme that explains how the observed distribution of exon sizes can be derived from a common ancestral genome by a random ( quasi)-evolutionary process . this formal model makes it possible to investigate the evolution of the genomes of particular organisms , and to estimate the number of evolutionary steps that separate them from the hypothetical ancestral genome . conceptually , our results support the opinion that at the initial stages of evolution , simple genomes had a lower fraction of introns ( introns late hypothesis ) . the model we propose explains the observed lognormal distribution of exon sizes , and suggests that introns may be inserted by a process that is independent of exon size . our findings also can be rationalized in relation to the phenomenon of alternative splicing ( 4 ) . in our study , we used genome data for 12 animal species provided by ensembl ( 5 ) ( http://www.ensembl.org/ ) , the joint project of the european molecular biology laboratory european bioinformatics institute and the sanger institute , and a plant genome sequence from the arabidopsis information resource ( tair , http://www.arabidopsis.org/ ) ( 6 ) . all duplicates of exons ( which originate primarily from multiple accessions ) starting and ending at the same points were excluded from the analysis . the natural logarithms of exon sizes were divided into bins of width  ln(exon size ) = 0.2 to obtain the exon size distributions presented in figure 1 and in supplementary data . the distributions were fit using the unweighted -measure normalized over the number of degrees of freedom ( 7 ) as the criterion of fitting quality . table 2 shows the corresponding values of  and the pearson correlation coefficient , r , which characterize the agreement between original data and the fitted model . figure 1.distributions of the natural logarithms of exon sizes for the homo sapiens and drosophila melanogaster genomes . both data sets can be approximated by a sum of two gaussian peaks with very high correlation between the data and the best - fit distributions : pearson correlation coefficients r = 0.997 and r = 0.998 for homo sapiens and drosophila melanogaster genomes , respectively . data are shown as gray bars ; thick lines represent the best fit approximation , while thin lines depict individual locations and shapes of two gaussian peaks . table 1.parameters of fitting for different models of exon size distributionspecies namescaenorhabditis elegansarabidopsis thalianaanopheles gambiaetetraodon nigroviridisbos taurusdrosophila melanogastercanis familiarisgallus gallusdanio reriomacaca mulattamus musculushomo sapienspan troglodytesexponential peakn  028  1.533  1.712  0.548  3.552  3.313  0.552  3.444  2.964  4.056  3.652  3.967  4.855  3.9  105.3  0.45.1  0.43.6  0.26.7  0.77.7  0.72.8  0.27.5  0.77.6  0.77.7  0.77.3  0.77.0  0.76.9  0.77.2  0.7lognormal peakn  1023  130  111  136  138  113  138  133  147  142  139  150  141 11.3  0.11.8  0.11.7  0.11.0  0.11.1  0.12.0  0.11.0  0.11.0  0.11.00.11.1  0.11.0  0.11.0  0.11.0  0.1m5.0  0.15.0  0.15.4  0.14.9  0.14.8  0.15.6  0.14.8  0.14.8  0.14.8  0.14.8  0.14.8  0.14.8  0.14.8  0.1exponential peak plus weibull peakn  1013  225  19.3  0.318  112  111  0.314  111  122  214  113  118  113  1  103.5  0.53.7  0.22.7  0.14.3  0.34.1  0.32.1  0.14.7  0.34.2  0.35.7  0.43.0  0.22.2  0.22.1  0.22.5  0.2nw  1032  440.8  2.20.1  0.23.5  2.762  160.1  0.133  1046  1123  1273  1931  1259  2453  17c2.3  0.33.2  0.63.1  0.32.9  0.22.4  0.13.1  0.42.5  0.12.4  0.12.6  0.12.4  0.12.5  0.12.5  0.12.4  0.1w  1010  200.3  0.90.1  0.20.5  0.48.6  2.30.1  0.24.5  1.37.3  1.72.8  1.59.2  2.53.8  1.55.7  2.46.3  2.1two lognormal peaksna 105  0.89  0.91.3  0.121  1.321  1.51.5  0.120  1.417  1.2261.2281.129  0.937  0.930  0.9nb  1018  0.921  1.010  0.118  1.322  1.511  0.123  1.420  1.2291.1211.216  1.123  1.119  1.0a0.7  0.10.8  0.10.6  0.10.8  0.12.1  0.10.5  0.12.1  0.12.2  0.12.50.10.80.10.8  0.10.9  0.10.8  0.1b1.5  0.11.8  0.11.9  0.11.9  0.10.8  0.12.1  0.10.8  0.10.8  0.10.80.12.70.23.1  0.23.1  0.13.1  0.2ma4.7  0.14.5  0.15.1  0.14.9  0.14.6  0.15.1  0.14.6  0.14.6  0.14.50.14.80.14.8  0.14.8  0.14.8  0.1mb5.2  0.05.4  0.05.5  0.05.0  0.04.8  0.05.8  0.04.8  0.04.8  0.04.80.04.90.15.4  0.15.5  0.15.1  0.1 
 table 2.fitting quality ,  and r , together with the number of data points and number of degrees of freedom , dfspecies namescaenorhabditis elegansarabidopsis thalianaanopheles gambiaetetraodon nigroviridisbos taurusdrosophila melanogastercanis familiarisgallus gallusdanio reriomacaca mulattamus musculushomo sapienspan troglodytesnumber of data points50455150515151505151515151exponential peak  103081357725914 94713 99131614 87610 87419 82616 19618 55929 37318 882r0.8310.8470.9100.7480.7770.9020.7710.7700.7860.7730.7250.7310.743df48434948494949484949494949lognormal peak 10301132469107115791131540140034002014199537492148r0.9840.9450.9770.9820.9750.9660.9770.9710.9640.9720.9710.9660.971df47424847484848474848484848exponential peak plus weibull peak  1075710434261216153145934592805481135476r0.9610.9580.9860.9900.9980.9850.9980.9980.9950.9960.9920.9900.994df45404645464646454646464646two lognormal peaks  10534881623086256244341277322342264r0.9970.9980.9980.9970.9960.9980.9960.9950.9970.9960.9960.9970.997df44394544454545444545454545 distributions of the natural logarithms of exon sizes for the homo sapiens and drosophila melanogaster genomes . both data sets can be approximated by a sum of two gaussian peaks with very high correlation between the data and the best - fit distributions : pearson correlation coefficients r = 0.997 and r = 0.998 for homo sapiens and drosophila melanogaster genomes , respectively . data are shown as gray bars ; thick lines represent the best fit approximation , while thin lines depict individual locations and shapes of two gaussian peaks . parameters of fitting for different models of exon size distribution fitting quality ,  and r , together with the number of data points and number of degrees of freedom , df three classes of fitting models were used : 
 a model based on the assumption that an exon can be split at any position with equal and constant probability ; this model leads to an exponential distribution of exon sizes : 

 where e is exon size , dn is the number of exons in a bin , n is the amplitude of the peak and  is the probability of splitting an exon at a particular place.two models that produce lognormal distributions of exon sizes . these models are based on a kolmogoroff process , which does not assume any relationship between exon size and probability of splitting an exon . particularly , to fit the data we used a model of single lognormal peak 

 and a mixture of two lognormal distributions , 

 where , n , na and nb are the amplitudes , m , ma and mb the mean positions , and  , a and b the variances of lognormal peaks observed in the data.a combination of a weibull distribution and exponential distribution 

 where nw , w and c are the amplitude , frequency parameter and shape parameter of the weibull distribution , respectively . 
  a model based on the assumption that an exon can be split at any position with equal and constant probability ; this model leads to an exponential distribution of exon sizes : 

 where e is exon size , dn is the number of exons in a bin , n is the amplitude of the peak and  is the probability of splitting an exon at a particular place . these models are based on a kolmogoroff process , which does not assume any relationship between exon size and probability of splitting an exon . particularly , to fit the data we used a model of single lognormal peak 

 and a mixture of two lognormal distributions , 

 where , n , na and nb are the amplitudes , m , ma and mb the mean positions , and  , a and b the variances of lognormal peaks observed in the data . a combination of a weibull distribution and exponential distribution 

 where nw , w and c are the amplitude , frequency parameter and shape parameter of the weibull distribution , respectively . sequences of pseudorandom numbers were obtained using the mersenne twister algorithm ( 8) implemented in the standard matlab 7.1 installation . the sequences of pseudorandom numbers were seeded with the values 14 ( sequence a ) and 16 ( sequence b ) , and were used to generate the multiplicative processes presented in figure 2 . each pseudorandom sequence consists of 10 random numbers , i.e. , 10 exon splitting events . the ratio between the initial lengths of the two ancestral exons used to generate the exon size distributions presented in figure 2 was 1:1000 . figure 2.distributions of the natural logarithm of exon sizes for two different realizations of a random multiplicative process starting from an ancestral genome comprising two exons . the diagrams at the top of the figure illustrate two different splitting patterns for sequences a and b ( see methods section ) . the peaks shown with gray bars in the plot originated from the short exon in ancestral genome ( also gray in diagram ) ; those shown in black bars originated from the long exon in ancestral genome ( also black in diagram ) . the probability , p , of a splitting event for a single exon at each step of the process is also indicated . for the three first steps of the pseudorandom number sequence shown in ( a ) , all splitting events happened to occur for the exons that originated from the short part of the model ancestral genome . thus , by the third step of the process , four out of five exons originated from this part of model ancestral genome . this subset of exons has largest cumulative probability ( p = 4/5 ) for the next splitting event . for the sequence in ( b ) , at the third step of the process the largest cumulative probability for the next splitting event ( p = 3/5 ) belongs to the subset of exons that originated from the long part of ancestral exon . this initial break of symmetry between the two subsets of exons in the model genome persists and produces the two different peak dispositions shown in panels ( a ) and ( b ) . distributions of the natural logarithm of exon sizes for two different realizations of a random multiplicative process starting from an ancestral genome comprising two exons . the diagrams at the top of the figure illustrate two different splitting patterns for sequences a and b ( see methods section ) . the peaks shown with gray bars in the plot originated from the short exon in ancestral genome ( also gray in diagram ) ; those shown in black bars originated from the long exon in ancestral genome ( also black in diagram ) . the probability , p , of a splitting event for a single exon at each step of the process is also indicated . for the three first steps of the pseudorandom number sequence shown in ( a ) , all splitting events happened to occur for the exons that originated from the short part of the model ancestral genome . thus , by the third step of the process , four out of five exons originated from this part of model ancestral genome . this subset of exons has largest cumulative probability ( p = 4/5 ) for the next splitting event . for the sequence in ( b ) , at the third step of the process the largest cumulative probability for the next splitting event ( p = 3/5 ) belongs to the subset of exons that originated from the long part of ancestral exon . this initial break of symmetry between the two subsets of exons in the model genome persists and produces the two different peak dispositions shown in panels ( a ) and ( b ) . the confidence intervals shown for the fitted parameters in figure 1 , error bars in figure 3 , tables 1 and 3 , and in the supplementary data are 95% confidence intervals estimated from the covariance matrix ( 7 ) . figure 3.widths ,  , of exon size distributions for several complete genomes ( see table 1 and in supplementary data ) . all the species were ordered according to the width of wider peak in the exon distributions . widths ,  , of exon size distributions for several complete genomes ( see table 1 and in supplementary data ) . all the species were ordered according to the width of wider peak in the exon distributions . these data were received from the third release of the alternative splicing database ( asd ) , ( http://www.ebi.ac.uk/asd/altsplice/ ) . this resource provides manually curated data on alternative splicing with all exons confirmed by est / mrna alignments ( 9,10 ) . these data were also divided into bins of width  ln ( exon size ) = 0.2 and fitted using the unweighted -measure to models of single exponential peak , single weibull peak , and single lognormal distribution ( figure 4 and table 3 ) . figure 4.top row presents comparison between statistical distributions of all exons ( gray bars ) with the distributions of alternatively spliced exons ( black bars ) from asd . the data for homo sapiens are in left panes ; for drosophila melanogaster in right panels . panels in bottom row show the fit of the alternatively spliced exon length distribution to a single lognormal peak ( full lines ) as well as indicate locations of narrow ( dashed lines ) and wide ( dotted lines ) peaks obtained for the statistical distributions of all exons . inserts in bottom panels show values of correlations coefficients between the statistical distributions of alternatively spliced exons ( ae ) and fitted curve ( fit ) ; narrow peak ( np ) and wide peak ( wp ) . 
 table 3.parameters of fitting for different models approximating size distributions of alternatively spliced exons ( see methods sections)species namesmus musculus 28 data pointshomo sapiens 35 data pointsexponential peakdf = 26df = 33n  1033.3  3.739.5  3.9  108.1  1.37.9  1.1  1016 88018 313r0.6480.677weibull peakdf = 25df = 32nw  1040  2251  22c2.4  0.12.4  0.1w  106.9  3.27.4  3.5 10522756r0.9900.987lognormal peakdf = 25df = 32n  1023.5  0.428.2  0.10.9  0.10.9  0.1m4.8  0.14.8  0.1  10378329r0.9920.994 top row presents comparison between statistical distributions of all exons ( gray bars ) with the distributions of alternatively spliced exons ( black bars ) from asd . the data for homo sapiens are in left panes ; for drosophila melanogaster in right panels . panels in bottom row show the fit of the alternatively spliced exon length distribution to a single lognormal peak ( full lines ) as well as indicate locations of narrow ( dashed lines ) and wide ( dotted lines ) peaks obtained for the statistical distributions of all exons . inserts in bottom panels show values of correlations coefficients between the statistical distributions of alternatively spliced exons ( ae ) and fitted curve ( fit ) ; narrow peak ( np ) and wide peak ( wp ) . parameters of fitting for different models approximating size distributions of alternatively spliced exons ( see methods sections ) in the case of fitting the model with two lognormal peaks , the fractions of exons contributing to each peak ( table 4 ) were estimated by formal integration of the peak areas , which , for a standard gaussian distribution , are proportional to the peak amplitudes . table 4.distributions of exons between narrow and wide peaks together with peak positions , exp m , and peak widths , exp speciesnarrow peakwide peakexp m ( bp)exp  ( bp)fraction ( % ) exp m ( bp)exp  ( bp)fraction ( % ) caenorhabditis elegans1152.03221854.3278arabidopsis thaliana902.21302335.7670anopheles gambiae1591.73122526.6388tetraodon nigroviridis1282.15531497.0047bos taurus1252.1652997.8248drosophila melanogaster1601.71123268.1888canis familiaris1262.15541018.2046gallus gallus1262.1454969.1346danio rerio1262.15539012.0547macaca mulatta1232.285712915.3443mus musculus1242.316421721.9836homo sapiens1222.336324522.0637pan troglodytes1232.286116122.8439 distributions of exons between narrow and wide peaks together with peak positions , exp m , and peak widths , exp  f - test critical values , f0.05 presented in table s1 in supplementary data , were calculated for significance level of 5% and number of degrees of freedom , df , given in table s1 and table 3 . for each pair of models under consideration when , the hypothesis that the variance i associated with   is greater than variance ii associated with , i.e. the assumption that model ii is a better description of the data than model i ( 11 ) , is verified at the p < 0.05 level . the homo sapiens ( 12 ) and drosophila melanogaster ( 13 ) genomes show ( figure 1 ) a striking similarity in the distributions of exon sizes ; in both , the distribution of the logarithm of exon size forms two distinctive peaks . other genomes [ see supplementary data in which we analyze all complete genome data provided by ensembl ( 5 ) and tair ( 6 ) ] show similar patterns . a simplistic model of intron insertion assumes that the probability of inserting an intron is equal at all positions of an exon ( making it more likely that a longer exon will be split ) . this type of process would lead to an exponential distribution of exon sizes showing a single peak on a logarithmic scale ( 14,15 ) . recently , the simplistic exponential model of exon size distribution was reconsidered by gudlaugsdottir and co - authors ( 15 ) who suggested treating distributions of exon sizes as a combination of weibull ( 16 ) and exponential distributions . however , they did not provide a model for an evolutionary process , which could lead to the weibull distribution but rather consider this distribution to be only an empirical approximation of the observed exon size data . our analysis of thirteen genomes ( table 2 ) , shows that the exponential distribution model always produces the largest -values , which signifies that it is the least adequate fit to the observed exon size distributions . the mixture of weibull and exponential distributions suggested by gudlaugsdottir and co - authors is a better approximation , but for most genomes ( except bos taurus , canis familiaris and gallus gallus ) a combination of two lognormal distributions produces the best agreement between the real data and a fitted model . the length distributions of alternatively spliced exons in homo sapiens and mus musculus , drawn from asd ( 9,10 ) , reveal a single peak pattern ( figure 4 ) . in this case , analysis of -values also shows that an exponential distribution is poorest fit to the data , and the lognormal distribution is the best ( table 3 ) . more sophisticated comparison between different models used in this work involves f - test criterion ( 11 ) and suggests similar hierarchy of the models ( see methods section and table s1 in supplementary data ) . the lognormal distribution is a normal ( gaussian ) distribution of the logarithm of some quantity . this kind of distribution commonly originates from a kolmogoroff random multiplicative process ( 17 ) , which was originally introduced to describe the distribution of ore particle sizes observed in geological samples ( 18,19 ) , and later was found to be useful as a paradigm for a whole universe of different breakage and splitting processes . a modified version of this process , in the context of our problem , can be described as follows . let us consider a single exon , which is split by a random mutational process into two parts ( equal in size , for the sake of illustration ) . then , in the next step of the process , let us assume that one randomly selected part of the ancestral exon undergoes the same splitting . subsequently , this process repeats for a large number of splitting events . the key assumption for this process , which is the same as the assumption of kolmogoroff , is independence of the probability of undergoing a splitting event from exon size . this version of a random multiplicative process is slightly different from that discussed by kolmogoroff , which assumes a constant frequency of splitting events for all parts of the fractionating set ( exons in a genome in our case ) at every time , while the process considered here assumes a single exon splitting event at each step of the process . thus , the probability of breaking a particular exon at the next step of the process is not constant , but constantly decreases with the increase in the number of exons . when started from a single ancestral exon , this model of the multiplicative process , similarly to kolmogoroff 's version , produces a lognormal distribution of exon sizes . the resulting distribution of exon sizes obtained in this manner is independent of the actual sequence of random splitting events . after a sufficiently large number of steps , the mean position of the peak in the distribution of ln ( exon size ) , m , shows an asymptotic linear shift with the logarithm of the number of splitting events , nspl ( for the process introduced here ) or on time , t ( for the kolmogoroff process ) , ln e0  m  ln nspl  t , where e0 is the length of the ancestral exon . similarly , the peak width ,  , is related to the same quantities , . the particular values of the coefficients for these proportionalities depend on the details of the process and can be ignored for our current purposes . however , we assume that the details of the splitting processes , i.e. the coefficients of proportionality , are the same for all species . the assumption that the splitting probability is independent of exon size is essential for a lognormal distribution . a similar hypothesis , that selection of exons from open reading frames is independent of exon size for large exons ( larger than a certain threshold of 105110 bp ) , was discussed earlier ( 20 ) . our model assumes that the probability of exon splitting is independent of exon size for all exons , irrespective of any threshold . we note , in this regard , that the currently most common point of view is that the evolutionary process splits exons by intron insertions . as mentioned in the introduction section , this generally acknowledged hypothesis of exon splitting by transposon insertion implicitly assumes a larger probability of splitting longer exons because they are larger targets for transposons . in contrast to this , the lognormal distribution of exon sizes suggests that the mechanism of intron insertion should be independent of exon size . a kolmogoroff process provides a conceptual background for understanding the lognormal nature of exon size distribution . however , the distribution of exon length in real genomes is somewhat more complicated and generally reveals two lognormal peaks ( figure 1 and supplementary data ) . to demonstrate how this two - peak pattern can be obtained for a random exon splitting process , let us consider a random process , similar to the one described before , initiated for a model ancestral genome that comprises two exons of unequal lengths . in this case , the two parts of the ancestral genome generate two distinct peaks in the exon size distribution . however , unlike previously , the resulting distribution of exon sizes is heavily dependent on the actual splitting pathway . figure 2 demonstrates this fact for two splitting patterns generated from two sequences of pseudorandom numbers ( see methods section ) . this figure reveals a ( spontaneous ) symmetry breaking between the two peaks in the exon size distribution . in a certain sense , this phenomenon is similar to bifurcation , but , in contrast to an instant bifurcation event , this type of behavior can be regarded as a kind of every subsequent splitting event has progressively less effect , and after 5000 steps , the ratio between numbers of exons contributing to the two different peaks remains nearly constant . we believe that this phenomenon of symmetry breaking plays a major role in the diversity of patterns of exon size distributions . the idealized model presented above can be reformulated in a different way where , instead of splitting exons on every step of the process , one randomly selected exon produces a descendant of larger size ( a doublet , for example ) . from a formal point of view , these two variants of the process correspond to two possible directions for the time variable . from a biological prospective , they can be viewed as processes producing orthologous and paralogous gene families : the exons produced by duplications could be considered as paralogs , while exons derived from splitting as orthologs . both variants of the model assume that the process starts from an ancestral genome with few exons , which are later split into smaller parts ( or duplicated to produce longer parts ) . thus , the proposed model , based on kolmogoroff 's ideas , supports the opinion that , at the initial stages of evolution , simple genomes had a lower fraction of introns ( introns late hypothesis ) . figure 3 shows the parameters fit to the exon size distributions for 13 different genomes , using a two component lognormal model . since peak width ,  , is a direct measure of the number of exon splitting events , our model suggests that exons in all of these genomes are distributed between two groups having different evolutionary histories . the width of wider peak varies among species . for simple organisms , like caenorhabditis elegans , its width is close to the width of narrow peak while for higher organisms , such as homo sapiens and pan troglodytes , it is about three times wider than the width of narrow peak . this may suggest both : that exons in this wider peak are older or are evolving more rapidly . we argue that the biological origins of the exons contributing to the narrow and wide peaks are related to the appearance of new splicing mechanisms , such as alternative splicing ( 4 ) . there recently has been great progress in the development of probabilistic methods for determination of exon boundaries ( 2125 ) . however , computational prediction of alternative splicing boundaries is still imperfect . thus , to test the hypothesis that the observed distribution of exons between the two classes correlates with alternative splicing , we used the sets of alternatively spliced exons , which were manually curated and confirmed by est / mrna aligning ( 9,10 ) . unfortunately , curated datasets are available only for human and mouse , but there is a remarkable correlation between these data and the narrow peaks in homo sapiens and mus musculus ( see data in tables 1 and 3 and figure 4 ) . this strongly suggests that alternatively spliced exons are major contributors to the narrow peaks in all of the species examined here . this observation favors the hypothesis that the narrow peak is younger , since alternative splicing is a comparatively advanced evolutionarily mechanism . if one assumes that the evolutionary process modifies all exons in all the species with the same rate , the observation of nearly equal width narrow peaks in all species leads one to conclude that those peaks appeared at approximately the same and quite recent time . at least hypothetically , one may conclude that existence of the narrow peak in all discussed genomes is a manifestation of a spontaneous symmetry break in genome evolution , which is correlated with the evolution of alternative splicing . the assumption that narrow peak is more recent does not completely eliminate the hypothesis that that this peak is more conserved . indeed , one may assume is that exons contributing in narrow peaks are simultaneously more recent , and more conserved ( with respect to the exon size distributions ) than those contributing to the wide peak . in support of this point of view , one may note that the fraction of exons contributing to the narrow peak , in general , correlates with the width of wide peak ( see data in table 4 ) . with two exceptions for genomes of anopheles gambiae and drosophila melanogaster , the increases in the fraction of exons in the narrow peak parallel the increases of width ,  , of the wider peak . in other words , its looks like the exons in the narrow peak are less frequently split by evolutionary processes ( narrow peak width ) , but more frequently duplicated ( larger fraction of exons ) than are exons in the wide peak . this point of view is in agreement with the general pattern of alternative splicing in which a single copy of an exon is replaced by two ( or several ) isoforms of similar lengths . however , we have shown that a very simple process , length independent splitting of exons , can produce the observed lognormal distribution of exon sizes . this suggests that it would be profitable to focus more attention on biological processes that are length independent , or on processes that could constrain the length of exons independently of exon length . processes involved in mrna splicing or associated with intron  exon boundaries are obvious candidates . furthermore , the observation that all eukaryotic organisms possess two exon length classes , one in common among all eukaryotes and one variable , suggests that exon splitting has played a key role in the evolution of eukaryotes . the difference in the peak widths of these two length classes suggests that the wider peak is older , or undergoes more rapid splitting , and mostly comprises nonalternatively spliced exons . the narrow peaks mostly comprise alternatively spliced exons , which are rather conserved in length but multiplied in number in more rapidly evolving genomes . the presence of a narrow peak in all genomes examined could be explained as a manifestation of a spontaneous symmetry break in genome evolution associated with the appearance of the mrna splicing mechanism . more detailed investigations , including characterization of conservation of exons in narrow and wide peaks between different species and between different subclasses within a particular genome , are extremely intriguing but beyond the scope of this article . the main goals of this work are to draw attention to statistical properties of exon size distribution and to highlight the utility of kolmogoroff process model in understanding of genome evolution background .
the quest for evolutionary mechanisms providing separation between the coding ( exons ) and noncoding ( introns ) parts of genomic dna remains an important focus of genetics . 
 this work combines an analysis of the most recent achievements of genomics and fundamental concepts of random processes to provide a novel point of view on genome evolution . 
 exon sizes in sequenced genomes show a lognormal distribution typical of a random kolmogoroff fractioning process . 
 this implies that the process of intron incretion may be independent of exon size , and therefore could be dependent on intron  exon boundaries . 
 all genomes examined have two distinctive classes of exons , each with different evolutionary histories . in the framework proposed in this article 
 , these two classes of exons can be derived from a hypothetical ancestral genome by ( spontaneous ) symmetry breaking . 
 we note that one of these exon classes comprises mostly alternatively spliced exons .