Abstract
Primula, well known for its heterostyly, is the largest genus in the family Primulaceae with more than 500 species. The considerable species number has introduced a huge challenge for taxonomy. The phylogenetic relationships among Primula still maintain unresolved due to frequent hybridization and introgression between closely related species. In this study, we sequenced and assembled the complete chloroplast genomes of Primula wilsonii Dunn, which is a PSESP (plant species with extremely small populations), using Illumina sequencing and compared its genomic sequences with those of four related Primula species. The chloroplast genomes of Primula species were similar in the basic structure, gene order, and GC content. The detected 38 SSRs (simple sequence repeats) loci and 17 hypervariable regions had many similarities in P. wilsonii, P. anisodora, P. miyabeana, and P. poissonii, but showed a significant difference compared with those in P. secundiflora. Slight variations were observed among Primula chloroplast genomes, in consideration of the relatively stable patterns of IR (inverted repeats) contraction and expansion. Phylogenetic analysis based on chloroplast genomes and protein-coding genes confirmed three major clades in Chinese Primula, but the infrageneric sections were not in accordance with morphological traits. The P. poissonii complex was confirmed here and P. anisodora was the most closely related species to P. wilsonii. Overall, the chloroplast genome sequences provided useful genetic and evolutionary information for phylogeny and population genetics on Chinese Primula species.