Testing substitution models within a phylogenetic tree

Phylogenetic tree reconstruction frequently assumes the homogeneity of the substitution process over the whole tree. To test this assumption statistically, we propose a test based on the sample covariance matrix of the set of substitution rate matrices estimated from pairwise sequence comparison. The sample covariance matrix is condensed into a one-dimensional test statistic Delta = sum ln(1 + delta(i)), where delta(i) are the eigenvalues of the sample covariance matrix. The test does not assume a specific mutational model. It analyses the variation in the estimated rate matrices. The distribution of this test statistic is determined by simulations based on the phylogeny estimated from the data. We study the power of the test under various scenarios and apply the test to X chromosome and mtDNA primate sequence data. Finally, we demonstrate how to include rate variation in the test.     

Constraining prior probabilities of phylogenetic trees

Although Bayesian methods are widely used in phylogenetic systematics today, the foundations of this methodology are still debated among both biologists and philosophers. The Bayesian approach to phylogenetic inference requires the assignment of prior probabilities to phylogenetic trees. As in other applications of Bayesian epistemology, the question of whether there is an objective way to assign these prior probabilities is a contested issue. This paper discusses the strategy of constraining the prior probabilities of phylogenetic trees by means of the Principal Principle. In particular, I discuss a proposal due to Velasco (Biol Philos 23:455–473, 2008) of assigning prior probabilities to tree topologies based on the Yule process. By invoking the Principal Principle I argue that prior probabilities of tree topologies should rather be assigned a weighted mixture of probability distributions based on Pinelis’ (P Roy Soc Lond B Bio 270:1425–1431, 2003) multi-rate branching process including both the Yule distribution and the uniform distribution. However, I argue that this solves the problem of the priors of phylogenetic trees only in a weak form.     

Cyanobacteria and the Great Oxidation Event: evidence from genes and fossils

Cyanobacteria are among the most ancient of evolutionary lineages, oxygenic photosynthesizers that may have originated before 3.0 Ga, as evidenced by free oxygen levels. Throughout the Precambrian, cyanobacteria were one of the most important drivers of biological innovations, strongly impacting early Earth's environments. At the end of the Archean Eon, they were responsible for the rapid oxygenation of Earth's atmosphere during an episode referred to as the Great Oxidation Event (GOE). However, little is known about the origin and diversity of early cyanobacterial taxa, due to: (1) the scarceness of Precambrian fossil deposits; (2) limited characteristics for the identification of taxa; and (3) the poor preservation of ancient microfossils. Previous studies based on 16S rRNA have suggested that the origin of multicellularity within cyanobacteria might have been associated with the GOE. However, single‐gene analyses have limitations, particularly for deep branches. We reconstructed the evolutionary history of cyanobacteria using genome scale data and re‐evaluated the Precambrian fossil record to get more precise calibrations for a relaxed clock analysis. For the phylogenomic reconstructions, we identified 756 conserved gene sequences in 65 cyanobacterial taxa, of which eight genomes have been sequenced in this study. Character state reconstructions based on maximum likelihood and Bayesian phylogenetic inference confirm previous findings, of an ancient multicellular cyanobacterial lineage ancestral to the majority of modern cyanobacteria. Relaxed clock analyses provide firm support for an origin of cyanobacteria in the Archean and a transition to multicellularity before the GOE. It is likely that multicellularity had a greater impact on cyanobacterial fitness and thus abundance, than previously assumed. Multicellularity, as a major evolutionary innovation, forming a novel unit for selection to act upon, may have served to overcome evolutionary constraints and enabled diversification of the variety of morphotypes seen in cyanobacteria today.     

Phylogenomic Evidence for the Origin of Obligate Anaerobic Anammox Bacteria Around the Great Oxidation Event

The anaerobic ammonium oxidation (anammox) bacteria can transform ammonium and nitrite to dinitrogen gas, and this obligate anaerobic process accounts for up to half of the global nitrogen loss in surface environments. Yet its origin and evolution, which may give important insights into the biogeochemistry of early Earth, remain enigmatic. Here, we performed a comprehensive phylogenomic and molecular clock analysis of anammox bacteria within the phylum Planctomycetes. After accommodating the uncertainties and factors influencing time estimates, which include implementing both a traditional cyanobacteria-based and a recently developed mitochondria-based molecular dating approach, we estimated a consistent origin of anammox bacteria at early Proterozoic and most likely around the so-called Great Oxidation Event (GOE; 2.32–2.5 Ga) which fundamentally changed global biogeochemical cycles. We further showed that during the origin of anammox bacteria, genes involved in oxidative stress adaptation, bioenergetics, and anammox granules formation were recruited, which might have contributed to their survival on an increasingly oxic Earth. Our findings suggest the rising levels of atmospheric oxygen, which made nitrite increasingly available, was a potential driving force for the emergence of anammox bacteria. This is one of the first studies that link the GOE to the evolution of obligate anaerobic bacteria.     

Inferring the root of a phylogenetic tree

Phylogenetic trees can be rooted by a number of criteria. Here, we introduce a Bayesian method for inferring the root of a phylogenetic tree by using one of several criteria: the outgroup, molecular clock, and nonreversible model of DNA substitution. We perform simulation analyses to examine the relative ability of these three criteria to correctly identify the root of the tree. The outgroup and molecular clock criteria were best able to identify the root of the tree, whereas the nonreversible model was able to identify the root only when the substitution process was highly nonreversible. We also examined the performance of the criteria for a tree of four species for which the topology and root position are well supported. Results of the analyses of these data are consistent with the simulation results.     

Characterization of a branch of the phylogenetic tree

We use a combination of analytic models and computer simulations to gain insight into the dynamics of evolution. Our results suggest that certain interesting phenomena should eventually emerge from the fossil record. For example, there should be a "tortoise and hare effect": those genera with the smallest species death rate are likely to survive much longer than genera with large species birth and death rates. A complete characterization of the behavior of a branch of the phylogenetic tree corresponding to a genus and accurate mathematical representations of the various stages are obtained. We apply our results to address certain controversial issues that have arisen in paleontology such as the importance of punctuated equilibrium and whether unique Cambrian phyla have survived to the present.     

Shaping the phylogenetic tree of influenza by cross-immunity

Cross-immunity among related strains can account for the selection producing the slender phylogenetic tree of influenza A and B in humans. Using a model of seasonal influenza epidemics with drift (Andreasen, 2003. Dynamics of annual influenza A epidemics with immuno-selection. J. Math. Biol. 46, 504-536), and assuming that two mutants arrive in the host population sequentially, we determine the threshold condition for the establishment of the second mutant in the presence of partial cross-protection caused by the first mutant and their common ancestors. For fixed levels of cross-protection, the chance that the second mutant establishes increases with rho the basic reproduction ratio and some temporary immunity may be necessary to explain the slenderness of flu's phylogenetic tree. In the presence of moderate levels of temporary immunity, an asymmetric situation can arise in the season after the two mutants were introduced and established: if the offspring of the new mutant arrives before the offspring of the resident type, then the mutant-line may produce a massive epidemic suppressing the original lineage. However, if the original lineage arrives first then both strains may establish and the phylogenetic tree may bifurcate.     

Rubisco genes indicate a close phylogenetic relation between the plastids of Chromophyta and Rhodophyta

The genes for both subunits of Rubisco (rbcL, rbcS) are located on the plastome of the brown alga Ectocarpus siliculosus (Chromophyta, Phaeophyceae). The organization of these genes in the form of an operon was similar to that found in rhodoplasts, cyanobacteria and the plastids of Cryptomonas . Sequence analysis of the complete operon revealed a high degree of homology and great structural similarities to corresponding genes from two red algae. In contrast, sequence homology to Rubisco genes from chloroplasts and cyanobacteria was much lower. This clearly indicated a close phylogenetic relationship between the plastids of Rhodophyta and Chromophyta which seem to have evolved independently from the chloroplasts (polyphyletic origin). Our data suggest that the plastids of Chromophyta and Cryptophyta have originated from endosymbiotic unicellular red algae. Surprisingly, red and brown algal Rubiscos show a significantly higher degree of homology to that from a hydrogen bacterium than to those from cyanobacteria.     

The probability of a gene tree topology within a phylogenetic network with applications to hybridization detection

Gene tree topologies have proven a powerful data source for various tasks, including species tree inference and species delimitation. Consequently, methods for computing probabilities of gene trees within species trees have been developed and widely used in probabilistic inference frameworks. All these methods assume an underlying multispecies coalescent model. However, when reticulate evolutionary events such as hybridization occur, these methods are inadequate, as they do not account for such events. Methods that account for both hybridization and deep coalescence in computing the probability of a gene tree topology currently exist for very limited cases. However, no such methods exist for general cases, owing primarily to the fact that it is currently unknown how to compute the probability of a gene tree topology within the branches of a phylogenetic network. Here we present a novel method for computing the probability of gene tree topologies on phylogenetic networks and demonstrate its application to the inference of hybridization in the presence of incomplete lineage sorting. We reanalyze a Saccharomyces species data set for which multiple analyses had converged on a species tree candidate. Using our method, though, we show that an evolutionary hypothesis involving hybridization in this group has better support than one of strict divergence. A similar reanalysis on a group of three Drosophila species shows that the data is consistent with hybridization. Further, using extensive simulation studies, we demonstrate the power of gene tree topologies at obtaining accurate estimates of branch lengths and hybridization probabilities of a given phylogenetic network. Finally, we discuss identifiability issues with detecting hybridization, particularly in cases that involve extinction or incomplete sampling of taxa.     

On the repeat-annotated phylogenetic tree reconstruction problem.

A new problem in phylogenetic inference is presented, based on recent biological findings indicating a strong association between reversals (i.e., inversions) and repeats. These biological findings are formalized here in a new mathematical model, called repeat-annotated phylogenetic trees (RAPT). We show that, under RAPT, the evolutionary process--including both the tree-topology as well as internal node genome orders--is uniquely determined, a property that is of major significance both in theory and in practice. Furthermore, the repeats are employed to provide linear-time algorithms for reconstructing both the genomic orders and the phylogeny, which are NP-hard problems under the classical model of sorting by reversals (SBR).     

分子系统发育树构建的简易方法

以系统发育树构建的原有距离方法为基础,吸取了NJ法和FM法中的部分理论,提出了以节点引入为手段的新的简易方法,通过该方法构建了分子系统发育树,结果表明这种方法更加快捷,而且所得结果与FM法完全一致。     

The root of the phylogenetic tree of human populations

Although African populations have been shown to be most divergent from any other human populations, it has been difficult to establish the root of the phylogenetic tree of human populations since the rate of evolutionary change may vary from population to population owing to the fluctuation of population size and other factors. However, the root can be determined by using the chimpanzee as an outgroup and by employing proper statistical methods. Using this strategy, we constructed phylogenetic trees of human populations for five different sets of gene frequency data. The data sets used were two sets of microsatellite loci data (25 and 8 loci, respectively), restriction fragment length polymorphism (RFLP) data (79 loci), protein polymorphism data (15 loci), and Alu insertion frequency data (4 loci). All these data sets showed that the root is located in the branch connecting African and non-African populations, and in the four data sets the root was established at a significant level. These results indicate that Africans are the first group of people that split from the rest of the human populations.      

树鼩进化分类地位的分子证据

树鼩隶属攀鼩目,广泛分布于东南亚、南亚和中国南部等地区。由于其独特的特点,如体型小、脑-体重比例高、生殖周期短、寿命短和饲养成本低等,在生物医学研究中被认为是可望替代灵长类动物的新型实验动物。然而,关于树鼩与灵长类动物的亲缘关系一直存在争议。明确树鼩的分类地位是创建实验动物的重要研究基础。该文介绍了近年来关于树鼩分类地位探讨的分子证据。在现有的研究中,大部分核DNA序列研究,包括近期树鼩全基因组序列分析,都支持树鼩是灵长动物的近缘旁系群,然而绝大部分基于线粒体DNA序列的研究却显示树鼩与啮齿动物的亲缘关系更为接近。这样的分歧主要是由于线粒体序列和核基因数据的差异以及不同的算法导致。综合现有不同DNA数据的研究结果,作者认为树鼩作为灵长类的近亲这一结论应该成为共识。     

The phylogenetic relationship within the genus Carcharhinus

The phylogenetic relationships of elasmobranchs, especially sharks, are unclear. All previously made research and hypotheses indicate that there are still unsolved relationships within and between the class Chondrichthyes. To find out the relationship and sister group of this genus, the ribosomal ITS1-2 regions sequence has been chosen to differentiate the genus Carcharhinus from the blue shark (genus Prionace) and from some other genus species as an outgroup. The results show that the blue shark is a member of the genus Carcharhinus, suggesting that maybe the blue shark belongs also to the genus Carcharhinus instead of Prionace, or that there is a misclassification, Prionace being not a separate genus.     

Conflicting phylogenetic signals at the base of the metazoan tree

A phylogenetic framework is essential for under-standing the origin and evolution of metazoan development. Despite a number of recent molecular studies and a rich fossil record of sponges and cnidarians, the evolutionary relationships of the early branching metazoan groups to each other and to a putative outgroup, the choanoflagellates, remain uncertain. This situation may be the result of the limited amount of phylogenetic information found in single genes and the small number of relevant taxa surveyed. To alleviate the effect of these analytical factors in the phylogenetic recons-truction of early branching metazoan lineages, we cloned multiple protein-coding genes from two choanoflagellates and diverse sponges, cnidarians, and a ctenophore. Comparisons of sequences for alpha-tubulin, beta-tubulin, elongation factor 2, HSP90, and HSP70 robustly support the hypothesis that choanoflagellates are closely affiliated with animals. However, analyses of single and concatenated amino acid sequences fail to resolve the relationships either between early branching metazoan groups or between Metazoa and choano-flagellates. We demonstrate that variable rates of evolution among lineages, sensitivity of the analyses to taxon selection, and conflicts in the phylogenetic signal contained in different amino acid sequences obscure the phylogenetic associations among the early branching Metazoa. These factors raise concerns about the ability to resolve the phylogenetic history of animals with molecular sequences. A consensus view of animal evolution may require investigations of genome-scale characters.     

Drought timing influences the legacy of tree growth recovery

Whether and how the timing of extreme events affects the direction and magnitude of legacy effects on tree growth is poorly understood. In this study, we use a global database of Ring‐Width Index (RWI) from 2,500 sites to examine the impact and legacy effects (the departure of observed RWI from expected RWI) of extreme drought events during 1948–2008, with a particular focus on the influence of drought timing. We assessed the recovery of stem radial growth in the years following severe drought events with separate groupings designed to characterize the timing of the drought. We found that legacies from extreme droughts during the dry season (DS droughts) lasted longer and had larger impacts in each of the 3 years post drought than those from extreme droughts during the wet season (WS droughts). At the global scale, the average integrated legacy from DS droughts (0.18) was about nine times that from WS droughts (0.02). Site‐level comparisons also suggest stronger negative impacts or weaker positive impacts of DS droughts on tree growth than WS droughts. Our results, therefore, highlight that the timing of drought is a crucial factor determining drought impacts on tree recovery. Further increases in baseline aridity could therefore exacerbate the impact of punctuated droughts on terrestrial ecosystems.     

Sampling phylogenetic tree space with the generalized Gibbs sampler

A recent article published in Cladistics is critical of a number of heuristic methods for phylogenetic inference based on parsimony scores. One of my papers is among those criticized, and I would appreciate the opportunity to make a public response. The specific criticism is that I have re‐invented an algorithm for economizing parsimony calculations on trees that differ by a subtree pruning and regrafting (SPR) rearrangement. This criticism is justified, and I apologize for incorrectly claiming originality for my presentation of this algorithm. However, I would like to clarify the intent of my paper, if I can do so without detracting from the sincerity of my apology. My paper is not about that algorithm, nor even primarily about parsimony. Rather, it is about a novel strategy for Markov chain Monte Carlo (MCMC) sampling in a state space consisting of trees. The sampler involves drawing from conditional distributions over sets of trees: a Gibbs‐like strategy that had not previously been used to sample tree‐space. I would like to see this technique incorporated into MCMC samplers for phylogenetics, as it may have advantages over commonly used Metropolis‐like strategies. I have recently used it to sample phylogenies of a biological invasion, and I am finding many applications for it in agent‐based Bayesian ecological modelling. It is thus my contention that my 2005 paper retains substantial value.