scholarly journals PRANC: ML species tree estimation from the ranked gene trees under coalescence

2020 ◽  
Vol 36 (18) ◽  
pp. 4819-4821
Author(s):  
Anastasiia Kim ◽  
James H Degnan
Keyword(s):  
Maximum Likelihood ◽  
Gene Tree ◽  
Species Tree ◽  
Supplementary Information ◽  
Supplementary Data ◽  
Gene Trees ◽  
Multispecies Coalescent ◽  
Branch Lengths ◽  
Tree Topologies ◽  
Tree Estimation

Abstract Summary PRANC computes the Probabilities of RANked gene tree topologies under the multispecies coalescent. A ranked gene tree is a gene tree accounting for the temporal ordering of internal nodes. PRANC can also estimate the maximum likelihood (ML) species tree from a sample of ranked or unranked gene tree topologies. It estimates the ML tree with estimated branch lengths in coalescent units. Availability and implementation PRANC is written in C++ and freely available at github.com/anastasiiakim/PRANC. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals Complexity of the simplest species tree problem

2021 ◽  
Author(s):  
Tianqi Zhu ◽  
Ziheng Yang
Keyword(s):  
Maximum Likelihood ◽  
Estimation Error ◽  
Gene Tree ◽  
Species Tree ◽  
Ancestral Population ◽  
Statistical Efficiency ◽  
Multispecies Coalescent ◽  
Full Likelihood ◽  
Tree Methods ◽  
Tree Estimation

Abstract The multispecies coalescent (MSC) model provides a natural framework for species tree estimation accounting for gene-tree conflicts. While a number of species tree methods under the MSC have been suggested and evaluated using simulation, their statistical properties remain poorly understood. Here we use mathematical analysis aided by computer simulation to examine the identifiability, consistency, and efficiency of different species tree methods in the case of three species and three sequences under the molecular clock. We consider four major species-tree methods including concatenation, two-step, independent-sites maximum likelihood (ISML) and maximum likelihood (ML). We develop approximations that predict that the probit transform of the species tree estimation error decreases linearly with the square root of the number of loci. Even in this simplest case major differences exist among the methods. Fulllikelihood methods are considerably more efficient than summary methods such as concatenation and two-step. They also provide estimates of important parameters such as species divergence times and ancestral population sizes while these parameters are not identifiable by summary methods. Our results highlight the need to improve the statistical efficiency of summary methods and the computational efficiency of full likelihood methods of species tree estimation.

scholarly journals Probabilities of Unranked and Ranked Anomaly Zones under Birth–Death Models

2019 ◽  
Vol 37 (5) ◽  
pp. 1480-1494 ◽  
Author(s):  
Anastasiia Kim ◽  
Noah A Rosenberg ◽  
James H Degnan
Keyword(s):  
Gene Tree ◽  
Species Tree ◽  
Tree Topology ◽  
Death Process ◽  
Gene Trees ◽  
Species Trees ◽  
Multispecies Coalescent ◽  
Speciation Rates ◽  
Tree Topologies ◽  
Birth Death

Abstract A labeled gene tree topology that is more probable than the labeled gene tree topology matching a species tree is called “anomalous.” Species trees that can generate such anomalous gene trees are said to be in the “anomaly zone.” Here, probabilities of “unranked” and “ranked” gene tree topologies under the multispecies coalescent are considered. A ranked tree depicts not only the topological relationship among gene lineages, as an unranked tree does, but also the sequence in which the lineages coalesce. In this article, we study how the parameters of a species tree simulated under a constant-rate birth–death process can affect the probability that the species tree lies in the anomaly zone. We find that with more than five taxa, it is possible for species trees to have both anomalous unranked and ranked gene trees. The probability of being in either type of anomaly zone increases with more taxa. The probability of anomalous gene trees also increases with higher speciation rates. We observe that the probabilities of unranked anomaly zones are higher and grow much faster than those of ranked anomaly zones as the speciation rate increases. Our simulation shows that the most probable ranked gene tree is likely to have the same unranked topology as the species tree. We design the software PRANC, which computes probabilities of ranked gene tree topologies given a species tree under the coalescent model.

scholarly journals Consistency of SVDQuartets and Maximum Likelihood for Coalescent-Based Species Tree Estimation

2020 ◽  
Vol 70 (1) ◽  
pp. 33-48 ◽  
Author(s):  
Matthew Wascher ◽  
Laura Kubatko
Keyword(s):  
Maximum Likelihood ◽  
Gene Tree ◽  
Species Tree ◽  
Gene Trees ◽  
Data Types ◽  
Multispecies Coalescent ◽  
Consistency Results ◽  
True Tree ◽  
Tree Inference ◽  
Statistical Consistency

Abstract Numerous methods for inferring species-level phylogenies under the coalescent model have been proposed within the last 20 years, and debates continue about the relative strengths and weaknesses of these methods. One desirable property of a phylogenetic estimator is that of statistical consistency, which means intuitively that as more data are collected, the probability that the estimated tree has the same topology as the true tree goes to 1. To date, consistency results for species tree inference under the multispecies coalescent (MSC) have been derived only for summary statistics methods, such as ASTRAL and MP-EST. These methods have been found to be consistent given true gene trees but may be inconsistent when gene trees are estimated from data for loci of finite length. Here, we consider the question of statistical consistency for four taxa for SVDQuartets for general data types, as well as for the maximum likelihood (ML) method in the case in which the data are a collection of sites generated under the MSC model such that the sites are conditionally independent given the species tree (we call these data coalescent independent sites [CIS] data). We show that SVDQuartets is statistically consistent for all data types (i.e., for both CIS data and for multilocus data), and we derive its rate of convergence. We additionally show that ML is consistent for CIS data under the JC69 model and discuss why a proof for the more general multilocus case is difficult. Finally, we compare the performance of ML and SDVQuartets using simulation for both data types. [Consistency; gene tree; maximum likelihood; multilocus data; hylogenetic inference; species tree; SVDQuartets.]

scholarly journals Retroposon Insertions within a Multispecies Coalescent Framework Suggest that Ratite Phylogeny is not in the ‘Anomaly Zone’

2019 ◽  
Author(s):  
Mark S. Springer ◽  
John Gatesy
Keyword(s):  
Phylogenetic Analysis ◽  
Sequence Data ◽  
Gene Tree ◽  
Reconstruction Error ◽  
Species Tree ◽  
Neutral Evolution ◽  
Gene Trees ◽  
Tree Reconstruction ◽  
Multispecies Coalescent ◽  
Branch Lengths

ABSTRACTSummary coalescence methods were developed to address the negative impacts of incomplete lineage sorting on species tree estimation with concatenation. Coalescence methods are statistically consistent if certain requirements are met including no intralocus recombination, neutral evolution, and no gene tree reconstruction error. However, the assumption of no intralocus recombination may not hold for many DNA sequence data sets, and neutral evolution is not the rule for genetic markers that are commonly employed in phylogenomic coalescence analyses. Most importantly, the assumption of no gene tree reconstruction error is routinely violated, especially for rapid radiations that are deep in the Tree of Life. With the sequencing of complete genomes and novel pipelines, phylogenetic analysis of retroposon insertions has emerged as a valuable alternative to sequence-based phylogenetic analysis. Retroposon insertions avoid or reduce several problems that beset analysis of sequence data with summary coalescence methods: 1) intralocus recombination is avoided because retroposon insertions are singular evolutionary events, 2) neutral evolution is approximated in many cases, and 3) gene tree reconstruction errors are rare because retroposons have low rates of homoplasy. However, the analysis of retroposons within a multispecies coalescent framework has not been realized. Here, we propose a simple workaround in which a retroposon insertion matrix is first transformed into a series of incompletely resolved gene trees. Next, the program ASTRAL is used to estimate a species tree in the statistically consistent framework of the multispecies coalescent. The inferred species tree includes support scores at all nodes and internal branch lengths in coalescent units. As a test case, we analyzed a retroposon dataset for palaeognath birds (ratites and tinamous) with ASTRAL and compared the resulting species tree to an MP-EST species tree for the same clade derived from thousands of sequence-based gene trees. The MP-EST species tree suggests an empirical case of the ‘anomaly zone’ with three very short internal branches at the base of Palaeognathae, and as predicted for anomaly zone conditions, the MP-EST species tree differs from the most common gene tree. Although identical in topology to the MP-EST tree, the ASTRAL species tree based on retroposons shows branch lengths that are much longer and incompatible with anomaly zone conditions. Simulation of gene trees from the retroposon-based species tree reveals that the most common gene tree matches the species tree. We contend that the wide discrepancies in branch lengths between sequence-based and retroposon-based species trees are explained by the greater accuracy of retroposon gene trees (bipartitions) relative to sequence-based gene trees. Coalescence analysis of retroposon data provides a promising alternative to the status quo by reducing gene tree reconstruction error that can have large impacts on both branch length estimates and evolutionary interpretations.

scholarly journals Phylogenomic terraces: presence and implication in species tree estimation from gene trees

2020 ◽  
Author(s):  
Ishrat Tanzila Farah ◽  
Md Muktadirul Islam ◽  
Kazi Tasnim Zinat ◽  
Atif Hasan Rahman ◽  
Md Shamsuzzoha Bayzid
Keyword(s):  
Incomplete Lineage Sorting ◽  
Gene Tree ◽  
Species Tree ◽  
Gene Trees ◽  
Species Trees ◽  
Lineage Sorting ◽  
Deep Coalescence ◽  
Tree Estimation ◽  
Tree Space ◽  
Multiple Species

AbstractSpecies tree estimation from multi-locus dataset is extremely challenging, especially in the presence of gene tree heterogeneity across the genome due to incomplete lineage sorting (ILS). Summary methods have been developed which estimate gene trees and then combine the gene trees to estimate a species tree by optimizing various optimization scores. In this study, we have formalized the concept of “phylogenomic terraces” in the species tree space, where multiple species trees with distinct topologies may have exactly the same optimization score (quartet score, extra lineage score, etc.) with respect to a collection of gene trees. We investigated the presence and implication of terraces in species tree estimation from multi-locus data by taking ILS into account. We analyzed two of the most popular ILS-aware optimization criteria: maximize quartet consistency (MQC) and minimize deep coalescence (MDC). Methods based on MQC are provably statistically consistent, whereas MDC is not a consistent criterion for species tree estimation. Our experiments, on a collection of dataset simulated under ILS, indicate that MDC-based methods may achieve competitive or identical quartet consistency score as MQC but could be significantly worse than MQC in terms of tree accuracy – demonstrating the presence and affect of phylogenomic terraces. This is the first known study that formalizes the concept of phylogenomic terraces in the context of species tree estimation from multi-locus data, and reports the presence and implications of terraces in species tree estimation under ILS.

scholarly journals STEM: species tree estimation using maximum likelihood for gene trees under coalescence

2009 ◽  
Vol 25 (7) ◽  
pp. 971-973 ◽  
Author(s):  
Laura S. Kubatko ◽  
Bryan C. Carstens ◽  
L. Lacey Knowles
Keyword(s):  
Maximum Likelihood ◽  
Species Tree ◽  
Gene Trees ◽  
Tree Estimation ◽  
Stem Species

The Multispecies Coalescent in Space and Time

2020 ◽  
Author(s):  
Patrick F. McKenzie ◽  
Deren A. R. Eaton
Keyword(s):  
Gene Tree ◽  
Species Tree ◽  
Gene Trees ◽  
Effective Population ◽  
Sliding Windows ◽  
Multispecies Coalescent ◽  
A Genome ◽  
Tree Inference ◽  
Rate Of Decay ◽  
Coalescent Time

AbstractA key distinction between species tree inference under the multi-species coalescent model (MSC), and the inference of gene trees in sliding windows along a genome, is in the effect of genetic linkage. Whereas the MSC explicitly assumes genealogies to be unlinked, i.e., statistically independent, genealogies located close together on genomes are spatially auto-correlated. Here we use tree sequence simulations with recombination to explore the effects of species tree parameters on spatial patterns of linkage among genealogies. We decompose coalescent time units to demonstrate differential effects of generation time and effective population size on spatial coalescent patterns, and we define a new metric, “phylogenetic linkage,” for measuring the rate of decay of phylogenetic similarity by comparison to distances among unlinked genealogies. Finally, we provide a simple example where accounting for phylogenetic linkage in sliding window analyses improves local gene tree inference.

scholarly journals DISSECT: an assignment-free Bayesian discovery method for species delimitation under the multispecies coalescent

2014 ◽  
Author(s):  
Graham Jones ◽  
Bengt Oxelman
Keyword(s):  
Simulated Data ◽  
Degree Of Approximation ◽  
Species Tree ◽  
Supplementary Information ◽  
Data Set ◽  
Multispecies Coalescent ◽  
Formal Framework ◽  
Beast Analysis ◽  
Tree Topologies ◽  
Discovery Method

Motivation: The multispecies coalescent model provides a formal framework for the assignment of individual organisms to species, where the species are modeled as the branches of the species tree. None of the available approaches so far have simultaneously co-estimated all the relevant parameters in the model, without restricting the parameter space by requiring a guide tree and/or prior assignment of individuals to clusters or species. Results: We present DISSECT, which explores the full space of possible clusterings of individuals and species tree topologies in a Bayesian framework. It uses an approximation to avoid the need for reversible-jump MCMC, in the form of a prior that is a modification of the birth-death prior for the species tree. It incorporates a spike near zero in the density for node heights. The model has two extra parameters: one controls the degree of approximation, and the second controls the prior distribution on the numbers of species. It is implemented as part of BEAST and requires only a few changes from a standard *BEAST analysis. The method is evaluated on simulated data and demonstrated on an empirical data set. The method is shown to be insensitive to the degree of approximation, but quite sensitive to the second parameter, suggesting that large numbers of sequences are needed to draw firm conclusions. Availability:http://code.google.com/p/beast-mcmc/, http://www.indriid.com/dissectinbeast.html Contact:[email protected], www.indriid.com Supplementary information: Supplementary material is available.

scholarly journals Consistency of SVDQuartets and Maximum Likelihood for Coalescent-based Species Tree Estimation

2019 ◽  
Author(s):  
Matthew Wascher ◽  
Laura Kubatko
Keyword(s):  
Maximum Likelihood ◽  
Species Tree ◽  
Gene Trees ◽  
Data Types ◽  
Coalescent Model ◽  
Multispecies Coalescent ◽  
Consistency Results ◽  
True Tree ◽  
Tree Inference ◽  
Statistical Consistency

AbtractNumerous methods for inferring species-level phylogenies under the coalescent model have been proposed within the last 20 years, and debates continue about the relative strengths and weaknesses of these methods. One desirable property of a phylogenetic estimator is that of statistical consistency, which means intuitively that as more data are collected, the probability that the estimated tree has the same topology as the true tree goes to 1. To date, consistency results for species tree inference under the multispecies coalescent have been derived only for summary statistics methods, such as ASTRAL and MP-EST. These methods have been found to be consistent given true gene trees, but may be inconsistent when gene trees are estimated from data for loci of finite length (Roch et al., 2019). Here we consider the question of statistical consistency for four taxa for SVDQuartets for general data types, as well as for the maximum likelihood (ML) method in the case in which the data are a collection of sites generated under the multispecies coalescent model such that the sites are conditionally independent given the species tree (we call these data Coalescent Independent Sites (CIS) data). We show that SVDQuartets is statistically consistent for all data types (i.e., for both CIS data and for multilocus data), and we derive its rate of convergence. We additionally show that ML is consistent for CIS data under the JC69 model, and discuss why a proof for the more general multilocus case is difficult. Finally, we compare the performance of maximum likelihood and SDVQuartets using simulation for both data types.

Sign in / Sign up

Export Citation Format

Share Document