Assessing ecological status with diatoms DNA metabarcoding: Scaling-up on a WFD monitoring network (Mayotte island,France)

Diatoms are excellent ecological indicators of water quality because they are broadly distributed, they show high species diversity and they respond rapidly to human pressures. In Europe, the Water Framework Directive (WFD) gives the legal basis for the use of this indicator for water quality assessment and its management. Several quality indices, like the Specific Polluosensitivity Index (SPI), were developed to assess the ecological quality status of rivers based on diatom communities. It is based on morphological identifications and count of diatom species present in natural biofilms using a microscope. This methodology requires high taxonomic skills and several hours of analysis per sample as 400 individuals must be identified to species level. Since several years, a molecular approach based on DNA metabarcoding combined to High-Throughput Sequencing (HTS) is developed to characterize species assemblages in environmental samples which is potentially faster and cheaper. The ability of this approach to provide reliable diatom inventories has been demonstrated and its application to water quality assessment is currently being improved. Despite optimization of the DNA metabarcoding process with diatoms, few studies had yet extended it at the scale of a freshwater monitoring network and evaluated the reliability of its quality assessment compared to the classical morphological approach.In the present study we applied DNA metabarcoding to the river monitoring network of the tropical Island Mayotte. This island is a French département since 2011 and the WFD has to be applied. This offered the opportunity to scale up the comparison of molecular and morphological approaches and their ability to produce comparable community inventories and water quality assessments. Benthic diatoms were sampled following WFD standards in 45 river sites in 2014 and 2015 (80 samples). All samples were submitted in parallel to the molecular and the morphological approaches. DNA metabarcoding was carried out using Genelute DNA extraction method, rbcL DNA barcode and PGM sequencing, while microscopic counts were carried out for the classical methodology. Diatom community structures in terms of molecular (OTUs) and of morphological (species) were significantly correlated. However, only 13% of the species was shared by both approaches, with qualitative and quantitative variation due to i) the incompleteness of the reference library (82% of morphological species are not represented in the database), ii) limits in taxonomic knowledge and iii) biases in the estimation of relative abundances linked to diatom cell biovolume. However, ecological quality status assessed with the molecular and morphological SPI values were congruent, and little affected by sequencing depth. DNA metabarcoding of diatom communities allowed a reliable estimation of the quality status for most of the rivers at the scale of the full biomonitoring network of Mayotte Island.     

Taxonomy‐free molecular diatom index for high‐throughput eDNA biomonitoring

Current biodiversity assessment and biomonitoring are largely based on the morphological identification of selected bioindicator taxa. Recently, several attempts have been made to use eDNA metabarcoding as an alternative tool. However, until now, most applied metabarcoding studies have been based on the taxonomic assignment of sequences that provides reference to morphospecies ecology. Usually, only a small portion of metabarcoding data can be used due to a limited reference database and a lack of phylogenetic resolution. Here, we investigate the possibility to overcome these limitations using a taxonomy‐free approach that allows the computing of a molecular index directly from eDNA data without any reference to morphotaxonomy. As a case study, we use the benthic diatoms index, commonly used for monitoring the biological quality of rivers and streams. We analysed 87 epilithic samples from Swiss rivers, the ecological status of which was established based on the microscopic identification of diatom species. We compared the diatom index derived from eDNA data obtained with or without taxonomic assignment. Our taxonomy‐free approach yields promising results by providing a correct assessment for 77% of examined sites. The main advantage of this method is that almost 95% of OTUs could be used for index calculation, compared to 35% in the case of the taxonomic assignment approach. Its main limitations are under‐sampling and the need to calibrate the index based on the microscopic assessment of diatoms communities. However, once calibrated, the taxonomy‐free molecular index can be easily standardized and applied in routine biomonitoring, as a complementary tool allowing fast and cost‐effective assessment of the biological quality of watercourses.     

Identifying error and accurately interpreting environmental DNA metabarcoding results: A case study to detect vertebrates at arid zone waterholes

Environmental DNA (eDNA) metabarcoding surveys enable rapid, noninvasive identification of taxa from trace samples with wide‐ranging applications from characterizing local biodiversity to identifying food‐web interactions. However, the technique is prone to error from two major sources: (a) contamination through foreign DNA entering the workflow, and (b) misidentification of DNA within the workflow. Both types of error have the potential to obscure true taxon presence or to increase taxonomic richness by incorrectly identifying taxa as present at sample sites, but multiple error sources can remain unaccounted for in metabarcoding studies. Here, we use data from an eDNA metabarcoding study designed to detect vertebrate species at waterholes in Australia's arid zone to illustrate where and how in the workflow errors can arise, and how to mitigate those errors. We detected the DNA of 36 taxa spanning 34 families, 19 orders and five vertebrate classes in water samples from waterholes, demonstrating the potential for eDNA metabarcoding surveys to provide rapid, noninvasive detection in remote locations, and to widely sample taxonomic diversity from aquatic through to terrestrial taxa. However, we initially identified 152 taxa in the samples, meaning there were many false positive detections. We identified the sources of these errors, allowing us to design a stepwise process to detect and remove error, and provide a template to minimize similar errors that are likely to arise in other metabarcoding studies. Our findings suggest eDNA metabarcoding surveys need to be carefully conducted and screened for errors to ensure their accuracy.     

Effect of marker choice and thermal cycling protocol on zooplankton DNA metabarcoding studies

DNA metabarcoding is a promising approach for rapidly surveying biodiversity and is likely to become an important tool for measuring ecosystem responses to environmental change. Metabarcoding markers need sufficient taxonomic coverage to detect groups of interest, sufficient sequence divergence to resolve species, and will ideally indicate relative abundance of taxa present. We characterized zooplankton assemblages with three different metabarcoding markers (nuclear 18S rDNA, mitochondrial COI, and mitochondrial 16S rDNA) to compare their performance in terms of taxonomic coverage, taxonomic resolution, and correspondence between morphology‐ and DNA‐based identification. COI amplicons sequenced on separate runs showed that operational taxonomic units representing >0.1% of reads per sample were highly reproducible, although slightly more taxa were detected using a lower annealing temperature. Mitochondrial COI and nuclear 18S showed similar taxonomic coverage across zooplankton phyla. However, mitochondrial COI resolved up to threefold more taxa to species compared to 18S. All markers revealed similar patterns of beta‐diversity, although different taxa were identified as the greatest contributors to these patterns for 18S. For calanoid copepod families, all markers displayed a positive relationship between biomass and sequence reads, although the relationship was typically strongest for 18S. The use of COI for metabarcoding has been questioned due to lack of conserved primer‐binding sites. However, our results show the taxonomic coverage and resolution provided by degenerate COI primers, combined with a comparatively well‐developed reference sequence database, make them valuable metabarcoding markers for biodiversity assessment.     

Environmental monitoring through protist next‐generation sequencing metabarcoding: assessing the impact of fish farming on benthic foraminifera communities

The measurement of species diversity represents a powerful tool for assessing the impacts of human activities on marine ecosystems. Traditionally, the impact of fish farming on the coastal environment is evaluated by monitoring the dynamics of macrobenthic infaunal populations. However, taxonomic sorting and morphology‐based identification of the macrobenthos demand highly trained specialists and are extremely time‐consuming and costly, making it unsuitable for large‐scale biomonitoring efforts involving numerous samples. Here, we propose to alleviate this laborious task by developing protist metabarcoding tools based on next‐generation sequencing (NGS) of environmental DNA and RNA extracted from sediment samples. In this study, we analysed the response of benthic foraminiferal communities to the variation of environmental gradients associated with salmon farms in Scotland. We investigated the foraminiferal diversity based on ribosomal minibarcode sequences generated by the Illumina NGS technology. We compared the molecular data with morphospecies counts and with environmental gradients, including distance to cages and redox used as a proxy for sediment oxygenation. Our study revealed high variations between foraminiferal communities collected in the vicinity of fish farms and at distant locations. We found evidence for species richness decrease in impacted sites, especially visible in the RNA data. We also detected some candidate bioindicator foraminiferal species. Based on this proof‐of‐concept study, we conclude that NGS metabarcoding using foraminifera and other protists has potential to become a new tool for surveying the impact of aquaculture and other industrial activities in the marine environment.     

Sorting specimen‐rich invertebrate samples with cost‐effective NGS barcodes: Validating a reverse workflow for specimen processing

Biologists frequently sort specimen‐rich samples to species. This process is daunting when based on morphology, and disadvantageous if performed using molecular methods that destroy vouchers (e.g., metabarcoding). An alternative is barcoding every specimen in a bulk sample and then presorting the specimens using DNA barcodes, thus mitigating downstream morphological work on presorted units. Such a “reverse workflow” is too expensive using Sanger sequencing, but we here demonstrate that is feasible with an next‐generation sequencing (NGS) barcoding pipeline that allows for cost‐effective high‐throughput generation of short specimen‐specific barcodes (313 bp of COI; laboratory cost <$0.50 per specimen) through next‐generation sequencing of tagged amplicons. We applied our approach to a large sample of tropical ants, obtaining barcodes for 3,290 of 4,032 specimens (82%). NGS barcodes and their corresponding specimens were then sorted into molecular operational taxonomic units (mOTUs) based on objective clustering and Automated Barcode Gap Discovery (ABGD). High diversity of 88–90 mOTUs (4% clustering) was found and morphologically validated based on preserved vouchers. The mOTUs were overwhelmingly in agreement with morphospecies (match ratio 0.95 at 4% clustering). Because of lack of coverage in existing barcode databases, only 18 could be accurately identified to named species, but our study yielded new barcodes for 48 species, including 28 that are potentially new to science. With its low cost and technical simplicity, the NGS barcoding pipeline can be implemented by a large range of laboratories. It accelerates invertebrate species discovery, facilitates downstream taxonomic work, helps with building comprehensive barcode databases and yields precise abundance information.     

Environmental metabarcodes for insects: in silico PCR reveals potential for taxonomic bias

Studies of insect assemblages are suited to the simultaneous DNA‐based identification of multiple taxa known as metabarcoding. To obtain accurate estimates of diversity, metabarcoding markers ideally possess appropriate taxonomic coverage to avoid PCR‐amplification bias, as well as sufficient sequence divergence to resolve species. We used in silico PCR to compare the taxonomic coverage and resolution of newly designed insect metabarcodes (targeting 16S) with that of existing markers [16S and cytochrome oxidase c subunit I (COI)] and then compared their efficiency in vitro. Existing metabarcoding primers amplified in silico <75% of insect species with complete mitochondrial genomes available, whereas new primers targeting 16S provided >90% coverage. Furthermore, metabarcodes targeting COI appeared to introduce taxonomic PCR‐amplification bias, typically amplifying a greater percentage of Lepidoptera and Diptera species, while failing to amplify certain orders in silico. To test whether bias predicted in silico was observed in vitro, we created an artificial DNA blend containing equal amounts of DNA from 14 species, representing 11 insect orders and one arachnid. We PCR‐amplified the blend using five primer sets, targeting either COI or 16S, with high‐throughput amplicon sequencing yielding more than 6 million reads. In vitro results typically corresponded to in silico PCR predictions, with newly designed 16S primers detecting 11 insect taxa present, thus providing equivalent or better taxonomic coverage than COI metabarcodes. Our results demonstrate that in silico PCR is a useful tool for predicting taxonomic bias in mixed template PCR and that researchers should be wary of potential bias when selecting metabarcoding markers.     

Evaluating next‐generation sequencing (NGS) methods for routine monitoring of wild bees: Metabarcoding,mitogenomics or NGS barcoding

Implementing cost‐effective monitoring programs for wild bees remains challenging due to the high costs of sampling and specimen identification. To reduce costs, next‐generation sequencing (NGS)‐based methods have lately been suggested as alternatives to morphology‐based identifications. To provide a comprehensive presentation of the advantages and weaknesses of different NGS‐based identification methods, we assessed three of the most promising ones, namely metabarcoding, mitogenomics and NGS barcoding. Using a regular monitoring data set (723 specimens identified using morphology), we found that NGS barcoding performed best for both species presence/absence and abundance data, producing only few false positives (3.4%) and no false negatives. In contrast, the proportion of false positives and false negatives was higher using metabarcoding and mitogenomics. Although strong correlations were found between biomass and read numbers, abundance estimates significantly skewed the communities' composition in these two techniques. NGS barcoding recovered the same ecological patterns as morphology. Ecological conclusions based on metabarcoding and mitogenomics were similar to those based on morphology when using presence/absence data, but different when using abundance data. In terms of workload and cost, we show that metabarcoding and NGS barcoding can compete with morphology, but not mitogenomics which was consistently more expensive. Based on these results, we advocate that NGS barcoding is currently the seemliest NGS method for monitoring of wild bees. Furthermore, this method has the advantage of potentially linking DNA sequences with preserved voucher specimens, which enable morphological re‐examination and will thus produce verifiable records which can be fed into faunistic databases.     

Comparing diversity levels in environmental samples: DNA sequence capture and metabarcoding approaches using 18S and COI genes

Environmental DNA studies targeting multiple taxa using metabarcoding provide remarkable insights into levels of species diversity in any habitat. The main drawbacks are the presence of primer bias and difficulty in identifying rare species. We tested a DNA sequence‐capture method in parallel with the metabarcoding approach to reveal possible advantages of one method over the other. Both approaches were performed using the same eDNA samples and the same 18S and COI regions, followed by high throughput sequencing. Metabarcoded eDNA libraries were PCR amplified with one primer pair from 18S and COI genes. DNA sequence‐capture libraries were enriched with 3,639 baits targeting the same gene regions. We tested amplicon sequence variants (ASVs) and operational taxonomic units (OTUs) in silico approaches for both markers and methods, using for this purpose the metabarcoding data set. ASVs methods uncovered more species for the COI gene, whereas the opposite occurred for the 18S gene, suggesting that clustering reads into OTUs could bias diversity richness especially using 18S with relaxed thresholds. Additionally, metabarcoding and DNA sequence‐capture recovered 80%–90% of the control sample species. DNA sequence‐capture was 8x more expensive, nonetheless it identified 1.5x more species for COI and 13x more genera for 18S than metabarcoding. Both approaches offer reliable results, sharing ca. 40% species and 72% families and retrieve more taxa when nuclear and mitochondrial markers are combined. eDNA metabarcoding is quite well established and low‐cost, whereas DNA‐sequence capture for biodiversity assessment is still in its infancy, is more time‐consuming but provides more taxonomic assignments.     

Toward an ecoregion scale evaluation of eDNA metabarcoding primers: A case study for the freshwater fish biodiversity of the Murray–Darling Basin (Australia)

High‐throughput sequencing of environmental DNA (i.e., eDNA metabarcoding) has become an increasingly popular method for monitoring aquatic biodiversity. At present, such analyses require target‐specific primers to amplify DNA barcodes from co‐occurring species, and this initial amplification can introduce biases. Understanding the performance of different primers is thus recommended prior to undertaking any metabarcoding initiative. While multiple software programs are available to evaluate metabarcoding primers, all programs have their own strengths and weaknesses. Therefore, a robust in silico workflow for the evaluation of metabarcoding primers will benefit from the use of multiple programs. Furthermore, geographic differences in species biodiversity are likely to influence the performance of metabarcoding primers and further complicate the evaluation process. Here, an in silico workflow is presented that can be used to evaluate the performance of metabarcoding primers on an ecoregion scale. This workflow was used to evaluate the performance of published and newly developed eDNA metabarcoding primers for the freshwater fish biodiversity of the Murray–Darling Basin (Australia). To validate the in silico workflow, a subset of the primers, including one newly designed primer pair, were used in metabarcoding analyses of an artificial DNA community and eDNA samples. The results show that the in silico workflow allows for a robust evaluation of metabarcoding primers and can reveal important trade‐offs that need to be considered when selecting the most suitable primer. Additionally, a new primer pair was described and validated that allows for more robust taxonomic assignments and is less influenced by primer biases compared to commonly used fish metabarcoding primers.     

Barcoding diatoms: evaluation of the V4 subregion on the 18S rRNA gene,including new primers and protocols

Diatoms are present in all types of water bodies and their species diversity is influenced greatly by environmental conditions. This means that diatom occurrence and abundances are suitable indicators of water quality. Furthermore, continuous screening of algal biodiversity can provide information about diversity changes in ecosystems. Thus, diatoms represent a desirable group for which to develop an easy to use, quick, efficient, and standardised organism identification tool to serve routine water quality assessments. Because conventional morphological identification of diatoms demands specialised in-depth knowledge, we have established standard laboratory procedures for DNA barcoding in diatoms. We (1) identified a short segment (about 400 bp) of the SSU (18S) rRNA gene which is applicable for the identification of diatom taxa, and (2) elaborated a routine protocol including standard primers for this group of microalgae. To test the universality of the primer binding sites and the discriminatory power of the proposed barcode region, 123 taxa, representing limnic diatom diversity, were included in the study and identified at species level. The effectiveness of the barcode was also scrutinised within a closely related species group, namely the Sellaphora pupula taxon complex and relatives.     

A from‐benchtop‐to‐desktop workflow for validating HTS data and for taxonomic identification in diet metabarcoding studies

The main objective of this work was to develop and validate a robust and reliable “from‐benchtop‐to‐desktop” metabarcoding workflow to investigate the diet of invertebrate‐eaters. We applied our workflow to faecal DNA samples of an invertebrate‐eating fish species. A fragment of the cytochrome c oxidase I (COI) gene was amplified by combining two minibarcoding primer sets to maximize the taxonomic coverage. Amplicons were sequenced by an Illumina MiSeq platform. We developed a filtering approach based on a series of nonarbitrary thresholds established from control samples and from molecular replicates to address the elimination of cross‐contamination, PCR/sequencing errors and mistagging artefacts. This resulted in a conservative and informative metabarcoding data set. We developed a taxonomic assignment procedure that combines different approaches and that allowed the identification of ~75% of invertebrate COI variants to the species level. Moreover, based on the diversity of the variants, we introduced a semiquantitative statistic in our diet study, the minimum number of individuals, which is based on the number of distinct variants in each sample. The metabarcoding approach described in this article may guide future diet studies that aim to produce robust data sets associated with a fine and accurate identification of prey items.     

Vector soup: high‐throughput identification of Neotropical phlebotomine sand flies using metabarcoding

Phlebotomine sand flies are haematophagous dipterans of primary medical importance. They represent the only proven vectors of leishmaniasis worldwide and are involved in the transmission of various other pathogens. Studying the ecology of sand flies is crucial to understand the epidemiology of leishmaniasis and further control this disease. A major limitation in this regard is that traditional morphological‐based methods for sand fly species identifications are time‐consuming and require taxonomic expertise. DNA metabarcoding holds great promise in overcoming this issue by allowing the identification of multiple species from a single bulk sample. Here, we assessed the reliability of a short insect metabarcode located in the mitochondrial 16S rRNA for the identification of Neotropical sand flies, and constructed a reference database for 40 species found in French Guiana. Then, we conducted a metabarcoding experiment on sand flies mixtures of known content and showed that the method allows an accurate identification of specimens in pools. Finally, we applied metabarcoding to field samples caught in a 1‐ha forest plot in French Guiana. Besides providing reliable molecular data for species‐level assignations of phlebotomine sand flies, our study proves the efficiency of metabarcoding based on the mitochondrial 16S rRNA for studying sand fly diversity from bulk samples. The application of this high‐throughput identification procedure to field samples can provide great opportunities for vector monitoring and eco‐epidemiological studies.     

Do we similarly assess diversity with microscopy and high-throughput sequencing? Case of microalgae in lakes

Diatoms are a species-diverse phylum of microalgae often presenting high biomass in aquatic habitats. This makes them excellent ecological indicators in rivers and lakes. They are routinely used to assess ecological quality of rivers and lakes using microscopy, which is time consuming. An alternative is to determine species in samples based on short DNA barcodes and high-throughput sequencing (HTS). Former studies showed that community structure and water quality assessments based on diatoms deliver similar results with both methods. But, none evaluated if diversities were assessed in the same way despite the importance of this ecological metric. Based on littoral benthic samplings carried out in 56 pristine alpine lakes, we compared different diversity indices measured with microscopy and metabarcoding. Each lake was sampled in three different places of its littoral. We showed that α (diversity measured in a single sampling site of a given lake) and ? (total diversity in a lake where three independent samples were considered) diversities obtained with HTS were higher than those obtained with microscopy. This may be explained by the capacity of HTS to detect morphologically cryptic species and to better detect rare taxa. On the other hand, β diversity obtained with HTS was smaller, which may be explained by the capacity of HTS to detect very rare species and free-floating extracellular DNA. Nevertheless, diversity indices obtained with both methodologies were well correlated each other. This study validates the possibility to assess diatom diversity with HTS in a comparable way to the classical microscopic analysis.     

Have the cake and eat it: Optimizing nondestructive DNA metabarcoding of macroinvertebrate samples for freshwater biomonitoring

DNA metabarcoding can contribute to improving cost‐effectiveness and accuracy of biological assessments of aquatic ecosystems, but significant optimization and standardization efforts are still required to mainstream its application into biomonitoring programmes. In assessments based on freshwater macroinvertebrates, a key challenge is that DNA is often extracted from cleaned, sorted and homogenized bulk samples, which is time‐consuming and may be incompatible with sample preservation requirements of regulatory agencies. Here, we optimize and evaluate metabarcoding procedures based on DNA recovered from 96% ethanol used to preserve field samples and thus including potential PCR inhibitors and nontarget organisms. We sampled macroinvertebrates at five sites and subsampled the preservative ethanol at 1 to 14 days thereafter. DNA was extracted using column‐based enzymatic (TISSUE) or mechanic (SOIL) protocols, or with a new magnetic‐based enzymatic protocol (BEAD), and a 313‐bp COI fragment was amplified. Metabarcoding detected at least 200 macroinvertebrate taxa, including most taxa detected through morphology and for which there was a reference barcode. Better results were obtained with BEAD than SOIL or TISSUE, and with subsamples taken 7–14 than 1–7 days after sampling, in terms of DNA concentration and integrity, taxa diversity and matching between metabarcoding and morphology. Most variation in community composition was explained by differences among sites, with small but significant contributions of subsampling day and extraction method, and negligible contributions of extraction and PCR replication. Our methods enhance reliability of preservative ethanol as a potential source of DNA for macroinvertebrate metabarcoding, with a strong potential application in freshwater biomonitoring.     

Taking eDNA underground: Factors affecting eDNA detection of subterranean fauna in groundwater

Stygofauna are aquatic fauna that have evolved to live underground. The impacts of anthropogenic climate change, extraction and pollution on groundwater pose major threats to groundwater health, prompting the need for efficient and reliable means to detect and monitor stygofaunal communities. Conventional survey techniques for these species rely on morphological identification and can be biased, labour-intensive and often indeterminate to lower taxonomic levels. By contrast, environmental DNA (eDNA)-based methods have the potential to dramatically improve on existing stygofaunal survey methods in a large range of habitats and for all life stages, reducing the need for the destructive manual collection of often critically endangered species or for specialized taxonomic expertise. We compared eDNA and haul-net samples collected in 2020 and 2021 from 19 groundwater bores and a cave on Barrow Island, northwest Western Australia, and assessed how sampling factors influenced the quality of eDNA detection of stygofauna. The two detection methods were complementary; eDNA metabarcoding was able to detect soft-bodied taxa and fish often missed by nets, but only detected seven of the nine stygofaunal crustacean orders identified from haul-net specimens. Our results also indicated that eDNA metabarcoding could detect 54%–100% of stygofauna from shallow-water samples and 82%–90% from sediment samples. However, there was significant variation in stygofaunal diversity between sample years and sampling types. The findings of this study demonstrate that haul-net sampling has a tendency to underestimate stygofaunal diversity and that eDNA metabarcoding of groundwater can substantially improve the efficiency of stygofaunal surveys.     

Considerations for metabarcoding‐based port biological baseline surveys aimed at marine nonindigenous species monitoring and risk assessments

1. Monitoring introduction and spread of nonindigenous species via maritime transport and performing risk assessments require port biological baseline surveys. Yet, the comprehensiveness of these surveys is often compromised by the large number of habitats present in a port, the seasonal variability, and the time‐consuming morphological approach used for taxonomic identification. Metabarcoding represents a promising alternative for rapid comprehensive port biological baseline surveys, but its application in this context requires further assessments.
2. We applied metabarcoding (based on barcodes of the cytochrome c oxidase subunit I and of the 18S ribosomal RNA gene) to 192 port samples collected (a) from diverse habitats (water column—including environmental DNA and zooplankton, sediment, and fouling structures), (b) at different sites (from inner to outer estuary), and iii) during the four seasons of the year.
3. By comparing the biodiversity metrics derived from each sample group, we show that each sampling method resulted in a distinct community profile and that environmental DNA alone cannot substitute for organismal sampling, and that, although sampling at different seasons and locations resulted in higher observed biodiversity, operational results can be obtained by sampling selected locations and seasons.
4. By assessing the taxonomic composition of the samples, we show that metabarcoding data allowed the detection of previously recorded nonindigenous species as well as to reveal presence of new ones, even if in low abundance.
5. Synthesis and application. Our comprehensive assessment of metabarcoding for port biological baseline surveys sets the basics for cost‐effective, standardized, and comprehensive monitoring of nonindigenous species and for performing risk assessments in ports. This development will contribute to the implementation of the recently entered into force International Convention for the Control and Management of Ships' Ballast Water and Sediments.

     

Towards routine DNA metabarcoding of macroinvertebrates using bulk samples for freshwater bioassessment: Effects of debris and storage conditions on the recovery of target taxa

1. Macroinvertebrates are commonly sampled for bioassessment of freshwater ecosystems. However, current bioassessment protocols involve laborious sorting of the animals from the debris (sample matrix) and morphological identification, where species level identifications are often difficult. DNA metabarcoding has the potential to improve bioassessment by reducing the time taken to process samples and improve the accuracy and speed of macroinvertebrate species identification.
2. In this study, we evaluated DNA metabarcoding of macroinvertebrate samples, which include macroinvertebrates and the debris collected in the sample nets, to test if bulk, unsorted samples can be used to assess macroinvertebrate diversity. First, we tested if the sample matrix prevented the detection of six target macroinvertebrate taxa when DNA metabarcoding. Second, we tested if sample storage influenced the detection of the same six target macroinvertebrates. We also explored different levels of replication at the sample, sub-sample, and polymerase chain reaction levels and compared the overall macroinvertebrate families detected using DNA metabarcoding to those identified morphologically.
3. We found that the presence of the sample matrix did not interfere with or inhibit the detection of the six target macroinvertebrate taxa. Furthermore, we found that the various sample storage methods did not affect target macroinvertebrate detection. The reliability of detection of the target macroinvertebrates improved as hierarchical levels of replication were combined. We found strong overlap between the detection of overall macroinvertebrate family diversity when comparing DNA metabarcoding to morphological identification.
4. Extracting DNA from the bulk macroinvertebrate samples that included the sample matrix and using this for DNA metabarcoding could improve bioassessment by removing the need for laborious sorting of samples. Furthermore, DNA metabarcoding detection of the six target taxa was not dependent on sample storage of up to 1 year in 95% ethanol, at room temperature or after heating. DNA metabarcoding had the advantage of identifying macroinvertebrate species, but good DNA barcode libraries are needed for widespread species identifications. Further investigation should focus on including multiple samples with different macroinvertebrate composition and densities to refine and standardise bulk sample processing protocols, and on building comprehensive DNA barcode libraries for aquatic macroinvertebrates.