Mock samples resolve biases in diversity estimates and quantitative interpretation of zooplankton metabarcoding data

Abstract

Metabarcoding is a rapidly developing tool in marine zooplankton ecology, although most zooplankton surveys continue to rely on visual identification for monitoring purposes. We attempted to resolve some of the biases associated with metabarcoding by sequencing a 313 b.p. fragment of the COI gene in 34 “mock” samples from the North Sea which were pre-sorted to species level, with biomass and abundance estimates obtained for each species and taxonomic group. The samples were preserved either in 97% ethanol or dried for 24 hours in a drying oven at 65° C (the routine way of preserving samples for dry weight measurements). The visual identification yielded a total of 59 unique holoplanktonic and 16 meroplanktonic species/taxa. Metabarcoding identified 86 holoplanktonic and 124 meroplanktonic species/taxa, which included all but 3 of the species identified visually as well as numerous species of hard-to-identify crustaceans, hydrozoan jellyfish and larvae of benthic animals. On a sample-to-sample basis, typically 90–95% of visually registered species were recovered, but the number of false positives was also high. We demonstrate robust correlations of relative sequence abundances to relative biomass for most taxonomic groups and develop conversion factors for different taxa to account for sequencing biases. We then combine the adjusted sequencing data with a single bulk biomass measurement for the entire sample to produce a quantitative parameter akin to species biomass. When examined with multivariate statistics, this parameter, which we call BWSR (Biomass weighed sequence reads) showed very similar trends to species biomass and comparable patterns to species abundance, highlighting the potential of metabarcoding not only for biodiversity estimation and mapping of presence/absence of species, but also for quantitative assessment of zooplankton communities.