Advances in metabarcoding techniques bring us closer to reliable monitoring of the marine benthos

Comparative study of traditional and DNA-based methods for environmental impact assessment: A case study of marine aggregate extraction in the North Sea

Environmental impact assessments of marine aggregate extraction are traditionally conducted based on morphological characteristics of macrobenthos, which is time-consuming, labour-intensive and requires specific taxonomic expert knowledge. Bulk DNA metabarcoding is suggested as a promising alternative. This study compares the traditional morphological and the bulk DNA metabarcoding method to assess the impact of sand extraction activities on three sandbanks in the Belgian North Sea. Substantial differences in the detected species were observed between methods: Abundant and/or large macrobenthos species were detected by both methods, while small species or species with an exoskeleton were usually only detected by the morphological method. Taxa uniquely detected by bulk DNA metabarcoding could be explained by specimens identified at a higher taxonomic level by morphology, or by specimens with very low read numbers, probably representing species missed in the morphological sorting process, DNA traces on the specimens or false positives during PCR amplification efficiency. Despite the difference in detected species, comparable alpha and beta diversity patterns were observed by both methods, indicating that bulk DNA metabarcoding can effectively detect the overall ecological changes associated with sand extraction. We further demonstrate that bulk DNA metabarcoding reduces sample processing both in time (44 % faster) and cost (26 % cheaper) compared to the morphology-based identification. However, biomass quantification remains challenging for bulk DNA metabarcoding since of the ten most abundant genera, only two genera (Echinocardium and Ophelia) showed a significant positive correlation between biomass and read numbers. Additionally, bulk DNA metabarcoding does not provide information on life stages or size of the identified specimens. As such, our results underpin the complementary nature of both methods, wherein DNA-based analyses allow for rapid detection of community changes (as similar patterns in alpha and beta diversity and biotic index were observed), while morphology-based analyses provide additional information on e.g. secondary production (biomass) and size composition. We show how the strengths of both methods can be combined to assess the impact of sand extraction.

Meiofauna at a tropical sandy beach in the SW Atlantic: the influence of seasonality on diversity

Background

Sandy beaches are dynamic environments housing a large diversity of organisms and providing important environmental services. Meiofaunal metazoan are small organisms that play a key role in the sediment. Their diversity, distribution and composition are driven by sedimentary and oceanographic parameters. Understanding the diversity patterns of marine meiofauna is critical in a changing world.

Methods

In this study, we investigate if there is seasonal difference in meiofaunal assemblage composition and diversity along 1 year and if the marine seascapes dynamics (water masses with particular biogeochemical features, characterized by temperature, salinity, absolute dynamic topography, chromophoric dissolved organic material, chlorophyll-a, and normalized fluorescent line height), rainfall, and sediment parameters (total organic matter, carbonate, carbohydrate, protein, lipids, protein-to-carbohydrate, carbohydrate-to-lipids, and biopolymeric carbon) affect significatively meiofaunal diversity at a tropical sandy beach. We tested two hypotheses here: (i) meiofaunal diversity is higher during warmer months and its composition changes significatively among seasons along a year at a tropical sandy beach, and (ii) meiofaunal diversity metrics are significantly explained by marine seascapes characteristics and sediment parameters. We used metabarcoding (V9 hypervariable region from 18S gene) from sediment samples to assess the meiofaunal assemblage composition and diversity (phylogenetic diversity and Shannon's diversity) over a period of 1 year.

Results

Meiofauna was dominated by Crustacea (46% of sequence reads), Annelida (28% of sequence reads) and Nematoda (12% of sequence reads) in periods of the year with high temperatures (>25 °C), high salinity (>31.5 ppt), and calm waters. Our data support our initial hypotheses revealing a higher meiofaunal diversity (phylogenetic and Shannon's Diversity) and different composition during warmer periods of the year. Meiofaunal diversity was driven by a set of multiple variables, including biological variables (biopolymeric carbon) and organic matter quality (protein content, lipid content, and carbohydrate-to-lipid ratio).

Transition of an estuarine benthic meiofauna assemblage 1.7 and 2.8 years after a mining disaster

Background

Estuaries are transitional coastal ecosystems that are threatened by multiple sources of human pollution. In 2015, mining tailings from an upstream dam failure caused massive metal contamination that impacted benthic assemblages on the Brazilian Rio Doce estuary.

Methods

In this study, we investigate and compare meiofaunal assemblages with eDNA metabarcoding 1.7 years (2017) and 2.8 years (2018) after the initial contamination by mine tailings in order to evaluate the continued impact of sediment mine tailing contaminants on the structure of benthic assemblages after the disaster.

Results

The community was dominated by Arthropoda and Nematoda 1.7 yr after the impacts (42 and 29% of meiofaunal sequence reads, respectively) but after 2.8 years Arthropoda (64.8% of meiofaunal sequence reads) and Rotifera (11.8%) were the most common taxa. This continued impact on meiofaunal assemblage revealed a lower phylogenetic diversity (7.8-fold) in 2018, despite overall decrease in metal concentration (Al, Ba, Cr, As, Fe, Zn, Mn, Pb, Cd, Co) in sediments. Our data suggests that differences in benthic assemblages and loss of diversity may be influenced by contaminants in sediments of this estuary, and indicate that broad eDNA assessments are greatly useful to understand the full range of biodiversity changes in dynamic estuarine ecosystems.

First evaluation of ballast water management systems on operational ships for minimizing introductions of nonindigenous zooplankton

Ballast water is a leading pathway for the global introduction of aquatic nonindigenous species. Most international ships are expected to install ballast water management systems (BWMS) by 2024 to treat ballast water before release. This study examines if ballast water discharges managed by BWMS are meeting standards for organisms ≥50 μm in minimum dimension (i.e., <10 organisms per m3; typically zooplankton). Representative samples of ballast water were collected from 29 ships (using 14 different BWMS) arriving to Canada during 2017-2018. Fourteen samples (48 %) had zooplankton concentrations clearly exceeding the standard (ranging from 18 to 3822 organisms per m3). Nonetheless, compared to earlier management strategies, BWMS appear to reduce the frequency of high-risk introduction events. BWMS filter mesh size was an important predictor of zooplankton concentration following treatment. Greater rates of compliance may be achieved as ship crews gain experience with operation and maintenance of BWMS.

Comparing sediment preservation methods for genomic biomonitoring of coastal marine ecosystems

To avoid loss of genetic information in environmental DNA (eDNA) field samples, the preservation of nucleic acids during field sampling is a critical step. In the development of standard operating procedures (SOPs) for eDNA-based compliance monitoring, the effect of different routinely used sediment preservations on biological community structures serving as bioindicators has gone untested. We compared eDNA metabarcoding results of marine bacterial communities from sample aliquots that were treated with a nucleic acid preservation solution (treated samples) and aliquots that were frozen without further treatment (non-treated samples). Sediment samples were obtained from coastal locations subjected to different stressors (aquaculture, urbanization, industry). DNA extraction efficiency, bacterial community profiles, and measures of alpha- and beta-diversity were highly congruent between treated and non-treated samples. As both preservation methods provide the same relevant information to environmental managers and regulators, we recommend the inclusion of both methods into SOPs for biomonitoring in marine coastal environments.

To denoise or to cluster, that is not the question: optimizing pipelines for COI metabarcoding and metaphylogeography

BACKGROUND

The recent blooming of metabarcoding applications to biodiversity studies comes with some relevant methodological debates. One such issue concerns the treatment of reads by denoising or by clustering methods, which have been wrongly presented as alternatives. It has also been suggested that denoised sequence variants should replace clusters as the basic unit of metabarcoding analyses, missing the fact that sequence clusters are a proxy for species-level entities, the basic unit in biodiversity studies. We argue here that methods developed and tested for ribosomal markers have been uncritically applied to highly variable markers such as cytochrome oxidase I (COI) without conceptual or operational (e.g., parameter setting) adjustment. COI has a naturally high intraspecies variability that should be assessed and reported, as it is a source of highly valuable information. We contend that denoising and clustering are not alternatives. Rather, they are complementary and both should be used together in COI metabarcoding pipelines.

RESULTS

Using a COI dataset from benthic marine communities, we compared two denoising procedures (based on the UNOISE3 and the DADA2 algorithms), set suitable parameters for denoising and clustering, and applied these steps in different orders. Our results indicated that the UNOISE3 algorithm preserved a higher intra-cluster variability. We introduce the program DnoisE to implement the UNOISE3 algorithm taking into account the natural variability (measured as entropy) of each codon position in protein-coding genes.  This correction increased the number of sequences retained by 88%. The order of the steps (denoising and clustering) had little influence on the final outcome.

CONCLUSIONS

We highlight the need for combining denoising and clustering, with adequate choice of stringency parameters, in COI metabarcoding. We present a program that uses the coding properties of this marker to improve the denoising step. We recommend researchers to report their results in terms of both denoised sequences (a proxy for haplotypes) and clusters formed (a proxy for species), and to avoid collapsing the sequences of the latter into a single representative. This will allow studies at the cluster (ideally equating species-level diversity) and at the intra-cluster level, and will ease additivity and comparability between studies.

eDNA metabarcoding of small plankton samples to detect fish larvae and their preys from Atlantic and Pacific waters

Zooplankton community inventories are the basis of fisheries management for containing fish larvae and their preys; however, the visual identification of early-stage larvae (the "missing biomass") is difficult and laborious. Here, eDNA metabarcoding was employed to detect zooplankton species of interest for fisheries from open and coastal waters. High-Throughput sequencing (HTS) from environmental samples using small water volumes has been proposed to detect species of interest whose DNA is the most abundant. We analyzed 6-L water samples taken from subtropical and tropical waters using Cytochrome oxidase I (COI) gene as metabarcode. In the open ocean, several commercial fish larvae and invertebrate species important in fish diet were found from metabarcodes and confirmed from individual barcoding. Comparing Atlantic, Mediterranean, Red Sea, and Pacific samples we found a lower taxonomic depth of OTU assignments in samples from tropical waters than in those from temperate ones, suggesting large gaps in reference databases for those areas; thus a higher effort of zooplankton barcoding in tropical oceans is highly recommended. This and similar simplified sampling protocols could be applied in early detection of species important for fisheries.

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.