Phylogenomic conflict resulting from ancient introgression following species diversification in Stewartia s.l. (Theaceae).

Phylogenetic analyses using diverse datasets can yield conflicting inference of evolutionary history. Phylogenetic conflicts observed in both animal and plant systems have often been explained by two competing (but not mutually exclusive) hypotheses, i.e., hybridization vs. incomplete lineage sorting (ILS). The likelihood of either process contributing to phylogenetic conflict in a given group is context-dependent, involving attributes of life history, distribution, and phylogeny, among others. Here we explore phylogenetic conflict in Stewartia s.l., a genus with ca. 20 species of trees and shrubs from the tea family (Theaceae) disjunctly distributed between eastern Asia (EAS) and eastern North America (ENA). We use both restriction-site associated DNA sequencing (RAD-seq) and complete plastome sequence data to reconstruct the phylogeny of the group using concatenation and coalescence approaches. Our results indicate strong conflicts between the topologies reconstructed using nuclear and plastid data. Four-taxon D-statistic (ABBA-BABA) tests detected prevailing signals of introgression. Bayesian Analysis of Macro-evolutionary Mixtures (BAMM) inferred that species diversification occurred in the middle to late Miocene. Ancestral range reconstructions indicated co-distribution of ancestral species (represented by internal nodes) for both the Hartia clade (in southern China) and the EAS Stewartia s.s. clade (Japan Archipelago and the Yangtze Valley of China). The latter clade experienced multiple events of dispersal and vicariance during its diversification history. Ancient introgressive hybridization following species diversification in the mid- to late-Miocene likely caused diverging histories in the nuclear and plastid genomes, leading to phylogenetic conflict in Stewartia s.l. Our study indicates that species diversification driven by both the intensification of the East Asian summer monsoon since the late Miocene and reduced risks of extinction due to frequent dispersal possibly via East China Sea Land Bridge impacted the anomalous species richness between EAS and ENA. Our study highlights the importance of using data from different genomes while reconstructing deep and shallow phylogenies of organisms.