Metabarcoding of zooplankton diversity within the Chukchi Borderland, Arctic Ocean: improved resolution from multi-gene markers and region-specific DNA databases

Abstract

DNA-based analyses have become powerful tools for characterizing the metazoan biodiversity of diverse marine ecosystems. Metabarcoding (i.e., large-scale taxonomic identification of complex samples via high-throughput sequencing of a DNA barcode region) frequently uses hypervariable regions of the nuclear eukaryotic 18S ribosomal RNA (rRNA) gene. However, species-level taxonomic identification is hampered by the conservative nature of the 18S gene in comparison to the mitochondrial cytochrome oxidase I (COI) barcode gene. Additionally, metabarcoding relies on reference DNA sequence databases for classification of millions of unknown sequence reads and molecular operational taxonomic units (OTUs); databases that are at present depauperate for marine zooplankton taxa. Here, we characterized the mesozooplankton community for the Chukchi Borderland (CBL) region, western Arctic Ocean, through metabarcoding analysis of the V4 and V9 hypervariable regions of 18S rRNA and a portion of COI. Characterization of zooplankton diversity for the epipelagic and upper mesopelagic layers (0–500 m) was based upon 17 metazoan taxonomic categories encompassing 24 orders in 14 classes. Taxonomic classification using the V4 and V9 markers was most reliable for orders, with copepods dominating OTU counts. To increase taxonomic resolution and allow detection of species, V4, V9, and COI OTUs were classified against DNA sequence databases for the Arctic Ocean for the subclass Copepoda. The geographic region-specific databases for 18S rRNA and COI resulted in the detection and identification of 6 genera and 49 species of copepods representing 23 families, a marked increase in the taxonomic classification of the 18S rRNA markers. The greatest copepod species diversity was captured with V4 (34 species) followed by COI (28 species) with the least copepod diversity detected by V9 (5 species). Our results demonstrate the power of using multiple gene markers, with DNA reference databases that are specific to the geographic region of interest, providing more accurate metabarcoding biodiversity measures for the copepod assemblages compared to universal metazoan sequence reference databases. Results from this study highlight the need for continued DNA barcode sequencing to increase species representation in the reference sequence databases that are crucial for accurate characterization of mesozooplankton communities.