Incomplete estimates of genetic diversity within species: Implications for DNA barcoding

DNA barcoding insufficiently identifies European wild bees (Hymenoptera, Anthophila) due to undefined species diversity, genus-specific barcoding gaps and database errors

Recent declines in insect abundances, especially populations of wild pollinators, pose a threat to many natural and agricultural ecosystems. Traditional species monitoring relies on morphological character identification and is inadequate for efficient and standardized surveys. DNA barcoding has become a standard approach for molecular identification of organisms, aiming to overcome the shortcomings of traditional biodiversity monitoring. However, its efficacy depends on the completeness of reference databases. Large DNA barcoding efforts are (almost entirely) lacking in many European countries and such patchy data limit Europe-wide analyses of precisely how to apply DNA barcoding in wild bee identification. Here, we advance towards an effective molecular identification of European wild bees. We conducted a high-effort survey of wild bees at the junction of central and southern Europe and DNA barcoded all collected morphospecies. For global analyses, we complemented our DNA barcode dataset with all relevant European species and conducted global analyses of species delimitation, general and genus-specific barcoding gaps and examined the error rate in DNA data repositories. We found that (i) a sixth of all specimens from Slovenia could not be reliably identified, (ii) species delimitation methods show numerous systematic discrepancies, (iii) there is no general barcoding gap across all bees and (iv) the barcoding gap is genus specific, but only after curating for errors in DNA data repositories. Intense sampling and barcoding efforts in underrepresented regions and strict curation of DNA barcode repositories are needed to enhance the use of DNA barcoding for the identification of wild bees.

Molecular evidence for a new endemic species of Acartia (Copepoda, Calanoida) from the Southeast Pacific coast

The loss of biodiversity in marine populations is one of the consequences of the increased events of extreme environmental conditions in the oceans, which can condition the persistence of populations to future scenarios of climate change. Therefore, it is extremely necessary to explore and monitor the genetic diversity of natural populations. In the Southeast Pacific Ocean (SEPO), specifically on the coast of Chile, the presence of the copepod Acartia tonsa has been indicated solely using morphological evidence, due to the absence of genetic information. In the present work, the genetic diversity, population structure and phylogenetic position within the genus Acartia, of populations identified morphologically as A. tonsa, was evaluated by amplification of the mitochondrial cytochrome c oxidase subunit I and nuclear marker 18 s. Our results showed that the populations identified as A. tonsa correspond to a new monophyletic group endemic to SEPO (GMYC = 1.00; PTP = 0.95). The populations showed moderate to high genetic diversity with an incipient structuring between populations and biogeographic zones. Our results suggest that despite the homogenizing effect of the Humboldt Current, isolation by distance and contrasting environmental conditions at different geographic scales have an important influence on the genetic diversity of zooplankton in the SEPO region.

Unique among high passes: Phylogenetic inferences from DNA barcoding of the butterfly fauna of Ladakh Trans-Himalaya, India

The butterfly assemblage of Ladakh Trans-Himalaya demands a thorough analysis of their population genetic structure owing to their typical biogeographic affinity and their adaptability to extreme cold-desert climates. No such effort has been taken till date, and in this backdrop, we created a barcode reference library of 60 specimens representing 23 species. Barcodes were generated from freshly collected leg samples using the Sanger sequencing method, followed by phylogenetic clade analyses and divergence calculation. Our data represents 22% of Ladakh's Rhopaloceran fauna with the novel barcode submission for six species, including one Schedule II species, Paralasa mani . Contrary to the 3% threshold rule, the interspecific divergence between two species pairs of typical mountain genus Hyponephele and Karanasa was found to be 2.3% and 2.2%, respectively. The addition of conspecific global barcodes revealed that most species showed little increase in divergence value, while a two-fold increase was noted in a few species. Bayesian clade clustering outcomes largely aligned with current morphological classifications, forming monophyletic clades of conspecific barcodes, with only minor exceptions observed for the taxonomically complicated genus Polyommatus and misidentified records of Aulocera in the database. We also observed variations within the same phylogenetic clades forming nested lineages, which may be attributed to the taxonomic intricacies present at the subspecies level globally, mostly among Eurasian species. Overall, our effort not only substantiated the effectiveness of DNA Barcoding for the identification and conservation of this climatically vulnerable assemblage but also highlighted the significance of deciphering the unique genetic composition among this geographically isolated population of Ladakh butterflies.

Mitochondrial genome comparison and phylogenetic position of Fannia pusio among the Calyptratae flies

Fannia pusio, the chicken dung fly species, remains unexplored despite its forensic, sanitary, and veterinary importance in the Nearctic and Neotropical regions. In this study, we obtained the complete mitochondrial genome of Fannia pusio for the first time using next-generation sequencing. We compared it with previously published mitogenomes of the genus from the Palearctic region, and its phylogenetic position was studied based on the concatenated protein-coding genes (PCGs) dataset of Calyptratae flies. The circular mitochondrial genome of F. pusio is 16,176 bp in length, with a high A + T content (78.3%), whose gene synteny, codon usage analysis, and amino acid frequency are similar to previously reported Fannia mitogenomes. All PCGs underwent purifying selection except the nad2 gene. Interspecific K2P distances of PCGs of Fannia yielded an average of 12.4% (8.1%-21.1%). The Fannia genus is monophyletic and closely related to Muscidae based on molecular data. Further taxonomic sampling is required to deep into the phylogenetic relationships of the originally proposed species-groups and subgroups within the genus. These results provide a valuable dataset for studying the mitochondrial genome evolution and a resource for the taxonomy and systematics of Fannia.

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.

Deep-sea meiofaunal communities in the south-eastern Levantine basin and their shaping factors - Morphological-taxonomy-free metabarcoding approach

The <3% dissimilar Amplicon Sequence Variant (ASV) clusters of the 18S-V4 barcode were used as species-proxies for the evaluation of ASV composition and ASV diversity indices characterizing the hitherto poorly investigated meiofaunal communities of the south-eastern part of the Levantine basin. Accompanied by abundance measurements, the relationships of these characteristics with sedimentary and bottom terrain parameters were interpreted. The construction of community composition profiles, namely ASVs' list and their estimated abundances, was done using our previously established procedure (Harbuzov et al., 2022, Marine Genomics 65, 100980), combining metabarcoding with sample reads normalization by the abundance of hard-bodied meiofaunal taxa. The study province included the 54-1418 m depth range, across vertical sub-bottom horizons ranging 0-17 cm. Oxygen, hydrogen sulfide and methane concentrations in the pore water, as well as sediment grain size spectra and sedimentary protein and carbohydrate levels, were measured, followed by an evaluation of their involvement in the shaping of the meiofaunal communities' characteristics. Community composition was generally site-and-horizon dependent and its abundance decreased with increasing bottom depth and across sub-bottom horizons, typical to benthic habitats which are nourished by organic carbon from the euphotic zone. The relatively sharply inclined continental slope bottom located in the northern part of the Israeli coast was an exception. Its meiofaunal community characteristics were speculated to be affected by intensive sediment mixing and lateral transport of food from the shelf, in addition to the effect of the euphotic zone-originated food sources.

A Measure of the DNA Barcode Gap for Applied and Basic Research

DNA barcoding has largely established itself as a mainstay for rapid molecular taxonomic identification in both academic and applied research. The use of DNA barcoding as a molecular identification method depends on a "DNA barcode gap"-the separation between the maximum within-species difference and the minimum between-species difference. Previous work indicates the presence of a gap hinges on sampling effort for focal taxa and their close relatives. Furthermore, both theory and empirical work indicate a gap may not occur for related pairs of biological species. Here, we present a novel evaluation approach in the form of an easily calculated set of nonparametric metrics to quantify the extent of proportional overlap in inter- and intraspecific distributions of pairwise differences among target species and their conspecifics. The metrics are based on a simple count of the number of overlapping records for a species falling within the bounds of maximum intraspecific distance and minimum interspecific distance. Our approach takes advantage of the asymmetric directionality inherent in pairwise genetic distance distributions, which has not been previously done in the DNA barcoding literature. We apply the metrics to the predatory diving beetle genus Agabus as a case study because this group poses significant identification challenges due to its morphological uniformity despite both relative sampling ease and well-established taxonomy. Results herein show that target species and their nearest neighbor species were found to be tightly clustered and therefore difficult to distinguish. Such findings demonstrate that DNA barcoding can fail to fully resolve species in certain cases. Moving forward, we suggest the implementation of the proposed metrics be integrated into a common framework to be reported in any study that uses DNA barcoding for identification. In so doing, the importance of the DNA barcode gap and its components for the success of DNA-based identification using DNA barcodes can be better appreciated.

taxalogue: a toolkit to create comprehensive CO1 reference databases

Background

Taxonomic identification through DNA barcodes gained considerable traction through the invention of next-generation sequencing and DNA metabarcoding. Metabarcoding allows for the simultaneous identification of thousands of organisms from bulk samples with high taxonomic resolution. However, reliable identifications can only be achieved with comprehensive and curated reference databases. Therefore, custom reference databases are often created to meet the needs of specific research questions. Due to taxonomic inconsistencies, formatting issues, and technical difficulties, building a custom reference database requires tremendous effort. Here, we present taxalogue, an easy-to-use software for creating comprehensive and customized reference databases that provide clean and taxonomically harmonized records. In combination with extensive geographical filtering options, taxalogue opens up new possibilities for generating and testing evolutionary hypotheses.

Methods

taxalogue collects DNA sequences from several online sources and combines them into a reference database. Taxonomic incongruencies between the different data sources can be harmonized according to available taxonomies. Dereplication and various filtering options are available regarding sequence quality or metadata information. taxalogue is implemented in the open-source Ruby programming language, and the source code is available at https://github.com/nwnoll/taxalogue. We benchmark four reference databases by sequence identity against eight queries from different localities and trapping devices. Subsamples from each reference database were used to compare how well another one is covered.

Results

taxalogue produces reference databases with the best coverage at high identities for most tested queries, enabling more accurate, reliable predictions with higher certainty than the other benchmarked reference databases. Additionally, the performance of taxalogue is more consistent while providing good coverage for a variety of habitats, regions, and sampling methods. taxalogue simplifies the creation of reference databases and makes the process reproducible and transparent. Multiple available output formats for commonly used downstream applications facilitate the easy adoption of taxalogue in many different software pipelines. The resulting reference databases improve the taxonomic classification accuracy through high coverage of the query sequences at high identities.

Morphological and genetic mechanisms underlying the plasticity of the coral Porites astreoides across depths in Bermuda

The widespread decline of shallow-water coral reefs has fueled interest in assessing whether mesophotic reefs can act as refugia replenishing deteriorated shallower reefs through larval exchange. Here we explore the morphological and molecular basis facilitating survival of planulae and adults of the coral Porites astreoides (Lamarck, 1816; Hexacorallia: Poritidae) along the vertical depth gradient in Bermuda. We found differences in micro-skeletal features such as bigger calyxes and coarser surface of the skeletal spines in shallow corals. Yet, tomographic reconstructions reveal an analogous mineral distribution between shallow and mesophotic adults, pointing to similar skeleton growth dynamics. Our study reveals patterns of host genetic connectivity and minimal symbiont depth-zonation across a broader depth range than previously known for this species in Bermuda. Transcriptional variations across life stages showed different regulation of metabolism and stress response functions, unraveling molecular responses to environmental conditions at different depths. Overall, these findings increase our understanding of coral acclimatory capability across broad vertical gradients, ultimately allowing better evaluation of the refugia potential of mesophotic reefs.

How we study cryptic species and their biological implications: A case study from marine shelled gastropods

Methodological and biological considerations are intertwined when studying cryptic species. A potentially large component of modern biodiversity, the frequency of cryptic species among taxonomic groups is not well documented. The term "cryptic species" is imprecisely used in scientific literature, causing ambiguity when interpreting their evolutionary and ecological significance. This study reviews how cryptic species have been defined, discussing implications for taxonomy and biology, and explores these implications with a case study based on recently published literature on extant shelled marine gastropods. Reviewed gastropods were recorded by species. Records of cryptic gastropods were presented by authors with variable levels of confidence but were difficult to disentangle from inherent biases in the study effort. These complexities notwithstanding, most gastropod species discussed were not cryptic. To the degree that this review's sample represents extinct taxa, the results suggest that a high proportion of shelled marine gastropod species are identifiable for study in the fossil record. Much additional work is needed to provide a more adequate understanding of the relative frequency of cryptic species in shelled marine gastropods, which should start with more explicit definitions and targeted case studies.

The role and risks of selective adaptation in extreme coral habitats

The alarming rate of climate change demands new management strategies to protect coral reefs. Environments such as mangrove lagoons, characterized by extreme variations in multiple abiotic factors, are viewed as potential sources of stress-tolerant corals for strategies such as assisted evolution and coral propagation. However, biological trade-offs for adaptation to such extremes are poorly known. Here, we investigate the reef-building coral Porites lutea thriving in both mangrove and reef sites and show that stress-tolerance comes with compromises in genetic and energetic mechanisms and skeletal characteristics. We observe reduced genetic diversity and gene expression variability in mangrove corals, a disadvantage under future harsher selective pressure. We find reduced density, thickness and higher porosity in coral skeletons from mangroves, symptoms of metabolic energy redirection to stress response functions. These findings demonstrate the need for caution when utilizing stress-tolerant corals in human interventions, as current survival in extremes may compromise future competitive fitness.

Bioinformatics approaches and big data analytics opportunities in improving fisheries and aquaculture

Aquaculture has witnessed an excellent growth rate during the last two decades and offers huge potential to provide nutritional as well as livelihood security. Genomic research has contributed significantly toward the development of beneficial technologies for aquaculture. The existing high throughput technologies like next-generation technologies generate oceanic data which requires extensive analysis using appropriate tools. Bioinformatics is a rapidly evolving science that involves integrating gene based information and computational technology to produce new knowledge for the benefit of aquaculture. Bioinformatics provides new opportunities as well as challenges for information and data processing in new generation aquaculture. Rapid technical advancements have opened up a world of possibilities for using current genomics to improve aquaculture performance. Understanding the genes that govern economically relevant characteristics, necessitates a significant amount of additional research. The various dimensions of data sources includes next-generation DNA sequencing, protein sequencing, RNA sequencing gene expression profiles, metabolic pathways, molecular markers, and so on. Appropriate bioinformatics tools are developed to mine the biologically relevant and commercially useful results. The purpose of this scoping review is to present various arms of diverse bioinformatics tools with special emphasis on practical translation to the aquaculture industry.

Similar Ones Are Not Related and Vice Versa—New Dendronotus Taxa (Nudibranchia: Dendronotidae) from the North Atlantic Ocean Provide a Platform for Discussion of Global Marine Biodiversity Patterns

One new species of the genus Dendronotus (Nudibranchia: Dendronotidae) is described from Norway and Northern Ireland, as well as from the adjacent North Sea, and one new subspecies of Dendronotus arcticus is described from Norway by applying a combination of fine-scale morphological and molecular phylogenetic data. The present case demonstrates multilevel morphological and molecular similarities and differences considering on the one hand a grouping of three similar looking sympatric taxa (D. yrjargul, D. arcticus gartensis n. subsp. and D. keatleyae n. sp.), and on the other hand two different looking apparently allopatric subspecies (D. arcticus arcticus and D. arcticus gartensis n. subsp.). The type species of the genus, D. frondosus, which is the commonest dendronotid in Norway and the United Kingdom, consistently demonstrates substantial molecular and fine-scale morphological differences from D. keatleyae n. sp. The present study, apart from providing purely taxonomic information, also provides new data for a broad discussion of global biodiversity patterns.

Extensive Cryptic Diversity and Ecological Associations Uncovered among Mexican and Global Collections of Naegleria and Vermamoeba Species by 18S Ribosomal DNA, Internal Transcribed Spacer, and Cytochrome Oxidase Subunit I Sequence Analysis

Free-living amoebae (FLA) are phagocytic protists that play crucial roles in microbial communities as significant microbial grazers. However, our current knowledge of their diversity, ecology, and population genetic structures is marginal due to the shallow and biased sampling of ecosystems and the use of few, poorly resolving molecular markers. Thirty-two FLA were isolated from soil and water samples collected across representative ecosystems of the State of Morelos in Central Mexico, including the drinking water distribution system (DWDS) from the state capital. We classified our isolates as members of Acanthamoeba, Vermamoeba, Naegleria, and Tetramitus by 18S ribosomal DNA (rDNA) sequencing. Vermamoeba isolates were recovered exclusively from the DWDS samples. In contrast, Naegleria strains displayed a broad distribution in soil and water samples across the natural ecosystems. We used a combination of phylogenetic and population genetic analyses of internal transcribed spacer (ITS) and cytochrome oxidase subunit I (COI) sequences from our isolates and a comprehensive set of reference sequences to analyze the currently known diversity of Naegleria spp. Significant associations were uncovered between the most prevalent lineages of Naegleria and Vermamoeba and broad ecological and geographical variables at regional and global levels. The population structure and cryptic diversity within the Naegleria galeacystis-Naegleria americana and Vermamoeba vermiformis species complexes were thoroughly analyzed. Our results prove that the genus Vermamoeba, which was previously thought to consist of only one species, actually encompasses at least seven widely distributed species, as indicated by consistent evidence from Bayesian phylogenetics, two species-delimitation programs, and population genetics analyses. IMPORTANCE Our study sheds new light on the population genetic structure of V. vermiformis and diverse Naegleria species. Using improved molecular markers and advanced analytical approaches, we discovered that N. americana, previously considered a single species, actually contains multiple distinct lineages, as revealed by COI sequencing. These lineages are highly differentiated, with little gene flow between them. Our findings demonstrate that the genus Vermamoeba holds multiple cryptic species, requiring a significant taxonomic revision in light of multilocus sequence analyses. These results advance our understanding of the ecology, molecular systematics, and biogeography of these genera and species complexes at both regional and global scales. This study has significant implications for diagnosing amoebal infections and evaluating health risks associated with FLA in domestic and recreational waters.

Demographic changes and life-history strategies predict the genetic diversity in crabs

Uncovering what predicts genetic diversity (GD) within species can help us access the status of populations and their evolutionary potential. Traits related to effective population size show a proportional association to GD, but evidence supports life-history strategies and habitat as the drivers of GD variation. Instead of investigating highly divergent taxa, focusing on one group could help to elucidate the factors influencing the GD. Additionally, most empirical data is based on vertebrate taxa; therefore, we might be missing novel patterns of GD found in neglected invertebrate groups. Here, we investigated the predictors of the GD in crabs (Brachyura) by compiling the most comprehensive cytochrome c oxidase subunit I (COI) available. Eight predictor variables were analysed across 150 species (16 992 sequences) using linear models (multiple linear regression) and comparative methods (PGLS). Our results indicate that population size fluctuation represents the most critical trait predicting GD, with species that have undergone bottlenecks followed by population expansion showing lower GD. Egg size, pelagic larval duration and habitat might play a role probably because of their association with how species respond to disturbances. Ultimately, K-strategists that have undergone bottlenecks are the species showing lower GD. Some variables do not show an association with GD as expected, most likely due to the taxon-specific role of some predictors, which should be considered in further investigations and generalizations. This work highlights the complexity underlying the predictors of GD and adds results from a marine invertebrate group to the current understanding of this topic.

DNA Barcoding of Fish Species Diversity in Guizhou, China

Guizhou is an important ecological barrier in the upper reaches of the Yangtze River and the Pearl River basins with abundant fish species. However, fish from these regions are threatened by anthropogenic activities, including overfishing and habitat destruction. Here, we assessed the fish diversity including more than half of the species from the region using DNA barcoding (partial sequence of cytochrome c oxidase subunit I (COI) gene). We obtained 800 mitochondrial COI barcode sequences from 82 genera, 18 families and 8 orders of fishes. The average Kimura two-parameter (K2P) distances within species and genera were 0.35% and 5.44%, respectively. The average interspecific distance was 15.54 times higher than the mean intraspecific distance. Moreover, DNA barcodes revealed 175 operational taxonomic units (OTUs) based on consensus demarcation schemes. Barcoding gaps were detected in 94.81% of morphospecies. Three fish species (Schistura fasciolata, Vanmanenia pingchowensis, and Misgurnus dabryanus) have considerable intraspecific variability, and each was divided into multiple molecular operational taxonomic units (MOTUs) using molecular definition methods (Automatic Barcode Gap Discovery, Refined Single Linkage, General Mixed Yule Coalescent, and Poisson Tree Processes), possibly indicating the occurrence of cryptic species. Altogether, our study reveals the complex diversity of fish species in Guizhou Province, serving as a reference for the conservation and monitoring of fish species in this region.

VLF: An R package for the analysis of very low frequency variants in DNA sequences

Here, we introduce VLF, an R package to determine the distribution of very low frequency variants (VLFs) in nucleotide and amino acid sequences for the analysis of errors in DNA sequence records. The package allows users to assess VLFs in aligned and trimmed protein-coding sequences by automatically calculating the frequency of nucleotides or amino acids in each sequence position and outputting those that occur under a user-specified frequency (default of p = 0.001). These results can then be used to explore fundamental population genetic and phylogeographic patterns, mechanisms and processes at the microevolutionary level, such as nucleotide and amino acid sequence conservation. Our package extends earlier work pertaining to an implementation of VLF analysis in Microsoft Excel, which was found to be both computationally slow and error prone. We compare those results to our own herein. Results between the two implementations are found to be highly consistent for a large DNA barcode dataset of bird species. Differences in results are readily explained by both manual human error and inadequate Linnean taxonomy (specifically, species synonymy). Here, VLF is also applied to a subset of avian barcodes to assess the extent of biological artifacts at the species level for Canada goose (Branta canadensis), as well as within a large dataset of DNA barcodes for fishes of forensic and regulatory importance. The novelty of VLF and its benefit over the previous implementation include its high level of automation, speed, scalability and ease-of-use, each desirable characteristics which will be extremely valuable as more sequence data are rapidly accumulated in popular reference databases, such as BOLD and GenBank.

Species discrimination in Schima (Theaceae): Next-generation super-barcodes meet evolutionary complexity

Plastid genome and nrDNA arrays, proposed recently as "super barcodes", might provide additional discriminatory power and overcome the limitations of traditional barcoding loci, yet super barcodes need to be tested for their effectiveness in more plant groups. Morphological homoplasy among Schima species makes the genus a model for testing the efficacy of super barcodes. In this study, we generated multiple datasets comprising standard DNA barcodes (matK, rbcL, trnH-psbA, nrITS) and super-barcodes (plastid genome, nrDNA arrays) across 58 individuals from 12 out of 13 species of Schima from China. No samples were correctly assigned to species using standard DNA barcodes and nrDNA arrays, while only 27.27% of species with multiple accessions were distinguished using the plastid genome and its partitioned datasets-the lowest estimated rate of super barcode success in the literature so far. For Schima and other taxa with similarly recently divergence and low levels of genetic variation, incomplete lineage sorting, hybridization, or taxonomic oversplitting are all possible causes of the failure. Taken together, our study suggests that by no means are super barcodes immune to the challenges imposed by evolutionary complexity. We therefore call for developing multi-locus nuclear markers for species discrimination in plant groups.

DNA Barcoding Medicinal Plant Species from Indonesia

Over the past decade, plant DNA barcoding has emerged as a scientific breakthrough and is often used to help with species identification or as a taxonomical tool. DNA barcoding is very important in medicinal plant use, not only for identification purposes but also for the authentication of medicinal products. Here, a total of 61 Indonesian medicinal plant species from 30 families and a pair of ITS2, matK, rbcL, and trnL primers were used for a DNA barcoding study consisting of molecular and sequence analyses. This study aimed to analyze how the four identified DNA barcoding regions (ITS2, matK, rbcL, and trnL) aid identification and conservation and to investigate their effectiveness for DNA barcoding for the studied species. This study resulted in 212 DNA barcoding sequences and identified new ones for the studied medicinal plant species. Though there is no ideal or perfect region for DNA barcoding of the target species, we recommend matK as the main region for Indonesian medicinal plant identification, with ITS2 and rbcL as alternative or complementary regions. These findings will be useful for forensic studies that support the conservation of medicinal plants and their national and global use.

Lack of Statistical Rigor in DNA Barcoding Likely Invalidates the Presence of a True Species' Barcode Gap

DNA barcoding has been largely successful in satisfactorily exposing levels of standing genetic diversity for a wide range of taxonomic groups through the employment of only one or a few universal gene markers. However, sufficient coverage of geographically-broad intra-specific haplotype variation within genomic databases like the Barcode of Life Data Systems (BOLD) and GenBank remains relatively sparse. As reference sequence libraries continue to grow exponentially in size, there is now the need to identify novel ways of meaningfully analyzing vast amounts of available DNA barcode data. This is an important issue to address promptly for the routine tasks of specimen identification and species discovery, which have seen broad adoption in areas as diverse as regulatory forensics and resource conservation. Here, it is demonstrated that the interpretation of DNA barcoding data is lacking in statistical rigor. To highlight this, focus is set specifically on one key concept that has become a household name in the field: the DNA barcode gap. Arguments outlined herein specifically center on DNA barcoding in animal taxa and stem from three angles: (1) the improper allocation of specimen sampling effort necessary to capture adequate levels of within-species genetic variation, (2) failing to properly visualize intra-specific and interspecific genetic distances, and (3) the inconsistent, inappropriate use, or absence of statistical inferential procedures in DNA barcoding gap analyses. Furthermore, simple statistical solutions are outlined which can greatly propel the use of DNA barcoding as a tool to irrefutably match unknowns to knowns on the basis of the barcoding gap with a high degree of confidence. Proposed methods examined herein are illustrated through application to DNA barcode sequence data from Canadian Pacific fish species as a case study.

DNA barcoding and phylogeography of the Hoplias malabaricus species complex

Hoplias malabaricus (Bloch, 1794) is a carnivorous fish species widely distributed from northern to southern South America. This taxon is believed to be a good model for the investigation of biogeographic events that shape the ichthyofauna evolution in the Neotropical freshwater systems. However, many studies have revealed that H. malabaricus hides a species complex that hampers its taxonomic identity and limit its practical value for evolutionary and biogeographic studies. In this paper, we used the mitochondrial gene cytochrome c oxidase subunit I (COI) to delimit cryptic species and explore the phylogeography of H. malabaricus sensu stricto. We found genetic evidence for putative new species in the genus Hoplias and showed that H. malabaricus (Bloch, 1794) is a major clade assigned to barcode index number (BIN) BOLD:ABZ3047. This species is structured in six subpopulations differentiated by high Fst values and restricts gene flow. The subpopulations of the São Francisco/East Atlantic/Eastern Northeast Atlantic/Parnaíba/Itapecuru River basins and Tapajós River Basin were the most differentiated and showed demographic fluctuations. The present distributional pattern is most likely explained through a scenario from the Pleistocene.

Population Structure and Genetic Diversity of the Pepper Weevil (Coleoptera: Curculionidae) Using the COI Barcoding Region

The pepper weevil Anthonomus eugenii Cano (Coleoptera: Curculionidae) is a pest of economic importance for Capsicum species pepper in North America that attacks the reproductive structures of the plant. The insect is distributed across Mexico, the United States, and the Caribbean, and is occasionally found during the pepper growing season in southern Ontario, Canada. Continuous spread of the insect to new areas is partially the result of global pepper trade. Here, we describe the genetic diversity of the pepper weevil using the mitochondrial COI barcoding region across most of its geographic range. In this study, 44 (H1-H44) highly similar haplotypes were identified, the greatest number of haplotypes and haplotype diversity were observed among specimens from its native Mexico, followed by specimens from the United States. Unlike Mexico, a low haplotype diversity was found among specimens from Canada, the Dominican Republic, Italy, and the Netherlands. Out of these 44 haplotypes, 29 are reported for the first time. Haplotype diversity in the Canadian population suggests either multiple and continuous introductions of the pepper weevil into this area or a single introduction of genetically diverse individuals. We discuss the importance of such population genetic data in tailoring pepper weevil management programs, using Canada as an example.

DNA barcode reference library construction and genetic diversity and structure analysis of Amomum villosum Lour. (Zingiberaceae) populations in Guangdong Province

Background

Amomum villosum Lour. is the plant that produces the famous traditional Chinese medicine Amomi Fructus. Frequent habitat destruction seriously threatens A. villosum germplasm resources. Genetic diversity is very important to the optimization of germplasm resources and population protection, but the range of inherited traits within A. villosum is unclear. In this study, we analyzed the genetic diversity and genetic structures of A. villosum populations in Guangdong and constructed a local reference DNA barcode library as a resource for conservation efforts.

Methods

DNA barcoding and Inter-Simple Sequence Repeat (ISSR) markers were used to investigate the population genetics of A. villosum. Five universal DNA barcodes were amplified and used in the construction of a DNA barcode reference library. Parameters including percentage of polymorphic sites (PPB), number of alleles (Na), effective number of alleles (Ne), Nei’s gene diversity index (H), and Shannon’s polymorphism information index (I) were calculated for the assessment of genetic diversity. Genetic structure was revealed by measuring Nei’s gene differentiation coefficient (Gst), total population genetic diversity (Ht), intra-group genetic diversity (Hs), and gene flow (Nm). Analysis of molecular variance (AMOVA), Mantel tests, unweighted pair-group method with arithmetic mean (UPGMA) dendrogram, and principal co-ordinates (PCoA) analysis were used to elucidate the genetic differentiation and relationship among populations.

Results

A total of 531 sequences were obtained from the five DNA barcodes with no variable sites from any of the barcode sequences. A total of 66 ISSR bands were generated from A. villosum populations using the selected six ISSR primers; 56 bands, 84.85% for all the seven A. villosum populations were polymorphic. The A. villosum populations showed high genetic diversity (H = 0.3281, I = 0.4895), whereas the gene flow was weak (Nm = 0.6143). Gst (0.4487) and AMOVA analysis indicated that there is obvious genetic differentiation amongA. villosum populations and more genetic variations existed within each population. The genetic relationship of each population was relatively close as the genetic distances were between 0.0844 and 0.3347.

Opportunities and challenges of macrogenetic studies

The rapidly emerging field of macrogenetics focuses on analysing publicly accessible genetic datasets from thousands of species to explore large-scale patterns and predictors of intraspecific genetic variation. Facilitated by advances in evolutionary biology, technology, data infrastructure, statistics and open science, macrogenetics addresses core evolutionary hypotheses (such as disentangling environmental and life-history effects on genetic variation) with a global focus. Yet, there are important, often overlooked, limitations to this approach and best practices need to be considered and adopted if macrogenetics is to continue its exciting trajectory and reach its full potential in fields such as biodiversity monitoring and conservation. Here, we review the history of this rapidly growing field, highlight knowledge gaps and future directions, and provide guidelines for further research.

The InBIO Barcoding Initiative Database: DNA barcodes of Portuguese Hemiptera 01

Background

The InBIO Barcoding Initiative (IBI) Hemiptera 01 dataset contains records of 131 specimens of Hemiptera. Most specimens have been morphologically identified to species or subspecies level and represent 88 species in total. The species of this dataset correspond to about 7.3% of continental Portuguese hemipteran species diversity. All specimens were collected in continental Portugal. Sampling took place from 2015 to 2019 and specimens are deposited in the IBI collection at CIBIO, Research Center in Biodiversity and Genetic Resources.

New information

This dataset increases the knowledge on the DNA barcodes and distribution of 88 species of Hemiptera from Portugal. Six species, from five different families, were new additions to the Barcode of Life Data System (BOLD), with another twenty five species barcodes' added from under-represented taxa in BOLD. All specimens have their DNA barcodes publicly accessible through BOLD online database and the distribution data can be accessed through the Global Biodiversity Information Facility (GBIF). Eutettixvariabilis and Fieberiellaflorii are recorded for the first time for Portugal and Siphantaacuta, an invasive species, previously reported from the Portuguese Azores archipelago, is recorded for the first time for continental Portugal.

Updating the bionomy and geographical distribution of Anopheles (Nyssorhynchus) albitarsis F: A vector of malaria parasites in northern South America

Anopheles albitarsis F is a putative species belonging to the Albitarsis Complex, recognized by rDNA, mtDNA, partial white gene, and microsatellites sequences. It has been reported from the island of Trinidad, Venezuela and Colombia, and incriminated as a vector of malaria parasites in the latter. This study examined mitochondrially encoded cytochrome c oxidase I (MT-CO1) sequences of An. albitarsis F from malaria-endemic areas in Colombia and Venezuela to understand its relations with other members of the Complex, revised and update the geographical distribution and bionomics of An. albitarsis F and explore hypotheses to explain its phylogenetic relationships and geographical expansion. Forty-five MT-CO1 sequences obtained in this study were analyzed to estimate genetic diversity and possible evolutionary relationships. Sequences generated 37 haplotypes clustered in a group where the genetic divergence of Venezuelan populations did not exceed 1.6% with respect to Colombian samples. Anopheles albitarsis F (π = 0.013) represented the most recent cluster located closer to An. albitarsis I (π = 0.009). Barcode gap was detected according to Albitarsis Complex lineages previously reported (threshold 0.014–0.021). Anopheles albitarsis F has a wide distribution in northern South America and might play an important role in the transmission dynamics of malaria due to its high expansion capacity. Future studies are required to establish the southern distribution of An. albitarsis F in Venezuela, and its occurrence in Guyana and Ecuador.

Genetic diversity of the Nubian ibex in Oman as revealed by mitochondrial DNA

The Nubian ibex (Capra nubiana) is patchily distributed across parts of Africa and Arabia. In Oman, it is one of the few free-ranging wild mammals found in the central and southern regions. Its population is declining due to habitat degradation, human expansion, poaching and fragmentation. Here, we investigated the population's genetic diversity using mitochondrial DNA (D-loop 186 bp and cytochrome b 487 bp). We found that the Nubian ibex in the southern region of Oman was more diverse (D-loop HD; 0.838) compared with the central region (0.511) and gene flow between them was restricted. We compared the genetic profiles of wild Nubian ibex from Oman with captive ibex. A Bayesian phylogenetic tree showed that wild Nubian ibex form a distinct clade independent from captive animals. This divergence was supported by high mean distances (D-loop 0.126, cytochrome b 0.0528) and high F _ST statistics (D-loop 0.725, cytochrome b 0.968). These results indicate that captive ibex are highly unlikely to have originated from the wild population in Oman and the considerable divergence suggests that the wild population in Oman should be treated as a distinct taxonomic unit. Further nuclear genetic work will be required to fully elucidate the degree of global taxonomic divergence of Nubian ibex populations.

Validation of a Triplex Quantitative Polymerase Chain Reaction Assay for Detection and Quantification of Traditional Protein Sources, Pisum sativum L. and Glycine max (L.) Merr., in Protein Powder Mixtures

Several botanicals have been traditionally used as protein sources, including the leguminous Pisum sativum L. and Glycine max (L.) Merr. While a rich history exists of cultivating these plants for their whole, protein-rich grain, modern use as powdered supplements present a new challenge in material authentication. The absence of clear morphological identifiers of an intact plant and the existence of long, complex supply chains behoove industry to create quick, reliable analytical tools to identify the botanical source of a protein product (many of which contain multiple sources). The utility of molecular tools for plant-based protein powder authentication is gaining traction, but few validated tools exist. Multiplex quantitative polymerase chain reaction (qPCR) can provide an economical means by which sources can be identified and relative proportions quantified. We followed established guidelines for the design, optimization, and validation of qPCR assay, and developed a triplex qPCR assay that can amplify and quantify pea and soy DNA targets, normalized by a calibrator. The assay was evaluated for analytical specificity, analytical sensitivity, efficiency, precision, dynamic range, repeatability, and reproducibility. We tested the quantitative ability of the assay using pea and soy DNA mixtures, finding exceptional quantitative linearity for both targets - 0.9983 (p < 0.0001) for soy and 0.9915 (p < 0.0001) for pea. Ratios based on mass of protein powder were also tested, resulting in non-linear patterns in data that suggested the requirement of further sample preparation optimization or algorithmic correction. Variation in fragment size within different lots of commercial protein powder samples was also analyzed, revealing low SD among lots. Ultimately, this study demonstrated the utility of qPCR in the context of protein powder mixtures and highlighted key considerations to take into account for commercial implementation.

Bioinformatic pipelines combining denoising and clustering tools allow for more comprehensive prokaryotic and eukaryotic metabarcoding

Environmental DNA metabarcoding is a powerful tool for studying biodiversity. However, bioinformatic approaches need to adjust to the diversity of taxonomic compartments targeted as well as to each barcode gene specificities. We built and tested a pipeline based on read correction with DADA2 allowing analysing metabarcoding data from prokaryotic (16S) and eukaryotic (18S, COI) life compartments. We implemented the option to cluster amplicon sequence variants (ASVs) into operational taxonomic units (OTUs) with swarm, a network-based clustering algorithm, and the option to curate ASVs/OTUs using LULU. Finally, taxonomic assignment was implemented via the Ribosomal Database Project Bayesian classifier (RDP) and BLAST. We validated this pipeline with ribosomal and mitochondrial markers using metazoan mock communities and 42 deep-sea sediment samples. The results show that ASVs and OTUs describe different levels of biotic diversity, the choice of which depends on the research questions. They underline the advantages and complementarity of clustering and LULU-curation for producing metazoan biodiversity inventories at a level approaching the one obtained using morphological criteria. While clustering removes intraspecific variation, LULU effectively removes spurious clusters, originating from errors or intragenomic variability. Swarm clustering affected alpha and beta diversity differently depending on genetic marker. Specifically, d-values > 1 appeared to be less appropriate with 18S for metazoans. Similarly, increasing LULU's minimum ratio level proved essential to avoid losing species in sample-poor data sets. Comparing BLAST and RDP underlined that accurate assignments of deep-sea species can be obtained with RDP, but highlighted the need for a concerted effort to build comprehensive, ecosystem-specific databases.

A Revised Phylogeny of the Mentha spicata Clade Reveals Cryptic Species

The genus Mentha is taxonomically and phylogenetically challenging due to complex genomes, polyploidization and an extensive historical nomenclature, potentially hiding cryptic taxa. A straightforward interpretation of phylogenetic relationships within the section Mentha is further hindered by dominant but outdated concepts on historically identified hybrid taxa. Mentha spicata is traditionally considered to be of hybrid origin, but the evidence for this is weak. Here, we aim to understand the phylogenetic relationships within the section Mentha using large sample sizes and to revisit the hybrid status and identity of M. spicata. We show that two of three traditional species in the subsection Spicatae are polyphyletic, as is the subsection as a whole, while the real number of cryptic species was underestimated. Compared to previous studies we present a fundamentally different phylogeny, with a basal split between M. spicata s.s. and M. longifolia s.s. Cluster analyses of morphological and genotypic data demonstrate that there is a dissociation between morphologically and genotypically defined groups of samples. We did not find any evidence that M. spicata is of hybrid origin, and we conclude its taxonomic status should be revised. The combination of genetic and phenotypic information is essential when evaluating hyperdiverse taxonomic groups.

Monitoring compliance of CITES lion bone exports from South Africa

From 2008 to 2018, South Africa permitted the export of captive-bred African lion (Panthera leo) skeletons to Southeast Asia under CITES Appendix II. Legal exports rose from approximately 50 individuals in 2008 to a maximum of 1,771 skeletons in 2016, and has led to ongoing concerns over possible laundering of non-lion, multiple-source and wild-sourced bones. South Africa is required under its obligations to CITES to employ mechanisms for monitoring and reporting trade, and to limit the potential for illegal trade and laundering of lion and other large felid bones. Monitoring tools for legal trade are critical to compliance with CITES. Here we evaluate the CITES-compliance procedure implemented by South Africa for export of lion bones and identify six essential general points for consideration in the implementation of animal export quota compliance protocols. We provide specific insight into the South African lion bone export monitoring system through: i) outlining the protocols followed; ii) assessing the utility of cranial morphology to identify species; iii) evaluating skeleton consignment weight as a monitoring tool; and iv) presenting molecular (DNA) species assignment and pairwise-comparative sample matching of individuals. We describe irregularities and illicit behaviour detected in the 2017 and 2018 lion bone quotas. Notably, we report that the compliance procedure successfully identified and prevented the attempted laundering of a tiger (P. tigris) skeleton in 2018. We emphasise the utility of mixed-method protocols for the monitoring of compliance in CITES Appendix II export quota systems.

High resolution DNA barcode library for European butterflies reveals continental patterns of mitochondrial genetic diversity

The study of global biodiversity will greatly benefit from access to comprehensive DNA barcode libraries at continental scale, but such datasets are still very rare. Here, we assemble the first high-resolution reference library for European butterflies that provides 97% taxon coverage (459 species) and 22,306 COI sequences. We estimate that we captured 62% of the total haplotype diversity and show that most species possess a few very common haplotypes and many rare ones. Specimens in the dataset have an average 95.3% probability of being correctly identified. Mitochondrial diversity displayed elevated haplotype richness in southern European refugia, establishing the generality of this key biogeographic pattern for an entire taxonomic group. Fifteen percent of the species are involved in barcode sharing, but two thirds of these cases may reflect the need for further taxonomic research. This dataset provides a unique resource for conservation and for studying evolutionary processes, cryptic species, phylogeography, and ecology.

Factors affecting the efficiency of molecular species delimitation in a species-rich insect family

In the context of global biodiversity loss, molecular species delimitation approaches can be very useful for accelerating species discovery through DNA taxonomy and inventory through DNA metabarcoding. In this study, the effect of some intrinsic factors on the efficiency of various single-marker species delimitation methods (fixed and variable nucleotide distance thresholds, ABGD, ASAP, GMYC, mPTP) was tested on more than 90 empirical data sets, derived from a set of 7,237 COI sequences attributed to 542 leaf beetles species (Coleoptera: Chrysomelidae). The considered factors were: (i) the number of haplotypes per species (as a proxy for genetic diversity), (ii) the geographic distance among conspecific collection localities (as a proxy of sampling width), (iii) the difficulty related to morphological identification of species, and (iv) the taxonomic rank. Distance-based methods, with on average more than 70% of match with morphological identification, outperformed those relying on phylogenetic trees, with less than 59%. A high number of haplotypes per species was found to have a negative effect on delimitation efficiency, whereas large geographic distances within species had a positive effect. All methods delimitations (except for GMYC) were significantly affected by the presence of species that are difficult to be identified, decreasing their efficiency. Finally, the only method influenced by the taxonomic rank of the data set was GMYC, showing lower efficiency in data sets at the genus than at higher levels. The observed biases we highlighted affecting efficiency could be accounted for when developing input data sets for species delimitation analyses to obtain a more reliable representation of biological diversity.

A DNA barcode-based survey of wild urban bees in the Loire Valley, France

The current decline of wild bees puts important ecosystem services such as pollination at risk. Both inventory and monitoring programs are needed to understand the causes of wild bee decline. Effective insect monitoring relies on both mass-trapping methods coupled with rapid and accurate identifications. Identifying wild bees using only morphology can be challenging, in particular, specimens from mass-trapped samples which are often in poor condition. We generated DNA barcodes for 2931 specimens representing 157 species (156 named and one unnamed species) and 28 genera. Automated cluster delineation reveals 172 BINs (Barcodes Index Numbers). A total of 36 species (22.93%) were found in highly urbanized areas. The majority of specimens, representing 96.17% of the species barcoded form reciprocally exclusive groups, allowing their unambiguous identification. This includes several closely related species notoriously difficult to identify. A total of 137 species (87.26%) show a "one-to-one" match between a named species and the BIN assignment. Fourteen species (8.92%) show deep conspecific lineages with no apparent morphological differentiation. Only two species pairs shared the same BIN making their identification with DNA barcodes alone uncertain. Therefore, our DNA barcoding reference library allows reliable identification by non-experts for the vast majority of wild bee species in the Loire Valley.

Evaluating DNA Barcoding for Species Identification and Discovery in European Gracillariid Moths

Gracillariidae is the most species-rich leaf-mining moth family with over 2,000 described species worldwide. In Europe, there are 263 valid named species recognized, many of which are difficult to identify using morphology only. Here we explore the use of DNA barcodes as a tool for identification and species discovery in European gracillariids. We present a barcode library including 6,791 COI sequences representing 242 of the 263 (92%) resident species. Our results indicate high congruence between morphology and barcodes with 91.3% (221/242) of European species forming monophyletic clades that can be identified accurately using barcodes alone. The remaining 8.7% represent cases of non-monophyly making their identification uncertain using barcodes. Species discrimination based on the Barcode Index Number system (BIN) was successful for 93% of species with 7% of species sharing BINs. We discovered as many as 21 undescribed candidate species, of which six were confirmed from an integrative approach; the other 15 require additional material and study to confirm preliminary evidence. Most of these new candidate species are found in mountainous regions of Mediterranean countries, the South-Eastern Alps and the Balkans, with nine candidate species found only on islands. In addition, 13 species were classified as deep conspecific lineages, comprising a total of 27 BINs with no intraspecific morphological differences found, and no known ecological differentiation. Double-digest restriction-site associated DNA sequencing (ddRAD) analysis showed strong mitonuclear discrepancy in four out of five species studied. This discordance is not explained by Wolbachia-mediated genetic sweeps. Finally, 26 species were classified as “unassessed species splits” containing 71 BINs and some involving geographical isolation or ecological specialization that will require further study to test whether they represent new cryptic species.

Fifty Years of Research on European Mink Mustela lutreola L., 1761 Genetics: Where Are We Now in Studies on One of the Most Endangered Mammals?

The purpose of this review is to present the current state of knowledge about the genetics of European mink Mustela lutreola L., 1761, which is one of the most endangered mammalian species in the world. This article provides a comprehensive description of the studies undertaken over the last 50 years in terms of cytogenetics, molecular genetics, genomics (including mitogenomics), population genetics of wild populations and captive stocks, phylogenetics, phylogeography, and applied genetics (including identification by genetic methods, molecular ecology, and conservation genetics). An extensive and up-to-date review and critical analysis of the available specialist literature on the topic is provided, with special reference to conservation genetics. Unresolved issues are also described, such as the standard karyotype, systematic position, and whole-genome sequencing, and hotly debated issues are addressed, like the origin of the Southwestern population of the European mink and management approaches of the most distinct populations of the species. Finally, the most urgent directions of future research, based on the research questions arising from completed studies and the implementation of conservation measures to save and restore M. lutreola populations, are outlined. The importance of the popularization of research topics related to European mink genetics among scientists is highlighted.

Gene flow and species delimitation in fishes of Western North America: Flannelmouth (Catostomus latipinnis) and Bluehead sucker (C. Pantosteus discobolus)

The delimitation of species boundaries, particularly those obscured by reticulation, is a critical step in contemporary biodiversity assessment. It is especially relevant for conservation and management of indigenous fishes in western North America, represented herein by two species with dissimilar life histories codistributed in the highly modified Colorado River (i.e., flannelmouth sucker, Catostomus latipinnis; bluehead sucker, C. (Pantosteus) discobolus). To quantify phylogenomic patterns and examine proposed taxonomic revisions, we first employed double-digest restriction site-associated DNA sequencing (ddRAD), yielding 39,755 unlinked SNPs across 139 samples. These were subsequently evaluated with multiple analytical approaches and by contrasting life history data. Three phylogenetic methods and a Bayesian assignment test highlighted similar phylogenomic patterns in each, but with considerable difference in presumed times of divergence. Three lineages were detected in bluehead sucker, supporting elevation of C. (P.) virescens to species status and recognizing C. (P.) discobolus yarrowi (Zuni bluehead sucker) as a discrete entity. Admixture in the latter necessitated a reevaluation of its contemporary and historic distributions, underscoring how biodiversity identification can be confounded by complex evolutionary histories. In addition, we defined three separate flannelmouth sucker lineages as ESUs (evolutionarily significant units), given limited phenotypic and genetic differentiation, contemporary isolation, and lack of concordance (per the genealogical concordance component of the phylogenetic species concept). Introgression was diagnosed in both species, with the Little Colorado and Virgin rivers in particular. Our diagnostic methods, and the agreement of our SNPs with previous morphological, enzymatic, and mitochondrial work, allowed us to partition complex evolutionary histories into requisite components, such as isolation versus secondary contact.

HACSim: an R package to estimate intraspecific sample sizes for genetic diversity assessment using haplotype accumulation curves

Assessing levels of standing genetic variation within species requires a robust sampling for the purpose of accurate specimen identification using molecular techniques such as DNA barcoding; however, statistical estimators for what constitutes a robust sample are currently lacking. Moreover, such estimates are needed because most species are currently represented by only one or a few sequences in existing databases, which can safely be assumed to be undersampled. Unfortunately, sample sizes of 5-10 specimens per species typically seen in DNA barcoding studies are often insufficient to adequately capture within-species genetic diversity. Here, we introduce a novel iterative extrapolation simulation algorithm of haplotype accumulation curves, called HACSim (Haplotype Accumulation Curve Simulator) that can be employed to calculate likely sample sizes needed to observe the full range of DNA barcode haplotype variation that exists for a species. Using uniform haplotype and non-uniform haplotype frequency distributions, the notion of sampling sufficiency (the sample size at which sampling accuracy is maximized and above which no new sampling information is likely to be gained) can be gleaned. HACSim can be employed in two primary ways to estimate specimen sample sizes: (1) to simulate haplotype sampling in hypothetical species, and (2) to simulate haplotype sampling in real species mined from public reference sequence databases like the Barcode of Life Data Systems (BOLD) or GenBank for any genomic marker of interest. While our algorithm is globally convergent, runtime is heavily dependent on initial sample sizes and skewness of the corresponding haplotype frequency distribution.

Using DNA barcoding to improve invasive pest identification at U.S. ports-of-entry

Interception of potential invasive species at ports-of-entry is essential for effective biosecurity and biosurveillance programs. However, taxonomic assessment of the immature stages of most arthropods is challenging; characters for identification are often dependent on adult morphology and reproductive structures. This study aims to strengthen the identification of such specimens through DNA barcoding, with a focus on microlepidoptera. A sample of 241 primarily immature microlepidoptera specimens intercepted at U.S. ports-of-entry from 2007 to 2011 were selected for analysis. From this sample, 201 COI-5P sequences were generated and analyzed for concordance between morphology-based and DNA-based identifications. The retrospective analysis of the data over 10 years (2009 to 2019) using the Barcode of Life Data (BOLD) system demonstrates the importance of establishing and growing DNA barcode reference libraries for use in specimen identification. Additionally, analysis of specimen identification using public data (43.3% specimens identified) vs. non-public data (78.6% specimens identified) highlights the need to encourage researchers to make data publicly accessible. DNA barcoding surpassed morphological identification with 42.3% (public) and 66.7% (non-public) of the sampled specimens achieving a species-level identification, compared to 38.3% species-level identification by morphology. Whilst DNA barcoding was not able to identify all specimens in our dataset, its incorporation into border security programs as an adjunct to morphological identification can provide secondary lines of evidence and lower taxonomic resolution in many cases. Furthermore, with increased globalization, database records need to be clearly annotated for suspected specimen origin versus interception location.