Estimating the barcoding gap in a global dataset of cox1 sequences for Odonata: close, but no cigar

DNA barcoding is currently unreliable for species identification in most crayfishes

DNA barcoding is commonly used for species identification. Despite this, there has not been a comprehensive assessment of the utility of DNA barcoding in crayfishes (Decapoda: Astacidea). Here we examined the extent to which local barcoding gaps (used for species identification) and global barcoding gaps (used for species discovery) exist among crayfishes, and whether global gaps met a previously suggested 10× threshold (mean interspecific difference being 10× larger than mean intra specific difference). We examined barcoding gaps using publicly available mitochondrial COI sequence data from the National Center for Biotechnology Information's nucleotide database. We created two versions of the COI datasets used for downstream analyses: one focused on the number of unique haplotypes (N H) per species, and another that focused on total number of sequences (N S; i.e., including redundant haplotypes) per species. A total of 81 species were included, with 58 species and five genera from the family Cambaridae and 23 species from three genera from the family Parastacidae. Local barcoding gaps were present in only 30 species (20 Cambaridae and 10 Parastacidae species). We detected global barcoding gaps in only four genera (Cambarus, Cherax, Euastacus, and Tenuibranchiurus), which were all below (4.2× to 5.2×) the previously suggested 10× threshold. We propose that a ~5× threshold would be a more appropriate working hypothesis for species discovery. While the N H and N S datasets yielded largely similar results, there were some discrepant inferences. To understand why some species lacked a local barcoding gap, we performed species delimitation analyses for each genus using the N H dataset. These results suggest that current taxonomy in crayfishes may be inadequate for the majority of examined species, and that even species with local barcoding gaps present may be in need of taxonomic revisions. Currently, the utility of DNA barcoding for species identification and discovery in crayfish is quite limited, and caution should be exercised when mitochondrial-based approaches are used in place of taxonomic expertise. Assessment of the evidence for local and global barcoding gaps is important for understanding the reliability of molecular species identification and discovery, but outcomes are dependent on the current state of taxonomy. As this improves (e.g., via resolving species complexes, possibly elevating some subspecies to the species-level status, and redressing specimen misidentifications in natural history and other collections), so too will the utility of DNA barcoding.

Species delimitation with limited sampling: An example from rare trapdoor spider genus Cyclocosmia (Mygalomorphae, Halonoproctidae)

The outcome of species delimitation depends on many factors, including conceptual framework, study design, data availability, methodology employed and subjective decision making. Obtaining sufficient taxon sampling in endangered or rare taxa might be difficult, particularly when non-lethal tissue collection cannot be utilized. The need to avoid overexploitation of the natural populations may thus limit methodological framework available for downstream data analyses and bias the results. We test species boundaries in rare North American trapdoor spider genus Cyclocosmia Ausserer (1871) inhabiting the Southern Coastal Plain biodiversity hotspot with the use of genomic data and two multispecies coalescent model methods. We evaluate the performance of each methodology within a limited sampling framework. To mitigate the risk of species over splitting, common in taxa with highly structured populations, we subsequently implement a species validation step via genealogical diversification index (gdi), which accounts for both genetic isolation and gene flow. We delimited eight geographically restricted lineages within sampled North American Cyclocosmia, suggesting that major river drainages in the region are likely barriers to dispersal. Our results suggest that utilizing BPP in the species discovery step might be a good option for datasets comprising hundreds of loci, but fewer individuals, which may be a common scenario for rare taxa. However, we also show that such results should be validated via gdi, in order to avoid over splitting.

Systematics of the Nososticta salomonis complex (Odonata: Zygoptera: Platycnemididae)

We examined the morphology, colour patterns and genetic relationships of Nososticta populations allied to N. salomonis (Selys) from across Melanesia. Seven species-level taxa are recognised in the N. salomonis 'complex': N. africana (Schmidt), N. boonei sp. nov., N. chrismulleri Theischinger & Richards, N. hedigeri sp. nov., N. salomonis (Selys), N. stueberi sp. nov., and N. tagula sp. nov. All of these species are black damselflies with blue markings, and they differ from all other Nososticta by having: 1) a prominent spike on the male superior appendage, 2) a prominent angular base of the male inferior appendage, and 3) a complex posterior lobe on the female pronotum bearing two pairs of processes in the rough shape of a chair when viewed laterally. A molecular phylogeny based on the DNA barcode fragment of the COI gene plus two nuclear genes indicates that these seven species are closely related, but more extensive sampling of Nososticta species is required to confirm that they form a monophyletic group.

Evaluation of partial 12S rRNA, 16S rRNA, COI and Cytb gene sequence datasets for potential single DNA barcode for hylids (Anura: Hylidae)

We evaluated the extent of intraspecific and interspecific genetic distances and the effectiveness of predefined threshold values using the main genes for estimates of biodiversity and specimens' identification in anurans. Partial sequences of the mitochondrial genes for small (12S) and large (16S) ribosomal subunits, cytochrome c oxidase subunit I (COI) and Cytochrome b (Cytb) of the family Hylidae were downloaded from GenBank and curated for length, coverage, and potential contaminations. We performed analyses for all sequences of each gene and the same species present in these datasets by distance and tree (monophyly)-based evaluations. We also evaluated the ability to identify specimens using these datasets applying "nearest neighbor" (NN), "best close match" (BCM) and "BOLD ID" tests. Genetic distance thresholds were generated by the function 'threshVal' and "localMinima" from SPIDER package and traditional threshold values (1%, 3%, 6% and 10%) were also evaluated. Coding genes, especially COI, had a better identification capacity than non-coding genes on barcoding gap and monophyly analysis and NN, BCM, BOLD ID tests. Considering the multiple factors involved in global DNA barcoding evaluations, we present a critical assessment of the use of these genes for biodiversity estimation and specimens' identification in anurans (e.g. hylids).

The Quality of Sequence Data Affects Biodiversity and Conservation Perspectives in the Neotropical Damselfly Megaloprepus caerulatus

Ideally, the footprint of the evolutionary history of a species is drawn from integrative studies including quantitative and qualitative taxonomy, biogeography, ecology, and molecular genetics. In today’s research, species delimitations and identification of conservation units is often accompanied by a set of—at minimum—two sequence markers appropriate for the systematic level under investigation. Two such studies re-evaluated the species status in the world’s largest Odonata, the Neotropical damselfly Megaloprepus caerulatus. The species status of the genus Megaloprepus has long been debated. Despite applying a highly similar set of sequence markers, the two studies reached different conclusions concerning species status and population genetic relationships. In this study, we took the unique opportunity to compare the two datasets and analyzed the reasons for those incongruences. The two DNA sequence markers used (16S rDNA and CO1) were re-aligned using a strict conservative approach and the analyses used in both studies were repeated. Going step by step back to the first line of data handling, we show that a high number of unresolved characters in the sequence alignments as well as internal gaps are responsible for the different outcomes in terms of species delimitations and population genetic relationships. Overall, this study shows that high quality raw sequence data are an indispensable requirement, not only in odonate research.

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.

A paradigm for the recognition of cryptic trematode species in tropical Indo-west Pacific fishes: the problematic genus Preptetos (Trematoda: Lepocreadiidae)

Molecular data have transformed approaches to trematode taxonomy by providing objective evidence for the delineation of species. However, although the data are objective, the interpretation of these data regarding species boundaries is subjective, especially when different markers conflict. Conserved markers can lead to an underestimation of richness and those used for finer species delineation have the capacity to inflate species recognition, perhaps unrealistically. Here we examine molecular and morphological evidence for species recognition in an especially confusing system, the lepocreadiid genus Preptetos Pritchard, 1960 in acanthuriform fishes of the tropical Indo-west Pacific. We consider species boundaries within this genus based on combined data (ITS2 and 28S rDNA; cox1 mtDNA and morphometrics) for substantial new collections. Delineation of species using only morphological data suggest fewer species than analysis of the sequence data; the latter suggests the presence of potential cryptic species and analysis of different markers suggests the presence of differing numbers of species. We conclude that an integrative interpretation creates the most satisfying taxonomic hypothesis. In the light of the new data, we have chosen and propose a model of trematode species recognition that demands reciprocal monophyly in the most discriminating available molecular marker plus distinction in morphology or host distribution. By invoking these criteria, we distinguish eight species in our new tropical Indo-west Pacific collections. Six of these are new (Preptetos allocaballeroi n. sp., Preptetos paracaballeroi n. sp., Preptetos pearsoni n. sp., Preptetos prudhoei n. sp., Preptetos quandamooka n. sp. and Preptetos zebravaranus n. sp.) and we continue to recognise Preptetos cannoni Barker, Bray & Cribb, 1993 and Preptetos laguncula Bray and Cribb, 1996. Notably; two of the new species, P. allocaballeroi n. sp. and P. paracaballeroi n. sp., are morphologically cryptic relative to each other. Our criteria lead us to recognise, as species, populations with unvarying morphology and similar host relationships but which may have a complex population structure over their range. In our view, this paradigm has the capacity to render tractable the interpretation of the species status of the huge trematode fauna of the tropical Indo-west Pacific.

Coverage and quality of DNA barcode references for Central and Northern European Odonata

Background

Dragonflies and damselflies (Odonata) are important components in biomonitoring due to their amphibiotic lifecycle and specific habitat requirements. They are charismatic and popular insects, but can be challenging to identify despite large size and often distinct coloration, especially the immature stages. DNA-based assessment tools rely on validated DNA barcode reference libraries evaluated in a supraregional context to minimize taxonomic incongruence and identification mismatches.

Methods

This study reports on findings from the analysis of the most comprehensive DNA barcode dataset for Central European Odonata to date, with 103 out of 145 recorded European species included and publicly deposited in the Barcode of Life Data System (BOLD). The complete dataset includes 697 specimens (548 adults, 108 larvae) from 274 localities in 16 countries with a geographic emphasis on Central Europe. We used BOLD to generate sequence divergence metrics and to examine the taxonomic composition of the DNA barcode clusters within the dataset and in comparison with all data on BOLD.

Results

Over 88% of the species included can be readily identified using their DNA barcodes and the reference dataset provided. Considering the complete European dataset, unambiguous identification is hampered in 12 species due to weak mitochondrial differentiation and partial haplotype sharing. However, considering the known species distributions only two groups of five species possibly co-occur, leading to an unambiguous identification of more than 95% of the analysed Odonata via DNA barcoding in real applications. The cases of small interspecific genetic distances and the observed deep intraspecific variation in Cordulia aenea (Linnaeus, 1758) are discussed in detail and the corresponding taxa in the public reference database are highlighted. They should be considered in future applications of DNA barcoding and metabarcoding and represent interesting evolutionary biological questions, which call for in depth analyses of the involved taxa throughout their distribution ranges.

DNA barcoding for identification of anuran species in the central region of South America

The use of COI barcodes for specimen identification and species discovery has been a useful molecular approach for the study of Anura. Here, we establish a comprehensive amphibian barcode reference database in a central area of South America, in particular for specimens collected in Mato Grosso do Sul state (Brazil), and to evaluate the applicability of the COI gene for species-level identification. Both distance- and tree-based methods were applied for assessing species boundaries and the accuracy of specimen identification was evaluated. A total of 204 mitochondrial COI barcode sequences were evaluated from 22 genera and 59 species (19 newly barcoded species). Our results indicate that morphological and molecular identifications converge for most species, however, some species may present cryptic species due to high intraspecific variation, and there is a high efficiency of specimen identification. Thus, we show that COI sequencing can be used to identify anuran species present in this region.

Italian odonates in the Pandora's box: A comprehensive DNA barcoding inventory shows taxonomic warnings at the Holarctic scale

The Odonata are considered among the most endangered freshwater faunal taxa. Their DNA-based monitoring relies on validated reference data sets that are often lacking or do not cover important biogeographical centres of diversification. This study presents the results of a DNA barcoding campaign on Odonata, based on the standard 658-bp 5' end region of the mitochondrial COI gene, involving the collection of 812 specimens (409 of which barcoded) from peninsular Italy and its main islands (328 localities), belonging to all the 88 species (31 Zygoptera and 57 Anisoptera) known from the country. Additional BOLD and GenBank data from Holarctic samples expanded the data set to 1,294 DNA barcodes. A multi-approach species delimitation analysis involving two distance (OT and ABGD) and four tree-based (PTP, MPTP, GMYC and bGMYC) methods was used to explore these data. Of the 88 investigated morphospecies, 75 (85%) unequivocally corresponded to distinct molecular operational units, whereas the remaining ones were classified as 'warnings' (i.e. showing a mismatch between morphospecies assignment and DNA-based species delimitation). These results are in contrast with other DNA barcoding studies on Odonata showing up to 95% of identification success. The species causing warnings were grouped into three categories depending on if they showed low, high or mixed genetic divergence patterns. The analysis of haplotype networks revealed unexpected intraspecific complexity at the Italian, Palearctic and Holarctic scale, possibly indicating the occurrence of cryptic species. Overall, this study provides new insights into the taxonomy of odonates and a valuable basis for future DNA and eDNA-based monitoring studies.

First records raise questions: DNA barcoding of Odonata in the middle of the Mediterranean

We present the results of the first-ever DNA barcoding study of odonates from the Maltese Islands. In total, 10 morphologically identified species were collected during a two-week long expedition in 2018. Eighty cytochrome c oxidase subunit I (COI) barcodes were obtained from the collected specimens. Intra- and interspecific distances ranged from 0.00% to 2.24% and 0.48% to 17.62%, respectively. Successful species identification based on ascribing a single morphological species to a single Barcode Index Number (BIN) was achieved for eight species (80%). In the case of two species, Ischnura genei and Anax parthenope, BINs were shared with other closely related species. The taxonomic status of I. genei is questionable and the phylogenetic relationship between A. imperator/parthenope is not clear. Further studies involving a series of adult specimens collected in a wide spatial range and nuclear markers are necessary to resolve these cases. Therefore, this dataset serves as an initial DNA barcode reference library for Maltese odonates, within a larger project: Aquatic Macroinvertebrates DNA Barcode Library of Malta.

A method to generate multilocus barcodes of pinned insect specimens using MiSeq

For molecular insect identification, amplicon sequencing methods are recommended because they offer a cost-effective approach for targeting small sets of informative genes from multiple samples. In this context, high-throughput multilocus amplicon sequencing has been achieved using the MiSeq Illumina sequencing platform. However, this approach generates short gene fragments of <500 bp, which then have to be overlapped using bioinformatics to achieve longer sequence lengths. This increases the risk of generating chimeric sequences or leads to the formation of incomplete loci. Here, we propose a modified nested amplicon sequencing method for targeting multiple loci from pinned insect specimens using the MiSeq Illumina platform. The modification exists in using a three-step nested PCR approach targeting near full-length loci in the initial PCR and subsequently amplifying short fragments of between 300 and 350 bp for high-throughput sequencing using Illumina chemistry. Using this method, we generated 407 sequences of three loci from 86% of all the specimens sequenced. Out of 103 pinned bee specimens of replicated species, 71% passed the 95% sequence similarity threshold between species replicates. This method worked best for pinned specimens aged between 0 and 5 years, with a limit of 10 years for pinned and 14 years for ethanol-preserved specimens. Hence, our method overcomes some of the challenges of amplicon sequencing using short read next generation sequencing and improves the possibility of creating high-quality multilocus barcodes from insect collections.

A further study on Franciscobasis Machado & Bedê, 2016 (Odonata: Coenagrionidae), a newly described genus from Minas Gerais, Brazil

The genus Franciscobasis Machado & Bedê, 2016 is endemic to the Serra da Canastra National Park in Minas Gerais state, Brazil. Two species of Franciscobasis were described simultaneously with the genus description: F. franciscoi and F. sonia, the latter described only from females. Through morphological and molecular analysis, we investigated if F. sonia may represent the young female of F. franciscoi. Resulting data did not present adequate differences between females to characterize them as different species. Therefore, we suggest that F. sonia is a junior synonym of F. franciscoi, and the female of F. franciscoi goes through a complex ontogenetic color change.

A DNA barcode library for ground beetles of Germany: the genus Amara Bonelli, 1810 (Insecta, Coleoptera, Carabidae)

The genus Amara Bonelli, 1810 is a very speciose and taxonomically difficult genus of the Carabidae. The identification of many of the species is accomplished with considerable difficulty, in particular for females and immature stages. In this study the effectiveness of DNA barcoding, the most popular method for molecular species identification, was examined to discriminate various species of this genus from Central Europe. DNA barcodes from 690 individuals and 47 species were analysed, including sequences from previous studies and more than 350 newly generated DNA barcodes. Our analysis revealed unique BINs for 38 species (81%). Interspecific K2P distances below 2.2% were found for three species pairs and one species trio, including haplotype sharing between Amara alpina/Amara torrida and Amara communis/Amara convexior/Amara makolskii. This study represents another step in generating an extensive reference library of DNA barcodes for carabids, highly valuable bioindicators for characterizing disturbances in various habitats.