Comprehensive Phylogenomic Inference of Eucalypts Reveals Taxonomic Limits through Gene Tree Concordance

Speciation can be a slow process for long-lived organisms with large populations due to persistent incomplete lineage sorting, while hybridisation can further obscure species boundaries through admixture. In eucalypts, there is tremendous morphological and environmental variation with almost a thousand described species, subspecies, and hybrid taxa falling into multiple subgenera and intermediate sections. Multilocus molecular studies have generated concatenated species trees supporting deeper divergence between the Corymbia + Blakella + Angophora clade with the rest of eucalypts. However, the proportion of gene trees that recover the concatenated tree topology, including the taxonomic relationships among individual subgenera and sections, remains underappreciated. In this study, we assessed the predictive power of current eucalypt taxonomy - represented as a concatenated tree - in explaining topological variation across gene trees. To do this, we first extracted 1,187 BUSCO loci from short-read data of 701 samples representing roughly 500 described species using CAPTUS. We then inferred individual gene trees and a concatenated species tree using IQ-TREE2. Finally, we calculated the concordance factors for every described subgenus and section on the concatenated tree, and developed a taxonomic consistency test that tracks the taxonomic placement of each sample across gene trees. We show that, even though most taxonomic groups are monophyletic on the concatenated tree with 100% bootstrap support, several groups exhibit very low concordance, indicating that most gene trees do not recover their relationships as represented by the concatenated tree. These results highlight important considerations for current eucalypt taxonomic classification and provide a phylogenomic inference pipeline that can be applied not only to future eucalypt studies, but also to other evolutionarily complex groups of taxa.