scholarly journals A Distance Method to Reconstruct Species Trees In the Presence of Gene Flow

2014 ◽  
Author(s):  
Lingfei Cui ◽  
Laura Kubatko
Keyword(s):  
Gene Flow ◽  
Evolutionary Biology ◽  
Sequence Data ◽  
Divergence Time ◽  
Species Tree ◽  
Estimation Accuracy ◽  
Species Trees ◽  
Distance Method ◽  
Data Set ◽  
Divergence Time Estimates

One of the central tasks in evolutionary biology is to reconstruct the evolutionary relationships among species from sequence data, particularly from multilocus data. In the last ten years, many methods have been proposed to use the variance in the gene histories to estimate species trees by explicitly modeling deep coalescence. However, gene flow, another process that may produce gene history variance, has been less studied. In this paper, we propose a simple yet innovative method for species trees estimation in the presence of gene flow. Our method, called STEST (Species Tree Estimation from Speciation Times), constructs species tree estimates from pairwise speciation time or species divergence time estimates. By using methods that estimate speciation times in the presence of gene flow, (for example, M1 (Yang 2010) or SIM3s (Zhu and Yang 2012)), STEST is able to estimate species trees from data subject to gene flow. We develop two methods, called STEST (M1) and STEST (SIM3s), for this purpose. Additionally, we consider the method STEST (M0), which instead uses the M0 method (Yang 2002), a coalescent-based method that does not assume gene flow, to estimate speciation times. It is therefore devised to estimate species trees in the absence of gene flow. Our simulation studies show that STEST (M0) outperforms STEST(M1), STEST (SIM3s) and STEM in terms of estimation accuracy and outperfroms *BEAST in terms of running time when the degree of gene flow is small. STEST (M1) outperforms STEST (M0), STEST (SIM3s), STEM and *BEAST in term of estimation accuracy when the degree of gene flow is large. An empirical data set analyzed by these methods gives species tree estimates that are consistent with the previous results.

scholarly journals Inference of species histories in the presence of gene flow

2018 ◽  
Author(s):  
Nicola F. Müller ◽  
Huw A. Ogilvie ◽  
Chi Zhang ◽  
Alexei Drummond ◽  
Tanja Stadler
Keyword(s):  
Gene Flow ◽  
Sequence Data ◽  
Real Data ◽  
Monte Carlo Sampling ◽  
Species Tree ◽  
Species Trees ◽  
Population Isolation ◽  
Data Set ◽  
Isolation With Migration ◽  
Over Time

AbstractWhen populations become isolated, members of these populations can diverge genetically over time. This leads to genetic differences between individuals of these populations that increase over time if the isolation persists. This process can be counteracted when genes are exchanged between populations. In order to study the speciation processes when gene flow is present, isolation-with-migration methods have been developed. These methods typically assume that the ranked topology of the species history is already known. However, this is often not the case and the species tree is therefore of interest itself. To infer it is currently only possible when assuming no gene flow. This assumption can lead to wrongly inferred speciation times and species tree topologies.Building on a recently introduced structured coalescent approach, we introduce a new method that allows inference of the species tree while explicitly modelling the flow of genes between coexisting species. By using Markov chain Monte Carlo sampling, we co-infer the species tree alongside evolutionary parameters of interest. By using simulations, we show that our newly introduced approach is able to reliably infer the species trees and parameters of the isolation-with-migration model from genetic sequence data. We then infer the species history of six great ape species including gene flow after population isolation. By using this dataset, we are able to show that our new methods is able to infer the correct species tree not only on simulated but also on a real data set where the species history has already been well studied. In line with previous results, we find some support for some gene flow between bonobos and common chimpanzees.

scholarly journals COAL_PHYRE: A Composite Likelihood Method for Estimating Species Tree Parameters from Genomic Data Using Coalescent Theory

2020 ◽  
Author(s):  
Geno Guerra ◽  
Rasmus Nielsen
Keyword(s):  
Population Size ◽  
Gene Tree ◽  
Divergence Time ◽  
Species Tree ◽  
Composite Likelihood ◽  
Likelihood Method ◽  
Computational Time ◽  
Species Trees ◽  
Divergence Time Estimates ◽  
Full Likelihood

2AbstractGenome-scale data are increasingly being used to infer phylogenetic trees. A major challenge in such inferences is that different regions of the genome may have local topologies that differ from the species tree due to incomplete lineage sorting (ILS). Another source of gene tree discrepancies is estimation errors arising from the randomness of the mutational process during sequence evolution. There are two major groups of methods for estimating species tree from whole-genome data: a set of full likelihood methods, which model both sources of variance, but do not scale to large numbers of independent loci, and a class of faster approximation methods which do not model the mutational variance.To bridge the gap between these two classes of methods, we present COAL_PHYRE (COmposite Approximate Likelihood for PHYlogenetic REconstruction), a composite likelihood based method for inferring population size and divergence time estimates of rooted species trees from aligned gene sequences. COAL_PHYRE jointly models coalescent variation across loci using the MSC and variation in local gene tree reconstruction using a normal approximation. To evaluate the accuracy and speed of the method, we compare against BPP, a powerful MCMC full-likelihood method, as well as ASTRAL-III, a fast approximate method. We show that COAL_PHYRE’s divergence time and population size estimates are more accurate than ASTRAL, and comparable to those obtained using BPP, with an order of magnitude decrease in computational time. We also present results on previously published data from a set of Gibbon species to evaluate the accuracy in topology and parameter inference on real data, and to illustrate the method’s ability to analyze data sets which are prohibitively large for MCMC methods.

scholarly journals Anchored Hybrid Enrichment-Based Phylogenomics of Leafhoppers and Treehoppers (Hemiptera: Cicadomorpha: Membracoidea)

2017 ◽  
Vol 1 (1) ◽  
pp. 57-72 ◽  
Author(s):  
Christopher H Dietrich ◽  
Julie M Allen ◽  
Alan R Lemmon ◽  
Emily Moriarty Lemmon ◽  
Daniela M Takiya ◽  
...  
Keyword(s):  
Dna Sequences ◽  
Sequence Data ◽  
Single Gene ◽  
Gene Tree ◽  
Divergence Time ◽  
Data Set ◽  
Nucleotide Sequence Data ◽  
Divergence Time Estimates ◽  
Morphological Criteria ◽  
Anchored Hybrid Enrichment

Abstract A data set comprising DNA sequences from 388 loci and >99,000 aligned nucleotide positions, generated using anchored hybrid enrichment, was used to estimate relationships among 138 leafhoppers and treehoppers representative of all major lineages of Membracoidea, the most diverse superfamily of hemipteran insects. Phylogenetic analysis of the concatenated nucleotide sequence data set using maximum likelihood produced a tree with most branches receiving high support. A separate coalescent gene tree analysis of the same data generally recovered the same strongly supported clades but was less well resolved overall. Several nodes pertaining to relationships among leafhopper subfamilies currently recognized based on morphological criteria were separated by short internodes and received low support. Although various higher taxa were corroborated with improved branch support, relationships among some major lineages of Membracoidea are only somewhat more resolved than previously published phylogenies based on single gene regions or morphology. In agreement with previous studies, the present results indicate that leafhoppers (Cicadellidae) are paraphyletic with respect to the three recognized families of treehoppers (Aetalionidae, Melizoderidae, and Membracidae). Divergence time estimates indicate that most of the poorly resolved divergence events among major leafhopper lineages occurred during the lower to middle Cretaceous and that most modern leafhopper subfamilies, as well as the lineage comprising the three recognized families of treehoppers, also arose during the Cretaceous.

scholarly journals Strategies for Partitioning Clock Models in Phylogenomic Dating: Application to the Angiosperm Evolutionary Timescale

2017 ◽  
Author(s):  
Charles S. P. Foster ◽  
Simon Y. W. Ho
Keyword(s):  
Sequence Data ◽  
Divergence Time ◽  
Molecular Dating ◽  
Sequence Length ◽  
Rate Heterogeneity ◽  
Data Set ◽  
Molecular Sequence Data ◽  
Molecular Sequence ◽  
Divergence Time Estimates ◽  
Time Estimates

AbstractEvolutionary timescales can be inferred from molecular sequence data using a Bayesian phylogenetic approach. In these methods, the molecular clock is often calibrated using fossil data. The uncertainty in these fossil calibrations is important because it determines the limiting posterior distribution for divergence-time estimates as the sequence length tends to infinity. Here we investigate how the accuracy and precision of Bayesian divergence-time estimates improve with the increased clock-partitioning of genome-scale data into clock-subsets. We focus on a data set comprising plastome-scale sequences of 52 angiosperm taxa. There was little difference among the Bayesian date estimates whether we chose clock-subsets based on patterns of among-lineage rate heterogeneity or relative rates across genes, or by random assignment. Increasing the degree of clock-partitioning usually led to an improvement in the precision of divergence-time estimates, but this increase was asymptotic to a limit presumably imposed by fossil calibrations. Our clock-partitioning approaches yielded highly precise age estimates for several key nodes in the angiosperm phylogeny. For example, when partitioning the data into 20 clock-subsets based on patterns of among-lineage rate heterogeneity, we inferred crown angiosperms to have arisen 198–178 Ma. This demonstrates that judicious clock-partitioning can improve the precision of molecular dating based on phylogenomic data, but the meaning of this increased precision should be considered critically.

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 Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l

2020 ◽  
Author(s):  
Diego F Morales-Briones ◽  
Gudrun Kadereit ◽  
Delphine T Tefarikis ◽  
Michael J Moore ◽  
Stephen A Smith ◽  
...  
Keyword(s):  
Network Inference ◽  
Incomplete Lineage Sorting ◽  
Gene Tree ◽  
Phylogenetic Network ◽  
Species Tree ◽  
Data Sets ◽  
Species Trees ◽  
Lineage Sorting ◽  
Data Set ◽  
Gene Tree Discordance

Abstract Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations. [Amaranthaceae; gene tree discordance; hybridization; incomplete lineage sorting; phylogenomics; species network; species tree; transcriptomics.]

scholarly journals Uprooting phylogenetic uncertainty in coalescent species delimitation: A meta-analysis of empirical studies

2015 ◽  
Vol 61 (5) ◽  
pp. 866-873 ◽  
Author(s):  
Itzue W. Caviedes-Solis ◽  
Nassima M. Bouzid ◽  
Barbara L. Banbury ◽  
Adam D. Leaché
Keyword(s):  
Species Delimitation ◽  
Sequence Data ◽  
Meta Analysis ◽  
Empirical Studies ◽  
Species Tree ◽  
Species Trees ◽  
Posterior Probabilities ◽  
Guide Tree ◽  
Tree Topologies ◽  
Accurate Quantification

Abstract Phylogenetic and phylogeographic studies rely on the accurate quantification of biodiversity. In recent studies of taxonomically ambiguous groups, species boundaries are often determined based on multi-locus sequence data. Bayesian Phylogenetics and Phylogeography (BPP) is a coalescent-based method frequently used to delimit species; however, empirical studies suggest that the requirement of a user-specified guide tree biases the range of possible outcomes. We evaluate fifteen multi-locus datasets using the most recent iteration of BPP, which eliminates the need for a user-specified guide tree and reconstructs the species tree in synchrony with species delimitation (= unguided species delimitation). We found that the number of species recovered with guided versus unguided species delimitation was the same except for two cases, and that posterior probabilities were generally lower for the unguided analyses as a result of searching across species trees in addition to species delimitation models. The guide trees used in previous studies were often discordant with the species tree topologies estimated by BPP. We also compared species trees estimated using BPP and *BEAST and found that when the topologies are the same, BPP tends to give higher posterior probabilities.

scholarly journals Terraces in Gene Tree Reconciliation-Based Species Tree Inference

2020 ◽  
Author(s):  
Michael J. Sanderson ◽  
Michelle M. McMahon ◽  
Mike Steel
Keyword(s):  
Incomplete Lineage Sorting ◽  
Gene Tree ◽  
Solution Space ◽  
Species Tree ◽  
Gene Trees ◽  
Species Trees ◽  
Lineage Sorting ◽  
Data Set ◽  
Tree Reconciliation ◽  
The Impact

AbstractTerraces in phylogenetic tree space are sets of trees with identical optimality scores for a given data set, arising from missing data. These were first described for multilocus phylogenetic data sets in the context of maximum parsimony inference and maximum likelihood inference under certain model assumptions. Here we show how the mathematical properties that lead to terraces extend to gene tree - species tree problems in which the gene trees are incomplete. Inference of species trees from either sets of gene family trees subject to duplication and loss, or allele trees subject to incomplete lineage sorting, can exhibit terraces in their solution space. First, we show conditions that lead to a new kind of terrace, which stems from subtree operations that appear in reconciliation problems for incomplete trees. Then we characterize when terraces of both types can occur when the optimality criterion for tree search is based on duplication, loss or deep coalescence scores. Finally, we examine the impact of assumptions about the causes of losses: whether they are due to imperfect sampling or true evolutionary deletion.

scholarly journals Using Phylogenomic Data to Explore the Effects of Relaxed Clocks and Calibration Strategies on Divergence Time Estimation: Primates as a Test Case

2017 ◽  
Author(s):  
Mario dos Reis ◽  
Gregg F. Gunnell ◽  
José Barba-Montoya ◽  
Alex Wilkins ◽  
Ziheng Yang ◽  
...  
Keyword(s):  
Sequence Data ◽  
Divergence Time ◽  
Time Estimation ◽  
Test Case ◽  
Base Pairs ◽  
Molecular Sequence Data ◽  
Divergence Time Estimation ◽  
Divergence Time Estimates ◽  
Time Estimates ◽  
Clock Models

AbstractPrimates have long been a test case for the development of phylogenetic methods for divergence time estimation. Despite a large number of studies, however, the timing of origination of crown Primates relative to the K-Pg boundary and the timing of diversification of the main crown groups remain controversial. Here we analysed a dataset of 372 taxa (367 Primates and 5 outgroups, 61 thousand base pairs) that includes nine complete primate genomes (3.4 million base pairs). We systematically explore the effect of different interpretations of fossil calibrations and molecular clock models on primate divergence time estimates. We find that even small differences in the construction of fossil calibrations can have a noticeable impact on estimated divergence times, especially for the oldest nodes in the tree. Notably, choice of molecular rate model (auto-correlated or independently distributed rates) has an especially strong effect on estimated times, with the independent rates model producing considerably more ancient estimates for the deeper nodes in the phylogeny. We implement thermodynamic integration, combined with Gaussian quadrature, in the program MCMCTree, and use it to calculate Bayes factors for clock models. Bayesian model selection indicates that the auto-correlated rates model fits the primate data substantially better, and we conclude that time estimates under this model should be preferred. We show that for eight core nodes in the phylogeny, uncertainty in time estimates is close to the theoretical limit imposed by fossil uncertainties. Thus, these estimates are unlikely to be improved by collecting additional molecular sequence data. All analyses place the origin of Primates close to the K-Pg boundary, either in the Cretaceous or straddling the boundary into the Palaeogene.

Sign in / Sign up

Export Citation Format

Share Document