Next generation sequencing reveals the hidden diversity of zooplankton assemblages

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.

Global mesozooplankton communities show lower connectivity in deep oceanic layers

Mesozooplankton is a key component of the ocean, regulating global processes such as the carbon pump, and ensuring energy transfer from lower to higher trophic levels. Yet, knowledge on mesozooplankton diversity, distribution and connectivity at global scale is still fragmented. To fill this gap, we applied DNA metabarcoding to mesozooplankton samples collected during the Malaspina-2010 circumnavigation expedition across the Atlantic, Indian and Pacific oceans from the surface to bathypelagic depths. We highlight the still scarce knowledge on global mesozooplankton diversity and identify the Indian Ocean and the deep sea as the oceanic regions with the highest proportion of hidden diversity. We report no consistent alpha-diversity patterns for mesozooplankton at a global scale, neither across vertical nor horizontal gradients. However, beta-diversity analysis suggests horizontal and vertical structuring of mesozooplankton communities mostly attributed to turnover and reveals an increase in mesozooplankton beta-diversity with depth, indicating reduced connectivity at deeper layers. Additionally, we identify a water mass type-mediated structuring of mesozooplankton bathypelagic communities instead of an oceanic basin-mediated as observed at upper layers. This suggests limited dispersal at deep ocean layers, most likely due to weaker currents and lower mixing of water mass types, thus reinforcing the importance of oceanic currents and barriers to dispersal in shaping global plankton communities.

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.

Potential PCR amplification bias in identifying complex ecological patterns: Higher species compositional homogeneity revealed in smaller-size coral reef zooplankton by metatranscriptomics

PCR-based high-throughput sequencing has permitted comprehensive resolution analyses of zooplankton diversity dynamics. However, significant methodological issues still surround analyses of complex bulk community samples, not least as in prevailing PCR-based approaches. Marine drifting animals-zooplankton-play essential ecological roles in the pelagic ecosystem, transferring energy and elements to higher trophic levels, such as fishes, cetaceans and others. In the present study, we collected 48 size-fractionated zooplankton samples in the vicinity of a coral reef island with environmental gradients. To investigate the spatiotemporal dynamics of zooplankton diversity patterns and the effect of PCR amplification biases across these complex communities, we first took metatranscriptomics approach. Comprehensive computational analyses revealed a clear pattern of higher/lower homogeneity in smaller/larger zooplankton compositions across samples respectively. Our study thus suggests changes in the role of dispersal across the sizes. Next, we applied in silico PCR to the metatranscriptomics datasets, in order to estimate the extent of PCR amplification bias. Irrespective of stringency criteria, we observed clear separations of size fraction sample clusters in both metatranscriptomics and in silico datasets. In contrast, the pattern-smaller-fractioned communities had higher compositional homogeneity than larger ones-was observed in the metatranscriptomics data but not in the in silico datasets. To investigate this discrepancy further, we analysed the mismatches of widely used mitochondrial CO1 primers and identified priming site mismatches likely driving PCR-based biases. Our results suggest the use of metatranscriptomics or, although less ideal, redesigning the CO1 primers is necessary to circumvent these issues.

Assessing the relevance of DNA metabarcoding compared to morphological identification for lake phytoplankton monitoring

Phytoplankton is a key biological group used to assess the ecological status of lakes. The classical monitoring approach relies on microscopic identification and counting of phytoplankton species, which is time-consuming and requires high taxonomic expertise. High-throughput sequencing, combined with metabarcoding, has recently demonstrated its potential as an alternative approach for plankton surveys. Several studies have confirmed the relevance of the diatom metabarcoding approach to calculate biotic indices based on species ecology. However, phytoplankton communities have not yet benefited from such validation. Here, by comparing the results obtained with the two methods (molecular and microscopic counting), we evaluated the relevance of metabarcoding approach for phytoplankton monitoring by considering different metrics: alpha diversity, taxonomic composition, community structure and a phytoplankton biotic index used to assess the trophic level of lakes. For this purpose, 55 samples were collected in four large alpine lakes (Aiguebelette, Annecy, Bourget, Geneva) during the year 2021. For each sample, a metabarcoding analysis based on two genetic markers (16S and 23S rRNA) was performed, in addition to the microscopic count. Regarding the trophic level of lakes, significant differences were found between index values obtained with the two approaches. The main hypothesis to explain these differences comes from the incompleteness, particularly at the species level, of the barcode reference library for the two genetic markers. It is therefore necessary to complete reference libraries for using such species-based biotic indices with metabarcoding data. Besides this, species richness and diversity were higher in the molecular inventories than in the microscopic ones. Moreover, despite differences in taxonomic composition of the floristic lists obtained by the two approaches, their community structures were similar. These results support the possibility of using metabarcoding for phytoplankton monitoring but in a different way. We suggest exploring alternative approaches to index development, such as a taxonomy-free approach.

Phytoplankton Community Integrity Index (PCII) - A potential supplementary tool for evaluating nutrient enrichment status of Hong Kong marine waters

Diagnosis of eutrophication requires evidence of disturbance to the balance of organisms. We describe a tool, the Plankton Community Integrity Index (PCII), derived from the Plankton Index (PI) for tracking change in the seasonal patterns of abundance of diatom and dinoflagellate lifeforms when plotted in state space. The tool uses a nutrient-minimum reference period to interpret PCII values as status indicators, with values close to 1 indicating "High" status and 0.6 a Biological Water Quality Criterion (BioWQC) target set at the "Fair"/"Good" status boundary. It has been applied to Hong Kong marine waters, using data from monthly samples from 1995 through 2021. A preliminary analysis, required for the PI method, confirmed monsoonal seasonality in the diatom lifeform. In 5 of the 9 water bodies examined, PCII time series correlated with those of Total Inorganic Nitrogen (TIN). Since 2020, all Water Control Zones met the operationally defined BioWQC target.

Assessing the Seasonal and Spatial Dynamics of Zooplankton through DNA Metabarcoding in a Temperate Estuary

Zooplankton are key components of estuarine trophic networks. However, routine monitoring is hindered by the difficulty of morphology-based identification. DNA-based methods allow us to circumvent some of these hurdles, providing precise species identifications regardless of the taxonomic expertise of the investigator or the developmental stage of the specimens. However, the process is dependent on the completeness of the reference libraries. In this study, we sought to evaluate the potential of DNA metabarcoding to assess the seasonal (summer, autumn, and early spring) and spatial dynamics of zooplankton (four locations spanning ca. 6 km) in the Lima estuary (NW Portugal). Two genetic markers were used: the cytochrome c oxidase subunit I and the V4 hypervariable region of the ribosomal 18S rRNA genes. Overall, 327 species were recovered, and both markers displayed minute overlap (7% were detected with both markers). Species richness, composition, and taxonomic distinctness were majorly influenced by the season, with a declining tendency from summer (highest number of exclusive species, n = 74) to spring. Second to season, the taxa composition was influenced by spatial variation where the most downstream site displayed the highest number of exclusive species, n = 53. A total of 16 non-indigenous species were detected using metabarcoding, but only one (Austrominus modestus) has been documented out in the estuary. In conclusion, both the seasonal and spatial gradients influenced the recovered richness, composition, and taxonomic distinctness, confirming the great aptitude of DNA metabarcoding for providing higher density monitoring and shedding new light on the composition and dynamics of complex zooplankton communities.

Metabarcoding analysis of marine zooplankton confirms the ecological role of a sheltered bight along an exposed continental shelf

Zooplankton plays an essential role in marine ecosystems as the link between primary producers (phytoplankton) and higher trophic levels in food webs, and as a dynamic pool of recruits for invertebrates and fish. Zooplankton communities are diverse with a patchy distribution at different spatial scales, influenced by oceanographic processes. The continental shelf of eastern South Africa is narrow and exposed to the western-boundary Agulhas Current, with some shelter against strong directional flow provided by the broader KwaZulu-Natal Bight, a coastal offset adjacent to an estuary. We compared zooplankton species richness, diversity and relative abundance of key taxa among sheltered and exposed shelf areas using metabarcoding and community analysis, to explore the ecological role of the bight in a highly dynamic ocean region. Metabarcoding recovered higher richness and diversity at a finer resolution than could previously be achieved with traditional microscopy. Of 271 operational taxonomic units (OTUs) recovered through metabarcoding, 63% could be matched with >95% sequence similarity to reference barcodes. OTUs were dominated by malacostracan crustaceans (161 spp.), ray-finned fishes (45 spp.) and copepods (28 spp.). Species richness, diversity and the relative abundance of key taxa differed between sheltered and exposed shelf areas. Lower species richness in the bight was partly attributed to structurally homogeneous benthic habitats, and an associated reduction of meroplanktonic species originating from local benthic-pelagic exchange. High relative abundance of a ray-finned fish in the bight, as observed based on fish eggs and read counts, confirmed that the bight is an important fish spawning area. Overall, zooplankton metabarcoding outputs were congruent with findings of previous ecological research using more traditional methods of observation.

Environmental DNA metabarcoding serves as a promising method for aquatic species monitoring and management: A review focused on its workflow, applications, challenges and prospects

Marine and freshwater biodiversity is under threat from both natural and manmade causes. Biological monitoring is currently a top priority for biodiversity protection. Given present limitations, traditional biological monitoring methods may not achieve the proposed monitoring aims. Environmental DNA metabarcoding technology reflects species information by capturing and extracting DNA from environmental samples, using molecular biology techniques to sequence and analyze the DNA, and comparing the obtained information with existing reference libraries to obtain species identification. However, its practical application has highlighted several limitations. This paper summarizes the main steps in the environmental application of eDNA metabarcoding technology in aquatic ecosystems, including the discovery of unknown species, the detection of invasive species, and evaluations of biodiversity. At present, with the rapid development of big data and artificial intelligence, certain advanced technologies and devices can be combined with environmental DNA metabarcoding technology to promote further development of aquatic species monitoring and management.

Improving taxonomic classification of marine zooplankton by molecular approach: registration of taxonomically verified 18S and 28S rRNA gene sequences

Background

Zooplankton plays an important role in the marine ecosystem. A high level of taxonomic expertise is necessary for accurate species identification based on morphological characteristics. As an alternative method to morphological classification, we focused on a molecular approach using 18S and 28S ribosomal RNA (rRNA) gene sequences. This study investigates how the accuracy of species identification by metabarcoding improves when taxonomically verified sequences of dominant zooplankton species are added to the public database. The improvement was tested by using natural zooplankton samples.

Methods

rRNA gene sequences were obtained from dominant zooplankton species from six sea areas around Japan and registered in the public database for improving the accuracy of taxonomic classifications. Two reference databases with and without newly registered sequences were created. Comparison of detected OTUs associated with single species between the two references was done using field-collected zooplankton samples from the Sea of Okhotsk for metabarcoding analysis to verify whether or not the newly registered sequences improved the accuracy of taxonomic classifications.

Results

A total of 166 sequences in 96 species based on the 18S marker and 165 sequences in 95 species based on the 28S marker belonging to Arthropoda (mostly Copepoda) and Chaetognatha were registered in the public database. The newly registered sequences were mainly composed of small non-calanoid copepods, such as species belonging to Oithona and Oncaea. Based on the metabarcoding analysis of field samples, a total of 18 out of 92 OTUs were identified at the species level based on newly registered sequences in the data obtained by the 18S marker. Based on the 28S marker, 42 out of 89 OTUs were classified at the species level based on taxonomically verified sequences. Thanks to the newly registered sequences, the number of OTUs associated with a single species based on the 18S marker increased by 16% in total and by 10% per sample. Based on the 28S marker, the number of OTUs associated with a single species increased by 39% in total and by 15% per sample. The improved accuracy of species identification was confirmed by comparing different sequences obtained from the same species. The newly registered sequences had higher similarity values (mean >0.003) than the pre-existing sequences based on both rRNA genes. These OTUs were identified at the species level based on sequences not only present in the Sea of Okhotsk but also in other areas.

Discussion

The results of the registration of new taxonomically verified sequences and the subsequent comparison of databases based on metabarcoding data of natural zooplankton samples clearly showed an increase in accuracy in species identification. Continuous registration of sequence data covering various environmental conditions is necessary for further improvement of metabarcoding analysis of zooplankton for monitoring marine ecosystems.

High-Throughput DNA Metabarcoding as an Approach for Ichthyoplankton Survey in Oujiang River Estuary, China

High-throughput DNA metabarcoding of mitochondrial 12S rRNA and Cyt b gene sequences was coupled with a morphology-based identification tool to assess ichthyoplankton community structure in Oujiang River Estuary, China. The performances of 12S and Cyt b barcoding markers were compared in terms of taxonomic resolution, detection and coverage, and their suitability was established for use as a quick and powerful ichthyoplankton assessment tool. A total of 30,138 ichthyoplankton (2462 eggs and 27,676 larvae) samples were collected from April to August 2015 and identified to 145 taxa belonging to 57 families and 105 genera. June and July were the main spawning months. Ichthyoplankton were more abundant around Lingkun and Qidu Islands and the upper parts of Oujiang River Estuary. The 12S gene marker presented higher species coverage and detection rate than Cyt b. DNA metabarcoding exhibited more representative species identification power than morphology. The findings reported in this study provided a key attempt towards the development of time-efficient and cost-effective ichthyoplankton identification and assessment tool.

Complementary molecular and visual sampling of fish on oil and gas platforms provides superior biodiversity characterisation

Offshore oil and gas platforms have the potential to provide complex refugia for fish and benthic colonisers. We compare two methods of biodiversity assessment for fish and elasmobranchs at seven decommissioned oil and gas platforms as well as five sediment sites, located 5 km from platforms, in the Gulf of Thailand. Using surveys from stereo-video ROV transects, and data from Environmental DNA (eDNA) water-column samples, we detected fish and elasmobranch taxa from 39 families and 66 genera across both platform and sediment sites with eDNA, compared with 18 families and 29 genera by stereo-ROV with platforms yielding significantly greater species richness. This study demonstrates that the combination of stereo-video ROV and eDNA provide effective, non-extractive and complementary methods to enhance data capture. This approach sets new benchmarks for evaluating fish assemblages surrounding platforms and will enhance measurements of biota to inform decisions on the fate of oil/gas infrastructure.

A triad of kicknet sampling, eDNA metabarcoding, and predictive modeling to assess richness of mayflies, stoneflies and caddisflies in rivers

Monitoring biodiversity is essential to understand the impacts of human activities and for effective management of ecosystems. Thereby, biodiversity can be assessed through direct collection of targeted organisms, through indirect evidence of their presence (e.g. signs, environmental DNA, camera trap, etc.), or through extrapolations from species distribution and species richness models. Differences in approaches used in biodiversity assessment, however, may come with individual challenges and hinder cross-study comparability. In the context of rapidly developing techniques, we compared three different approaches in order to better understand assessments of aquatic macroinvertebrate diversity. Specifically, we compared the community composition and species richness of three orders of aquatic macroinvertebrates (mayflies, stoneflies, and caddisflies, hereafter EPT) obtained via eDNA metabarcoding and via traditional in situ kicknet sampling to catchment-level based predictions of a species richness model. We used kicknet data from 24 sites in Switzerland and compared taxonomic lists to those obtained using eDNA amplified with two different primer sets. Richness detected by these methods was compared to the independent predictions made by a statistical species richness model, that is, a generalized linear model using landscape-level features to estimate EPT diversity. Despite the ability of eDNA to consistently detect some EPT species found by traditional sampling, we found important discrepancies in community composition between the kicknet and eDNA approaches, particularly at a local scale. We found the EPT-specific primer set fwhF2/EPTDr2n, detected a greater number of targeted EPT species compared to the more general primer set mlCOIintF/HCO2198. Moreover, we found that the species richness measured by eDNA from either primer set was poorly correlated to the richness measured by kicknet sampling (Pearson correlation = 0.27) and that the richness estimated by eDNA and kicknet were poorly correlated with the prediction of the species richness model (Pearson correlation = 0.30 and 0.44, respectively). The weak relationships between the traditional kicknet sampling and eDNA with this model indicates inherent limitations in upscaling species richness estimates, and possibly a limited ability of the model to meet real world expectations. It is also possible that the number of replicates was not sufficient to detect ambiguous correlations. Future challenges include improving the accuracy and sensitivity of each approach individually, yet also acknowledging their respective limitations, in order to best meet stakeholder demands and address the biodiversity crisis we are facing.

Can a key boreal Calanus copepod species now complete its life-cycle in the Arctic? Evidence and implications for Arctic food-webs

The changing Arctic environment is affecting zooplankton that support its abundant wildlife. We examined how these changes are influencing a key zooplankton species, Calanus finmarchicus, principally found in the North Atlantic but expatriated to the Arctic. Close to the ice-edge in the Fram Strait, we identified areas that, since the 1980s, are increasingly favourable to C. finmarchicus. Field-sampling revealed part of the population there to be capable of amassing enough reserves to overwinter. Early developmental stages were also present in early summer, suggesting successful local recruitment. This extension to suitable C. finmarchicus habitat is most likely facilitated by the long-term retreat of the ice-edge, allowing phytoplankton to bloom earlier and for longer and through higher temperatures increasing copepod developmental rates. The increased capacity for this species to complete its life-cycle and prosper in the Fram Strait can change community structure, with large consequences to regional food-webs.

Meta-analysis shows both congruence and complementarity of DNA and eDNA metabarcoding to traditional methods for biological community assessment

DNA metabarcoding is increasingly used for the assessment of aquatic communities, and numerous studies have investigated the consistency of this technique with traditional morpho‐taxonomic approaches. These individual studies have used DNA metabarcoding to assess diversity and community structure of aquatic organisms both in marine and freshwater systems globally over the last decade. However, a systematic analysis of the comparability and effectiveness of DNA‐based community assessment across all of these studies has hitherto been lacking. Here, we performed the first meta‐analysis of available studies comparing traditional methods and DNA metabarcoding to measure and assess biological diversity of key aquatic groups, including plankton, microphytobentos, macroinvertebrates, and fish. Across 215 data sets, we found that DNA metabarcoding provides richness estimates that are globally consistent to those obtained using traditional methods, both at local and regional scale. DNA metabarcoding also generates species inventories that are highly congruent with traditional methods for fish. Contrastingly, species inventories of plankton, microphytobenthos and macroinvertebrates obtained by DNA metabarcoding showed pronounced differences to traditional methods, missing some taxa but at the same time detecting otherwise overseen diversity. The method is generally sufficiently advanced to study the composition of fish communities and replace more invasive traditional methods. For smaller organisms, like macroinvertebrates, plankton and microphytobenthos, DNA metabarcoding may continue to give complementary rather than identical estimates compared to traditional approaches. Systematic and comparable data collection will increase the understanding of different aspects of this complementarity, and increase the effectiveness of the method and adequate interpretation of the results.

Experimental validation of taxon-specific mini-barcode primers for metabarcoding of zooplankton

Metabarcoding to determine the species composition and diversity of marine zooplankton communities is a fast-developing field in which the standardization of methods is yet to be fully achieved. The selection of genetic markers and primer choice are particularly important because they substantially influence species detection rates and accuracy. Validation is therefore an important step in the design of metabarcoding protocols. We developed taxon-specific mini-barcode primers for the cytochrome c oxidase subunit I (COI) gene region and used an experimental approach to test species detection rates and primer accuracy of the newly designed primers for prawns, shrimps and crabs and published primers for marine lobsters and fish. Artificially assembled mock communities (with known species ratios) and unsorted coastal tow-net zooplankton samples were sequenced and the detected species were compared with those seeded in mock communities to test detection rates. Taxon-specific primers increased detection rates of target taxa compared with a universal primer set. Primer cocktails (multiple primer sets) significantly increased species detection rates compared with single primer pairs and could detect up to 100% of underrepresented target taxa in mock communities. Taxon-specific primers recovered fewer false-positive or false-negative results than the universal primer. The methods used to design taxon-specific mini-barcodes and the experimental mock community validation protocols shown here can easily be applied to studies on other groups and will allow for a level of standardization among studies undertaken in different ecosystems or geographic locations.

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.

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.

Environmental conditions drive zooplankton community structure in the epipelagic oceanic water of the southern Gulf of Mexico: A molecular approach

Zooplankton plays a pivotal role in sustaining the majority of marine ecosystems. The distribution patterns and diversity of zooplankton provide key information for understanding the functioning of these ecosystems. Nevertheless, due to the numerous cryptic and sibling species and the lack of diagnostic characteristics for early developmental stages, the identification of the global-to-local patterns of zooplankton biodiversity and biogeography remains challenging in different research fields. The spatial and temporal changes in the zooplankton community in the open waters of the southern Gulf of Mexico were assessed using metabarcoding analysis of the V9 region of 18S rRNA and mitochondrial cytochrome oxidase c subunit I (COI). Additionally, a multiscale analysis was implemented to evaluate which environmental predictors may explain the variability in the structure of the zooplankton community. Our findings suggest that the synergistic effects of dissolved oxygen concentration, temperature, and longitude (intended as a proxy for still unidentified predictors) may explain both spatial and temporal zooplankton variability even with low contribution. Furthermore, the zooplankton distribution probably reflects the coexistence of three heterogeneous ecoregions and a bio-physical partitioning of the studied area. Finally, some taxa were either exclusive or predominant with either 18S or COI markers. This may suggest that comprehensive assessments of the zooplankton community may be more accurately met by the use of multilocus approaches.

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.

Evaluating the performance of DNA metabarcoding for assessment of zooplankton communities in Western Lake Superior using multiple markers

For DNA metabarcoding to attain its potential as a community assessment tool, we need to better understand its performance versus traditional morphological identification and work to address any remaining performance gaps in incorporating DNA metabarcoding into community assessments. Using fragments of the 18S nuclear and 16S mitochondrial rRNA genes and two fragments of the mitochondrial COI marker, we examined the use of DNA metabarcoding and traditional morphological identification for understanding the diversity and composition of crustacean zooplankton at 42 sites across western Lake Superior. We identified 51 zooplankton taxa (genus or species, depending on the finest resolution of the taxon across all identification methods), of which 17 were identified using only morphological traits, 13 using only DNA and 21 using both methods. The taxa found using only DNA metabarcoding included four species and one genus-level identification not previously known to occur in Lake Superior, the presence of which still needs to be confirmed. A substantial portion of taxa that were identified to genus or species by morphological identification, but not identified using DNA metabarcoding, had zero ("no record") or ≤ 2 ("underrepresented records") reference barcodes in the BOLD or NCBI databases (63% for COI, 80% for 16S, 74% for 18S). The two COI marker fragments identified the most genus- and species-level taxa, whereas 18S was the only marker whose family-level percent sequence abundance patterns showed high correlation to composition patterns from morphological identification, based on a NMDS analysis of Bray-Curtis similarities. Multiple replicates were collected at a subset of sites and an occupancy analysis was performed, which indicated that rare taxa were more likely to be detected using DNA metabarcoding than traditional morphology. Our results support that DNA metabarcoding can augment morphological identification for estimating zooplankton diversity and composition of zooplankton over space and time, but may require use of multiple markers. Further addition of taxa to reference DNA databases will improve our ability to use DNA metabarcoding to identify zooplankton and other invertebrates in aquatic surveys.

Building on gAMBI in ports for a challenging biological invasions scenario: Blue-gNIS as a proof of concept

The status of aquatic ecosystems has historically been monitored by the use of biotic indices. However, few biotic measures consider the presence of non-indigenous species as a sign of anthropogenic pollution and habitat disturbance even when this may seriously affect the metric scores and ecological status classifications of an environment. Today, biological invasions are currently one of the greatest threats to biodiversity and sustainable blue economies around the world. In this work, environmental assessments were conducted in the Port of Gijon, Northern Spain, using eDNA metabarcoding, and the gAMBI (genetics based AZTI Marine Biotic Index) was estimated. Results indicate a high/good ecological status within the port. However, nine non-indigenous species and five invasive species were found, and a modification of the gAMBI that includes species invasiveness was proposed: Blue-gNIS. The index was preliminary tested against existing validated indices such as gAMBI, BENTIX (based on the ecology of macroinvertebrates) and ALEX (based on the invasiveness of the species). Blue-gNIS classified the port in a good ecological status and showed its potential usefulness to achieve more complete water quality assessments of ports.

DNA metabarcoding of zooplankton communities: species diversity and seasonal variation revealed by 18S rRNA and COI

Background

Zooplankton is an important component of aquatic organisms and has important biological and economical significance in freshwater ecosystems. However, traditional methods that rely on morphology to classify zooplankton require expert taxonomic skills. Moreover, traditional classification methods are time-consuming and labor-intensive, which is not practical for the design of conservation measures and ecological management tools based on zooplankton diversity assessment.

Methods

We used DNA metabarcoding technology with two different markers: the nuclear small subunit ribosomal RNA (18S rRNA) and mitochondrial cytochrome c oxidase (COI), to analyze 72 zooplankton samples collected in 4 seasons and 9 locations from the Sanmenxia Reservoir. We investigated seasonal changes in the zooplankton community and their relationship with water environmental factors.

Results

A total of 190 species of zooplankton were found, belonging to 12 phyla, 24 classes, 61 orders, 111 families, and 174 genera. Protozoa, especially ciliates, were the most diverse taxa. Richness and relative abundance of zooplankton showed significant seasonal changes. Both alpha and beta diversity showed seasonal trends: the diversity in summer and autumn was higher than that in winter and spring. The zooplankton diversity was most similar in winter and spring. By correlating metabarcoding data and water environmental factors, we proved that water temperature, chemical oxygen demand, total nitrogen and ammoniacal nitrogen were the main environmental factors driving the seasonal changes in zooplankton in the Sanmenxia Reservoir. Water temperature, followed by total nitrogen, were the most influential factors. This study highlights the advantages and some limitations of zooplankton molecular biodiversity assessment using two molecular markers.

Exploring the Use of Environmental DNA (eDNA) to Detect Animal Taxa in the Mesopelagic Zone

Animal biodiversity in the ocean’s vast mesopelagic zone is relatively poorly studied due to technological and logistical challenges. Environmental DNA (eDNA) analyses show great promise for efficiently characterizing biodiversity and could provide new insight into the presence of mesopelagic species, including those that are missed by traditional net sampling. Here, we explore the utility of eDNA for identifying animal taxa. We describe the results from an August 2018 cruise in Slope Water off the northeast United States. Samples for eDNA analysis were collected using Niskin bottles during five CTD casts. Sampling depths along each cast were selected based on the presence of biomass as indicated by the shipboard Simrad EK60 echosounder. Metabarcoding of the 18S V9 gene region was used to assess taxonomic diversity. eDNA metabarcoding results were compared with those from net-collected (MOCNESS) plankton samples. We found that the MOCNESS sampling recovered more animal taxa, but the number of taxa detected per liter of water sampled was significantly higher in the eDNA samples. eDNA was especially useful for detecting delicate gelatinous animals which are undersampled by nets. We also detected eDNA changes in community composition with depth, but not with sample collection time (day vs. night). We provide recommendations for applying eDNA-based methods in the mesopelagic including the need for studies enabling interpretation of eDNA signals and improvement of barcode reference databases.

The crossover from microscopy to genes in marine diversity: from species to assemblages in marine pelagic copepods

An accurate identification of species and communities is a prerequisite for analysing and recording biodiversity and community shifts. In the context of marine biodiversity conservation and management, this review outlines past, present and forward-looking perspectives on identifying and recording planktonic diversity by illustrating the transition from traditional species identification based on morphological diagnostic characters to full molecular genetic identification of marine assemblages. In this process, the article presents the methodological advancements by discussing progress and critical aspects of the crossover from traditional to novel and future molecular genetic identifications and it outlines the advantages of integrative approaches using the strengths of both morphological and molecular techniques to identify species and assemblages. We demonstrate this process of identifying and recording marine biodiversity on pelagic copepods as model taxon. Copepods are known for their high taxonomic and ecological diversity and comprise a huge variety of behaviours, forms and life histories, making them a highly interesting and well-studied group in terms of biodiversity and ecosystem functioning. Furthermore, their short life cycles and rapid responses to changing environments make them good indicators and core research components for ecosystem health and status in the light of environmental change. This article is part of the theme issue 'Integrative research perspectives on marine conservation'.

Protocols for Monitoring Harmful Algal Blooms for Sustainable Aquaculture and Coastal Fisheries in Chile

Harmful algae blooms (HABs) cause acute effects on marine ecosystems due to their production of endogenous toxins or their enormous biomass, leading to significant impacts on local economies and public health. Although HAB monitoring has been intensively performed at spatiotemporal scales in coastal areas of the world over the last decades, procedures have not yet been standardized. HAB monitoring procedures are complicated and consist of many methodologies, including physical, chemical, and biological water sample measurements. Each monitoring program currently uses different combinations of methodologies depending on site specific purposes, and many prior programs refer to the procedures in quotations. HAB monitoring programs in Chile have adopted the traditional microscopic and toxin analyses but not molecular biology and bacterial assemblage approaches. Here we select and optimize the HAB monitoring methodologies suitable for Chilean geography, emphasizing on metabarcoding analyses accompanied by the classical tools with considerations including cost, materials and instrument availability, and easiness and efficiency of performance. We present results from a pilot study using the standardized stepwise protocols, demonstrating feasibility and plausibility for sampling and analysis for the HAB monitoring. Such specific instructions in the standardized protocol are critical obtaining quality data under various research environments involving multiple stations, different analysts, various time-points, and long HAB monitoring duration.

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.

Sensitivity of Modeled CO2 Air–Sea Flux in a Coastal Environment to Surface Temperature Gradients, Surfactants, and Satellite Data Assimilation

This work evaluates the sensitivity of CO2 air–sea gas exchange in a coastal site to four different model system configurations of the 1D coupled hydrodynamic–ecosystem model GOTM–ERSEM, towards identifying critical dynamics of relevance when specifically addressing quantification of air–sea CO2 exchange. The European Sea Regional Ecosystem Model (ERSEM) is a biomass and functional group-based biogeochemical model that includes a comprehensive carbonate system and explicitly simulates the production of dissolved organic carbon, dissolved inorganic carbon and organic matter. The model was implemented at the coastal station L4 (4 nm south of Plymouth, 50°15.00’N, 4°13.02’W, depth of 51 m). The model performance was evaluated using more than 1500 hydrological and biochemical observations routinely collected at L4 through the Western Coastal Observatory activities of 2008–2009. In addition to a reference simulation (A), we ran three distinct experiments to investigate the sensitivity of the carbonate system and modeled air–sea fluxes to (B) the sea-surface temperature (SST) diurnal cycle and thus also the near-surface vertical gradients, (C) biological suppression of gas exchange and (D) data assimilation using satellite Earth observation data. The reference simulation captures well the physical environment (simulated SST has a correlation with observations equal to 0.94 with a p > 0.95). Overall, the model captures the seasonal signal in most biogeochemical variables including the air–sea flux of CO2 and primary production and can capture some of the intra-seasonal variability and short-lived blooms. The model correctly reproduces the seasonality of nutrients (correlation > 0.80 for silicate, nitrate and phosphate), surface chlorophyll-a (correlation > 0.43) and total biomass (correlation > 0.7) in a two year run for 2008–2009. The model simulates well the concentration of DIC, pH and in-water partial pressure of CO2 (pCO2) with correlations between 0.4–0.5. The model result suggest that L4 is a weak net source of CO2 (0.3–1.8 molCm−2 year−1). The results of the three sensitivity experiments indicate that both resolving the temperature profile near the surface and assimilation of surface chlorophyll-a significantly impact the skill of simulating the biogeochemistry at L4 and all of the carbonate chemistry related variables. These results indicate that our forecasting ability of CO2 air–sea flux in shelf seas environments and their impact in climate modeling should consider both model refinements as means of reducing uncertainties and errors in any future climate projections.

Biodiversity between sand grains: Meiofauna composition across southern and western Sweden assessed by metabarcoding

The meiofauna is an important part of the marine ecosystem, but its composition and distribution patterns are relatively unexplored. Here we assessed the biodiversity and community structure of meiofauna from five locations on the Swedish western and southern coasts using a high-throughput DNA sequencing (metabarcoding) approach. The mitochondrial cytochrome oxidase 1 (COI) mini-barcode and nuclear 18S small ribosomal subunit (18S) V1-V2 region were amplified and sequenced using Illumina MiSeq technology. Our analyses revealed a higher number of species than previously found in other areas: thirteen samples comprising 6.5 dm3 sediment revealed 708 COI and 1,639 18S metazoan OTUs. Across all sites, the majority of the metazoan biodiversity was assigned to Arthropoda, Nematoda and Platyhelminthes. Alpha and beta diversity measurements showed that community composition differed significantly amongst sites. OTUs initially assigned to Acoela, Gastrotricha and the two Platyhelminthes sub-groups Macrostomorpha and Rhabdocoela were further investigated and assigned to species using a phylogeny-based taxonomy approach. Our results demonstrate that there is great potential for discovery of new meiofauna species even in some of the most extensively studied locations.

Zooplankton biodiversity monitoring in polluted freshwater ecosystems: A technical review

Freshwater ecosystems harbor a vast diversity of micro-eukaryotes (rotifers, crustaceans and protists), and such diverse taxonomic groups play important roles in ecosystem functioning and services. Unfortunately, freshwater ecosystems and biodiversity therein are threatened by many environmental stressors, particularly those derived from intensive human activities such as chemical pollution. In the past several decades, significant efforts have been devoted to halting biodiversity loss to recover services and functioning of freshwater ecosystems. Biodiversity monitoring is the first and a crucial step towards diagnosing pollution impacts on ecosystems and making conservation plans. Yet, bio-monitoring of ubiquitous micro-eukaryotes is extremely challenging, owing to many technical issues associated with micro-zooplankton such as microscopic size, fuzzy morphological features, and extremely high biodiversity. Here, we review current methods used for monitoring zooplankton biodiversity to advance management of impaired freshwater ecosystems. We discuss the development of traditional morphology-based identification methods such as scanning electron microscope (SEM) and ZOOSCAN and FlowCAM automatic systems, and DNA-based strategies such as metabarcoding and real-time quantitative PCR. In addition, we summarize advantages and disadvantages of these methods when applied for monitoring impacted ecosystems, and we propose practical DNA-based monitoring workflows for studying biological consequences of environmental pollution in freshwater ecosystems. Finally, we propose possible solutions for existing technical issues to improve accuracy and efficiency of DNA-based biodiversity monitoring.

•

Freshwater ecosystems and associated biodiversity have been highly degraded.

•

Biodiversity monitoring is crucial for diagnosing degradation degrees.

•

Here we review available methods for monitoring zooplankton biodiversity.

•

We propose possible solutions for existing technical issues.

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.

Environmental DNA Sequencing Reveals a Highly Complex Eukaryote Community in Sansha Yongle Blue Hole, Xisha, South China Sea

We report an Illumina high-throughput sequencing protocol of eukaryotic microbes in the world’s deepest marine blue hole, Sansha Yongle Blue Hole, Xisha, South China Sea. The variable V9 region of small subunit (SSU) rDNA, was sequenced using this approach from the waters of blue hole and outer reef slope. 917,771 unique eukaryotic 18S rRNA gene sequences and 6093 operational taxonomic units (OTUs) were identified. Significant differences in the eukaryotic composition were observed between the blue hole and outer reef slope, and the richness in the blue hole was much higher than that in the outer reef slope. The richness and diversity of eukaryotes in the blue hole were both lowest at 60 m and highest at 100 m depth. Eukaryotic microalgae assemblages dominated by Dinophyceae were the most abundant in the 10–20 m water column in the hole. Fauna was the main group at and below a depth of 60 m, where Araneae and Cyclopoida were dominant in the 60 m and 80 m water layer, respectively. There was a large number of Entoprocta at a depth of 180 m in the hole, where little oxygen was detected. Turbidity and nitrite concentration had a significant effect on the eukaryote community structure (p < 0.01).

Prokaryotic and eukaryotic microbiomes associated with blooms of the ichthyotoxic dinoflagellate Cochlodinium (Margalefidinium) polykrikoides in New York, USA, estuaries

While harmful algal blooms caused by the ichthyotoxic dinoflagellate, Cochlodinium (Margalefidinium) polykrikoides, are allelopathic and may have unique associations with bacteria, a comprehensive assessment of the planktonic communities associated with these blooms has been lacking. Here, we used high-throughput amplicon sequencing to assess size fractionated (0.2 and 5 μm) bacterial (16S) and phytoplankton assemblages (18S) associated with blooms of C. polykrikoides during recurrent blooms in NY, USA. Over a three-year period, samples were collected inside (‘patch’) and outside (‘non-patch’) dense accumulations of C. polykrikoides to assess the microbiome associated with these blooms. Eukaryotic plankton communities of blooms had significantly lower diversity than non-bloom samples, and non-bloom samples hosted 30 eukaryotic operational taxonomic units (OTUs) not found within blooms, suggesting they may have been allelopathically excluded from blooms. Differential abundance analyses revealed that C. polykrikoides blooms were significantly enriched in dinoflagellates (p<0.001) and the experimental enrichment of C. polykrikoides led to a significant increase in the relative abundance of eight genera of dinoflagellates but a significant decline in other eukaryotic plankton. Amoebophrya co-dominated both within- and near- C. polykrikoides blooms and was more abundant in bloom patches. The core bacterial microbiome of the >0.2μm fraction of blooms was dominated by an uncultured bacterium from the SAR11 clade, while the >5μm size fraction was co-dominated by an uncultured bacterium from Rhodobacteraceae and Coraliomargarita. Two bacterial lineages within the >0.2μm fraction, as well as the Gammaproteobacterium, Halioglobus, from the >5μm fraction were unique to the microbiome of blooms, while there were 154 bacterial OTUs only found in non-bloom waters. Collectively, these findings reveal the unique composition and potential function of eukaryotic and prokaryotic communities associated with C. polykrikoides blooms.

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.

Comparison of Water Sampling between Environmental DNA Metabarcoding and Conventional Microscopic Identification: A Case Study in Gwangyang Bay, South Korea

Our study focuses on methodological comparison of plankton community composition in relation to ecological monitoring and assessment with data sampling. Recently, along with the advancement of monitoring techniques, metabarcoding has been widely used in the context of environmental DNA (eDNA). We examine the applicability of eDNA metabarcoding for effective monitoring and assessment of community composition, compared with conventional observation using microscopic identification in a coastal ecosystem, Gwangynag Bay in South Korea. Our analysis is based primarily on two surveys at a total of 15 study sites in early and late summer (June and September) of the year 2018. The results of our study demonstrate the similarity and dissimilarity of biological communities in composition, richness and diversity between eDNA metabarcoding and conventional microscopic identification. It is found that, overall, eDNA metabarcoding appears to provide a wider variety of species composition, while conventional microscopic identification depicts more distinct plankton communities in sites. Finally, we suggest that eDNA metabarcoding is a practically useful method and can be potentially considered as a valuable alternative for biological monitoring and diversity assessments.

Phylotypic Characterization of Mycobionts and Photobionts of Rock Tripe Lichen in East Antarctica

Saxicolous rock ripe lichens that grow on rocks in the East Antarctic fellfields were sampled for phylotypic characterization of its constituent mycobionts (fungi) and photobionts (algae and cyanobacteria). The rock tripe lichen-forming fungal and algal phylotypes were classified under the common lichen-forming genera of ascomycetes, namely, Umbilicaria, and green algae, namely, Trebouxia and Coccomyxa. However, phylotypes of the green algal chloroplasts and the lichen-associated cyanobacteria showed unexpectedly high diversity. The detected chloroplast phylotypes were not fully affiliated with the green algal genera Trebouxia or Coccomyxa. The predominant chloroplast phylotype demonstrated maximum resemblance to Neglectella solitaria, which is neither a known Antarctic species nor a typical lichen photobiont. Another dominant chloroplast phylotype belonged to the atypical Antarctic green algae family. Cyanobacterial phylotypes were dominated by those affiliated with the Microcoleus species rather than the well-known lichen-associates, Nostoc species. The occurrences of these Microcoleus-affiliated cyanobacterial phylotypes were specifically abundant within the Yukidori Valley site, one of the Antarctic Specially Protected Areas (ASPA). The ASPA site, along with another 50 km-distant site, yielded most of the cryptic diversity in the phylotypes of chloroplasts and cyanobacteria, which may contribute to the phenotypic variability within the rock tripe lichen photobionts.

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.

High diversity and potential translocation of DNA viruses in ballast water

Ballast water is a common vector for the transport of invasive species to new marine and aquatic environments. We used a metagenomics approach to examine the composition and diversity of viral communities in ballast water from ships originating in Mexico, Saudi Arabia, New York, and Panama, and in water from the port of their destination in Busan, Korea. Myoviridae was the most abundant virus family in ballast water, followed Podoviridae and Siphoviridae. We also identified viruses that infect invertebrates, amoebas, and algae in ballast water and in the Busan port water. Interestingly, there were several viruses that infect humans or other animals (Swinepox virus, Raccoonpox virus, Suid herpesvirus, and Human endogenous retrovirus) in the samples from New York and Panama. In addition, there were giant viruses in all the ballast water samples, especially, identified Megavirus chilensis in New York and Panama, and Pandoravirus salinus in Mexico and Saudi Arabia. These results provide detailed descriptions of the characteristics of the viruses present in ballast water, document significant viral diversity, and indicate the potential translocation of viruses via ballast water.

Multi-marker metabarcoding approach to study mesozooplankton at basin scale

Zooplankton plays a pivotal role in marine ecosystems and the characterisation of its biodiversity still represents a challenge for marine ecologists. In this study, mesozooplankton composition from 46 samples collected in summer along the western Adriatic Sea, was retrieved by DNA metabarcoding analysis. For the first time, the highly variable fragments of the mtDNA COI and the V9 region of 18S rRNA genes were used in a combined matrix to compile an inventory of mesozooplankton at basin scale. The number of sequences retrieved after quality filtering were 824,148 and 223,273 for COI and 18S (V9), respectively. The taxonomical assignment against reference sequences, using 95% (for COI) and 97% (for 18S) similarity thresholds, recovered 234 taxa. NMDS plots and cluster analysis divided coastal from offshore samples and the most representative species of these clusters were distributed according to the dominant surface current pattern of the Adriatic for the summer period. For selected sampling sites, mesozooplankton species were also identified under a stereo microscope providing insights on the strength and weakness of the two approaches. In addition, DNA metabarcoding was shown to be helpful for the monitoring of non-indigenous marine metazoans and spawning areas of commercial fish species. We defined pros and cons of applying this approach at basin scale and the benefits of combining the datasets from two genetic markers.

Metabarcoding analysis on European coastal samples reveals new molecular metazoan diversity

Although animals are among the best studied organisms, we still lack a full description of their diversity, especially for microscopic taxa. This is partly due to the time-consuming and costly nature of surveying animal diversity through morphological and molecular studies of individual taxa. A powerful alternative is the use of high-throughput environmental sequencing, providing molecular data from all organisms sampled. We here address the unknown diversity of animal phyla in marine environments using an extensive dataset designed to assess eukaryotic ribosomal diversity among European coastal locations. A multi-phylum assessment of marine animal diversity that includes water column and sediments, oxic and anoxic environments, and both DNA and RNA templates, revealed a high percentage of novel 18S rRNA sequences in most phyla, suggesting that marine environments have not yet been fully sampled at a molecular level. This novelty is especially high among Platyhelminthes, Acoelomorpha, and Nematoda, which are well studied from a morphological perspective and abundant in benthic environments. We also identified, based on molecular data, a potentially novel group of widespread tunicates. Moreover, we recovered a high number of reads for Ctenophora and Cnidaria in the smaller fractions suggesting their gametes might play a greater ecological role than previously suspected.

Omics Advances in Ecotoxicology

Toxic substances in the environment generate adverse effects at all levels of biological organization from the molecular level to community and ecosystem. Given this complexity, it is not surprising that ecotoxicologists have struggled to address the full consequences of toxic substance release at ecosystem level, due to the limits of observational and experimental tools to reveal the changes in deep structure at different levels of organization. -Omics technologies, consisting of genomics and ecogenomics, have the power to reveal, in unprecedented detail, the cellular processes of an individual or biodiversity of a community in response to environmental change with high sample/observation throughput. This represents a historic opportunity to transform the way we study toxic substances in ecosystems, through direct linkage of ecological effects with the systems biology of organisms. Three recent examples of -omics advance in the assessment of toxic substances are explored here: (1) the use of functional genomics in the discovery of novel molecular mechanisms of toxicity of chemicals in the environment; (2) the development of laboratory pipelines of dose-dependent, reduced transcriptomics to support high-throughput chemical testing at the biological pathway level; and (3) the use of eDNA metabarcoding approaches for assessing chemical effects on biological communities in mesocosm experiments and through direct observation in field monitoring. -Omics advances in ecotoxicological studies not only generate new knowledge regarding mechanisms of toxicity and environmental effect, improving the relevance and immediacy of laboratory toxicological assessment, but can provide a wholly new paradigm for ecotoxicology by linking ecological models to mechanism-based, systems biology approaches.

Evaluation of marine zooplankton community structure through environmental DNA metabarcoding

Zooplankton dominate the abundance and biomass of multicellular animals in pelagic marine environments; however, traditional methods to characterize zooplankton communities are invasive and laborious. This study compares zooplankton taxonomic composition revealed through metabarcoding of the cytochrome oxidase I (COI) and 18S rRNA genes to traditional morphological identification by microscopy. Triplicates of three different sample types were collected from three coral reef sites in the Florida Keys National Marine Sanctuary: (1) 1 L surface seawater samples prefiltered through 3 μm filters and subsequently collected on 0.22 μm filters for eDNA (PF-eDNA); (2) 1 L surface seawater samples filtered on 0.22 μm pore-size filters (environmental DNA; eDNA), and (3) zooplankton tissue samples from 64 μm, 200 μm, and 500 μm mesh size net tows. The zooplankton tissue samples were split, with half identified morphologically and tissue DNA (T-DNA) extracted from the other half. The COI and 18S rRNA gene metabarcoding of PF-eDNA, eDNA, and T-DNA samples was performed using Illumina MiSeq. Of the families detected with COI and 18S rRNA gene metabarcoding, 40% and 32%, respectively, were also identified through morphological assessments. Significant differences in taxonomic composition were observed between PF-DNA, eDNA, and T-DNA with both genetic markers. PF-eDNA resulted in detection of fewer taxa than the other two sample types; thus, prefiltering is not recommended. All dominant copepod taxa (> 5% of total abundance) were detected with eDNA, T-DNA, and morphological assessments, demonstrating that eDNA metabarcoding is a promising technique for future biodiversity assessments of pelagic zooplankton in marine systems.

Nucleic acids-based tools for ballast water surveillance, monitoring, and research

Understanding the risks of biological invasion posed by ballast water-whether in the context of compliance testing, routine monitoring, or basic research-is fundamentally an exercise in biodiversity assessment, and as such should take advantage of the best tools available for tackling that problem. The past several decades have seen growing application of genetic methods for the study of biodiversity, driven in large part by dramatic technological advances in nucleic acids analysis. Monitoring approaches based on such methods have the potential to increase dramatically sampling throughput for biodiversity assessments, and to improve on the sensitivity, specificity, and taxonomic accuracy of traditional approaches. The application of targeted detection tools (largely focused on PCR but increasingly incorporating novel probe-based methodologies) has led to a paradigm shift in rare species monitoring, and such tools have already been applied for early detection in the context of ballast water surveillance. Rapid improvements in community profiling approaches based on high throughput sequencing (HTS) could similarly impact broader efforts to catalogue biodiversity present in ballast tanks, and could provide novel opportunities to better understand the risks of biotic exchange posed by ballast water transport-and the effectiveness of attempts to mitigate those risks. These various approaches still face considerable challenges to effective implementation, depending on particular management or research needs. Compliance testing, for instance, remains dependent on accurate quantification of viable target organisms; while tools based on RNA detection show promise in this context, the demands of such testing require considerable additional investment in methods development. In general surveillance and research contexts, both targeted and community-based approaches are still limited by various factors: quantification remains a challenge (especially for taxa in larger size classes), gaps in nucleic acids reference databases are still considerable, uncertainties in taxonomic assignment methods persist, and many applications have not yet matured sufficiently to offer standardized methods capable of meeting rigorous quality assurance standards. Nevertheless, the potential value of these tools, their growing utilization in biodiversity monitoring, and the rapid methodological advances over the past decade all suggest that they should be seriously considered for inclusion in the ballast water surveillance toolkit.

Censusing marine eukaryotic diversity in the twenty-first century

The ocean constitutes one of the vastest and richest biomes on our planet. Most recent estimations, all based on indirect approaches, suggest that there are millions of marine eukaryotic species. Moreover, a large majority of these are small (less than 1 mm), cryptic and still unknown to science. However, this knowledge gap, caused by the lack of diagnostic morphological features in small organisms and the limited sampling of the global ocean, is currently being filled, thanks to new DNA-based approaches. The molecular technique of PCR amplification of homologous gene regions combined with high-throughput sequencing, routinely used to census unculturable prokaryotes, is now also being used to characterize whole communities of marine eukaryotes. Here, we review how this methodological advancement has helped to better quantify the magnitude and patterns of marine eukaryotic diversity, with an emphasis on taxonomic groups previously largely overlooked. We then discuss obstacles remaining to achieve a global understanding of marine eukaryotic diversity. In particular, we argue that 18S variable regions do not provide sufficient taxonomic resolution to census marine life, and suggest combining broad eukaryotic surveys targeting the 18S rRNA region with more taxon-focused analyses of hypervariable regions to improve our understanding of the diversity of species, the functional units of marine ecosystems.

This article is part of the themed issue ‘From DNA barcodes to biomes’.

Indigenous species barcode database improves the identification of zooplankton

Incompleteness and inaccuracy of DNA barcode databases is considered an important hindrance to the use of metabarcoding in biodiversity analysis of zooplankton at the species-level. Species barcoding by Sanger sequencing is inefficient for organisms with small body sizes, such as zooplankton. Here mitochondrial cytochrome c oxidase I (COI) fragment barcodes from 910 freshwater zooplankton specimens (87 morphospecies) were recovered by a high-throughput sequencing platform, Ion Torrent PGM. Intraspecific divergence of most zooplanktons was < 5%, except Branchionus leydign (Rotifer, 14.3%), Trichocerca elongate (Rotifer, 11.5%), Lecane bulla (Rotifer, 15.9%), Synchaeta oblonga (Rotifer, 5.95%) and Schmackeria forbesi (Copepod, 6.5%). Metabarcoding data of 28 environmental samples from Lake Tai were annotated by both an indigenous database and NCBI Genbank database. The indigenous database improved the taxonomic assignment of metabarcoding of zooplankton. Most zooplankton (81%) with barcode sequences in the indigenous database were identified by metabarcoding monitoring. Furthermore, the frequency and distribution of zooplankton were also consistent between metabarcoding and morphology identification. Overall, the indigenous database improved the taxonomic assignment of zooplankton.

Vertical gradients in species richness and community composition across the twilight zone in the North Pacific Subtropical Gyre

Although metazoan animals in the mesopelagic zone play critical roles in deep pelagic food webs and in the attenuation of carbon in midwaters, the diversity of these assemblages is not fully known. A metabarcoding survey of mesozooplankton diversity across the epipelagic, mesopelagic and upper bathypelagic zones (0-1500 m) in the North Pacific Subtropical Gyre revealed far higher estimates of species richness than expected given prior morphology-based studies in the region (4,024 OTUs, 10-fold increase), despite conservative bioinformatic processing. Operational taxonomic unit (OTU) richness of the full assemblage peaked at lower epipelagic-upper mesopelagic depths (100-300 m), with slight shoaling of maximal richness at night due to diel vertical migration, in contrast to expectations of a deep mesopelagic diversity maximum as reported for several plankton groups in early systematic and zoogeographic studies. Four distinct depth-stratified species assemblages were identified, with faunal transitions occurring at 100 m, 300 m and 500 m. Highest diversity occurred in the smallest zooplankton size fractions (0.2-0.5 mm), which had significantly lower % OTUs classified due to poor representation in reference databases, suggesting a deep reservoir of poorly understood diversity in the smallest metazoan animals. A diverse meroplankton assemblage also was detected (350 OTUs), including larvae of both shallow and deep living benthic species. Our results provide some of the first insights into the hidden diversity present in zooplankton assemblages in midwaters, and a molecular reappraisal of vertical gradients in species richness, depth distributions and community composition for the full zooplankton assemblage across the epipelagic, mesopelagic and upper bathypelagic zones.