<title>Introduction</title> <p>The Tree of Life is a synoptic depiction of the pathways of evolutionary differentiation between Earth life forms ##UREF##0##[1]##, and contains valuable clues on the key issue of understanding the diversification of life in the planet ##REF##10821281##[2]##. The branching pattern of the Tree of Life, which is being captured at increasing resolution by the advent of molecular tools ##REF##17008515##[3]##, can be examined to investigate fundamental questions, such as whether it follows universal rules, and at what extent random differentiation mechanisms explain the shape of phylogenetic trees. The examination of the structure of the Tree of Life can also help to infer whether evolution acts at intraspecific scales in a way different from the action of evolution at the interspecific scale. Here we address these fundamental questions on the basis of a comprehensive comparative analysis of phylogenetic trees representing different fractions and domains of the Tree of Life, from interspecific to intraspecific scales. We draw from previous analyses of the geometry of the Tree of Life ##REF##16969944##[4]##, the characterization of other branching systems ##UREF##1##[5]##, ##REF##16023678##[6]##, and using tools derived from modern network theory ##REF##12736684##[7]##–##REF##16197119##[10]## to examine the scaling of the branching in the Tree of Life ##UREF##2##[11]##, ##UREF##3##[12]##. Our analysis is based on a thorough data set of more than 5000 interspecific phylogenies and a sample of 67 intraspecific phylogenies (see ##SUPPL##0##Text S1##), thereby testing the universality of the results derived across scales.</p> <p>A phylogenetic tree is a set of nodes, each node representing a diversification event, connected by branches (links). For each node <italic>i</italic>, a subtree <italic>S<sub>i</sub></italic> is made up of a root at node <italic>i</italic> and all the descendant nodes stemming from this root. The subtree size <italic>A<sub>i</sub></italic> gives the number of subtaxa that diversify from node <italic>i</italic> (including itself). Beyond this measure of the diversity degree, the characterization of how the diversity is arranged through the phylogenies can be achieved through the cumulative branch size, <italic>C<sub>i</sub></italic>, a measure of the subtree shape. It is defined ##REF##10335841##[13]## as the sum of the branch sizes associated to all the nodes in the subtree <italic>S<sub>i</sub></italic>, <italic>C<sub>i</sub> = </italic>Σ<italic>A<sub>j</sub></italic>. For the same tree size, and restricting to binary branching events, the smallest value of the cumulative branch size is obtained for a completely symmetric, balanced tree, whereas the most asymmetric, the pectinate or comb-like tree in which all branches split successively from a single one, yields the largest <italic>C<sub>i</sub></italic> value ##REF##10335841##[13]##. To be clearer, we show in ##FIG##0##Figure 1## the analysis of <italic>A<sub>i</sub></italic> and <italic>C<sub>i</sub></italic> for a completely balanced tree (##FIG##0##Figure 1A##) and for a completely imbalanced tree (##FIG##0##Figure 1B##). A portion of a real phylogenetic tree is also shown (##FIG##0##Figure 1C##). How the shape of the tree (i.e., the distribution of the biological diversification) does change with tree size (i.e., with the number of taxa it contains) is given by the scaling of the subtree shape <italic>C</italic> vs. the subtree size <italic>A</italic>, as described by the allometric scaling relation <italic>C∼A<sup>η</sup></italic>. We quantitatively characterize the shape of each tree in our data set by calculating the functions <italic>F</italic>(<italic>A</italic>) and <italic>F</italic>(<italic>C</italic>), which are the complementary cumulative distribution functions (CCDF) of <italic>A<sub>i</sub></italic> and <italic>C<sub>i</sub></italic> values in the tree, respectively, and the value of the allometric scaling exponent, <italic>η</italic>. We compare the results derived from the analyses of inter- and intra-specific phylogenetic trees among them, to test for the preservation of branching patterns across evolutionary scales, and against those derived from the analyses of randomly-generated trees to test whether the allometric scaling derived can be modeled using simple, random branching rules.</p>
<title>Materials and Methods</title> <title>Phylogenies databases</title> <p>On June 30th 2007 we downloaded the 5,212 phylogenetic trees available at that time in the database TreeBASE (<ext-link ext-link-type="uri" xlink:href="http://www.treebase.org">http://www.treebase.org</ext-link>). TreeBASE constitutes a large database of interspecific phylogenies, which were collected from previously published research papers. The size of trees oscillates from 10 to 600 tips. Most of the bifurcations in these trees are binary, as confirmed by the fact that the ratio between the number of tips and the total number of nodes gives 0.52 when averaged over all the trees (for perfect binary trees, the ratio is 0.50).</p> <p>As a comprehensive database comparable to TreeBASE does not exist for intraspecific phylogenies, we constructed an intraspecific data set by manually compiling 67 intraspecific phylogenies from several published phylogenetic analysis [S1–S45]. We compiled this data set in such a way that it contains: 1) Organisms from the main different environments (terrestrial, marine and fresh water), climatic regions (from polar to desert), and branches of life (##SUPPL##1##Table S1##). 2) Phylogenetic trees reconstructed with the main phylogenetic tree estimation methods, i.e., neighbor-joining, maximum parsimony and maximum likelihood methods.</p> <p>In order to test whether the results derived from the examination of the relatively small (67 phylogenies) intraspecific data base can be compared with the much larger (5212) set of interspecific phylogenies extracted from TreeBASE, we sampled the literature to construct a dataset of 67 interspecific phylogenies drawn from the literature [S46–S85] using the same criteria as those to derived the intraspecific phylogeny data base (##SUPPL##1##Table S1##), obtaining full agreement (##SUPPL##2##Figure S1##). The intra- and interspecific phylogenies derived from the literature ranged between 30 and 170 tips, and they contained mainly binary branching events. An example for each kind of phylogenies is shown in ##SUPPL##3##Figures S2A## and ##SUPPL##4##S3A##.</p> <title>Branch size and cumulative branch size distributions</title> <p>We associate to each node <italic>i</italic> of a phylogenetic tree two quantities, the size <italic>A<sub>i</sub></italic> (number of nodes) of the subtree <italic>S<sub>i</sub></italic> made up of node <italic>i</italic> and all the descendant nodes below it, that is, the subtree which does not contain the global root of the original tree, and the cumulative branch size, <italic>C<sub>i</sub></italic>, defined as the sum of the branch sizes associated to all the nodes in the subtree <italic>S<sub>i</sub></italic>, <italic>C<sub>i</sub></italic> = Σ<italic>A<sub>j</sub></italic>. To characterize the probability distributions of the <italic>A<sub>i</sub></italic> and <italic>C<sub>i</sub></italic> values on a particular phylogenetic tree we compute the respective complementary cumulative distribution functions (CCDF): <italic>F</italic>(<italic>A</italic>)<italic> = </italic>probability(<italic>A<sub>i</sub>&gt;A</italic>), and <italic>F</italic>(<italic>C</italic>)<italic> = </italic>probability(<italic>C<sub>i</sub>&gt;C</italic>). We observe that these quantities scale, for large values of <italic>A</italic> and <italic>C</italic>, as power laws:  and . The exponents <italic>τ<sub>A</sub></italic> and <italic>τ<sub>C</sub></italic>, thus, characterize the probabilities of {<italic>A<sub>i</sub></italic>} and {<italic>C<sub>i</sub></italic>}: and, respectively.</p> <title>Allometric scaling relationship</title> <p>We observe that a functional relationship among the values of <italic>C</italic> and <italic>A</italic>, i.e. among shape and size, exists and also follows a power law, <italic>C∼A<sup>η</sup></italic>, characterized by an exponent <italic>η</italic>. Since this relationship encodes the variation of a system property as size is varied, we can call this an <italic>allometric scaling relationship</italic>, to stress its connections with other functional relationships relating function and size ##UREF##2##[11]##, ##REF##10335841##[13]##, ##UREF##14##[27]##. We note that introduction of the change of variables <italic>C∼A<sup>η</sup></italic> into  leads to, from which <italic>η</italic> = (1−<italic>τ<sub>A</sub></italic>)/(1−<italic>τ<sub>C.</sub></italic>). Thus, only two out of the three exponents are independent. As simple examples for which the above exponents can be computed by direct counting, we mention the pectinate or fully unbalanced tree, i.e. a tree in which all branching occurs successively along a single branch, characterized by the exponents <italic>τ<sub>A</sub></italic> = 0, <italic>τ<sub>C</sub></italic> = 1/2, <italic>η</italic> = 2, or the fully symmetric or Cayley tree, characterized by <italic>τ<sub>A</sub></italic> = 2, and <italic>C∼A</italic>ln<italic>A</italic>, which except for the weak logarithmic correction corresponds to <italic>η</italic> = 1 and <italic>τ<sub>C</sub></italic> = 2. ##SUPPL##3##Figures S2B## and ##SUPPL##4##S3B## show, in contrast, the allometric scaling relationship for the particular examples of intra- and inter-specific phylogenies displayed in ##SUPPL##3##Figures S2A## and ##SUPPL##4##S3A##.</p> <p>In order to investigate whether observations differ from random expectations, we have compared the allometric scaling found here with the prediction of a null model ##UREF##16##[29]##, the Equal-rates Markov (ERM) model. The ERM model was attributed to Harding ##UREF##17##[30]##, and to Cavalli-Sforza and Edwards ##REF##28563688##[31]##, although it is based on models of the diversification process that date back at least to Yule ##UREF##10##[23]##. The main assumption of the ERM model is that the phylogeny is the product of random branching. This is the result when the “effective speciation rate” (the difference between extinction and speciation rate) is equal for all species. The effective speciation rate may change chronologically, provided that it is the same for all lineages at a given time ##UREF##10##[23]##. For this model we obtain <italic>C∼A</italic> ln <italic>A</italic>, or <italic>η</italic> = 1, and also <italic>τ<sub>A</sub></italic> = <italic>τ<sub>C</sub></italic> = 2. The random asymmetries introduced by the ERM are not strong enough to change the scaling behavior from the symmetric tree result.</p> <p>The quantity <italic>C<sub>i</sub>/A<sub>i</sub></italic> can be thought as a measure of the average <italic>depth</italic> or <italic>distance</italic> of the phylogenetic tree leaves to the node <italic>i</italic>. This can be seen taking into account that C<italic><sub>i</sub></italic> = Σ(<italic>d<sub>ij</sub></italic>+1), where <italic>d<sub>ij</sub></italic> corresponds to the distance of each of the nodes <italic>j</italic> of the subtree <italic>S<sub>i</sub></italic> to the root <italic>i</italic>. Thus, the relationship between <italic>C</italic> and <italic>A</italic> can be written as <italic>C<sub>i</sub></italic> = <italic>A<sub>i</sub></italic>+〈<italic>d</italic>〉<italic><sub>i</sub>A<sub>i</sub></italic>, where 〈<italic>d</italic>〉<italic><sub>i</sub></italic> is the average depth of the nodes in the subtree <italic>S<sub>i</sub></italic>. The relationship between <italic>C<sub>i</sub>/A<sub>i</sub></italic> and the depth is obtained: C<italic><sub>i</sub></italic>/A<italic><sub>i</sub></italic> = 〈<italic>d</italic>〉<italic><sub>i</sub></italic>+1. This quantity is closely related to the Sackin's index defined as the distance of the leaves to the root: <italic>S</italic> = Σ<italic><sub>l</sub></italic>
<sub>∈<italic>leaves</italic></sub>
<italic>d<sub>l</sub></italic>
<sub>,<italic>root</italic></sub>
##UREF##18##[32]##, ##UREF##19##[33]##. It can be shown that for binary trees <italic>C = 2S+1</italic>, where <italic>C</italic> = Σ<sub>∀<italic>i</italic></sub>
<italic>d<sub>i</sub></italic>
<sub>,<italic>root</italic></sub>. Since the scaling law relating the increase of the depth or Sackin's index with three size is known to be the same as the scaling of the Colless' index, measuring the symmetry or balance of a phylogenetic tree ##UREF##20##[34]##, our results for <italic>η</italic> can be put in the context of the numerous studies available on the unbalance of phylogenetic trees ##REF##16969944##[4]##, ##UREF##7##[17]##, ##REF##16282168##[35]##. Thus, connections between several methodologies previously used to analyze the topology of trees, such as size distributions ##REF##16197119##[10]##, ##UREF##10##[23]##, unbalance and depth ##REF##16969944##[4]##, ##REF##15959469##[8]##, ##UREF##18##[32]##–##REF##16282168##[35]##, and transport efficiency ##REF##12736684##[7]##, ##REF##10335841##[13]##, ##UREF##14##[27]##, ##UREF##15##[28]##, are revealed within the framework presented here.</p>
<title>Results</title> <p>The branch-size CCDF displays power-law tails of the form  for large branch size <italic>A</italic> (##FIG##1##Figure 2A##). The power-law exponents <italic>τ<sub>A</sub></italic> are remarkably similar for the data sets analyzed: <italic>τ<sub>A</sub></italic> = 1.76±0.03, and 1.74±0.02 for intra- and interspecific phylogenies, respectively. Similarly, the cumulative-branch-size CCDF also displays a power-law tail of the form  at large <italic>C</italic>, with a similar agreement between the exponents of the intra- and interspecific data sets: <italic>τ<sub>C</sub></italic> = 1.53±0.02 and 1.53±0.02, respectively (##FIG##1##Figure 2B##). The discrepancy observed between the two data sets at the tail of the distributions can be explained by the different sizes of the typical trees on them: each tree contributes a natural cutoff to the overall distribution, and since the intraspecific trees are smaller in average, their cutoff appears at smaller tree sizes.</p> <p>The allometric exponent, <italic>η</italic>, that characterizes the scaling of tree shape with tree size (##FIG##2##Figure 3A##), is also remarkably similar for the intraspecific (<italic>η</italic> = 1.43±0.01) and the interspecific (<italic>η</italic> = 1.44±0.01) phylogenies. This constancy of the exponents is still more remarkable when realizing (inset of ##FIG##2##Figure 3A##) that it does not only apply to average properties of sets of intraspecific and interspecific trees, but also to individual phylogenies of groups of organisms pertaining to different kingdoms and living across widely contrasting environments, as it is reflected by the very narrow range of <italic>η</italic> obtained from different phylogenies (〈<italic>η</italic>〉 = 1.47, <italic>σ</italic> = 0.03, ##FIG##2##Figure 3A##). The scaling exponents for our large interspecific data set are also matched almost perfectly (##SUPPL##2##Figure S1##) by those derived from a set of 67 interspecific phylogenies randomly drawn from the published literature thereby validating the uniformity of the scaling rules of the broad interspecific phylogenies and the smaller set of intraspecific ones used here. The later was also derived from a similar random sample taken from the published literature (see ##SUPPL##0##Text S1##).</p> <p>The allometric scaling of <italic>C∼A<sup>1.44</sup></italic> derived from our analysis falls somehow in between those obtained by simulated phylogenies derived from two extreme topologies: The symmetric tree gives <italic>C∼A</italic> ln <italic>A</italic>, which corresponds to <italic>η</italic> = 1 with a logarithmic correction, while the pectinate tree has <italic>η</italic> = 2. The natural null model for tree construction, the Equal-Rates Markov (ERM) model ##UREF##4##[14]##, ##UREF##5##[15]##, yields a scaling <italic>C∼A</italic> ln <italic>A</italic> similar to the symmetric tree with <italic>η</italic> = 1 but different from the scaling displayed by empirical inter- and intraspecific phylogenies, particularly for large ones (##FIG##2##Figure 3B##). Therefore some topological aspects of phylogenetic trees are not adequately reproduced by the ERM model. Our results imply that successful lineages diversify more profusely than expected under random branching, generating the large imbalances that characterize emerging depictions of the Tree of Life ##REF##16969944##[4]##. Alternative models introducing correlations, such as the proportional-to-distinguishable-arrangements (PDA) model ##REF##16969944##[4]##, ##UREF##6##[16]## or the beta splitting model ##UREF##7##[17]##, could generate more realistic phylogenies. Guided by previous biological allometric scaling analysis, we have assumed a power-law scaling of the form <italic>C</italic>∼<italic>A<sup>η</sup></italic>. However, other ansatz could also fit the data. The important point, however, is that these modeling approaches should give similar scaling properties for intra- as for interspecific branching.</p>
<title>Discussion</title> <p>Traditionally, microevolutionary and macroevolutionary processes have been studied independently by population geneticists and evolutionary biologists, respectively ##UREF##8##[18]##. The divide between these two levels of generation of biological diversity is an old one, rooted in the controversy between Darwinian gradualism and the saltationism proposed by others, prominently paleontologists, to explain macroevolutionary processes ##REF##28563569##[19]##. The debate as to whether macroevolution is more than the accumulation of microevolutionary events remains active ##UREF##8##[18]##, ##UREF##9##[20]##, ##REF##11258393##[21]##, although refined paleontological evidence supports the continuum between micro- and macroevolution for some lineages ##REF##15241603##[22]##. The results presented here show that the branching and scaling patterns in intraspecific and interspecific phylogenies do not differ significantly for the topological properties we have calculated. Thus, shall saltation processes be a factor at the macroevolutive level, this is not reflected in the topology of phylogenetic branching as examined here. Evidence for possible differences in phylogenetic topologies between the inter- and intraspecific levels may require a detailed analysis of branching times, which we have not attempted.</p> <p>Processes leading to scaling laws in size distributions in natural systems have been formulated as growth models ##UREF##10##[23]##, ##UREF##11##[24]##. Many of the findings carry over to scaling properties found in networks ##UREF##12##[25]## and their description in terms of branching processes ##UREF##13##[26]##. But most of these models predict branching topologies similar to the ERM model. An alternative approach to understand the observed exponent would be to trace analogies with scaling laws in different branching systems ##UREF##1##[5]##, ##REF##16023678##[6]##, ##UREF##14##[27]## which have been explained by invoking a natural optimization criterion based in the fact that the observed trees contain the largest possible number of apices within the smallest number of branching levels. For binary trees of size <italic>A</italic>, where nodes are restricted to occupy uniformly a <italic>D</italic> dimensional Euclidean space, the minimum value of <italic>C</italic> scales as <italic>A<sup>η</sup></italic>, with <italic>η</italic> = <italic>(D+1)/D</italic>. This scaling also describes the <italic>D</italic>-dimensional tree with the maximum size for a given depth (the average distance between root and leaves). The value of <italic>η</italic> obtained in our phylogeny analysis, <italic>η</italic>≅1.44, is achieved only for optimal trees restricted to spaces of <italic>D</italic>≅2.27 dimensions. Given the apparently unlimited number of variables that may yield differences among taxa, restricting their representation to a space with such a small number of dimensions seems unreasonable. This interpretation suggests that the evolutionary process yielding the observed phylogenies is not the most parsimonious one, which could potentially yield a similar biodiversity with fewer branching levels. In fact, the natural choice <italic>D</italic> = ∞ gives an optimal exponent <italic>η</italic> = 1, which correspond to the ERM value and departs from observed scaling. Optimal traffic networks ##UREF##15##[28]## also led to the exponent <italic>τ<sub>A</sub></italic> = 2 which departs from the empirical scaling exponent reported here for phylogenetic trees.</p> <p>In summary, the remarkably similar allometric exponents reported here to characterize universally the scaling properties of intra- and inter-specific phylogenies across kingdoms, reproductive systems and environments, strongly suggests the conservation of branching rules, and hence of the evolutionary processes that drive biological diversification, across the entire history of life. Although at short branch sizes the topology of observed phylogenies cannot differ much from that expected under random and symmetric trees, due to the restriction of binary bifurcations in phylogenetic tree reconstruction, significant departures become universally evident as trees become larger, where the null ERM model and real phylogenies differ (##FIG##1##Figure 2B##). These deviations suggest (a) that the evolution of life leads to less biodiversity than an optimal tree can possibly generate; and (b) the operation of a mechanism generating a correlated branching, where some memory of past evolutionary events is maintained along each branch. This correlated branching pattern implies that entities that diversify faster than average lead to new biological forms that diversify more than average themselves. Invariance across the broad scales considered here indicates that relatively simple rules govern the phylogenetic branching and the unfolding of biodiversity. Their deviation from random models indicates that evolutionary success is a correlated trait within lineages, yielding present asymmetries in the structure of the Tree of Life.</p>
<p>Conceived and designed the experiments: VME. Analyzed the data: EAH CJT KK VME. Contributed reagents/materials/analysis tools: VME. Wrote the paper: EAH VME EHG CMD. Discussed and interpreted the results: CD EH KK CT EH.</p> <p>Understanding the patterns and processes of diversification of life in the planet is a key challenge of science. The Tree of Life represents such diversification processes through the evolutionary relationships among the different taxa, and can be extended down to intra-specific relationships. Here we examine the topological properties of a large set of interspecific and intraspecific phylogenies and show that the branching patterns follow allometric rules conserved across the different levels in the Tree of Life, all significantly departing from those expected from the standard null models. The finding of non-random universal patterns of phylogenetic differentiation suggests that similar evolutionary forces drive diversification across the broad range of scales, from macro-evolutionary to micro-evolutionary processes, shaping the diversity of life on the planet.</p>
<title>Supporting Information</title>
<fig id="pone-0002757-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002757.g001</object-id><label>Figure 1</label><caption><title>Branch size and cumulative branch size examples.</title><p>The values of the branch size (<italic>A</italic>) and of the cumulative branch size (<italic>C</italic>) are shown (in brackets, as (<italic>A</italic>,<italic>C</italic>)) at each node of three small example trees. A: a completely balanced tree of 15 nodes; B: a completely imbalanced tree of 15 nodes; C: a subtree of 15 nodes of a real phylogenetic tree, the intraspecific <italic>Vibrio vulnificus</italic> phylogeny presented in full in ##SUPPL##3##Fig. S2A##. Note that the value of <italic>C</italic> at the root is maximum for the fully imbalanced tree, and minimum for the balanced one.</p></caption></fig> <fig id="pone-0002757-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002757.g002</object-id><label>Figure 2</label><caption><title>Average distributions.</title><p>Cumulative complementary distribution functions (CCDFs) averaged and logarithmically binned over all phylogenetic trees in the interspecific (empty squares) and intraspecific (solid circles) data sets. A: CCDF of branch size, <italic>F</italic>(<italic>A</italic>). Solid line corresponds to a power law  with the exponent given by the best fit to the interspecific data set <italic>τ<sub>A</sub></italic> = 1.74. B: CCDF of the cumulative branch size, <italic>F</italic>(<italic>C</italic>). The line corresponds to a power law with the exponent given by the best fit to the interspecific data set <italic>τ<sub>C</sub></italic> = 1.53.</p></caption></fig> <fig id="pone-0002757-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002757.g003</object-id><label>Figure 3</label><caption><title>Allometric scaling.</title><p>A: Plot of the logarithmically binned set of values of branch size, <italic>A</italic>, and cumulative branch size, <italic>C</italic>, for the interspecific (empty squares) and intraspecific (solid circles) data sets considered. The line corresponds to a power law <italic>C∼A<sup>η</sup></italic>, with the exponent given by the best fit through all data, <italic>η</italic> = 1.44. The inset shows probability distributions of the values of <italic>η</italic> fitted to each individual tree (left: interspecific, right: intraspecific data sets) illustrating the small dispersion in the values. B: Plot of the logarithmically binned set of values of <italic>C</italic> as a function of <italic>A</italic> for the interspecific data, normalized by the prediction from the ERM model (the horizontal line). Data systematically deviate from ERM, especially for large size <italic>A</italic>.</p></caption></fig>
<inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula>
<supplementary-material content-type="local-data" id="pone.0002757.s001"><label>Text S1</label><caption><p>Scaling of branch size and cumulative branch size: TreeBASE vs. manually selected data sets. We provide the list of references corresponding to the selected intraspecific and interspecific phylogenetic trees; the statistics of all data sets with two specific examples; and a summary table of taxa in the data sets.</p><p>(0.06 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002757.s002"><label>Table S1</label><caption><p>Break-down of the number of analyzed inter- and intra-species trees with respect to taxa.</p><p>(0.03 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002757.s003"><label>Figure S1</label><caption><p>Cumulative complementary distribution functions (CCDFs) for branch size (<italic>F</italic>(<italic>A</italic>), panel A) and cumulative branch size (<italic>F</italic>(<italic>C</italic>), panel B), and the allometric scaling relation (<italic>C</italic> {similar, tilde operator } <italic>A</italic>
<sup>η</sup>, panel B) averaged and logarithmically binned over all phylogenetic trees. Empty squares are for the interspecific TreeBASE data set, solid circles are for the manually compiled intraspecific data set, and triangles are for the new manually compiled interspecific data set of reduced size. Solid lines are power laws fitted to the TreeBASE behavior, as in ##FIG##1##Figs. 2## and ##FIG##2##3## of the main text.</p><p>(1.22 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002757.s004"><label>Figure S2</label><caption><p>A: An example of an intraspecific phylogenetic tree: different strains of the bacteria Vibrio vulnificus [S19]. Most of the branchings are binary, but there are some 3rd order branchings. B: The allometric scaling plot showing the relationship of cumulative branch size (<italic>C</italic>) to branch size (<italic>A</italic>) from each node of that tree. The solid line corresponds to the fitting <italic>C</italic> {similar, tilde operator } <italic>A</italic>
<sup>1.43</sup> to this intraspecific dataset.</p><p>(2.66 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002757.s005"><label>Figure S3</label><caption><p>A: An example of an interspecific phylogenetic tree: the catfish species (order Siluriformes) [S80]. Most of the branchings are binary, but there are some 3rd order branchings. B: The allometric scaling plot showing the relationship of cumulative branch size (<italic>C</italic>) to branch size (<italic>A</italic>) from each node of that tree. The solid line corresponds to the fitting <italic>C</italic> {similar, tilde operator } <italic>A</italic>
<sup>1.44</sup> to this intraspecific dataset.</p><p>(2.58 MB TIF)</p></caption></supplementary-material>
<fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>We acknowledge financial support from MEC (Spain) and FEDER, project FISICOS, from CSIC (Spain) project PIE 200750I016, from SBF (Switzerland) through project C05.0148 (Physics of Risk), and from the European Commission through the NEST-Complexity project EDEN. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group>
<graphic xlink:href="pone.0002757.g001"/> <inline-graphic xlink:href="pone.0002757.e001.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002757.e002.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002757.e003.jpg" mimetype="image"/> <graphic xlink:href="pone.0002757.g002"/> <graphic xlink:href="pone.0002757.g003"/> <inline-graphic xlink:href="pone.0002757.e004.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002757.e005.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002757.e006.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002757.e007.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002757.e008.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0002757.e009.jpg" mimetype="image"/>
<media xlink:href="pone.0002757.s001.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002757.s002.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002757.s003.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002757.s004.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002757.s005.tif"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["1"], "element-citation": ["\n"], "surname": ["Cracraft", "Donoghue"], "given-names": ["J", "MJ"], "year": ["2004"], "source": ["Assembling the Tree of Life"], "publisher-loc": ["Oxford"], "publisher-name": ["Oxford University Press"]}, {"label": ["5"], "element-citation": ["\n"], "surname": ["Rodriguez-Iturbe", "Rinaldo"], "given-names": ["I", "A"], "year": ["1997"], "source": ["Fractal river basins: chance and self-organization"], "publisher-loc": ["New York"], "publisher-name": ["Cambridge University Press"]}, {"label": ["11"], "element-citation": ["\n"], "surname": ["LaBarbera"], "given-names": ["M"], "year": ["1989"], "article-title": ["Analyzing Body Size as a Factor in Ecology and Evolution."], "source": ["Annu Rev Ecol Syst"], "volume": ["20"], "fpage": ["97"], "lpage": ["117"]}, {"label": ["12"], "element-citation": ["\n"], "surname": ["Webb", "Brook", "Shine"], "given-names": ["JK", "BW", "R"], "year": ["2002"], "article-title": ["What makes a species vulnerable to extinction? Comparative life-history traits of two sympatric snakes."], "source": ["Ecol Res"], "volume": ["17"], "fpage": ["59"], "lpage": ["67"]}, {"label": ["14"], "element-citation": ["\n"], "surname": ["Mooers", "Heard"], "given-names": ["AO", "SB"], "year": ["1997"], "article-title": ["Inferring evolutionary process from phylogenetic tree shape."], "source": ["Q Rev Biol"], "volume": ["72"], "fpage": ["31"], "lpage": ["54"]}, {"label": ["15"], "element-citation": ["\n"], "surname": ["Caldarelli", "Cartozo", "De Los Rios", "Servedio"], "given-names": ["G", "CC", "P", "VDP"], "year": ["2004"], "article-title": ["Widespread occurrence of the inverse square distribution in social sciences and taxonomy."], "source": ["Phys Rev E"], "volume": ["69"], "fpage": ["035101(1\u20133)"]}, {"label": ["16"], "element-citation": ["\n"], "surname": ["Pinelis"], "given-names": ["I"], "year": ["2003"], "article-title": ["Evolutionary models of phylogenetic trees."], "source": ["Proc R Soc Lond B"], "volume": ["270"], "fpage": ["1425"], "lpage": ["1431"]}, {"label": ["17"], "element-citation": ["\n"], "surname": ["Aldous"], "given-names": ["DJ"], "year": ["2001"], "article-title": ["Stochastic models and descriptive statistics for phylogenetic trees from Yule to today."], "source": ["Stat Sci"], "volume": ["16"], "fpage": ["23"], "lpage": ["34"]}, {"label": ["18"], "element-citation": ["\n"], "surname": ["Simons"], "given-names": ["A"], "year": ["2002"], "article-title": ["The continuity of microevolution and macroevolution."], "source": ["J Evol Biol"], "volume": ["15"], "fpage": ["688"], "lpage": ["701"]}, {"label": ["20"], "element-citation": ["\n"], "surname": ["Grantham"], "given-names": ["T"], "year": ["2007"], "article-title": ["Is macroevolution more than succesive rounds of microevolution?"], "source": ["Paleontology"], "volume": ["50"], "fpage": ["75"], "lpage": ["85"]}, {"label": ["23"], "element-citation": ["\n"], "surname": ["Yule"], "given-names": ["GU"], "year": ["1924"], "article-title": ["A mathematical theory of evolution, based on the conclusions of Dr. J. C. Willis."], "source": ["Philos Trans R Soc Lond A"], "volume": ["213"], "fpage": ["21"], "lpage": ["87"]}, {"label": ["24"], "element-citation": ["\n"], "surname": ["Simon"], "given-names": ["HA"], "year": ["1995"], "article-title": ["On a class of skew distribution functions."], "source": ["Biometrika"], "volume": ["42"], "fpage": ["425"], "lpage": ["440"]}, {"label": ["25"], "element-citation": ["\n"], "surname": ["Bornholdt", "Ebel"], "given-names": ["S", "H"], "year": ["2001"], "article-title": ["World Wide Web Scaling Exponent from Simon's 1955 Model."], "source": ["Phys Rev E"], "volume": ["64"], "fpage": ["035104(R)"]}, {"label": ["26"], "element-citation": ["\n"], "surname": ["Durrett"], "given-names": ["R"], "year": ["2007"], "source": ["Random Graph Dynamics"], "publisher-loc": ["Cambridge"], "publisher-name": ["Cambridge University Press"]}, {"label": ["27"], "element-citation": ["\n"], "surname": ["Brown", "Gillooly", "Allen", "Savage", "West"], "given-names": ["JH", "JF", "AP", "VM", "GB"], "year": ["2004"], "article-title": ["Toward a Metabolic Theory of Ecology."], "source": ["Ecology"], "volume": ["85"], "fpage": ["1771"], "lpage": ["1789"]}, {"label": ["28"], "element-citation": ["\n"], "surname": ["Barth\u00e9lemy", "Flammini"], "given-names": ["M", "A"], "year": ["2006"], "article-title": ["Optimal traffic networks."], "source": ["J Stat Mech"], "volume": ["07"], "fpage": ["L07002"]}, {"label": ["29"], "element-citation": ["\n"], "surname": ["Harvey", "Colwell", "Silvertown", "May"], "given-names": ["PH", "RK", "JW", "RM"], "year": ["1983"], "article-title": ["Null models in ecology."], "source": ["Ann Rev Ecol Syst"], "volume": ["14"], "fpage": ["189"], "lpage": ["211"]}, {"label": ["30"], "element-citation": ["\n"], "surname": ["Harding"], "given-names": ["EF"], "year": ["1971"], "article-title": ["The probabilities of rooted tree-shapes generated by random bifurcation."], "source": ["Adv Appl Prob"], "volume": ["3"], "fpage": ["44"], "lpage": ["77"]}, {"label": ["32"], "element-citation": ["\n"], "surname": ["Sackin"], "given-names": ["MJ"], "year": ["1972"], "article-title": ["\u201cGood\u201d and \u201cbad\u201d phenograms."], "source": ["Sys Zool"], "volume": ["21"], "fpage": ["225"], "lpage": ["226"]}, {"label": ["33"], "element-citation": ["\n"], "surname": ["Shao", "Sokal"], "given-names": ["KT", "R"], "year": ["1990"], "article-title": ["Tree balance."], "source": ["Sys Zool"], "volume": ["39"], "fpage": ["226"], "lpage": ["276"]}, {"label": ["34"], "element-citation": ["\n"], "surname": ["Ford"], "given-names": ["DJ"], "year": ["2006"], "article-title": ["Probabilities on cladograms: introduction to the alpha model (PhD Thesis, Stanford University)"]}]
PLoS One. 2008 Jul 23; 3(7):e2757