A curated benchmark dataset for molecular identification based on genome skimming

Genome skimming is a promising sequencing strategy for DNA-based taxonomic identification. However, the lack of standardized datasets for benchmarking genome skimming tools presents a challenge in comparing new methods to existing ones. As part of the development of varKoder, a new tool for DNA-based identification, we curated four datasets designed for comparing molecular identification tools using low-coverage genomes. These datasets comprise vast phylogenetic and taxonomic diversity from closely related species to all taxa currently represented on NCBI SRA. One of them consists of novel sequences from taxonomically verified samples in the plant clade Malpighiales, while the other three datasets compile publicly available data. All include raw genome skim sequences to enable comprehensive testing and validation of a variety molecular species identification methods. We also provide the two-dimensional graphical representations of genomic data used in varKoder. These datasets represent a reliable resource for researchers to assess the accuracy, efficiency, and robustness of new tools to varKoder and other methods in a consistent and reproducible manner.

An Accessible Metagenomic Strategy Allows for Better Characterisation of Invertebrate Bulk Samples

DNA-based techniques are a popular approach for assessing biodiversity in ecological research, especially for organisms which are difficult to detect or identify morphologically. Metabarcoding, the most established method for determining species composition and relative abundance in bulk samples, can be more sensitive and time- and cost-effective than traditional morphological approaches. However, one drawback of this method is PCR bias caused by between-species variation in the amplification efficiency of a marker gene. Metagenomics, bypassing PCR amplification, has been proposed as an alternative to overcome this bias. Several studies have already shown the promising potential of metagenomics, but they all indicate the unavailability of reference genomes for most species in any ecosystem as one of the primary bottlenecks preventing its wider implementation. In this study, we present a strategy that combines unassembled reads of low-coverage whole genome sequencing and publicly available reference genomes to construct a genomic reference database, thus circumventing high sequencing costs and intensive bioinformatic processing. We show that this approach is superior to metabarcoding for approximating relative biomass of macrobenthos species from bulk samples. Furthermore, these results can be obtained with a sequencing effort comparable to metabarcoding. The strategy presented here can thus accelerate the implementation of metagenomics in biodiversity assessments, as it should be relatively easy to adopt by laboratories familiar with metabarcoding and can be used as an accessible alternative.

Exploring Mitonuclear Discordance: Ghost Introgression From an Ancient Extinction Lineage in the Odorrana swinhoana Complex

Mitonuclear discordance, the incongruence between mitochondrial DNA (mtDNA) and nuclear DNA (nuDNA), is a well-documented phenomenon with various potential explanations. One emerging hypothesis, ghost introgression, refers to the genetic contribution of an ancient, extinct or unsampled lineage and can now be tested using modern genomic data and demographic models. In this study, we investigated the evolutionary history of the Odorrana swinhoana complex (Anura: Ranidae), which includes O. swinhoana, O. utsunomiyaorum and an unidentified population with highly divergent mtDNA. While mitochondrial phylogeny suggested this population as a basal lineage, nuclear data from ddRADseq revealed it as a mixture of the most derived O. swinhoana nuclear sequences combined with ancient mtDNA. Demographic modelling further supported ghost introgression, as all models incorporating a ghost population outperformed those without it. These findings suggest that an eastward expansion of western O. swinhoana replaced an ancient Odorrana lineage, leaving only its mtDNA and fragments of its nuclear genome in the hybrid population. Our results provide one of the first documented cases of ghost introgression in amphibians and highlight its potential as a widespread evolutionary process. This study also underscores the risks of relying solely on mtDNA for phylogenetic reconstruction and species delimitation.

Specimen Identification Through Multilocus Species Tree Constructed From Single-Copy Orthologs (SCOs): A Case Study in Cymbidium Subgenus Jensoa

Standard barcodes and ultra-barcode encounter significant challenges when delimiting and discriminating closely related species characterized by deep coalescence, hybrid speciation, gene flow, or low sequence variation. Single-copy orthologs (SCOs) have been widely recognized as standardized nuclear markers in metazoan DNA taxonomy, yet their application in plant taxonomy remains unexplored. This study evaluates the efficacy of SCOs for identifying recently diverged species within the Cymbidium subgenus Jensoa, where ultra-barcodes have previously shown limited resolution. Remarkably, over 90% of the 9094 targeted reference SCOs, inferred from three Cymbidium genomes, were successfully retrieved for all 11 representative species in subg. Jensoa using ALiBaSeq at a minimal 5× depth from whole genome shotgun sequences. The species tree, reconstructed from multiple refined SCO matrices under the coalescent model, effectively distinguished all species and identified mislabeled or misidentified specimens. The comprehensive and refined SCO matrices produced by our pipeline not only enhance phylogenetic analysis but also improve the precision of species diagnosis. Additionally, biparentally inherited SCOs, serving as multi-locus markers, not only augment the effectiveness of DNA barcoding but also support a transition to multi-locus, species-tree-based specimen assignment strategies.

Complete plastomes serve as desirable molecular makers for precise identification of Asparagus cochinchinensis (Asparagaceae) and nine other congeneric species frequently utilized as its adulterants

BACKGROUD

The processed tuberous roots of Asparagus cochinchinensis (Asparagaceae), known as Asparagi Radix, have long been used in East Asia (particularly in China) as traditional medicines and play an indispensable role in the pharmaceutical industry. However, the frequent adulteration of Asparagi Radix with processed tuberous roots obtained from nine other congeneric species could potentially compromise the quality control measures for related pharmaceutical products, while also posing challenges to the conservation and rational exploitation of the nine adulterant congeneric species that are also used as traditional ethnomedicines. Given this issue, this study aims to develop a molecular authentication method for the accurate identification of A. cochinchinensis and the nine congeneric adulterants, employing the genome skimming approach to generate complete plastid genomes (plastomes) and nuclear ribosomal DNA (nrDNA) arrays as the candidate molecular markers.

RESULTS

Through comprehensive phylogenetic and genetic distance analyses based on extensive sampling at both inter- and intra-specific levels, the efficacy of the two candidate molecular markers was assessed by investigating whether their inter-specific genetic divergences align with the taxonomically delineated species boundaries.

CONCLUSION

The results indicated that complete plastomes exhibit superior performance for accurately identifying A. cochinchinensis (the botanical source of Asparagi Radix) and the nine congeneric adulterants, thus can serve as the optimal molecular markers for effective authentication of Asparagi Radix. The desirable discriminative power demonstrated by complete plastomes suggests that the PCR-free molecular authentication method developed in this study will not only contribute to the quality control of pharmaceutical products derived from Asparagi Radix but also facilitate the conservation efforts and rational exploitation of the nine Asparagus species commonly used as adulterants.

An integrative approach clarifies species delimitation and biogeographic history of Debregeasia (Urticaceae)

Integrative data from plastid and nuclear loci are increasingly utilized to resolve species boundaries and phylogenetic relationships within major angiosperm clades. Debregeasia (Urticaceae), an economically important genus, presents challenges in species delimitation due to its overlapping morphological traits and unstable taxonomic assignments. Here, we analyzed 14 morphological traits and generated 12 data matrices from the plastomes and nrDNA using genome skimming from the nine recognized morphospecies to clarify species boundaries and assess barcode performance in Debregeasia. We also used a universal set of 353 nuclear genes to explore reticulate evolution and biogeographic history of Debregeasia. Plastomes of Debregeasia exhibited the typical quadripartite structure with conserved gene content and marginal independent variations in the SC/IR boundary at inter- and intra-specific levels. Three Debregeasia species were non-monophyletic and could not be discerned by any barcode; however, ultra-barcodes identified the remaining six (67%), outperforming standard barcodes (56%). Our phylogenetic analyses placed Debregeasia wallichiana outside the genus and suggested six monophyletic clades in Debregeasia, although the placement between Debregeasia hekouensis and Debregeasia libera varied. There was extensive trait overlap in key morphologically diagnostic characters, with reticulation analysis showing potentially pervasive hybridization, likely influenced by speciation patterns and overlaps between species ranges. We inferred that Debregeasia crown diversification began at ca. 12.82 Ma (95% HPD: 11.54-14.63 Ma) in the mid-Miocene within Australia, followed by vicariance and later long-distance dispersal, mainly out of southern China. Our findings highlight the utility of genomic data with integrative lines of evidence to refine species delimitation and explore evolutionary relationships in complex plant lineages.

Assembly-free reads accurate identification (AFRAID) approach outperforms other methods of DNA barcoding in the walnut family (Juglandaceae)

DNA barcoding has been extensively used for species identification. However, species identification of mixed samples or degraded DNA is limited by current DNA barcoding methods. In this study, we use plant species in Juglandaceae to evaluate an assembly-free reads accurate identification (AFRAID) method of species identification, a novel approach for precise species identification in plants. Specifically, we determined (1) the accuracy of DNA barcoding approaches in delimiting species in Juglandaceae, (2) the minimum size of chloroplast dataset for species discrimination, and (3) minimum amount of next generation sequencing (NGS) data required for species identification. We found that species identification rates were highest when whole chloroplast genomes were used, followed by taxon-specific DNA barcodes, and then universal DNA barcodes. Species identification of 100% was achieved when chloroplast genome sequence coverage reached 20% and the original sequencing data reached 500,000 reads. AFRAID accurately identified species for all samples tested after 500,000 clean reads, with far less computing time than common approaches. These results provide a new approach to accurately identify species, overcoming limitations of traditional DNA barcodes. Our method, which uses next generation sequencing to generate partial chloroplast genomes, reveals that DNA barcode regions are not necessarily fixed, accelerating the process of species identification.

From light into shadow: comparative plastomes in Petrocosmea and implications for low light adaptation

BACKGROUND

Plastids originated from an ancient endosymbiotic event and evolved into the photosynthetic organelles in plant cells. They absorb light energy and carbon dioxide, converting them into chemical energy and oxygen, which are crucial for plant development and adaptation. However, little is known about the plastid genome to light adaptation. Petrocosmea, a member of the Gesneriaceae family, comprises approximately 70 species with diverse light environment, serve as an ideal subject for studying plastomes adapt to light.

RESULTS

In this study, we selected ten representative species of Petrocosmea from diverse light environments, assembled their plastid genomes, and conducted a comparative genomic analysis. We found that the plastid genome of Petrocosmea is highly conserved in both structure and gene content. The phylogenetic relationships reconstructed based on the plastid genes were divided into five clades, which is consistent with the results of previous studies. The vast majority of plastid protein-coding genes were under purifying selection, with only the rps8 and rps16 genes identified under positive selection in different light environments. Notably, significant differences of evolutionary rate were observed in NADH dehydrogenase, ATPase ribosome, and RNA polymerase between Clade A and the other clades. Additionally, we identified ycf1 and several intergenic regions (trnH-psbA, trnK-rps16, rpoB-trnC, petA-psbJ, ccsA-trnL, rps16-trnQ, and trnS-trnG) as candidate barcodes for this emerging ornamental horticulture.

CONCLUSION

We newly assembled ten plastid genomes of Petrocosmea and identified several hypervariable regions, providing genetic resources and candidate markers for this promising emerging ornamental horticulture. Furthermore, our study suggested that rps8 and rps16 were under positive selection and that the evolutionary patterns of NADH dehydrogenase, ATPase ribosome, and RNA polymerase were related to the diversity light environment in Petrocosmea. This revealed an evolutionary scenario for light adaptation of the plastid genome in plants.

Molecular identification and studies on genetic diversity and structure-related GC heterogeneity of Spatholobus Suberectus based on ITS2

To determine the role of internal transcribed spacer 2 (ITS2) in the identification of Spatholobus suberectus and explore the genetic diversity of S. suberectus. A total of 292 ITS2s from S. suberectus and 17 other plant species were analysed. S. suberectus was clustered separately in the phylogenetic tree. The genetic distance between species was greater than that within S. suberectus. Synonymous substitution rate (Ks) analysis revealed that ITS2 diverged the most recently within S. suberectus (Ks = 0.0022). These findings suggested that ITS2 is suitable for the identification of S. suberectus. The ITS2s were divided into 8 haplotypes and 4 evolutionary branches on the basis of secondary structure, indicating that there was variation within S. suberectus. Evolutionary analysis revealed that the GC content of paired regions (pGC) was greater than that of unpaired regions (upGC), and the pGC showed a decreasing trend, whereas the upGC remained unchanged. Single-base mutation was the main cause of base pair substitution. In both the initial state and the equilibrium state, the substitution rate of GC was higher than that of AU. The increase in the GC content was partly attributed to GC-biased gene conversion (gBGC). High GC content reflected the high recombination and mutation rates of ITS2, which is the basis for species identification and genetic diversity. We characterized the sequence and structural characteristics of S. suberectus ITS2 in detail, providing a reference and basis for the identification of S. suberectus and its products, as well as the protection and utilization of wild resources.

Phylogeny, adaptive evolution, and taxonomy of Acronema (Apiaceae): evidence from plastid phylogenomics and morphological data

Introduction

The genus Acronema, belonging to Apiaceae, includes approximately 25 species distributed in the high-altitude Sino-Himalayan region from E Nepal to SW China. This genus is a taxonomically complex genus with often indistinct species boundaries and problematic generic delimitation with Sinocarum and other close genera, largely due to the varied morphological characteristics.

Methods

To explore the phylogenetic relationships and clarify the limits of the genus Acronema and its related genera, we reconstructed a reliable phylogenetic framework with high support and resolution based on two molecular datasets (plastome data and ITS sequences) and performed morphological analyses.

Results

Both phylogenetic analyses robustly supported that Acronema was a non-monophyletic group that fell into two clades: Acronema Clade and East-Asia Clade. We also newly sequenced and assembled sixteen Acronema complete plastomes and performed comprehensively comparative analyses for this genus. The comparative results showed that the plastome structure, gene number, GC content, codon bias patterns were high similarity, but varied in borders of SC/IR and we identified six different types of SC/IR border. The SC/IR boundaries of Acronema chienii were significantly different from the other Acronema members which was consistent with the type VI pattern in the genus Tongoloa. We also identified twelve potential DNA barcode regions (ccsA, matK, ndhF, ndhG, psaI, psbI, rpl32, rps15, ycf1, ycf3, psaI-ycf4 and psbM-trnD) for species identification in Acronema. The molecular evolution of Acronema was relatively conservative that only one gene (petG) was found to be under positive selection (ω = 1.02489).

Discussion

The gene petG is one of the genes involved in the transmission of photosynthetic electron chains during photosynthesis, which plays a crucial role in the process of photosynthesis in plants. This is also a manifestation of the adaptive evolution of plants in high-altitude areas to the environment. In conclusion, our study provides novel insights into the plastome adaptive evolution, phylogeny, and taxonomy of genus Acronema.

Testing plastomes and nuclear ribosomal DNA sequences as the next-generation DNA barcodes for species identification and phylogenetic analysis in Acer

BACKGROUND

Acer is a taxonomically intractable and speciose genus that contains over 150 species. It is challenging to distinguish Acer species only by morphological method due to their abundant variations. Plastome and nuclear ribosomal DNA (nrDNA) sequences are recommended as powerful next-generation DNA barcodes for species discrimination. However, their efficacies were still poorly studied. The current study will evaluate the application of plastome and nrDNA in species identification and perform phylogenetic analyses for Acer.

RESULT

Based on a collection of 83 individuals representing 55 species (c. 55% of Chinese species) from 13 sections, our barcoding analyses demonstrated that plastomes exhibited the highest (90.47%) species discriminatory power among all plastid DNA markers, such as the standard plastid barcodes matK + rbcL + trnH-psbA (61.90%) and ycf1 (76.19%). And the nrDNA (80.95%) revealed higher species resolution than ITS (71.43%). Acer plastomes show abundant interspecific variations, however, species identification failure may be due to the incomplete lineage sorting (ILS) and chloroplast capture resulting from hybridization. We found that the usage of nrDNA contributed to identifying those species that were unidentified by plastomes, implying its capability to some extent to mitigate the impact of hybridization and ILS on species discrimination. However, combining plastome and nrDNA is not recommended given the cytonuclear conflict caused by potential hybridization. Our phylogenetic analysis covering 19 sections (95% sections of Acer) and 128 species (over 80% species of this genus) revealed pervasive inter- and intra-section cytonuclear discordances, hinting that hybridization has played an important role in the evolution of Acer.

CONCLUSION

Plastomes and nrDNA can significantly improve the species resolution in Acer. Our phylogenetic analysis uncovered the scope and depth of cytonuclear conflict in Acer, providing important insights into its evolution.

A comprehensive DNA barcoding reference database for Plecoptera of Switzerland

DNA barcoding is an essential tool in modern biodiversity sciences. Despite considerable work to barcode the tree of life, many groups, including insects, remain partially or totally unreferenced, preventing barcoding from reaching its full potential. Aquatic insects, especially the three orders Ephemeroptera, Plecoptera, and Trichoptera (EPT), are key freshwater quality indicators worldwide. Among them, Plecoptera (stoneflies), which are among the most sensitive aquatic insects to habitat modification, play a central role in river monitoring surveys. Here, we present an update of the Plecoptera reference database for (meta)barcoding in Switzerland, now covering all 118 species known from this country. Fresh specimens, mostly from rare or localized species, were collected, and 151 new CO1 barcodes were generated. These were merged with the 422 previously published sequences, resulting in a dataset of 573 barcoded specimens. Our CO1 dataset was delimited in 115 CO1 clusters based on a priori morphological identifications, of which 17% are newly reported for Switzerland, and 4% are newly reported globally. Among the 115 CO1 clusters, 85% showed complete congruence with morphology. Distance-based analysis indicated local barcoding gaps in 97% of the CO1 clusters. This study significantly improves the Swiss reference database for stoneflies, enhancing future species identification accuracy and biodiversity monitoring. Additionally, this work reveals cryptic diversity and incongruence between morphology and barcodes, both presenting valuable opportunities for future integrative taxonomic studies. Voucher specimens, DNA extractions and reference barcodes are available for future developments, including metabarcoding and environmental DNA surveys.

Transcriptome and comparative chloroplast genome analysis of Taxus yunnanensis individuals with high and low paclitaxel yield

Paclitaxel is a potent anti-cancer drug that is mainly produced through semi-synthesis, which still requires plant materials as precursors. The content of paclitaxel and 10-deacetyl baccatin III (10-DAB) in Taxus yunnanensis has been found to differ from that of other Taxus species, but there is little research on the mechanism underlying the variation in paclitaxel content in T. yunnanensis of different provenances. In this experiment, the contents of taxoids and precursors in twigs between a high paclitaxel-yielding individual (TG) and a low paclitaxel-yielding individual (TD) of T. yunnanensis were compared, and comparative analyses of transcriptomes as well as chloroplast genomes were performed. High-performance liquid chromatography (HPLC) detection showed that 10-DAB and baccatin III contents in TG were 18 and 47 times those in TD, respectively. Transcriptomic analysis results indicated that genes encoding key enzymes in the paclitaxel biosynthesis pathway, such as taxane 10-β-hydroxylase (T10βH), 10-deacetylbaccatin III 10-O-acetyltransferase (DBAT), and debenzoyl paclitaxel N-benzoyl transferase (DBTNBT), exhibited higher expression levels in TG. Additionally, qRT-PCR showed that the relative expression level of T10βH and DBAT in TG were 29 and 13 times those in TD, respectively. In addition, six putative transcription factors were identified that may be involved in paclitaxel biosynthesis from transcriptome data. Comparative analysis of plastid genomes showed that the TD chloroplast contained a duplicate of rps12, leading to a longer plastid genome length in TD relative to TG. Fifteen mutation hotspot regions were identified between the two plastid genomes that can serve as candidate DNA barcodes for identifying high-paclitaxel-yield individuals. This experiment provides insight into the difference in paclitaxel accumulation among different provenances of T. yunnanensis individuals.

Nuclear Magnetic Resonance Fingerprints and Mini DNA Markers for the Authentication of Cinnamon Species Ingredients Used in Food and Natural Health Products

Cinnamomum verum (syn C. zeylanicum) is considered 'true' cinnamon. However, it is reported that less expensive sources of cinnamon from C. cassia (syn C. aromaticum), C. loureiroi, and C. burmannii (toxic coumarin) may be used in the place of C. verum. We lack the quality assurance tools that are required to differentiate C. verum from other cinnamon species when verifying that the correct species is sourced from ingredient suppliers. The current research on cinnamon species authentication using DNA tools is limited to a few species and the use of high-quality DNA extracted from raw leaf materials. The cinnamon bark traded in the supply chain contains much less DNA and poorer-quality DNA than leaves. Our research advances DNA methods to authenticate cinnamon, as we utilized full-length chloroplast genomes via a genome skimming approach for C. burmannii and C. cassia to facilitate the design of optimal mini DNA markers. Furthermore, we developed and validated the use of NMR fingerprints for several commercial cinnamon species, including the quantification of 16 molecules. NMR fingerprints provided additional data that were useful for quality assessment in cinnamon extract powders and product consistency. Both the new mini DNA markers and NMR fingerprints were tested on commercial cinnamon products.

LocoGSE, a sequence-based genome size estimator for plants

Extensive research has focused on exploring the range of genome sizes in eukaryotes, with a particular emphasis on land plants, where significant variability has been observed. Accurate estimation of genome size is essential for various research purposes, but existing sequence-based methods have limitations, particularly for low-coverage datasets. In this study, we introduce LocoGSE, a novel genome size estimator designed specifically for low-coverage datasets generated by genome skimming approaches. LocoGSE relies on mapping the reads on single copy consensus proteins without the need for a reference genome assembly. We calibrated LocoGSE using 430 low-coverage Angiosperm genome skimming datasets and compared its performance against other estimators. Our results demonstrate that LocoGSE accurately predicts monoploid genome size even at very low depth of coverage (<1X) and on highly heterozygous samples. Additionally, LocoGSE provides stable estimates across individuals with varying ploidy levels. LocoGSE fills a gap in sequence-based plant genome size estimation by offering a user-friendly and reliable tool that does not rely on high coverage or reference assemblies. We anticipate that LocoGSE will facilitate plant genome size analysis and contribute to evolutionary and ecological studies in the field. Furthermore, at the cost of an initial calibration, LocoGSE can be used in other lineages.

Species delimitation of tea plants (Camellia sect. Thea) based on super-barcodes

BACKGROUND

The era of high throughput sequencing offers new paths to identifying species boundaries that are complementary to traditional morphology-based delimitations. De novo species delimitation using traditional or DNA super-barcodes serve as efficient approaches to recognizing putative species (molecular operational taxonomic units, MOTUs). Tea plants (Camellia sect. Thea) form a group of morphologically similar species with significant economic value, providing the raw material for tea, which is the most popular nonalcoholic caffeine-containing beverage in the world. Taxonomic challenges have arisen from vague species boundaries in this group.

RESULTS

Based on the most comprehensive sampling of C. sect. Thea by far (165 individuals of 39 morphospecies), we applied three de novo species delimitation methods (ASAP, PTP, and mPTP) using plastome data to provide an independent evaluation of morphology-based species boundaries in tea plants. Comparing MOTU partitions with morphospecies, we particularly tested the congruence of MOTUs resulting from different methods. We recognized 28 consensus MOTUs within C. sect. Thea, while tentatively suggesting that 11 morphospecies be discarded. Ten of the 28 consensus MOTUs were uncovered as morphospecies complexes in need of further study integrating other evidence. Our results also showed a strong imbalance among the analyzed MOTUs in terms of the number of molecular diagnostic characters.

CONCLUSION

This study serves as a solid step forward for recognizing the underlying species boundaries of tea plants, providing a needed evidence-based framework for the utilization and conservation of this economically important plant group.

Retrieving complete plastid genomes of endangered Guibourtia timber using hybridization capture for forensic identification and phylogenetic analysis

The high economic value and increased demand for timber have led to illegal logging and overexploitation, threatening wild populations. In this context, there is an urgent need to develop effective and accurate forensic tools for identifying endangered Guibourtia timber species to protect forest ecosystem resources and regulate their trade. In this study, a hybridization capture method was developed and applied to explore the feasibility of retrieving complete plastid genomes from Guibourtia sapwood and heartwood specimens stored in a xylarium (wood collection). We then carried out forensic identification and phylogenetic analyses of Guibourtia within the subfamily Detarioideae. This study is the first to successfully retrieve high-quality plastid genomes from xylarium specimens, with 76.95-99.97% coverage. The enrichment efficiency, sequence depth, and coverage of plastid genomes from sapwood were 16.73 times, 70.47 times and 1.14 times higher, respectively, than those from heartwood. Although the DNA capture efficiency of heartwood was lower than that of sapwood, the hybridization capture method used in this study is still suitable for heartwood DNA analysis. Based on the complete plastid genome, we identified six endangered or commonly traded Guibourtia woods at the species level. This technique also has the potential for geographical traceability, especially for Guibourtia demeusei and Guibourtia ehie. Meanwhile, Bayesian phylogenetic analysis suggested that these six Guibourtia species diverged from closely related species within the subfamily Detarioideae ca. 18 Ma during the Miocene. The DNA reference database established based on the xylarium specimens provides admissible evidence for diversity conservation and evolutionary analyses of endangered Guibourtia species.

Global exact optimisations for chloroplast structural haplotype scaffolding

BACKGROUND

Scaffolding is an intermediate stage of fragment assembly. It consists in orienting and ordering the contigs obtained by the assembly of the sequencing reads. In the general case, the problem has been largely studied with the use of distances data between the contigs. Here we focus on a dedicated scaffolding for the chloroplast genomes. As these genomes are small, circular and with few specific repeats, numerous approaches have been proposed to assemble them. However, their specificities have not been sufficiently exploited.

RESULTS

We give a new formulation for the scaffolding in the case of chloroplast genomes as a discrete optimisation problem, that we prove the decision version to be [Formula: see text]-Complete. We take advantage of the knowledge of chloroplast genomes and succeed in expressing the relationships between a few specific genomic repeats in mathematical constraints. Our approach is independent of the distances and adopts a genomic regions view, with the priority on scaffolding the repeats first. In this way, we encode the structural haplotype issue in order to retrieve several genome forms that coexist in the same chloroplast cell. To solve exactly the optimisation problem, we develop an integer linear program that we implement in Python3 package khloraascaf. We test it on synthetic data to investigate its performance behaviour and its robustness against several chosen difficulties.

CONCLUSIONS

We succeed to model biological knowledge on genomic structures to scaffold chloroplast genomes. Our results suggest that modelling genomic regions is sufficient for scaffolding repeats and is suitable for finding several solutions corresponding to several genome forms.

Analyses of Nuclear Reads Obtained Using Genome Skimming

In this protocol paper, we review a set of methods developed in recent years for analyzing nuclear reads obtained from genome skimming. As the cost of sequencing drops, genome skimming (low-coverage shotgun sequencing of a sample) becomes increasingly a cost-effective method of measuring biodiversity at high resolution. While most practitioners only use assembled over-represented organelle reads from a genome skim, the vast majority of the reads are nuclear. Using assembly-free and alignment-free methods described in this protocol, we can compare samples to each other and reference genomes to compute distances, characterize underlying genomes, and infer evolutionary relationships.

Comparative Analysis of Luisia (Aeridinae, Orchidaceae) Plastomes Shed Light on Plastomes Evolution and Barcodes Investigation

Luisia, a genus of the subtribe Aeridinae of Orchidaceae, comprises ca. 40 species. Members of Luisia exhibit unique morphological characteristics and represent a valuable ornamental orchid genus. However, due to the scarcity of distinct morphological characters, species identification within this genus is ambiguous and controversial. In the present study, next-generation sequencing (NGS) methods were used to assemble the plastomes of five Luisia species and compare them with one publicly available Luisia plastid genome data. The plastomes of Luisia possessed a quadripartite structure, with sizes ranging from 146,243 bp to 147,430 bp. The plastomes of six Luisia species contained a total of 120 genes, comprising 74 protein-coding genes, 38 tRNA genes and eight rRNA genes. Notably, all ndh genes were pseudogenized or lost. An analysis of codon usage bias showed that leucine (Leu) exhibited the highest frequency, while cysteine (Cys) exhibited the lowest frequency. A total of 57 to 64 SSRs and 42 to 49 long repeats were identified. Five regions and five coding sequences were identified for DNA barcodes, based on the nucleotide diversity (Pi) analysis. The species of Luisia constituted a monophyletic group and were sister to Paraphalaenopsis with strong support. Our study deepens the understanding of species identification, plastome evolution and the phylogenetic positions of Luisia.

A Review on the Morphology, Cultivation, Identification, Phytochemistry, and Pharmacology of Kitagawia praeruptora (Dunn) Pimenov

Kitagawia praeruptora (Dunn) Pimenov, commonly known as Qianhu in China, is a widely used folk Chinese herbal medicine. This article reviews its botanical traits, ethnopharmacology, cultivation techniques, identification, phytochemical compositions, and pharmacological effects. Over 70 coumarin compounds, including simple coumarins, pyranocoumarins, and furanocoumarins, have been isolated within this plant. Additionally, K. praeruptora contains other components such as flavonoids, fatty acids, benzoic acids, and sterols. This information highlights the importance of utilizing active ingredients and excavating pharmacological effects. With its remarkable versatility, K. praeruptora exhibits a wide range of pharmacological effects. It has been found to possess expectorant and bronchodilator properties, cardiovascular protection, antimicrobial and antioxidant activities, anti-tumor effects, and even antidiabetic properties. It is recommended to focus on the development of new drugs that leverage the active ingredients of K. praeruptora and explore its potential for new clinical applications and holistic utilization.

Comparative and phylogenetic analysis of Chiloschista (Orchidaceae) species and DNA barcoding investigation based on plastid genomes

BACKGROUND

Chiloschista (Orchidaceae, Aeridinae) is an epiphytic leafless orchid that is mainly distributed in tropical or subtropical forest canopies. This rare and threatened orchid lacks molecular resources for phylogenetic and barcoding analysis. Therefore, we sequenced and assembled seven complete plastomes of Chiloschista to analyse the plastome characteristics and phylogenetic relationships and conduct a barcoding investigation.

RESULTS

We are the first to publish seven Chiloschista plastomes, which possessed the typical quadripartite structure and ranged from 143,233 bp to 145,463 bp in size. The plastomes all contained 120 genes, consisting of 74 protein-coding genes, 38 tRNA genes and eight rRNA genes. The ndh genes were pseudogenes or lost in the genus, and the genes petG and psbF were under positive selection. The seven Chiloschista plastomes displayed stable plastome structures with no large inversions or rearrangements. A total of 14 small inversions (SIs) were identified in the seven Chiloschista plastomes but were all similar within the genus. Six noncoding mutational hotspots (trnNGUU-rpl32 > rpoB-trnCGCA > psbK-psbI > psaC-rps15 > trnEUUC-trnTGGU > accD-psaI) and five coding sequences (ycf1 > rps15 > matK > psbK > ccsA) were selected as potential barcodes based on nucleotide diversity and species discrimination analysis, which suggested that the potential barcode ycf1 was most suitable for species discrimination. A total of 47-56 SSRs and 11-14 long repeats (> 20 bp) were identified in Chiloschista plastomes, and they were mostly located in the large single copy intergenic region. Phylogenetic analysis indicated that Chiloschista was monophyletic. It was clustered with Phalaenopsis and formed the basic clade of the subtribe Aeridinae with a moderate support value. The results also showed that seven Chiloschista species were divided into three major clades with full support.

CONCLUSION

This study was the first to analyse the plastome characteristics of the genus Chiloschista in Orchidaceae, and the results showed that Chiloschista plastomes have conserved plastome structures. Based on the plastome hotspots of nucleotide diversity, several genes and noncoding regions are suitable for phylogenetic and population studies. Chiloschista may provide an ideal system to investigate the dynamics of plastome evolution and DNA barcoding investigation for orchid studies.

Towards exhaustive community ecology via DNA metabarcoding

Exhaustive biodiversity data, covering all the taxa in an environment, would be fundamental to understand how global changes influence organisms living at different trophic levels, and to evaluate impacts on interspecific interactions. Molecular approaches such as DNA metabarcoding are boosting our ability to perform biodiversity inventories. Nevertheless, even though a few studies have recently attempted exhaustive reconstructions of communities, holistic assessments remain rare. The majority of metabarcoding studies published in the last years used just one or two markers and analysed a limited number of taxonomic groups. Here, we provide an overview of emerging approaches that can allow all-taxa biological inventories. Exhaustive biodiversity assessments can be attempted by combining a large number of specific primers, by exploiting the power of universal primers, or by combining specific and universal primers to obtain good information on key taxa while limiting the overlooked biodiversity. Multiplexes of primers, shotgun sequencing and capture enrichment may provide a better coverage of biodiversity compared to standard metabarcoding, but still require major methodological advances. Here, we identify the strengths and limitations of different approaches, and suggest new development lines that might improve broad scale biodiversity analyses in the near future. More holistic reconstructions of ecological communities can greatly increase the value of metabarcoding studies, improving understanding of the consequences of ongoing environmental changes on the multiple components of biodiversity.

A DNA barcode library for woody plants in tropical and subtropical China

The application of DNA barcoding has been significantly limited by the scarcity of reliable specimens and inadequate coverage and replication across all species. The deficiency of DNA barcode reference coverage is particularly striking for highly biodiverse subtropical and tropical regions. In this study, we present a comprehensive barcode library for woody plants in tropical and subtropical China. Our dataset includes a standard barcode library comprising the four most widely used barcodes (rbcL, matK, ITS, and ITS2) for 2,520 species from 4,654 samples across 49 orders, 144 families, and 693 genera, along with 79 samples identified at the genus level. This dataset also provides a super-barcode library consisting of 1,239 samples from 1,139 species, 411 genera, 113 families, and 40 orders. This newly developed library will serve as a valuable resource for DNA barcoding research in tropical and subtropical China and bordering countries, enable more accurate species identification, and contribute to the conservation and management of tropical and subtropical forests.

The role of seasonality in shaping the interactions of honeybees with other taxa

The Eltonian niche of a species is defined as its set of interactions with other taxa. How this set varies with biotic, abiotic and human influences is a core question of modern ecology. In seasonal environments, the realized Eltonian niche is likely to vary due to periodic changes in the occurrence and abundance of interaction partners and changes in species behavior and preferences. Also, human management decisions may leave strong imprints on species interactions. To compare the impact of seasonality to that of management effects, honeybees provide an excellent model system. Based on DNA traces of interaction partners archived in honey, we can infer honeybee interactions with floral resources and microbes in the surrounding habitats, their hives, and themselves. Here, we resolved seasonal and management-based impacts on honeybee interactions by sampling beehives repeatedly during the honey-storing period of honeybees in Finland. We then use a genome-skimming approach to identify the taxonomic contents of the DNA in the samples. To compare the effects of the season to the effects of location, management, and the colony itself in shaping honeybee interactions, we used joint species distribution modeling. We found that honeybee interactions with other taxa varied greatly among taxonomic and functional groups. Against a backdrop of wide variation in the interactions documented in the DNA content of honey from bees from different hives, regions, and beekeepers, the imprint of the season remained relatively small. Overall, a honey-based approach offers unique insights into seasonal variation in the identity and abundance of interaction partners among honeybees. During the summer, the availability and use of different interaction partners changed substantially, but hive- and taxon-specific patterns were largely idiosyncratic as modified by hive management. Thus, the beekeeper and colony identity are as important determinants of the honeybee's realized Eltonian niche as is seasonality.

Chloroplast genome, nuclear ITS regions, mitogenome regions, and Skmer analysis resolved the genetic relationship among Cinnamomum species in Sri Lanka

Cinnamomum species have gained worldwide attention because of their economic benefits. Among them, C. verum (synonymous with C. zeylanicum Blume), commonly known as Ceylon Cinnamon or True Cinnamon is mainly produced in Sri Lanka. In addition, Sri Lanka is home to seven endemic wild cinnamon species, C. capparu-coronde, C. citriodorum, C. dubium, C. litseifolium, C. ovalifolium, C. rivulorum and C. sinharajaense. Proper identification and genetic characterization are fundamental for the conservation and commercialization of these species. While some species can be identified based on distinct morphological or chemical traits, others cannot be identified easily morphologically or chemically. The DNA barcoding using rbcL, matK, and trnH-psbA regions could not also resolve the identification of Cinnamomum species in Sri Lanka. Therefore, we generated Illumina Hiseq data of about 20x coverage for each identified species and a C. verum sample (India) and assembled the chloroplast genome, nuclear ITS regions, and several mitochondrial genes, and conducted Skmer analysis. Chloroplast genomes of all eight species were assembled using a seed-based method.According to the Bayesian phylogenomic tree constructed with the complete chloroplast genomes, the C. verum (Sri Lanka) is sister to previously sequenced C. verum (NC_035236.1, KY635878.1), C. dubium and C. rivulorum. The C. verum sample from India is sister to C. litseifolium and C. ovalifolium. According to the ITS regions studied, C. verum (Sri Lanka) is sister to C. verum (NC_035236.1), C. dubium and C. rivulorum. Cinnamomum verum (India) shares an identical ITS region with C. ovalifolium, C. litseifolium, C. citriodorum, and C. capparu-coronde. According to the Skmer analysis C. verum (Sri Lanka) is sister to C. dubium and C. rivulorum, whereas C. verum (India) is sister to C. ovalifolium, and C. litseifolium. The chloroplast gene ycf1 was identified as a chloroplast barcode for the identification of Cinnamomum species. We identified an 18 bp indel region in the ycf1 gene, that could differentiate C. verum (India) and C. verum (Sri Lanka) samples tested.

Cryptic Species in Ecotoxicology

The advent of genetic methods has led to the discovery of an increasing number of species that previously could not be distinguished from each other on the basis of morphological characteristics. Even though there has been an exponential growth of publications on cryptic species, such species are rarely considered in ecotoxicology. Thus, the particular question of ecological differentiation and the sensitivity of closely related cryptic species is rarely addressed. Tackling this question, however, is of key importance for evolutionary ecology, conservation biology, and, in particular, regulatory ecotoxicology. At the same time, the use of species with (known or unknown) cryptic diversity might be a reason for the lack of reproducibility of ecotoxicological experiments and implies a false extrapolation of the findings. Our critical review includes a database and literature search through which we investigated how many of the species most frequently used in ecotoxicological assessments show evidence of cryptic diversity. We found a high proportion of reports indicating overlooked species diversity, especially in invertebrates. In terrestrial and aquatic realms, at least 67% and 54% of commonly used species, respectively, were identified as cryptic species complexes. The issue is less prominent in vertebrates, in which we found evidence for cryptic species complexes in 27% of aquatic and 6.7% of terrestrial vertebrates. We further exemplified why different evolutionary histories may significantly determine cryptic species' ecology and sensitivity to pollutants. This in turn may have a major impact on the results of ecotoxicological tests and, consequently, the outcome of environmental risk assessments. Finally, we provide a brief guideline on how to deal practically with cryptic diversity in ecotoxicological studies in general and its implementation in risk assessment procedures in particular. Environ Toxicol Chem 2023;42:1889-1914. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Target capture and genome skimming for plant diversity studies

Recent technological advances in long-read high-throughput sequencing and assembly methods have facilitated the generation of annotated chromosome-scale whole-genome sequence data for evolutionary studies; however, generating such data can still be difficult for many plant species. For example, obtaining high-molecular-weight DNA is typically impossible for samples in historical herbarium collections, which often have degraded DNA. The need to fast-freeze newly collected living samples to conserve high-quality DNA can be complicated when plants are only found in remote areas. Therefore, short-read reduced-genome representations, such as target capture and genome skimming, remain important for evolutionary studies. Here, we review the pros and cons of each technique for non-model plant taxa. We provide guidance related to logistics, budget, the genomic resources previously available for the target clade, and the nature of the study. Furthermore, we assess the available bioinformatic analyses, detailing best practices and pitfalls, and suggest pathways to combine newly generated data with legacy data. Finally, we explore the possible downstream analyses allowed by the type of data generated using each technique. We provide a practical guide to help researchers make the best-informed choice regarding reduced genome representation for evolutionary studies of non-model plants in cases where whole-genome sequencing remains impractical.

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.

Molecular Barcoding: A Tool to Guarantee Correct Seafood Labelling and Quality and Preserve the Conservation of Endangered Species

The recent increase in international fish trade leads to the need for improving the traceability of fishery products. In relation to this, consistent monitoring of the production chain focusing on technological developments, handling, processing and distribution via global networks is necessary. Molecular barcoding has therefore been suggested as the gold standard in seafood species traceability and labelling. This review describes the DNA barcoding methodology for preventing food fraud and adulteration in fish. In particular, attention has been focused on the application of molecular techniques to determine the identity and authenticity of fish products, to discriminate the presence of different species in processed seafood and to characterize raw materials undergoing food industry processes. In this regard, we herein present a large number of studies performed in different countries, showing the most reliable DNA barcodes for species identification based on both mitochondrial (COI, cytb, 16S rDNA and 12S rDNA) and nuclear genes. Results are discussed considering the advantages and disadvantages of the different techniques in relation to different scientific issues. Special regard has been dedicated to a dual approach referring to both the consumer's health and the conservation of threatened species, with a special focus on the feasibility of the different genetic and genomic approaches in relation to both scientific objectives and permissible costs to obtain reliable traceability.

Hidden species diversity and mito-nuclear discordance within the Mediterranean cone snail, Lautoconus ventricosus

The Mediterranean cone snail, Lautoconus ventricosus, is currently considered a single species inhabiting the whole Mediterranean basin and the adjacent Atlantic coasts. Yet, no population genetic study has assessed its taxonomic status. Here, we collected 245 individuals from 75 localities throughout the Mediterranean Sea and used cox1 barcodes, complete mitochondrial genomes, and genome skims to test whether L. ventricosus represents a complex of cryptic species. The maximum likelihood phylogeny based on complete mitochondrial genomes recovered six main clades (hereby named blue, brown, green, orange, red, and violet) with sufficient sequence divergence to be considered putative species. On the other hand, phylogenomic analyses based on 437 nuclear genes only recovered four out of the six clades: blue and orange clades were thoroughly mixed and the brown one was not recovered. This mito-nuclear discordance revealed instances of incomplete lineage sorting and introgression, and may have caused important differences in the dating of main cladogenetic events. Species delimitation tests proposed the existence of at least three species: green, violet, and red+blue+orange (i.e., cyan). Green plus cyan (with sympatric distributions) and violet, had West and East Mediterranean distributions, respectively, mostly separated by the Siculo-Tunisian biogeographical barrier. Morphometric analyses of the shell using species hypotheses as factor and shell length as covariate showed that the discrimination power of the studied parameters was only 70.2%, reinforcing the cryptic nature of the uncovered species, and the importance of integrative taxonomic approaches considering morphology, ecology, biogeography, and mitochondrial and nuclear population genetic variation.

Taxonomy and Biological Control: New Challenges in an Old Relationship

Biological control and taxonomy are continuously developing fields with remarkable impacts on society. At least 80 years of literature have documented this relationship, which remains essentially the same in its mutualistic nature, as well as in its major challenges. From the perspective of Brazilian taxonomists, we discuss the impacts of important scientific and social developments that directly affect research in these areas, posing new challenges for this lasting relationship. The increasing restrictions and concerns regarding the international transit of organisms require improvements in research related to risk assessment for exotic biological control agents and also stimulate prospecting within the native biota. In our view, this is a positive situation that can foster a closer relationship between taxonomists and applied entomologists, as well as local surveys and taxonomic studies that are necessary before new programs and agents can be implemented. We discuss the essential role of molecular biology in this context, as an iconic example of the synergy between applied sciences and natural history. As our society comes to need safer and more sustainable solutions for food security and the biodiversity crisis, scientific progress will build upon this integration, where biological control and taxonomy play an essential role.

Chloroplast genome characterization of Uncaria guianensis and Uncaria tomentosa and evolutive dynamics of the Cinchonoideae subfamily

Uncaria species are used in traditional medicine and are considered of high therapeutic value and economic importance. This work describes the assembly and annotation of the chloroplast genomes of U. guianensis and U. tomentosa, as well as a comparative analysis. The genomes were sequenced on MiSeq Illumina, assembled with NovoPlasty, and annotated using CHLOROBOX GeSeq. Addictionaly, comparative analysis were performed with six species from NCBI databases and primers were designed in Primer3 for hypervariable regions based on the consensus sequence of 16 species of the Rubiaceae family and validated on an in-silico PCR in OpenPrimeR. The genome size of U. guianensis and U. tomentosa was 155,505 bp and 156,390 bp, respectively. Both Species have 131 genes and GC content of 37.50%. The regions rpl32-ccsA, ycf1, and ndhF-ccsA showed the three highest values of nucleotide diversity within the species of the Rubiaceae family and within the Uncaria genus, these regions were trnH-psbA, psbM-trnY, and rps16-psbK. Our results indicates that the primer of the region ndhA had an amplification success for all species tested and can be promising for usage in the Rubiaceae family. The phylogenetic analysis recovered a congruent topology to APG IV. The gene content and the chloroplast genome structure of the analyzed species are conserved and most of the genes are under negative selection. We provide the cpDNA of Neotropical Uncaria species, an important genomic resource for evolutionary studies of the group.

Horizon scan of DNA-based methods for quality control and monitoring of herbal preparations

Herbal medicines and preparations are widely used in healthcare systems globally, but concerns remain about their quality and safety. New herbal products are constantly being introduced to the market under varying regulatory frameworks, with no global consensus on their definition or characterization. These biologically active mixtures are sold through complex globalized value chains, which create concerns around contamination and profit-driven adulteration. Industry, academia, and regulatory bodies must collaborate to develop innovative strategies for the identification and authentication of botanicals and their preparations to ensure quality control. High-throughput sequencing (HTS) has significantly improved our understanding of the total species diversity within DNA mixtures. The standard concept of DNA barcoding has evolved over the last two decades to encompass genomic data more broadly. Recent research in DNA metabarcoding has focused on developing methods for quantifying herbal product ingredients, yielding meaningful results in a regulatory framework. Techniques, such as loop-mediated isothermal amplification (LAMP), DNA barcode-based Recombinase Polymerase Amplification (BAR-RPA), DNA barcoding coupled with High-Resolution Melting (Bar-HRM), and microfluidics-based methods, offer more affordable tests for the detection of target species. While target capture sequencing and genome skimming are considerably increasing the species identification resolution in challenging plant clades, ddPCR enables the quantification of DNA in samples and could be used to detect intended and unwanted ingredients in herbal medicines. Here, we explore the latest advances in emerging DNA-based technologies and the opportunities they provide as taxa detection tools for evaluating the safety and quality of dietary supplements and herbal medicines.

Multi-marker DNA metabarcoding reveals spatial and sexual variation in the diet of a scarce woodland bird

Avian diet can be affected by site-specific variables, such as habitat, as well as intrinsic factors such as sex. This can lead to dietary niche separation, which reduces competition between individuals, as well as impacting how well avian species can adapt to environmental variation. Estimating dietary niche separation is challenging, due largely to difficulties in accurately identifying food taxa consumed. Consequently, there is limited knowledge of the diets of woodland bird species, many of which are undergoing serious population declines. Here, we show the effectiveness of multi-marker fecal metabarcoding to provide in-depth dietary analysis of a declining passerine in the UK, the Hawfinch (Coccothraustes coccothraustes). We collected fecal samples from (n = 262) UK Hawfinches prior to, and during, the breeding seasons in 2016-2019. We detected 49 and 90 plant and invertebrate taxa, respectively. We found Hawfinch diet varied spatially, as well as between sexes, indicating broad dietary plasticity and the ability of Hawfinches to utilize multiple resources within their foraging environments.

DNA barcoding in herbal medicine: Retrospective and prospective

DNA barcoding has been widely used for herb identification in recent decades, enabling safety and innovation in the field of herbal medicine. In this article, we summarize recent progress in DNA barcoding for herbal medicine to provide ideas for the further development and application of this technology. Most importantly, the standard DNA barcode has been extended in two ways. First, while conventional DNA barcodes have been widely promoted for their versatility in the identification of fresh or well-preserved samples, super-barcodes based on plastid genomes have rapidly developed and have shown advantages in species identification at low taxonomic levels. Second, mini-barcodes are attractive because they perform better in cases of degraded DNA from herbal materials. In addition, some molecular techniques, such as high-throughput sequencing and isothermal amplification, are combined with DNA barcodes for species identification, which has expanded the applications of herb identification based on DNA barcoding and brought about the post-DNA-barcoding era. Furthermore, standard and high-species coverage DNA barcode reference libraries have been constructed to provide reference sequences for species identification, which increases the accuracy and credibility of species discrimination based on DNA barcodes. In summary, DNA barcoding should play a key role in the quality control of traditional herbal medicine and in the international herb trade.

Diversity of dung beetle-associated yeasts from pristine environments of Botswana

Yeast-insect interactions are increasingly becoming an attractive source of discovery for previously unknown, unique, diverse, and industrially relevant yeast species. Despite a wealth of studies that have recently focused on yeasts in symbiotic association with Hymenopteran insects, yeasts associated with Coleopteran insects, such as lignocellulosic-rich dung-dependent beetles, remain poorly studied. Trends in yeast discovery suggest that species richness and diversity can be attributed to the ecological niche of the insect. Here, we considered the potential of dung beetles inhabiting the extreme environments of Botswana, characterized by desert-like conditions (semi-arid to arid and hot) as well as protected pristine environments, as possible attribute niches that can shape the extremophilic and diverse life history strategies of yeasts. We obtained a total of 97 phylogenetically diverse yeast isolates from six species of dung beetles from Botswana's unexplored environments, representing 19 species belonging to 11 genera. The findings suggest that the guts of dung beetles are a rich niche for non-Saccharomyces yeast species. Meyerozyma and Pichia were the most dominant genera associated with dung beetles, representing 55% (53 out of 97) of the yeast isolates in our study. Trichosporon and Cutaneotrichosporon genera represented 32% (31 out of 97) of the isolates. The remaining isolates belonged to Apiotrichum, Candida, Diutina, Naganishia, Rhodotorula, and Wickerhamiella genera (12 out of 97). We found out that about 62% (60 out of 97) of the isolates were potentially new species because of their low internal transcribed spacer (ITS) sequence similarity when compared to the most recent optimal species delineation threshold. A single isolate was unidentifiable using the ITS sequences. Using an in silico polymerase chain reaction-restriction fragment length polymorphism approach, we revealed that there was genetic diversity within isolates of the same species. Our results contribute to the knowledge and understanding of the diversity of dung beetle-associated yeasts.

Enhancing DNA barcode reference libraries by harvesting terrestrial arthropods at the Smithsonian's National Museum of Natural History

The use of DNA barcoding has revolutionised biodiversity science, but its application depends on the existence of comprehensive and reliable reference libraries. For many poorly known taxa, such reference sequences are missing even at higher-level taxonomic scales. We harvested the collections of the Smithsonian's National Museum of Natural History (USNM) to generate DNA barcoding sequences for genera of terrestrial arthropods previously not recorded in one or more major public sequence databases. Our workflow used a mix of Sanger and Next-Generation Sequencing (NGS) approaches to maximise sequence recovery while ensuring affordable cost. In total, COI sequences were obtained for 5,686 specimens belonging to 3,737 determined species in 3,886 genera and 205 families distributed in 137 countries. Success rates varied widely according to collection data and focal taxon. NGS helped recover sequences of specimens that failed a previous run of Sanger sequencing. Success rates and the optimal balance between Sanger and NGS are the most important drivers to maximise output and minimise cost in future projects. The corresponding sequence and taxonomic data can be accessed through the Barcode of Life Data System, GenBank, the Global Biodiversity Information Facility, the Global Genome Biodiversity Network Data Portal and the NMNH data portal.

Highly divergent regions in the complete plastome sequences of Aquilaria are suitable for DNA barcoding applications including identifying species origin of agarwood products

Members of Aquilaria Lam. (Thymelaeaceae) are evergreen trees that are widely distributed in the Indomalesia region. Aquilaria is highly prized for its unique scented resin, agarwood, which is often the subject of unlawful trade activities. Survival of the tree is heavily threatened by destructive harvesting and agarwood poaching, leading to its protection under the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES). Unfortunately, an efficient species identification method, which is crucial to aid in the conservation efforts of Aquilaria is lacking. Here, we described our search for a suitable specific DNA barcode for Aquilaria species using eight complete plastome sequences. We identified five highly variable regions (HVR) (matK-rps16, ndhF-rpl32, psbJ-petA, trnD, and trnT-trnL) in the plastomes. These regions were further analyzed using the neighbor-joining (NJ) method to assess their ability at discriminating the eight species. Coupled with in silico primer design, two potential barcoding regions, psbJ-petA and trnT-trnL, were identified. Their strengths in species delimitation were evaluated individually and in combination, via DNA barcoding analysis. Our findings showed that the combined dataset, psbJ-petA + trnT-trnL, effectively resolved members of the genus Aquilaria by clustering all species into their respective clades. In addition, we demonstrated that the newly proposed DNA barcode was capable at identifying the species of origin of six commercial agarwood samples that were included as unknown samples. Such achievement offers a new technical advancement, useful in the combat against illicit agarwood trades and in assisting the conservation of these valuable species in natural populations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03479-1.

Testing complete plastomes and nuclear ribosomal DNA sequences for species identification in a taxonomically difficult bamboo genus Fargesia

Fargesia, the largest genus within the temperate bamboo tribe Arundinarieae, has more than 90 species mainly distributed in the mountains of Southwest China. The Fargesia bamboos are important components of the subalpine forest ecosystems that provide food and habitat for many endangered animals, including the giant panda. However, species-level identification of Fargesia is difficult. Moreover, the rapid radiation and slow molecular evolutionary rate of Fargesia pose a significant challenge to using DNA barcoding with standard plant barcodes (rbcL, matK, and ITS) in bamboos. With progress in the sequencing technologies, complete plastid genomes (plastomes) and nuclear ribosomal DNA (nrDNA) sequences have been proposed as organelle barcodes for species identification; however, these have not been tested in bamboos. We collected 196 individuals representing 62 species of Fargesia to comprehensively evaluate the discriminatory power of plastomes and nrDNA sequences compared to standard barcodes. Our analysis indicates that complete plastomes have substantially higher discriminatory power (28.6%) than standard barcodes (5.7%), whereas nrDNA sequences show a moderate improvement (65.4%) compared to ITS (47.2%). We also found that nuclear markers performed better than plastid markers, and ITS alone had higher discriminatory power than complete plastomes. The study also demonstrated that plastomes and nrDNA sequences can contribute to intrageneric phylogenetic resolution in Fargesia. However, neither of these sequences were able to discriminate all the sampled species, and therefore, more nuclear markers need to be identified.

Poaceae Chloroplast Genome Sequencing: Great Leap Forward in Recent Ten Years

The first complete chloroplast genome of rice (Oryza sativa) was published in 1989, ushering in a new era of studies of chloroplast genomics in Poaceae. Progresses in Next-Generation Sequencing (NGS) and Third-Generation Sequencing (TGS) technologiesand in the development of genome assembly software, have significantly advanced chloroplast genomics research. Poaceae is one of the most targeted families in chloroplast genome research because of its agricultural, ecological, and economic importance. Over the last 30 years, 2,050 complete chloroplast genome sequences from 40 tribes and 282 genera have been generated, most (97%) of them in the recent ten years. The wealth of data provides the groundwork for studies on species evolution, phylogeny, genetic transformation, and other aspects of Poaceae chloroplast genomes. As a result, we have gained a deeper understanding of the properties of Poaceae chloroplast genomes. Here, we summarize the achievements of the studies of the Poaceae chloroplast genomes and envision the challenges for moving the area ahead.

DNA barcoding of Cymbidium by genome skimming: Call for next-generation nuclear barcodes

Cymbidium is an orchid genus that has undergone rapid radiation and has high ornamental, economic, ecological and cultural importance, but its classification based on morphology is controversial. The plastid genome (plastome), as an extension of plant standard DNA barcodes, has been widely used as a potential molecular marker for identifying recently diverged species or complicated plant groups. In this study, we newly generated 237 plastomes of 50 species (at least two individuals per species) by genome skimming, covering 71.4% of members of the genus Cymbidium. Sequence-based analyses (barcoding gaps and automatic barcode gap discovery) and tree-based analyses (maximum likelihood, Bayesian inference and multirate Poisson tree processes model) were conducted for species identification of Cymbidium. Our work provides a comprehensive DNA barcode reference library for Cymbidium species identification. The results show that compared with standard DNA barcodes (rbcL + matK) as well as the plastid trnH-psbA, the species identification rate of the plastome increased moderately from 58% to 68%. At the same time, we propose an optimized identification strategy for Cymbidium species. The plastome cannot completely resolve the species identification of Cymbidium, the main reasons being incomplete lineage sorting, artificial cultivation, natural hybridization and chloroplast capture. To further explore the potential use of nuclear data in identifying species, the Skmer method was adopted and the identification rate increased to 72%. It appears that nuclear genome data have a vital role in species identification and are expected to be used as next-generation nuclear barcodes.

Be positive: customized reference databases and new, local barcodes balance false taxonomic assignments in metabarcoding studies

Background

In metabarcoding analyses, the taxonomic assignment is crucial to place sequencing data in biological and ecological contexts. This fundamental step depends on a reference database, which should have a good taxonomic coverage to avoid unassigned sequences. However, this goal is rarely achieved in many geographic regions and for several taxonomic groups. On the other hand, more is not necessarily better, as sequences in reference databases belonging to taxonomic groups out of the studied region/environment context might lead to false assignments.

Methods

We investigated the effect of using several subsets of a cytochrome c oxidase subunit I (COI) reference database on taxonomic assignment. Published metabarcoding sequences from the Mediterranean Sea were assigned to taxa using COInr, which is a comprehensive, non-redundant and recent database of COI sequences obtained both from BOLD and NCBI, and two of its subsets: (i) all sequences except insects (COInr-WO-Insecta), which represent the overwhelming majority of COInr database, but are irrelevant for marine samples, and (ii) all sequences from taxonomic families present in the Mediterranean Sea (COInr-Med). Four different algorithms for taxonomic assignment were employed in parallel to evaluate differences in their output and data consistency.

Results

The reduction of the database to more specific custom subsets increased the number of unassigned sequences. Nevertheless, since most of them were incorrectly assigned by the less specific databases, this is a positive outcome. Moreover, the taxonomic resolution (the lowest taxonomic level to which a sequence is attributed) of several sequences tended to increase when using customized databases. These findings clearly indicated the need for customized databases adapted to each study. However, the very high proportion of unassigned sequences points to the need to enrich the local database with new barcodes specifically obtained from the studied region and/or taxonomic group. Including novel local barcodes to the COI database proved to be very profitable: by adding only 116 new barcodes sequenced in our laboratory, thus increasing the reference database by only 0.04%, we were able to improve the resolution for ca. 0.6-1% of the Amplicon Sequence Variants (ASVs).

The plastid genome of twenty-two species from Ferula, Talassia, and Soranthus: comparative analysis, phylogenetic implications, and adaptive evolution

BACKGROUND

The Ferula genus encompasses 180-185 species and is one of the largest genera in Apiaceae, with many of Ferula species possessing important medical value. The previous studies provided more information for Ferula, but its infrageneric relationships are still confusing. In addition, its genetic basis of its adaptive evolution remains poorly understood. Plastid genomes with more variable sites have the potential to reconstruct robust phylogeny in plants and investigate the adaptive evolution of plants. Although chloroplast genomes have been reported within the Ferula genus, few studies have been conducted using chloroplast genomes, especially for endemic species in China.

RESULTS

Comprehensively comparative analyses of 22 newly sequenced and assembled plastomes indicated that these plastomes had highly conserved genome structure, gene number, codon usage, and repeats type and distribution, but varied in plastomes size, GC content, and the SC/IR boundaries. Thirteen mutation hotspot regions were detected and they would serve as the promising DNA barcodes candidates for species identification in Ferula and related genera. Phylogenomic analyses with high supports and resolutions showed that Talassia transiliensis and Soranthus meyeri were nested in the Ferula genus, and thus they should be transferred into the Ferula genus. Our phylogenies also indicated the monophyly of subgenera Sinoferula and subgenera Narthex in Ferula genus. Twelve genes with significant posterior probabilities for codon sites were identified in the positively selective analysis, and their function may relate to the photosystem II, ATP subunit, and NADH dehydrogenase. Most of them might play an important role to help Ferula species adapt to high-temperatures, strong-light, and drought habitats.

CONCLUSION

Plastome data is powerful and efficient to improve the support and resolution of the complicated Ferula phylogeny. Twelve genes with significant posterior probabilities for codon sites were helpful for Ferula to adapt to the harsh environment. Overall, our study supplies a new perspective for comprehending the phylogeny and evolution of Ferula.

A comparative study of the chloroplast genomes of five Lepidium species with high medicinal value

Plantgenomics is a rapidly developing field in medicinal plant research. This study analysed the relevant information of chloroplasts genome sequences of five medicinal plants from the genus Lepidium . We sequenced the complete chloroplast (cp) genomes of Lepidium apetalum Willd. and Lepidium perfoliatum Linnaeus., and assessed their genetic profiles against the reported profiles of Lepidium sativum Linnaeus., Lepidium meyenii Walp., and Lepidium virginicum Linn. We found that L. apetalum and L. perfoliatum possessed 130 distinct genes that included 85 protein-coding, 37 transfer RNA (tRNA), and eight ribosomal RNA (rRNA) genes. Our repeat analyses revealed that L. apetalum harboured 20 direct repeats, 16 palindrome repeats, 30 tandem repeats, and 87 simple sequence repeats, whereas, L. perfoliatum had 15 direct repeats, 20 palindrome repeats, four reverse repeats, 21 tandem repeats, and 98 simple sequence repeats. Using syntenic analysis, we also revealed a high degree of sequence similarity within the coding regions of Lepidium medicinal plant cp genomes, and a high degree of divergence among the intergenic spacers. Pairwise alignment and single-nucleotide polymorphism (SNP) examinations further revealed certain Lepidium -specific gene fragments. Codon usage analysis showed that codon 14 was the most frequently used codon in the Lepidium coding sequences. Further, correlation investigations suggest that L. apetalum and L. perfoliatum originate from similar genetic backgrounds. Analysis of codon usage bias of Lepidium cp genome was strongly influenced by mutation and natural selection. We showed that L. apetalum and L. perfoliatum will likely enhance breeding, species recognition, phylogenetic evolution, and cp genetic engineering of the Lepidium medicinal plants.

Deoxyribonucleic Acid Barcoding for the Identification of Botanicals

The Natural Herbal Products industry uses botanicals or herbs as raw materials for production of herbal products or dietary supplements. Recently, the demand for natural herbal products has increased tremendously and this has led to adulteration and to counterfeit herbal products. The present chapter deals with currently used molecular methods from "simple" single genomic regions to high-throughput whole genome or transcriptome sequencing methods used in the identification of botanicals.

New distribution record and DNA barcoding of Sapria himalayana Griff. (Rafflesiaceae), a rare and endangered holoparasitic plant from Mizoram, India

Sapria himalayana Griff. is a rare and endangered holoparasitic plant that prefers a specific host (Tetrastigma sp.). It is one of the lesser-known and poorly understood plants facing threats of extinction owing to human interference in the evergreen forests of Mizoram. The flower is the only visible part of this endophyte and blooms from November to December. The plant was encountered for the first time in the evergreen forest near Rullam village in the Serchhip District of Mizoram, India. In the present study, DNA barcoding was used to identify the plants, and the internal transcribed spacer 2 (ITS2) region of S. himalayana was amplified and sequenced. The ITS2 sequence could accurately identify up to the species level for this endangered species. The absence of the ribulose-biphosphate carboxylase gene (rbcL) region in the genome supports its holoparasitic nature. Hence, DNA barcoding can help in taxonomic and biodiversity research and aid in selecting taxa for various molecular ecology and population genetics studies. The phylogenetic tree was analyzed using the maximum-likelihood method, and our findings showed that species from different families were clearly discriminated in a phylogenetic tree. To the best of our knowledge, this is the first report of DNA barcoding using ITS2 region of S. himalayana from Mizoram, India.

Phylogenomics, plastome structure and species identification in Mahonia (Berberidaceae)

BACKGROUND

Elucidating the phylogenetic relationships within species-rich genera is essential but challenging, especially when lineages are assumed to have been going through radiation events. Mahonia Nutt. (Berberidaceae) is a genus with cosmopolitan distribution, comprising approximately 100 species, two of which are known as Caulis Mahoniae (M. bealei and M. fortunei) with crucial pharmacological significance in Chinese herbal medicine. Mahonia is a taxonomically challenging genus, and intrageneric phylogenetic relationships still need to be explored using genome data. Universal DNA barcodes and floral morphological attributes have limited discriminatory power in Mahonia.

RESULTS

We sequenced 17 representative plastomes and integrated three published plastome data together to conduct comparative and phylogenetic analyses. We found that Mahonia and Berberis share a large IR expansion (~ 12 kb), which is recognized as a typical character of Berberideae. Repeated sequences are revealed in the species of Mahonia, which are valuable for further population genetic studies. Using a comparative plastome analysis, we determined eight hypervariable regions whose discriminative power is comparable to that of the whole plastid genomes. The incongruence of the ITS and the plastome tree topologies may be ascribed to ancestral hybridization events and/or to incomplete lineage sorting. In addition, we suggest that leaf epidermal characters could help to distinguish closely related species in Mahonia.

CONCLUSIONS

We propose an integrative approach combining special barcodes and micromorphological traits to circumscribe Mahonia species. The results cast a new light on the development of an integrative method for accurate species circumscription and provide abundant genetic resources for further research on Mahonia.

A manager's guide to using eDNA metabarcoding in marine ecosystems

Environmental DNA (eDNA) metabarcoding is a powerful tool that can enhance marine ecosystem/biodiversity monitoring programs. Here we outline five important steps managers and researchers should consider when developing eDNA monitoring program: (1) select genes and primers to target taxa; (2) assemble or develop comprehensive barcode reference databases; (3) apply rigorous site occupancy based decontamination pipelines; (4) conduct pilot studies to define spatial and temporal variance of eDNA; and (5) archive samples, extracts, and raw sequence data. We demonstrate the importance of each of these considerations using a case study of eDNA metabarcoding in the Ports of Los Angeles and Long Beach. eDNA metabarcoding approaches detected 94.1% (16/17) of species observed in paired trawl surveys while identifying an additional 55 native fishes, providing more comprehensive biodiversity inventories. Rigorous benchmarking of eDNA metabarcoding results improved ecological interpretation and confidence in species detections while providing archived genetic resources for future analyses. Well designed and validated eDNA metabarcoding approaches are ideally suited for biomonitoring applications that rely on the detection of species, including mapping invasive species fronts and endangered species habitats as well as tracking range shifts in response to climate change. Incorporating these considerations will enhance the utility and efficacy of eDNA metabarcoding for routine biomonitoring applications.

Quantifying the uncertainty of assembly-free genome-wide distance estimates and phylogenetic relationships using subsampling

Computing distance between two genomes without alignments or even access to assemblies has many downstream analyses. However, alignment-free methods, including in the fast-growing field of genome skimming, are hampered by a significant methodological gap. While accurate methods (many k-mer-based) for assembly-free distance calculation exist, measuring the uncertainty of estimated distances has not been sufficiently studied. In this paper, we show that bootstrapping, the standard non-parametric method of measuring estimator uncertainty, is not accurate for k-mer-based methods that rely on k-mer frequency profiles. Instead, we propose using subsampling (with no replacement) in combination with a correction step to reduce the variance of the inferred distribution. We show that the distribution of distances using our procedure matches the true uncertainty of the estimator. The resulting phylogenetic support values effectively differentiate between correct and incorrect branches and identify controversial branches that change across alignment-free and alignment-based phylogenies reported in the literature.