Abstract:Abstract Phylogenies inferred using both concatenation- and coalescent-based analyses typically render highly congruent trees. However, when they disagree, they often differ with respect to historically contentious and evolutionarily important relationships. These relationships may also involve etiolated lineages where increased sampling is not possible. Recently, methods aimed at interrogating single relationships or trees have emerged as promising investigative tools to examine these cases. Although recent methods such as “Edge-based Phylogenomic Support analYsis” (EPSY) led to insights into both systematic error and real biological signal, whether they are consistent with the coalescent in cases with high Incomplete Lineage Sorting (ILS) has yet to be characterized. Here, we use simulations and an empirical dataset to test the performance of EPSY, concatenation, and coalescent-based summary analyses under high levels of ILS. We focused on high-ILS scenarios because these represent the typical difficult cases that researchers often face due to the prevalence of ILS in phylogenomic datasets. ILS is known to be a major cause of phylogenomic conflict, which confounds many biological conclusions that depend on a resolved phylogeny, such as inferring ancestral character states, biogeographic reconstructions, and domestication histories. Our study found that EPSY was consistent with the coalescent in a high-ILS empirical dataset. In high-ILS simulations EPSY infers the correct edge more than half the time, whereas coalescent based methods and concatenation methods inferred the actual tree 37.8% and 25% of the time, respectively. All methods have conditions under which they generate the most accurate inferences. Given the levels of ILS in simulations, 26.2% of the time no method recovered the true tree. This zone where no current method can infer the true topology is likely due to properties of the species tree, such as the length of internal edges adjacent to a conflict and/or the length of the shortest branch. Nevertheless, the EPSY approach proves to be a valuable complement to phylogenomic analyses for interrogating regions of the tree with conflicting hypotheses generated from past studies or alternative inference methods. Our analyses highlight that robust phylogenetic trees may not be possible under some scenarios regardless of method and data source.

Species Tree Inference Methods Intended to Deal with Incomplete Lineage Sorting Are Robust to the Presence of Paralogs

Species Tree Estimation and the Impact of Gene Loss Following Whole-Genome Duplication

Coalescent Methods Are Robust to the Simultaneous Effects of Long Branches and Incomplete Lineage Sorting

Using paralogs for phylogenetic inference mitigates the effects of long-branch attraction

When species trees disagree: an approach consistent with the coalescent that quantifies phylogenomic support for contentious relationships

Phylogenomics with Paralogs

A Distance-Based Method for Inferring Phylogenetic Networks in the Presence of Incomplete Lineage Sorting

ASTRAL-Pro: Quartet-Based Species-Tree Inference despite Paralogy

The Impact of Missing Data on Species Tree Estimation.

reconcILS: A gene tree-species tree reconciliation algorithm that allows for incomplete lineage sorting

The Multilocus Multispecies Coalescent: A Flexible New Model of Gene Family Evolution

Bayesian Inference of Species Trees from Multilocus Data

Genes with Minimal Phylogenetic Information Are Problematic for Coalescent Analyses when Gene Tree Estimation is Biased

Defining Loci in Restriction-Based Reduced Representation Genomic Data from Nonmodel Species: Sources of Bias and Diagnostics for Optimal Clustering

From gene trees to species trees II: Species tree inference in the deep coalescence model

DupLoss-2: Improved Phylogenomic Species Tree Inference under Gene Duplication and Loss

Large-scale Species Tree Estimation

The Inference of Gene Trees with Species Trees

Estimating phylogenetic trees from genome-scale data

Statistical Inference of Allopolyploid Species Networks in the Presence of Incomplete Lineage Sorting

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