Historical biogeography and evolutionary diversification of Lilium (Liliaceae): New insights from plastome phylogenomics.

Here, we infer the historical biogeography and evolutionary diversification of the genus Lilium. For this purpose, we used the complete plastomes of 64 currently accepted species in the genus Lilium (14 plastomes were newly sequenced) to recover the phylogenetic backbone of the genus and a time-calibrated phylogenetic framework to estimate biogeographical history scenarios and evolutionary diversification rates of Lilium. Our results suggest that ancient climatic changes and geological tectonic activities jointly shaped the distribution range and drove evolutionary radiation of Lilium, including the Middle Miocene Climate Optimum (MMCO), the late Miocene global cooling, as well as the successive uplift of the Qinghai-Tibet Plateau (QTP) and the strengthening of the monsoon climate in East Asia during the late Miocene and the Pliocene. This case study suggests that the unique geological and climatic events in the Neogene of East Asia, in particular the uplift of QTP and the enhancement of monsoonal climate, may have played an essential role in formation of uneven distribution of plant diversity in the Northern Hemisphere.

Comparative and phylogenetic analyses of Loranthaceae plastomes provide insights into the evolutionary trajectories of plastome degradation in hemiparasitic plants.

BACKGROUND: The lifestyle transition from autotrophy to heterotrophy often leads to extensive degradation of plastomes in parasitic plants, while the evolutionary trajectories of plastome degradation associated with parasitism in hemiparasitic plants remain poorly understood. In this study, phylogeny-oriented comparative analyses were conducted to investigate whether obligate Loranthaceae stem-parasites experienced higher degrees of plastome degradation than closely related facultative root-parasites and to explore the potential evolutionary events that triggered the 'domino effect' in plastome degradation of hemiparasitic plants. RESULTS: Through phylogeny-oriented comparative analyses, the results indicate that Loranthaceae hemiparasites have undergone varying degrees of plastome degradation as they evolved towards a heterotrophic lifestyle. Compared to closely related facultative root-parasites, all obligate stem-parasites exhibited an elevated degree plastome degradation, characterized by increased downsizing, gene loss, and pseudogenization, thereby providing empirical evidence supporting the theoretical expectation that evolution from facultative parasitism to obligate parasitism may result in a higher degree of plastome degradation in hemiparasites. Along with infra-familial divergence in Loranthaceae, several lineage-specific gene loss/pseudogenization events occurred at deep nodes, whereas further independent gene loss/pseudogenization events were observed in shallow branches. CONCLUSIONS: The findings suggest that in addition to the increasing levels of nutritional reliance on host plants, cladogenesis can be considered as another pivotal evolutionary event triggering the 'domino effect' in plastome degradation of hemiparasitic plants. These findings provide new insights into the evolutionary trajectory of plastome degradation in hemiparasitic plants.

Plastid phylogenomics and fossil evidence provide new insights into the evolutionary complexity of the 'woody clade' in Saxifragales.

BACKGROUND: The "woody clade" in Saxifragales (WCS), encompassing four woody families (Altingiaceae, Cercidiphyllaceae, Daphniphyllaceae, and Hamamelidaceae), is a phylogenetically recalcitrant node in the angiosperm tree of life, as the interfamilial relationships of the WCS remain contentious. Based on a comprehensive sampling of WCS genera, this study aims to recover a robust maternal backbone phylogeny of the WCS by analyzing plastid genome (plastome) sequence data using Bayesian inference (BI), maximum likelihood (ML), and maximum parsimony (MP) methods, and to explore the possible causes of the phylogenetic recalcitrance with respect to deep relationships within the WCS, in combination with molecular and fossil evidence. RESULTS: Although the four WCS families were identically resolved as monophyletic, the MP analysis recovered different tree topologies for the relationships among Altingiaceae, Cercidiphyllaceae, and Daphniphyllaceae from the ML and BI phylogenies. The fossil-calibrated plastome phylogeny showed that the WCS underwent a rapid divergence of crown groups in the early Cretaceous (between 104.79 and 100.23 Ma), leading to the origin of the stem lineage ancestors of Altingiaceae, Cercidiphyllaceae, Daphniphyllaceae, and Hamamelidaceae within a very short time span (∼4.56 Ma). Compared with the tree topology recovered in a previous study based on nuclear genome data, cytonuclear discordance regarding the interfamilial relationships of the WCS was detected. CONCLUSIONS: Molecular and fossil evidence imply that the early divergence of the WCS might have experienced radiative diversification of crown groups, extensive extinctions at the genus and species levels around the Cretaceous/Paleocene boundary, and ancient hybridization. Such evolutionarily complex events may introduce biases in topological estimations within the WCS due to incomplete lineage sorting, cytonuclear discordance, and long-branch attraction, potentially impacting the accurate reconstruction of deep relationships.

Unraveling the serial glycosylation in the biosynthesis of steroidal saponins in the medicinal plant Paris polyphylla and their antifungal action.

Sugar-sugar glycosyltransferases play important roles in constructing complex and bioactive saponins. Here, we characterized a series of UDP-glycosyltransferases responsible for biosynthesizing the branched sugar chain of bioactive steroidal saponins from a widely known medicinal plant Paris polyphylla var. yunnanensis. Among them, a 2'-O-rhamnosyltransferase and three 6'-O-glucosyltrasferases catalyzed a cascade of glycosylation to produce steroidal diglycosides and triglycosides, respectively. These UDP-glycosyltransferases showed astonishing substrate promiscuity, resulting in the generation of a panel of 24 terpenoid glycosides including 15 previously undescribed compounds. A mutant library containing 44 variants was constructed based on the identification of critical residues by molecular docking simulations and protein model alignments, and a mutant UGT91AH1Y187A with increased catalytic efficiency was obtained. The steroidal saponins exhibited remarkable antifungal activity against four widespread strains of human pathogenic fungi attributed to ergosterol-dependent damage of fungal cell membranes, and 2'-O-rhamnosylation appeared to correlate with strong antifungal effects. The findings elucidated the biosynthetic machinery for their production of steroidal saponins and revealed their potential as new antifungal agents.

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

BACKGROUND: Paris yunnanensis (Melanthiaceae) is a traditional Chinese medicinal plant of significant pharmaceutical importance. Due to previous taxonomic confusion, a congeneric species, Paris liiana, has been mistaken for P. yunnanensis and cultivated on a large scale, leading to the mixing of commercial products (i.e., seedlings and processed rhizomes) of P. yunnanensis with those of P. liiana. This may have adverse effects on quality control in the standardization of P. yunnanensis productions. As the lack of PCR amplifiable genomic DNA within processed rhizomes is an intractable obstacle to the authentication of P. yunnanensis products using PCR-based diagnostic tools, this study aimed to develop a PCR-free method to authenticate commercial P. yunnanensis products, by applying genome skimming to generate complete plastomes and nrDNA arrays for use as the molecular tags. RESULTS: Based on a dense intraspecies sampling of P. liiana and P. yunnanensis, the robustness of the proposed authentication systems was evaluated by phylogenetic inferences and experimental authentication of commercial seedling and processed rhizome samples. The results indicate that the genetic criteria of both complete plastomes and nrDNA arrays were consistent with the species boundaries to achieve accurate discrimination of P. yunnanensis and P. liinna. Owing to its desirable accuracy and sensitivity, genome skimming can serve as an effective and sensitive tool for monitoring and controlling the trade of P. yunnanensis products. CONCLUSION: This study provides a new way to solve the long-standing problem of the molecular authentication of processed plant products due to the lack of PCR amplifiable genomic DNA. The proposed authentication system will support quality control in the standardization of P. yunnanensis products in cultivation and drug production. This study also provides molecular evidence to clarify the long-standing taxonomic confusion regarding the species delimitation of P. yunnanensis, which will contribute to the rational exploration and conservation of the species.

Miocene Dipteronia (Sapindaceae) samaras from South Korea and their biogeographical implications.

Dipteronia, now endemic to East Asia, was widely distributed in North America during the Paleogene; however, its fossil records in Asia are scarce and none are of the Neogene. Here, we report the first Neogene Dipteronia samaras from South Korea. The more complete fossil records suggest that Dipteronia possibly originated in either Asia or North America and that its two known lineages have different geographical histories. The Dipteronia sinensis lineage was established in Asia and North America in the Paleocene and reached its maximum range in the Eocene, followed by stepwise range contraction and extirpation in North America, South Korea, and southwestern China, finally becoming endemic to central China. In contrast, the Dipteronia dyeriana lineage might have been restricted to southwestern China, where it originated, indicating historical confinement. The current restricted distribution of Dipteronia possibly resulted from its evolutionary deceleration in a constantly changing environment.

Phylogeny and evolution of Asparagaceae subfamily Nolinoideae: new insights from plastid phylogenomics.

BACKGROUND AND AIMS: Asparagaceae subfamily Nolinoideae is an economically important plant group, but the deep relationships and evolutionary history of the lineage remain poorly understood. Based on a large data set including 37 newly sequenced samples and publicly available plastomes, this study aims to better resolve the inter-tribal relationships of Nolinoideae, and to rigorously examine the tribe-level monophyly of Convallarieae, Ophiopogoneae and Polygonateae. METHODS: Maximum likelihood (ML) and Bayesian inference (BI) methods were used to infer phylogenetic relationships of Nolinoideae at the genus level and above. The diversification history of Nolinoideae was explored using molecular dating. KEY RESULTS: Both ML and BI analyses identically recovered five clades within Nolinoideae, respectively corresponding to Dracaeneae + Rusceae, Polygonateae + Theropogon, Ophiopogoneae, Nolineae, and Convallarieae excluding Theropogon, and most deep nodes were well supported. As Theropogon was embedded in Polygonateae, the plastome phylogeny failed to resolve Convallarieae and Polygonateae as reciprocally monophyletic. Divergence time estimation showed that the origins of most Nolinoideae genera were dated to the Miocene and Pliocene. The youthfulness of Nolinoideae genera is well represented in the three herbaceous tribes (Convallarieae, Ophiopogoneae and Polygonateae) chiefly distributed in temperate areas of the Northern Hemisphere, as the median stem ages of all 14 genera currently belonging to them were estimated at <12.37 Ma. CONCLUSIONS: This study recovered a robust backbone phylogeny, providing new insights for better understanding the evolution and classification of Nolinoideae. Compared with the deep relationships recovered by a previous study based on transcriptomic data, our data suggest that ancient hybridization or incomplete lineage sorting may have occurred in the early diversification of Nolinoideae. Our findings will provide important reference for further study of the evolutionary complexity of Nolinoideae using nuclear genomic data. The recent origin of these herbaceous genera currently belonging to Convallarieae, Ophiopogoneae and Polygonateae provides new evidence to support the hypothesis that the global expansion of temperate habitats caused by the climate cooling over the past 15 million years may have dramatically driven lineage diversification and speciation in the Northern Hemisphere temperate flora.

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

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.

Parisfargosides A-E, five new cholestane glycosides from the rhizomes of Paris fargesii.

Five new cholestane glycosides, named parisfargosides A-E (1-5), were isolated from the rhizomes of Paris fargesii. Their structures were elucidated on the basis of UV, HR-ESI-MS, 1D and 2D NMR data as well as chemical methods. The structures of all compounds contained α, ß-unsaturated ketone unit. Compounds 3-5 possessed a 16,23-cyclocholest skeleton with 6/6/6/5/5 condensed ring, and the absolute configurations of C-16 and C-23 were confirmed according to ROESY spectra with pyridine­d5 and DMSO­d6 as solvents. In addition, the platelet aggregation activity and cytotoxic activity against five human cancer cell lines (HL-60, A549, SMMC-7721, MDA-MB-231, and SW480) of compounds 1-5 were evaluated.

Parisvaniosides A-E, five new steroidal saponins from Paris vaniotii.

The species of Paris genus is a prolific source of structurally diverse steroidal saponins responsible for multivarious biological properties. The first phytochemical investigation on the steroidal saponin constituents from the rhizomes of Paris vaniotii Lévl. led to the discovery and structural characterization of four new spirostanol saponins, named parisvaniosides A-D (1-4), and one new furostanol glycoside, named parisvanioside E (5), along with eleven known analogues (6-16). Their structures were unambiguously established on the basis of extensive spectroscopic analysis and comparison with the reported spectroscopic data. Compound 1 is a rare spirostanol saponin sharing with a C-9/C-11 double bond and a peroxy group located between C-5 and C-8 of the aglycone, whereas 3 and 4 are unusual C-27 steroidal sapoins with hydroxyl/methoxyl at both C-5 and C-6. Furthermore, 5 is the first furostanol saponin with a unique aglycone featuring two trisubstituted double bonds in ring B. All isolated saponins were evaluated for their anti-inflammatory effects on a lipopolysaccharide (LPS)-stimulated nitric oxide (NO) production model in RAW 264.7 macrophages.

Deciphering the Taxonomic Delimitation of Ottelia acuminata (Hydrocharitaceae) Using Complete Plastomes as Super-Barcodes.

Accurate species delimitation and identification, which is a challenging task in traditional morphology-based taxonomy, is crucial to species conservation. Ottelia acuminata (Hydrocharitaceae) is a severely threatened submerged macrophyte endemic to southwestern China. The taxonomy of O. acuminata, which has long been in dispute, remains unresolved, impeding effective conservation and management practices. Here, we aim to address the long-standing issues concerning species boundary and intraspecific subdivision of O. acuminata using complete plastome sequences as super-barcodes. The taxonomic delimitation of O. acuminata was explored using phylogenetic inference and two independent sequence-based species delimitation schemes: automatic barcode gap discovery (ABGD) and multi-rate Poisson tree processes (mPTP). The reciprocally reinforcing results support the reduction of the closely related congeneric species, O. balansae and O. guanyangensis, as two conspecific varieties of O. acuminata. Within the newly defined O. acuminata, accurate varietal identification can be achieved using plastome super-barcodes. These findings will help inform future decisions regarding conservation, management and restoration of O. acuminata. This case study suggests that the use of plastome super-barcodes can provide a solution for species delimitation and identification in taxonomically difficult plant taxa, thus providing great potential to lessen the challenges of inventorying biodiversity, as well as biologically monitoring and assessing threatened species.

Highly degenerate plastomes in two hemiparasitic dwarf mistletoes: Arceuthobium chinense and A. pini (Viscaceae).

MAIN CONCLUSION: The leafless and endophytic habitat may significantly relax the selection pressure on photosynthesis, and plastid transcription and translation, causing the loss/pseudogenization of several essential plastid-encoding genes in dwarf mistletoes. Dwarf mistletoes (Arceuthobium spp., Viscaceae) are the most destructive plant parasites to numerous conifer species worldwide. In this study, the plastid genomes (plastomes) of Arceuthobium chinense Lecomte and A. pini Hawksworth and Wiens were sequenced and characterized. Although dwarf mistletoes are hemiparasites capable of photosynthesis, their plastomes were highly degenerated, as indicated by the smallest plastome size, the lowest GC content, and relatively very few intact genes among the Santalales hemiparasites. Unexpectedly, several essential housekeeping genes (rpoA, rpoB, rpoC1, and rpoC2) and some core photosynthetic genes (psbZ and petL), as well as the rpl33 gene, that is indispensable for plants under stress conditions, were deleted or pseudogenized in the Arceuthobium plastomes. Our data suggest that the leafless and endophytic habit, which heavily relies on the coniferous hosts for nutrients and carbon requirement, may largely relax the selection pressure on photosynthesis, as well as plastid transcription and translation, thus resulting in the loss/pseudogenization of such essential plastid-encoding genes in dwarf mistletoes. Therefore, the higher level of plastome degradation in Arceuthobium species than other Santalales hemiparasites is likely correlated with the evolution of leafless and endophytic habit. A higher degree of plastome degradation in Arceuthobium. These findings provide new insights into the plastome degeneration associated with parasitism in Santalales and deepen our understanding of the biology of dwarf mistletoes.

Plastome phylogenomics of the East Asian endemic genus Dobinea.

Dobinea is a dioecious genus endemic to East Asia that consists of two extant species: Dobinea delavayi and Dobinea vulgaris. Although the genus is morphologically distinct, its phylogenetic position remains controversial. In this study, we investigated the phylogenetic relationships between Dobinea and related taxa by sequencing the whole plastome DNA sequences for both extant species of Dobinea and comparing them to published plastomes within Sapindales. The complete plastomes of D. vulgaris and D. delavayi were 160,683 and 160, 154 base pairs (bp) in length, including a pair of inverted repeat regions (IRs, 26,889 and 26,759 bp) divided by the large single-copy region (LSC, 87,962 and 87,555 bp) and small single-copy region (SSC, 18,943 and 19,081 bp), and identically encoded 113 unique genes (79 protein-coding genes, 30 tRNAs, and 4 rRNA genes). Plastid phylogenomic analyses showed that Dobinea was a well-supported monophyletic unit and sister to the clade including tribes Anacardieae and Rhoideae, which suggests that Dobinea is a member of Anacardiaceae. In addition, molecular dating inferred D. delavayi and D. vulgaris diverged approximately 10.76 Ma, suggesting the divergence between these two species may have been driven by the intensification of the Asian summer monsoon and the establishment of distinct monsoon regimes in East Asia.

Plastome phylogenomics of Cephalotaxus (Cephalotaxaceae) and allied genera.

BACKGROUND AND AIMS: Cephalotaxus is a paleo-endemic genus in East Asia that consists of about 7-9 conifer species. Despite its great economic and ecological importance, the relationships between Cephalotaxus and related genera, as well as the interspecific relationships within Cephalotaxus, have long been controversial, resulting in contrasting taxonomic proposals in delimitation of Cephalotaxaceae and Taxaceae. Based on plastome data, this study aims to reconstruct a robust phylogeny to infer the systematic placement and the evolutionary history of Cephalotaxus. METHODS: A total of 11 plastomes, representing all species currently recognized in Cephalotaxus and two Torreya species, were sequenced and assembled. Combining these with previously published plastomes, we reconstructed a phylogeny of Cephalotaxaceae and Taxaceae with nearly full taxonomic sampling. Under a phylogenetic framework and molecular dating, the diversification history of Cephalotaxus and allied genera was explored. KEY RESULTS: Phylogenetic analyses of 81 plastid protein-coding genes recovered robust relationships between Cephalotaxus and related genera, as well as providing a well-supported resolution of interspecific relationships within Cephalotaxus, Taxus, Torreya and Amentotaxus. Divergence time estimation indicated that most extant species of these genera are relatively young, although fossil and other molecular evidence consistently show that these genera are ancient plant lineages. CONCLUSIONS: Our results justify the taxonomic proposal that recognizes Cephalotaxaceae as a monotypic family, and contribute to a clear-cut delineation between Cephalotaxaceae and Taxaceae. Given that extant species of Cephalotaxus are derived from recent divergence events associated with the establishment of monsoonal climates in East Asia and Pleistocene climatic fluctuations, they are not evolutionary relics.

The complete plastomes of two flowering epiparasites (Phacellaria glomerata and P. compressa): Gene content, organization, and plastome degradation.

A plant parasite obligately parasitizing another plant parasite is referred to as epiparasite, which is extremely rare in angiosperms, and their complete plastome sequences have not been characterized to date. In this study, the complete plastomes of two flowering epiparasites: Phacellaria compressa and P. glomerata (Amphorogynaceae, Santalales) were sequenced. The plastomes of both species are of similar size, structure, gene content, and arrangement of genes to other hemiparasites in Santalales. Their plastomes were characterized by the functional loss of plastid-encoded NAD(P)H-dehydrogenase and infA genes, which strongly coincides with the general pattern of plastome degradation observed in Santalales hemiparasites. Our study demonstrates that the relatively higher level of nutritional reliance on the host plants and the reduced vegetative bodies of P. compressa and P. glomerata do not appear to cause any unique plastome degradation compared with their closely related hemiparasites.

The Complete Plastomes of Five Hemiparasitic Plants (Osyris wightiana, Pyrularia edulis, Santalum album, Viscum liquidambaricolum, and V. ovalifolium): Comparative and Evolutionary Analyses Within Santalales.

Most species of Santalales (the sandalwood order) are hemiparasites, including both facultative and obligate hemiparasites. Despite its rich diversity, only a small fraction of the species in the sandalwood order have sequenced plastomes. The evolution of parasitism-associated plastome reduction in Santalales remains under-studied. Here, we report the complete plastomes of three facultative hemiparasites (Pyrularia edulis, Cervantesiaceae; Osyris wightiana, and Santalum album, Santalaceae), and two obligate hemiparasites (Viscum liquidambaricolum and Viscum ovalifolium, Viscaceae). Coupled with publicly available data, we investigated the dynamics of plastome degradation in Santalales hemiparasites. Our results indicate that these hemiparasites can be characterized by various degrees of plastome downsizing, structural rearrangement, and gene loss. The loss or pseudogenization of ndh genes was commonly observed in Santalales hemiparasites, which may be correlated to the lifestyle shift from photoautotroph to hemiparasitism. However, the obligate hemiparasites did not exhibit a consistently higher level of gene loss or pseudogenization compared to facultative hemiparasites, which suggests that the degree of plastome reduction is not correlated with the trophic level facultative or obligate hemiparasitism. Instead, closely related taxa tend to possess highly similar plastome size, structure, and gene content. This implies the parasitism-associated plastome degradation in Santalales may evolve in a lineage-specific manner.

Molecular networking uncovers steroidal saponins of Paris tengchongensis.

Based on a method combining the LC-MS/MS molecular networking strategy with the conventional means of phytochemical research, the chemical constituents and the availability of Paris tengchongensis, a new species found in 2017 from Yunnan Province, were investigated for the first time. The molecular networking showed that this species contained the characteristic steroidal glycosides of the genus Paris by comparison of those of Paris polyphylla var. yunnanensis. Furthermore, the detailed investigation on the 80% EtOH extract of its rhizomes resulted to the isolation of twenty steroidal glycosides including three new spirostane-type saponins, named paristengosides A-C (1-3). Their structures were confirmed by spectroscopic analyses (HRMS and NMR) and chemical methods. The new isolates were evaluated for their cytotoxicities against two human cancer cell lines (HEL and MDA-MB-231), anti-inflammatory effects on a lipopolysaccharide (LPS)-stimulated NO production model in RAW264.7 macrophages, anti-AChE, and antimicrobial activities. The results from the molecular networking and the investigation on the chemical constituents suggested that P. tengchongensis can be used as a potential resource of Rhizoma Paridis.

Ultra-Barcoding Discovers a Cryptic Species in Paris yunnanensis (Melanthiaceae), a Medicinally Important Plant.

Ultra-barcoding is a technique using whole plastomes and nuclear ribosomal DNA (nrDNA) sequences for plant species identification. Paris yunnanensis is a medicinal plant of great economic importance for the pharmaceutical industry. However, the alpha taxonomy of P. yunnanensis is still uncertain, hindering effective conservation and management of the germplasm. To resolve long-standing taxonomic disputes regarding this species, we newly generated the complete plastomes and nrDNA sequences from 22 P. yunnanensis accessions. Ultra-barcoding analyses suggest that P. yunnanensis as currently circumscribed is made up of two distinct genetic lineages, corresponding to the two phenotypes ("typical" and "high stem" form) identified early in our study. With distinct morphologies and distribution, the "high stem" form should be recognized as a previously unrecognized species; here it is described as a new species, P. liiana sp. nov. Moreover, the ultra-barcoding data do not support treatment of P. yunnanensis as a conspecific variety under Paris polyphylla. Our study represents a guiding practical application of ultra-barcoding for discovery of cryptic species in taxonomically challenging plant taxa. The findings highlight the great potential of ultra-barcoding as an effective tool for resolving perplexing problems in plant taxonomy.