Abstract:The genomic revolution offers renewed hope of resolving rapid radiations in the Tree of Life. The development of the multispecies coalescent model and improved gene tree estimation methods can better accommodate gene tree heterogeneity caused by incomplete lineage sorting (ILS) and gene tree estimation error stemming from the short internal branches. However, the relative influence of these factors in species tree inference is not well understood. Using anchored hybrid enrichment, we generated a data set including 423 single-copy loci from 64 taxa representing 39 families to infer the species tree of the flowering plant order Malpighiales. This order includes 9 of the top 10 most unstable nodes in angiosperms, which have been hypothesized to arise from the rapid radiation during the Cretaceous. Here, we show that coalescent-based methods do not resolve the backbone of Malpighiales and concatenation methods yield inconsistent estimations, providing evidence that gene tree heterogeneity is high in this clade. Despite high levels of ILS and gene tree estimation error, our simulations demonstrate that these two factors alone are insufficient to explain the lack of resolution in this order. To explore this further, we examined triplet frequencies among empirical gene trees and discovered some of them deviated significantly from those attributed to ILS and estimation error, suggesting gene flow as an additional and previously unappreciated phenomenon promoting gene tree variation in Malpighiales. Finally, we applied a novel method to quantify the relative contribution of these three primary sources of gene tree heterogeneity and demonstrated that ILS, gene tree estimation error, and gene flow contributed to 10.0$\%$, 34.8$\%$, and 21.4$\%$ of the variation, respectively. Together, our results suggest that a perfect storm of factors likely influence this lack of resolution, and further indicate that recalcitrant phylogenetic relationships like the backbone of Malpighiales may be better represented as phylogenetic networks. Thus, reducing such groups solely to existing models that adhere strictly to bifurcating trees greatly oversimplifies reality, and obscures our ability to more clearly discern the process of evolution. [Coalescent; concatenation; flanking region; hybrid enrichment, introgression; phylogenomics; rapid radiation, triplet frequency.].

Asymmetric Distribution of Gene Trees Can Arise under Purifying Selection If Differences in Population Size Exist

Effects of discordance between species and gene trees on phylogenetic diversity conservation

The Structure of Genealogies in the Presence of Purifying Selection: A "Fitness-Class Coalescent"

Implications of gene tree heterogeneity on downstream phylogenetic analyses: A case study employing the Fair Proportion index

Evolutionary Rate Variation among Lineages in Gene Trees has a Negative Impact on Species-Tree Inference

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

Genealogies under purifying selection

A tale of too many trees: a conundrum for phylogenetic regression

Probabilities of unranked and ranked anomaly zones under birth-death models

Anomaly zones for uniformly sampled gene trees under the gene duplication and loss model

Natural selection and repeated patterns of molecular evolution following allopatric divergence.

The Impact of Missing Data on Species Tree Estimation.

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

The Perfect Storm: Gene Tree Estimation Error, Incomplete Lineage Sorting, and Ancient Gene Flow Explain the Most Recalcitrant Ancient Angiosperm Clade, Malpighiales.

Complex Models of Sequence Evolution Improve Fit, but not Gene Tree Discordance, for Tetrapod Mitogenomes

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

Distributions of topological tree metrics between a species tree and a gene tree

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

Maximum-frequency gene tree: a simplified genome-scale approach to overcoming incongruence in molecular phylogenies

Genetic Diversity and the Efficacy of Purifying Selection across Plant and Animal Species

Pervasive incomplete lineage sorting illuminates speciation and selection in primates