Genetic monitoring of open ocean biodiversity: An evaluation of DNA metabarcoding for processing continuous plankton recorder samples

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.

Saline mine water influences eukaryote life in shallow groundwater of a tropical sandy stream

Eukaryotic communities in groundwater may be particularly sensitive to disturbance because they are adapted to stable environmental conditions and often have narrow spatial distributions. Traditional methods for characterising these communities, focussing on groundwater-inhabiting macro- and meiofauna (stygofauna), are challenging because of limited taxonomic knowledge and expertise (particularly in less-explored regions), and the time and expense of morphological identification. The primary objective of this study was to evaluate the vulnerability of eukaryote communities in shallow groundwater to mine water discharge containing elevated concentrations of magnesium (Mg) and sulfate (SO4). The study was undertaken in a shallow sand bed aquifer within a wet-dry tropical setting. The aquifer, featuring a saline mine water gradient primarily composed of elevated Mg and SO4, was sampled from piezometers in the creek channel upstream and downstream of the mine water influence during the dry season when only subsurface water flow was present. Groundwater communities were characterised using both morphological assessments of stygofauna from net samples and environmental DNA (eDNA) targeting the 18S rDNA and COI mtDNA genes. eDNA data revealed significant shifts in community composition in response to mine waters, contrasting with findings from traditional morphological composition data. Changes in communities determined using eDNA data were notably associated with concentrations of SO42-, Mg2+ and Na+, and water levels in the piezometers. This underscores the importance of incorporating molecular approaches in impact assessments, as relying solely on traditional stygofauna sampling methods in similar environments may lead to inaccurate conclusions about the responses of the assemblage to studied impacts.

Do red tide events promote an increase in zooplankton biodiversity?

Red tides occurring off the southern coast of Korea impact the marine ecosystem and aquaculture industries. Zooplankton are crucial in the food web, connecting primary producers to higher predators and interact diversely with red tide organisms. This study explores dynamics of the zooplankton community over seven years including three red tide and four non-red tide years in Tongyeong using metabarcoding. In non-red tide years, zooplankton diversity showed typical seasonal patterns, increasing from June to early October. However, during red tide years, diversity remained high, with a shift in species composition-decreased Copepoda and increased Branchiopoda, Echinodermata, Malacostraca, and Annelida. Diversity indices were significantly higher in red tide years across all periods except for the richness in 'after' that showed an insignificant higher value. The differences in zooplankton assemblages across periods were influenced by surface temperatures and the density of the red tide-causing alga Margalefidinium polykrikoides. Eight species emerged as indicator species and showed direct correlations with M. polykrikoides and among them, seven species were indicator species for red tide occurrence years. The ecological characteristics of M. polykrikoides blooms and their recurrent occurrences over several decades suggest that zooplankton may adapt to the toxins and use these blooms as spawning cues. Overall, this study provides comprehensive understanding on changes in zooplankton communities during red tide events, offering novel insights into their ecology.

A circumpolar study of surface zooplankton biodiversity of the Southern Ocean based on eDNA metabarcoding

Under pressure from climate change and fishing, the Southern Ocean ecosystems have been changing. Zooplankton plays a vital role in the food web of the Southern Ocean and is crucial for maintaining ecosystem stability. Investigating the circumpolar-scale species composition and biodiversity of zooplankton is crucial for ensuring ecosystem-based conservation and management of the Southern Ocean in a changing climate. Here, we utilized eDNA metabarcoding to assess the biodiversity of zooplankton in the surface seawater surrounding the Antarctica based on samples collected during two expeditions spanning from 2021 to 2022. The main purpose of this paper is to provide more baseline information about circumpolar zooplankton biodiversity based on the emerging eDNA metabarcoding tool. This comprehensive approach led to the identification of over 300 distinct zooplankton species, forming a diverse community dominated by Jellyfish, Mollusca and Polychaete. Surprisingly, common dominant taxonomic groups such as krill and copepods in the Southern Ocean did not show high relative abundance (reads) in surface seawater. The results of redundancy analysis (RDA) and correlation analysis highlighted that water temperature and chlorophyll a had the most significant impact on the reads and diversity of zooplankton. Notably, the influence of water temperature on zooplankton seemed to be primarily indirect, potentially mediated by its effects on primary productivity. Increasing in primary production might lead to lower zooplankton biodiversity in the Southern Ocean in future. This research underscores the effectiveness of eDNA metabarcoding as a valuable tool for monitoring zooplankton diversity in open seas. Given the ongoing changes in temperature, sea ice extent and their impact on primary production, our findings lay a crucial foundation for using eDNA techniques to establish long-term biodiversity monitoring programs across extensive marine ecosystems in the future.

Are genomic indices effective alternatives to morphology based benthic indices in biomonitoring studies? Perspectives from a major harbour and marine protected area

Ecological assessments are currently being conducted by traditional morpho-taxonomical identification techniques that are time-consuming and often inaccurate. Biomonitoring programs are increasingly being complemented by the more rapid and efficient DNA barcoding approach. We compared the congruency of morpho-taxonomic (AMBI - AZTI's Marine Biotic Index) and genomic (gAMBI) benthic indices in ecological quality status (EcoQS) assignation in Mumbai harbour and Malvan Marine Protected area (MPA). The study, first of its kind to adopt the gAMBI tool in the selected milieu, contributed substantial number of macrobenthic cytochrome c oxidase subunit I gene (COI) sequences that were previously unavailable in the reference library, adding sufficient genetic resources for establishing ecostatus. AMBI and gAMBI values based on presence/absence data related significantly with those derived from abundance data matrices. Taxonomic and genomic indices derived ecostatus corresponded sufficiently well despite minor discrepancies, underscoring the viability of gAMBI as a superior alternative to AMBI in monitoring studies.

Metabarcoding reveals high genetic diversity of harmful algae in the coastal waters of Texas, Gulf of Mexico

Environmental-DNA (eDNA) for metabarcoding is a rapid and effective means to investigate microplankton community composition and species diversity. The objective of this study was to examine the genetic diversity of the phytoplankton community in the Gulf of Mexico, with particular emphasis on harmful algal bloom species. Samples were collected at stations along the coast of Texas in September-October 2017 that were inundated by low salinity waters in the aftermath of Hurricane Harvey. Metabarcodes were generated from the eDNA targeting both the V4 and V8-V9 regions of the 18S rDNA gene. Evaluation of the metabarcodes revealed an unexpectedly high number of harmful algal species during this short period, including five that had not been documented in this region previously. A total of 36 harmful algal species could be differentiated based on V4 and V8-V9 metabarcode markers. Using a phylogenetic approach, the taxonomic resolution of each marker differed and not all species could be differentiated using solely one marker. The V4 region resolved species within some genera (e.g., Heterocapsa), while the V8-V9 marker was necessary to resolve species within other genera (e.g., Chattonella). In other cases, species differentiation within a genus required a combination of both markers (e.g., Prorocentrum, Karenia), or another marker will be needed to resolve all species (e.g., Alexandrium, Dinophysis). We conclude that no single marker can delineate all species, so it is recommended HAB monitoring programs use more than one marker. Overall, the observed diversity of HAB species along the Texas coast using metabarcoding exceeded reports from other parts of the world.

Large mesopelagic fish biomass in the Southern Ocean resolved by acoustic properties

The oceanic mesopelagic zone, 200-1000 m below sea level, holds abundant small fishes that play central roles in ecosystem function. Global mesopelagic fish biomass estimates are increasingly derived using active acoustics, where echosounder-generated signals are emitted, reflected by pelagic organisms and detected by transducers on vessels. Previous studies have interpreted a ubiquitous decline in acoustic reflectance towards the Antarctic continent as a reduction in mesopelagic fish biomass. Here, we use empirical data to estimate species-specific acoustic target strength for the dominant mesopelagic fish of the Scotia Sea in the Southern Ocean. We use these data, alongside estimates of fish relative abundance from net surveys, to interpret signals received in acoustic surveys and calculate mesopelagic biomass of the broader Southern Ocean. We estimate the Southern Ocean mesopelagic fish biomass to be approximately 274 million tonnes if Antarctic krill contribute to the acoustic signal, or 570 million tonnes if mesopelagic fish alone are responsible. These quantities are approximately 1.8 and 3.8 times greater than previous net-based biomass estimates. We also show a peak in fish biomass towards the seasonal ice-edge, corresponding to the preferred feeding grounds of penguins and seals, which may be at risk under future climate change scenarios. Our study provides new insights into the abundance and distributions of ecologically significant mesopelagic fish stocks across the Southern Ocean ecosystem.

Measuring mass: variation among 3,161 species of Canadian Coleoptera and the prospects of a mass registry for all insects

Although biomass values are critical for diverse ecological and evolutionary analyses, they are unavailable for most insect species. Museum specimens have the potential to address this gap, but the variation introduced by sampling and preservation methods is uncertain. This study quantifies species-level variation in the body mass of Canadian Coleoptera based on the analysis of 3,744 specimens representing 3,161 Barcode Index Number (BIN) clusters. Employing the BIN system as a proxy for species allows the inclusion of groups where the taxonomic impediment prevents the assignment of specimens to a Linnaean species. By validating the reproducibility of measurements and evaluating the error introduced by operational complexities such as curatorial practice and the loss of body parts, this study demonstrates that museum specimens can speed the assembly of a mass registry. The results further indicate that congeneric species of Coleoptera generally have limited variation in mass, so a genus-level identification allows prediction of the body mass of species that have not been weighed or measured. Building on the present results, the construction of a mass registry for all insects is feasible.

Pragmatic Applications and Universality of DNA Barcoding for Substantial Organisms at Species Level: A Review to Explore a Way Forward

DNA barcodes are regarded as hereditary succession codes that serve as a recognition marker to address several queries relating to the identification, classification, community ecology, and evolution of certain functional traits in organisms. The mitochondrial cytochrome c oxidase 1 (CO1) gene as a DNA barcode is highly efficient for discriminating vertebrate and invertebrate animal species. Similarly, different specific markers are used for other organisms, including ribulose bisphosphate carboxylase (rbcL), maturase kinase (matK), transfer RNA-H and photosystem II D1-ApbsArabidopsis thaliana (trnH-psbA), and internal transcribed spacer (ITS) for plant species; 16S ribosomal RNA (16S rRNA), elongation factor Tu gene (Tuf gene), and chaperonin for bacterial strains; and nuclear ITS for fungal strains. Nevertheless, the taxon coverage of reference sequences is far from complete for genus or species-level identification. Applying the next-generation sequencing approach to the parallel acquisition of DNA barcode sequences could greatly expand the potential for library preparation or accurate identification in biodiversity research. Overall, this review articulates on the DNA barcoding technology as applied to different organisms, its universality, applicability, and innovative approach to handling DNA-based species identification.

Towards harmonization of DNA metabarcoding for monitoring marine macrobenthos: the effect of technical replicates and pooled DNA extractions on species detection

DNA-based monitoring methods are potentially faster and cheaper compared to traditional morphological benthic identification. DNA metabarcoding involves various methodological choices which can introduce bias leading to a different outcome in biodiversity patterns. Therefore, it is important to harmonize DNA metabarcoding protocols to allow comparison across studies and this requires a good understanding of the effect of methodological choices on diversity estimates. This study investigated the impact of DNA and PCR replicates on the detection of macrobenthos species in locations with high, medium and low diversity. Our results show that two to three DNA replicates were needed in locations with a high and medium diversity to detect at least 80% of the species found in the six DNA replicates, while three to four replicates were needed in the location with low diversity. In contrast to general belief, larger body size or higher abundance of the species in a sample did not increase its detection prevalence among DNA replicates. However, rare species were less consistently detected across all DNA replicates of the location with high diversity compared to locations with less diversity. Our results further show that pooling of DNA replicates did not significantly alter diversity patterns, although a small number of rare species was lost. Finally, our results confirm high variation in species detection between PCR replicates, especially for the detection of rare species. These results contribute to create reliable, time and cost efficient metabarcoding protocols for the characterization of macrobenthos.

Zooplankton diversity monitoring strategy for the urban coastal region using metabarcoding analysis

Marine ecosystems in urban coastal areas are exposed to many risks due to human activity. Thus, long-term and continuous monitoring of zooplankton diversity is necessary. High-throughput DNA metabarcoding has gained recognition as an efficient and highly sensitive approach to accurately describing the species diversity of marine zooplankton assemblages. In this study, we collected 30 zooplankton samples at about 2-week intervals for 1 year. Zooplankton diversity showing a typical four season pattern. Of the “total” and “common” zooplankton, we assigned 267 and 64 taxa. The cluster structure and seasonal diversity pattern were rough when only the “common” zooplankton was used. Our study examined how to maximize the benefits of metabarcoding for monitoring zooplankton diversity in urban coastal areas. The results suggest that to take full advantage of metabarcoding when monitoring a zooplankton community, it is necessary to carefully investigate potential ecosystem threats (non-indigenous species) through sufficient curation rather than disregarding low-abundance operational taxonomic units.

Environmental DNA metabarcoding for monitoring metazoan biodiversity in Antarctic nearshore ecosystems

Antarctic benthic ecosystems support high biodiversity but their characterization is limited to a few well-studied areas, due to the extreme environment and remoteness making access and sampling difficult. Our aim was to compare water and sediment as sources of environmental DNA (eDNA) to better characterise Antarctic benthic communities and further develop practical approaches for DNA-based biodiversity assessment in remote environments. We used a cytochrome c oxidase subunit I (COI) metabarcoding approach to characterise metazoan communities in 26 nearshore sites across 12 locations in the Vestfold Hills (East Antarctica) based on DNA extracted from either sediment cores or filtered seawater. We detected a total of 99 metazoan species from 12 phyla across 26 sites, with similar numbers of species detected in sediment and water eDNA samples. However, significantly different communities were detected in the two sample types at sites where both were collected (i.e., where paired samples were available). For example, nematodes and echinoderms were more likely to be detected exclusively in sediment and water eDNA samples, respectively. eDNA from water and sediment core samples are complementary sample types, with epifauna more likely to be detected in water column samples and infauna in sediment. More reference DNA sequences are needed for infauna/meiofauna to increase the proportion of sequences and number of taxa that can be identified. Developing a better understanding of the temporal and spatial dynamics of eDNA at low temperatures would also aid interpretation of eDNA signals from polar environments. Our results provide a preliminary scan of benthic metazoan communities in the Vestfold Hills, with additional markers required to provide a comprehensive biodiversity survey. However, our study demonstrates the choice of sample type for eDNA studies of benthic ecosystems (sediment, water or both) needs to be carefully considered in light of the research or monitoring question of interest.

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.

Pooling size sorted Malaise trap fractions to maximize taxon recovery with metabarcoding

Background

Small and rare specimens can remain undetected when metabarcoding is applied on bulk samples with a high specimen size heterogeneity. This is especially critical for Malaise trap samples, where most of the biodiversity is contributed by small taxa with low biomass. The separation of samples in different size fractions for downstream analysis is one possibility to increase detection of small and rare taxa. However, experiments systematically testing different size sorting approaches and subsequent proportional pooling of fractions are lacking, but would provide important information for the optimization of metabarcoding protocols. We set out to find a size sorting strategy for Malaise trap samples that maximizes taxonomic recovery but remains scalable and time efficient.

Methods

Three Malaise trap samples were sorted into four size classes using dry sieving. Each fraction was homogenized and lysed. The corresponding lysates were pooled to simulate unsorted samples. Pooling was additionally conducted in equal proportions and in four different proportions enriching the small size fraction of samples. DNA from the individual size classes as well as the pooled fractions was extracted and metabarcoded using the FwhF2 and Fol-degen-rev primer set. Additionally, alternative wet sieving strategies were explored.

Results

The small size fractions harboured the highest diversity and were best represented when pooling in favour of small specimens. Metabarcoding of unsorted samples decreases taxon recovery compared to size sorted samples. A size separation into only two fractions (below 4 mm and above) can double taxon recovery compared to not size sorting. However, increasing the sequencing depth 3- to 4-fold can also increase taxon recovery to levels comparable with size sorting, but remains biased towards biomass rich taxa in the sample.

Conclusion

We demonstrate that size fractionation of Malaise trap bulk samples can increase taxon recovery. While results show distinct patterns, the lack of statistical support due to the limited number of samples processed is a limitation. Due to increased speed and lower risk of cross-contamination as well as specimen damage we recommend wet sieving and proportional pooling of the lysates in favour of the small size fraction (80–90% volume). However, for large-scale projects with time constraints, increasing sequencing depth is an alternative solution.

Measuring the absolute abundance of the microbiome by adding yeast containing 16S rRNA gene from a hyperthermophile

High-throughput sequencing (HTS) of 16S rRNA gene amplicons provides compositional information regarding the microbial community, but not the absolute abundance of the bacteria. We aimed to develop a standardized method for quantifying the absolute abundance of bacteria in microbiome studies. To demonstrate the utility of our approach, we quantified the number of bacteria from the compositional data of the fecal and cecal microbiomes. The 16S rRNA gene of a hyperthermophile, Thermus aquaticus, was cloned into Pichia pastoris (yeast) genome, and an equivalent amount of the yeast was added to the stool and cecal samples of mice before DNA extraction. 16S rRNA gene library construction and HTS were performed after DNA extraction. The absolute abundances of bacteria were calculated using T. aquaticus reads. The average relative abundances of T. aquaticus in the five stool and five cecal samples were 0.95% and 0.33%, respectively, indicating that the number of bacteria in a cecum sample is 2.9 times higher than that in a stool sample. The method proposed for quantifying the absolute abundance of the bacterial population in this study is expected to overcome the limitation of showing only compositional data in most microbiome studies.

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.

DNA metabarcoding illustrates biological pollution threats of Red Sea - Dead Sea water conveyance to Dead Sea biodiversity

The Dead Sea has a hypersaline environment where only extremophile species like Archaea, Bacteria, and fungi can survive. The Red Sea-Dead Sea Water Conveyance (RDSC) is constructing a pipeline of 180 km to import water from the Red Sea in the shrinking Dead Sea. Both seas exhibit highly different hydrographic features that determine their biodiversity. Using environmental DNA (eDNA) metabarcoding and amplifying a fragment of the cytochrome oxidase subunit I gene (COI) from water samples, we compared the communities of the Red Sea and the Dead Sea to understand the potential impact of the water conveyance project on biodiversity, following by an identification of potential biopollutants able to enter the Dead Sea for their small size. The results suggest a high likelihood of acquiring harmful algae into the Dead Sea. This study alerts about the real risk of losing the unique Dead Sea biota when the conveyance is actually undertaken.

A Comparison of DNA Metabarcoding and Microscopy Methodologies for the Study of Aquatic Microbial Eukaryotes

The procedures and methodologies employed to study microbial eukaryotic plankton have been thoroughly discussed. Two main schools exist—one insisting on classic microscopy methodologies and the other supporting modern high-throughput sequencing (DNA metabarcoding). However, few studies have attempted to combine both these approaches; most studies implement one method while ignoring the other. This work aims to contribute to this discussion and examine the advantages and disadvantages of each methodology by comparing marine plankton community results from microscopy and DNA metabarcoding. The results obtained by the two methodologies do not vary significantly for Bacillariophyta, although they do for Dinoflagellata and Ciliophora. The lower the taxonomic level, the higher the inconsistency between the two methodologies for all the studied groups. Considering the different characteristics of microscopy-based identification and DNA metabarcoding, this work underlines that each method should be chosen depending on the aims of the study. DNA metabarcoding provides a better estimate of the taxonomic richness of an ecosystem while microscopy provides more accurate quantitative results regarding abundance and biomass. In any case, the combined use of the two methods, if properly standardized, can provide much more reliable and accurate results for the study of marine microbial eukaryotes.

Diversity of dinoflagellate assemblages in coastal temperate and offshore tropical waters of Australia

BACKGROUND

Dinoflagellates are a ubiquitous and ecologically important component of marine phytoplankton communities, with particularly notable species including those associated with harmful algal blooms (HABs) and those that bioluminesce. High-throughput sequencing offers a novel approach compared to traditional microscopy for determining species assemblages and distributions of dinoflagellates, which are poorly known especially in Australian waters.

RESULTS

We assessed the composition of dinoflagellate assemblages in two Australian locations: coastal temperate Port Phillip Bay and offshore tropical waters of Davies Reef (Great Barrier Reef). These locations differ in certain environmental parameters reflecting latitude as well as possible anthropogenic influences. Molecular taxonomic assessment revealed more species than traditional microscopy, and it showed statistically significant differences in dinoflagellate assemblages between locations. Bioluminescent species and known associates of HABs were present at both sites. Dinoflagellates in both areas were mainly represented by the order Gymnodiniales (66%-82% of total sequence reads). In the warm waters of Davies Reef, Gymnodiniales were equally represented by the two superclades, Gymnodiniales sensu stricto (33%) and Gyrodinium (34%). In contrast, in cooler waters of Port Phillip Bay, Gymnodiniales was mainly represented by Gyrodinium (82%). In both locations, bioluminescent dinoflagellates represented up to 0.24% of the total sequence reads, with Protoperidinium the most abundant genus. HAB-related species, mainly represented by Gyrodinium, were more abundant in Port Phillip Bay (up to 47%) than at Davies Reef (28%), potentially reflecting anthropogenic influence from highly populated and industrial areas surrounding the bay. The entire assemblage of dinoflagellates, as well as the subsets of HAB and bioluminescent species, were strongly correlated with water quality parameters (R2 = 0.56-0.92). Significant predictors differed between the subsets: HAB assemblages were explained by salinity, temperature, dissolved oxygen, and total dissolved solids; whereas, bioluminescent assemblages were explained only by salinity and dissolved oxygen, and had greater variability.

CONCLUSION

High-throughput sequencing and genotyping revealed greater diversity of dinoflagellate assemblages than previously known in both subtropical and temperate Australian waters. Significant correlations of assemblage structure with environmental variables suggest the potential for explaining the distribution and composition of both HAB species and bioluminescent species.

The quest for absolute abundance: The use of internal standards for DNA-based community ecology

To characterize microbiomes and other ecological assemblages, ecologists routinely sequence and compare loci that differ among focal taxa. Counts of these sequences convey information regarding the occurrence and relative abundances of taxa, but provide no direct measure of their absolute abundances, due to the technical limitations of the sequencing process. The relative abundances in compositional data are inherently constrained and difficult to interpret. The incorporation of internal standards (ISDs; colloquially referred to as 'spike-ins') into DNA pools can ameliorate the problems posed by relative abundance data and allow absolute abundances to be approximated. Unfortunately, many laboratory and sampling biases cause ISDs to underperform or fail. Here, we discuss how careful deployment of ISDs can avoid these complications and be an integral component of well-designed studies seeking to characterize ecological assemblages via sequencing of DNA.

Capturing open ocean biodiversity: Comparing environmental DNA metabarcoding to the continuous plankton recorder

Environmental DNA (eDNA) metabarcoding is emerging as a novel, objective tool for monitoring marine metazoan biodiversity. Zooplankton biodiversity in the vast open ocean is currently monitored through continuous plankton recorder (CPR) surveys, using ship-based bulk plankton sampling and morphological identification. We assessed whether eDNA metabarcoding (2 L filtered seawater) could capture similar Southern Ocean zooplankton biodiversity as conventional CPR bulk sampling (~1,500 L filtered seawater per CPR sample). We directly compared eDNA metabarcoding with (a) conventional morphological CPR sampling and (b) bulk DNA metabarcoding of CPR collected plankton (two transects for each comparison, 40 and 44 paired samples, respectively). A metazoan-targeted cytochrome c oxidase I (COI) marker was used to characterize species-level diversity. In the 2 L seawater eDNA samples, this marker amplified large amounts of non-metazoan picoplanktonic algae, but eDNA metabarcoding still detected up to 1.6 times more zooplankton species than morphologically analysed bulk CPR samples. COI metabarcoding of bulk DNA samples mostly avoided nonmetazoan amplifications and recovered more zooplankton species than eDNA metabarcoding. However, eDNA metabarcoding detected roughly two thirds of metazoan species and identified similar taxa contributing to community differentiation across the subtropical front separating transects. We observed a diurnal pattern in eDNA data for copepods which perform diel vertical migrations, indicating a surprisingly short temporal eDNA signal. Compared to COI, a eukaryote-targeted 18S ribosomal RNA marker detected a higher proportion, but lower diversity, of metazoans in eDNA. With refinement and standardization of methodology, eDNA metabarcoding could become an efficient tool for monitoring open ocean biodiversity.

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).

Environmental DNA can act as a biodiversity barometer of anthropogenic pressures in coastal ecosystems

Loss of biodiversity from lower to upper trophic levels reduces overall productivity and stability of coastal ecosystems in our oceans, but rarely are these changes documented across both time and space. The characterisation of environmental DNA (eDNA) from sediment and seawater using metabarcoding offers a powerful molecular lens to observe marine biota and provides a series of ‘snapshots’ across a broad spectrum of eukaryotic organisms. Using these next-generation tools and downstream analytical innovations including machine learning sequence assignment algorithms and co-occurrence network analyses, we examined how anthropogenic pressures may have impacted marine biodiversity on subtropical coral reefs in Okinawa, Japan. Based on 18 S ribosomal RNA, but not ITS2 sequence data due to inconsistent amplification for this marker, as well as proxies for anthropogenic disturbance, we show that eukaryotic richness at the family level significantly increases with medium and high levels of disturbance. This change in richness coincides with compositional changes, a decrease in connectedness among taxa, an increase in fragmentation of taxon co-occurrence networks, and a shift in indicator taxa. Taken together, these findings demonstrate the ability of eDNA to act as a barometer of disturbance and provide an exemplar of how biotic networks and coral reefs may be impacted by anthropogenic activities.

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.

Larval assemblages over the abyssal plain in the Pacific are highly diverse and spatially patchy

Abyssal plains are among the most biodiverse yet least explored marine ecosystems on our planet, and they are increasingly threatened by human impacts, including future deep seafloor mining. Recovery of abyssal populations from the impacts of polymetallic nodule mining will be partially determined by the availability and dispersal of pelagic larvae leading to benthic recolonization of disturbed areas of the seafloor. Here we use a tree-of-life (TOL) metabarcoding approach to investigate the species richness, diversity, and spatial variability of the larval assemblage at mesoscales across the abyssal seafloor in two mining-claim areas in the eastern Clarion Clipperton Fracture Zone (CCZ; abyssal Pacific). Our approach revealed a previously unknown taxonomic richness within the meroplankton assemblage, detecting larvae from 12 phyla, 23 Classes, 46 Orders, and 65 Families, including a number of taxa not previously reported at abyssal depths or within the Pacific Ocean. A novel suite of parasitic copepods and worms were sampled, from families that are known to associate with other benthic invertebrates or demersal fishes as hosts. Larval assemblages were patchily distributed at the mesoscale, with little similarity in OTUs detected among deployments even within the same 30 × 30 km study area. Our results provide baseline observations on larval diversity prior to polymetallic nodule mining in this region, and emphasize our overwhelming lack of knowledge regarding larvae of the benthic boundary layer in abyssal plain ecosystems.

SPIKEPIPE: A metagenomic pipeline for the accurate quantification of eukaryotic species occurrences and intraspecific abundance change using DNA barcodes or mitogenomes

The accurate quantification of eukaryotic species abundances from bulk samples remains a key challenge for community ecology and environmental biomonitoring. We resolve this challenge by combining shotgun sequencing, mapping to reference DNA barcodes or to mitogenomes, and three correction factors: (a) a percent-coverage threshold to filter out false positives, (b) an internal-standard DNA spike-in to correct for stochasticity during sequencing, and (c) technical replicates to correct for stochasticity across sequencing runs. The SPIKEPIPE pipeline achieves a strikingly high accuracy of intraspecific abundance estimates (in terms of DNA mass) from samples of known composition (mapping to barcodes R²  = .93, mitogenomes R²  = .95) and a high repeatability across environmental-sample replicates (barcodes R²  = .94, mitogenomes R²  = .93). As proof of concept, we sequence arthropod samples from the High Arctic, systematically collected over 17 years, detecting changes in species richness, species-specific abundances, and phenology. SPIKEPIPE provides cost-efficient and reliable quantification of eukaryotic communities.

New light in the dark - a proposed multidisciplinary framework for studying functional ecology of groundwater fauna

Groundwaters provide the vast majority of unfrozen freshwater resources on the planet, but our knowledge of subsurface ecosystems is surprisingly limited. Stygofauna, or stygobionts -subterranean obligate aquatic animals - provide ecosystem services such as grazing biofilms and maintaining water quality, but we know little about how their ecosystems function. The cryptic nature of groundwaters, together with the high degree of local endemism and stygofaunal site-specific adaptations, represent major obstacles for the field. To overcome these challenges, and integrate biodiversity and ecosystem function, requires a holistic design drawing on classical ecology, taxonomy, molecular ecology and geochemistry. This study presents an approach based on the integration of existing concepts in groundwater ecology with three more novel scientific techniques: compound specific stable isotope analysis (CSIA) of amino acids, radiocarbon analysis (¹⁴C) and DNA analyses of environmental samples, stygofauna and gut contents. The combination of these techniques allows elucidation of aspects of ecosystem function that are often obscured in small invertebrates and cryptic systems. Carbon (δ¹³C) and nitrogen (δ¹⁵N) CSIA provides a linkage between biogeochemical patterns and ecological dynamics. It allows the identification of stygofaunal food web structures and energy flows based on the metabolic pathway of specific amino groups. Concurrently, ¹⁴C provides complementary data on the carbon recycling and incorporation within the stygobiotic trophic webs. Changes in groundwater environmental conditions (e.g. aquifer recharge), and subsequent community adaptations, can be pinpointed via the measurementof the radiocarbon fingerprint of water, sediment and specimens. DNA analyses are a rapidly expanding approach in ecology. eDNA is mainly employed as a biomonitoring tool, while metabarcoding of individuals and/or gut contents provides insight into diet regimes. In all cases, the application of the approaches in combination provides more powerful data than any one alone. By combining quantitative (CSIA and ¹⁴C) and qualitative (eDNA and DNA metabarcoding) approaches via Bayesian Mixing Models (BMM), linkages can be made between community composition, energy and nutrient sources in the system, and trophic function. This suggested multidisciplinary design will contribute to a more thorough comprehension of the biogeochemical and ecological patterns within these undervalued but essential ecosystems.

Molecular analysis indicates high levels of carabid weed seed consumption in cereal fields across Central Europe

Carabid beetles are abundant in temperate agroecosystems and can play a pivotal role as biocontrol agents. While there is good knowledge regarding their effects on invertebrate pests in some systems, comparably little is known on the rate of seed feeding under field conditions. Molecular approaches are ideally suited for investigating carabid feeding interactions; to date, however, they have only been applied to animal prey. We sampled adult carabid beetles in organic cereal fields in three regions along a Central European transect. Regurgitates from populations of the three most common species, Poecilus cupreus, Pseudoophonus rufipes and Pterostichus melanarius, were screened for plant DNA, cereal aphids, collembolans and earthworms. The frequency of carabid individuals positive for plant DNA was high (> 70%) and independent of carabid species, sex, region and the time point of sampling. Detections for non-pest and pest prey were comparably lower, with 21.6% for collembolans, 18.1% for earthworms and 4.2% for aphids, respectively. Despite the prolonged detection period of plant DNA in carabid guts, as compared to animal prey, these first results suggest that weed seeds form an important part of the adult carabid diet. It would also lend support to the hypothesis that seed-feeding carabids are biocontrol agents of weeds, with effects of regulation on the weed seedbank that depend on behavioural and contextual factors including carabid species preferences for weed seed species, their life stage and tillage practices.

Marine environmental DNA biomonitoring reveals seasonal patterns in biodiversity and identifies ecosystem responses to anomalous climatic events

Marine ecosystems are changing rapidly as the oceans warm and become more acidic. The physical factors and the changes to ocean chemistry that they drive can all be measured with great precision. Changes in the biological composition of communities in different ocean regions are far more challenging to measure because most biological monitoring methods focus on a limited taxonomic or size range. Environmental DNA (eDNA) analysis has the potential to solve this problem in biological oceanography, as it is capable of identifying a huge phylogenetic range of organisms to species level. Here we develop and apply a novel multi-gene molecular toolkit to eDNA isolated from bulk plankton samples collected over a five-year period from a single site. This temporal scale and level of detail is unprecedented in eDNA studies. We identified consistent seasonal assemblages of zooplankton species, which demonstrates the ability of our toolkit to audit community composition. We were also able to detect clear departures from the regular seasonal patterns that occurred during an extreme marine heatwave. The integration of eDNA analyses with existing biotic and abiotic surveys delivers a powerful new long-term approach to monitoring the health of our world's oceans in the context of a rapidly changing climate.

Assessing eukaryotic biodiversity in the Florida Keys National Marine Sanctuary through environmental DNA metabarcoding

Environmental DNA (eDNA) is the DNA suspended in the environment (e.g., water column), which includes cells, gametes, and other material derived from but not limited to shedding of tissue, scales, mucus, and fecal matter. Amplifying and sequencing marker genes (i.e., metabarcoding) from eDNA can reveal the wide range of taxa present in an ecosystem through analysis of a single water sample. Metabarcoding of eDNA provides higher resolution data than visual surveys, aiding in assessments of ecosystem health. This study conducted eDNA metabarcoding of two molecular markers (cytochrome c oxidase I (COI) and 18S ribosomal RNA (rRNA) genes) to survey eukaryotic diversity across multiple trophic levels in surface water samples collected at three sites along the coral reef tract within the Florida Keys National Marine Sanctuary (FKNMS) during four research cruises in 2015. The 18S rRNA gene sequences recovered 785 genera while the COI gene sequences recovered 115 genera, with only 33 genera shared between the two datasets, emphasizing the complementarity of these marker genes. Community composition for both genetic markers clustered by month of sample collection, suggesting that temporal variation has a larger effect on biodiversity than spatial variability in the FKNMS surface waters. Sequences from both marker genes were dominated by copepods, but each marker recovered distinct phytoplankton groups, with 18S rRNA gene sequences dominated by dinoflagellates and COI sequences dominated by coccolithophores. Although eDNA samples were collected from surface waters, many benthic species such as sponges, crustaceans, and corals were identified. These results show the utility of eDNA metabarcoding for cataloging biodiversity to establish an ecosystem baseline against which future samples can be compared in order to monitor community changes.

Optimising the detection of marine taxonomic richness using environmental DNA metabarcoding: the effects of filter material, pore size and extraction method

The analysis of environmental DNA (eDNA) using metabarcoding has increased in use as a method for tracking biodiversity of ecosystems. Little is known about eDNA in marine human-modified environments, such as commercial ports, which are key sites to monitor for anthropogenic impacts on coastal ecosystems. To optimise an eDNA metabarcoding protocol in these environments, seawater samples were collected in a commercial port and methodologies for concentrating and purifying eDNA were tested for their effect on eukaryotic DNA yield and subsequent richness of Operational Taxonomic Units (OTUs). Different filter materials [Cellulose Nitrate (CN) and Glass Fibre (GF)], with different pore sizes (0.5 µm, 0.7 µm and 1.2 µm) and three previously published liquid phase extraction methods were tested. The number of eukaryotic OTUs detected differed by a factor of three amongst the method combinations. The combination of CN filters with phenol-chloroform-isoamyl alcohol extractions recovered a higher amount of eukaryotic DNA and OTUs compared to GF filters and the chloroform-isoamyl alcohol extraction method. Pore size was not independent of filter material but did affect the yield of eukaryotic DNA. For the OTUs assigned to a highly successful non-indigenous species, Styelaclava, the two extraction methods with phenol significantly outperformed the extraction method without phenol; other experimental treatments did not contribute significantly to detection. These results highlight that careful consideration of methods is warranted because choice of filter material and extraction method create false negative detections of marine eukaryotic OTUs and underestimate taxonomic richness from environmental samples.