Assessment of species gaps in DNA barcode libraries of non-indigenous species (NIS) occurring in European coastal regions

50%, not great, not terrible: Pan-European gap-analysis shows the real status of the DNA barcode reference libraries in two aquatic invertebrate groups and points the way ahead

The essential key to routine molecular species identification (DNA barcoding/metabarcoding) is the existence of an error-free DNA barcode reference library providing full coverage of all species. Published studies generally state the need to produce more barcodes, and control their quality, but unfortunately, the number of barcoded species is still low. However, to initiate real progress, we need to know where the gaps lie, how big they are and why they persist. Our aims were to draw and understand the current state of knowledge regarding species diversity, distribution, and barcode coverage, and offer solutions for improvement. In this study, we used two groups of aquatic insects, beetles and true bugs. We have compiled and critically evaluated an essentially complete and up-to-date European list, containing 1527 species. The list served as a basis for the barcode gap analyses in the Barcode-of-Life-Data-System (BOLD) conducted in three subsequent years (2020-2022). The overall barcode coverage of the pan-European fauna was around 50 % in both groups. The lowest coverage was in the Mediterranean, the Balkans and South-eastern Europe. The coverage in each country depended significantly on the local diversity, the number of rare, endemic species and the similarity of its fauna to that of the most active barcoding European countries. Gap analyses showed a very small increase in species coverage (<1 % in European aquatic beetles) despite an ~25 % increase in the number of barcodes. Hence, it is clear that future barcoding campaigns must prioritise quality over quantity. To visibly improve reference libraries, we need to increase the involvement of taxonomic experts and focus on targeted studies and underexplored but biodiversity-rich areas.

Compilation, Revision, and Annotation of DNA Barcodes of Marine Invertebrate Non-Indigenous Species (NIS) Occurring in European Coastal Regions

The introduction of non-indigenous species (NIS) is one of the major threats to the integrity of European coastal ecosystems. DNA-based assessments have been increasingly adopted for monitoring NIS. However, the accuracy of DNA-based taxonomic assignments is largely dependent on the completion and reliability of DNA barcode reference libraries. As such, we aimed to compile and audit a DNA barcode reference library for marine invertebrate NIS occurring in Europe. To do so, we compiled a list of NIS using three databases: the European Alien Species Information Network (EASIN), the Information System on Aquatic Non-indigenous and Cryptogenic Species (AquaNIS), and the World Register of Introduced Marine Species (WRiMS). For each species, we retrieved the available cytochrome c oxidase subunit I (COI) mitochondrial gene sequences from the Barcode of Life Data System (BOLD) and used the Barcode, Audit & Grade System (BAGS) to check congruence between morphospecies names and Barcode Index Numbers (BINs). From the 1249 species compiled, approximately 42% had records on BOLD, among which 56% were discordant. We further analyzed these cases to determine the causes of the discordances and attributed additional annotation tags. Of the 622 discordant BINs, after revision, 35% were successfully solved, which increased the number of NIS detected in metabarcoding datasets from 12 to 16. However, a fair number of BINs remained discordant. Reliability of reference barcode records is particularly critical in the case of NIS, where erroneous identification may trigger action or inaction when not required.

DNA sequence and taxonomic gap analyses to quantify the coverage of aquatic cyanobacteria and eukaryotic microalgae in reference databases: Results of a survey in the Alpine region

The taxonomic identification of organisms based on the amplification of specific genetic markers (metabarcoding) implicitly requires adequate discriminatory information and taxonomic coverage of environmental DNA sequences in taxonomic databases. These requirements were quantitatively examined by comparing the determination of cyanobacteria and microalgae obtained by metabarcoding and light microscopy. We used planktic and biofilm samples collected in 37 lakes and 22 rivers across the Alpine region. We focused on two of the most used and best represented genetic markers in the reference databases, namely the 16S rRNA and 18S rRNA genes. A sequence gap analysis using blastn showed that, in the identity range of 99-100%, approximately 30% (plankton) and 60% (biofilm) of the sequences did not find any close counterpart in the reference databases (NCBI GenBank). Similarly, a taxonomic gap analysis showed that approximately 50% of the cyanobacterial and eukaryotic microalgal species identified by light microscopy were not represented in the reference databases. In both cases, the magnitude of the gaps differed between the major taxonomic groups. Even considering the species determined under the microscope and represented in the reference databases, 22% and 26% were still not included in the results obtained by the blastn at percentage levels of identity ≥95% and ≥97%, respectively. The main causes were the absence of matching sequences due to amplification and/or sequencing failure and potential misidentification in the microscopy step. Our results quantitatively demonstrated that in metabarcoding the main obstacles in the classification of 16S rRNA and 18S rRNA sequences and interpretation of high-throughput sequencing biomonitoring data were due to the existence of important gaps in the taxonomic completeness of the reference databases and the short length of reads. The study focused on the Alpine region, but the extent of the gaps could be much greater in other less investigated geographic areas.

Environmental DNA metabarcoding for benthic monitoring: A review of sediment sampling and DNA extraction methods

Environmental DNA (eDNA) metabarcoding (parallel sequencing of DNA/RNA for identification of whole communities within a targeted group) is revolutionizing the field of aquatic biomonitoring. To date, most metabarcoding studies aiming to assess the ecological status of aquatic ecosystems have focused on water eDNA and macroinvertebrate bulk samples. However, the eDNA metabarcoding has also been applied to soft sediment samples, mainly for assessing microbial or meiofaunal biota. Compared to classical methodologies based on manual sorting and morphological identification of benthic taxa, eDNA metabarcoding offers potentially important advantages for assessing the environmental quality of sediments. The methods and protocols utilized for sediment eDNA metabarcoding can vary considerably among studies, and standardization efforts are needed to improve their robustness, comparability and use within regulatory frameworks. Here, we review the available information on eDNA metabarcoding applied to sediment samples, with a focus on sampling, preservation, and DNA extraction steps. We discuss challenges specific to sediment eDNA analysis, including the variety of different sources and states of eDNA and its persistence in the sediment. This paper aims to identify good-practice strategies and facilitate method harmonization for routine use of sediment eDNA in future benthic monitoring.

Water eDNA metabarcoding is effective in detecting non-native species in marinas, but detection errors still hinder its use for passive monitoring

Marinas are high-priority targets for marine non-indigenous species (NIS), where they compose a large portion of the biofouling communities. The practicality of water samples collection makes environmental DNA (eDNA) metabarcoding an interesting tool for routine NIS surveys. Here the effectiveness of water-eDNA-metabarcoding to identify biofouling NIS, in 10 marinas from western France, was examined. Morphological identification of specimens collected in quadrats brought out 18 sessile benthic NIS beneath floating pontoons. Water-eDNA-metabarcoding detected two thirds of them, failing to detect important NIS. However, sampling and bioinformatics filtering steps can be optimized to identify more species. In addition, this method allowed the detection of additional NIS from neighboring micro-habitats. Caution should, however, be taken when reporting putative novel NIS, because of errors in species assignment. This work highlights that water-eDNA-metabarcoding is effective for active (targeted) NIS surveys and could be significantly improved for its further use in marine NIS passive surveys.

Global gap-analysis of amphipod barcode library

In the age of global climate change and biodiversity loss there is an urgent need to provide effective and robust tools for diversity monitoring. One of the promising techniques for species identification is the use of DNA barcoding, that in Metazoa utilizes the so called 'gold-standard' gene of cytochrome c oxidase (COI). However, the success of this method relies on the existence of trustworthy barcode libraries of the species. The Barcode of Life Data System (BOLD) aims to provide barcodes for all existing organisms, and is complemented by the Barcode Index Number (BIN) system serving as a tool for potential species recognition. Here we provide an analysis of all public COI sequences available in BOLD of the diverse and ubiquitous crustacean order Amphipoda, to identify the barcode library gaps and provide recommendations for future barcoding studies. Our gap analysis of 25,702 records has shown that although 3,835 BINs (indicating putative species) were recognised by BOLD, only 10% of known amphipod species are represented by barcodes. We have identified almost equal contribution of both records (sequences) and BINs associated with freshwater and with marine realms. Three quarters of records have a complete species-level identification provided, while BINs have just 50%. Large disproportions between identification levels of BINs coming from freshwaters and the marine environment were observed, with three quarters of the former possessing a species name, and less than 40% for the latter. Moreover, the majority of BINs are represented by a very low number of sequences rendering them unreliable according to the quality control system. The geographical coverage is poor with vast areas of Africa, South America and the open ocean acting as "white gaps". Several, of the most species rich and highly abundant families of Amphipoda (e.g., Phoxocephalidae, Ampeliscidae, Caprellidae), have very poor representation in the BOLD barcode library. As a result of our study we recommend stronger effort in identification of already recognised BINs, prioritising the studies of families that are known to be important and abundant components of particular communities, and targeted sampling programs for taxa coming from geographical regions with the least knowledge.

On the Diversity of Phyllodocida (Annelida: Errantia), with a Focus on Glyceridae, Goniadidae, Nephtyidae, Polynoidae, Sphaerodoridae, Syllidae, and the Holoplanktonic Families

Phyllodocida is a clade of errantiate annelids characterized by having ventral sensory palps, anterior enlarged cirri, axial muscular proboscis, compound chaetae (if present) with a single ligament, and of lacking dorsolateral folds. Members of most families date back to the Carboniferous, although the earliest fossil was dated from the Devonian. Phyllodocida holds 27 well-established and morphologically homogenous clades ranked as families, gathering more than 4600 currently accepted nominal species. Among them, Syllidae and Polynoidae are the most specious polychaete groups. Species of Phyllodocida are mainly found in the marine benthos, although a few inhabit freshwater, terrestrial and planktonic environments, and occur from intertidal to deep waters in all oceans. In this review, we (1) explore the current knowledge on species diversity trends (based on traditional species concept and molecular data), phylogeny, ecology, and geographic distribution for the whole group, (2) try to identify the main knowledge gaps, and (3) focus on selected families: Alciopidae, Goniadidae, Glyceridae, Iospilidae, Lopadorrhynchidae, Polynoidae, Pontodoridae, Nephtyidae, Sphaerodoridae, Syllidae, Tomopteridae, Typhloscolecidae, and Yndolaciidae. The highest species richness is concentrated in European, North American, and Australian continental shelves (reflecting a strong sampling bias). While most data come from shallow coastal and surface environments most world oceans are clearly under-studied. The overall trends indicate that new descriptions are constantly added through time and that less than 10% of the known species have molecular barcode information available.

Integration of DNA-Based Approaches in Aquatic Ecological Assessment Using Benthic Macroinvertebrates

Benthic macroinvertebrates are among the most used biological quality elements for assessing the condition of all types of aquatic ecosystems worldwide (i.e., fresh water, transitional, and marine). Current morphology-based assessments have several limitations that may be circumvented by using DNA-based approaches. Here, we present a comprehensive review of 90 publications on the use of DNA metabarcoding of benthic macroinvertebrates in aquatic ecosystems bioassessments. Metabarcoding of bulk macrozoobenthos has been preferentially used in fresh waters, whereas in marine waters, environmental DNA (eDNA) from sediment and bulk communities from deployed artificial structures has been favored. DNA extraction has been done predominantly through commercial kits, and cytochrome c oxidase subunit I (COI) has been, by far, the most used marker, occasionally combined with others, namely, the 18S rRNA gene. Current limitations include the lack of standardized protocols and broad-coverage primers, the incompleteness of reference libraries, and the inability to reliably extrapolate abundance data. In addition, morphology versus DNA benchmarking of ecological status and biotic indexes are required to allow general worldwide implementation and higher end-user confidence. The increased sensitivity, high throughput, and faster execution of DNA metabarcoding can provide much higher spatial and temporal data resolution on aquatic ecological status, thereby being more responsive to immediate management needs.

Status and prospects of marine NIS detection and monitoring through (e)DNA metabarcoding

In coastal ecosystems, non-indigenous species (NIS) are recognized as a major threat to biodiversity, ecosystem functioning and socio-economic activities. Here we present a systematic review on the use of metabarcoding for NIS surveillance in marine and coastal ecosystems, through the analysis of 42 publications. Metabarcoding has been mainly applied to environmental DNA (eDNA) from water samples, but also to DNA extracted from bulk organismal samples. DNA extraction kits have been widely used and the 18S rRNA and the COI genes the most employed markers, but less than half of the studies targeted more than one marker loci. The Illumina MiSeq platform has been used in >50% of the publications. Current weaknesses include potential occurrence of false negatives due to the primer-biased or faulty DNA amplification and the incompleteness of reference libraries. This is particularly concerning in the case of NIS surveillance, where proficiency in species level detection is critical. Until these weaknesses are resolved, ideally NIS metabarcoding should be supported by complementary approaches, such as morphological analysis or more targeted molecular approaches (e.g. qPCR, ddPCR). Even so, metabarcoding has already proved to be a highly sensitive tool to detect small organisms or undifferentiated life stages across a wide taxonomic range. In addition, it also seems to be very effective in ballast water management and to improve the spatial and temporal sampling frequency of NIS surveillance in marine and coastal ecosystems. Although specific protocols may be required for species-specific NIS detection, for general monitoring it would be vital to settle on a standard protocol able to generate comparable results among surveillance campaigns and regions of the globe, seeking the best approach for detecting the broadest range of species, while minimizing the chances of a false positive or negative detection.