- Research
- Open access
- Published:
Genome skimming as an efficient tool for authenticating commercial products of the pharmaceutically important Paris yunnanensis (Melanthiaceae)
BMC Plant Biology volume 23, Article number: 344 (2023)
Abstract
Background
Paris yunnanensis (Melanthiaceae) is a traditional Chinese medicinal plant of significant pharmaceutical importance. Due to previous taxonomic confusion, a congeneric species, Paris liiana, has been mistaken for P. yunnanensis and cultivated on a large scale, leading to the mixing of commercial products (i.e., seedlings and processed rhizomes) of P. yunnanensis with those of P. liiana. This may have adverse effects on quality control in the standardization of P. yunnanensis productions. As the lack of PCR amplifiable genomic DNA within processed rhizomes is an intractable obstacle to the authentication of P. yunnanensis products using PCR-based diagnostic tools, this study aimed to develop a PCR-free method to authenticate commercial P. yunnanensis products, by applying genome skimming to generate complete plastomes and nrDNA arrays for use as the molecular tags.
Results
Based on a dense intraspecies sampling of P. liiana and P. yunnanensis, the robustness of the proposed authentication systems was evaluated by phylogenetic inferences and experimental authentication of commercial seedling and processed rhizome samples. The results indicate that the genetic criteria of both complete plastomes and nrDNA arrays were consistent with the species boundaries to achieve accurate discrimination of P. yunnanensis and P. liinna. Owing to its desirable accuracy and sensitivity, genome skimming can serve as an effective and sensitive tool for monitoring and controlling the trade of P. yunnanensis products.
Conclusion
This study provides a new way to solve the long-standing problem of the molecular authentication of processed plant products due to the lack of PCR amplifiable genomic DNA. The proposed authentication system will support quality control in the standardization of P. yunnanensis products in cultivation and drug production. This study also provides molecular evidence to clarify the long-standing taxonomic confusion regarding the species delimitation of P. yunnanensis, which will contribute to the rational exploration and conservation of the species.
Background
Paris yunnanensis (Melanthiaceae), a rhizomatous herbaceous perennial plant that is found in southwestern China and northern Myanmar, is a medicinal plant of great economic importance. The well-known traditional Chinese medicine, Rhizoma Paridis, is made from the rhizomes of P. yunnanensis and its congeneric species Paris chinensis [1] (hereafter referred to as ‘medicinal Paris’). For more than 2,000 years, Rhizoma Paridis has been used to treat more than 30 diseases, such as asthma, bleeding, fever, flu, infection, injuries, insect bites, mumps, rheumatic pain syndrome, pruritus, snakebites, swelling, tonsillitis, and ulcers [2, 3]. Recently, modern pharmacological investigations have successively discovered its anti-inflammatory, analgesic, antipyretic, antitumor, immunostimulatory, antimicrobial, and other therapeutic properties [2, 4]. On this basis, more than 90 commercial drugs and health products using Rhizoma Paridis as a raw material have been developed in China [2]. The annual revenue generated by these pharmaceutical products was estimated to be approximately 10 billion CNY (ca. 1.6 billion USD) [5].
In contrast to the substantial demand and consumption of raw materials for drug production (approximately 3 × 106 kg per year), the natural propagation and growth of medicinal Paris and congeneric species are very slow [2, 3, 6, 7]. Species with slow growth rates are generally more prone to overexploitation [8, 9], and the large-scale harvesting of rhizomes from the wild in recent decades has inevitably led to sharp decreases in the natural populations of both P. chinensis and P. yunnanensis [2,3,4,5,6, 10]. With the continuous depletion of natural resources, the pharmaceutical industry that uses Rhizoma Paridis as a raw material has encountered a serious supply crisis [3, 6, 10]. Under this circumstance, the commercial cultivation of medicinal Paris has rapidly expanded in China since the early 2000s, with P. yunnanensis being cultivated on the largest scale to produce raw material for drug production [2, 3].
Notably, P. yunnanensis is highly variable in terms of petal color, folia size and shape, and pedicel length; therefore, the taxonomic delimitation of the species has frequently changed [2, 11,12,13,14,15]. Recently, a total of 4 nominal taxa (i.e., Paris birmanica, Paris daliensis, Paris polyphylla var. emeiensis, and Paris polyphylla var. nana) have been synonymized to P. yunnanensis based on morphological analyses [2, 14]. Inferred from morphological and molecular evidence, a previous study revealed that P. yunnanensis is a collective species consisting of two morphologically distinct and genetically disparate lineages that have little overlap between their distribution ranges; as a result, a new species, Paris liiana, was discovered and described [15]. Due to the previous ambiguity in species delineation, P. liiana has been mistaken for and cultivated as P. yunnanensis on a large scale (Fig. S1). This inevitably leads to the commercial products of P. yunnanensis (i.e., seedlings and processed rhizomes) being mixed with those of P. liiana. In view of the differences in the rhizome chemical compositions and contents between the two species [2, 4, 16], the mixing with P. liiana may have adverse impacts on quality control in either the standardized cultivation of P. yunnanensis or drug production. Given the great importance of P. yunnanensis to the pharmaceutical industry, there is an urgent need for an authentication method that can accurately and sensitively detect contaminating P. liiana among the commercial products of P. yunnanensis.
Although P. liiana and P. yunnanensis are easily distinguished based on the differences in flower and fruit morphologies [2, 15], their seedling and processed rhizomes are morphologically indistinguishable (Fig. 1), making the morphological identification of such commercial products difficult. Since the beginning of this century, interspecific DNA sequence variations have been commonly used for species discrimination [17,18,19,20,21,22,23]. To date, many PCR-based diagnostic tools, such as standard (or lineage-specific) DNA barcodes [10, 24,25,26,27], RAPD fingerprints [28], and sequence-characterized amplified region (SCAR) markers [29,30,31,32,33,34], have been developed for authenticating the commercial products of medicinal plants. Given that genomic DNA of some commercial plant products have been highly degraded during processing, the lack of PCR amplifiable genomic DNA within the processed botanic materials becomes an intractable obstacle to the application of these diagnostic tools that exclusively rely on PCR to recovered target DNA regions for species identification [35,36,37,38,39].
The introduction and development of high-throughput sequencing technologies has made the identification of plant species no longer dependent on PCR [35, 40, 41]. Genome skimming, a PCR-free approach that recovers complete plastid genomes (plastomes) and nuclear tandemly repeated ribosomal DNA (nrDNA) arrays (18 S rDNA–ITS1–5.8 S rDNA–ITS2–26 S rDNA) as the molecular tags from shallow-depth genome sequencing data, has been recommended as an alternative method for plant species identification [21, 35, 42, 43]. Compared with single or multiple sequence regions generated by Sanger sequencing, either complete plastomes or nrDNA arrays possess far more variable loci; therefore, the application of both datasets in plant species discrimination is considered to have great potential for improving resolution [21, 23, 35, 40,41,42,43]. Consistent with this theoretical expectation, empirical studies have shown that although the application of genome skimming is not yet robust enough to distinguish all species, especially in evolutionarily young and complex plant taxa, it can substantially improve the efficacy of species discrimination in most of the studied taxa [38, 41, 44,45,46,47,48,49,50,51,52]. Better yet, previous studies have shown that even using trace and highly degraded genomic DNA extracted from processed plant materials to build shotgun libraries, complete plastomes and nrDNA arrays can be recovered by application of high-throughput sequencing technologies [35, 37, 38], and the use of genome skimming has shown desirable performance in accurate and sensitive identification of species in mixed samples [53, 54].
In view of the favorable advantages of high-throughput sequencing technologies in terms of species identification, this study aimed to develop a PCR-free method for effective authenticating commercial seedlings and processed rhizomes of P. yunnanensis by employing genome skimming to generate complete plastomes and nrDNA arrays for use as the molecular tags. To achieve this goal, individual plants representing the phenotypic diversity of P. liiana and P. yunnanensis were densely sampled and subjected to genome skimming to generate complete plastome and nrDNA array for each sample. On this basis, a reference dataset for the discrimination of P. liiana and P. yunnanensis was built, and the efficacy of genome skimming in identifying the source species of commercial seedling and processed rhizome samples was empirically verified.
Results
Illumina sequencing and recovery of complete plastomes and nrDNA arrays
Illumina sequencing yielded approximately 12.10‒40.61 million paired-end (2 × 150 bp) clean reads per sample, and based on these reads, the reference-guild assembly recovered the complete plastome and entire nrDNA arrays for each sample. These newly generated plastomes (varied from 157,674 to 158,384 bp in length, with average sequencing coverage ranging from 286 to 6,016 times; Table S4) possessed an identical quadripartite structure that consisted of a large single-copy (LSC) and a small single-copy (SSC) separated by two copies of inverted repeat (IR) regions (Table S4 & Table S5). The assembly of the nrDNA arrays completely recovered 18 S rDNA, ITS1, 5.8 S rDNA, ITS2, and 26 rDNA regions for each sample, with the sequence length varying from 5,836 to 5,852 bp and the average sequencing coverage ranging from 375 to 1,487 times (Table S4 & Table S5). The plastomes and nrDNA arrays generated in this study were deposited in the NCBI GenBank database, and their accession numbers are shown in Table S1.
Complete plastomes and nrDNA arrays for discriminating between Paris liiana and Paris yunnanensis
Based on intensive intraspecific sampling, the efficacy of the complete platomes and nrDNA arrays in discriminating between P. liiana and P. yunnanensis was evaluated by examining whether and to what extent each dataset matches the species boundaries. Briefly, ML analysis of complete plastomes (Fig. 2) resolved both species as fully supported monophyletic entities (BS = 100%), and the two species were phylogenetically distinct from each other in the tree topology. Similar to previous studies [49, 55], ML analysis of the nrDNA arrays generated tree topologies that were completely different from the plastome phylogeny (Fig. 3). Even so, both P. liiana and P. yunnanensis were identically resolved as monophyletic entities with full branch support (BS = 100%) by the nrDNA phylogeny. Additionally, neither plastome nor nrDNA phylogeny resolved the four synonymized taxa, i.e., P. birmanica, P. daliensis, P. polyphylla var. emeiensis, and P. polyphylla var. nana, as monophyletic entities but embedded them within accessions of “Typical” P. yunnanensis in the ML trees (Fig. 2 & Fig. 3). The well-supported species-level monophyly of P. liiana and P. yunnanensis showed that the genetic criteria of either the plastome or nrDNA dataset are consistent with the species boundaries of these two species.
Experimental authentication of commercial seedling and processed rhizome samples
Based on the reference plastome and nrDNA datasets, ML phylogeny was determined in order to identify the species to which the commercial seedling and rhizome samples belong. For the commercial seedlings marketed as P. yunnanensis, the ML phylogeny assigned these commercial products to two distinct species. Among them, only three seedlings (Fig. 4) and three process rhizomes (Fig. 5) were clustered with P. yunnanensis (BS = 100%), while seven seedling (Fig. 4) and five rhizome (Fig. 5) samples were clustered with P. liiana (BS = 100%). As the aggregation of query sequences and reference sequences was well supported, the species assignment scheme proposed by phylogenetic analysis can be considered credible [56].
Discussion
Genome skimming for authenticating the commercial products of Paris yunnanensis
The primary goal of the standardized cultivation of medicinal plants is to produce high-quality raw materials that meet therapeutic requirements, which means that quality control should be prioritized in the production of medicinal plants [57, 58]. Due to previous taxonomic ambiguity, commercial seedlings and processed rhizomes of P. yunnanensis are likely mixed with those of P. liiana, which has adverse impacts on quality control in either the standardized cultivation of P. yunnanensis or drug production [58]. Therefore, the detection of P. liiana among the commercial products of P. yunnanensis is crucial for the standardization of P. yunnanensis productions. Given the great importance of P. yunnanensis to the pharmaceutical industry, an accurate and sensitive authentication method is urgently needed to distinguish the commercial seedlings and processed rhizomes of P. yunnanensis from those of P. liiana. In view of the intractable obstacle of PCR-based diagnostic tools in authenticating such commercial plant products (i.e., the processed rhizomes of Paris species) whose genomic DNA have been highly degraded during processing [35,36,37,38,39], this study aimed to develop a PCR-free method for accurate and effective authentication of the commercial seedlings and processed rhizomes of P. yunnanensis, using genome skimming to generate complete plastomes and nrDNA arrays for use as the molecular tags.
The availability of a high-quality reference dataset is an indispensable prerequisite for the application of genome skimming for species discrimination [35, 59,60,61,62]. Accordingly, the creation of robust reference complete plastome and nrDNA datasets is the essential first step in the practical application of genome skimming for the authentication of P. yunnanenis and its commercial products. In this study, the complete plastomes and nrDNA arrays of 25 P. yunnanensis and 16 P. liinna individuals, representing the intraspecific phenotypic diversity of both species, were sampled to construct the reference dataset. These sequences were obtained from correctly identified specimens, and the sequences constituted comprehensive and high-quality complete reference plastome and nrDNA datasets. Additionally, phylogenetic analysis of both datasets resolved P. yunnanensis and P. liinna as monophyletic entities with high statistical supports, indicating that the genetic criteria of complete plastomes and nrDNA arrays were consistent with the species boundaries of P. yunnanensis and P. liinna. This suggests that accurate discrimination of P. yunnanensis and P. liinna can be achieved using either complete plastomes or nrDNA arrays as the molecular tags.
To validate the efficacy of the proposed authentication system, this study experimentally applied genome skimming to identify the source species of commercial seedling and processed rhizome samples. Although the processed rhizomes of Paris species represent botanic materials with highly degraded genomic DNA, the complete plastomes and nrDNA arrays of all samples were recovered from Illumina reads. The desirable success rate of recovering complete plastomes and nrDNA arrays suggests that they are ideal molecular tags for authentication the processed rhizomes of Paris yunnanensis. Based on the high-quality reference plastome and nrDNA datasets, ML phylogeny unambiguously detected P. liiana among commercial seedling and processed rhizomes marketed as P. yunnanensis, validating the desired accuracy and practicality of the proposed authentication system.
To date, numerous PCR-based diagnostic tools have been developed for authenticating the commercial products of medicinal plants [10, 24,25,26,27,28,29,30,31,32,33,34]. Although these methods exhibit the merits of time and cost effectiveness, the success of their authentication depends largely on the DNA quality and integrity of examined samples. This should limit their practical application in the authentication of some processed plant products, where genomic DNA have been highly degraded during processing [35,36,37,38,39].
The current study represents the first case to develop a PCR-free method to authenticate commercial products of P. yunnanensis, a pharmaceutically important plant species. Although considerable bioinformatic skills are needed to recover complete plastomes and nrDNA arrays from low coverage genome sequencing data and to perform phylogenetic analyses, the method developed in this study provides a good solution for authentication of processed plant products that lack PCR amplifiable genomic DNA. Recently, the advances in high-throughput sequencing technology have largely decreased the cost of low coverage genome sequencing, and the development of new approaches have made it easier to recover complete plastomes and nrDNA arrays from shotgun sequencing data [63,64,65,66]. Accordingly, the authentication system developed in this study can serve as an effective and sensitive tool for monitoring and controlling the trade of seedlings and processed rhizomes, which will contribute to quality control in the standardized cultivation of P. yunnanensis and drug production.
Implications for conservation and exploitation of Paris yunnanensis
This study provides further evidence supporting the taxonomic revision that considered P. birmanica, P. daliensis, P. polyphylla var. emeiensis, and P. polyphylla var. nana to be synonyms for P. yunnanensis [2, 14]. Specifically, neither plastome nor nrDNA phylogeny resolved the four synonymized taxa, i.e., P. birmanica, P. daliensis, P. polyphylla var. emeiensis, and P. polyphylla var. nana, as monophyletic entities but embedded them within accessions of “Typical” P. yunnanensis in the phylogenetic trees, suggesting that these taxa have not genetically diverged from “Typical” P. yunnanensis. Accordingly, either recognizing P. birmanica and P. daliensis as segregated species [12, 67] or treating P. polyphylla var. emeiensis and P. polyphylla var. nana as conspecific varieties of P. polyphylla [68, 69] is likely a taxonomic error. This phylogenetic inference is in good agreement with morphological evidence: except for the purple petals of P. birmanica [12, 68, 70], the thickened free portion of the connectivum in P. daliensis [67], a nearly absent pedicel in P. polyphylla var. emeiensis [69], and an extremely shortened (less than 5 cm) pedicel in P. polyphylla var. nana [68], the four taxa are highly similar to “typical” P. yunnanensis in terms of their leaf and flower morphologies (Fig. 6). This suggests that the diagnostic characteristics that were used to define the four taxa may represent intraspecific morphological variations or phenotypic plasticity [2, 14], justifying to synonymize P. birmanica, P. daliensis, P. polyphylla var. emeiensis, and P. polyphylla var. nana with P. yunnanensis [2, 14].
On the other hand, P. yunnanensis is a medicinal plant whose conservation has attracted substantial concerns, as overexploitation of the species for medicinal purposes has resulted in the rapid depletion of its natural populations [2, 3]. Given that accurate delineation of species boundaries is the essential first step to properly address issues regarding species conservation and exploitation [71,72,73], this study clarifies the long-standing taxonomic confusion and provides robust phylogenetic evidence to delineate a clear-cut species boundary of P. yunnanensis, thus has great potential to improve strategies for the conservation and exploitation of the economically important species. Specifically, this study determines that previous taxonomic confusion has resulted in the proliferation of as many as four synonyms in P. yunnanensis. Among them, P. daliensis is categorized as critically endangered [74], P. polyphylla var. emeiensis and P. polyphylla var. nana are classified as endangered [69, 74], and P. birmanica is classified as vulnerable [75]. This suggests that previous taxonomic errors may have resulted in the inappropriate allocation of limited conservation resources to these taxa that are not at high risk of extinction. Additionally, previous taxonomic errors in species delimitation must have led to the omission of these synonymized taxa from the introduction and cultivation of P. yunnanensis and its germplasm conservation. Therefore, individual plants and populations of both “typical” P. yunnanensis and these synonymized taxa should be simultaneously collected to preserve the germplasm resources, from which elite germplasm can be excavated for breeding needs.
The collection and preservation of propagules is the preferred strategy for the conservation of plant germplasm resources [76, 77]. For the propagules (i.e., rhizomes, seeds, and seedlings) of P. yunnanensis, traditional morphology-based species identification is a fairly difficult task due to the absence of diagnostic morphologies in these organs. Nevertheless, based on the reference datasets of complete plastome and nrDNA arrays, determining whether propagules that will be collected and preserved belong to P. yunnanesis can be easy: genome skimming can be applied to recover their complete plastomes or nrDNA arrays, and then, these sequences can be queried against the corresponding reference dataset.
Conclusions
For the first time, this study developed a PCR-free authentication system for detecting whether the commercial seedlings and processed rhizomes of Paris yunnanensisis are adulterated with the congeneric species, Paris liiana, by applying genome skimming to generate complete plastomes and nrDNA arrays for use as the molecular tags. Although the authentication system relies heavily on bioinformatic skills and is not time- and cost-efficient, this study provides a new way to solve the long-standing problem of the molecular authentication of processed plant products due to their lack of PCR amplifiable genomic DNA. Owing to its desirable accuracy and sensitivity, the authentication system will be conductive to the standardization of P. yunnanensisis products to guarantee their quality and effectiveness.
Methods
Plant sampling, Illumina sequencing, and recovery of complete plastomes and nrDNA arrays
In this study, the complete plastomes and nrDNA arrays of 25 P. yunnanensis and 16 P. liinna individuals were sampled to establish the reference dataset for the use of genome skimming to distinguish between P. liinna and P. yunnanensis. Among them, 12 and 10 individual plants, whose morphological characteristics represent the phenotypic diversity of P. yunnanensis and P. liinna, respectively, were sequenced for the first time in this study. The plant materials sampled in this study was identified by Dr. Yunheng Ji. The voucher specimens were deposited at herbarium of Kunming Institute of Botany (Chinese Academy of Sciences), and the original sources of the plant samples and voucher information are presented in Table S1. The remaining complete plastomes and nrDNA arrays were obtained from the NCBI GenBank database (Table S2). Notably, the intraspecific sampling of P. yunnanensis included not only individual plants with the “typical” morphological characteristics of the species but also individuals that morphologically represented the four synonymized taxa (Table S1), i.e., Paris birmanica (3 individuals), Paris daliensis (3 individuals), Paris polyphylla var. emeiensis (2 individuals), and Paris polyphylla var. nana (2 individuals), which were reduced to P. yunnanensis by Liang and Soukup [14] and Ji [2]. In addition, to experimentally validate the efficacy of genome skimming for the detection of P. liiana among the commercial seedlings and rhizomes of P. yunnanensis, a total of 10 commercial seedlings and eight processed rhizome samples marketed as P. yunnanensis (voucher information is presented in Table S3) were sampled for shotgun sequencing.
The genomic DNA of these samples was extracted with the CTAB method [78]. Illumina sequencing libraries with an average insert size of approximately 400 bp were built using a TruSeq DNA PCR-free Prep Kit (Illumina Inc., San Diego, CA, USA) according to the manufacturer’s protocol and then subjected to paired-end sequencing on an Illumina Novaseq 6000 platform to generate approximately two Gbp of raw reads for each sample. The software Trimmomatic v0.40 [79] was used to remove adaptors and to filter low-quality reads with default parameters. A reference-guided method was used for the assembly of complete plastomes and nrDNA arrays for each sample; the published complete plastome (GenBank accession: MN175247) and nrDNA sequence (GenBank accession: MN647572) of P. liiana were used as the reference, and the GetOrganelle v1.7.5.0 pipeline [66] was use to recover complete plastomes and nrDNA arrays from the trimmed Illumina reads, with default parameters and preset options. The obtained plastomes were annotated with the online program Geseq v2.03 [80] and further validated by performing a BLAST search against the NCBI protein database with Geneious v10.2.3 [81]. For the nrDNA arrays, the 28 S, 18 S, and 5.8 S ribosomal RNA genes and their boundaries with the intergenic transcribed spacer (ITS1 and ITS2) regions were annotated and defined by comparison with the reference sequence in Geneious v10.2.3 [81].
Data analysis
According to phylogenetic tree topologies, the robustness of the complete plastomes and nrDNA arrays for the discrimination of P. liiana and P. yunnanensis was assessed by examining the monophyly of each species. In addition to the reference sequences of P. liiana and P. yunnanensis, 24 complete plastomes and nrDNA arrays that represent the remaining congeneric species were included in phylogenetic analysis (Table S2), with Trillium tschonoskii as the outgroup, according to the findings of previous studies [55, 82,83,84]. The complete plastomes and nrDNA arrays were independently aligned with MAFFT v7.402 [85] using the default parameters. For each dataset, phylogenetic relationships were reconstructed using the maximum likelihood (ML) method. The ML analyses were conducted with the online software IQ-TREE v2.2.0 [63, 64], and 1,000 ultrafast bootstrap (BS) replicates were performed to generate the percentage of BS support for each node.
To investigate whether genome skimming enables the accurate and sensitive detection of P. liiana among the commercial seedlings and rhizomes of P. yunnanensis, complete plastomes and nrDNA arrays from 10 commercial seedling and eight processed rhizome samples marketed as P. yunnanensis were incorporated into the plastome and nr DNA datasets to reconstruct the ML phylogeny. For each dataset, ML analyses were conducted using the previously described protocol. According to tree topology, the phylogenetic placement of each sample was investigated to determine whether the ML phylogeny accurately assigned these samples to their source species. In general, if the query sequence was clustered with the reference sequences from a certain taxon with high branch support, credible species assignment was considered to have been achieved [56].
Data Availability
The complete plastomes and nrDNA arrays generated in this study are available at the NCBI GenBank database, with the following accession numbers: OP021412, OP021413, OP021414, OP021415, OP021474, OP021475, OP021476, OP021477, OP021478, OP021493, OP021494, OP021495, OP021496, OP021793, OP021797, OP021802, OP021804, OP021806, OP021817, OP021820, OP021823, OP021830, OP021831, OP021835, OP021843, OP021845, OP806312, OP806313, OP806314, OP806315, OP806316, OP806317, OP837990, OP837991, OP837992, OP837993, OP837994, and OP837995.
Abbreviations
- bp:
-
Base pair
- BP:
-
Bootstrap percentage
- BS:
-
Bootstrap
- CTAB:
-
Cetyl trimethylammonium bromide
- DNA:
-
Deoxyribonucleic acid
- Gb:
-
Giga base pairs
- IR:
-
Inverted repeat
- ITS:
-
internal transcribed spacer of nuclear ribosomal DNA
- LSC:
-
Large single-copy
- MCMC:
-
Markov Chain Monte Carlo
- ML:
-
Maximum Likelihood
- PP:
-
Posterior probability
- rRNA:
-
Ribosomal RNA
- SSC:
-
Small single copy
- tRNA:
-
Transfer RNA
References
China Pharmacopoeia Committee. Pharmacopoeia of the people’s Republic of China. Sec I, the 2020 edition. Beijing:China Medical Science Press; 2020.
Ji Y. A monograph of Paris. Beijing and Singapore: Science Press and Springer; 2021.
Li H, Su B, Zhang ZY, et al. An assessment on the rarely medical Paris plants in China with exploring the future development of its plantation. J West China For Sci. 2015a;44:1–6. https://doi.org/10.16473/j.cnki.xblykx1972.2015.03.001.
Ding YG, Zhao YL, Zhang J, et al. The traditional uses, phytochemistry, and pharmacological properties of Paris L. (Liliaceae): a review. J Ethnopharmacol. 2021;278:114293. https://doi.org/10.1016/j.jep.2021.114293.
Huang LQ, Xiao PG, Wnag YY. Investigation on resources of rare and endangered medicinal plants in China. Shanghai: Shanghai Science and Technology Press; 2012.
He J, Zhang S, Wang H, et al. Advances in studies on and uses of Paris polyphylla var. Yunnanensis (Trilliaceae). Acta Bot Yunnan. 2006;28:271–6.
Cunningham AB, Brinckmann JA, Bi YF, et al. Paris in the spring: a review of the trade, conservation and opportunities in the shift from wild harvest to cultivation of Paris polyphylla (Trilliaceae). J Ethnopharmacol. 2018;222:208–16. https://doi.org/10.1016/j.jep.2018.04.048.
Turner IM. Species loss in fragments of tropical rain forest: a review of the evidence. J Appl Ecol. 1996;33:200–9. https://doi.org/10.2307/2404743.
Stork NE, Coddington JA, Colwell RK, et al. Vulnerability and resilience of tropical forest species to land–use change. Conserv Biol. 2009;23:1438–47. https://doi.org/10.1111/j.1523-1739.2009.01335.x.
Yang Y, Zhai YH, Liu T, et al. Detection of Valeriana jatamansi as an adulterant of medicinal Paris by length variation of chloroplast psb A -trnH region. Planta Med. 2011;77:87–91. https://doi.org/10.1055/s-0030-1250072.
Hara H. Variations in Paris polyphylla Smith with reference to other Asiatic species. J Fac Sci Univ Tokyo Sect III. 1969;10:141–80.
Takhtajan A. A revision of Daiswa (Trilliaceae). Britton. 1983;35:255–70. https://doi.org/10.2307/2806025.
Li HThe. Genus Paris. Beijing: Science Press; 1998.
Liang SY, Soukup VG. Paris L. In: Wu ZY, Raven PH,editor. Flora of China. Beijing: Science Press and St Louis, MO: Botanical Garden Press; 2000. p. 88–95.
Ji Y, Liu C, Yang J, et al. Ultra-barcoding discovers a cryptic species in Paris yunnanensis (Melanthiaceae), a medicinally important plant. Front Plant Sci. 2020;11:411. https://doi.org/10.3389/fpls.2020.00411.
Wang JL, Yang QS, Jiang Y, et al. Geo-herbalism study of Paridis Rhizoma. Chin Trad Herb Drugs. 2022;53:2572–81.
Brower AVZ, DeSalle R, Vogler A. Gene trees, species trees and systematics: a cladistic perspective. Annu Rev Ecol Syst. 1996;27:423–50. https://doi.org/10.1146/annurev.ecolsys.27.1.423.
Hebert PDN, Cywinska A, Ball SL, et al. Biological identifications through DNA barcodes. Proc R Soc Lond B Biol Sci. 2003;270:313–21. https://doi.org/10.1098/rspb.2002.2218.
Hollingsworth PM, Forrest LL, Spouge JL, et al. A DNA barcode for land plants. Proc Natl Acad Sci. 2009;106:12794–7. https://doi.org/10.1073/pnas.0905845106.
Hollingsworth PM, Li DZ, van der Bank M, et al. Telling plant species apart with DNA: from barcodes to genomes. Philos Trans R Soc B Biol Sci. 2016;371:20150338. https://doi.org/10.1098/rstb.2015.0338.
Hollingsworth PM, Graham SW, Little DP. Choosing and using a plant DNA barcode. PLoS ONE. 2011;6:e19254. https://doi.org/10.1371/journal.pone.0019254.
Hollingsworth PM. Refining the DNA barcode for land plants. Proc Natl Acad Sci U S A. 2011;108:19451–2. https://doi.org/10.1073/pnas.1116812108.
Gostel MR, Kress WJ. The expanding role of DNA barcodes: indispensable tools for ecology, evolution, and conservation. Diversity. 2022;14:213. https://doi.org/10.3390/d14030213.
Chen S, Pang X, Song J, et al. A renaissance in herbal medicine identification: from morphology to DNA. Biotechnol Adv. 2014;32:1237–44. https://doi.org/10.1016/j.biotechadv.2014.07.004.
Techen N, Parveen I, Pan Z, et al. DNA barcoding of medicinal plant material for identification. Curr Opin Biotechnol. 2014;25:103–10. https://doi.org/10.1016/j.copbio.2013.09.010.
Jiao J, Huang WL, Bai ZQ, et al. DNA barcoding for the efficient and accurate identification of medicinal polygonate rhizoma in China. PLoS ONE. 2018;13:e0201015. https://doi.org/10.1371/journal.pone.0201015.
Hua Z, Jiang C, Song S, et al. Accurate identification of taxon-specific molecular markers in plants based on DNA signature sequence. Mol Ecol Resour. 2023;23:106–17. https://doi.org/10.1111/1755-0998.13697.
Zhang KY, Leung HW, Yeung HW, et al. Differentiation of Lycium barbarum from its related Lycium species using random amplified polymorphic DNA. Planta Med. 2001;67:379–81. https://doi.org/10.1055/s-2001-14310.
Kiran U, Khan S, Mirza KJ, et al. SCAR markers: a potential tool for authentication of herbal drugs. Fitoterapia. 2010;81:969–76. https://doi.org/10.1016/j.fitote.2010.08.002.
Ganie SH, Upadhyay P, Das S, et al. Authentication of medicinal plants by DNA markers. Plant Gene. 2015;4:83–99. https://doi.org/10.1016/j.plgene.2015.10.002.
Kim WJ, Moon BC, Yang S, et al. Rapid authentication of the herbal medicine plant species Aralia continentalis Kitag. And Angelica biserrata CQ Yuan and RH Shan using ITS2 sequences and multiplex-SCAR markers. Molecules. 2016;21:270. https://doi.org/10.3390/molecules21030270.
Wetters S, Horn T, Nick P. Goji who? Morphological and DNA based authentication of a “Superfood. Front Plant Sci. 2018;9:1859. https://doi.org/10.3389/fpls.2018.01859.
Shen Z, Lu T, Zhang Z, et al. Authentication of traditional chinese medicinal herb “Gusuibu” by DNA-based molecular methods. Ind Crops Prod. 2019;141:111756. https://doi.org/10.1016/j.indcrop.2019.111756.
Nob P, Kim WJ, Yang S, et al. PCR-based rapid diagnostic tools for the authentication of medicinal. Phytomed. 2021;91:153667. https://doi.org/10.1016/j.phymed.2021.153667.
Li XW, Yang Y, Henry RJ, et al. Plant DNA barcoding: from gene to genome: plant identification using DNA barcodes. Biol Rev. 2015b;90:157–66. https://doi.org/10.1111/brv.12104.
Raclariu AC, Heinrich M, Ichim MC, et al. Benefits and limitations of DNA barcoding and metabarcoding in herbal product authentication. Phytochem Anal. 2018;29:123–8. https://doi.org/10.1002/pca.273.
Zeng C, Hollingsworth PM, Yang J, et al. Genome skimming herbarium specimens for DNA barcoding and phylogenomics. Plant Methods. 2018;14:43. https://doi.org/10.1186/s13007-018-0300-0.
Guo XR, Shi NX, Xie PX, et al. Plastome sequencing for accurate and effective authentication of Polygonatum kingianum (Asparagaceae). Ind Crops Prod. 2022;184:115056. https://doi.org/10.1016/j.indcrop.2022.115056.
Mahima K, Sunil Kumar KN, Rakhesh KV, et al. Advancements and future prospective of DNA barcodes in the herbal drug industry. Front Pharmacol. 2022;13:947512. https://doi.org/10.3389/fphar.2022.947512.
Nock CJ, Water DLE, Edwards MA, et al. Chloroplast genome sequences from total DNA for plant identification: chloroplast genome sequences for plant identification. Plant Biotechnol J. 2011;9:328–33. https://doi.org/10.1111/j.1467-7652.2010.00558.x.
Kane N, Sveinsson S, Dempewolf H, et al. Ultra-barcoding in cacao (Theobroma spp.; Malvaceae) using whole chloroplast genomes and nuclear ribosomal DNA. Am J Bot. 2012;99:320–9. https://doi.org/10.3732/ajb.1100570.
Dodsworth S. Genome skimming for next-generation biodiversity analysis. Trends Plant Sci. 2015;20:525–7. https://doi.org/10.1016/j.tplants.2015.06.012.
Coissac E, Hollingsworth PM, Lavergne S, et al. From barcodes to genomes: extending the concept of DNA barcoding. Mol Ecol. 2016;25:1423–8. https://doi.org/10.1111/mec.13549.
Ruhsam M, Rai HS, Mathews S, et al. Does complete plastid genome sequencing improve species discrimination and phylogenetic resolution in Araucaria ? Mol Ecol Resour. 2015;15:1067–78. https://doi.org/10.1111/1755-0998.12375.
Turner B, Paun O, Munzinger J, et al. Sequencing of whole plastid genomes and nuclear ribosomal DNA of Diospyros species (Ebenaceae) endemic to New Caledonia: many species, little divergence. Ann Bot. 2016;117:1175–85. https://doi.org/10.1093/aob/mcw060.
Zhu SY, Niu ZT, Xue QY, et al. Accurate authentication of Dendrobium officinale and its closely related species by comparative analysis of complete plastomes. Acta Pharm Sin B. 2018;8:969–80. https://doi.org/10.1016/j.apsb.2018.05.009.
Ji Y, Liu CK, Yang ZY, et al. Testing and using complete plastomes and ribosomal DNA sequences as the next generation DNA barcodes in Panax (Araliaceae). Mol Ecol Resour. 2019a;19:1333–45. https://doi.org/10.1111/1755-0998.13050.
Ji Y, Yang J, Landis J, et al. Deciphering the taxonomic delimitation of Ottelia acuminata (Hydrocharitaceae) using complete plastomes as super-barcodes. Front Plant Sci. 2021;12:681270. https://doi.org/10.3389/fpls.2021.681270.
Ji Y, Yang J, Landis JB, et al. Genome skimming contributes to clarifying species limits in Paris section Axiparis (Melanthiaceae). Front Plant Sci. 2022;13:832034. https://doi.org/10.3389/fpls.2022.832034.
Fu C, Mo Z, Yang J, et al. Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendron. Mol Ecol Resour. 2022;22:404–14. https://doi.org/10.1111/1755-0998.13479.
Yu XQ, Jiang YZ, Folk RA, et al. Species discrimination in Schima (Theaceae): next-generation super‐barcodes meet evolutionary complexity. Mol Ecol Resour. 2022;22:3161–75. https://doi.org/10.1111/1755-0998.13683.
Wang SY, Zhou N, Shi NX, et al. Testing and using complete plastomes for authentication of medicinal Polygonatum species (Asparagaceae). Ind Crops Prod. 2023;197:116557. https://doi.org/10.1016/j.indcrop.2023.116557.
Lang D, Tang M, Hu J, et al. Genome-skimming provides accurate quantification for pollen mixtures. Mol Ecol Resour. 2019;19:1433–46. https://doi.org/10.1111/1755-0998.1306.
Bohmann K, Mirarab S, Bafna V, et al. Beyond DNA barcoding: the unrealized potential of genome skim data in sample identification. Mol Ecol. 2020;29:2521–34. https://doi.org/10.1111/mec.15507.
Ji Y, Yang L, Chase M, et al. Plastome phylogenomics, biogeography, and clade diversification of Paris (Melanthiaceae). BMC Plant Biol. 2019b;19:543. https://doi.org/10.1186/s12870-019-2147-6.
Elias M, Hill RI, Willmott KR, et al. Limited performance of DNA barcoding in a diverse community of tropical butterflies. Proc R Soc B Biol Sci. 2007;274:2881–9. https://doi.org/10.1098/rspb.2007.1035.
Canter PH, Thomas H, Ernst E. Bringing medicinal plants into cultivation: opportunities and challenges for biotechnology. Trends Biotechnol. 2005;23:180–5. https://doi.org/10.1016/j.tibtech.2005.02.002.
Zhang BG, Peng Y, Zhang Z, et al. GAP production of TCM herbs in China. Planta Med. 2010;76:1948–55. https://doi.org/10.1055/s-0030-1250527.
Ekrem T, Willassen E, Stur E. A comprehensive DNA sequence library is essential for identification with DNA barcodes. Mol Phylogenet Evol. 2007;43:530–42. https://doi.org/10.1016/j.ympev.2006.11.021.
Miller SE, Hausmann A, Hallwachs W, et al. Advancing taxonomy and bioinventories with DNA barcodes. Philos Trans R Soc B Biol Sci. 2016;371:20150339. https://doi.org/10.1098/rstb.2015.0339.
Morinière J, Hendrich L, Balke M, et al. A DNA barcode library for Germany′s mayflies, stoneflies and caddisflies (Ephemeroptera, Plecoptera and Trichoptera). Mol Ecol Resour. 2017;17:1293–307. https://doi.org/10.1111/1755-0998.12683.
Zhang N, Ramachandran P, Wen J, et al. Development of a reference standard library of chloroplast genome sequences, GenomeTrakrCP. Planta Med. 2017;83:1420–30. https://doi.org/10.1055/s-0043-113449.
Nguyen LT, Schmidt HA, von Haeseler A, et al. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32:268–74. https://doi.org/10.1093/molbev/msu300.
Trifinopoulos J, Nguyen L-T, von Haeseler A, et al. W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res. 2016;44:W232–5. https://doi.org/10.1093/nar/gkw256.
Kumar S, Stecher G, Li M, et al. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9. https://doi.org/10.1093/molbev/msy096.
Jin JJ, Yu WB, Yang JB, et al. GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol. 2020;21:241. https://doi.org/10.1186/s13059-020-02154-5.
Li H. New taxa of the genus Paris L. Acta Bot 1992;Suppl 5:13–8.
Li H. A study on the taxonomy of the genus Paris L. Bull Bot Res. 1986;6:111–46.
Yin HX, Zhang H, Xue D. Paris polyphylla var. Emeiensis H. X. Yin, H. Zhang and D. Xue, a new variety of Trilliaceae from Sichuan, China. Acta Phytotaxon Sin. 2007;45:823–7. https://doi.org/10.1360/aps06160.
Li H, Noltie H. Miscellaneous notes on the genus Paris. Edinb J Bot. 1997;54:351–3.
Hopkins GW, Freckleton RP. Declines in the numbers of amateur and professional taxonomists: implications for conservation. Anim Conserv. 2002;5:245–9. https://doi.org/10.1017/S1367943002002299.
Sites JW, Marshall JC. Delimiting species: a Renaissance issue in systematic biology. Trends Ecol Evol. 2003;18:462–70. https://doi.org/10.1016/S0169-5347(03)00184-8.
Mace GM. The role of taxonomy in species conservation. Philos Trans R Soc Lond B Biol Sci. 2004;359:711–9. https://doi.org/10.1098/rstb.2003.1454.
Qin HN, Yang Y, Dong SY, et al. Threatened species list of China’s higher plants. Biodivers Sci. 2017;25:696–744.
Yang J, Yang X, Oo TN, et al. Rediscovery and supplemental description of Paris birmanica (Melanthiaceae), a species endemic to Myanmar. Phytotaxa. 2019;391:240–2. https://doi.org/10.11646/phytotaxa.391.3.8.
Margules CR, Pressey RL. Systematic conservation planning. Nature. 2000;405:243–53. https://doi.org/10.1038/35012251.
Milligan AL, Putwain PD, Cox ES, et al. Developing an integrated land management strategy for the restoration of moorland vegetation on Molinia caerulea-dominated vegetation for conservation purposes in upland Britain. Biol Conserv. 2004;119:371–85. https://doi.org/10.1016/j.biocon.2003.12.002.
Doyle JJ, Doyle JL. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 1987;19:11–5.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20. https://doi.org/10.1093/bioinformatics/btu170.
Tillich M, Lehwark P, Pellizzer T, et al. GeSeq–versatile and accurate annotation of organelle genomes. Nucleic Acids Res. 2017;45:W6–W11. https://doi.org/10.1093/nar/gkx391.
Kearse M, Moir R, Wilson A, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9. https://doi.org/10.1093/bioinformatics/bts199.
Ji Y, Fritsch PW, Li H, et al. Phylogeny and classification of Paris (Melanthiaceae) inferred from DNA sequence data. Ann Bot. 2006;98:245–56. https://doi.org/10.1093/aob/mcl095.
Huang YL, Li XJ. Analysis of complete chloroplast genome sequences improves phylogenetic resolution in Paris (Melanthiaceae). Front Plant Sci. 2016;7. https://doi.org/10.3389/fpls.2016.01797.
Yang L, Yang Z, Liu C, et al. Chloroplast phylogenomic analysis provides insights into the evolution of the largest eukaryotic genome holder, Paris japonica (Melanthiaceae). BMC Plant Biol. 2019;19:293. https://doi.org/10.1186/s12870-019-1879-7.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80. https://doi.org/10.1093/molbev/mst010.
Acknowledgements
The authors are grateful to Shuying Wang, Zhangming Wang, Lifang Yang, Tingzhou Zhao, Guohua Zhou for their helps in sample collection.
Funding
The present study is supported by the National Natural Science Foundation of China (31872673), the NSFC-Joint Foundation of Yunnan Province (U1802287), the Key R & D program of Yunnan Province (202103AC100003), and the Academicians and Experts Workstation of Yunnan Province (202105AF150074).
Author information
Authors and Affiliations
Contributions
Yunheng Ji conceived and designed the research framework. Nian Zhou, Lilei Tang, Pingxuan Xie, and Ke Miao collected and analyzed the data. Nian Zhou, Lilei Tang, Pingxuan Xie, and Ke Miao wrote the original draft manuscript. Haiyang Liu, Chengjin Yang, and Yunheng Ji revised and edited the final manuscript. All authors have read and agreed to the published version of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Ethics approval and consent to participate
Collection of all plant samples in this study completely followed the Regulations on the Protection of Wild Plants of the People’s Republic of China, the IUCN Policy Statement on Research Involving Species at Risk of Extinction and the Convention on the Trade in Endangered Species of Wild Fauna and Flora. All the methods included in this study are in accordance with the relevant guidelines.
Consent for publication
Not applicable.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
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/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Zhou, N., Tang, L., Xie, P. et al. Genome skimming as an efficient tool for authenticating commercial products of the pharmaceutically important Paris yunnanensis (Melanthiaceae). BMC Plant Biol 23, 344 (2023). https://doi.org/10.1186/s12870-023-04365-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12870-023-04365-x