Yunheng Ji1,2*
Jin Yang1,3
Jacob B. Landis4,5Shuying Wang1,3Lei Jin1,6Pingxuan Xie1,6Haiyang Liu7
Jun-Bo Yang8*
Ting-Shuang Yi8*
- 1CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
- 2Yunnan Key Laboratory for Integrative Conservation of Plant Species With Extremely Small Population, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
- 3School of Life Sciences, Yunnan University, Kunming, China
- 4Section of Plant Biology and the L. H. Bailey Hortorium, School of Integrative Plant Science, Cornell University, Ithaca, NY, United States
- 5BTI Computational Biology Center, Boyce Thompson Institute, Ithaca, NY, United States
- 6School of Traditional Chinese Medicine, Guangdong Pharmaceutical University, Guangzhou, China
- 7State Key Laboratory of Phytochemistry and Plant Resources in West China, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
- 8Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China
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. Nine nominal species within Paris section Axiparis, showing their leaf size and shape variations. Paris axialis (A), P. dulongensis (B), P. forrestii (C), P. guizhouensis (D), P. lihengiana (E), P. rugosa (F), P. tenchongensis (G), P. vaniotii (H), and P. variabilis (I).
Figure 2. Distribution of the nine nominal species within Paris section Axiparis.
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.
Table 1. Plant samples used in this study with voucher information and GenBank accession numbers.
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,1 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 server2 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 server3 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. Phylogenetic tree of Paris section Axiparis based on Maximum likelihood (ML, A) and Bayesian inference (BI, B) analyses of complete plastomes. Numbers above branches indicate ML bootstrap (BS) percentages and Bayesian posterior probabilities (PP). Species delimitation schemes proposed by automatic barcode gap discovery (ABGD), multi-rate Poisson tree processes model (mPTP), and species delimitation plugin (SDP) of Geneious are reflected on the tree topology.
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. Phylogenetic tree of Paris section Axiparis based on Maximum likelihood (ML, A) and Bayesian inference (BI, B) analyses of nuclear ribosomal DNA (nrDNA). Numbers above branches indicate ML bootstrap (BS) percentages and Bayesian posterior probabilities (PP). Species delimitation schemes proposed by automatic barcode gap discovery (ABGD), multi-rate Poisson tree processes model (mPTP), and species delimitation plugin (SDP) of Geneious are reflected on the tree topology.
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).
Table 2. The number of putative species recognized by Automatic Barcode Gap Discovery (ABGD) analyses among 51 complete plastomes using three distance metrics.
Table 3. The number of putative species recognized by Automatic Barcode Gap Discovery (ABGD) analyses among 51 nrDNA sequences using three distance metrics.
Table 4. Posterior delimitation probability of two putative species proposed by mPTP analyses of complete plastomes and nuclear ribosomal DNA (nrDNA) sequences.
Table 5. Species delimitation plugin (SDP) analyses show the distinctiveness [Rodrigo’s P(RD)] and significance [Rosenberg’s P(AB)] of the two clades recovered by phylogenetic analyses of complete plastomes and nuclear ribosomal DNA (nrDNA) as species-level taxonomic units.
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.
Data Availability Statement
The datasets presented in this study can be found in online repositories. All the relevant accession numbers are included in Supplementary Table 1.
Author Contributions
YJ, J-BY, and T-SY conceived and designed the research framework. YJ, JY, and LJ collected samples. JY, J-BY, SW, LJ, and PX collected and analyzed the data. YJ wrote the original draft manuscript. JL, T-SY, and HL revised and edited the final manuscript. All authors have read and agreed to the published version of the manuscript.
Funding
The present study was granted 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 a grant of the Large-scale Scientific Facilities of the Chinese Academy of Sciences (No. 2017-LSF-GBOWS-02).
Conflict of Interest
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.
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.
Acknowledgments
We are grateful to Chengjin Yang, Guohua Zhou, Ren Zhao, Tingzhou Zhao, Zhangming Wang, Changming Yang, Yongxing Xiong, Guohai He, Yulong Li, Fengtian Lin, and Qiang Zhao for their help in sample collection.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.832034/full#supplementary-material
Footnotes
- ^ http://purl.org/phylo/treebase/phylows/study/TB2:S29422
- ^ https://bioinfo.mnhn.fr/abi/public/abgd/abgdweb.html; the latest access date: July 9, 2020.
- ^ http://species.h-its.org/ptp/; the latest access date: August 14, 2020.
References
Aguilar, C., Wood, P. L., Cusi, J. C., Guzmán, A., Huari, F., Lundberg, M., 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. Zookeys 364, 47–91. doi: 10.3897/zookeys.364.6109
An, Z. S., Kutzbach, J. E., Prell, W. L., and Porter, S. C. (2001). Evolution of Asian monsoons and phased uplift of the Himalaya–Tibetan plateau since Late Miocene times. Nature 411, 62–66. doi: 10.1038/35075035
Baird, N. A., Etter, P. D., Atwood, T. S., Currey, M. C., Shiver, A. L., Lewis, Z. A., et al. (2008). Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS One 3:e3376. doi: 10.1371/journal.pone.0003376
Bickford, D., Lohman, D. J., Sodhi, N. S., Ng, P. K., Meier, R., Winker, K., et al. (2007). Cryptic species as a window on diversity and conservation. Trends Ecol. Evol. 22, 148–155. doi: 10.1016/j.tree.2006.11.004
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170
Brower, A. V. Z., DeSalle, R., and Vogler, A. (1996). Gene trees, species trees and systematics: a cladistic perspective. Annu. Rev. Ecol. Syst. 27, 423–450. doi: 10.1146/annurev.ecolsys.27.1.423
Camargo, A., Morando, M., Avila, L. J., and Sites, J. 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). Evolution 66, 2834–2849. doi: 10.1111/j.1558-5646.2012.01640.x
Carstens, B. C., and Satler, J. D. (2013). The carnivorous plant described as Sarracenia alata contains two cryptic species. Biol. J. Linn. Soc. 109, 737–746. doi: 10.1111/bij.12093
Cheng, Z., Shu, H., Zhang, S., Luo, B., Gu, R., Zhang, R., et al. (2020). From folk taxonomy to species confirmation of Acorus (Acoraceae): evidences based on phylogenetic and metabolomic analyses. Front. Plant. Sci. 11:965. doi: 10.3389/fpls.2020.00965
Coissac, E., Hollingsworth, P. M., Lavergne, S., and Taberlet, P. (2016). From barcodes to genomes: extending the concept of DNA barcoding. Mol. Ecol. 25, 1423–1428. doi: 10.1111/mec.13549
De Queiroz, K. (2007). Species concepts and species delimitation. Syst. Bot. 56, 879–886. doi: 10.1080/10635150701701083
Dodsworth, S. (2015). Genome skimming for next-generation biodiversity analysis. Trends Plant Sci. 20, 525–527. doi: 10.1016/j.tplants.2015.06.012
Donkpegan, A. S. L., Doucet, J. L., Hardy, O. J., Heuertz, M., and Piñeiro, R. (2020). Miocene diversification in the Savannahs precedes tetraploid rainforest radiation in the African tree genus Afzelia (Detarioideae. Fabaceae). Front. Plant Sci. 11:798. doi: 10.3389/fpls.2020.00798
Doyle, J. J., and Doyle, J. (1987). A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem. Bull. 19, 11–15.
Duminil, J., Kenfack, D., Viscosi, V., Grumiau, L., and Hardy, O. J. (2012). Testing species delimitation in sympatric species complexes: the case of an African tropical tree. Carapa spp. (Meliaceae). Mol. Phylogenet. Evol. 62, 275–285. doi: 10.1016/j.ympev.2011.09.020
Duminil, J., and Michele, M. D. (2009). Plant species delimitation: a comparison of morphological and molecular markers. Plant Biosyst. 143, 528–542. doi: 10.1080/11263500902722964
Eisenring, M., Altermatt, F., Westram, A. M., and Jokela, J. (2016). Habitat requirements and ecological niche of two cryptic amphipod species at landscape and local scales. Ecosphere 7:e01319. doi: 10.1002/ecs2.1319
Ellis, J. S., Knight, M. E., Carvell, C., and Goulson, D. (2006). Cryptic species identification: a simple diagnostic tool for discriminating between two problematic bumblebee species. Mol. Ecol. Notes 6, 540–542. doi: 10.1111/j.1471-8286.2006.01231.x
Ence, D. D., and Carstens, B. C. (2011). SpedeSTEM: a rapid and accurate method for species delimitation. Mol. Ecol. Resour. 11, 473–480. doi: 10.1111/j.1755-0998.2010.02947.x
Farris, J.S., Källersjö, M., Kluge, A.C., and Bult, C. (1994). Testing significance of incongruence. Cladistics 10, 315–319. doi: 10.1111/j.1096-0031.1994.tb00181.x
Firetti, F., Zuntini, A. R., Gaiarsa, J. W., Oliveira, R. S., Lohmann, L. G., and Van Sluys, M. A. (2017). Complete chloroplast genome sequences contribute to plant species delimitation: a case study of the Anemopaegma species complex. Am. J. Bot. 104, 1493–1509. doi: 10.3732/ajb.1700302
Fu, C. N., Mo, Z. Q., Yang, J. B., Cai, J., Ye, L. J., Zou, J. Y., et al. (2021). Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendron. Mol. Ecol. Resour. 22, 404–414. doi: 10.1111/1755-0998.13479
Fu, C. N., Wu, C. S., Ye, L. J., Mo, Z. Q., Liu, J., Chang, Y. W., et al. (2019). Prevalence of isomeric plastomes and effectiveness of plastome super-barcodes in yews (Taxus) worldwide. Sci. Rep. 9:2773. doi: 10.1038/s41598-019-39161-x
Fujita, M. K., Leaché, A. D., Burbrink, F. T., McGuire, J. A., and Moritz, C. (2012). Coalescent-based species delimitation in an integrative taxonomy. Trends Ecol. Evol. 27, 480–488. doi: 10.1016/j.tree.2012.04.012
Hebert, P. D., Cywinska, A., Ball, S. L., and De-Waard, J. R. (2003). Biological identifications through DNA barcodes. Proc. Biol. Sci. 270, 313–321. doi: 10.1098/rspb.2002.2218
Hebert, P. D., Penton, E. H., Burns, J. M., Janzen, D. H., and Hallwachs, W. (2004). Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proc. Natl. Acad. Sci. U.S.A. 101, 14812–14817. doi: 10.1073/pnas.0406166101
Hollingsworth, P. M. (2011). Refining the DNA barcode for land plants. Proc. Natl. Acad. Sci. U.S.A. 108, 19451–19452. doi: 10.1073/pnas.1116812108
Hollingsworth, P. M., Forrest, L. L., Spouge, J. L., Hajibabaei, M., Ratnasingham, S., Van der Bank, M., et al. (2009). A DNA barcode for land plants. Proc. Natl. Acad. Sci. U.S.A. 106, 12794–12797. doi: 10.1073/pnas.0905845106
Hollingsworth, P. M., Graham, S. W., and Little, D. P. (2011). Choosing and using a plant DNA barcode. PLoS One 6:e19254. doi: 10.1371/journal.pone.0019254
Hollingsworth, P. M., Li, D. Z., Van der Bank, M., and Twyford, A. D. (2016). Telling plant species apart with DNA: from barcodes to genomes. Philos. Trans. R. Soc. B Biol. Sci. 371:20150338. doi: 10.1098/rstb.2015.0338
Huang, W., Xie, X., Huo, L., Liang, X., Wang, X., and Chen, X. (2020). An integrative DNA barcoding framework of ladybird beetles (Coleoptera: Coccinellidae). Sci. Rep. 10:10063. doi: 10.1038/s41598-020-66874-1
Huang, Y., Li, X., Yang, Z., Yang, C., Yang, J., and Ji, Y. (2016). Analysis of complete chloroplast genome sequences improves phylogenetic resolution in Paris (Melanthiaceae). Front. Plant Sci. 7:1797. doi: 10.3389/fpls.2016.01797
Jacobs, S. J., Kristofferson, C., Uribe-Convers, S., Latvis, M., and Tank, D. C. (2018). Incongruence in molecular species delimitation schemes: What to do when adding more data is difficult. Mol. Ecol. 27, 2397–2413. doi: 10.1111/mec.14590
Ji, Y., Fritsch, P. W., Li, H., Xiao, T., and Zhou, Z. (2006). Phylogeny and classification of Paris (Melanthiaceae) inferred from DNA sequence data. Ann. Bot. 98, 245–256. doi: 10.1093/aob/mcl095
Ji, Y., Liu, C., Yang, J., Jin, L., Yang, Z., and Yang, J. B. (2020). Ultra-barcoding discovers a cryptic species in Paris yunnanensis (Melanthiaceae), a medicinally important plant. Front. Plant Sci. 11:411. doi: 10.3389/fpls.2020.00411
Ji, Y., Liu, C., Yang, Z., Yang, L., He, Z., Wang, H., et al. (2019a). Testing and using complete plastomes and ribosomal DNA sequences as the next generation DNA barcodes in Panax (Araliaceae). Mol. Ecol. Resour. 19, 1333–1345. doi: 10.1111/1755-0998.13050
Ji, Y., Yang, L., Chase, M. W., Liu, C., Yang, Z., Yang, J., et al. (2019b). Plastome phylogenomics, biogeography, and clade diversification of Paris (Melanthiaceae). BMC Plant Biol. 19:543. doi: 10.1186/s12870-019-2147-6
Ji, Y., Yang, C., and Huang, Y. (2017). A new species of Paris sect. Axiparis (Melanthiaceae) from Yunnan, China. Phytotaxa 306, 234–236. doi: 10.11646/phytotaxa.306.3.6
Ji, Y., Yang, J., Landis, J. B., Wang, S., Yang, Z., and Zhang, Y. (2021). Deciphering the taxonomic delimitation of Ottelia acuminata (Hydrocharitaceae) using complete plastomes as super-barcodes. Front. Plant Sci. 12:681270. doi: 10.3389/fpls.2021.681270
Jin, J. J., Yu, W. B., Yang, J. B., Song, Y., dePamphilis, C. W., Yi, T. S., et al. (2020). GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol. 21:241. doi: 10.1186/s13059-020-02154-5
Jirapatrasilp, P., Backeljau, T., Prasankok, P., Chanabun, R., and Panha, S. (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. doi: 10.1016/j.ympev.2019.106531
Kane, N., Sveinsson, S., Dempewolf, H., Yang, J. Y., Zhang, D., Engels, J. M., et al. (2012). Ultra-barcoding in cacao (Theobroma spp.; Malvaceae) using whole chloroplast genomes and nuclear ribosomal DNA. Am. J. Bot. 99, 320–329. doi: 10.3732/ajb.1100570
Kapli, P., Lutteropp, S., Zhang, J., Kobert, K., Pavlidis, P., Stamatakis, A., et al. (2017). Multi-rate Poisson tree processes for single-locus species delimitation under maximum likelihood and Markov chain Monte Carlo. Bioinformatics 33, 1630–1638. doi: 10.1093/bioinformatics/btx025
Katoh, K., and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010
Kearse, M., Moir, R., Wilson, A., Stones-Havas, S., Cheung, M., Sturrock, S., et al. (2012). Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649. doi: 10.1093/bioinformatics/bts199
Kekkonen, M., and Hebert, P. D. (2014). DNA barcode-based delineation of putative species: efficient start for taxonomic workflows. Mol. Ecol. Resour. 14, 706–715. doi: 10.1111/1755-0998.12233
Knope, M. L., Bellinger, M. R., Datlof, E. M., Gallaher, T. J., and Johnson, M. A. (2020). Insights into the Evolutionary History of the Hawaiian Bidens (Asteraceae) Adaptive Radiation Revealed Through Phylogenomics. J. Hered. 111, 119–137. doi: 10.1093/jhered/esz066
Kress, W. J., Wurdack, K. J., Zimmer, E. A., Weigt, L. A., and Janzen, D. H. (2005). Use of DNA barcodes to identify flowering plants. Proc. Natl. Acad. Sci. U.S.A. 102, 8369–8374. doi: 10.1073/pnas.0503123102
Lambert, F., Clark, J. L., and Joly, S. (2017). Species delimitation in the Caribbean Gesneria viridiflora complex (Gesneriaceae) reveals unsuspected endemism. Taxon 66, 1171–1183. doi: 10.12705/665.10
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 Lincei 24, 355–356.
Li, F. J., Rousseau, D. D., Wu, N. Q., Hao, Q. Z., and Pei, Y. 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. 274, 439–447. doi: 10.1016/j.epsl.2008.07.038
Li, L., Zhang, J., Lu, Z. Q., Zhao, J. L., and Li, Q. J. (2021). Genomic data reveal two distinct species from the widespread alpine ginger Roscoea tibetica Batalin (Zingiberaceae). J. Syst. Evol. 59, 1232–1243. doi: 10.1111/jse.12596
Li, X., Yang, Y., Henry, R. J., Rossetto, M., Wang, Y., and Chen, S. (2015). Plant DNA barcoding: from gene to genome. Biol. Rev. Camb. Philos. Soc. 90, 157–166. doi: 10.1111/brv.12104
Lowe, T. M., and Chan, P. P. (2016). tRNAscan-SE On-line: integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 44, W54–W57. doi: 10.1093/nar/gkw413
Ma, Z. Y., Wen, J., Tian, J. P., Gui, L. L., and Liu, X. Q. (2020). Testing morphological trait evolution and assessing species delimitations in the grape genus using a phylogenomic framework. Mol. Phylogenet. Evol. 148:106809. doi: 10.1016/j.ympev.2020.106809
Mace, G. M. (2004). The role of taxonomy in species conservation. Philos. Trans. R. Soc. Lond. B Biol. Sci. 359, 711–719. doi: 10.1098/rstb.2003.1454
Masters, B. C., Fan, V., and Ross, H. A. (2011). Species Delimitation-a Geneious plugin for the exploration of species boundaries. Mol. Ecol. Resour. 11, 154–157. doi: 10.1111/j.1755-0998.2010.02896.x
Miller, M. R., Dunham, J. P., Amores, A., Cresko, W. A., and Johnson, E. A. (2007). Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers. Genome Res. 17, 240–248. doi: 10.1101/gr.5681207
Mutanen, M., Kekkonen, M., Prosser, S. W., Hebert, P. D., and Kaila, L. (2015). One species in eight: DNA barcodes from type specimens resolve a taxonomic quagmire. Mol. Ecol. Resour. 15, 967–984. doi: 10.1111/1755-0998.12361
Nock, C. J., Waters, D. L., Edwards, M. A., Bowen, S. G., Rice, N., Cordeiro, G. M., et al. (2011). Chloroplast genome sequences from total DNA for plant identification. Plant Biotechnol. J. 9, 328–333. doi: 10.1111/j.1467-7652.2010.00558.x
Pellicer, J., Kelly, L. J., Leitch, I. J., Zomlefer, W. B., and Fay, M. F. (2014). A universe of dwarfs and giants: genome size and chromosome evolution in the monocot family Melanthiaceae. New Phytol. 201, 1484–1497. doi: 10.1111/nph.12617
Philippe, H., Brinkmann, H., Lavrov, D. V., Littlewood, D. T. J., Manuel, M., Wörheide, G., et al. (2011). Resolving difficult phylogenetic questions: Why more sequences are not enough. PLoS Biol. 9:e1000602. doi: 10.1371/journal.pbio.1000602
Pons, J., Barraclough, T. G., Gomez-Zurita, J., Cardoso, A., Duran, D. P., Hazell, S., et al. (2006). Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst. Biol. 55, 595–609. doi: 10.1080/10635150600852011
Posada, D., and Buckley, T. R. (2004). Model selection and model averaging in phylogenetics: advantages of akaike information criterion and Bayesian approaches over likelihood ratio tests. Syst. Biol. 53, 793–808. doi: 10.1080/10635150490522304
Posada, D., and Crandall, K. A. (1998). MODELTEST: testing the model of DNA substitution. Bioinformatics 14, 817–818. doi: 10.1093/bioinformatics/14.9.817
Puillandre, N., Lambert, A., Brouillet, S., and Achaz, G. (2011). ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol. Ecol. 21, 1864–1877. doi: 10.1111/j.1365-294X.2011.05239.x
Puillandre, N., Modica, M. V., Zhang, Y., Sirovich, L., Boisselier, M. C., Cruaud, C., et al. (2012). Large-scale species delimitation method for hyperdiverse groups. Mol. Ecol. 21, 2671–2691. doi: 10.1111/j.1365-294X.2012.05559.x
Rannala, B., and Yang, Z. (2013). Improved reversible jump algorithms for Bayesian species delimitation. Genetics 194, 245–253. doi: 10.1534/genetics.112.149039
Rodrigo, A., Bertels, F., Heled, J., Noder, R., Shearman, H., and Tsai, P. (2008). The perils of plenty: What are we going to do with all these genes? Philos. Trans. R. Soc. Lond. B Biol. Sci. 363, 3893–3902. doi: 10.1098/rstb.2008.0173
Ronquist, F., Teslenko, M., Mark Vander, P., Ayres, D. L., Darling, A., Höhna, S., et al. (2012). MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 61, 539–542. doi: 10.1093/sysbio/sys029
Rosenberg, N. A. (2007). Statistical tests for taxonomic distinctiveness from observations of monophyly. Evolution 61, 317–323. doi: 10.1111/j.1558-5646.2007.00023.x
Ruhsam, M., Rai, H. S., Mathews, S., Ross, T. G., Graham, S. W., Raubeson, L. A., et al. (2015). Does complete plastid genome sequencing improve species discrimination and phylogenetic resolution in Araucaria? Mol. Ecol. Resour. 15, 1067–1078. doi: 10.1111/1755-0998.12375
Sites, J. W., and Crandall, K. A. (1997). Testing Species Boundaries in Biodiversity Studies. Conserv. Biol. 11, 1289–1297. doi: 10.1046/j.1523-1739.1997.96254.x
Sites, J. W., and Marshall, J. C. (2003). Delimiting species: a Renaissance issue in systematic biology. Trends Ecol. 18, 462–470. doi: 10.1016/S0169-5347(03)00184-8
Ślipiko, M., Myszczyñski, K., Buczkowska, K., Ba̧czkiewicz, A., Szczeciñska, M., and Sawicki, J. (2020). Molecular delimitation of European leafy liverworts of the genus Calypogeia based on plastid super-barcodes. BMC Plant Biol. 20:243. doi: 10.1186/s12870-020-02435-y
Song, C., Lin, X. L., Wang, Q., and Wang, X. H. (2018). DNA barcodes successfully delimit morphospecies in a superdiverse insect genus. Zool. Scr. 47, 311–324. doi: 10.1111/zsc.12284
Stamatakis, A. (2006). RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22, 2688–2690. doi: 10.1093/bioinformatics/btl446
Straub, S. C., Parks, M., Weitemier, K., Fishbein, M., Cronn, R. C., and Liston, A. (2012). Navigating the tip of the genomic iceberg: next-generation sequencing for plant systematics. Am. J. Bot. 99, 349–364. doi: 10.3732/ajb.1100335
Su, N., Liu, B. B., Wang, J. R., Tong, R. C., Ren, C., Chang, Z. 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. doi: 10.3389/fpls.2021.743643
Su, X., Wu, G., Li, L., and Liu, J. (2015). Species delimitation in plants using the Qinghai-Tibet Plateau endemic Orinus (Poaceae: Tridentinae) as an example. Ann. Bot. 116, 35–48. doi: 10.1093/aob/mcv062
Sukumaran, J., and Knowles, L. L. (2017). Multispecies coalescent delimits structure, not species. Proc. Natl. Acad. Sci. U.S.A. 114, 1607–1612. doi: 10.1073/pnas.1607921114
Swofford, D. L. (2002). PAUP: Phylogenetic Analysis Using Parsimony (and other methods), 4.0 Beta. Sunderland, MA: Sinauer Associates.
Takhtajan, A. (1983). A revision of Daiswa (Trilliaceae). Brittonia 35, 255–270. doi: 10.2307/2806025
Tillich, M., Lehwark, P., Pellizzer, T., Ulbricht-Jones, E., Fischer, A., Bock, R., et al. (2017). GeSeq- Versatile and accurate annotation of organelle genomes. Nucleic Acids Res. 45, W6–W11. doi: 10.1093/nar/gkx391
Wiley, E. O., and Mayden, R. (2000). “The evolutionary species concept,” in Species Concepts and Phylogenetic Theory: A Debate, eds Q. D. Wheeler and R. Meier (New York, NY: Columbia University Press), 70–89.
Wu, W., Ng, W. L., Yang, J. X., Li, W. M., and Ge, X. 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. doi: 10.1186/s12870-018-1410-6
Xu, Z., Wei, N., Tan, Y., Peng, S., Ngumbau, V. M., Hu, G., et al. (2019). Paris lihengiana (Melanthiaceae: Parideae), a new species from Yunnan, China. Phytotaxa 392, 045–053. doi: 10.11646/phytotaxa.392.1.4
Yang, Z., and Rannala, B. (2010). Bayesian species delimitation using multilocus sequence data. Proc. Natl. Acad. Sci. U.S.A. 107, 9264–9269. doi: 10.1073/pnas.0913022107
Yang, Z., Yang, C., and Ji, Y. (2019). Paris variabilis (Melanthiaceae), A new species from southwestern China. Phytotaxa 401, 190–198. doi: 10.11646/phytotaxa.401.3.4
Yao, Y. F., Bruch, A. A., Mosbrugger, V., and Li, C. S. (2011). Quantitative reconstruction of Miocene climate patterns and evolution in southern China based on plant fossils. Palaeogeogr. Palaeoclimatol. Palaeoecol. 304, 291–307. doi: 10.1016/j.palaeo.2010.04.012
Zhang, J., Kapli, P., Pavlidis, P., and Stamatakis, A. (2013). A general species delimitation method with applications to phylogenetic placements. Bioinformatics 29, 2869–2876. doi: 10.1093/bioinformatics/btt499
Keywords: plastome, ribosomal DNA, taxonomy, Paris Linn., species delimitation
Citation: Ji Y, Yang J, Landis JB, Wang S, Jin L, Xie P, Liu H, Yang J-B and Yi T-S (2022) Genome Skimming Contributes to Clarifying Species Limits in Paris Section Axiparis (Melanthiaceae). Front. Plant Sci. 13:832034. doi: 10.3389/fpls.2022.832034
Received: 09 December 2021; Accepted: 22 February 2022;
Published: 04 April 2022.
Edited by:
Carl J. Rothfels, University of California, Berkeley, United StatesReviewed by:
Ran Wei, Institute of Botany (CAS), ChinaMaría Fernanda Torres Jimenez, Vilnius University, Lithuania
Copyright © 2022 Ji, Yang, Landis, Wang, Jin, Xie, Liu, Yang and Yi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yunheng Ji, jiyh@mail.kib.ac.cn; Jun-Bo Yang, jbyang@mail.kib.ac.cn; Ting-Shuang Yi, tingshuangyi@mail.kib.ac.cn