Metabarcoding approach for nonindigenous species surveillance in marine coastal waters

Primed and ready: nanopore metabarcoding can now recover highly accurate consensus barcodes that are generally indel-free

BACKGROUND

DNA metabarcoding applies high-throughput sequencing approaches to generate numerous DNA barcodes from mixed sample pools for mass species identification and community characterisation. To date, however, most metabarcoding studies employ second-generation sequencing platforms like Illumina, which are limited by short read lengths and longer turnaround times. While third-generation platforms such as the MinION (Oxford Nanopore Technologies) can sequence longer reads and even in real-time, application of these platforms for metabarcoding has remained limited possibly due to the relatively high read error rates as well as the paucity of specialised software for processing such reads.

RESULTS

We show that this is no longer the case by performing nanopore-based, cytochrome c oxidase subunit I (COI) metabarcoding on 34 zooplankton bulk samples, and benchmarking the results against conventional Illumina MiSeq sequencing. Nanopore R10.3 sequencing chemistry and super accurate (SUP) basecalling model reduced raw read error rates to ~ 4%, and consensus calling with amplicon_sorter (without further error correction) generated metabarcodes that were ≤ 1% erroneous. Although Illumina recovered a higher number of molecular operational taxonomic units (MOTUs) than nanopore sequencing (589 vs. 471), we found no significant differences in the zooplankton communities inferred between the sequencing platforms. Importantly, 406 of 444 (91.4%) shared MOTUs between Illumina and nanopore were also found to be free of indel errors, and 85% of the zooplankton richness could be recovered after just 12-15 h of sequencing.

CONCLUSION

Our results demonstrate that nanopore sequencing can generate metabarcodes with Illumina-like accuracy, and we are the first study to show that nanopore metabarcodes are almost always indel-free. We also show that nanopore metabarcoding is viable for characterising species-rich communities rapidly, and that the same ecological conclusions can be obtained regardless of the sequencing platform used. Collectively, our study inspires confidence in nanopore sequencing and paves the way for greater utilisation of nanopore technology in various metabarcoding applications.

Assessing the potential of nematode metabarcoding for benthic monitoring of offshore oil platforms

Environmental DNA metabarcoding is gaining momentum as a time and cost-effective tool for biomonitoring and environmental impact assessment. Yet, its use as a replacement for the conventional marine benthic monitoring based on morphological analysis of macrofauna is still challenging. Here we propose to study the meiofauna, which is much better represented in sediment DNA samples. We focus on nematodes, which are the most numerous and diverse group of meiofauna. Our aim is to assess the potential of nematode metabarcoding to monitor impacts associated with offshore oil platform activities. To achieve this goal, we used nematode-optimized marker (18S V1V2-Nema) and universal eukaryotic marker (18S V9) region to analyse 252 sediment DNA samples collected near three offshore oil platforms in the North Sea. For both markers, we analysed changes in alpha and beta diversity in relation to distance from the platforms and environmental variables. We also defined three impact classes based on selected environmental variables that are associated with oil extraction activities and used random forest classifiers to compare the predictive performance of both datasets. Our results show that alpha- and beta-diversity of nematodes varies with the increasing distance from the platforms. The variables directly related to platform activity, such as Ba and THC, strongly influence the nematode community. The nematode metabarcoding data provide more robust predictive models than eukaryotic data. Furthermore, the nematode community appears more stable in time and space, as illustrated by the overlap of nematode datasets obtained from the same platform three years apart. A significative negative correlation between distance and Shannon diversity also advocates for higher performance of the V1V2-Nema over the V9. Overall, these results suggest that the sensitivity of nematodes is higher compared to the eukaryotic community. Hence, nematode metabarcoding has the potential to become an effective tool for benthic monitoring in marine environment.

DNA metabarcoding reveals spatial and temporal variation of fish eye fluke communities in lake ecosystems

Eye flukes (Diplostomidae) are diverse and abundant trematode parasites that form multi-species communities in fish with negative effects on host fitness and survival. However, the environmental factors and host-related characteristics that determine species diversity, composition, and coexistence in such communities remain poorly understood. Here, we developed a cost-effective cox1 region-specific DNA metabarcoding approach to characterize parasitic diplostomid communities in two common fish species (Eurasian perch and common roach) collected from seven temperate lakes in Estonia. We found considerable inter- and intra-lake, as well as inter-host species, variation in diplostomid communities. Sympatric host species characterization revealed that parasite communities were typically more diverse in roach than perch. Additionally, we detected five positive and two negative diplostomid species associations in roach, whereas only a single negative association was observed in perch. These results indicate that diplostomid communities in temperate lakes are complex and dynamic systems exhibiting both spatial and temporal heterogeneity. They are influenced by various environmental factors and by host-parasite and inter-parasite interactions. We expect that the described methodology facilitates ecological and biodiversity research of diplostomid parasites. It is also adaptable to other parasite groups where it could serve to improve current understanding of diversity, distribution, and interspecies interactions of other understudied taxa.

Aquatic environmental DNA: A review of the macro-organismal biomonitoring revolution

Environmental DNA (eDNA) is the fastest growing biomonitoring tool fuelled by two key features: time efficiency and sensitivity. Technological advancements allow rapid biodiversity detection at both species and community levels with increasing accuracy. Concurrently, there has been a global demand to standardise eDNA methods, but this is only possible with an in-depth overview of the technological advancements and a discussion of the pros and cons of available methods. We therefore conducted a systematic literature review of 407 peer-reviewed papers on aquatic eDNA published between 2012 and 2021. We observed a gradual increase in the annual number of publications from four (2012) to 28 (2018), followed by a rapid growth to 124 publications in 2021. This was mirrored by a tremendous diversification of methods in all aspects of the eDNA workflow. For example, in 2012 only freezing was applied to preserve filter samples, whereas we recorded 12 different preservation methods in the 2021 literature. Despite an ongoing standardisation debate in the eDNA community, the field is seemingly moving fast in the opposite direction and we discuss the reasons and implications. Moreover, by compiling the largest PCR-primer database to date, we provide information on 522 and 141 published species-specific and metabarcoding primers targeting a wide range of aquatic organisms. This works as a user-friendly 'distillation' of primer information that was hitherto scattered across hundreds of papers, but the list also reflects which taxa are commonly studied with eDNA technology in aquatic environments such as fish and amphibians, and reveals that groups such as corals, plankton and algae are under-studied. Efforts to improve sampling and extraction methods, primer specificity and reference databases are crucial to capture these ecologically important taxa in future eDNA biomonitoring surveys. In a rapidly diversifying field, this review synthetises aquatic eDNA procedures and can guide eDNA users towards best practice.

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.

Plankton biodiversity and species co-occurrence based on environmental DNA – a multiple marker study

Metabarcoding in combination with high-throughput sequencing (HTS) allows simultaneous detection of multiple taxa by targeting single or several taxonomically informative gene regions from environmental DNA samples. In this study, a multiple-marker HTS approach was applied to investigate the plankton diversity and seasonal succession in the Baltic Sea from winter to autumn. Four different markers targeting the 16S, 18S, and 28S ribosomal RNA genes were employed, including a marker for more efficient dinoflagellate detection. Typical seasonal changes were observed in phyto- and bacterioplankton communities. In phytoplankton, the appearance patterns of selected common, dominant, or harmful species followed the patterns also confirmed based on 20 years of phytoplankton monitoring data. In the case of zooplankton, both macro- and microzooplankton species were detected. However, no seasonal patterns were detected in their appearance. In total, 15 and 2 new zoo- and phytoplankton species were detected from the Baltic Sea. HTS approach was especially useful for detecting microzooplankton species as well as for investigating the co-occurrence and potential interactions of different taxa. The results of this study further exemplify the efficiency of metabarcoding for biodiversity monitoring and the advantage of employing multiple markers through the detection of species not identifiable based on a single marker survey and/or by traditional morphology-based methods.

Towards reproducible metabarcoding data: Lessons from an international cross-laboratory experiment

Advances in high-throughput sequencing (HTS) are revolutionizing monitoring in marine environments by enabling rapid, accurate and holistic detection of species within complex biological samples. Research institutions worldwide increasingly employ HTS methods for biodiversity assessments. However, variance in laboratory procedures, analytical workflows and bioinformatic pipelines impede the transferability and comparability of results across research groups. An international experiment was conducted to assess the consistency of metabarcoding results derived from identical samples and primer sets using varying laboratory procedures. Homogenized biofouling samples collected from four coastal locations (Australia, Canada, New Zealand and the USA) were distributed to 12 independent laboratories. Participants were asked to follow one of two HTS library preparation workflows. While DNA extraction, primers and bioinformatic analyses were purposefully standardized to allow comparison, many other technical variables were allowed to vary among laboratories (amplification protocols, type of instrument used, etc.). Despite substantial variation observed in raw results, the primary signal in the data was consistent, with the samples grouping strongly by geographical origin for all data sets. Simple post hoc data clean-up by removing low-quality samples gave the best improvement in sample classification for nuclear 18S rRNA gene data, with an overall 92.81% correct group attribution. For mitochondrial COI gene data, the best classification result (95.58%) was achieved after correction for contamination errors. The identified critical methodological factors that introduced the greatest variability (preservation buffer, sample defrosting, template concentration, DNA polymerase, PCR enhancer) should be of great assistance in standardizing future biodiversity studies using metabarcoding.

Suitability of a dual COI marker for marine zooplankton DNA metabarcoding

As DNA metabarcoding has become an emerging tool for surveying biodiversity, including its application in legally binding assessments, reliable and efficient barcodes are requested, especially for the highly diverse group of zooplankton. This study focuses on comparing the efficiency of two mitochondrial COI barcodes based on the internal primers mlCOIintF and mlCOIintR utilizing mesozooplankton samples collected in a Mediterranean lagoon. Our results indicate that after a slight adjustment, the mlCOIintR primer performs in combination with jdgLCO1490 (herein) very comparably to the much more widely used primer system mlCOIintF/jgHCO2198+dgHCO2198, in terms of level of taxonomic resolution, species detection and their relative abundance in terms of numbers of reads. As for some groups, like Ctenophora, this barcode is not suitable; a combination of them may be the best option to rely on the Folmer region in its entirety without the risk of losing information for a limited primer match.

Ecosystems monitoring powered by environmental genomics: A review of current strategies with an implementation roadmap

A decade after environmental scientists integrated high-throughput sequencing technologies in their toolbox, the genomics-based monitoring of anthropogenic impacts on the biodiversity and functioning of ecosystems is yet to be implemented by regulatory frameworks. Despite the broadly acknowledged potential of environmental genomics to this end, technical limitations and conceptual issues still stand in the way of its broad application by end-users. In addition, the multiplicity of potential implementation strategies may contribute to a perception that the routine application of this methodology is premature or "in development", hence restraining regulators from binding these tools into legal frameworks. Here, we review recent implementations of environmental genomics-based methods, applied to the biomonitoring of ecosystems. By taking a general overview, without narrowing our perspective to particular habitats or groups of organisms, this paper aims to compare, review and discuss the strengths and limitations of four general implementation strategies of environmental genomics for monitoring: (a) Taxonomy-based analyses focused on identification of known bioindicators or described taxa; (b) De novo bioindicator analyses; (c) Structural community metrics including inferred ecological networks; and (d) Functional community metrics (metagenomics or metatranscriptomics). We emphasise the utility of the three latter strategies to integrate meiofauna and microorganisms that are not traditionally utilised in biomonitoring because of difficult taxonomic identification. Finally, we propose a roadmap for the implementation of environmental genomics into routine monitoring programmes that leverage recent analytical advancements, while pointing out current limitations and future research needs.

Biases in bulk: DNA metabarcoding of marine communities and the methodology involved

With the growing anthropogenic pressure on marine ecosystems, the need for efficient monitoring of biodiversity grows stronger. DNA metabarcoding of bulk samples is increasingly being implemented in ecosystem assessments and is more cost-efficient and less time-consuming than monitoring based on morphology. However, before raw sequences are obtained from bulk samples, a profound number of methodological choices must be made. Here, we critically review the recent methods used for metabarcoding of marine bulk samples (including benthic, plankton and diet samples) and indicate how potential biases can be introduced throughout sampling, preprocessing, DNA extraction, marker and primer selection, PCR amplification and sequencing. From a total of 64 studies evaluated, our recommendations for best practices include to (a) consider DESS as a fixative instead of ethanol, (b) use the DNeasy PowerSoil kit for any samples containing traces of sediment, (c) not limit the marker selection to COI only, but preferably include multiple markers for higher taxonomic resolution, (d) avoid touchdown PCR profiles, (e) use a fixed annealing temperature for each primer pair when comparing across studies or institutes, (f) use a minimum of three PCR replicates, and (g) include both negative and positive controls. Although the implementation of DNA metabarcoding still faces several technical complexities, we foresee wide-ranging advances in the near future, including improved bioinformatics for taxonomic assignment, sequencing of longer fragments and the use of whole-genome information. Despite the bulk of biases involved in metabarcoding of bulk samples, if appropriate controls are included along the data generation process, it is clear that DNA metabarcoding provides a valuable tool in ecosystem assessments.

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.

Comparison of Larval Fish Detections Using Morphology-Based Taxonomy versus High-Throughput Sequencing for Invasive Species Early Detection

When first introduced, invasive species typically evade detection; DNA barcoding coupled with high-throughput sequencing (HTS) may be more sensitive and accurate than morphology-based taxonomy, and thereby improve invasive (or rare) species detection. We quantified the relative error of species detection between morphology-based and HTS-based taxonomic identification of ichthyoplankton collections from the Port of Duluth, Minnesota, an aquatic non-native species introduction 'hot-spot' in the Laurentian Great Lakes. We found HTS-based taxonomy identified 28 species and morphology-based taxonomy 30 species, of which 27 were common to both. Among samples, 76% of family-level taxonomic assignments agreed; however, only 42% of species assignments agreed. Most errors were attributed to morphology-based taxonomy, whereas HTS-based taxonomy error was low. For this study system, for most non-native fishes, the detection probability by randomized survey for larvae was similar to that by a survey that is optimized for non-native species early detection of juveniles and adults. We conclude that classifying taxonomic errors by comparing HTS results against morphology-based taxonomy is an important step toward incorporating HTS-based taxonomy into biodiversity surveys.

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.

Monitoring the ecological status of rivers with diatom eDNA metabarcoding: A comparison of taxonomic markers and analytical approaches for the inference of a molecular diatom index

Recently, several studies demonstrated the usefulness of diatom eDNA metabarcoding as an alternative to assess the ecological quality of rivers and streams. However, the choice of the taxonomic marker as well as the methodology for data analysis differ between these studies, hampering the comparison of their results and effectiveness. The aim of this study was to compare two taxonomic markers commonly used in diatom metabarcoding and three distinct analytical approaches to infer a molecular diatom index. We used the values of classical morphological diatom index as a benchmark for this comparison. We amplified and sequenced both a fragment of the rbcL gene and the V4 region of the 18S rRNA gene for 112 epilithic samples from Swiss and French rivers. We inferred index values using three analytical approaches: by computing it directly from taxonomically assigned sequences, by calibrating de novo the ecovalues of all metabarcodes, and by using a supervised machine learning algorithm to train predictive models. In general, the values of index obtained using the two "taxonomy-free" approaches, encompassing molecular assignment and machine learning, were closer correlated to the values of the morphological index than the values based on taxonomically assigned sequences. The correlations of the three analytical approaches were higher in the case of rbcL compared to the 18S marker, highlighting the importance of the reference database which is more complete for the rbcL marker. Our study confirms the effectiveness of diatom metabarcoding as an operational tool for rivers ecological quality assessment and shows that the analytical approaches by-passing the taxonomic assignments are particularly efficient when reference databases are incomplete.

DNA metabarcoding and morphological analysis - Assessment of zooplankton biodiversity in transitional waters

Zooplankton biodiversity assessment is a crucial element in monitoring marine ecosystem processes and community responses to environmental alterations. In order to evaluate the suitability of metabarcoding for zooplankton biodiversity assessment and biomonitoring as a fast and more cost-effective method, seasonal zooplankton sampling was carried out in the Venice Lagoon and the nearby coastal area (Northern Adriatic Sea). The molecular analysis showed higher taxa richness compared to the classical morphological method (224 vs. 88 taxa), discriminating better the meroplanktonic component, morphologically identified only up to order level. Both methods revealed a similar spatio-temporal distribution pattern and the sequence abundances and individual counts were significantly correlated for various taxonomic groups. These results indicate that DNA metabarcoding is an efficient tool for biodiversity assessments in ecosystems with high spatial and temporal variability, where high sampling effort is required as well as fast alert systems for non-native species (NIS).

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

DNA metabarcoding has the capacity to bolster current biodiversity assessment techniques, including the early detection and monitoring of non-indigenous species (NIS). However, the success of this approach is greatly dependent on the availability, taxonomic coverage and reliability of reference sequences in genetic databases, whose deficiencies can potentially compromise species identifications at the taxonomic assignment step. In this study we assessed lacunae in availability of DNA sequence data from four barcodes (COI, 18S, rbcL and matK) for NIS occurring in European marine and coastal environments. NIS checklists were based on EASIN and AquaNIS databases. The highest coverage was found for COI for Animalia and rbcL for Plantae (up to 63%, for both) and 18S for Chromista (up to 51%), that greatly increased when only high impact species were taken into account (up to 82 to 89%). Results show that different markers have unbalanced representations in genetic databases, implying that the parallel use of more than one marker can act complimentarily and may greatly increase NIS identification rates through DNA-based tools. Furthermore, based on the COI marker, data for approximately 30% of the species had maximum intra-specific distances higher than 3%, suggesting that many NIS may have undescribed or cryptic diversity. Although completing the gaps in reference libraries is essential to make the most of the potential of the DNA-based tools, a careful compilation, verification and annotation of available sequences is fundamental to assemble large curated and reliable reference libraries that provide support for rigorous species identifications.

Environmental DNA and RNA as Records of Human Exposome, Including Biotic/Abiotic Exposures and Its Implications in the Assessment of the Role of Environment in Chronic Diseases

Most of environment-related diseases often result from multiple exposures of abiotic and/or biotic stressors across various life stages. The application of environmental DNA/RNA (eDNA/eRNA) to advance ecological understanding has been very successfully used. However, the eminent extension of eDNA/eRNA-based approaches to estimate human exposure to biotic and/or abiotic environmental stressors to understand the environmental causes of chronic diseases has yet to start. Here, we introduce the potential of eDNA/eRNA for bio-monitoring of human exposome and health effects in the real environmental or occupational settings. This review is the first of its kind to discuss how eDNA/eRNA-based approaches can be applied for assessing the human exposome. eDNA-based exposome assessment is expected to rely on our ability to capture the genome- and epigenome-wide signatures left behind by individuals in the indoor and outdoor physical spaces through shedding, excreting, etc. Records of eDNA/eRNA exposome may reflect the early appearance, persistence, and presence of biotic and/or abiotic-exposure-mediated modifications in these nucleic acid molecules. Functional genome- and epigenome-wide mapping of eDNA offer great promise to help elucidate the human exposome. Assessment of longitudinal exposure to physical, biological, and chemical agents present in the environment through eDNA/eRNA may enable the building of an integrative causal dynamic stochastic model to estimate environmental causes of human health deficits. This model is expected to incorporate key biological pathways and gene networks linking individuals, their geographic locations, and random multi-hits of environmental factors. Development and validation of monitoring of eDNA/eRNA exposome should seriously be considered to introduce into safety and risk assessment and as surrogates of chronic exposure to environmental stressors. Here we highlight that eDNA/eRNA reflecting longitudinal exposure of both biotic and abiotic environmental stressors may serve as records of human exposome and discuss its application as molecular tools for understanding the toxicogenomics basis of environment-related health deficits.

Elucidating Biodiversity Shifts in Ballast Water Tanks during a Cross-Latitudinal Transfer: Complementary Insights from Molecular Analyses

In this study, the evolution of ballast water (BW) assemblages across different trophic levels was characterized over a 21 day cross-latitudinal vessel transit using a combination of molecular methods. Triplicate BW samples were collected every second day and size-fractionated (<2.7, 10, >50 μm). Measurements of adenosine triphosphate (ATP) and metabarcoding of environmental nucleic acid (DNA and RNA) analyses, complemented by microscopy and flow cytometry, were performed on each sample. Measured ATP concentrations exhibited high variance between replicates and a strong negative trend in the large (≥50 μm) fraction over the voyage. In concert with microscopy, the metabarcoding data indicated a die-off of larger metazoans during the first week of study and gradual reductions in dinoflagellates and ochrophytes. The ATP and metabarcoding data signaled persistent or increased cellular activity of heterotrophic bacteria and protists in the BW, which was supported by flow cytometry. The metabarcoding showed the presence of active bacteria in all size fractions, suggesting that the sequential filtration approach does not ensure taxonomical differentiation, which has implications for BW quality assessment. Although our data show that ATP and metabarcoding have potential for indicative BW screening for BW compliance monitoring, further research and technological development is needed to improve representativeness of sampling and deliver the unequivocal response criteria required by the international Ballast Water Management Convention.

Biodiversity assessment of tropical shelf eukaryotic communities via pelagic eDNA metabarcoding

Our understanding of marine communities and their functions in an ecosystem relies on the ability to detect and monitor species distributions and abundances. Currently, the use of environmental DNA (eDNA) metabarcoding is increasingly being applied for the rapid assessment and monitoring of aquatic species. Most eDNA metabarcoding studies have either focussed on the simultaneous identification of a few specific taxa/groups or have been limited in geographical scope. Here, we employed eDNA metabarcoding to compare beta diversity patterns of complex pelagic marine communities in tropical coastal shelf habitats spanning the whole Caribbean Sea. We screened 68 water samples using a universal eukaryotic COI barcode region and detected highly diverse communities, which varied significantly among locations, and proved good descriptors of habitat type and environmental conditions. Less than 15% of eukaryotic taxa were assigned to metazoans, most DNA sequences belonged to a variety of planktonic "protists," with over 50% of taxa unassigned at the phylum level, suggesting that the sampled communities host an astonishing amount of micro-eukaryotic diversity yet undescribed or absent from COI reference databases. Although such a predominance of micro-eukaryotes severely reduces the efficiency of universal COI markers to investigate vertebrate and other metazoans from aqueous eDNA, the study contributes to the advancement of rapid biomonitoring methods and brings us closer to a full inventory of extant marine biodiversity.

Introductions and transfers of species by ballast water in the Adriatic Sea

Following the Editorial addressing the BALMAS project, we open the ballast water management special issue for the Adriatic Sea by providing background information on non-indigenous species and the mechanisms (vectors) of transport. Problems allocating introduction mechanisms for various species with certainty are described; in general, key introduction mechanisms are shipping, with ballast water and biofouling as dominant vectors, and aquaculture activities. The dominant mechanisms for introduction may differ through time, between regions and across species. We highlight ballast water as the focus of an international convention to prevent future introductions, reviewing management options and suggesting future research needs. This assessment is not restricted in application to the Adriatic Sea, but is applicable to other coastal waters. Results of such future work may contribute to the experience building phase planned by the International Maritime Organization for a harmonised implementation of the Ballast Water Management Convention.

Discrimination of Spatial Distribution of Aquatic Organisms in a Coastal Ecosystem Using eDNA

The nonlinearity and complexity of coastal ecosystems often cause difficulties when analyzing spatial and temporal patterns of ecological traits. Environmental DNA (eDNA) monitoring has provided an alternative to overcoming the aforementioned issues associated with classical monitoring. We determined aquatic community taxonomic composition using eDNA based on a meta-barcoding approach that characterizes the general ecological features in the Gwangyang Bay coastal ecosystem. We selected the V9 region of the 18S rDNA gene (18S V9), primarily because of its broad range among eukaryotes. Our results produced more detailed spatial patterns in the study area previously categorized (inner bay, main channel of the bay and outer bay) by Kim et al. (2019). Specifically, the outer bay zone was clearly identified by CCA using genus-level identification of aquatic organisms based on meta-barcoding data. We also found significant relationships between environmental factors. Therefore, eDNA monitoring based on meta-barcoding approach holds great potential as a complemental monitoring tool to identify spatial taxonomic distribution patterns in coastal areas.

Macrobenthic fauna from an upwelling coastal area of Peru (Warm Temperate South-eastern Pacific province -Humboldtian ecoregion)

Background

A total of 162 species and subspecies of marine macroinvertebrates were recorded in the submerged soft and hard substrates around the PERU LNG marine terminal and surrounding area, in the central coast of Peru, 167 km south of Lima, Peru. The collection of specimens was carried out from June 2011 to June 2015 as part of the research studies conducted by the Biodiversity Monitoring and Assessment Program (BMAP) around the marine terminal. The area is part of the Humboldt Current Large Marine Ecosystem, one of the most important upwelling systems in the world.

New information

We identified specimens belonging to 83 families and seven phyla. The list was assembled from the taxonomic identifications made by the BMAP. We identified species and subspecies belonging to phyla Annelida, Arthropoda, Brachiopoda, Cnidaria, Echinodermata, Bryozoa and Mollusca. Phyla Annelida (60 spp.), Arthropoda (47 spp.)and Mollusca (45 spp.) exhibited the largest number of species.

The future of biotic indices in the ecogenomic era: Integrating (e)DNA metabarcoding in biological assessment of aquatic ecosystems

The bioassessment of aquatic ecosystems is currently based on various biotic indices that use the occurrence and/or abundance of selected taxonomic groups to define ecological status. These conventional indices have some limitations, often related to difficulties in morphological identification of bioindicator taxa. Recent development of DNA barcoding and metabarcoding could potentially alleviate some of these limitations, by using DNA sequences instead of morphology to identify organisms and to characterize a given ecosystem. In this paper, we review the structure of conventional biotic indices, and we present the results of pilot metabarcoding studies using environmental DNA to infer biotic indices. We discuss the main advantages and pitfalls of metabarcoding approaches to assess parameters such as richness, abundance, taxonomic composition and species ecological values, to be used for calculation of biotic indices. We present some future developments to fully exploit the potential of metabarcoding data and improve the accuracy and precision of their analysis. We also propose some recommendations for the future integration of DNA metabarcoding to routine biomonitoring programs.

Historical baselines in marine bioinvasions: Implications for policy and management

The human-mediated introduction of marine non-indigenous species is a centuries- if not millennia-old phenomenon, but was only recently acknowledged as a potent driver of change in the sea. We provide a synopsis of key historical milestones for marine bioinvasions, including timelines of (a) discovery and understanding of the invasion process, focusing on transfer mechanisms and outcomes, (b) methodologies used for detection and monitoring, (c) approaches to ecological impacts research, and (d) management and policy responses. Early (until the mid-1900s) marine bioinvasions were given little attention, and in a number of cases actively and routinely facilitated. Beginning in the second half of the 20th century, several conspicuous non-indigenous species outbreaks with strong environmental, economic, and public health impacts raised widespread concerns and initiated shifts in public and scientific perceptions. These high-profile invasions led to policy documents and strategies to reduce the introduction and spread of non-indigenous species, although with significant time lags and limited success and focused on only a subset of transfer mechanisms. Integrated, multi-vector management within an ecosystem-based marine management context is urgently needed to address the complex interactions of natural and human pressures that drive invasions in marine ecosystems.

The Hidden Diversity of Flagellated Protists in Soil

Protists are among the most diverse and abundant eukaryotes in soil. However, gaps between described and sequenced protist morphospecies still present a pending problem when surveying environmental samples for known species using molecular methods. The number of sequences in the molecular PR² database (∼130,000) is limited compared to the species richness expected (>1 million protist species) - limiting the recovery rate. This is important, since high throughput sequencing (HTS) methods are used to find associative patterns between functional traits, taxa and environmental parameters. We performed HTS to survey soil flagellates in 150 grasslands of central Europe, and tested the recovery rate of ten previously isolated and cultivated cercomonad species, among locally found diversity. We recovered sequences for reference soil flagellate species, but also a great number of their phylogenetically evaluated genetic variants, among rare and dominant taxa with presumably own biogeography. This was recorded among dominant (cercozoans, Sandona), rare (apusozoans) and a large hidden diversity of predominantly aquatic protists in soil (choanoflagellates, bicosoecids) often forming novel clades associated with uncultured environmental sequences. Evaluating the reads, instead of the OTUs that individual reads are usually clustered into, we discovered that much of this hidden diversity may be lost due to clustering.

Under the canopy: Community-wide effects of invasive algae in Marine Protected Areas revealed by metabarcoding

We analysed with multigene (18S and COI) metabarcoding the effects of the proliferation of invasive seaweeds on rocky littoral communities in two Spanish Marine Protected Areas. The invasive algae studied were Caulerpa cylindracea, Lophocladia lallemandii and Asparagopsis armata. They are canopy-forming, landscape-dominant seaweeds, and we were interested in their effects on the underlying communities of meiobenthos and macrobenthos, separated in two size fractions through sieving. A new semiquantitative treatment of metabarcoding data is introduced. The results for both markers showed that the presence of the invasive seaweed had a significant effect on the understory communities for Lophocladia lallemandii and Asparagopsis armata but not for Caulerpa cylindracea. Likewise, changes in MOTU richness and diversity with invasion status varied in magnitude and direction depending on the alga considered. Our results showed that metabarcoding allows monitoring of the less conspicuous, but not least important, effects of the presence of dominant invasive seaweeds.

Zooplankton Community Profiling in a Eutrophic Freshwater Ecosystem-Lake Tai Basin by DNA Metabarcoding

Communities of zooplankton, a critical portion of aquatic ecosystems, can be adversely affected by contamination resulting from human activities. Understanding the influence of environmental change on zooplankton communities under field-conditions is hindered by traditional labor-intensive approaches that are prone to taxonomic and enumeration mistakes. Here, metabarcoding of cytochrome c oxidase I (COI) region of mitochondrial DNA was used to characterize the genetic diversity of zooplankton. The species composition of zooplankton communities determined by metabarcoding was consistent with the results based on the traditional morphological approach. The spatial distribution of common species (frequency of occurrence >10 samples) by metabarcoding exhibited good agreement with morphological data. Furthermore, metabarcoding can clearly distinguish the composition of the zooplankton community between lake and river ecosystems. In general, rotifers were more abundant in riverine environments than lakes and reservoirs. Finally, the sequence read number of different taxonomic groups using metabarcoding was positively correlated with the zooplankton biomass inferred by density and body length of zooplankton. Overall, the utility of metabarcoding for taxonomic profiling of zooplankton communities was validated by the morphology-based method on a large ecological scale. Metabarcoding of COI could be a powerful and efficient biomonitoring tool to protect local aquatic ecosystems.

Community barcoding reveals little effect of ocean acidification on the composition of coastal plankton communities: Evidence from a long-term mesocosm study in the Gullmar Fjord, Skagerrak

The acidification of the oceans could potentially alter marine plankton communities with consequences for ecosystem functioning. While several studies have investigated effects of ocean acidification on communities using traditional methods, few have used genetic analyses. Here, we use community barcoding to assess the impact of ocean acidification on the composition of a coastal plankton community in a large scale, in situ, long-term mesocosm experiment. High-throughput sequencing resulted in the identification of a wide range of planktonic taxa (Alveolata, Cryptophyta, Haptophyceae, Fungi, Metazoa, Hydrozoa, Rhizaria, Straminipila, Chlorophyta). Analyses based on predicted operational taxonomical units as well as taxonomical compositions revealed no differences between communities in high CO2 mesocosms (~ 760 μatm) and those exposed to present-day CO2 conditions. Observed shifts in the planktonic community composition were mainly related to seasonal changes in temperature and nutrients. Furthermore, based on our investigations, the elevated CO2 did not affect the intraspecific diversity of the most common mesozooplankter, the calanoid copepod Pseudocalanus acuspes. Nevertheless, accompanying studies found temporary effects attributed to a raise in CO2. Differences in taxa composition between the CO2 treatments could, however, only be observed in a specific period of the experiment. Based on our genetic investigations, no compositional long-term shifts of the plankton communities exposed to elevated CO2 conditions were observed. Thus, we conclude that the compositions of planktonic communities, especially those in coastal areas, remain rather unaffected by increased CO2.

Marine invasions enter the genomic era: three lessons from the past, and the way forward

The expanding scale and increasing rate of marine biological invasions have been documented since the early 20th century. Besides their global ecological and economic impacts, non-indigenous species (NIS) also have attracted much attention as opportunities to explore important eco-evolutionary processes such as rapid adaptation, long-distance dispersal and range expansion, and secondary contacts between divergent evolutionary lineages. In this context, genetic tools have been extensively used in the past 20 years. Three important issues appear to have emerged from such studies. First, the study of NIS has revealed unexpected cryptic diversity in what had previously been assumed homogeneous entities. Second, there has been surprisingly little evidence of strong founder events accompanying marine introductions, a pattern possibly driven by large propagule loads. Third, the evolutionary processes leading to successful invasion have been difficult to ascertain due to faint genetic signals. Here we explore the potential of novel tools associated with high-throughput sequencing (HTS) to address these still pressing issues. Dramatic increase in the number of loci accessible via HTS has the potential to radically increase the power of analyses aimed at species delineation, exploring the population genomic consequences of range expansions, and examining evolutionary processes such as admixture, introgression, and adaptation. Nevertheless, the value of this new wealth of genomic data will ultimately depend on the ability to couple it with expanded "traditional" efforts, including exhaustive sampling of marine populations over large geographic scales, integrated taxonomic analyses, and population level exploration of quantitative trait differentiation through common-garden and other laboratory experiments.

Real-time PCR detection of Didemnum perlucidum (Monniot, 1983) and Didemnum vexillum (Kott, 2002) in an applied routine marine biosecurity context

Prevention and early detection are well recognized as the best strategies for minimizing the risks posed by nonindigenous species (NIS) that have the potential to become marine pests. Central to this is the ability to rapidly and accurately identify the presence of NIS, often from complex environmental samples like biofouling and ballast water. Molecular tools have been increasingly applied to assist with the identification of NIS and can prove particularly useful for taxonomically difficult groups like ascidians. In this study, we have developed real-time PCR assays suited to the specific identification of the ascidians Didemnum perlucidum and Didemnum vexillum. Despite being recognized as important global pests, this is the first time specific molecular detection methods have been developed that can support the early identification and detection of these species from a broad range of environmental sample types. These fast, robust and high-throughput assays represent powerful tools for routine marine biosecurity surveillance, as detection and confirmation of the early presence of species could assist in the timely establishment of emergency responses and control strategies. This study applied the developed assays to confirm the ability to detect Didemnid eDNA in water samples. While previous work has focused on detection of marine larvae from water samples, the development of real-time PCR assays specifically aimed at detecting eDNA of sessile invertebrate species in the marine environment represents a world first and a significant step forwards in applied marine biosecurity surveillance. Demonstrated success in the detection of D. perlucidum eDNA from water samples at sites where it could not be visually identified suggests value in incorporating such assays into biosecurity survey designs targeting Didemnid species.

Detection and characterisation of the biopollutant Xenostrobus securis (Lamarck 1819) Asturian population from DNA Barcoding and eBarcoding

DNA efficiently contributes to detect and understand marine invasions. In 2014 the potential biological pollutant pygmy mussel (Xenostrobus securis) was observed for the first time in the Avilés estuary (Asturias, Bay of Biscay). The goal of this study was to assess the stage of invasion, based on demographic and genetic (DNA Barcoding) characteristics, and to develop a molecular tool for surveying the species in environmental DNA. A total of 130 individuals were analysed for the DNA Barcode cytochrome oxidase I gene in order to determine genetic diversity, population structure, expansion trends, and to inferring introduction hits. Reproduction was evidenced by bimodal size distributions of 1597 mussels. High population genetic variation and genetically distinct clades might suggest multiple introductions from several source populations. Finally, species-specific primers were developed within the DNA barcode for PCR amplification from water samples in order to enabling rapid detection of the species in initial expansion stages.