Skip to main content

Genome skimming as an efficient tool for authenticating commercial products of the pharmaceutically important Paris yunnanensis (Melanthiaceae)

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.

Peer Review reports

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].

Fig. 1
figure 1

Morphologies of arial shoot (a), flower (b), fruit (c), seedling (d), rhizome (e), and processed rhizomes (f) of Paris yunnanensis and Paris liiana

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.

Fig. 2
figure 2

Cladogram of the maximum-likelihood (ML) phylogeny resulted from analyzing the complete plastomes of Paris species. The numbers at nodes represents ML bootstrap (BS) percentages

Fig. 3
figure 3

Cladogram of the maximum-likelihood (ML) phylogeny resulted from analyzing the nuclear ribosome DNA (nrDNA) arrays of Paris species. The numbers at nodes represents ML bootstrap (BS) percentages

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].

Fig. 4
figure 4

Species assignment of the commercial seedlings marketed as Paris yunnanensis (highlighted in red) via maximum-likelihood (ML) analyses of complete plastomes (a) and nuclear ribosome DNA (nrDNA) arrays (b). The numbers at nodes represents ML bootstrap (BS) percentages

Fig. 5
figure 5

Species assignment of the processed rhizomes marketed as Paris yunnanensis (highlighted in red) via maximum-likelihood (ML) analyses of complete plastomes (a) and nuclear ribosome DNA (nrDNA) arrays (b). The numbers at nodes represents ML bootstrap (BS) percentages

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].

Fig. 6
figure 6

Morphology of “typical” Paris yunnanensis (a) and synonymized taxa: Paris birmanica (b), Paris daliensis (c), Paris polyphylla var. nana (d), and Paris polyphylla var. emeiensis (e)

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

  1. China Pharmacopoeia Committee. Pharmacopoeia of the people’s Republic of China. Sec I, the 2020 edition. Beijing:China Medical Science Press; 2020.

  2. Ji Y. A monograph of Paris. Beijing and Singapore: Science Press and Springer; 2021.

    Google Scholar 

  3. 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.

    Article  Google Scholar 

  4. 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.

    Article  CAS  PubMed  Google Scholar 

  5. Huang LQ, Xiao PG, Wnag YY. Investigation on resources of rare and endangered medicinal plants in China. Shanghai: Shanghai Science and Technology Press; 2012.

    Google Scholar 

  6. 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.

    CAS  Google Scholar 

  7. 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.

    Article  CAS  PubMed  Google Scholar 

  8. 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.

    Article  Google Scholar 

  9. 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.

    Article  PubMed  Google Scholar 

  10. 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.

    Article  CAS  PubMed  Google Scholar 

  11. Hara H. Variations in Paris polyphylla Smith with reference to other Asiatic species. J Fac Sci Univ Tokyo Sect III. 1969;10:141–80.

    Google Scholar 

  12. Takhtajan A. A revision of Daiswa (Trilliaceae). Britton. 1983;35:255–70. https://doi.org/10.2307/2806025.

    Article  Google Scholar 

  13. Li HThe. Genus Paris. Beijing: Science Press; 1998.

    Google Scholar 

  14. 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.

  15. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Wang JL, Yang QS, Jiang Y, et al. Geo-herbalism study of Paridis Rhizoma. Chin Trad Herb Drugs. 2022;53:2572–81.

    Google Scholar 

  17. 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.

    Article  Google Scholar 

  18. 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.

    Article  CAS  Google Scholar 

  19. 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.

    Article  PubMed Central  Google Scholar 

  20. 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.

    Article  CAS  Google Scholar 

  21. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 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.

    Article  CAS  Google Scholar 

  24. 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.

    Article  CAS  PubMed  Google Scholar 

  25. 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.

    Article  CAS  PubMed  Google Scholar 

  26. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 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.

    Article  CAS  PubMed  Google Scholar 

  28. 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.

    Article  CAS  PubMed  Google Scholar 

  29. 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.

    Article  CAS  PubMed  Google Scholar 

  30. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 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.

    Article  CAS  Google Scholar 

  34. 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.

    Article  CAS  Google Scholar 

  35. 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.

    Article  PubMed  Google Scholar 

  36. 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.

    Article  CAS  PubMed  Google Scholar 

  37. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. 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.

    Article  CAS  Google Scholar 

  39. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. 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.

    Article  CAS  PubMed  Google Scholar 

  41. 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.

    Article  CAS  PubMed  Google Scholar 

  42. 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.

    Article  CAS  PubMed  Google Scholar 

  43. 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.

    Article  CAS  PubMed  Google Scholar 

  44. 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.

    Article  CAS  PubMed  Google Scholar 

  45. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  47. 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.

    Article  CAS  PubMed  Google Scholar 

  48. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  49. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  50. 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.

    Article  CAS  PubMed  Google Scholar 

  51. 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.

    Article  CAS  PubMed  Google Scholar 

  52. 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.

    Article  CAS  Google Scholar 

  53. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  56. 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.

    Article  CAS  Google Scholar 

  57. 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.

    Article  CAS  PubMed  Google Scholar 

  58. 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.

    Article  CAS  PubMed  Google Scholar 

  59. 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.

    Article  CAS  PubMed  Google Scholar 

  60. 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.

    Article  Google Scholar 

  61. 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.

    Article  CAS  PubMed  Google Scholar 

  62. 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.

    Article  CAS  PubMed  Google Scholar 

  63. 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.

    Article  CAS  PubMed  Google Scholar 

  64. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Li H. New taxa of the genus Paris L. Acta Bot 1992;Suppl 5:13–8.

  68. Li H. A study on the taxonomy of the genus Paris L. Bull Bot Res. 1986;6:111–46.

    Google Scholar 

  69. 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.

    Article  Google Scholar 

  70. Li H, Noltie H. Miscellaneous notes on the genus Paris. Edinb J Bot. 1997;54:351–3.

    Google Scholar 

  71. 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.

    Article  Google Scholar 

  72. 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.

    Article  Google Scholar 

  73. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Qin HN, Yang Y, Dong SY, et al. Threatened species list of China’s higher plants. Biodivers Sci. 2017;25:696–744.

    Article  Google Scholar 

  75. 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.

    Article  Google Scholar 

  76. Margules CR, Pressey RL. Systematic conservation planning. Nature. 2000;405:243–53. https://doi.org/10.1038/35012251.

    Article  CAS  PubMed  Google Scholar 

  77. 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.

    Article  Google Scholar 

  78. Doyle JJ, Doyle JL. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 1987;19:11–5.

    Google Scholar 

  79. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  82. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. 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.

  84. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

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

Authors

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

Correspondence to Yunheng Ji.

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

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-023-04365-x

Keywords