Complete chloroplast of four Sanicula taxa (Apiaceae) endemic to China: lights into genome structure, comparative analysis, and phylogenetic relationships

Abstract Background The genus Sanicula comprises ca. 45 taxa, widely distributed from East Asia to North America, which is a taxonomically difficult genus with high medicinal value in Apiaceae. The systematic classification of the genus has been controversial for a long time due to varied characters in key morphological traits. China is one of the most important distributed centers, with ca. 18 species and two varieties. At present, chloroplast genomes are generally considered to be conservative and play an important role in evolutionary relationship study. To investigate the plastome evolution and phylogenetic relationships of Chinese Sanicula, we comprehensively analyzed the structural characteristics of 13 Chinese Sanicula chloroplasts and reconstructed their phylogenetic relationships. Results In present study, four newly complete chloroplast genome of Sanicula taxa by using Illumina sequencing were reported, with the typical quadripartite structure and 155,396–155,757 bp in size. They encoded 126 genes, including 86 protein-coding genes, 32 tRNA genes and 8 rRNA genes. Genome structure, distributions of SDRs and SSRs, gene content, among Sanicula taxa, were similar. The nineteen intergenic spacers regions, including atpH-atpI, ndhC-trnM, petB-petD, petD-rpoA, petN-psbM, psaJ-rpl33, rbcL-accD, rpoB-trnC, rps16-trnQ, trnE-psbD, trnF-ndhJ, trnH-psbA, trnN-ndhF, trnS-psbZ, trnS-trnR, trnT-trnF, trnV-rps12, ycf3-trnS and ycf4-cemA, and one coding region (ycf1 gene) were the most variable. Results of maximum likelihood analysis based on 79 unique coding genes of 13 Chinese Sanicula samples and two Eryngium (Apiaceae-Saniculoideae) species as outgroup taxa revealed that they divided into four subclades belonged to two clades, and one subclade was consistent with previously traditional Sanicula section of its system. The current classification based on morphology at sect. Sanicla and Sect. Tuberculatae in Chinese Sanicula was not supported by analysis of cp genome phylogeny. Conclusions The chloroplast genome structure of Sanicula was similar to other angiosperms and possessed the typical quadripartite structure with the conserved genome arrangement and gene features. However, their size varied owing to expansion/contraction of IR/SC boundaries. The variation of non-coding regions was larger than coding regions of the chloroplast genome. Phylogenetic analysis within these Chinese Sanicula were determined using the 79 unique coding genes. These results could provide important data for systematic, phylogenomic and evolutionary research in the genus for the future studies.