Abstract
Guizhou Province is one of the most important karst regions of southwest China, with 22 Primulina species. These species are highly habitat-specialized and dependent on the soils of the karst region, and many inter-species classifications are unclear. Therefore, studying the chloroplast genomes and estimating the divergence times of there species can not only provide a better understanding of interspecific relationships but also help to know the species speciation and divergence in the karst environment. So, we sequenced and assembled the chloroplast genomes of nine Primulina species (including six endemic species of Guizhou) and conducted chloroplast genome comparison analysis and phylogenetic study. The chloroplast genome structures of the nine Primulina species were quadripartite with total lengths of 152,869-153,364 base pairs (bp) and GC content of 37.55-37.64%. There are 132 functional genes annotated, respectively. A total of 375 simple sequence repeats and 375 interspersed nuclear elements were identified. The 30 highly preferred codons identified were used at similar frequencies in different species, respectively. A phylogenetic tree constructed on the basis of the 38 chloroplast genomes showed that Primulina species form a monophyletic group. Eleven mutational hotspot regions that could serve as potential molecular markers were identified, of which two regions near the 3' and 5' ends of the ycf1 gene were of appropriate size and can serve as molecular markers for phylogenetic studies of Primulina. The results of molecular clock analyses indicate that the three major branches of Primulina begin to diverge in the Miocene, and the number of species proliferated in the Pliocene and Pleistocene. Most of the species of Primulina in Guizhou Province were formed in the Pleistocene and rapidly diverged within a short period of time. This research study enriches the genetic resource information of Primulina and deepens the understanding of the phylogenetic relationships of the genus.