Complete Chloroplast Genomes and the Phylogenetic Analysis of Three Native Species of Paeoniaceae from the Sino-Himalayan Flora Subkingdom

Paeonia delavayi var. lutea, Paeonia delavayi var. angustiloba, and Paeonia ludlowii are Chinese endemics that belong to the Paeoniaceae family and have vital medicinal and ornamental value. It is often difficult to classify Paeoniaceae plants based on their morphological characteristics, and the limited genomic information has strongly hindered molecular evolution and phylogenetic studies of Paeoniaceae. In this study, we sequenced, assembled, and annotated the chloroplast genomes of P. delavayi var. lutea, P. delavayi var. angustiloba, and P. ludlowii. The chloroplast genomes of these strains were comparatively analyzed, and their phylogenetic relationships and divergence times were inferred. These three chloroplast genomes exhibited a typical quadripartite structure and were 152,687–152,759 bp in length. Each genome contains 126–132 genes, including 81–87 protein-coding genes, 37 transfer RNAs, and 8 ribosomal RNAs. In addition, the genomes had 61–64 SSRs, with mononucleotide repeats being the most abundant. The codon bias patterns of the three species tend to use codons ending in A/U. Six regions of high variability were identified (psbK-psbL, trnG-UCC, petN-psbM, psbC, rps8-rpl14, and ycf1) that can be used as DNA molecular markers for phylogenetic and taxonomic analysis. The Ka/Ks ratio indicates positive selection for the rps18 gene associated with self-replication. The phylogenetic analysis of 99 chloroplast genomes from Saxifragales clarified the phylogenetic relationships of Paeoniaceae and revealed that P. delavayi var. lutea, P. delavayi var. angustiloba, and P. ludlowii are monophyletic groups and sisters to P. delavayi. Divergence time estimation revealed two evolutionary divergences of Paeoniaceae species in the early Oligocene and Miocene. Afterward, they underwent rapid adaptive radiation from the Pliocene to the early Pleistocene when P. delavayi var. lutea, P. delavayi var. angustiloba, and P. ludlowii formed. The results of this study enrich the chloroplast genomic information of Paeoniaceae and reveal new insights into the phylogeny of Paeoniaceae.