Abstract
Cherries (Prunus Subgenus Cerasus) have economic value and ecological significance, yet their phylogeny, geographic origin, timing, and dispersal patterns remain challenging to understand. To fill this gap, we conducted a comprehensive analysis of the complete chloroplast genomes of 54 subg. Cerasus individuals, along with 36 additional genomes from the NCBI database, resulting in a total of 90 genomes for comparative analysis. The chloroplast genomes of subg. Cerasus exhibited varying sizes and consisted of 129 genes, including protein-coding, transfer RNA, and ribosomIal RNA genes. Genomic variation was investigated through InDels and SNPs, showcasing distribution patterns and impact levels. A comparative analysis of chloroplast genome boundaries highlighted variations in inverted repeat (IR) regions among Cerasus and other Prunus species. Phylogeny based on whole-chloroplast genome sequences supported the division of Prunus into three subgenera, I subg. Padus, II subg. Prunus and III subg. Cerasus. The subg. Cerasus was subdivided into seven lineages (IIIa to IIIg), which matched roughly to taxonomic sections. The subg. Padus first diverged 51.42 Mya, followed by the separation of subg. Cerasus from subg. Prunus 39.27 Mya. The subg. Cerasus started diversification at 15.01 Mya, coinciding with geological and climatic changes, including the uplift of the Qinghai-Tibet Plateau and global cooling. The Himalayans were the refuge of cherries, from which a few species reached Europe through westward migration and another species reached North America through northeastward migration. The mainstage of cherry evolution was on the Qing-Tibet Plateau and later East China and Japan as well. These findings strengthen our understanding of the evolution of cherry and provide valuable insights into the conservation and sustainable utilization of cherry's genetic resources.