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.].

Reconstructing the Basal Angiosperm Phylogeny: Evaluating Information Content of Mitochondrial Genes

Angiosperm Phylogeny Inferred from Sequences of Four Mitochondrial Genes

Mitochondrial Genes from 18 Angiosperms Fill Sampling Gaps for Phylogenomic Inferences of the Early Diversification of Flowering Plants

Resolution of Deep Angiosperm Phylogeny Using Conserved Nuclear Genes and Estimates of Early Divergence Times

Angiosperm Phylogeny: 17 Genes, 640 Taxa

Expanded Taxonomic Sampling Coupled with Gene Genealogy Interrogation Provides Unambiguous Resolution for the Evolutionary Root of Angiosperms

Phylogenomic Insights into Deep Phylogeny of Angiosperms Based on Broad Nuclear Gene Sampling

Coalescent Versus Concatenation Methods and the Placement of Amborella As Sister to Water Lilies.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Molecular phylogeny of extant gymnosperms and seed plant evolution: analysis of nuclear 18S rRNA sequences.

The Evolutionary Root of Flowering Plants

Genome-scale Angiosperm Phylogenies Based on Nuclear,plastome,and Mitochondrial Datasets

Analysis of 81 genes from 64 plastid genomes resolves relationships in angiosperms and identifies genome-scale evolutionary patterns

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

Phylogenomics and Coalescent Analyses Resolve Extant Seed Plant Relationships

Phylogenomic and Syntenic Data Demonstrate Complex Evolutionary Processes in Early Radiation of the Rosids

Modeling compositional heterogeneity resolves deep phylogeny of flowering plants

Highly Conserved Low-Copy Nuclear Genes As Effective Markers for Phylogenetic Analyses in Angiosperms

The chloroplast genome of Nymphaea alba: whole-genome analyses and the problem of identifying the most basal angiosperm

Early Diversification and Karyotype Evolution of Flowering Plants

Novel phylogeny of angiosperms inferred from whole-genome microsynteny analysis