Evaluating Metabarcoding Markers for Identifying Zooplankton and Ichthyoplankton Communities to Species in the Salish Sea: Morphological Comparisons and Rare, Threatened or Invasive Species

Zooplankton and ichthyoplankton community assessments depend on species diagnostics, yet morphological identifications are time consuming, require taxonomic expertise, and are hampered by lack of diagnostic characters –particularly for larval stages. Metabarcoding can identify multiple species in communities from short DNA sequences in comparison to reference databases. To evaluate species resolution across phylogenetic groups and food webs of zooplankton and ichthyoplankton, we compare five metabarcode mitochondrial (mt)DNA markers from gene regions of (a)cytochrome c oxidase subunit I, (b) cytochrome b, (c) 16S ribosomal RNA, and (d) 12S ribosomal RNA for DNA extracted from net tows in the Northeastern Pacific Ocean’s Salish Sea across seven sites and two seasons. Species resolved by metabarcoding are compared to invertebrate morphological identifications and biomass estimates. Results indicate that species resolution for different zooplankton and ichthyoplankton taxa can markedly vary among gene regions and markers, and in comparison to morphological identifications. Thus, researchers seeking “universal” metabarcoding should take caution that several markers and gene regions likely will be needed; all will miss some taxa and yield incomplete overlap. Species resolution requires careful attention to taxon marker selection and coverage in reference sequence repositories. In summary, combined multi-marker metabarcode and morphological approaches improve broadscale zooplankton diagnostics.