A metagenetic approach for revealing community structure of marine planktonic copepods

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.

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 detects biodiversity and ecological features of phytoplankton communities in Mediterranean transitional waters

Climate changes and anthropogenic pressures are causing a biodiversity decline in terms of species number and genetic diversity, reducing the adaptability and evolvability of natural communities. Transitional water ecosystems are more sensitive to habitat reduction and degradation and, thus, are more exposed to biodiversity declines requiring biodiversity monitoring programs for their conservation. Environmental DNA (eDNA) metabarcoding represents a high-throughput tool for biodiversity assessment that is facilitating data collection for biodiversity monitoring. In this study, we applied, for the first time, eDNA metabarcoding in a Mediterranean coastal lagoon to assess the ecological features of eukaryotic phytoplankton communities. We sampled water in seven different lagoon sites and amplified the extracted DNA with primers targeting the variable region 4 (V4) of the 18S rRNA gene marker. The results demonstrated the validity of eDNA studies to provide insights into lagoon phytoplankton composition, establish the structure and spatial variation of phytoplankton communities, and evaluate its correlation to abiotic factors. Finally, the genetic distances analysis suggests that the different spatial distribution of OTUs, at least for the Tetraselmis genus, reflects the genetic background.

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.

A powerful long metabarcoding method for the determination of complex diets from faecal analysis of the European pond turtle (Emys orbicularis, L. 1758)

High‐throughput sequencing has become an accurate method for the identification of species present in soil, water, faeces, gut or stomach contents. However, information at the species level is limited due to the choice of short barcodes and based on the idea that DNA is too degraded to allow longer sequences to be amplified. We have therefore developed a long DNA metabarcoding method based on the sequencing of short reads followed by de novo assembly, which can precisely identify the taxonomic groups of organisms associated with complex diets, such as omnivorous individuals. The procedure includes 11 different primer pairs targeting the COI gene, the large subunit of the ribulose‐1,5‐bisphosphate carboxylase gene, the maturase K gene, the 28S rRNA and the trnL‐trnF chloroplastic region. We validated this approach using 32 faeces samples from an omnivorous reptile, the European pond turtle (Emys orbicularis, L. 1758). This metabarcoding approach was assessed using controlled experiments including mock communities and faecal samples from captive feeding trials. The method allowed us to accurately identify prey DNA present in the diet of the European pond turtles to the species level in most of the cases (82.4%), based on the amplicon lengths of multiple markers (168–1,379 bp, average 546 bp), and produced by de novo assembly. The proposed approach can be adapted to analyse various diets, in numerous conservation and ecological applications. It is consequently appropriate for detecting fine dietary variations among individuals, populations and species as well as for the identification of rare food items.

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

High-throughput sequencing of 5S-IGS in oaks: Exploring intragenomic variation and algorithms to recognize target species in pure and mixed samples

Measuring biological diversity is a crucial but difficult undertaking, as exemplified in oaks where complex patterns of morphological, ecological, biogeographical and genetic differentiation collide with traditional taxonomy, which measures biodiversity in number of species (or higher taxa). In this pilot study, we generated high-throughput sequencing amplicon data of the intergenic spacer of the 5S nuclear ribosomal DNA cistron (5S-IGS) in oaks, using six mock samples that differ in geographical origin, species composition and pool complexity. The potential of the marker for automated genotaxonomy applications was assessed using a reference data set of 1,770 5S-IGS cloned sequences, covering the entire taxonomic breadth and distribution range of western Eurasian Quercus, and applying similarity (blast) and evolutionary approaches (maximum-likelihood trees and Evolutionary Placement Algorithm). Both methods performed equally well, allowing correct identification of species in sections Ilex and Cerris in the pure and mixed samples, and main lineages shared by species of sect. Quercus. Application of different cut-off thresholds revealed that medium- to high-abundance (>10 or 25) sequences suffice for a net species identification of samples containing one or a few individuals. Lower thresholds identify phylogenetic correspondence with all target species in highly mixed samples (analogous to environmental bulk samples) and include rare variants pointing towards reticulation, incomplete lineage sorting, pseudogenic 5S units and in situ (natural) contamination. Our pipeline is highly promising for future assessments of intraspecific and interpopulation diversity, and of the genetic resources of natural ecosystems, which are fundamental to empower fast and solid biodiversity conservation programmes worldwide.

Nematodes as bioindicators of polluted sediments using metabarcoding and microscopic taxonomy

The use of bioindicator species is a widely applied approach to evaluate ecological conditions, and several indices have been designed for this purpose. To assess the impact of pollution, especially in sediments, a pollution-sensitive index based on nematodes, one of the most abundant and species-rich groups of metazoa, was developed. The NemaSPEAR[%] index in its original form relies on the morphological inspection of nematode species. The application of a morphologically based NemaSPEAR[%] at the genus-level was previously validated. The present study evaluated a NemaSPEAR[%] index based on metabarcoding of nematode communities and tested the potential of fragments from the 28S rDNA, 18S rDNA and cytochrome c oxidase subunit I (COI) genes. In general, molecular-based results tended to show a poorer condition than morphology-based results for the investigated sites. At the genus level, NemaSPEAR[%] values based on morphological data strongly correlated with those based on molecular data for both the 28S rDNA and the 18S rDNA gene fragments (R² = 0.86 and R² = 0.74, respectively). Within the dominant genera (>3%) identified by morphology, 68% were detected by at least one of the two ribosomal markers. At the species level, however, concordance was less pronounced, as there were several deviations of the molecular from the morphological data. These differences could mostly be attributed to shortcomings in the reference database used in the molecular-based assignments. Our pilot study shows that a molecularly based, genus-level NemaSPEAR[%] can be successfully applied to evaluate polluted sediment. Future studies need to validate this approach further, e.g. with bulk extractions of whole meiofaunal communities in order to circumvent time-consuming nematode isolation. Further database curation with abundant NemaSPEAR[%] species will also increase the applicability of this approach.

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

Large-scale metabarcoding analysis of epipelagic and mesopelagic copepods in the Pacific

A clear insight into the large-scale community structure of planktonic copepods is critical to understanding the mechanisms controlling diversity and biogeography of marine taxa in terms of their high abundance, ubiquity, and sensitivity to environmental changes. Here, we applied a 28S metabarcoding approach to large-scale communities of epipelagic and mesopelagic copepods at 70 stations across the Pacific Ocean and three stations in the Arctic Ocean. Major patterns of community structure and diversity, influenced by water mass structures, agreed with results from previous morphology-based studies. However, a large-scale metabarcoding approach could detect community changes even under stable environmental conditions, including changes in the north/south subtropical gyres and east/west areas within each subtropical gyre. There were strong effects of the epipelagic environment on mesopelagic communities, and community subdivisions were observed in the environmentally stable mesopelagic layer. In each sampling station, higher operational taxonomic unit (OTU) numbers and lower phylogenetic diversity were observed in the mesopelagic layer than in the epipelagic layer, indicating a recent rapid increase in species numbers in the mesopelagic layer. The phylogenetic analysis utilizing representative sequences of OTUs revealed trends of recent emergence of cold-water OTUs, which are mainly distributed at high latitudes with low water temperatures. Conversely, the high diversity of copepods at low latitudes was suggested to have been formed through long evolution under high water temperature conditions. The metabarcoding results suggest that evolutionary processes have strong impacts on current patterns of copepod diversity, and support the “out of the tropics” theory explaining latitudinal diversity gradients of copepods. Diversity patterns in both epipelagic and mesopelagic copepods was highly correlated to sea surface temperature; thus, predicted global warming may have a significant impact on copepod diversity in both layers.

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.

Metabarcoding data allow for reliable biomass estimates in the most abundant animals on earth

Microscopic organisms are the dominant and most diverse organisms on Earth. Nematodes, as part of this microscopic diversity, are by far the most abundant animals and their diversity is equally high. Molecular metabarcoding is often applied to study the diversity of microorganisms, but has yet to become the standard to determine nematode communities. As such, the information metabarcoding provides, such as in terms of species coverage, taxonomic resolution and especially if sequence reads can be linked to the abundance or biomass of nematodes in a sample, has yet to be determined. Here, we applied metabarcoding using three primer sets located within ribosomal rRNA gene regions to target assembled mock-communities consisting of 18 different nematode species that we established in 9 different compositions. We determined abundances and biomass of all species added to examine if relative sequence abundance or biomass can be linked to relative sequence reads. We found that nematode communities are not equally represented by the three different primer sets and we found that relative read abundances almost perfectly correlated positively with relative species biomass for two of the primer sets. This strong biomass-read number correlation suggests that metabarcoding reads can reveal biomass information even amongst more complex nematode communities as present in the environment and possibly can be transferred to better study other groups of organisms. This biomass-read link is of particular importance for more reliably assessing nutrient flow through food-webs, as well as adjusting biogeochemical models through user-friendly and easily obtainable metabarcoding data.

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.

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.

How quantitative is metabarcoding: A meta-analytical approach

Metabarcoding has been used in a range of ecological applications such as taxonomic assignment, dietary analysis and the analysis of environmental DNA. However, after a decade of use in these applications there is little consensus on the extent to which proportions of reads generated corresponds to the original proportions of species in a community. To quantify our current understanding, we conducted a structured review and meta‐analysis. The analysis suggests that a weak quantitative relationship may exist between the biomass and sequences produced (slope = 0.52 ± 0.34, p < 0.01), albeit with a large degree of uncertainty. None of the tested moderators, sequencing platform type, the number of species used in a trial or the source of DNA, were able to explain the variance. Our current understanding of the factors affecting the quantitative performance of metabarcoding is still limited: additional research is required before metabarcoding can be confidently utilized for quantitative applications. Until then, we advocate the inclusion of mock communities when metabarcoding as this facilitates direct assessment of the quantitative ability of any given study.

Why the COI barcode should be the community DNA metabarcode for the metazoa

Metabarcoding of complex metazoan communities is increasingly being used to measure biodiversity in terrestrial, freshwater and marine ecosystems, revolutionizing our ability to observe patterns and infer processes regarding the origin and conservation of biodiversity. A fundamentally important question is which genetic marker to amplify, and although the mitochondrial cytochrome oxidase subunit I (COI) gene is one of the more widely used markers in metabarcoding for the Metazoa, doubts have recently been raised about its suitability. We argue that (a) the extensive coverage of reference sequence databases for COI; (b) the variation it presents; (c) the comparative advantages for denoising protein-coding genes; and (d) recent advances in DNA sequencing protocols argue in favour of standardizing for the use of COI for metazoan community samples. We also highlight where research efforts should focus to maximize the utility of metabarcoding.

Biodiversity estimates and ecological interpretations of meiofaunal communities are biased by the taxonomic approach

Accurate assessments of biodiversity are crucial to advising ecosystem-monitoring programs and understanding ecosystem function. Nevertheless, a standard operating procedure to assess biodiversity accurately and consistently has not been established. This is especially true for meiofauna, a diverse community (>20 phyla) of small benthic invertebrates that have fundamental ecological roles. Recent studies show that metabarcoding is a cost-effective and time-effective method to estimate meiofauna biodiversity, in contrast to morphological-based taxonomy. Here, we compare biodiversity assessments of a diverse meiofaunal community derived by applying multiple taxonomic methods based on comparative morphology, molecular phylogenetic analysis, DNA barcoding of individual specimens, and metabarcoding of environmental DNA. We show that biodiversity estimates are strongly biased across taxonomic methods and phyla. Such biases affect understanding of community structures and ecological interpretations. This study supports the urgency of improving aspects of environmental high-throughput sequencing and the value of taxonomists in correctly understanding biodiversity estimates.

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.

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

DNA metabarcoding is an efficient method for measuring biodiversity, but the process of initiating long-term DNA-based monitoring programmes, or integrating with conventional programs, is only starting. In marine ecosystems, plankton surveys using the continuous plankton recorder (CPR) have characterized biodiversity along transects covering millions of kilometres with time-series spanning decades. We investigated the potential for use of metabarcoding in CPR surveys. Samples (n = 53) were collected in two Southern Ocean transects and metazoans identified using standard microscopic methods and by high-throughput sequencing of a cytochrome c oxidase subunit I marker. DNA increased the number of metazoan species identified and provided high-resolution taxonomy of groups problematic in conventional surveys (e.g., larval echinoderms and hydrozoans). Metabarcoding also generally produced more detections than microscopy, but this sensitivity may make cross-contamination during sampling a problem. In some samples, the prevalence of DNA from large plankton such as krill masked the presence of smaller species. We investigated adding a fixed amount of exogenous DNA to samples as an internal control to allow determination of relative plankton biomass. Overall, the metabarcoding data represent a substantial shift in perspective, making direct integration into current long-term time-series challenging. We discuss a number of hurdles that exist for progressing DNA metabarcoding from the current snapshot studies to the requirements of a long-term monitoring programme. Given the power and continually increasing efficiency of metabarcoding, it is almost certain this approach will play an important role in future plankton monitoring.

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

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.

Evaluation of metagenetic community analysis of planktonic copepods using Illumina MiSeq: Comparisons with morphological classification and metagenetic analysis using Roche 454

Metagenetics is a rapid and taxonomically comprehensive method for revealing community structures within environmental samples, based on large amounts of sequence data produced by high-throughput sequencers. Because community structures of planktonic copepods are important in the ocean owing to their high diversity and abundance, a metagenetic analysis of the 28S D2 region using Roche 454 was previously developed. However, the Illumina MiSeq platform with a high sequence output is being used more frequently in metagenetics, and non-calanoid copepods have not previously been fully evaluated. Here, we evaluated an Illumina MiSeq-based metagenetic analysis using a mock community and field-collected samples that were examined in a previous study using Roche 454, and the community structure, including non-calanoid copepods, was compared among morphological and metagenetic analyses. We removed a singleton read and applied an appropriate abundance threshold to remove erroneous Molecular Operational Taxonomic Units (MOTUs) with low-abundance sequences in the MiSeq-based analysis. Results showed that the copepod community was successfully characterized using Illumina MiSeq. Higher-quality sequences were obtained using MiSeq than by Roche 454, which reduced the overestimation of diversity, especially at a strict 99% similarity threshold for MOTU clustering. Taxonomic compositions in terms of both biomass and presence/absence of species, including non-calanoids, were more appropriately represented in the MiSeq- than in Roche 454-based analysis. Our data showed that metagenetic analysis using Illumina MiSeq is more useful for revealing copepod communities than Roche 454, owing to the lower cost and higher quality.

A quantitative protocol for DNA metabarcoding of springtails (Collembola)

We developed a novel protocol with superior quantitative analysis results for DNA metabarcoding of Collembola, a major soil microarthropod order. Degenerate PCR primers were designed for conserved regions in the mitochondrial cytochrome c oxidase subunit I (mtCOI) and 16S ribosomal RNA (mt16S) genes based on published collembolan mitogenomes. The best primer pair was selected based on its ability to amplify each gene, irrespective of the species. DNA was extracted from 10 natural communities sampled in a temperate forest (with typically 25-30 collembolan species per 10 soil samples) and 10 mock communities (with seven cultured collembolan species). The two gene regions were then amplified using the selected primers, ligated with adapters for 454 technology, and sequenced. Examination of the natural community samples showed that 32 and 36 operational taxonomic units (defined at a 90% sequence similarity threshold) were recovered from the mtCOI and mt16S data, respectively, which were comparable to the results of the microscopic identification of 25 morphospecies. Further, sequence abundances for each collembolan species from the mtCOI and mt16S data of the mock communities, after normalization by using a species as the internal control, showed good correlation with the number of individuals in the samples (R = 0.91-0.99), although relative species abundances within a mock community sample estimated from sequences were skewed from community composition in terms of the number of individuals or biomass of the species. Thus, this protocol enables the comparison of collembolan communities in a quantitative manner by metabarcoding.

Effect of marker choice and thermal cycling protocol on zooplankton DNA metabarcoding studies

DNA metabarcoding is a promising approach for rapidly surveying biodiversity and is likely to become an important tool for measuring ecosystem responses to environmental change. Metabarcoding markers need sufficient taxonomic coverage to detect groups of interest, sufficient sequence divergence to resolve species, and will ideally indicate relative abundance of taxa present. We characterized zooplankton assemblages with three different metabarcoding markers (nuclear 18S rDNA, mitochondrial COI, and mitochondrial 16S rDNA) to compare their performance in terms of taxonomic coverage, taxonomic resolution, and correspondence between morphology‐ and DNA‐based identification. COI amplicons sequenced on separate runs showed that operational taxonomic units representing >0.1% of reads per sample were highly reproducible, although slightly more taxa were detected using a lower annealing temperature. Mitochondrial COI and nuclear 18S showed similar taxonomic coverage across zooplankton phyla. However, mitochondrial COI resolved up to threefold more taxa to species compared to 18S. All markers revealed similar patterns of beta‐diversity, although different taxa were identified as the greatest contributors to these patterns for 18S. For calanoid copepod families, all markers displayed a positive relationship between biomass and sequence reads, although the relationship was typically strongest for 18S. The use of COI for metabarcoding has been questioned due to lack of conserved primer‐binding sites. However, our results show the taxonomic coverage and resolution provided by degenerate COI primers, combined with a comparatively well‐developed reference sequence database, make them valuable metabarcoding markers for biodiversity assessment.

Metagenomic analysis of the marine coastal invertebrates of South Korea as assessed by Ilumina MiSeq

Metagenomic analysis was carried out for the first time on the marine coastal invertebrates of South Korea. Samples collected from coastal areas of Korea were filtered through a 63 µm mesh, then their 18S rDNA V4 regions were amplified. High-throughput sequencing of PCR amplicons using Illumina MiSeq and BLAST against the SILVA database showed that a total of 319 eukaryotic Operational Taxonomic Units were identified at the species level. Annelida, Arthropoda, Mollusca, Nematoda, and Platyhelminthes and for 92.23% of the total 103 metazoan species belonging to 101 genera, 75 families, and 10 phyla. Of these, several taxa previously unreported to exist in Korea were detected at the family level compared with the integrated database from three main Korean biodiversity DBs (MABIK, KOMBIS, and MRBR). Analysis of beta diversity of the community structure of invertebrates indicated that the composition of marine invertebrates is likely to be affected by habitat type rather than geographical distance. The present study showed that metagenomic high-throughput technology can be used to unravel species diversity and for various studies regarding marine invertebrate community structure.

Please mind the gap - Visual census and cryptic biodiversity assessment at central Red Sea coral reefs

Coral reefs harbor the most diverse assemblages in the ocean, however, a large proportion of the diversity is cryptic and, therefore, undetected by standard visual census techniques. Cryptic and exposed communities differ considerably in species composition and ecological function. This study compares three different coral reef assessment protocols: i) visual benthic reef surveys: ii) visual census of Autonomous Reef Monitoring Structures (ARMS) plates; and iii) metabarcoding techniques of the ARMS (including sessile, 106-500 μm and 500-2000 μm size fractions), that target the cryptic and exposed communities of three reefs in the central Red Sea. Visual census showed a dominance of Cnidaria (Anthozoa) and Rhodophyta on the reef substrate, while Porifera, Bryozoa and Rhodophyta were the most abundant groups on the ARMS plates. Metabarcoding, targeting the 18S rRNA gene, significantly increased estimates of the species diversity (p < 0.001); revealing that Annelida were generally the dominant phyla (in terms of reads) of all fractions and reefs. Furthermore, metabarcoding detected microbial eukaryotic groups such as Syndiniophyceae, Mamiellophyceae and Bacillariophyceae as relevant components of the sessile fraction. ANOSIM analysis showed that the three reef sites showed no differences based on the visual census data. Metabarcoding showed a higher sensitivity for identifying differences between reef communities at smaller geographic scales than standard visual census techniques as significant differences in the assemblages were observed amongst the reefs. Comparison of the techniques showed no similar patterns for the visual techniques while the metabarcoding of the ARMS showed similar patterns amongst fractions. Establishing ARMS as a standard tool in reef monitoring will not only advance our understanding of local processes and ecological community response to environmental changes, as different faunal components will provide complementary information but also improve the estimates of biodiversity in coral reef benthic communities. This study lays the foundations for further studies looking at integrating traditional reef survey methodologies with complementary approaches, such as metabarcoding, which investigate other components of the reef community.

Response of Spring Diatoms to CO2 Availability in the Western North Pacific as Determined by Next-Generation Sequencing

Next-generation sequencing (NGS) technologies have enabled us to determine phytoplankton community compositions at high resolution. However, few studies have adopted this approach to assess the responses of natural phytoplankton communities to environmental change. Here, we report the impact of different CO₂ levels on spring diatoms in the Oyashio region of the western North Pacific as estimated by NGS of the diatom-specific rbcL gene (DNA), which encodes the large subunit of RubisCO. We also examined the abundance and composition of rbcL transcripts (cDNA) in diatoms to assess their physiological responses to changing CO₂ levels. A short-term (3-day) incubation experiment was carried out on-deck using surface Oyashio waters under different pCO₂ levels (180, 350, 750, and 1000 μatm) in May 2011. During the incubation, the transcript abundance of the diatom-specific rbcL gene decreased with an increase in seawater pCO₂ levels. These results suggest that CO₂ fixation capacity of diatoms decreased rapidly under elevated CO₂ levels. In the high CO₂ treatments (750 and 1000 μatm), diversity of diatom-specific rbcL gene and its transcripts decreased relative to the control treatment (350 μatm), as well as contributions of Chaetocerataceae, Thalassiosiraceae, and Fragilariaceae to the total population, but the contributions of Bacillariaceae increased. In the low CO₂ treatment, contributions of Bacillariaceae also increased together with other eukaryotes. These suggest that changes in CO₂ levels can alter the community composition of spring diatoms in the Oyashio region. Overall, the NGS technology provided us a deeper understanding of the response of diatoms to changes in CO₂ levels in terms of their community composition, diversity, and photosynthetic physiology.

In silico assessment of primers for eDNA studies using PrimerTree and application to characterize the biodiversity surrounding the Cuyahoga River

Analysis of environmental DNA (eDNA) enables the detection of species of interest from water and soil samples, typically using species-specific PCR. Here, we describe a method to characterize the biodiversity of a given environment by amplifying eDNA using primer pairs targeting a wide range of taxa and high-throughput sequencing for species identification. We tested this approach on 91 water samples of 40 mL collected along the Cuyahoga River (Ohio, USA). We amplified eDNA using 12 primer pairs targeting mammals, fish, amphibians, birds, bryophytes, arthropods, copepods, plants and several microorganism taxa and sequenced all PCR products simultaneously by high-throughput sequencing. Overall, we identified DNA sequences from 15 species of fish, 17 species of mammals, 8 species of birds, 15 species of arthropods, one turtle and one salamander. Interestingly, in addition to aquatic and semi-aquatic animals, we identified DNA from terrestrial species that live near the Cuyahoga River. We also identified DNA from one Asian carp species invasive to the Great Lakes but that had not been previously reported in the Cuyahoga River. Our study shows that analysis of eDNA extracted from small water samples using wide-range PCR amplification combined with high-throughput sequencing can provide a broad perspective on biological diversity.

18S rRNA V9 metabarcoding for diet characterization: a critical evaluation with two sympatric zooplanktivorous fish species

The potential of the 18S rRNA V9 metabarcoding approach for diet assessment was explored using MiSeq paired‐end (PE; 2 × 150 bp) technology. To critically evaluate the method′s performance with degraded/digested DNA, the diets of two zooplanktivorous fish species from the Bay of Biscay, European sardine (Sardina pilchardus) and European sprat (Sprattus sprattus), were analysed. The taxonomic resolution and quantitative potential of the 18S V9 metabarcoding was first assessed both in silico and with mock and field plankton samples. Our method was capable of discriminating species within the reference database in a reliable way providing there was at least one variable position in the 18S V9 region. Furthermore, it successfully discriminated diet between both fish species, including habitat and diel differences among sardines, overcoming some of the limitations of traditional visual‐based diet analysis methods. The high sensitivity and semi‐quantitative nature of the 18S V9 metabarcoding approach was supported by both visual microscopy and qPCR‐based results. This molecular approach provides an alternative cost and time effective tool for food‐web analysis.

Assessment of Zooplankton Community Composition along a Depth Profile in the Central Red Sea

The composition of zooplankton in the water column has received limited attention in the main body of the Red Sea and this study investigates the change in the community both spatially and temporally across 11 stations in the central Red Sea. Using molecular methods to target the v9 region of the 18S rRNA gene a total of approximately 11.5 million reads were sequenced resulting in 2528 operational taxonomic units (OTUs) at 97% similarity. The phylum Arthropoda dominated in terms of reads accounting for on average 86.2% and 65.3% for neuston nets and vertical multinets respectively. A reduction in the number of OTUs was noticed with depth for both total metazoa and Maxillopoda whilst there was also a significant change in the composition of the Maxillopoda community. The genus Corycaeus had a higher proportion of reads in the epipelagic zone with Pleuromamma becoming increasingly dominant with depth. No significant difference was observed in the community between night and day sampling however there was a significant difference in the zooplankton community between two sampling periods separated by 10 days.

Metagenetic tools for the census of marine meiofaunal biodiversity: An overview

Marine organisms belonging to meiofauna (size range: 20-500 μm) are amongst the most abundant and highly diversified metazoans on Earth including 22 over 35 known animal Phyla and accounting for more than 2/3 of the abundance of metazoan organisms. In any marine system, meiofauna play a key role in the functioning of the food webs and sustain important ecological processes. Estimates of meiofaunal biodiversity have been so far almost exclusively based on morphological analyses, but the very small size of these organisms and, in some cases, the insufficient morphological distinctive features limit considerably the census of the biodiversity of this component. Molecular approaches recently applied also to small invertebrates (including meiofauna) can offer a new momentum for the census of meiofaunal biodiversity. Here, we provide an overview on the application of metagenetic approaches based on the use of next generation sequencing platforms to study meiofaunal biodiversity, with a special focus on marine nematodes. Our overview shows that, although such approaches can represent a useful tool for the census of meiofaunal biodiversity, there are still different shortcomings and pitfalls that prevent their extensive use without the support of the classical taxonomic identification. Future investigations are needed to address these problems and to provide a good match between the contrasting findings emerging from classical taxonomic and molecular/bioinformatic tools.