scholarly journals A Bayesian Implementation of the Multispecies Coalescent Model with Introgression for Phylogenomic Analysis

2019 ◽  
Vol 37 (4) ◽  
pp. 1211-1223 ◽  
Author(s):  
Tomáš Flouri ◽  
Xiyun Jiao ◽  
Bruce Rannala ◽  
Ziheng Yang
Keyword(s):  
Gene Flow ◽  
Genomic Sequence ◽  
Sequence Data ◽  
Incomplete Lineage Sorting ◽  
Mosquito Species ◽  
Phylogenomic Analysis ◽  
Data Sets ◽  
Lineage Sorting ◽  
Coalescent Model ◽  
Multispecies Coalescent

Abstract Recent analyses suggest that cross-species gene flow or introgression is common in nature, especially during species divergences. Genomic sequence data can be used to infer introgression events and to estimate the timing and intensity of introgression, providing an important means to advance our understanding of the role of gene flow in speciation. Here, we implement the multispecies-coalescent-with-introgression model, an extension of the multispecies-coalescent model to incorporate introgression, in our Bayesian Markov chain Monte Carlo program Bpp. The multispecies-coalescent-with-introgression model accommodates deep coalescence (or incomplete lineage sorting) and introgression and provides a natural framework for inference using genomic sequence data. Computer simulation confirms the good statistical properties of the method, although hundreds or thousands of loci are typically needed to estimate introgression probabilities reliably. Reanalysis of data sets from the purple cone spruce confirms the hypothesis of homoploid hybrid speciation. We estimated the introgression probability using the genomic sequence data from six mosquito species in the Anopheles gambiae species complex, which varies considerably across the genome, likely driven by differential selection against introgressed alleles.

scholarly journals A Bayesian implementation of the multispecies coalescent model with introgression for comparative genomic analysis

2019 ◽  
Author(s):  
Thomas Flouris ◽  
Xiyun Jiao ◽  
Bruce Rannala ◽  
Ziheng Yang
Keyword(s):  
Gene Flow ◽  
Genomic Sequence ◽  
Sequence Data ◽  
Incomplete Lineage Sorting ◽  
Mosquito Species ◽  
Genomic Analysis ◽  
Comparative Genomic ◽  
Lineage Sorting ◽  
Multispecies Coalescent ◽  
Important Means

AbstractRecent analyses suggest that cross-species gene flow or introgression is common in nature, especially during species divergences. Genomic sequence data can be used to infer introgression events and to estimate the timing and intensity of introgression, providing an important means to advance our understanding of the role of gene flow in speciation. Here we implement the multispecies-coalescent-with-introgression (MSci) model, an extension of the multispecies-coalescent (MSC) model to incorporate introgression, in our Bayesian Markov chain Monte Carlo (MCMC) program BPP. The MSci model accommodates deep coalescence (or incomplete lineage sorting) and introgression and provides a natural framework for inference using genomic sequence data. Computer simulation confirms the good statistical properties of the method, although hundreds or thousands of loci are typically needed to estimate introgression probabilities reliably. Re-analysis of datasets from the purple cone spruce confirms the hypothesis of homoploid hybrid speciation. We estimated the introgression probability using the genomic sequence data from six mosquito species in the Anopheles gambiae species complex, which varies considerably across the genome, likely driven by differential selection against introgressed alleles.

scholarly journals MSCquartets 1.0: Quartet methods for species trees and networks under the multispecies coalescent model in R

2020 ◽  
Author(s):  
John A Rhodes ◽  
Hector Baños ◽  
Jonathan D Mitchell ◽  
Elizabeth S Allman
Keyword(s):  
Network Inference ◽  
Incomplete Lineage Sorting ◽  
R Package ◽  
Species Tree ◽  
Supplementary Information ◽  
Species Trees ◽  
Lineage Sorting ◽  
Coalescent Model ◽  
Multispecies Coalescent ◽  
Tree Inference

Abstract Summary MSCquartets is an R package for species tree hypothesis testing, inference of species trees, and inference of species networks under the Multispecies Coalescent model of incomplete lineage sorting and its network analog. Input for these analyses are collections of metric or topological locus trees which are then summarized by the quartets displayed on them. Results of hypothesis tests at user-supplied levels are displayed in a simplex plot by color-coded points. The package implements the QDC and WQDC algorithms for topological and metric species tree inference, and the NANUQ algorithm for level-1 topological species network inference, all of which give statistically consistent estimators under the model. Availability MSCquartets is available through the Comprehensive R Archive Network: https://CRAN.R-project.org/package=MSCquartets. Supplementary information Supplementary materials, including example data and analyses, are incorporated into the package.

scholarly journals Detecting destabilizing species in the phylogenetic backbone of Potentilla (Rosaceae) using low-copy nuclear markers

AoB Plants ◽  
2020 ◽  
Vol 12 (3) ◽  
Author(s):  
Nannie L Persson ◽  
Ingrid Toresen ◽  
Heidi Lie Andersen ◽  
Jenny E E Smedmark ◽  
Torsten Eriksson
Keyword(s):  
Incomplete Lineage Sorting ◽  
Gene Tree ◽  
Data Sets ◽  
Gene Trees ◽  
Lineage Sorting ◽  
Nuclear Markers ◽  
Data Set ◽  
Multispecies Coalescent ◽  
Single Marker ◽  
Named Groups

Abstract The genus Potentilla (Rosaceae) has been subjected to several phylogenetic studies, but resolving its evolutionary history has proven challenging. Previous analyses recovered six, informally named, groups: the Argentea, Ivesioid, Fragarioides, Reptans, Alba and Anserina clades, but the relationships among some of these clades differ between data sets. The Reptans clade, which includes the type species of Potentilla, has been noticed to shift position between plastid and nuclear ribosomal data sets. We studied this incongruence by analysing four low-copy nuclear markers, in addition to chloroplast and nuclear ribosomal data, with a set of Bayesian phylogenetic and Multispecies Coalescent (MSC) analyses. A selective taxon removal strategy demonstrated that the included representatives from the Fragarioides clade, P. dickinsii and P. fragarioides, were the main sources of the instability seen in the trees. The Fragarioides species showed different relationships in each gene tree, and were only supported as a monophyletic group in a single marker when the Reptans clade was excluded from the analysis. The incongruences could not be explained by allopolyploidy, but rather by homoploid hybridization, incomplete lineage sorting or taxon sampling effects. When P. dickinsii and P. fragarioides were removed from the data set, a fully resolved, supported backbone phylogeny of Potentilla was obtained in the MSC analysis. Additionally, indications of autopolyploid origins of the Reptans and Ivesioid clades were discovered in the low-copy gene trees.

scholarly journals The BPP program for species tree estimation and species delimitation

2015 ◽  
Vol 61 (5) ◽  
pp. 854-865 ◽  
Author(s):  
Ziheng Yang
Keyword(s):  
Species Delimitation ◽  
Genomic Sequence ◽  
Sequence Data ◽  
Species Tree ◽  
Joint Analysis ◽  
Coalescent Model ◽  
Species Phylogeny ◽  
Multispecies Coalescent ◽  
Tree Estimation ◽  
Nuclear Loci

Abstract This paper provides an overview and a tutorial of the BPP program, which is a Bayesian MCMC program for analyzing multi-locus genomic sequence data under the multispecies coalescent model. An example dataset of five nuclear loci from the East Asian brown frogs is used to illustrate four different analyses, including estimation of species divergence times and population size parameters under the multispecies coalescent model on a fixed species phylogeny (A00), species tree estimation when the assignment and species delimitation are fixed (A01), species delimitation using a fixed guide tree (A10), and joint species delimitation and species-tree estimation or unguided species delimitation (A11). For the joint analysis (A11), two new priors are introduced, which assign uniform probabilities for the different numbers of delimited species, which may be useful when assignment, species delimitation, and species phylogeny are all inferred in one joint analysis. The paper ends with a discussion of the assumptions, the strengths and weaknesses of the BPP analysis.

scholarly journals Identifying model violations under the multispecies coalescent model using P2C2M.SNAPP

PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e8271
Author(s):  
Drew J. Duckett ◽  
Tara A. Pelletier ◽  
Bryan C. Carstens
Keyword(s):  
Posterior Distribution ◽  
Incomplete Lineage Sorting ◽  
R Package ◽  
Predictive Distribution ◽  
Species Trees ◽  
Lineage Sorting ◽  
Posterior Predictive Distribution ◽  
Coalescent Model ◽  
Multispecies Coalescent ◽  
Phylogenetic Estimation

Phylogenetic estimation under the multispecies coalescent model (MSCM) assumes all incongruence among loci is caused by incomplete lineage sorting. Therefore, applying the MSCM to datasets that contain incongruence that is caused by other processes, such as gene flow, can lead to biased phylogeny estimates. To identify possible bias when using the MSCM, we present P2C2M.SNAPP. P2C2M.SNAPP is an R package that identifies model violations using posterior predictive simulation. P2C2M.SNAPP uses the posterior distribution of species trees output by the software package SNAPP to simulate posterior predictive datasets under the MSCM, and then uses summary statistics to compare either the empirical data or the posterior distribution to the posterior predictive distribution to identify model violations. In simulation testing, P2C2M.SNAPP correctly classified up to 83% of datasets (depending on the summary statistic used) as to whether or not they violated the MSCM model. P2C2M.SNAPP represents a user-friendly way for researchers to perform posterior predictive model checks when using the popular SNAPP phylogenetic estimation program. It is freely available as an R package, along with additional program details and tutorials.

Phylogenomics Uncovers Confidence and Conflict in the Rapid Radiation of Australo-Papuan Rodents

2019 ◽  
Vol 69 (3) ◽  
pp. 431-444 ◽  
Author(s):  
Emily J Roycroft ◽  
Adnan Moussalli ◽  
Kevin C Rowe
Keyword(s):  
Sequence Data ◽  
Incomplete Lineage Sorting ◽  
Phylogenetic Signal ◽  
Strong Support ◽  
Bootstrap Support ◽  
Data Sets ◽  
Lineage Sorting ◽  
Nucleotide Sequence Data ◽  
Branch Support ◽  
Partitioning Strategy

Abstract The estimation of robust and accurate measures of branch support has proven challenging in the era of phylogenomics. In data sets of potentially millions of sites, bootstrap support for bifurcating relationships around very short internal branches can be inappropriately inflated. Such overestimation of branch support may be particularly problematic in rapid radiations, where phylogenetic signal is low and incomplete lineage sorting severe. Here, we explore this issue by comparing various branch support estimates under both concatenated and coalescent frameworks, in the recent radiation Australo-Papuan murine rodents (Muridae: Hydromyini). Using nucleotide sequence data from 1245 independent loci and several phylogenomic inference methods, we unequivocally resolve the majority of genus-level relationships within Hydromyini. However, at four nodes we recover inconsistency in branch support estimates both within and among concatenated and coalescent approaches. In most cases, concatenated likelihood approaches using standard fast bootstrap algorithms did not detect any uncertainty at these four nodes, regardless of partitioning strategy. However, we found this could be overcome with two-stage resampling, that is, across genes and sites within genes (using -bsam GENESITE in IQ-TREE). In addition, low confidence at recalcitrant nodes was recovered using UFBoot2, a recent revision to the bootstrap protocol in IQ-TREE, but this depended on partitioning strategy. Summary coalescent approaches also failed to detect uncertainty under some circumstances. For each of four recalcitrant nodes, an equivalent (or close to equivalent) number of genes were in strong support ($>$ 75% bootstrap) of both the primary and at least one alternative topological hypothesis, suggesting notable phylogenetic conflict among loci not detected using some standard branch support metrics. Recent debate has focused on the appropriateness of concatenated versus multigenealogical approaches to resolving species relationships, but less so on accurately estimating uncertainty in large data sets. Our results demonstrate the importance of employing multiple approaches when assessing confidence and highlight the need for greater attention to the development of robust measures of uncertainty in the era of phylogenomics.

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.

scholarly journals Reticulate Evolutionary History in a Recent Radiation of Montane Grasshoppers Revealed by Genomic Data

2021 ◽  
Author(s):  
Vanina Tonzo ◽  
AdriÀ Bellvert ◽  
Joaquín Ortego
Keyword(s):  
Gene Flow ◽  
Biological Diversity ◽  
Sequence Data ◽  
Incomplete Lineage Sorting ◽  
Phylogenetic Analyses ◽  
Secondary Contact ◽  
Lineage Sorting ◽  
Glacial Cycles ◽  
Reciprocal Monophyly ◽  
Geometric Morphometric

AbstractInferring the ecological and evolutionary processes underlying lineage and phenotypic diversification is of paramount importance to shed light on the origin of contemporary patterns of biological diversity. However, reconstructing phylogenetic relationships in recent evolutionary radiations represents a major challenge due to the frequent co-occurrence of incomplete lineage sorting and introgression. In this study, we combined high throughput sequence data (ddRADseq), geometric morphometric information, and novel phylogenetic inference methods that explicitly account for gene flow to infer the evolutionary relationships and the timing and mode of diversification in a complex of Ibero-Maghrebian montane grasshoppers of the subgenus Dreuxius (genus Omocestus). Our analyses supported the phenotypic distinctiveness of most sister taxa, two events of historical introgression involving lineages at different stages of the diversification continuum, and the recent Pleistocene origin (< 1 Ma) of the complex. Phylogenetic analyses did not recover the reciprocal monophyly of taxa from Iberia and northwestern Africa, supporting overseas migration between the two continents during the Pleistocene. Collectively, these results indicate that periods of isolation and secondary contact linked to Pleistocene glacial cycles likely contributed to both allopatric speciation and post divergence gene flow in the complex. This study exemplifies how the integration of multiple lines of evidence can help to reconstruct complex histories of reticulated evolution and highlights the important role of Quaternary climatic oscillations as a diversification engine in the Ibero-Maghrebian biodiversity hotspot.

scholarly journals The genomic view of diversification

2018 ◽  
Author(s):  
Julie Marin ◽  
Guillaume Achaz ◽  
Anton Crombach ◽  
Amaury Lambert
Keyword(s):  
Gene Flow ◽  
Incomplete Lineage Sorting ◽  
Molecular Data ◽  
Species Tree ◽  
Data Sets ◽  
Gene Trees ◽  
Lineage Sorting ◽  
Parameters Tuning ◽  
History Of ◽  
New Framework

AbstractEvolutionary relationships between species are traditionally represented in the form of a tree, the species tree. Its reconstruction from molecular data is hindered by frequent conflicts between gene genealogies. Usually, these disagreements are explained by incomplete lineage sorting (ILS) due to random coalescences of gene lineages inside the edges of the species tree. This paradigm, the multi-species coalescent (MSC), is constantly violated by the ubiquitous presence of gene flow, leading to incongruences between gene trees that cannot be explained by ILS alone. Here we argue instead in favor of a vision acknowledging the importance of gene flow and where gene histories shape the species tree rather than the opposite. We propose a new framework for modeling the joint evolution of gene and species lineages relaxing the hierarchy between the species tree and gene trees. We implement this framework in two mathematical models called the gene-based diversification models (GBD): 1) GBD-forward following all evolving genomes and 2) GBD-backward based on coalescent theory. They feature four parameters tuning colonization, gene flow, genetic drift and genetic differentiation. We propose a quick inference method based on differences between gene trees. Applied to two empirical data-sets prone to gene flow, we find a better support for the GBD model than for the MSC model. Along with the increasing awareness of the extent of gene flow, this work shows the importance of considering the richer signal contained in genomic histories, rather than in the mere species tree, to better apprehend the complex evolutionary history of species.

scholarly journals Gene Tree Discord, Simplex Plots, and Statistical Tests under the Coalescent

2021 ◽  
Author(s):  
Elizabeth S Allman ◽  
Jonathan D Mitchell ◽  
John A Rhodes
Keyword(s):  
Incomplete Lineage Sorting ◽  
Hypothesis Test ◽  
Statistical Tests ◽  
Gene Tree ◽  
Complex Model ◽  
Species Tree ◽  
Gene Trees ◽  
Lineage Sorting ◽  
Coalescent Model ◽  
Multispecies Coalescent

Abstract A 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 are not in accord with the MSC, and thus that 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 is in accord with the MSC, the plots reveal when substantial incomplete lineage sorting is present. Applications to both simulated and empirical multilocus data sets illustrate the insights provided. [Gene tree discordance; hypothesis test; multispecies coalescent model; quartet concordance factor; simplex plot; species tree].

Sign in / Sign up

Export Citation Format

Share Document