Scaling properties of protein family phylogenies
© Herrada et al; licensee BioMed Central Ltd. 2011
Received: 14 March 2011
Accepted: 6 June 2011
Published: 6 June 2011
Skip to main content
© Herrada et al; licensee BioMed Central Ltd. 2011
Received: 14 March 2011
Accepted: 6 June 2011
Published: 6 June 2011
One of the classical questions in evolutionary biology is how evolutionary processes are coupled at the gene and species level. With this motivation, we compare the topological properties (mainly the depth scaling, as a characterization of balance) of a large set of protein phylogenies with those of a set of species phylogenies.
The comparative analysis between protein and species phylogenies shows that both sets of phylogenies share a remarkably similar scaling behavior, suggesting the universality of branching rules and of the evolutionary processes that drive biological diversification from gene to species level. In order to explain such generality, we propose a simple model which allows us to estimate the proportion of evolvability/robustness needed to approximate the scaling behavior observed in the phylogenies, highlighting the relevance of the robustness of a biological system (species or protein) in the scaling properties of the phylogenetic trees.
The invariance of the scaling properties at levels spanning from genes to species suggests that rules that govern the incapability of a biological system to diversify are equally relevant both at the gene and at the species level.
During the last century, an important effort has been devoted to the understanding of diversification patterns and processes in terms of branching evolutionary trees [1–7]. Tempo and mode of genetic change, and their connections with tempo and mode of speciation is an important issue in this context. In that sense, we address the question of whether similar forces act across the gene level and species-level evolution [8–10], through a comparative analysis of the topological behavior of protein and species phylogenies.
Previous analyses of the topological properties of phylogenies have revealed universal patterns of phylogenetic differentiation [3, 6, 7, 11, 12]. This means that the impact of evolutionary forces shaping the diversity of life on Earth on the shape of phylogenetic trees is, at least to the level of detail captured by the descriptors used, similar across a broad range of scales, from macro-evolution to speciation and population differentiation, and across diverse organisms such as eukaryotes, eubacteria, archaea or viruses, thereby. This together with the fact that evolutionary forces work at molecular level motivates the study of the topology of evolutionary relationships among molecular entities, looking for patterns of differentiation at such molecular level, thereby extending the examination of the universality of the scaling of branching laws in phylogenies all the way from molecular- to macro-evolution.
A protein family phylogeny is represented as a tree, i.e., as an acyclic graph of nodes connected by branches (links), where each node represents a diversification event. For each node in a phylogeny, a subtree (or subfamily) S is made of the root at the selected node and all of its descendant nodes. The subtree size A is the total number of subfamily members that diversify from the root (including itself). The characterization of how protein diversity is arranged through the phylogenies can be achieved in a variety of ways [19–26]. We focus here on the mean depth, d, of the subtree S (see Methods) [6, 27, 28] defined as: , where, for a given node j, d root,j is its topological distance to the root of the subtree S, that is, the number of nodes one has to go through so as to go from that node to the root (including the root in the counting), and the sum is over all nodes in the subtree S. Note that we use here the mean depth over all subtree nodes and not just the leaves, which gives a different but related measure [4, 29, 30]. In the remainder, when no subindex is indicated, we understand that mean depth and other quantities refer to a whole tree or a subtree depending on the context.
How the shape of a phylogenetic tree, i.e., the distribution of protein diversification, changes with tree size, i.e., with the number of proteins it contains, can be analyzed by examining the dependence of the mean depth on subfamily size d = d(A). This gives information on the balance characteristics of the tree. To be clearer, in the additional file 1 we show the analysis of A and d for a fully balanced and a fully imbalanced 15-tip tree, as well as for a 15-tip subtree of a real phylogenetic tree. For a given tree size, the smallest value of the mean depth corresponds to the fully polytomic tree. The mean depth d as a function of tree size A is given in this case by
The leading contribution at large sizes is logarithmic: d ~ ln A. This logarithmic scaling is not exclusive of fully balanced trees, it is also the behavior of the Equal-Rates Markov (ERM) model [28, 31, 32], the natural null model for stochastic tree construction, in which, at each time step, one of the existing leaves of the tree is chosen at random and bifurcated into two new leaves.
We report here the patterns of mean depth for protein families, and compare the branching patterns derived for protein families, from the PANDIT database with those of species phylogenies, reported previously from the TreeBASE database . This comparison shows that branching patterns are mostly preserved across evolutionary scales spanning from genes to species.
The depth scaling behavior shared by protein and species phylogenies can be explained by different branching mechanisms. In this direction, during the last decade, several models have been published proposing different mechanisms to capture the topology of phylogenetic trees [6, 27, 28, 33, 35, 36]. Most of the models proposed yield a logarithmic scaling of the mean depth, i.e., ERM-type for large sizes [31, 32, 37], which is not a good description of our data (see Figure 2 and additional file 2), at least at the tree sizes available; the AB model proposed in Ref.  is one of the few models that deviate from the ERM-like scaling leading to a squared logarithmic d ~ (ln A)2 (see also ); models with power law scaling of the mean depth d ~ A η have also been defined in terms of statistical rules assigning probabilities to different splittings or types of trees  or in terms of (simplified) evolutionary events (in the sense specified in Ref. ) occurring in time [27, 28].
An alternative explanation of the scaling properties of the phylogenetic trees  suggests that the non-ERM behavior is a small-size transient behavior, which would cross-over to the ERM scaling d ~ ln A as larger tree sizes become available.
The process conducive to trees that deviate from ERM behavior is the presence of temporal correlations, which leads to asymptotic or just finite-size deviations with respect to the ERM behavior depending on whether these correlations are permanent or restricted to finite but large times. We, thus, explored the role of such correlations through a simple model based on the inheritability of the evolvability, i.e., the ability to evolve [38, 39], as a biological characteristic which is itself inherited by sister species in speciation events. The process starts with the root, which we consider capable to speciate. At each time step, all present species capable to speciate branch simultaneously. Each branching event yields two new daughter species, for which we allow two possible outcomes:
with probability p, the new species inherit the evolvability of the mother species, i.e., they have the same capacity as the mother species to speciate again;
with probability 1- p, one of the daughter species is unable to speciate again, that is, only one of the two daughter species preserves the ability to evolve. Stemming from the definition of robustness as the property of a system to remain invariant in the presence of genetic or environmental perturbations , we consider a species' inability to speciate its robustness.
The trees generated with this algorithm yield a scaling very close to those observed for phylogenetic trees in both PANDIT and TreeBASE for p = 0.24 (see Figure 6, and additional file 3). This result identifies the prevalence of imbalanced branching events (occurring with frequency 1- p = 0.76) relative to balanced ones (p = 0.24), which is consistent with earlier reports [5, 6, 33].
At large n, the leading contributions are A ~ z n and C ~ nz n (we do not write explicitly prefactors which may depend on z but not on n). Taking into account Eq. (7) in Methods (i.e. d = (C/A) - 1) and inverting the relationship between A and n (n ~ ln A), we obtain that for large sizes the leading order of the mean depth is d ~ ln A, which indicates that what we observe in the simulations is a long transient behavior. This transient behavior leads to the fact that our model ts the proper behavior of the data at the sizes in the databases, but the asymptotic scaling at the larger sizes will finally be d ~ ln A, as in the ERM.
The development of high-throughput "-omics" has provided the data required to address the traditional debate on how gene-level evolution shapes the species-level evolution [8–10]. This debate connects with that on the (dis)continuity between micro- and macro-evolution, and gradualism versus saltationism [41–43]. In the context of these debates, the universal scaling of phylogenetic trees at intra and inter-specific levels shown earlier  suggested the conservation of the evolutionary processes that drive biological diversification across the entire history of life. Here we extend this observation further to demonstrate that the universality of the scaling properties can also be extrapolated to the gene-level. The results presented here show that the branching and scaling patterns in protein families do not differ significantly from the patterns observed in species phylogenies, at least for the topological properties we have calculated. We do not observe any discrepancy between the shape of protein phylogenies and species phylogenies. Moreover, the results presented here shows no evidence for possible differences in phylogenetic trees among protein families with different biological functions, further providing evidence of universal, conserved evolutionary processes from genes to species.
In 2006, Cotton and Page published a comparative analysis between human gene phylogenies and species phylogenies . They found quantitative differences between human paralogous gene and orthologous gene phylogenies. Their research focused on the comparison between (small) paralogous and orthologous gene families, while here we have analyzed complete protein families, which included both paralogous and orthologous protein members, focusing on the comparison between protein and species phylogenies. Our approach is based on a scaling analysis, examining how variables change with tree size, whereas the Cotton-Page's approach is based on a quantitative analysis of small sizes. This implies that despite their finding of quantitative differences between paralogous and orthologous gene phylogenies, we expect that both phylogenies would display scaling behavior similar to that we described here for complete protein phylogenies and organism phylogenies .
Different evolutionary models and mechanisms have been proposed to explain the branching patterns arising in evolution [6, 27, 28, 33, 36, 37]. Here we have introduced a simple model accounting for differences in the degree of evolvability, which is emerging as a key trait constraint as important as robustness in evolution [44–47]. The model we proposed can be interpreted in the framework of the balance between evolvability as the potential of a biological system for future adaptive mutation and evolution , and robustness as the property of a system to produce relatively invariant output in the presence of a perturbation . Indeed, the symmetric diversification event should correspond to the biological context in which the biological system is evolvable, while the asymmetric diversification process should correspond to a biological context where the new biological system, which has just appeared from the diversification process, is robust and unable of unlimited diversification.
The asymptotic behavior of our model at long tree sizes recovers the logarithmic behavior of the ERM scaling, so that, as in the models by , the non-ERM behavior occurs as a transient for the relatively small tree sizes present in the databases. Despite this, the local (i.e. present for finite sizes) imbalance in real trees can be interpreted in terms of the evolvability concept. The prevalence of the unbalanced branching found is consistent with previous works [6, 33, 48–51], and has been traditionally explained by the presence of variations in the speciation and/or extinction rates throughout the Tree of Life [4, 5].
Different biological explanations for these variations in the speciation and/or extinction rates have been proposed, such as: refractory period , mass extinctions , specialization  or environment effects . The consideration of an evolutionary scenario based on the evolvability/robustness interplay has led us to postulate the presence of asymmetric diversification events over the depth scaling during evolutionary processes giving rise to a new biological system which is unable to undergo a new diversification event. An incapability to diversify may occur at different levels of evolution, and can be found at the macroevolutionary level with taxa that require very long refractory periods or with random massive extinctions of taxa, as well as at the microevolutionary or gene level, where the elements unable to diversify are individuals from a population or genetic variants from a cell, embryo or individual.
In summary, the finding of universal scaling properties at gene and species level, characterized by the similar scaling laws, strongly suggest the universality of branching rules, and hence of the evolutionary processes that drive biological diversification across the entire history of life, from genes to species. The topological characterization of phylogenetic trees has proven helpful to analyze the relevance of the robustness of a biological system (species or protein) in the scaling properties of the phylogenetic trees. Thus, the invariance of the scaling properties at levels spanning from genes to species suggests that the mechanisms leading to the incapability of a biological system to diversify for a very long period of time act at both the gene- and species-level.
We analyzed the 7,738 protein families available in the PANDIT database (http://www.ebi.ac.uk/goldman-srv/pandit/ accession date May 27, 2008) . PANDIT is based upon Pfam http://pfam.sanger.ac.uk/, and constitutes a large collection of protein family phylogenies from different signalling pathways, cellular organelles and biological functions, reconstructed with five different methods: NJ , BioNJ , Weighbor , FastME  and Phyml . The size of each of the protein phylogenies, T, ranges from 2 to more than 2000 tips (i.e. proteins within families) and, in agreement with previous reports [20, 61–64], shows a power law distribution P(T) ~T -γ (see Figure 1). Most of the bifurcations in these phylogenies are binary, with only 22% of polytomic bifurcations.
Accounting for the initial condition, that is, the root, with C = 1 and S = 0, yields C = 2S + 1 for binary trees. Thus, at large sizes, both quantities, C and S, become proportional and scale in the same way with size.
We acknowledge financial support from the European Commission through the NEST-Complexity project EDEN (043251) and from MICINN (Spain) and FEDER through project FISICOS (FIS2007-60327).
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.