Phylogenomic analyses based on genome-skimming data reveal cyto-nuclear discordance in the evolutionary history of Cotoneaster (Rosaceae)

[Display omitted] •This is the first phylogenomic study of Cotoneaster with extensive taxon sampling.•We assembled 69 complete chloroplast genomes and extracted sequences of 204 low-copy nuclear genes.•Both complete chloroplast genome and nuclear gene sequences yielded robust phylogenetic trees of Cotoneaster, although discordance existed in some nodes.•The species tree based on nuclear gene sequences largely agrees with the traditional classification proposed based on morphological evidence.•Strong cyto-nuclear discordance discovered may be caused by incomplete lineage sorting (ILS) and/or hybridization. As a consequence of hybridization, polyploidization, and apomixis, the genus Cotoneaster (Rosaceae) represents one of the most complicated and controversial lineages in Rosaceae, with ca. 370 species which have been classified into two subgenera and several sections, and is notorious for its taxonomic difficulty. The infrageneric relationships and taxonomy of Cotoneaster have remained poorly understood. Previous studies have focused mainly on natural hybridization involving only several species, and phylogeny based on very limited markers. In the present study, the sequences of complete chloroplast genomes and 204 low-copy nuclear genes of 72 accessions, representing 69 species as ingroups, were used to conduct the most comprehensive phylogenetic analysis so far for Cotoneaster. Based on the sequences of complete chloroplast genomes and many nuclear genes, our analyses yield two robust phylogenetic trees respectively. Chloroplast genome and nuclear data confidently resolved relationships of this genus into two major clades which largely supported current classification based on morphological evidence. However, conflicts between the chloroplast genome and low-copy nuclear phylogenies were observed in both the species level and clade level. Cyto-nuclear discordance in the phylogeny could be caused by frequent hybridization events and incomplete sorting lineage (ILS). In addition, our divergence-time analysis revealed an evolutionary radiation of the genus from late Miocene to date.