scholarly journals Computing the probability of gene trees concordant with the species tree in the multispecies coalescent

2020 ◽  
Author(s):  
Jakub Truszkowski ◽  
Celine Scornavacca ◽  
Fabio Pardi
Keyword(s):  
Species Tree ◽  
Gene Trees ◽  
Multispecies Coalescent

scholarly journals Impact of Ghost Introgression on Coalescent-based Species Tree Inference and Estimation of Divergence Time

2022 ◽  
Author(s):  
XiaoXu Pang ◽  
Da-Yong Zhang
Keyword(s):  
Incomplete Lineage Sorting ◽  
Divergence Time ◽  
Species Tree ◽  
Gene Trees ◽  
Species Trees ◽  
Lineage Sorting ◽  
Multispecies Coalescent ◽  
Tree Inference ◽  
Tree Methods ◽  
The Impact

The species studied in any evolutionary investigation generally constitute a very small proportion of all the species currently existing or that have gone extinct. It is therefore likely that introgression, which is widespread across the tree of life, involves "ghosts," i.e., unsampled, unknown, or extinct lineages. However, the impact of ghost introgression on estimations of species trees has been rarely studied and is thus poorly understood. In this study, we use mathematical analysis and simulations to examine the robustness of species tree methods based on a multispecies coalescent model under gene flow sourcing from an extant or ghost lineage. We found that very low levels of extant or ghost introgression can result in anomalous gene trees (AGTs) on three-taxon rooted trees if accompanied by strong incomplete lineage sorting (ILS). In contrast, even massive introgression, with more than half of the recipient genome descending from the donor lineage, may not necessarily lead to AGTs. In cases involving an ingroup lineage (defined as one that diverged no earlier than the most basal species under investigation) acting as the donor of introgression, the time of root divergence among the investigated species was either underestimated or remained unaffected, but for the cases of outgroup ghost lineages acting as donors, the divergence time was generally overestimated. Under many conditions of ingroup introgression, the stronger the ILS was, the higher was the accuracy of estimating the time of root divergence, although the topology of the species tree is more prone to be biased by the effect of introgression.

scholarly journals Practical Speedup of Bayesian Inference of Species Phylogenies by Restricting the Space of Gene Trees

2019 ◽  
Author(s):  
Yaxuan Wang ◽  
Huw A. Ogilvie ◽  
Luay Nakhleh
Keyword(s):  
Bayesian Inference ◽  
Species Tree ◽  
Biological Data ◽  
Small Data ◽  
Data Sets ◽  
Gene Trees ◽  
Multispecies Coalescent ◽  
Tree Inference ◽  
The Individual ◽  
Species Tree Inference

AbstractSpecies tree inference from multi-locus data has emerged as a powerful paradigm in the post-genomic era, both in terms of the accuracy of the species tree it produces as well as in terms of elucidating the processes that shaped the evolutionary history. Bayesian methods for species tree inference are desirable in this area as they have been shown to yield accurate estimates, but also to naturally provide measures of confidence in those estimates. However, the heavy computational requirements of Bayesian inference have limited the applicability of such methods to very small data sets.In this paper, we show that the computational efficiency of Bayesian inference under the multispecies coalescent can be improved in practice by restricting the space of the gene trees explored during the random walk, without sacrificing accuracy as measured by various metrics. The idea is to first infer constraints on the trees of the individual loci in the form of unresolved gene trees, and then to restrict the sampler to consider only resolutions of the constrained trees. We demonstrate the improvements gained by such an approach on both simulated and biological data.

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 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 Practical Speedup of Bayesian Inference of Species Phylogenies by Restricting the Space of Gene Trees

2020 ◽  
Vol 37 (6) ◽  
pp. 1809-1818
Author(s):  
Yaxuan Wang ◽  
Huw A Ogilvie ◽  
Luay Nakhleh
Keyword(s):  
Bayesian Inference ◽  
Species Tree ◽  
Biological Data ◽  
Small Data ◽  
Data Sets ◽  
Gene Trees ◽  
Multispecies Coalescent ◽  
Tree Inference ◽  
The Individual ◽  
Species Tree Inference

Abstract Species tree inference from multilocus data has emerged as a powerful paradigm in the postgenomic era, both in terms of the accuracy of the species tree it produces as well as in terms of elucidating the processes that shaped the evolutionary history. Bayesian methods for species tree inference are desirable in this area as they have been shown not only to yield accurate estimates, but also to naturally provide measures of confidence in those estimates. However, the heavy computational requirements of Bayesian inference have limited the applicability of such methods to very small data sets. In this article, we show that the computational efficiency of Bayesian inference under the multispecies coalescent can be improved in practice by restricting the space of the gene trees explored during the random walk, without sacrificing accuracy as measured by various metrics. The idea is to first infer constraints on the trees of the individual loci in the form of unresolved gene trees, and then to restrict the sampler to consider only resolutions of the constrained trees. We demonstrate the improvements gained by such an approach on both simulated and biological data.

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

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.

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 Gene tree discord, simplex plots, and statistical tests under the coalescent

2020 ◽  
Author(s):  
Elizabeth S. Allman ◽  
Jonathan D. Mitchell ◽  
John A. Rhodes
Keyword(s):  
Incomplete Lineage Sorting ◽  
Statistical Tests ◽  
Gene Tree ◽  
Complex Model ◽  
Species Tree ◽  
Data Sets ◽  
Gene Trees ◽  
Lineage Sorting ◽  
Multispecies Coalescent ◽  
Tree Inference

AbstractA simple graphical device, the simplex plot of quartet concordance factors, is introduced to aid in the exploration of a collection of gene trees on a common set of taxa. A single plot summarizes all gene tree discord, and allows for visual comparison to the expected discord from the multispecies coalescent model (MSC) of incomplete lineage sorting on a species tree. A formal statistical procedure is described that can quantify the deviation from expectation for each subset of four taxa, suggesting when the data is not in accord with the MSC, and thus either gene tree inference error is substantial or a more complex model such as that on a network may be required. If the collection of gene trees appears to be in accord with the MSC, the plots may reveal when substantial incomplete lineage sorting is present and coalescent based species tree inference is preferred over concatenation approaches. Applications to both simulated and empirical multilocus data sets illustrate the insights provided.

Sign in / Sign up

Export Citation Format

Share Document