Abstract
Pulsatilla (Ranunculaceae) consists of about 40 species, and many of them have horticultural and/or medicinal value. However, it is difficult to recognize and identify wild Pulsatilla species. Universal molecular markers have been used to identify these species, but insufficient phylogenetic signal was available. Here, we compared the complete chloroplast genomes of seven Pulsatilla species. The chloroplast genomes of Pulsatilla were very similar and their length ranges from 161,501 to 162,669 bp. Eight highly variable regions and potential sources of molecular markers such as simple sequence repeats, large repeat sequences, and single nucleotide polymorphisms were identified, which are valuable for studies of infra- and inter-specific genetic diversity. The SNP number differentiating any two Pulsatilla chloroplast genomes ranged from 112 to 1214, and provided sufficient data for species delimitation. Phylogenetic trees based on different data sets were consistent with one another, with the IR, SSC regions and the barcode combination rbcL + matK + trnH-psbA produced slightly different results. Phylogenetic relationships within Pulsatilla were certainly resolved using the complete cp genome sequences. Overall, this study provides plentiful chloroplast genomic resources, which will be helpful to identify members of this taxonomically challenging group in further investigation.
Similar content being viewed by others
Introduction
DNA barcoding is an effective tool to identify many plant species rapidly and accurately1,2,3,4. However, there is no single universal barcode that can be successfully used to identify all plants to the species level5. Consequently, two alternative strategies have been proposed to distinguish among plant species: the first one is the use of complete chloroplast genomes6, named ‘super-barcoding’, and the second one is an approach that involves searching for mutational hotspots7, or using comparative plastid analyses to find loci with suitable species-level divergence8,9. Analyses of entire chloroplast genome sequences provide an effective way to develop both of these strategies.
In most angiosperms, the chloroplast genomes are inherited maternally and have a consistent structure, including two inverted repeats (IR), one large (LSC) and one small (SSC) single copy region. The chloroplast genome always contains 110–130 genes that exhibit a range of levels of polymorphism4,5. Thus, chloroplast genome sequence data are extremely valuable for studies of plant population genetics, phylogeny reconstruction, species identification, and genome evolution2,4,5,6.
The Ranunculaceae is a large family, which includes approximately 59 genera and 2500 species. Many plants of Ranunculaceae are pharmaceutically important10. The genus Pulsatilla Adans. consists of about 40 species which are distributed in temperate subarctic and mountainous areas of the Northern Hemisphere10. There are always long, soft hairs covering plants of Pulsatilla species. Most of the flowers of Pulsatilla are large and showy, and therefore the genus has horticultural importance11,12. The flowers are solitary and bisexual. In one flower, there are always six tepals, numerous stamens and carpels, with the outermost stamens resembling degenerated petals, excluding P. kostyczewii11,13,14,15.
In China, there are eleven species of Pulsatilla. Some species of Pulsatilla have been used in traditional Chinese medicine for many years, such as for “detoxification” or “blood-cooling”, because Pulsatilla species contain numerous secondary metabolites, including phytosterols, triterpenoid saponins and anthocyanins16. At the same time, all members of Pulsatilla produced the lactone protoanemonin17,18,19.
In Europe, some species of Pulsatilla are rare, endangered and endemic. Those taxa are protected due to their small populations and disappearing localities, and those species have been placed on the Red Lists of Endangered Species20.
Taxonomically, Pulsatilla is an especially complex and challenging group. In all treatments published before, three subgenera have been recognized: subgenus Kostyczewianae (only one species), subgenus Preonanthus, and the largest subgenus Pulsatilla. However, the intragenic morphological variability of Pulsatilla was especially complicated12. The recognition and identification of wild Pulsatilla species is particularly difficult based on traditional approaches.
Molecular markers are significant to explore the phylogenetic relationships of the genus Pulsatilla. Phylogenetic relationships between Pulsatilla and closely related genera have been dedicated during the past years21,22,23,24,25. Previous studies have attempted to identify these species among Pulsatilla with universal molecular markers, but the species resolution was relatively low15.
In this study, we present seven complete cp genomes from two subgenera of Pulsatilla obtained through next-generation sequencing (NGS) and genomic comparative analyses with four previously published cp genome sequences of Pulsatilla from NCBI, with Anemoclema glaucifolium as the outgroup. We identify microsatellites (SSRs), larger repeat sequences, and highly variable regions, with the aim of developing DNA barcodes and testing the feasibility of phylogenetic analyses of Pulsatilla using the chloroplast genome.
Results and discussion
Chloroplast genome features
We have obtained 1.95 Gb of average NGS clean data for each species, with minimum and maximum values of 1.14 Gb (P. dahurica), and 3.56 Gb (Pulsatilla alpina), respectively. The read number for each species ranged from 6,468,944 (P. dahurica) to 15,816,765 (P. alpina). The average length of the reads was 150 bp on the Illumina Sequencing System. The seven new Pulsatilla cp genomes ranged from 161,501 bp (P. grandis) to 162,669 bp (P. alpina) in length and 151.5 × to 503.4 × coverage. These seven novel Pulsatilla cp genome sequences were submitted to GenBank (Tables 1, 2). Their quadripartite structure is similar to the majority of cp genomes of land plants, which are composed of a pair of IRs (31,184–31,416 bp), separated by the LSC (81,615–82,149 bp) and SSC (of 17,431–17,908 bp) regions26,27 (Fig. 1; Table 2). Previous studies of other angiosperm groups have found that chloroplast genomes are conserved28 or highly polymorphic29,30. These genomes which we reported are highly conserved in gene order, gene content and intron number, which is in accordance with the results from many other taxa26,27. However, in some taxa, e.g. Amorphophallus of Araceae, some genes were lost31, and in others, e.g. Pelargonium32, the structure and gene order diverges from what is reported here and in most other angiosperms. The cp genomes of P. alpina, P. grandis, P. hirsutissima, P. ludoviciana, P. multifida and P. occidentalis had the same GC content of 37.6%, while P. dahurica had a subtle difference (37.5%) compared with the others.
Chloroplast genome comparison
In most angiosperms, the IR regions of cp genomes of angiosperms are highly conserved, but the expansion and contraction of IR region boundaries are ever present33,34. At the same time, several lineages of land plant chloroplast genomes show great structural rearrangement, even loss of IR regions or some gene families35. The expansion and contraction in IRs are significant evolutionary events, because they can change gene content and chloroplast genome size30,36. Expansion of the IRs has been reported in Araceae36,37. Sometimes, the size of LSC increases and that of SSC decreases, becoming only 7000 bp in Pothos38. At the same time, a linear chloroplast genome was also reported in some groups, e.g. maize35,39. Expansion and contraction of the IR regions can also lead to duplication of certain genes or conversion of duplicate genes to single copy, respectively30,36. Changes in the size of the IRs can also cause rearrangement of the genes in the SSC as recently observed in Zantedeschia36.
The Pulsatilla chloroplast genomes were compared to previously published data and showed typical Anemoneae (Ranunculaceae) genome structure37,39. As reported for Anemoclema, Anemone, Clematis and Hepatica, the IR regions of genus Pulsatilla are roughly 4.4 kb longer than those of other genera of the family Ranunculaceae, such as Aconitum, Coptis, Thalictrum, Megaleranthis, Ranunculus, and Trollius37,39. The gene orders located within the IR-SSC and IR-LSC boundaries are similar among tribe Anemoneae but different from those of other genera of Ranunculaceae (Fig. 2, Fig. S1).
We compared the IR/SC boundary regions of Pulsatilla, and the junction positions are very similar and conserved within genus Pulsatilla. In the four boundary regions (LSC, IRa, SSC, IRb) of seven Pulsatilla cp genomes, the LSC/IRa and IRb/LSC border was in the intergenic region, and the adjacent genes is rps36, rps8 and rps4, respectively. The genes ycf1 andψycf1 have crossed the SSC/IRb and IRa/SSC boundary, respectively, which was also found in Monsteroideae (Araceae)28. The pseudogene ycf1 has been found in other groups30,36. The IR regions were highly conserved, with nucleotide diversity values in those regions less than 2%.
Sequence divergence
Multiple alignments of plastid genomes were performed to investigate levels of genome divergence. Based on MAFFT analysis, there are three inversions in LSC of Pulsatilla, same as the tribe Anemoneae (Fig. 1, Fig. S1)37. The mVISTA analysis has revealed high sequence similarity across the coding region and there exists more variability in non-coding regions. Sequence identity among the seven species was 96.68–98.66%. The number of nucleotide substitutions and sequence distance (Pi) were the highest (1214, 0.0063) between P. alpine and P. dahurica, with the lowest (112, 0.0005) between P. vernalis and P. patens (Fig. 3; Table 3).
Identification of highly variable regions
Chloroplast genome markers, especially several universal chloroplast regions, have been widely used in plant systematics and identification at multiple taxonomic levels. Highly suitable polymorphic chloroplast loci have been identified and designed as unique markers in different groups28,40. However, relationships within the genus Pulsatilla have not been well resolved because of the low polymorphism of these universal markers15. In order to facilitate identification of closely related species of Pulsatilla, we sought to identify highly variable regions of the chloroplast genome, as previously described9,27,41,42,43,44. As a result, we identified nine divergent hotspot regions, including six intergenic spacer regions (rps4-rps16, rps16-matK, ndhC-trnV, psbE-petL, ndhD-ccsA, ccsA-ndhF) and four protein-coding regions (ycf1, ndhF, ndhI) (Fig. 4; Table 4). Most commonly employed loci, e.g. trnL-trnF, trnH-psbA were not selected in our finding. The nine highly variable regions included 684 variable sites, including 181 indels. However, these indels are not suitable for the phylogenetic inference because Maximum likelihood model used only substitutions not indels28. Their nucleotide diversity values ranged from 0.00802 to 0.02212. The region of ccsA-ndhF showed the highest variability, the next most variable regions were rps4-rps16, ndhC-trnV, and psbE-petL. The diversity level of two protein-coding regions (ycf1, ndhF) was the lowest.
Among the nine divergent hotspot regions, the ndhI is difficult to align. There are large numbers of indels in ndhI and the intergenic spacer between ndhI and ndhG, these regions were not considered suitable for the phylogenetic inference of the Pulsatilla. Thus, we selected eight regions, four (rps4-rps16, rps16-matK, ndhC-trnV, psbE-petL) in the LSC and four (ndhD-ccsA, ccsA-ndhF, ycf1, ndhF) in the SSC, with relatively high variability as potential molecular markers for the study of species identification and phylogeny in Pulsatilla. Five hotpots were found in chloroplast genome of Veroniceae (Plantaginaceae), and two universal marker, trnH-psbA and matK were identified, respectively45,46. Ten highly variable regions were selected as potential molecular markers for Fritillaria, including ycf144,47, which was also selected in this study. Sequences of these variable regions founded in this study could be regarded as potential molecular markers for species identification and evolutionary studies and have been shown to be valuable for studies in other groups (e.g., Fritillaria)44.
SSRs and large repeat sequences
Oligonucleotide repeats play an important role for generating indels, inversion and substitutions29. Repeat sequences in the chloroplast genome could provide valuable information for understanding not only the sequence divergence but the evolutionary history of the plants48,49,50. We have detected five types of large repeats (forward, reverse, palindromic, complement and tandem repeats) in the seven Pulsatilla cp genomes. Among them, the most common repeat types are forward and palindromic repeats, followed by reverse repeats, and only little complement repeats were found in Pulsatilla cp genomes (Fig. 5A). Most of the repeats were short, ranging from 30–49 bp (Fig. 5B).
We also identified multiple microsatellite repeats, also known as simple sequence repeats (SSR) or short tandem repeats (STR)49. Due to their codominant inheritance and high variability, SSRs are robust and effective markers for species identification and population genetic analyses49,50,51,52,53. Most of the mononucleotide repeats were composed of A/T. The other microsatellites types were also dominated by AT/TA, with very little G/C (Fig. 5C). In this study, plentiful microsatellite loci were found through the comparative analysis of Pulsatilla cp genome sequences. In total, we detected six types of microsatellite (mononucleotide, dinucleotide, trinucleotide, tetranucleotide, pentanucleotide and hexanucleotide repeats) based on the comparison of seven Pulsatilla cp genomes (Fig. 5D). Each Pulsatilla cp genome had 69–87 microsatellites. The lengths of repeat motifs of these microsatellites ranged from 10 to 21 bp (Fig. 5E). Among the four structural regions in the cp genomes, most of the repeats and microsatellites were distributed in LSC, followed by SSC, and fewest in IRa/IRb (Fig. 5F), which were also reported in other studies in angiosperms29,30. These SSRs and repeat sequences are uncorrelated with genome size and phylogenetic position of the species36, but will provide important information for further studies of phylogenetic reconstruction and infra- and inter-specific genetic diversity54,55.
Phylogenetic analyses
Chloroplast genomes have been widely used and have made significant contributions to phylogeny reconstruction at different taxonomic levels in plants7,9,24,35. To better clarify the evolutionary relationships within Pulsatilla, we used each data set to construct phylogenetic trees using the ML analytical methods. We also construct phylogenetic trees with those eight highly variable regions using the ML, MP analytical methods. All tree topology structures were identical. Therefore, here we presented the phylogenetic studies using the ML tree with the support values from the MP analyses recorded at the corresponding nodes (Fig. 6). The phylogenetic tree based on all data sets (except the IR and SSC regions) from the complete plastid genome sequences yielded the same topology. The phylogenetic tree based on chloroplast genome differed from that of the DNA barcode combination rbcL + matK + trnH-psbA, but with higher support values. The phylogenetic trees based on data from complete plastid genome sequences showed that the species of Pulsatilla formed a monophyletic group which in turn includes two strongly supported (bootstrap = 100) clades. One clade comprised P. alpina and P. occidentalis, members of subg. Preonanthus. The other comprised two subclades: (1) members of P. hirsutissima, P. ludoviciana, P. multifidi, P. patens and P. vernalis, and (2) species of P. chinensis, P. dahurica, P. grandis and P. pratensis. All the species of the two subclades are members of the subg. Pulsatilla. These results were congruent with our former results based on universal markers15.
In phylogenetic analyses, compared to the combination of barcodes, the full chloroplast genome sequence data formed distinct clades with high bootstrap support, improving the inadequate resolution of barcodes combination. The LSC regions and coding regions have the same topology structures with robust support. However, sequencing of genomic DNA is still expensive. It is necessary to utilize variation within chloroplast regions for rapid species-specific assay5,9,31,33,42. Here we found that phylogenetic inference based on highly variable regions yielded a tree with the same topology as the one recovered based on complete chloroplast genome sequences, demonstrating the high utility of hotspots of variability for species identification and phylogenetic analysis. More samples and laboratory works are needed in the future to increase the number of these variable regions available for study.
Conclusions
In this study, we generated complete chloroplast genomes of seven species of Pulsatilla and compared them to four previously published cp genome sequences of Pulsatilla. Chloroplast genomes of Pulsatilla share many features with those of other angiosperms. Informative differences between cp genomes of Pulsatilla were most evident in inversions in the large single copy region and expansion of the inverted repeat region. We identified multiple potentially valuable genetic markers, including large repeat sequences, numerous SSRs, and eight highly variable regions. Genetic markers have provided a reference for the improvement of plants fingerprints and the identification of similar Pulsatilla. In addition, it will be better to construct the phylogenetic relationships of Pulsatilla species using these highly variable regions and genetic markers in the future. Totally, this study provides a basis for future studies of horticultural cultivation, conservation, population genetics, phyletic evolution, development of DNA barcodes, and diverse research in Pulsatilla.
Materials and methods
Plant material, DNA extraction and sequencing
Fresh leaves of P. dahurica were collected from Jilin province of China and dried with silica gel. Dry leaves of other six Pulsatilla species were taken from herbarium specimens. We extracted total genomic DNA with the DNeasy Plant mini kits (QIAGEN, Guangzhou, China). The genomic DNA was sequenced using the Illumina Miseq platform (Illumina, San Diego, CA, USA).
Chloroplast genome assembly and annotation
Whole chloroplast genome sequencing was done for the seven species of Pulsatilla. For each species, high-quality Illumina sequencing reads were assembled into scaffolds with de novo sequence assembly software Spades, SOAPdenovo and CLC Genomics Workbench v.6.5 (CLC Bio), respectively56. We checked the contigs against the reference genome of P. chinensis (MG001341), using BLAST (https://blast.ncbi.nlm.nih.gov/) and oriented aligned contigs according to the reference genome. We mapped all the raw reads back to assembled sequences to check the assembly and then constructed the complete cp genomes using Geneious v.9.057. We submitted all the newly sequence data in raw format (fastq) and obtained SRA accesions (Table 1).
Annotations of cp genome sequences were performed using Plastid Genome Annotator (https://github.com/quxiaojian/PGA)58 and adjusted in Geneious v.9.0. To verify the accuracy of the annotations, we also used GeSeq59 to annotate each chloroplast genome in this study. We checked every boundary of tRNAs using tRNAscan-SE v.2.060. The circular genome maps were generated in OGDRAW (https://ogdraw.mpimp-golm.mpg.de/)61.
Genome comparisons
We aligned the cp genomes of Pulsatilla using multiple alignment of MAFFT v762 and manually edited in Geneious v.9.0. The contraction and expansion of inverted repeat regions were also examined among the seven species (excluded P. chinensis, MG001341) of the genus Pulsatilla using Irscope63. Then, we performed multiple alignments of the eight genomes of Pulsatilla in the mVISTA program64 under Shuffle-LAGAN mode, with default parameters for other options, using the annotation genome of P. chinensis as a reference, with the aim of comparing and visualized the similarities and differences among different plastid genomes.
To analyse chloroplast genome organisation and gene arrangement, we perform the analyses of collinear blocks with Mauve v 2.3.165 plugin in Geneious v.9.0, including only one copy of the IR, assuming collinear genomes for the full alignment. Detailed gene inversions were identified by comparing the gene order of Pulsatilla samples and Anemoclema to Berberis.
To observe the plastid genome divergence and determine parsimony informative sites, we conducted sliding window analysis after alignment to determine the nucleotide diversity (Pi) of the cp genome using DnaSP v5, with 200 bp of step size and 600 bp window length66. We defined hotspots as those regions with a higher value of Pi. We computed the variable sites across the complete cp genomes and the sequence characteristics of hotspots by DnaSP v5.0.
Repeated sequences identification
We identified repeat sequences, including palindromic, reverse and forward repeats, using the online software REPuter (https://bibiserv.cebitec.uni-bielefeld.de/reputer), with the following settings: Hamming distance of 3 and minimum repeat size of 30 bp67. We used the online program Tandem Repeats Finder (https://tandem.bu.edu/trf/trf.html) to find the tandem repeat sequences, in which the similarity percentage of two repeat copies was at least 90% and the minimal repeat size was 10 bp. The alignment parameters for match, mismatch, and indels were set at 2, 7, and 7, respectively. We identified microsatellites (SSRs) by MISA68 with thresholds of 10, 5, 4, 3, 3, and 3 for mono-, di-, tri-, tetra-, penta-, and hexa-nucleotide, respectively.
Phylogenetic analyses
For the purpose of reconstructing the phylogenetic relationships, four published complete cp genome sequences from the genus Pulsatilla and Anemoclema were also included in our analyses. The monotypic genus Anemoclema (MH205609) was selected as the outgroup. Because molecular evolutionary rates among the different cp genome regions are diverse, analyses of phylogenetic relationships were performed based on the following seven datasets: (a) the complete cp genome sequences; (b) coding genes (CDS); (c) one inverted repeat (IR) region (IRb); (d) the large single copy region (LSC); (e) the small single copy region (SSC), (f) the consensus sequences of eight highly variable regions; and (g) the DNA barcodes combination (rbcL + matK + trnH-psbA). We applied Maximum Likelihood (ML) analysis for each of the seven datasets to construct tree-sets. Maximum Parsimony (MP) analyses were also applied for the consensus sequences of eight highly variable regions and the DNA barcodes combination.
We conducted ML analyses with RAxMLHPC2 v.8.0.969 on the Cyberinfrastructure for Phylogenetic Research (CIPRES) Science Gateway v.3.370. Then the analysis of 1000 rapid bootstrap replicates (-x) was followed by a search for the best-scoring ML tree in one program (-f a). The best-fit model for nucleotide and amino acid sequences were evaluated using jModelTest 271. We applied the GTR + G model to nucleotide data for both bootstrapping and best-tree searching phases, with other parameters as the default settings. We performed the maximum Parsimony (MP) analysis on PAUP* v.4.0b1072. All the characters were treated as unordered and equally weighted. The heuristic search specified 1000 random sequence addition replicates with TBR branch swapping, saving only 10 trees per replicate. We obtained the strict consensus tree from all the most-parsimonious trees (MPTs) detected during the search. We calculated bootstrap percentages (BP) from 10,000 rapid bootstrap replicates, each comprising 10 random sequence addition replicates and saving only one tree per replicate.
References
Hebert, P. D. N., Cywinska, A., Ball, S. L. & Waard, J. R. Biological identifications through DNA barcodes. Proc. R. Soc. Lond. B Biol. Sci. 270, 313–321 (2003).
Kress, W. J., Wurdack, K. J., Zimmer, E. A., Weigt, L. A. & Janzen, D. H. Use of DNA barcodes to identify flowering plants. Proc. Natl. Acad. Sci. U.S.A. 102, 8369–8374 (2005).
Kress, W. J. et al. Plant DNA barcodes and a community phylogeny of a tropical forest dynamics plot in Panama. Proc. Natl. Acad. Sci. U.S.A. 106, 18621–18626 (2009).
CBOL Plant Working Group. A DNA barcode for land plants. Proc. Natl. Acad. Sci. U.S.A. 106, 12794–12797 (2009).
Li, X. W. et al. Plant DNA barcoding: from gene to genome. Biol. Rev. Camb. Philos. Soc. 90, 157–166 (2015).
Hollingsworth, P. M., Graham, S. W. & Little, D. P. Choosing and using a plant DNA barcode. PLoS ONE 6, e19254 (2011).
Ahmed, I. et al. Identification of chloroplast genome loci suitable for highresolution phylogeographic studies of Colocasiaesculenta (L.) Schott (Araceae) and closely related taxa. Mol. Ecol. Resour. 13, 929–937 (2013).
Kuang, D. Y. et al. Complete chloroplast genome sequence of Magnolia kwangsiensis (Magnoliaceae): implication for DNA barcoding and population genetics. Genome 54, 663–673 (2011).
Dong, W., Liu, J., Yu, J., Wang, L. & Zhou, S. Highly variable chloroplast markers for evaluating plant phylogeny at low taxonomic levels and for DNA barcoding. PLoS ONE 7, e35071 (2012).
Tamura, M. Ranunculaceae. in Die natürlichen pflanzenfamilien, Aufl. II. 2nd ed. 17a IV (ed Hiepko, P.) (Duncker und Humblot, Berlin, 1995).
Wang, W. C. et al. Flora of China Vol. 6, 133–438 (Science Press, Beijing; Missouri Botanical Garden Press, St. Louis, 2001).
Grey-Wilson, C. Pasque-flowers. The Genus Pulsatilla (. The Charlotte Louise Press, Kenning Hall, 2014).
Ren, Y., Gu, T. Q. & Chang, H. L. Floral development of Dichocarpum, Thalictrum, and Aquilegia (Thalictroideae, Ranunculaceae). Plant Syst. Evol. 292, 203–213 (2011).
Ren, Y., Chang, H. L. & Endress, P. K. Floral development in Anemoneae (Ranunculaceae). Bot. J. Linn. Soc. 162, 77–100 (2015).
Li, Q. J. et al. Efficient identification of Pulsatilla (Ranunculaceae) using DNA barcodes and micro-morphological characters. Front. Plant Sci. 10, 1196 (2019).
China Pharmacopoeia Committee. Pharmacopoeia of the People’s Republic of China (Chemical Industry Press, Beijing, 2015).
Xu, Q. et al. Antitumor activity of Pulsatilla chinensis (Bunge) Regel saponins in human liver tumor 7402 cells in vitro and in vivo. Phytomedicine 19, 293–300 (2012).
Wang, X. W., Fan, F. G. & Cao, Q. Modified Pulsatilla decoction attenuates oxazolone-induced colitis in mice through suppression of inflammation and epithelial barrier disruption. Mol. Med. Rep. 14, 1173–1179 (2016).
Suh, S. Y. & An, W. G. Systems pharmacological approach of Pulsatillae radix on treating Crohn’s disease. Evid. Based Complement. Altern. Med. 2017, 1–21 (2017).
Szczecińska, M. & Sawicki, J. Genomic resources of three Pulsatilla species reveal evolutionary hotspots, species-specific sites and variable plastid structure in the family Ranunculaceae. Int. J. Mol. Sci. 16, 22258–22279 (2015).
Hoot, S. B., Reznicek, A. A. & Palmer, J. D. Phylogenetic relationship in Anemone (Ranunculaceae) based on morphology and chloroplast DNA. Syst. Bot. 19, 169–200 (1994).
Hoot, S. B. Phylogenetic relationships in Anemone (Ranunculaceae) based on DNA restriction site variation and morphology. In Systematics and Evolution of the Ranunculiflorae (eds Jensen, U. & Kadereit, J. W.) 295–300 (Springer, Vienn, 1995).
Szczecińska, M., Gabor, S., Katarzyna, W., Jakub, S. & Dusan, G. Genetic diversity and population structure of the rare and endangered plant species Pulsatillapatens (L.) Mill in east central Europe. PLoS ONE 11, e0151730 (2016).
Jiang, N. et al. Phylogenetic reassessment of tribe Anemoneae (Ranunculaceae): non-monophyly of Anemone s.l. revealed by plastid datasets. PLoS ONE 12, e0174792 (2017).
Sramkó, G., Laczkó, L., Volkova, P. A., Bateman, R. M. & Mlinarec, J. Evolutionary history of the Pasque-flowers (Pulsatilla, Ranunculaceae): molecular phylogenetics, systematics and rDNA evolution. Mol. Phylogenet. Evol. 135, 45–61 (2019).
Yu, X. Q., Yang, D., Guo, C. & Gao, L. M. Plant phylogenomics based on genome-partitioning strategies: progress and prospects. Plant Divers. 40, 158–164 (2018).
Tang, Y., Yukawa, T., Bateman, R. M., Jiang, H. & Peng, H. Phylogeny and classification of the East Asian Amitostigma alliance (Orchidaceae: Orchideae) based on six DNA markers. BMC Evol. Biol. 15, 96 (2015).
Henriquez, C. L. et al. Molecular evolution of chloroplast genomes in Monsteroideae (Araceae). Planta 251, 72 (2020).
Abdullah, et al. Correlations among oligonucleotide repeats, nucleotide substitutions, and insertion–deletion mutations in chloroplast genomes of plant family Malvaceae. J. Syst. Evol. https://doi.org/10.1111/jse.12585 (2020).
Abdullah, et al. Chloroplast genome of Hibiscus rosa-sinensis (Malvaceae): comparative analyses and identification of mutational hotspots. Genomics 112, 581–591 (2020).
Liu, E. et al. Comparative analysis of complete chloroplast genome sequences of four major Amorphophallus species. Sci. Rep. 9, 809 (2019).
Chumley, T. W. et al. The complete chloroplast genome sequence of Pelargonium × hortorum: organization and evolution of the largest and most highly rearranged chloroplast genome of land plants. Mol. Biol. Evol. 23, 2175–2190.
Ogihara, Y. et al. Structural features of a wheat plastome as revealed by complete sequencing of chloroplast DNA. Mol. Genet. Genomics 266, 740–746 (2002).
Raman, G. & Park, S. The complete chloroplast genome sequence of Ampelopsis: gene organization, comparative analysis, and phylogenetic relationships to other angiosperms. Front. Plant Sci. 7, 341 (2016).
Daniell, H., Lin, C. S., Yu, M. & Chang, W. J. Chloroplast genomes: diversity, evolution, and applications in genetic engineering. Genome Biol. 17, 134 (2016).
Henriquez, C. L. et al. Evolutionary dynamics in chloroplast genomes of subfamily Aroideae (Araceae). Genomics 112, 2349–2360 (2020).
Wang, W. & Messing, J. High-throughput sequencing of three Lemnoideae (duckweeds) chloroplast genomes from total DNA. PLoS ONE 6, e24670 (2011).
Adbulch, et al. Complete chloroplast genomes of Anthurium huixtlense and Pothos scandens (Pothoideae, Araceae): unique inverted repeat expansion and contraction affect rate of evolution. J. Mol. Evol. 88, 562–574 (2020).
Oldenburg, D. J. & Bendich, A. J. The linear plastid chromosomes of maize: terminal sequences, structures, and implications for DNA replication. Curr. Genet. 62, 431–442 (2016).
Sun, J. H. et al. Evolutionary and phylogenetic aspects of the chloroplast genome of Chaenomeles species. Sci. Rep. 10, 11466 (2020).
Liu, H. J. et al. Comparative analysis of complete chloroplast genomes of Anemoclema, Anemone, Pulsatilla, and Hepatica revealing structural variations among genera in tribe Anemoneae (Ranunculaceae). Front. Plant Sci. 9, 1097 (2018).
Meng, J. et al. Comparative analysis of the complete chloroplast genomes of four Aconitum medicinal species. Molecules 23, 1015 (2018).
Du, Y. P. et al. Complete chloroplast genome sequences of Lilium: insights into evolutionary dynamics and phylogenetic analyses. Sci. Rep. 7, 5751 (2017).
Bi, Y. et al. Chloroplast genomic resources for phylogeny and DNA barcoding: a case study on Fritillaria. Sci. Rep. 8, 1184 (2018).
Choi, K. S., Chung, M. G. & Park, S. The complete chloroplast genome sequences of three Veroniceae species (Plantaginaceae): comparative analysis and highly divergent regions. Front. Plant Sci. 7, 1–8 (2016).
Menezes, A. P. A. et al. Chloroplast genomes of Byrsonima species (Malpighiaceae): Comparative analysis and screening of high divergence sequences. Sci. Rep. 8, 1–12 (2018).
Li, Y., Zhang, Z., Yang, J. & Lv, G. Complete chloroplast genome of seven Fritillaria species, variable DNA markers identification and phylogenetic relationships within the genus. PLoS ONE 13, e0194613 (2018).
Grassi, F., Labra, M., Scienza, A. & Imazio, S. Chloroplast SSR markers to assess DNA diversity in wild and cultivated grapevines. Vitis 41, 157–158 (2002).
Chen, C., Zhou, P., Choi, Y. A., Huang, S. & Gmitter, F. G. Jr. Mining and characterizing microsatellites from citrus ESTs. Theor. Appl. Genet. 112, 1248–1257 (2006).
Powell, W., Morgante, M., McDevitt, R., Vendramin, G. G. & Rafalski, J. A. Polymorphic simple sequence repeat regions in chloroplast genomes: applications to the population genetics of pines. Proc. Natl. Acad. Sci. U.S.A. 92, 7759–7763 (1995).
Kaundun, S. S. & Matsumoto, S. Heterologous nuclear and chloroplast microsatellite amplification and variation in tea, Camellia sinensis. Genome 45, 1041–1048 (2002).
Doorduin, L. et al. The complete chloroplast genome of 17 individuals of pest species Jacobaea vulgaris: SNPs, microsatellites and barcoding markers for population and phylogenetic studies. DNA Res. 18, 93–105 (2011).
Jiao, Y. et al. Development of simple sequence repeat (SSR) markers from a genome survey of Chinese bayberry (Myrica rubra). BMC Genomics 13, 201 (2012).
He, S., Wang, Y., Volis, S., Li, D. & Yi, T. Genetic diversity and population structure: implications for conservation of wild soybean (Glycinesoja Sieb. et Zucc) based on nuclear and chloroplast microsatellite variation. Int. J. Mol. Sci. 13, 12608–12628 (2012).
Zhang, N. et al. An analysis of Echinacea chloroplast genomes: implications for future botanical identification. Sci. Rep. 7, 216 (2017).
Li, H. T. et al. Origin of angiosperms and the puzzle of the Jurassic gap. Nat. Plants 5, 461–470 (2019).
Kearse, M. et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649 (2012).
Qu, X. J., Moore, M. J., Li, D. Z. & Yi, T. S. PGA: a software package for rapid, accurate, and flexible batch annotation of plastomes. Plant Methods 15, 12 (2019).
Tillich, M. et al. GeSeq-versatile and accurate annotation of organelle genomes. Nucleic Acids Res. 45, W6–W11 (2017).
Lowe, T. M. & Chen, P. P. tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 44, W54–W57 (2016).
Lohse, M., Drechsel, O. & Bock, R. OrganellarGenomeDRAW (OGDRAW): a tool for the easy generation of high-quality custom graphical maps of plastid and mitochondrial genomes. Curr. Genet. 52, 267–274 (2007).
Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780 (2013).
Amiryousefi, A., Hyvönen, J. & Poczai, P. IRscope: an online program to visualize the junction sites of chloroplast genomes. Bioinformatics 34, 3030–3031 (2018).
Frazer, K. A., Pachter, L., Poliakov, A., Rubin, E. M. & Dubchak, I. VISTA: computational tools for comparative genomics. Nucleic Acids. Res. 32, W273–W279 (2004).
Darling, A. C., Mau, B., Blattner, F. R. & Perna, N. T. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 14, 1394–1403 (2004).
Librado, P. & Rozas, J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451–1452 (2009).
Kurtz, S. et al. REPuter: the manifold applications of repeat analysis on a genomic scale. Nucleic Acids. Res. 29, 4633–4642 (2001).
Thiel, T., Michalek, W., Varshney, R. & Graner, A. Exploiting EST databases for the development and characterization of gene-derived SSR markers in barley (Hordeumvulgare L.). Theor. Appl. Genet. 106, 411–422 (2003).
Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313 (2014).
Miller, M. A., Pfeiffer, W. & Schwartz, T. Creating the CIPRES science gateway for inference of large phylogenetic trees. New Orleans, LA. Proceedings of the Gateway Computing Environments Workshop (GCE), pp. 1–8 (2010).
Darriba, D., Taboada, G. L., Doallo, R. & Posada, D. jModelTest 2: more models, new heuristics and parallel computing. Nat. Methods 9, 772 (2012).
Swofford, D. L. PAUP*: phylogenetic analysis using Parsimony (*and other methods), version 4.0b10 (Sinauer Associates, Sunderland, 2003).
Acknowledgements
We sincerely thank Prof. Zhong-hu Li, Dr. Jian-jun Jin, De-long Guan and Molecular Biology Experiment Center, Germplasm Bank of Wild Species in Southwest China for help with data analysis. This project was supported by the National Nature Science Foundation of China (No. 31770200, 31300158, 31770203, 31660048) and the Fundamental Research for the Central Universities (No. 2452020179).
Author information
Authors and Affiliations
Contributions
L.Zhao designed the research. Q.J.L., L.Zhang, X.H.Z, Z.Y.C., and L.Zhao. collected the samples, Q.J.L., N.S., L.Zhang, R.C.T., and J.R.W. conceived the experiments, Q.J.L., N.S., L.Z., R.C.T. and J.R.W. did computational analysis and deposited sequences. Q.J.L., N.S., L.Zhang, X.H.Z., Z.Y.C., L.Zhao and D.P. wrote the manuscript. All authors have read and approved the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Li, Qj., Su, N., Zhang, L. et al. Chloroplast genomes elucidate diversity, phylogeny, and taxonomy of Pulsatilla (Ranunculaceae). Sci Rep 10, 19781 (2020). https://doi.org/10.1038/s41598-020-76699-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-76699-7
This article is cited by
-
Twelve newly assembled jasmine chloroplast genomes: unveiling genomic diversity, phylogenetic relationships and evolutionary patterns among Oleaceae and Jasminum species
BMC Plant Biology (2024)
-
Comparative analysis of chloroplast genomes of Pulsatilla species reveals evolutionary and taxonomic status of newly discovered endangered species Pulsatilla saxatilis
BMC Plant Biology (2024)
-
Comparative analysis of complete chloroplast genomes of five Anemone species and phylogenetic analysis within Tribe Anemoneae (Ranunculaceae)
Journal of Plant Biochemistry and Biotechnology (2024)
-
Complete chloroplast genome sequences of the medicinal plant Aconitum transsectum (Ranunculaceae): comparative analysis and phylogenetic relationships
BMC Genomics (2023)
-
Assembly, annotation and analysis of the chloroplast genome of the Algarrobo tree Neltuma pallida (subfamily: Caesalpinioideae)
BMC Plant Biology (2023)