Genome Skimming Contributes to Clarifying Species Limits in Paris Section Axiparis (Melanthiaceae)

Abstract

Paris L. section Axiparis H. Li (Melanthiaceae) is a taxonomically perplexing taxon with considerable confusion regarding species delimitation. Based on the analyses of morphology and geographic distribution of each species currently recognized in the taxon, we propose a revision scheme that reduces the number of species in P. sect. Axiparis from nine to two. To verify this taxonomic proposal, we employed a genome skimming approach to recover the plastid genomes (plastomes) and nuclear ribosomal DNA (nrDNA) regions of 51 individual plants across the nine described species of P. sect. Axiparis by sampling multiple accessions per species. The species boundaries within P. sect. Axiparis were explored using phylogenetic inference and three different sequence-based species delimitation methods (ABGD, mPTP, and SDP). The mutually reinforcing results indicate that there are two species-level taxonomic units in P. sect. Axiparis (Paris forrestii s.l. and P. vaniotii s.l.) that exhibit morphological uniqueness, non-overlapping distribution, genetic distinctiveness, and potential reproductive isolation, providing strong support to the proposed species delimitation scheme. This study confirms that previous morphology-based taxonomy overemphasized intraspecific and minor morphological differences to delineate species boundaries, therefore resulting in an overestimation of the true species diversity of P. sect. Axiparis. The findings clarify species limits and will facilitate robust taxonomic revision in P. sect. Axiparis.

Introduction

Species delimitation is the crucially important first step for designing research in many fields of biology (Mace, 2004). Traditionally, species boundaries are defined by taxonomists based on the analysis of morphological variation (Lewin, 1981; Henderson, 2005), which has resulted in massive disagreements over species identification and delimitation due to either phenotypic plasticity or lack of taxonomically robust morphological characters at the species level (Hebert et al., 2003; De Queiroz, 2007). To compensate, additional data types, such as molecular and ecological profiles, are needed to explicitly decipher species boundaries (Sites and Crandall, 1997; Sites and Marshall, 2003; Ellis et al., 2006; Bickford et al., 2007; Duminil and Michele, 2009; Su et al., 2015; Eisenring et al., 2016; Lambert et al., 2017; Sukumaran and Knowles, 2017; Cheng et al., 2020; Ma et al., 2020).

Analysis of DNA sequence variation can provide useful information for identifying and delineating species (Brower et al., 1996; Hebert et al., 2003; Kress et al., 2005; Pons et al., 2006; Hollingsworth et al., 2009, 2011, 2016; Hollingsworth, 2011; Duminil et al., 2012; Puillandre et al., 2012). With next-generation DNA sequencing (NGS) technologies, genome-wide sequence variation has begun to replace one or a few sequence loci for the identification and delimitation of plant species (Li et al., 2015; Coissac et al., 2016; Hollingsworth et al., 2016). The genome skimming approach, which uses NGS technologies to generate multi-copy and highly repetitive genome components, such as whole plastid genomes (plastomes) and nuclear ribosomal DNA (nrDNA) clusters via relatively low coverage genome sequencing (Straub et al., 2012), has been increasingly used for species identification and delimitation in recent years (Nock et al., 2011; Kane et al., 2012; Dodsworth, 2015; Li et al., 2015; Ruhsam et al., 2015; Firetti et al., 2017; Fu et al., 2019, 2021; Ji et al., 2019a,2020, 2021; Knope et al., 2020; Ślipiko et al., 2020; Su et al., 2021). Compared with restriction-site associated DNA sequencing (RAD-seq; Miller et al., 2007; Baird et al., 2008), another NSG-based technique that is extensively used to generate genomic data for plant species identification and delimitation (e.g., Wu et al., 2018; Donkpegan et al., 2020; Ma et al., 2020; Zhou et al., 2020; Li et al., 2021), the promising advantage of using genome skimming for species identification and delimitation is the avoidance of problems encountered with RAD-seq (Kane et al., 2012; Hollingsworth et al., 2016), such as only applying to diploids and generating asymmetric data between distinctly related taxa.

This study focuses on Paris Linn. section Axiparis H. Li (Melanthiaceae), a taxonomically perplexing plant group that includes nine described species distributed from Central China to the Himalayas (Li, 1998; Ji et al., 2006, 2019b; Huang et al., 2016). Since the description of the first two species (Paris vaniotii and P. forrestii) of the section (Léveillé, 1906; Takhtajan, 1983), a total of four (P. axialis, P. guizhouensis, P. lihengiana, and P. variabilis) and three (P. dulongensis, P. rugosa, and P. tengchongensis) species whose morphologies are similar to P. vaniotii and P. forrestii respectively have been described (Li, 1984, 1992; He, 1990; Ji et al., 2017; Xu et al., 2019; Yang et al., 2019). The rapid accumulation in the number of species over the last 40 years has led to considerable taxonomic confusion in P. sect. Axiparis. After critical examination of types and specimens assigned to these species, it was found that, except for slight differences in leaf shape and size (Figure 1), these recently described species exhibit high levels of similarity in flower, fruit, and seed morphology with P. vaniotii and P. forrestii (Ji, 2021). Remarkably, leaf shape and size of foliar displays have high levels of intraspecific variation that sometimes far exceeds the divergence used for species diagnosis (Ji, 2021). Also noteworthy, we found that P. dulongensis, P. rugosa, and P. tengchongensis share the morphological similarity of stamens numbering twice as many as the sepal number (two-whorled stamens), and being geographically distributed from the Hengduan Mountains (Southwestern China) to the Himalayas (Figure 2). Comparatively, P. axialis, P. guizhouensis, P. lihengiana, P. vaniotii, and P. variabilis, have stamens numbering three times that of the sepal number (three-whorled stamen), and these species are distributed from Central China to the Wumeng Mountains in Southwestern China (Figure 2). The overlap of morphological features and species ranges implies that (1) P. dulongensis, P. rugose, and P. tengchongensis are likely conspecific with P. forrestii, and (2) P. axialis, P. guizhouensis, P. lihengiana, and P. variabilis may belong to P. vaniotii. Accordingly, previous morphological-based taxonomic studies (Li, 1984; He, 1990; Li, 1992; Ji et al., 2017; Xu et al., 2019; Yang et al., 2019) may have overemphasized intraspecific and minor morphological variations to establish species, therefore led to the proliferation of synonyms in P. sect. Axiparis.

FIGURE 1

FIGURE 2

Due to the above-mentioned aspects, P. sect. Axiparis is in great need of taxonomic revision. In this study, we attempt to use genome skimming to establish a basic understanding of species boundaries in this taxonomically perplexing taxon. Specifically, we generated complete plastomes and nrDNA sequences by sampling multiple accessions per species within P. sect. Axiparis via low coverage genome sequencing. Based on phylogenetic inference and multiple sequence-based species delimitation methods, we aim to test our taxonomic proposal that reduces the number of species in P. sect. Axiparis from nine to two (Paris forrestii s.l. and P. vaniotii s.l.).

Materials and Methods

Plant Samples, Illumina Sequencing, Assembly, and Annotation

In total, we sampled 51 individual plants including representatives of each of the nine described species of P. sect. Axiparis, in which 46 accessions were newly sequenced in the present study (Table 1). Voucher specimens and leaf tissue were collected from the field. For each taxon, multiple individuals within a species representing different localities were included, and at least one accession was harvested from the type locality. The sampling strategy entirely covers the nine species’ distribution ranges and allows robust exploring of species-level monophyly and species boundaries within P. sect. Axiparis to critically test for our taxonomic proposal.

Total genomic DNA for each accession was isolated from ∼10 mg silica gel dried leaf tissues using cetyltrimethylammonium bromide (CTAB; Doyle and Doyle, 1987). Approximately 5 μg of purified genomic DNA was used to construct PCR-free shotgun libraries with a TruSeq DNA Sample Prep Kit (Illumina, Inc., San Diego, CA, United States) according to the manufacturer’s instructions. Paired-end (150 bp) sequencing was performed on the Illumina HiSeq 2500 platform (Illumina, Inc., San Diego, CA, United States) to generate approximately two giga base pairs (Gbp) of raw data for each sample. Remarkably, Paris is a fairly distinctive angiosperm genus in possessing giant genomes (Pellicer et al., 2014; Ji, 2021), and the minimum documented genome size in the genus (P. bashanensis, 1C = 28.73 Gbp; Ji, 2021) is much larger than the mean genome size of angiosperms (1C = 5.7 Gbp; Pellicer et al., 2014). Due to the relatively low sequencing coverage and lack of nuclear genome reference, it is difficult to recover a sufficient number of unlinked and single-copy nuclear loci from the genome skimming data. Therefore, only complete plastomes and nrDNA sequences were assembled in this study for phylogenetic reconstruction and species delimitation analyses.

Shotgun reads of each sample were deposited in NCBI short-read archive (SRA), with accession numbers being shown in Supplementary Table 1. Trimmomatic v0.40 (Bolger et al., 2014) was used to remove low-quality reads and adaptors from the Illumina raw reads with default parameters. The filtered reads were assembled into complete plastomes and nrDNA clusters with the GetOrganelle pipeline v1.7.5.0 (Jin et al., 2020) on a Linux system, using the previously published complete plastome (MN125565) and nrDNA sequence (MN174877) of P. forrestii (Ji et al., 2019b) as the reference. The assembled plastomes were annotated using the online software Geseq v2.03 (Tillich et al., 2017), with default parameters. The transfer RNA (tRNA) genes were further verified using tRNAscan-SE v2.0 (Lowe and Chan, 2016) with default parameters. For nrDNA annotation, the ribosomal RNA genes (26S, 18S, and 5.8S ribosomal RNA genes) and boundaries with the intergenic transcribed spacer (ITS) regions were annotated and defined by comparison with the reference in Geneious v10.2.3 (Kearse et al., 2012).

Data Analyses

The complete plastome and entire nrDNA clusters (including 18S rRNA, ITS1, 5.8S rRNA, ITS2, and 26S rRNA) were aligned using MAFFT v7.450 (Katoh and Standley, 2013). The alignments of sequences were deposited the online database Treebase,¹ and only single nucleotide polymorphisms (SNPs) were included in phylogenetic and species delimitation analyses. Partition homogeneity tests (Farris et al., 1994) were performed with PAUP* 4.0b10 (Swofford, 2002) to determine the degree of congruence between the concatenated plastome and nrDNA matrices, as well as between coding (18S rRNA, 5.8S rRNA, and 26S rRNA) and non-coding regions (ITS1 and ITS2) of nrDNA sequences, using a heuristic tree search algorithm for 500 replicates.

Phylogenetic relationships were inferred using both maximum likelihood (ML) and Bayesian inference (BI) methods, based on which we tested whether the nine described species within P. sect. Axiparis are monophyletic units. Previous studies revealed that ancient intergeneric hybridization may have occurred between Trillium and P. sect. Paris (resulting in the speciation of P. japonica: P. sect. Kinugasa), and past intersectional hybridization likely took place between P. sect. Axiparis and P. sect. Euthyra (Ji et al., 2019b; Ji, 2021). Given that such reticulate relationships may result in phylogenetic errors (Philippe et al., 2011), we selected P. thibetica (complete plastome: MN125569; nrDNA: MN174890; Ji et al., 2019b) as the outgroup, which represents P. sect. Thibeticae, the closest relative of P. sect. Axiparis (Ji et al., 2006, 2019b; Ji, 2021). The best-fit sequence substitution model (GTR + G for concatenated plastome matrix, and GTR + I + G for concatenated nrDNA matrix) were determined by Modeltest v3.7 (Posada and Crandall, 1998) using the Akaike information criterion (Posada and Buckley, 2004). The ML analysis was performed with RAxML-HPC BlackBox v8.1.24 (Stamatakis, 2006). The best-scoring ML tree for each dataset was produced with 1,000 bootstrap (BS) replicates to provide support values for each node. The BI tree was inferred using MrBayes v3.2 (Ronquist et al., 2012). Two independent Markov chain Monte Carlo (MCMC) runs were performed with one million generations, sampling every 100 generations, with the initial 25% of the sampled trees as burn-in. Posterior probability (PP) values were computed based on the remaining trees.

Relying on reciprocal monophyly alone to delineate species boundaries likely under-represents true species-level diversity (Ence and Carstens, 2011; Rannala and Yang, 2013). The development of some coalescence-based species delimitation methods, such as Bayesian Phylogenetics and Phylogeography (BPP; Rannala and Yang, 2013), and spedeSTEM (Ence and Carstens, 2011), provide the solution to this issue since these methods do not require reciprocal monophyly of any species in a given gene tree. Both BPP and spedeSTEM are only suitable for processing genomic segments no longer than 500 or 1,000 bp (Ence and Carstens, 2011; Rannala and Yang, 2013), and thus cannot be used to analyze either complete plastomes or entire nrDNA sequences. Accordingly, we used the following approaches to explore species boundaries within P. sect. Axiparis based on each dataset (concatenated plastome and nrDNA matrices): (1) the distance-based method automatic barcode gap discovery (ABGD; Puillandre et al., 2011), (2) the tree-based method multi-rate Poisson tree processes model (mPTP; Kapli et al., 2017), and (3) the coalescence-based method species delimitation plugin (SDP; Masters et al., 2011) in Geneious v.10.2.3 (Kearse et al., 2012). The outgroup (P. thibetica) sequences were removed from the species delimitation analyses.

The ABGD analyses were conducted with the online server² with default settings (Pmin = 0.001, Pmax = 0.1, Steps = 10, X = 1.5, Nb = 20). All three genetic distance models (JC69, K2P, and uncorrected P-distances) specified by the program were used. Next, we used the mPTP v0.2.3 algorithm (Kapli et al., 2017), an improvement to PTP (Zhang et al., 2013), to delineate species boundaries and to estimate the posterior probability (PP) values for the putative species. The mPTP analyses were performed on the web server³ with standard default settings (The MCMC algorithm was conducted for 0.1 million generations, sampling every 100 generations and a 10% burn-in), using the inferred ML trees of plastome and nrDNA as inputs since the branch lengths of ML tree represent expected numbers of substitutions per site. Based on ML and BI trees of the concatenated plastome and nrDNA matrices, we performed SDP analyses to generate species delimitation scheme in P. sect. Axiparis. The distinctiveness of these candidate species proposed by SDP analyses was estimated using Rosenberg’s P_(AB): the probability of reciprocal monophyly under a random coalescent model (Rosenberg, 2007), and Rodrigo’s P_(RD): the probability that a clade has the observed degree of distinctiveness due to a random coalescent process (Rodrigo et al., 2008).

Results

Genome Skimming

The summary of low coverage genome sequencing and assembly of plastome and nrDNA sequences is presented in Supplementary Table 1. Plastome assembly generated the complete plastome of all samples, which possess a typical quadripartite structure, with the sequence length varying from 156,061 to 157,653 bp. The plastomes identically contain 114 genes, including 80 protein-coding genes, 30 tRNA genes, and four plastid rRNA genes (Supplementary Table 2). In addition, the assembly of nrDNA entirely covered 18S, ITS1, 5.8S, ITS2, and 26S regions in all accessions, with the sequence length being either 5,851 or 5,852 bp. DNA sequences of the newly generated plastomes and nrDNA sequences in this study were deposited in GenBank, with accession numbers being shown in Table 1.

Phylogenetic Inferences

Alignment of the plastome sequences yielded a matrix of 160,681 positions, in which we identified 1,724 variable sites (1.07%) with 1,192 (0.74%) being parsimoniously informative (Supplementary Table 3). The ML (Figure 3A) and BI (Figure 3B) analyses of the concatenated plastome matrix produced similar tree topologies except that several shallow nodes with low BS support collapse in BI tree. Two diverging clades were recovered in P. sect. Axiparis: Clade I (BS = 100%, PP = 1.00) comprises accessions of P. dulongensis, P. forrestii, P. rugosa, and P. tengchongensis, and Clade II (BS = 100%, PP = 1.00) includes accessions of P. axialis, P. guizhouensis, P. lihengiana, P. vaniotii, and P. variabilis. Strikingly, none of the nine nominal species within P. sect. Axiparis were recovered as a monophyletic unit in either ML or BI phylogeny.

FIGURE 3

Although partition homogeneity tests detected significant conflicts (p < 0.01) between the nrDNA and plastome phylogenies, the coding (18S rRNA, 5.8S rRNA, and 26S rRNA) and non-coding regions (ITS1 and ITS2) of nrDNA sequences are congruent with each other in the tree topology (p > 0.05). Therefore, we concatenated both coding and non-coding regions of nrDNA sequences for phylogenetic and species delimitation analyses. The concatenated nrDNA matrix possessed 193 variable sites (3.30%) with 112 (1.91%) being parsimoniously informative (Supplementary Table 3). ML (Figure 4A) and BI (Figure 4B) analyses generated congruent tree topologies despite that some weakly supported shallow nodes in ML tree collapse in BI phylogeny. The nrDNA trees recovered two major clades within P. sect. Axiparis, corresponding to those inferred from the plastomes but with slightly lower support values (clade I: BS = 98%, PP = 1.00; clade II: BS = 97%, PP = 1.00). Similar to the plastome phylogeny, analyses of nrDNA sequences failed to resolve any of the nine nominal species within P. sect. Axiparis as a monophyletic unit.

FIGURE 4

Species Delimitation Scenarios

Overall, ABGD, mPTP, and SDP analyses of the concatenated plastome and nrDNA matrices produced highly congruent results that are reflected in the phylogenetic trees (Figures 3, 4). The ABGD analyses of the concatenated plastome (Table 2) and nrDNA (Table 3) matrices resulted in a stable count (n = 2) of species division with a range of prior intraspecific values (plastomes: P = 0.001000–0.004642; nrDNA: P = 0.002000–0.007573) with JC69, K2P, and P-distances initial and recursive partitions. One proposed species includes individuals of P. dulongensis, P. forrestii, P. rugosa, and P. tengchongensis, while the other comprises individuals of P. axialis, P. guizhouensis, P. lihengiana, P. vaniotii, and P. variabilis (Figures 3, 4). The mPTP analyses yielded the same delimitation proposal as ABGD: all individuals were grouped into two species-like entities that coincide with the two clades recovered by our phylogenetic analyses (Figures 3, 4), and both received strong posterior support (Table 4). Similarly, the SDP analyses showed that either plastome or nrDNA matrix, excluding the outgroup, comprised two putative species (corresponding to Clade I and Clade II in the inferred phylogenetic tree; Figures 3, 4), which were distinctive with P_(RD) < 0.05, and significant with P_(AB) < 10^–5 (Table 5).

Discussion

The advancement of molecular-based approaches has brought about great progress in species delimitation (Wiley and Mayden, 2000; Sites and Marshall, 2003; De Queiroz, 2007). Concomitantly, a wide variety of species delimitation methods have been developed, especially over the past 20 years (Yang and Rannala, 2010; Ence and Carstens, 2011; Masters et al., 2011; Puillandre et al., 2011; Rannala and Yang, 2013; Zhang et al., 2013). However, no method is perfect and each has specific weaknesses. For instance, ABGD is problematic when species are represented with only a few specimens (Puillandre et al., 2011, 2012). In some cases, mPTP may lead to an overestimation of species numbers (Song et al., 2018; Huang et al., 2020). Given that factors such as recent diversification, radiative speciation, and restricted intraspecific gene flow may result in the absence of reciprocal monophyly among closely related species (Ruhsam et al., 2015; Hollingsworth et al., 2016); SDP may yield biased delimitation schemes since it primarily relies on reciprocal monophyly to explore species boundaries (Masters et al., 2011). Accordingly, multiple methods should be simultaneously employed to develop a robust species delimitation framework, given that using different delimitation approaches allows accommodations for the weaknesses of each approach (Hebert et al., 2004; Fujita et al., 2012; Aguilar et al., 2013; Kekkonen and Hebert, 2014; Mutanen et al., 2015). However, previous studies using a single or few sequence regions have indicated that employing multiple methods on the same dataset always produced incongruent delimitation proposals (Camargo et al., 2012; Mutanen et al., 2015; Jacobs et al., 2018; Jirapatrasilp et al., 2019) likely due to the inadequacies of available genetic information to properly delineate species boundaries (Carstens and Satler, 2013; Jacobs et al., 2018).

High variability of key diagnostic morphological characters makes species delimitation a particularly difficult task, and thus has led to considerable taxonomic confusion in P. sect. Axiparis. This study aims to critically test our taxonomical proposal that only recognizes two broadly defined species (vs. nine narrowly defined species) in P. sect. Axiparis. To achieve this goal, we employed genome skimming to recover complete plastomes and nrDNA sequences that contain more evolutionarily informative variation than a single or few sequence regions for phylogenetic and species delimitation analyses. Based on the concatenated plastome and nrDNA matrices, we not only investigated the species-level monophyly of the nine narrowly defined species by sampling multiple individuals per species, but also explored species boundaries using multiple sequence-based species delimitation methods (ABGD, mPTP, and SDP). To determine the robustness of the delimitation proposals generated by ABGD, mPTP, and SDP analyses, we also examined to what extent these methods generate congruent results.

Phylogenetic analyses recovered none of the nine narrowly defined species within P. sect. Axiparis as a monophyletic unit, suggesting that the genetic differentiation among them is low. By contrast, both plastome and nrDNA phylogenies recovered two well-supported clades that possess fairly distinct morphological traits (Figures 3, 4) and distribution ranges (Figure 2). Although members of P. sect. Axiparis represent recently diverged entities with their origins estimated no earlier than the late Miocene (Ji et al., 2019b), ABGD, mPTP, and SDP analyses based on different models generated consistent proposals in delineating species boundaries. Briefly, the ABGD analyses (Tables 2, 3) partitioned all the samples into two clusters comprised of individuals having two-whorled and (Clade I) three-whorled stamens (Clade II). This implies that there are likely two distinct species within P. sect. Axiparis with significant genetic gaps between them (Puillandre et al., 2012). Moreover, the mPTP analyses grouped all accessions into two putative species with high posterior delimitation probability, coinciding with the results found in the ABGD analyses (Figures 3, 4, and Table 4). The putative species boundaries proposed by ABGD and mPTP analyses are further validated by the SDP analyses: with P_(RD) value < 0.05 and the P_(AB) value < 10^–5 (Table 5), both of the two candidate species can be recognized as evolutionarily distinctive entities (Rosenberg, 2007; Rodrigo et al., 2008). These reciprocally reinforcing results suggest that only two species-level taxonomic units can be recognized in P. sect. Axiparis, providing robust support to our taxonomic proposal to reduce P. axialis, P. guizhouensis, P. lihengiana, and P. variabilis as the synonyms of P. vaniotii, and to expand the species boundary of P. forrestii to accommodate P. dulongensis, P. rugose, and P. tengchongensis.

Importantly, sequence-based delimitation results can only be considered a hypothesis that needs to be further validated with multiple data types, such as phenotypic and ecological information (Sukumaran and Knowles, 2017). In addition to morphological distinctiveness (two-whorled stamens vs. three-whorled stamens), there is a range of evidence robustly supporting the species delimitation scheme proposed by our data. Within Paris, significant conflict between the plastome and nrDNA trees was also detected in a previous study: although the plastome phylogeny recovered P. sect. Axiparis as a well-supported clade, it was resolved as non-monophyletic by nrDNA phylogeny (Ji et al., 2019b). By comparing differences in plastome and nrDNA tree topologies, ancient hybridization was proposed to have occurred between the common ancestor of P. luquanensis and P. mairei (P. sect. Euthyra) and that of P. dulongensis, P. forrestii, P. rugose, and P. tengchongensis (Ji et al., 2019b). This implies that the two broadly defined species, namely P. forrestii s.l. (syn. P. dulongensis, P. rugose, and P. tengchongensis) and P. vaniotii s.l. (syn. P. axialis, P. guizhouensis, P. lihengiana, and P. variabilis), originated separately. Also noteworthy, P. forrestii s.l. and P. vaniotii s.l. possess allopatric geographic distributions, between which there are great climatic differences. From the late Miocene onward, driven by the intensification of the East Asian monsoon in the summer (An et al., 2001), the species range of P. vaniotii s.l. has been primarily governed by Pacific monsoon, whereas that of P. forrestii s.l. has been mainly affected by Indian monsoon (Li et al., 2008; Yao et al., 2011; Figure 2). The long-term environmental and ecological heterogeneity resulted from climate differences between the two distinct monsoon regimes would block the regional gene flow and boosted vicariance in Paris (Ji et al., 2019b). Based on phenology recorded in herbarium specimens and observation of plants cultivated in common gardens, it was found that the flowering of P. vaniotii s.l. is approximately 30–40 days earlier than that of P. forrestii s.l., and there is little overlap between their flowering periods (Ji, 2021). This implies that outcrossing between P. forrestii s.l. and P. vaniotii s.l. is nearly impossible under natural conditions, suggesting that reproductive isolation between them may have been formed under the combined effect of geographic isolation and habitat heterogeneity. Therefore, recognizing P. forrestii s.l. and P. vaniotii s.l. as distinct species under the unified species concept (De Queiroz, 2007) is reasonable, given that the species boundary proposed by our data reflects the unity of morphological uniqueness, non-overlapping species range, genetic distinctiveness, and reproductive isolation.

Conclusion

This study represents a guiding practical application of genome skimming for exploring species boundaries in taxonomically perplexing plant taxa, using not only phylogenetic inferences but also multiple sequence-based species delimitation methods. The analyses of concatenated plastome and nrDNA matrices yielded identical schemes in delineating species boundaries in Paris sect. Axiparis, which are highly congruent with morphological characteristics and geographic distributions. The results robustly support our revision proposal that reduces the number of species in P. sect. Axiparis from nine to two, and confirm that previous morphology-based taxonomy (Li, 1984, 1992; He, 1990; Ji et al., 2017; Xu et al., 2019; Yang et al., 2019) overemphasized intraspecific and minor morphological differences to establish species, therefore resulted in over-splitting of species in P. sect. Axiparis.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1

  AguilarC.WoodP. L.CusiJ. C.GuzmánA.HuariF.LundbergM.et al (2013). Integrative taxonomy and preliminary assessment of species limits in the Liolaemus walkeri complex (Squamata, Liolaemidae) with descriptions of three new species from Peru.Zookeys36447–91. 10.3897/zookeys.364.6109

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AnZ. S.KutzbachJ. E.PrellW. L.PorterS. C. (2001). Evolution of Asian monsoons and phased uplift of the Himalaya–Tibetan plateau since Late Miocene times.Nature41162–66. 10.1038/35075035

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BairdN. A.EtterP. D.AtwoodT. S.CurreyM. C.ShiverA. L.LewisZ. A.et al (2008). Rapid SNP discovery and genetic mapping using sequenced RAD markers.PLoS One3:e3376. 10.1371/journal.pone.0003376

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  BickfordD.LohmanD. J.SodhiN. S.NgP. K.MeierR.WinkerK.et al (2007). Cryptic species as a window on diversity and conservation.Trends Ecol. Evol.22148–155. 10.1016/j.tree.2006.11.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data.Bioinformatics302114–2120. 10.1093/bioinformatics/btu170

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  BrowerA. V. Z.DeSalleR.VoglerA. (1996). Gene trees, species trees and systematics: a cladistic perspective.Annu. Rev. Ecol. Syst.27423–450. 10.1146/annurev.ecolsys.27.1.423

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  CamargoA.MorandoM.AvilaL. J.SitesJ. W. (2012). Species delimitation with ABC and other coalescent-based methods: a test of accuracy with simulations and an empirical example with lizards of the Liolaemus darwinii complex (Squamata: Liolaemidae).Evolution662834–2849. 10.1111/j.1558-5646.2012.01640.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  CarstensB. C.SatlerJ. D. (2013). The carnivorous plant described as Sarracenia alata contains two cryptic species.Biol. J. Linn. Soc.109737–746. 10.1111/bij.12093

  • CrossRef
  • Google Scholar

• 9

  ChengZ.ShuH.ZhangS.LuoB.GuR.ZhangR.et al (2020). From folk taxonomy to species confirmation of Acorus (Acoraceae): evidences based on phylogenetic and metabolomic analyses.Front. Plant. Sci.11:965. 10.3389/fpls.2020.00965

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  CoissacE.HollingsworthP. M.LavergneS.TaberletP. (2016). From barcodes to genomes: extending the concept of DNA barcoding.Mol. Ecol.251423–1428. 10.1111/mec.13549

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  De QueirozK. (2007). Species concepts and species delimitation.Syst. Bot.56879–886. 10.1080/10635150701701083

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  DodsworthS. (2015). Genome skimming for next-generation biodiversity analysis.Trends Plant Sci.20525–527. 10.1016/j.tplants.2015.06.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  DonkpeganA. S. L.DoucetJ. L.HardyO. J.HeuertzM.PiñeiroR. (2020). Miocene diversification in the Savannahs precedes tetraploid rainforest radiation in the African tree genus Afzelia (Detarioideae. Fabaceae).Front. Plant Sci.11:798. 10.3389/fpls.2020.00798

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  DoyleJ. J.DoyleJ. (1987). A rapid DNA isolation procedure for small quantities of fresh leaf tissue.Phytochem. Bull.1911–15.

  • Google Scholar
  • View reference in article

• 15

  DuminilJ.KenfackD.ViscosiV.GrumiauL.HardyO. J. (2012). Testing species delimitation in sympatric species complexes: the case of an African tropical tree. Carapa spp. (Meliaceae).Mol. Phylogenet. Evol.62275–285. 10.1016/j.ympev.2011.09.020

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  DuminilJ.MicheleM. D. (2009). Plant species delimitation: a comparison of morphological and molecular markers.Plant Biosyst.143528–542. 10.1080/11263500902722964

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  EisenringM.AltermattF.WestramA. M.JokelaJ. (2016). Habitat requirements and ecological niche of two cryptic amphipod species at landscape and local scales.Ecosphere7:e01319. 10.1002/ecs2.1319

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  EllisJ. S.KnightM. E.CarvellC.GoulsonD. (2006). Cryptic species identification: a simple diagnostic tool for discriminating between two problematic bumblebee species.Mol. Ecol. Notes6540–542. 10.1111/j.1471-8286.2006.01231.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  EnceD. D.CarstensB. C. (2011). SpedeSTEM: a rapid and accurate method for species delimitation.Mol. Ecol. Resour.11473–480. 10.1111/j.1755-0998.2010.02947.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  FarrisJ.S.KällersjöM.KlugeA.C.BultC. (1994). Testing significance of incongruence.Cladistics10315–319. 10.1111/j.1096-0031.1994.tb00181.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  FirettiF.ZuntiniA. R.GaiarsaJ. W.OliveiraR. S.LohmannL. G.Van SluysM. A. (2017). Complete chloroplast genome sequences contribute to plant species delimitation: a case study of the Anemopaegma species complex.Am. J. Bot.1041493–1509. 10.3732/ajb.1700302

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  FuC. N.MoZ. Q.YangJ. B.CaiJ.YeL. J.ZouJ. Y.et al (2021). Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendron.Mol. Ecol. Resour.22404–414. 10.1111/1755-0998.13479

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  FuC. N.WuC. S.YeL. J.MoZ. Q.LiuJ.ChangY. W.et al (2019). Prevalence of isomeric plastomes and effectiveness of plastome super-barcodes in yews (Taxus) worldwide.Sci. Rep.9:2773. 10.1038/s41598-019-39161-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  FujitaM. K.LeachéA. D.BurbrinkF. T.McGuireJ. A.MoritzC. (2012). Coalescent-based species delimitation in an integrative taxonomy.Trends Ecol. Evol.27480–488. 10.1016/j.tree.2012.04.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  HeS. (1990). A new medicinal species of Paris from Guizhou province, China.Guizhou Sci.3:16.

  • Google Scholar
  • View reference in article

• 26

  HebertP. D.CywinskaA.BallS. L.De-WaardJ. R. (2003). Biological identifications through DNA barcodes.Proc. Biol. Sci.270313–321. 10.1098/rspb.2002.2218

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  HebertP. D.PentonE. H.BurnsJ. M.JanzenD. H.HallwachsW. (2004). Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator.Proc. Natl. Acad. Sci. U.S.A.10114812–14817. 10.1073/pnas.0406166101

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  HendersonA. (2005). The methods of herbarium taxonomy.Syst. Bot.30456–459.

  • Google Scholar
  • View reference in article

• 29

  HollingsworthP. M. (2011). Refining the DNA barcode for land plants.Proc. Natl. Acad. Sci. U.S.A.10819451–19452. 10.1073/pnas.1116812108

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  HollingsworthP. M.ForrestL. L.SpougeJ. L.HajibabaeiM.RatnasinghamS.Van der BankM.et al (2009). A DNA barcode for land plants.Proc. Natl. Acad. Sci. U.S.A.10612794–12797. 10.1073/pnas.0905845106

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  HollingsworthP. M.GrahamS. W.LittleD. P. (2011). Choosing and using a plant DNA barcode.PLoS One6:e19254. 10.1371/journal.pone.0019254

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  HollingsworthP. M.LiD. Z.Van der BankM.TwyfordA. D. (2016). Telling plant species apart with DNA: from barcodes to genomes.Philos. Trans. R. Soc. B Biol. Sci.371:20150338. 10.1098/rstb.2015.0338

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  HuangW.XieX.HuoL.LiangX.WangX.ChenX. (2020). An integrative DNA barcoding framework of ladybird beetles (Coleoptera: Coccinellidae).Sci. Rep.10:10063. 10.1038/s41598-020-66874-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  HuangY.LiX.YangZ.YangC.YangJ.JiY. (2016). Analysis of complete chloroplast genome sequences improves phylogenetic resolution in Paris (Melanthiaceae).Front. Plant Sci.7:1797. 10.3389/fpls.2016.01797

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  JacobsS. J.KristoffersonC.Uribe-ConversS.LatvisM.TankD. C. (2018). Incongruence in molecular species delimitation schemes: What to do when adding more data is difficult.Mol. Ecol.272397–2413. 10.1111/mec.14590

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  JiY. (2021). A Monograph of Paris (Melanthaceae).Beijing: Science Press.

  • Google Scholar
  • View reference in article

• 37

  JiY.FritschP. W.LiH.XiaoT.ZhouZ. (2006). Phylogeny and classification of Paris (Melanthiaceae) inferred from DNA sequence data.Ann. Bot.98245–256. 10.1093/aob/mcl095

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  JiY.LiuC.YangJ.JinL.YangZ.YangJ. B. (2020). Ultra-barcoding discovers a cryptic species in Paris yunnanensis (Melanthiaceae), a medicinally important plant.Front. Plant Sci.11:411. 10.3389/fpls.2020.00411

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  JiY.LiuC.YangZ.YangL.HeZ.WangH.et al (2019a). Testing and using complete plastomes and ribosomal DNA sequences as the next generation DNA barcodes in Panax (Araliaceae).Mol. Ecol. Resour.191333–1345. 10.1111/1755-0998.13050

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  JiY.YangL.ChaseM. W.LiuC.YangZ.YangJ.et al (2019b). Plastome phylogenomics, biogeography, and clade diversification of Paris (Melanthiaceae).BMC Plant Biol.19:543. 10.1186/s12870-019-2147-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  JiY.YangC.HuangY. (2017). A new species of Paris sect. Axiparis (Melanthiaceae) from Yunnan, China.Phytotaxa306234–236. 10.11646/phytotaxa.306.3.6

  • CrossRef
  • Google Scholar

• 42

  JiY.YangJ.LandisJ. B.WangS.YangZ.ZhangY. (2021). Deciphering the taxonomic delimitation of Ottelia acuminata (Hydrocharitaceae) using complete plastomes as super-barcodes.Front. Plant Sci.12:681270. 10.3389/fpls.2021.681270

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  JinJ. J.YuW. B.YangJ. B.SongY.dePamphilisC. W.YiT. S.et al (2020). GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes.Genome Biol.21:241. 10.1186/s13059-020-02154-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  JirapatrasilpP.BackeljauT.PrasankokP.ChanabunR.PanhaS. (2019). Untangling a mess of worms: species delimitations reveal morphological crypsis and variability in Southeast Asian semi-aquatic earthworms (Almidae, Glyphidrilus).Mol. Phylogenet. Evol.139:106531. 10.1016/j.ympev.2019.106531

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  KaneN.SveinssonS.DempewolfH.YangJ. Y.ZhangD.EngelsJ. M.et al (2012). Ultra-barcoding in cacao (Theobroma spp.; Malvaceae) using whole chloroplast genomes and nuclear ribosomal DNA.Am. J. Bot.99320–329. 10.3732/ajb.1100570

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  KapliP.LutteroppS.ZhangJ.KobertK.PavlidisP.StamatakisA.et al (2017). Multi-rate Poisson tree processes for single-locus species delimitation under maximum likelihood and Markov chain Monte Carlo.Bioinformatics331630–1638. 10.1093/bioinformatics/btx025

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  KatohK.StandleyD. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability.Mol. Biol. Evol.30772–780. 10.1093/molbev/mst010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  KearseM.MoirR.WilsonA.Stones-HavasS.CheungM.SturrockS.et al (2012). Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data.Bioinformatics281647–1649. 10.1093/bioinformatics/bts199

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  KekkonenM.HebertP. D. (2014). DNA barcode-based delineation of putative species: efficient start for taxonomic workflows.Mol. Ecol. Resour.14706–715. 10.1111/1755-0998.12233

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  KnopeM. L.BellingerM. R.DatlofE. M.GallaherT. J.JohnsonM. A. (2020). Insights into the Evolutionary History of the Hawaiian Bidens (Asteraceae) Adaptive Radiation Revealed Through Phylogenomics.J. Hered.111119–137. 10.1093/jhered/esz066

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  KressW. J.WurdackK. J.ZimmerE. A.WeigtL. A.JanzenD. H. (2005). Use of DNA barcodes to identify flowering plants.Proc. Natl. Acad. Sci. U.S.A.1028369–8374. 10.1073/pnas.0503123102

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  LambertF.ClarkJ. L.JolyS. (2017). Species delimitation in the Caribbean Gesneria viridiflora complex (Gesneriaceae) reveals unsuspected endemism.Taxon661171–1183. 10.12705/665.10

  • CrossRef
  • Google Scholar

• 53

  LéveilléH. (1906). Nouvelles contributions à la connaissance des Liliacées, Amaryllidacées, Iridacées et Hemodoracées de Chine.Mem. Pontif. Accad. Rom. Nuovi Lincei24355–356.

  • Google Scholar

• 54

  LewinR. A. (1981). Three species concepts.Taxon30609–613. 10.2307/1219942

  • CrossRef
  • Google Scholar

• 55

  LiF. J.RousseauD. D.WuN. Q.HaoQ. Z.PeiY. P. (2008). Late Neogene evolution of the East Asian monsoon revealed by terrestrial mollusk record in western Chinese loess plateau: from winter to summer dominated subregime.Earth Planet Sci. Lett.274439–447. 10.1016/j.epsl.2008.07.038

  • CrossRef
  • Google Scholar

• 56

  LiH. (1984). Two new species of the genus Paris.Acta Bot. Yunnan.6273–276.

  • Google Scholar

• 57

  LiH. (1992). New taxa of the genus Paris L.Plant Divers.141–3.

  • Google Scholar

• 58

  LiH. (1998). The Genus Paris (Trilliaceae).Beijing: Science Press.

  • Google Scholar

• 59

  LiL.ZhangJ.LuZ. Q.ZhaoJ. L.LiQ. J. (2021). Genomic data reveal two distinct species from the widespread alpine ginger Roscoea tibetica Batalin (Zingiberaceae).J. Syst. Evol.591232–1243. 10.1111/jse.12596

  • CrossRef
  • Google Scholar

• 60

  LiX.YangY.HenryR. J.RossettoM.WangY.ChenS. (2015). Plant DNA barcoding: from gene to genome.Biol. Rev. Camb. Philos. Soc.90157–166. 10.1111/brv.12104

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  LoweT. M.ChanP. P. (2016). tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes.Nucleic Acids Res.44W54–W57. 10.1093/nar/gkw413

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  MaZ. Y.WenJ.TianJ. P.GuiL. L.LiuX. Q. (2020). Testing morphological trait evolution and assessing species delimitations in the grape genus using a phylogenomic framework.Mol. Phylogenet. Evol.148:106809. 10.1016/j.ympev.2020.106809

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63

  MaceG. M. (2004). The role of taxonomy in species conservation.Philos. Trans. R. Soc. Lond. B Biol. Sci.359711–719. 10.1098/rstb.2003.1454

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  MastersB. C.FanV.RossH. A. (2011). Species Delimitation-a Geneious plugin for the exploration of species boundaries.Mol. Ecol. Resour.11154–157. 10.1111/j.1755-0998.2010.02896.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 65

  MillerM. R.DunhamJ. P.AmoresA.CreskoW. A.JohnsonE. A. (2007). Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers.Genome Res.17240–248. 10.1101/gr.5681207

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  MutanenM.KekkonenM.ProsserS. W.HebertP. D.KailaL. (2015). One species in eight: DNA barcodes from type specimens resolve a taxonomic quagmire.Mol. Ecol. Resour.15967–984. 10.1111/1755-0998.12361

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67

  NockC. J.WatersD. L.EdwardsM. A.BowenS. G.RiceN.CordeiroG. M.et al (2011). Chloroplast genome sequences from total DNA for plant identification.Plant Biotechnol. J.9328–333. 10.1111/j.1467-7652.2010.00558.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  PellicerJ.KellyL. J.LeitchI. J.ZomleferW. B.FayM. F. (2014). A universe of dwarfs and giants: genome size and chromosome evolution in the monocot family Melanthiaceae.New Phytol.2011484–1497. 10.1111/nph.12617

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  PhilippeH.BrinkmannH.LavrovD. V.LittlewoodD. T. J.ManuelM.WörheideG.et al (2011). Resolving difficult phylogenetic questions: Why more sequences are not enough.PLoS Biol.9:e1000602. 10.1371/journal.pbio.1000602

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

  PonsJ.BarracloughT. G.Gomez-ZuritaJ.CardosoA.DuranD. P.HazellS.et al (2006). Sequence-based species delimitation for the DNA taxonomy of undescribed insects.Syst. Biol.55595–609. 10.1080/10635150600852011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71

  PosadaD.BuckleyT. R. (2004). Model selection and model averaging in phylogenetics: advantages of akaike information criterion and Bayesian approaches over likelihood ratio tests.Syst. Biol.53793–808. 10.1080/10635150490522304

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  PosadaD.CrandallK. A. (1998). MODELTEST: testing the model of DNA substitution.Bioinformatics14817–818. 10.1093/bioinformatics/14.9.817

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  PuillandreN.LambertA.BrouilletS.AchazG. (2011). ABGD, Automatic Barcode Gap Discovery for primary species delimitation.Mol. Ecol.211864–1877. 10.1111/j.1365-294X.2011.05239.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74

  PuillandreN.ModicaM. V.ZhangY.SirovichL.BoisselierM. C.CruaudC.et al (2012). Large-scale species delimitation method for hyperdiverse groups.Mol. Ecol.212671–2691. 10.1111/j.1365-294X.2012.05559.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75

  RannalaB.YangZ. (2013). Improved reversible jump algorithms for Bayesian species delimitation.Genetics194245–253. 10.1534/genetics.112.149039

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  RodrigoA.BertelsF.HeledJ.NoderR.ShearmanH.TsaiP. (2008). The perils of plenty: What are we going to do with all these genes?Philos. Trans. R. Soc. Lond. B Biol. Sci.3633893–3902. 10.1098/rstb.2008.0173

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 77

  RonquistF.TeslenkoM.Mark VanderP.AyresD. L.DarlingA.HöhnaS.et al (2012). MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space.Syst. Biol.61539–542. 10.1093/sysbio/sys029

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78

  RosenbergN. A. (2007). Statistical tests for taxonomic distinctiveness from observations of monophyly.Evolution61317–323. 10.1111/j.1558-5646.2007.00023.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 79

  RuhsamM.RaiH. S.MathewsS.RossT. G.GrahamS. W.RaubesonL. A.et al (2015). Does complete plastid genome sequencing improve species discrimination and phylogenetic resolution in Araucaria?Mol. Ecol. Resour.151067–1078. 10.1111/1755-0998.12375

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 80

  SitesJ. W.CrandallK. A. (1997). Testing Species Boundaries in Biodiversity Studies.Conserv. Biol.111289–1297. 10.1046/j.1523-1739.1997.96254.x

  • CrossRef
  • Google Scholar

• 81

  SitesJ. W.MarshallJ. C. (2003). Delimiting species: a Renaissance issue in systematic biology.Trends Ecol.18462–470. 10.1016/S0169-5347(03)00184-8

  • CrossRef
  • Google Scholar

• 82

  ŚlipikoM.MyszczyñskiK.BuczkowskaK.Ba̧czkiewiczA.SzczeciñskaM.SawickiJ. (2020). Molecular delimitation of European leafy liverworts of the genus Calypogeia based on plastid super-barcodes.BMC Plant Biol.20:243. 10.1186/s12870-020-02435-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 83

  SongC.LinX. L.WangQ.WangX. H. (2018). DNA barcodes successfully delimit morphospecies in a superdiverse insect genus.Zool. Scr.47311–324. 10.1111/zsc.12284

  • CrossRef
  • Google Scholar

• 84

  StamatakisA. (2006). RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models.Bioinformatics222688–2690. 10.1093/bioinformatics/btl446

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 85

  StraubS. C.ParksM.WeitemierK.FishbeinM.CronnR. C.ListonA. (2012). Navigating the tip of the genomic iceberg: next-generation sequencing for plant systematics.Am. J. Bot.99349–364. 10.3732/ajb.1100335

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 86

  SuN.LiuB. B.WangJ. R.TongR. C.RenC.ChangZ. Y.et al (2021). On the species delimitation of the Maddenia group of Prunus (Rosaceae): evidence from plastome and nuclear sequences and morphology.Front. Plant Sci.12:743643. 10.3389/fpls.2021.743643

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 87

  SuX.WuG.LiL.LiuJ. (2015). Species delimitation in plants using the Qinghai-Tibet Plateau endemic Orinus (Poaceae: Tridentinae) as an example.Ann. Bot.11635–48. 10.1093/aob/mcv062

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 88

  SukumaranJ.KnowlesL. L. (2017). Multispecies coalescent delimits structure, not species.Proc. Natl. Acad. Sci. U.S.A.1141607–1612. 10.1073/pnas.1607921114

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 89

  SwoffordD. L. (2002). PAUP: Phylogenetic Analysis Using Parsimony (and other methods), 4.0 Beta.Sunderland, MA: Sinauer Associates.

  • Google Scholar

• 90

  TakhtajanA. (1983). A revision of Daiswa (Trilliaceae).Brittonia35255–270. 10.2307/2806025

  • CrossRef
  • Google Scholar

• 91

  TillichM.LehwarkP.PellizzerT.Ulbricht-JonesE.FischerA.BockR.et al (2017). GeSeq- Versatile and accurate annotation of organelle genomes.Nucleic Acids Res.45W6–W11. 10.1093/nar/gkx391

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 92

  WileyE. O.MaydenR. (2000). “The evolutionary species concept,” in Species Concepts and Phylogenetic Theory: A Debate, edsWheelerQ. D.MeierR. (New York, NY: Columbia University Press), 70–89.

  • Google Scholar

• 93

  WuW.NgW. L.YangJ. X.LiW. M.GeX. J. (2018). High cryptic species diversity is revealed by genome-wide polymorphisms in a wild relative of banana, Musa itinerans, and implications for its conservation in subtropical China.BMC Plant Biol.18:194. 10.1186/s12870-018-1410-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 94

  XuZ.WeiN.TanY.PengS.NgumbauV. M.HuG.et al (2019). Paris lihengiana (Melanthiaceae: Parideae), a new species from Yunnan, China.Phytotaxa392045–053. 10.11646/phytotaxa.392.1.4

  • CrossRef
  • Google Scholar

• 95

  YangZ.RannalaB. (2010). Bayesian species delimitation using multilocus sequence data.Proc. Natl. Acad. Sci. U.S.A.1079264–9269. 10.1073/pnas.0913022107

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 96

  YangZ.YangC.JiY. (2019). Paris variabilis (Melanthiaceae), A new species from southwestern China.Phytotaxa401190–198. 10.11646/phytotaxa.401.3.4

  • CrossRef
  • Google Scholar

• 97

  YaoY. F.BruchA. A.MosbruggerV.LiC. S. (2011). Quantitative reconstruction of Miocene climate patterns and evolution in southern China based on plant fossils.Palaeogeogr. Palaeoclimatol. Palaeoecol.304291–307. 10.1016/j.palaeo.2010.04.012

  • CrossRef
  • Google Scholar

• 98

  ZhangJ.KapliP.PavlidisP.StamatakisA. (2013). A general species delimitation method with applications to phylogenetic placements.Bioinformatics292869–2876. 10.1093/bioinformatics/btt499

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 99

  ZhouM.YangG.SunG.GuoZ.GongX.PanY. (2020). Resolving complicated relationships of the Panax bipinnatifidus complex in southwestern China by RAD-seq data.Mol. Phylogenet. Evol.149:106851. 10.1016/j.ympev.2020.106851

  • Pubmed Abstract
  • CrossRef
  • Google Scholar