Complex evolutionary history of Rosales mediated by extensive incomplete lineage sorting and hybridization

The angiosperm order Rosales still represents a major challenge for phylogenetic reconstruction. Although its circumscription is now well-defined, phylogenetic relationships among families are still uncertain. Here, we used nuclear, plastid, and mitochondrial genomic data from 33 species representing all nine families to further clarify interfamilial relationships and the group's evolutionary history. We detected significant phylogenetic conflict among the three datasets. Further analyses at the nuclear level identified incomplete lineage sorting (ILS) as the main cause of unstable phylogenetic positions among families. The discordant placements of Rhamnaceae and Elaeagnaceae based on plastid and mitochondrial data are caused by ancient hybridization events, potentially involving differences in organellar inheritance. Our molecular dating confirms earlier suggestions that the ancient rapid diversification of the three Rosaceae subfamilies could be the main reason for the difficulties in resolving their phylogenetic relationships. Our findings provide new insights into the interfamilial relationships of Rosales and demonstrate that the evolutionary history of this order was shaped by ancient and rapid radiation as well as extensive ILS and reticulate evolution. They also suggest that previous attempts to clarify interfamilial relationships in this order were hampered by combining nuclear and organellar sequence data, leading to inconsistent topologies observed across earlier studies.