Censusing marine eukaryotic diversity in the twenty-first century

A comparison of two gene regions for assessing community composition of eukaryotic marine microalgae from coastal ecosystems

Two gene regions commonly used to characterise the diversity of eukaryotic communities using metabarcoding are the 18S ribosomal DNA V4 and V9 gene regions. We assessed the effectiveness of these two regions for characterising diverisity of coastal eukaryotic microalgae communities (EMCs) from tropical and temperate sites. We binned amplicon sequence variants (ASVs) into the high level taxonomic groups: dinoflagellates, pennate diatoms, radial centric diatoms, polar centric diatoms, chlorophytes, haptophytes and 'other microalgae'. When V4 and V9 generated ASV abundances were compared, the V9 region generated a higher number of raw reads, captured more diversity from all high level taxonomic groups and was more closely aligned with the community composition determined using light microscopy. The V4 region did resolve more ASVs to a deeper taxonomic resolution within the dinoflagellates, but did not effectively resolve other major taxonomic divisions. When characterising these communities via metabarcoding, the use of multiple gene regions is recommended, but the V9 gene region can be used in isolation to provide high-level community biodiversity to reflect relative abundances within groups. This approach reduces the cost of sequencing multiple gene regions whilst still providing important baseline ecosystem function information.

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.

VESPA: an optimized protocol for accurate metabarcoding-based characterization of vertebrate eukaryotic endosymbiont and parasite assemblages

Protocols for characterizing taxonomic assemblages by deep sequencing of short DNA barcode regions (metabarcoding) have revolutionized our understanding of microbial communities and are standardized for bacteria, archaea, and fungi. Unfortunately, comparable methods for host-associated eukaryotes have lagged due to technical challenges. Despite 54 published studies, issues remain with primer complementarity, off-target amplification, and lack of external validation. Here, we present VESPA (Vertebrate Eukaryotic endoSymbiont and Parasite Analysis) primers and optimized metabarcoding protocol for host-associated eukaryotic community analysis. Using in silico prediction, panel PCR, engineered mock community standards, and clinical samples, we demonstrate VESPA to be more effective at resolving host-associated eukaryotic assemblages than previously published methods and to minimize off-target amplification. When applied to human and non-human primate samples, VESPA enables reconstruction of host-associated eukaryotic endosymbiont communities more accurately and at finer taxonomic resolution than microscopy. VESPA has the potential to advance basic and translational science on vertebrate eukaryotic endosymbiont communities, similar to achievements made for bacterial, archaeal, and fungal microbiomes.

Interrelationships and biogeography of the New World pufferfish genus Sphoeroides (Tetraodontiformes: Tetraodontidae) inferred using ultra-conserved DNA elements

Colonization of the New World by marine taxa has been hypothesized to have occurred through the Tethys Sea or by crossing the East Pacific Barrier. To better understand patterns and timing of diversification, geological events can be coupled with time calibrated phylogenetic hypotheses to infer major drivers of diversification. Phylogenetic relationships among members of Sphoeroides, a genus of four toothed pufferfishes (Tetraodontiformes: Tetraodontidae) which are found nearly exclusively in the New World (eastern Pacific and western Atlantic), were reconstructed using sequences from ultra-conserved DNA elements, nuclear markers with clear homology among many vertebrate taxa. Hypotheses derived from concatenated maximum-likelihood and species tree summary methods support a paraphyletic Sphoeroides, with Colomesus deeply nested within the genus. Analyses also revealed S. pachygaster, a pelagic species with a cosmopolitan distribution, as the sister taxon to the remainder of Sphoeroides and recovered distinct lineages within S. pachygaster, indicating that this cosmopolitan species may represent a species complex. Ancestral range reconstruction may suggest the genus colonized the New World through the eastern Pacific before diversifying in the western Atlantic, though date estimates for these events are uncertain due to the lack of reliable fossil record for the genus.

How we study cryptic species and their biological implications: A case study from marine shelled gastropods

Methodological and biological considerations are intertwined when studying cryptic species. A potentially large component of modern biodiversity, the frequency of cryptic species among taxonomic groups is not well documented. The term "cryptic species" is imprecisely used in scientific literature, causing ambiguity when interpreting their evolutionary and ecological significance. This study reviews how cryptic species have been defined, discussing implications for taxonomy and biology, and explores these implications with a case study based on recently published literature on extant shelled marine gastropods. Reviewed gastropods were recorded by species. Records of cryptic gastropods were presented by authors with variable levels of confidence but were difficult to disentangle from inherent biases in the study effort. These complexities notwithstanding, most gastropod species discussed were not cryptic. To the degree that this review's sample represents extinct taxa, the results suggest that a high proportion of shelled marine gastropod species are identifiable for study in the fossil record. Much additional work is needed to provide a more adequate understanding of the relative frequency of cryptic species in shelled marine gastropods, which should start with more explicit definitions and targeted case studies.

Community dynamics and co-occurrence relationships of pelagic ciliates and their potential prey at a coastal and an offshore station in the ultra-oligotrophic Eastern Mediterranean Sea

Ciliates have been recognized as one of the major components of the microbial food web, especially in ultra-oligotrophic waters, such as the Eastern Mediterranean Sea, where nutrients are scarce and the microbial community is dominated by pico- and nano-sized organisms. For this reason, ciliates play an important role in these ecosystems since they are the main planktonic grazers. Regardless the importance of these organisms, little is known about the community structure of heterotrophic and mixotrophic ciliates and how they are associated to their potential prey. In this study, we used 18S V4 rRNA gene metabarcoding to analyze ciliate community dynamics and how the relationship with potential prey changes according to different seasons and depths. Samples were collected seasonally at two stations of the Eastern Mediterranean Sea (HCB: coastal, M3A: offshore) from the surface and deep chlorophyll maximum (DCM) layers. The ciliate community structure varied across depths in HCB and across seasons in M3A, and the network analysis showed that in both stations, mixotrophic oligotrichs were positively associated with diatoms and showed few negative associations with ASVs annotated as marine Stramenopiles (MAST). On the other hand, heterotrophic tintinnids showed negative relationships in both HCB and M3A stations, mostly with Ochrophyta and Chlorophyta. These results showed, in first place that, although the two stations are close to each other, the ciliate dynamics differed between them. Moreover, mixotrophic and heterotrophic ciliates may have different ecological niches since mixotrophic ciliates may be more selective compared to heterotrophic species regarding their prey. These findings are the first glimpse into an understanding of the dynamics between heterotrophic and mixotrophic ciliates and their role in microbial assemblages and dynamics of ultra-oligotrophic environments.

High taxonomic diversity and miniaturization in benthic communities under persistent natural CO2 disturbances

Metabarcoding techniques have revolutionized ecological research in recent years, facilitating the differentiation of cryptic species and revealing previously hidden diversity. In the current scenario of climate change and ocean acidification, biodiversity loss is one of the main threats to marine ecosystems. Here, we explored the effects of ocean acidification on marine benthic communities using DNA metabarcoding to assess the diversity of algae and metazoans. Specifically, we examined the natural pH gradient generated by the Fuencaliente CO2 vent system, located near La Palma Island (Canary Islands). High-resolution COI metabarcoding analyses revealed high levels of taxonomic diversity in an acidified natural area for the first time. This high number of species arises from the detection of small and cryptic species that were previously undetectable by other techniques. Such species are apparently tolerant to the acidification levels expected in future oceans. Hence and following our results, future subtropical communities are expected to keep high biodiversity values under an acidification scenario, although they will tend toward overall miniaturization due to the dominance of small algal and invertebrate species, leading to changes in ecosystem functions.

Advanced approaches for the diagnosis and chemoprevention of canine vector-borne pathogens and parasites-Implications for the Asia-Pacific region and beyond

Vector-borne pathogens (VBPs) of canines are a diverse range of infectious agents, including viruses, bacteria, protozoa and multicellular parasites, that are pernicious and potentially lethal to their hosts. Dogs across the globe are afflicted by canine VBPs, but the range of different ectoparasites and the VBPs that they transmit predominate in tropical regions. Countries within the Asia-Pacific have had limited prior research dedicated to exploring the epidemiology of canine VBPs, whilst the few studies that have been conducted show VBP prevalence to be high, with significant impacts on dog health. Moreover, such impacts are not restricted to dogs, as some canine VBPs are zoonotic. We reviewed the status of canine VBPs in the Asia-Pacific, with particular focus on nations in the tropics, whilst also investigating the history of VBP diagnosis and examining recent progress in the field, including advanced molecular methods, such as next-generation sequencing (NGS). These tools are rapidly changing the way parasites are detected and discovered, demonstrating a sensitivity equal to, or exceeding that of, conventional molecular diagnostics. We also provide a background to the armoury of chemopreventive products available for protecting dogs from VBP. Here, field-based research within high VBP pressure environments has underscored the importance of ectoparasiticide mode of action on their overall efficacy. The future of canine VBP diagnosis and prevention at a global level is also explored, highlighting how evolving portable sequencing technologies may permit diagnosis at point-of-care, whilst further research into chemopreventives will be essential if VBP transmission is to be effectively controlled.

Evaluating biodiversity for coral reef reformation and monitoring on complex 3D structures using environmental DNA (eDNA) metabarcoding

Quantifying coral reef biodiversity is challenging for cryptofauna and organisms in early life stages. We demonstrate the utility of eDNA metabarcoding as a tool for comprehensively evaluating invertebrate communities on complex 3D structures for reef reformation, and the role these structures play in provisioning habitat for organisms. 3D design and printing were used to create 18 complex tiles, which were used to form artificial reef structures. eDNA was collected from scraping tile surfaces for organismal biomass and from seawater samples around the artificial reefs in the Gulf of Eilat/Aqaba, Red Sea. Metabarcoding targeted the mitochondrial COI gene with specific primers for marine biodiversity. We provide the first eDNA biodiversity baseline for the Gulf of Eilat/Aqaba, capturing extensive information on species abundance, richness, and diversity. Tile tops had higher phylogenetic diversity and richness, despite a higher abundance of organisms on tile bottoms, highlighting the detection of cryptic organisms with eDNA. We recommend eDNA metabarcoding for reef restoration initiatives, especially for complex marine structures, to improve success and evaluation of biodiversity.

Benthic invertebrates in Svalbard fjords-when metabarcoding does not outperform traditional biodiversity assessment

To protect and restore ecosystems and biodiversity is one of the 10 challenges identified by the United Nations's Decade of the Ocean Science. In this study we used eDNA from sediments collected in two fjords of the Svalbard archipelago and compared the taxonomic composition with traditional methods through metabarcoding, targeting mitochondrial CO1, to survey benthos. Clustering of 21.6 mill sequence reads with a d value of 13 in swarm, returned about 25 K OTU reads. An identification search with the BOLD database returned 12,000 taxonomy annotated sequences spanning a similarity range of 50% to 100%. Using an acceptance filter of minimum 90% similarity to the CO1 reference sequence, we found that 74% of the ca 100 taxon identified sequence reads were Polychaeta and 22% Nematoda. Relatively few other benthic invertebrate species were detected. Many of the identified sequence reads were extra-organismal DNA from terrestrial, planktonic, and photic zone sources. For the species rich Polychaeta, we found that, on average, only 20.6% of the species identified from morphology were also detected with DNA. This discrepancy was not due to missing reference sequences in the search database, because 90-100% (mean 96.7%) of the visually identified species at each station were represented with barcodes in Boldsystems. The volume of DNA samples is small compared with the volume searched in visual sorting, and the replicate DNA-samples in sum covered only about 2% of the surface area of a grab. This may considerably reduce the detection rate of species that are not uniformly distributed in the sediments. Along with PCR amplification bias and primer mismatch, this may be an important reason for the limited congruence of species identified with the two approaches. However, metabarcoding also identified 69 additional species that are usually overlooked in visual sample sorting, demonstrating how metabarcoding can complement traditional methodology by detecting additional, less conspicuous groups of organisms.

Profiling Analysis of Filter Feeder Polypedilum (Chironomidae) Gut Contents Using eDNA Metabarcoding Following Contrasting Habitat Types-Weir and Stream

We analyzed the dietary composition of Polypedilum larvae among two contrasting habitats (river and weir). Our approach was (i) to apply eDNA-based sampling to reveal the gut content of the chironomid larvae, (ii) the diversity of gut contents in the two aquatic habitats, and (iii) assessment of habitat sediment condition with the food sources in the gut. The most abundant food was Chlorophyta in the gut of the river (20%) and weir (39%) chironomids. The average ratio of fungi, protozoa, and zooplankton in river chironomids gut was 5.9%, 7.2%, and 3.8%, while it was found decreased to 1.2%, 2.5%, and 0.1% in weir chironomids. Aerobic fungi in river midge guts were 3.6% and 10.34% in SC and IS, while they were in the range of 0.34-2.58% in weir midges. The hierarchical clustering analysis showed a relationship of environmental factors with food contents. Abiotic factors (e.g., pH) in the river and weir habitats correlated the clustered pattern with phytoplankton and minor groups of fungi. This study could help understand the food source diversity in the chironomid and habitat environmental conditions by using eDNA metabarcoding as an effective tool to determine dietary composition.

Life barcoded by DNA barcodes

The modern concept of DNA-based barcoding for cataloguing biodiversity was proposed in 2003 by first adopting an approximately 600 bp fragment of the mitochondrial COI gene to compare via nucleotide alignments with known sequences from specimens previously identified by taxonomists. Other standardized regions meeting barcoding criteria then are also evolving as DNA barcodes for fast, reliable and inexpensive assessment of species composition across all forms of life, including animals, plants, fungi, bacteria and other microorganisms. Consequently, global DNA barcoding campaigns have resulted in the formation of many online workbenches and databases, such as BOLD system, as barcode references, and facilitated the development of mini-barcodes and metabarcoding strategies as important extensions of barcode techniques. Here we intend to give an overview of the characteristics and features of these barcode markers and major reference libraries existing for barcoding the planet’s life, as well as to address the limitations and opportunities of DNA barcodes to an increasingly broader community of science and society.

Combining multi-marker metabarcoding and digital holography to describe eukaryotic plankton across the Newfoundland Shelf

The planktonic diversity throughout the oceans is vital to ecosystem functioning and linked to environmental change. Plankton monitoring tools have advanced considerably with high-throughput in-situ digital cameras and genomic sequencing, opening new challenges for high-frequency observations of community composition, structure, and species discovery. Here, we combine multi-marker metabarcoding based on nuclear 18S (V4) and plastidial 16S (V4–V5) rRNA gene amplicons with a digital in-line holographic microscope to provide a synoptic diversity survey of eukaryotic plankton along the Newfoundland Shelf (Canada) during the winter transition phase of the North Atlantic bloom phenomenon. Metabarcoding revealed a rich eukaryotic diversity unidentifiable in the imaging samples, confirming the presence of ecologically important saprophytic protists which were unclassifiable in matching images, and detecting important groups unobserved or taxonomically unresolved during similar sequencing campaigns in the Northwest Atlantic Ocean. In turn, imaging analysis provided quantitative observations of widely prevalent plankton from every trophic level. Despite contrasting plankton compositions portrayed by each sampling method, both capture broad spatial differences between the northern and southern sectors of the Newfoundland Shelf and suggest complementary estimations of important features in eukaryotic assemblages. Future tasks will involve standardizing digital imaging and metabarcoding for wider use and consistent, comparable ocean observations.

Ecosystem-specific microbiota and microbiome databases in the era of big data

The rapid development of sequencing methods over the past decades has accelerated both the potential scope and depth of microbiota and microbiome studies. Recent developments in the field have been marked by an expansion away from purely categorical studies towards a greater investigation of community functionality. As in-depth genomic and environmental coverage is often distributed unequally across major taxa and ecosystems, it can be difficult to identify or substantiate relationships within microbial communities. Generic databases containing datasets from diverse ecosystems have opened a new era of data accessibility despite costs in terms of data quality and heterogeneity. This challenge is readily embodied in the integration of meta-omics data alongside habitat-specific standards which help contextualise datasets both in terms of sample processing and background within the ecosystem. A special case of large genomic repositories, ecosystem-specific databases (ES-DB's), have emerged to consolidate and better standardise sample processing and analysis protocols around individual ecosystems under study, allowing independent studies to produce comparable datasets. Here, we provide a comprehensive review of this emerging tool for microbial community analysis in relation to current trends in the field. We focus on the factors leading to the formation of ES-DB's, their comparison to traditional microbial databases, the potential for ES-DB integration with meta-omics platforms, as well as inherent limitations in the applicability of ES-DB's.

Adapting mark-recapture methods to estimating accepted species-level diversity: a case study with terrestrial Gastropoda

We introduce a new method of estimating accepted species diversity by adapting mark-recapture methods to comparisons of taxonomic databases. A taxonomic database should become more complete over time, so the error bar on an estimate of its completeness and the known diversity of the taxon it treats will decrease. Independent databases can be correlated, so we use the time course of estimates comparing them to understand the effect of correlation. If a later estimate is significantly larger than an earlier one, the databases are positively correlated, if it is significantly smaller, they are negatively correlated, and if the estimate remains roughly constant, then the correlations have averaged out. We tested this method by estimating how complete MolluscaBase is for accepted names of terrestrial gastropods. Using random samples of names from an independent database, we determined whether each name led to a name accepted in MolluscaBase. A sample tested in August 2020 found that 16.7% of tested names were missing; one in July 2021 found 5.3% missing. MolluscaBase grew by almost 3,000 accepted species during this period, reaching 27,050 species. The estimates ranged from 28,409 ± 365 in 2021 to 29,063 ± 771 in 2020. All estimates had overlapping 95% confidence intervals, indicating that correlations between the databases did not cause significant problems. Uncertainty beyond sampling error added 475 ± 430 species, so our estimate for accepted terrestrial gastropods species at the end of 2021 is 28,895 ± 630 species. This estimate is more than 4,000 species higher than previous ones. The estimate does not account for ongoing flux of species into and out of synonymy, new discoveries, or changing taxonomic methods and concepts. The species naming curve for terrestrial gastropods is still far from reaching an asymptote, and combined with the additional uncertainties, this means that predicting how many more species might ultimately be recognized is presently not feasible. Our methods can be applied to estimate the total number of names of Recent mollusks (as opposed to names currently accepted), the known diversity of fossil mollusks, and known diversity in other phyla.

Do pseudogenes pose a problem for metabarcoding marine animal communities?

Because DNA metabarcoding typically employs sequence diversity among mitochondrial amplicons to estimate species composition, nuclear mitochondrial pseudogenes (NUMTs) can inflate diversity. This study quantifies the incidence and attributes of NUMTs derived from the 658 bp barcode region of cytochrome c oxidase I (COI) in 156 marine animal genomes. NUMTs were examined to ascertain if they could be recognized by their possession of indels or stop codons. In total, 309 NUMTs ≥ 150 bp were detected, with an average of 1.98 per species (range = 0-33) and a mean length of 391 bp ± 200 bp. Among this total, 75 (24.3 %) lacked indels or stop codons. NUMTs appear to pose the greatest interpretational risk when short (< 313 bp) amplicons are used, such as in eDNA studies, dietary analyses, or processed fish identification. Employing the standard amplicon length (313 bp) for marine metabarcoding, NUMTs could potentially inflate the OTU count by 21% above the true species count while also raising intraspecific variation at COI by 15%. However, when both amplicon length and position are considered, inflation in OTU counts and in barcode variation were just 9% and 10%, respectively, suggesting NUMTs will not seriously distort biodiversity assessments. There was a weak positive correlation between genome size and NUMT count but no variation among phyla or trophic groups. Until bioinformatic advances improve NUMT detection, the best defense involves targeting long amplicons and developing reference databases that include both mitochondrial sequences and their NUMT derivatives.

It’s time for a new type of type to facilitate naming the microbial world

Since January 1, 2001, the only acceptable nomenclatural type for species under the International Code of Nomenclature of Prokaryotes (ICNP) has been pure cultures. Here, we argue that this requirement is discordant with the more inclusive nature of nomenclatural types accepted under other codes of nomenclature and posit that the unique rigidity of the ICNP has failed to serve the broad research community and has stifled progress. This case is based on the axiom that many archaea and bacteria are interdependent in nature and therefore difficult, if not impossible, to grow, preserve, and distribute as pure cultures. As such, a large proportion of Earth's biodiversity cannot be named under the current system, which limits our ability to communicate about microbial diversity within and beyond the microbiology research community. Genome sequence data are now encouraged for valid publication of new taxa in microbial systematics journals, and metagenome-assembled genomes and single cell-amplified genomes are being generated rapidly from every biome on Earth. Thus, genome sequences are available for both cultivated and uncultivated microorganisms and can readily serve as a new category of nomenclatural type, allowing for a unified nomenclature for all archaea and bacteria, whether or not they are available as pure cultures. Ideally this would be under a single code of nomenclature but, as we review here, the newly established SeqCode will operate in parallel with the ICNP as a first step toward this goal.

Functional repertoire convergence of distantly related eukaryotic plankton lineages abundant in the sunlit ocean

Marine planktonic eukaryotes play critical roles in global biogeochemical cycles and climate. However, their poor representation in culture collections limits our understanding of the evolutionary history and genomic underpinnings of planktonic ecosystems. Here, we used 280 billion Tara Oceans metagenomic reads from polar, temperate, and tropical sunlit oceans to reconstruct and manually curate more than 700 abundant and widespread eukaryotic environmental genomes ranging from 10 Mbp to 1.3 Gbp. This genomic resource covers a wide range of poorly characterized eukaryotic lineages that complement long-standing contributions from culture collections while better representing plankton in the upper layer of the oceans. We performed the first, to our knowledge, comprehensive genome-wide functional classification of abundant unicellular eukaryotic plankton, revealing four major groups connecting distantly related lineages. Neither trophic modes of plankton nor its vertical evolutionary history could completely explain the functional repertoire convergence of major eukaryotic lineages that coexisted within oceanic currents for millions of years.

High-Resolution Screening for Marine Prokaryotes and Eukaryotes With Selective Preference for Polyethylene and Polyethylene Terephthalate Surfaces

Marine plastic debris serve as substrates for the colonization of a variety of prokaryote and eukaryote organisms. Of particular interest are the microorganisms that have adapted to thrive on plastic as they may contain genes, enzymes or pathways involved in the adhesion or metabolism of plastics. We implemented DNA metabarcoding with nanopore MinION sequencing to compare the 1-month-old biomes of hydrolyzable (polyethylene terephthalate) and non-hydrolyzable (polyethylene) plastics surfaces vs. those of glass and the surrounding water in a Mediterranean Sea marina. We sequenced longer 16S rRNA, 18S rRNA, and ITS barcode loci for a more comprehensive taxonomic profiling of the bacterial, protist, and fungal communities, respectively. Long read sequencing enabled high-resolution mapping to genera and species. Using previously established methods we performed differential abundance screening and identified 30 bacteria and five eukaryotic species, that were differentially abundant on plastic compared to glass. This approach will allow future studies to characterize the plastisphere communities and to screen for microorganisms with a plastic-metabolism potential.

18S rDNA amplicon sequence data (V1–V3) of the Palmyra Atoll National Wildlife Refuge, Central Pacific

To address the global biodiversity crisis, standardized data that are rapidly obtainable through minimally invasive means are needed for documenting change and informing conservation within threatened and diverse systems, such as coral reefs. In this data paper, we describe 18S rRNA gene amplicon data (V1–V3 region) generated from samples collected to begin characterizing coral reef eukaryotic community composition at the Palmyra Atoll National Wildlife Refuge in the Central Pacific Ocean. Sixteen samples were obtained across four sample types: sediments from two sieved fractions (100–500 μm, n = 3; 500 μm-2 mm, n = 3) and sessile material scrapings (n = 3) from Autonomous Reef Monitoring Structures (ARMS) sampled in 2015, as well as seawater from 2012 (n = 7). After filtering and contaminant removal, 3,861 Amplicon Sequence Variants (ASVs) were produced from 1,062,238 reads. The rarefaction curves demonstrated adequate sampling depth, and communities grouped by sample type. The dominant orders across samples were polychaete worms (Eunicida), demosponges (Poecilosclerida), and bryozoans (Cheilostomatida). The ten most common orders in terms of relative abundance comprised ~60% of all sequences and 23% of ASVs, and included reef-building crustose coralline algae (CCA; Corallinophycidae) and stony corals (Scleractinia), two taxa associated with healthy reefs. Highlighting the need for further study, ~21% of the ASVs were identified as uncultured, incertae sedis, or not assigned to phylum or order. This data paper presents the first 18S rDNA survey at Palmyra Atoll and serves as a baseline for biodiversity assessment, monitoring, and conservation of this remote and pristine ecosystem.

Seasonal Diversity of Microeukaryotes in the Han River, Korea Through 18S rRNA Gene Metabarcoding

Freshwater ecosystems contain a large diversity of microeukaryotes that play important roles in maintaining their structure. Microeukaryote communities vary in composition and abundance on the basis of temporal and environmental variables and may serve as useful bioindicators of environmental changes. In the present study, 18S rRNA metabarcoding was employed to investigate the seasonal diversity of microeukaryote communities during four seasons in the Han River, Korea. In total, 882 unique operational taxonomic units (OTUs) were detected, including various diatoms, metazoans (e.g., arthropods and rotifers), chlorophytes, and fungi. Although alpha diversity revealed insignificant differences based on seasons, beta diversity exhibited a prominent variation in the community composition as per seasons. The analysis revealed that the diversity of microeukaryotes was primarily driven by seasonal changes in the prevailing conditions of environmental water temperature and dissolved oxygen. Moreover, potential indicator OTUs belonging to diatoms and chlorophytes were associated with seasonal and environmental factors. This analysis was a preliminary study that established a continuous monitoring system using metabarcoding. This approach could be an effective tool to manage the Han River along with other freshwater ecosystems in Korea.

DNA barcoding and morphological identification of spiny lobsters in South Korean waters: a new record of Panulirus longipes and Panulirus homarus homarus

To date, 19 species of spiny lobsters from the genus Panulirus have been discovered, of which only P. japonicus, P. penicilatus, P. stimpsoni, and P. versicolor have been documented in South Korean waters. In this study, we aimed to identify and update the current list of spiny lobster species that inhabit South Korean waters based on the morphological features and the phylogenetic profile of cytochrome oxidase I (COI) of mitochondrial DNA (mtDNA). Spiny lobsters were collected from the southern and eastern coasts of Jeju Island, South Korea. Phylogenetic analyses were performed using neighbor-joining (NJ), maximum likelihood (ML), and Bayesian inference (BI) methods. The ML tree was used to determine the spiny lobster lineages, thereby clustering the 17 specimens collected in this study into clades A, B, C, and D, which were reciprocally monophyletic with P. japonicus, P. homarus homarus, P. longipes, and P. stimpsoni, respectively. These clades were also supported by morphological examinations. Interestingly, morphological variations, including the connected pleural and transverse groove at the third abdominal somite, were observed in four specimens that were genetically confirmed as P. japonicus. This finding is novel within the P. japonicus taxonomical reports. Additionally, this study updates the documentation of spiny lobsters inhabiting South Korean waters as P. longipes and P. homarus homarus were recorded for the first time in this region.

Comparing species detection success between molecular markers in DNA metabarcoding of coastal macroinvertebrates

DNA metabarcoding has great potential to improve marine biomonitoring programs by providing a rapid and accurate assessment of species composition in zoobenthic communities. However, some methodological improvements are still required, especially regarding failed detections, primers efficiency and incompleteness of databases. Here we assessed the efficiency of two different marker loci (COI and 18S) and three primer pairs in marine species detection through DNA metabarcoding of the macrozoobenthic communities colonizing three types of artificial substrates (slate, PVC and granite), sampled between 3 and 15 months of deployment. To accurately compare detection success between markers, we also compared the representativeness of the detected species in public databases and revised the reliability of the taxonomic assignments. Globally, we recorded extensive complementarity in the species detected by each marker, with 69% of the species exclusively detected by either 18S or COI. Individually, each of the three primer pairs recovered, at most, 52% of all species detected on the samples, showing also different abilities to amplify specific taxonomic groups. Most of the detected species have reliable reference sequences in their respective databases (82% for COI and 72% for 18S), meaning that when a species was detected by one marker and not by the other, it was most likely due to faulty amplification, and not by lack of matching sequences in the database. Overall, results showed the impact of marker and primer applied on species detection ability and indicated that, currently, if only a single marker or primer pair is employed in marine zoobenthos metabarcoding, a fair portion of the diversity may be overlooked.

Seasonal and spatial patterns of eukaryotic phytoplankton communities in an urban river based on marker gene

The seasonal and spatial eukaryotic phytoplankton composition in the Fenhe River was investigated based on the 18S rDNA V4 region. The relationship between phytoplankton functional groups and environmental factors was explored to effectively capture the responses of these taxa to environmental gradients and their effects on ecosystem function. Our results indicated that the Chlorophyta and Bacillariophyta had higher relative abundance than other taxa, and their diversity and richness indices in spring were higher than those in other seasons. The linear discriminant analysis effect size (LEfSe) analyses detected that the potential seasonal biomarkers included Desmodesmus, Cyclotella, Pseudoschroederia, Discostella, Scenedesmus, Monoraphidium, and Nannochloropsis; the spatial biomarkers included Amphora, Neochloris, Hindakia, Pseudomuriella, Coccomyxa, Chloroidium, Scherffelia, Chromochloris, and Scotinosphaera. The systemic evolution and distribution characteristics of the first 50 representative sequences showed that the dominant genus included Desmodesmus in spring, Pseudopediastrum in summer, Mychonastes in autumn, and Monoraphidium in winter. Main seasonal variation of phytoplankton functional groups was as follows: spring (J + F + C + X1) → summer (J + F + X1 + X2) → autumn (J + F + X1 + C) → winter (X1 + J + B + X2). Pearson correlation, redundancy analysis, and variance partitioning analysis showed temperature and phosphate were the determining factors causing the changes of phytoplankton functional groups and community composition in the Fenhe River.

High-Throughput Sequencing of Environmental DNA as a Tool for Monitoring Eukaryotic Communities and Potential Pathogens in a Coastal Upwelling Ecosystem

The marine environment includes diverse microeukaryotic organisms that play important functional roles in the ecosystem. With molecular approaches, eukaryotic taxonomy has been improved, complementing classical analysis. In this study, DNA metabarcoding was performed to describe putative pathogenic eukaryotic microorganisms in sediment and marine water fractions collected in Galicia (NW Spain) from 2016 to 2018. The composition of eukaryotic communities was distinct between sediment and water fractions. Protists were the most diverse group, with the clade TSAR (Stramenopiles, Alveolata, Rhizaria, and Telonemida) as the primary representative organisms in the environment. Harmful algae and invasive species were frequently detected. Potential pathogens, invasive pathogenic organisms as well as the causative agents of harmful phytoplanktonic blooms were identified in this marine ecosystem. Most of the identified pathogens have a crucial impact on the aquacultural sector or affect to relevant species in the marine ecosystem, such as diatoms. Moreover, pathogens with medical and veterinary importance worldwide were also found, as well as pathogens that affect diatoms. The evaluation of the health of a marine ecosystem that directly affects the aquacultural sector with a zoonotic concern was performed with the metabarcoding assay.

Quantification of marine benthic communities with metabarcoding

DNA metabarcoding methods have been implemented in studies aimed at detecting and quantifying marine benthic biodiversity. In such surveys, universal barcodes are amplified and sequenced from environmental DNA. To quantify biodiversity with DNA metabarcoding, a relation between the number of DNA sequences of a species and its biomass and/or the abundance is required. However, this relationship is complicated by many factors, and it is often unknown. In this study, we validate estimates of biomass and abundance from molecular approaches with those from the traditional morphological approach. Abundance and biomass were quantified from 126 samples of benthic intertidal mudflat using traditional morphological approaches and compared with frequency of occurrence and relative read abundance estimates from a molecular approach. A relationship between biomass and relative read abundance was found for two widely dispersed annelid taxa (Pygospio and Scoloplos). None of the other taxons, however, showed such a relationship. We discuss how quantification of abundance and biomass using molecular approaches are hampered by the ecology of DNA i.e. all the processes that determine the amount of DNA in the environment, including the ecology of the benthic species as well as the compositional nature of sequencing data.

18S rRNA amplicon sequence data (V1–V3) of the Bronx river estuary, New York

Characterising and monitoring biological diversity to foster sustainable ecosystems is highly recommended as urban centres rapidly expand. However, much of New York City’s biodiversity remains undescribed, including in the historically degraded, but recovering Bronx River Estuary. In a pilot study to identify organisms and characterise biodiversity patterns there, 18S rRNA gene amplicons (V1–V3 region), obtained from river sediments and surface waters of Hunts Point Riverside and Soundview Parks, were sequenced. Across 48 environmental samples collected over three seasons in 2015 and 2016, following quality control and contaminant removal, 2,763 Amplicon Sequence Variants (ASVs) were identified from 1,918,463 sequences. Rarefaction analysis showed sufficient sampling depth, and community composition varied over time and by substrate at the study sites over the sampling period. Protists, plants, fungi and animals, including organisms of management concern, such as Eastern oysters (Crassostrea virginica), wildlife pathogens and groups related to Harmful Algal Blooms, were detected. The most common taxa identified in river sediments were annelid worms, nematodes and diatoms. In the water column, the most commonly observed organisms were diatoms, algae of the phylum Cryptophyceae, ciliates and dinoflagellates. The presented dataset demonstrates the reach of 18S rRNA metabarcoding for characterising biodiversity in an urban estuary.

Towards reproducible metabarcoding data: Lessons from an international cross-laboratory experiment

Advances in high-throughput sequencing (HTS) are revolutionizing monitoring in marine environments by enabling rapid, accurate and holistic detection of species within complex biological samples. Research institutions worldwide increasingly employ HTS methods for biodiversity assessments. However, variance in laboratory procedures, analytical workflows and bioinformatic pipelines impede the transferability and comparability of results across research groups. An international experiment was conducted to assess the consistency of metabarcoding results derived from identical samples and primer sets using varying laboratory procedures. Homogenized biofouling samples collected from four coastal locations (Australia, Canada, New Zealand and the USA) were distributed to 12 independent laboratories. Participants were asked to follow one of two HTS library preparation workflows. While DNA extraction, primers and bioinformatic analyses were purposefully standardized to allow comparison, many other technical variables were allowed to vary among laboratories (amplification protocols, type of instrument used, etc.). Despite substantial variation observed in raw results, the primary signal in the data was consistent, with the samples grouping strongly by geographical origin for all data sets. Simple post hoc data clean-up by removing low-quality samples gave the best improvement in sample classification for nuclear 18S rRNA gene data, with an overall 92.81% correct group attribution. For mitochondrial COI gene data, the best classification result (95.58%) was achieved after correction for contamination errors. The identified critical methodological factors that introduced the greatest variability (preservation buffer, sample defrosting, template concentration, DNA polymerase, PCR enhancer) should be of great assistance in standardizing future biodiversity studies using metabarcoding.

Inventory of the benthic eukaryotic diversity in the oldest European lake

We have profound knowledge on biodiversity on Earth including plants and animals. In the recent decade, we have also increased our understanding on microorganisms in different hosts and the environment. However, biodiversity is not equally well studied among different biodiversity groups and Earth's systems with eukaryotes in freshwater sediments being among the least known. In this study, we used high-throughput sequencing of the 18S rRNA gene to investigate the entire diversity of benthic eukaryotes in three distinct habitats (littoral sediment and hard substrate, profundal sediment) of Lake Ohrid, the oldest European lake. Eukaryotic sequences were dominated by annelid and arthropod animals (54% of all eukaryotic reads) and protists (Ochrophyta and Ciliophora; together 40% of all reads). Eukaryotic diversity was 15% higher in the deep profundal than on either near-surface hard substrates or littoral sediments. The three habitats differed in their taxonomic and functional community composition. Specifically, heterotrophic organisms accounted for 92% of the reads in the profundal, whereas phototrophs accounted for 43% on the littoral hard substrate. The profundal community was the most homogeneous, and its network was the most complex, suggesting its highest stability among the sampled habitats.

Biases in bulk: DNA metabarcoding of marine communities and the methodology involved

With the growing anthropogenic pressure on marine ecosystems, the need for efficient monitoring of biodiversity grows stronger. DNA metabarcoding of bulk samples is increasingly being implemented in ecosystem assessments and is more cost-efficient and less time-consuming than monitoring based on morphology. However, before raw sequences are obtained from bulk samples, a profound number of methodological choices must be made. Here, we critically review the recent methods used for metabarcoding of marine bulk samples (including benthic, plankton and diet samples) and indicate how potential biases can be introduced throughout sampling, preprocessing, DNA extraction, marker and primer selection, PCR amplification and sequencing. From a total of 64 studies evaluated, our recommendations for best practices include to (a) consider DESS as a fixative instead of ethanol, (b) use the DNeasy PowerSoil kit for any samples containing traces of sediment, (c) not limit the marker selection to COI only, but preferably include multiple markers for higher taxonomic resolution, (d) avoid touchdown PCR profiles, (e) use a fixed annealing temperature for each primer pair when comparing across studies or institutes, (f) use a minimum of three PCR replicates, and (g) include both negative and positive controls. Although the implementation of DNA metabarcoding still faces several technical complexities, we foresee wide-ranging advances in the near future, including improved bioinformatics for taxonomic assignment, sequencing of longer fragments and the use of whole-genome information. Despite the bulk of biases involved in metabarcoding of bulk samples, if appropriate controls are included along the data generation process, it is clear that DNA metabarcoding provides a valuable tool in ecosystem assessments.

A Comparison of DNA Metabarcoding and Microscopy Methodologies for the Study of Aquatic Microbial Eukaryotes

The procedures and methodologies employed to study microbial eukaryotic plankton have been thoroughly discussed. Two main schools exist—one insisting on classic microscopy methodologies and the other supporting modern high-throughput sequencing (DNA metabarcoding). However, few studies have attempted to combine both these approaches; most studies implement one method while ignoring the other. This work aims to contribute to this discussion and examine the advantages and disadvantages of each methodology by comparing marine plankton community results from microscopy and DNA metabarcoding. The results obtained by the two methodologies do not vary significantly for Bacillariophyta, although they do for Dinoflagellata and Ciliophora. The lower the taxonomic level, the higher the inconsistency between the two methodologies for all the studied groups. Considering the different characteristics of microscopy-based identification and DNA metabarcoding, this work underlines that each method should be chosen depending on the aims of the study. DNA metabarcoding provides a better estimate of the taxonomic richness of an ecosystem while microscopy provides more accurate quantitative results regarding abundance and biomass. In any case, the combined use of the two methods, if properly standardized, can provide much more reliable and accurate results for the study of marine microbial eukaryotes.

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.

Expanding ecological assessment by integrating microorganisms into routine freshwater biomonitoring

Bioindication has become an indispensable part of water quality monitoring in most countries of the world, with the presence and abundance of bioindicator taxa, mostly multicellular eukaryotes, used for biotic indices. In contrast, microbes (bacteria, archaea and protists) are seldom used as bioindicators in routine assessments, although they have been recognized for their importance in environmental processes. Recently, the use of molecular methods has revealed unexpected diversity within known functional groups and novel metabolic pathways that are particularly important in energy and nutrient cycling. In various habitats, microbial communities respond to eutrophication, metals, and natural or anthropogenic organic pollutants through changes in diversity and function. In this review, we evaluated the common trends in these changes, documenting that they have value as bioindicators and can be used not only for monitoring but also for improving our understanding of the major processes in lotic and lentic environments. Current knowledge provides a solid foundation for exploiting microbial taxa, community structures and diversity, as well as functional genes, in novel monitoring programs. These microbial community measures can also be combined into biotic indices, improving the resolution of individual bioindicators. Here, we assess particular molecular approaches complemented by advanced bioinformatic analysis, as these are the most promising with respect to detailed bioindication value. We conclude that microbial community dynamics are a missing link important for our understanding of rapid changes in the structure and function of aquatic ecosystems, and should be addressed in the future environmental monitoring of freshwater ecosystems.

Identification of plastic-associated species in the Mediterranean Sea using DNA metabarcoding with Nanopore MinION

Plastic debris in the ocean form a new ecosystem, termed ‘plastisphere’, which hosts a variety of marine organisms. Recent studies implemented DNA metabarcoding to characterize the taxonomic composition of the plastisphere in different areas of the world. In this study, we used a modified metabarcoding approach which was based on longer barcode sequences for the characterization of the plastisphere biota. We compared the microbiome of polyethylene food bags after 1 month at sea to the free-living biome in two proximal but environmentally different locations on the Mediterranean coast of Israel. We targeted the full 1.5 kb-long 16S rRNA gene for bacteria and 0.4–0.8 kb-long regions within the 18S rRNA, ITS, tufA and COI loci for eukaryotes. The taxonomic barcodes were sequenced using Oxford Nanopore Technology with multiplexing on a single MinION flow cell. We identified between 1249 and 2141 species in each of the plastic samples, of which 61 species (34 bacteria and 27 eukaryotes) were categorized as plastic-specific, including species that belong to known hydrocarbon-degrading genera. In addition to a large prokaryotes repertoire, our results, supported by scanning electron microscopy, depict a surprisingly high biodiversity of eukaryotes within the plastisphere with a dominant presence of diatoms as well as other protists, algae and fungi.

Marine biomonitoring with eDNA: Can metabarcoding of water samples cut it as a tool for surveying benthic communities?

In the marine realm, biomonitoring using environmental DNA (eDNA) of benthic communities requires destructive direct sampling or the setting-up of settlement structures. Comparatively much less effort is required to sample the water column, which can be accessed remotely. In this study we assess the feasibility of obtaining information from the eukaryotic benthic communities by sampling the adjacent water layer. We studied two different rocky-substrate benthic communities with a technique based on quadrat sampling. We also took replicate water samples at four distances (0, 0.5, 1.5, and 20 m) from the benthic habitat. Using broad range primers to amplify a ca. 313 bp fragment of the cytochrome oxidase subunit I gene, we obtained a total of 3,543 molecular operational taxonomic units (MOTUs). The structure obtained in the two environments was markedly different, with Metazoa, Archaeplastida and Stramenopiles being the most diverse groups in benthic samples, and Hacrobia, Metazoa and Alveolata in the water. Only 265 MOTUs (7.5%) were shared between benthos and water samples and, of these, 180 (5.1%) were identified as benthic taxa that left their DNA in the water. Most of them were found immediately adjacent to the benthos, and their number decreased as we moved apart from the benthic habitat. It was concluded that water eDNA, even in the close vicinity of the benthos, was a poor proxy for the analysis of benthic structure, and that direct sampling methods are required for monitoring these complex communities via metabarcoding.

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.

MARES, a replicable pipeline and curated reference database for marine eukaryote metabarcoding

The use of DNA metabarcoding to characterise the biodiversity of environmental and community samples has exploded in recent years. However, taxonomic inferences from these studies are contingent on the quality and completeness of the sequence reference database used to characterise sample species-composition. In response, studies often develop custom reference databases to improve species assignment. The disadvantage of this approach is that it limits the potential for database re-use, and the transferability of inferences across studies. Here, we present the MARine Eukaryote Species (MARES) reference database for use in marine metabarcoding studies, created using a transparent and reproducible pipeline. MARES includes all COI sequences available in GenBank and BOLD for marine taxa, unified into a single taxonomy. Our pipeline facilitates the curation of sequences, synonymization of taxonomic identifiers used by different repositories, and formatting these data for use in taxonomic assignment tools. Overall, MARES provides a benchmark COI reference database for marine eukaryotes, and a standardised pipeline for (re)producing reference databases enabling integration and fair comparison of marine DNA metabarcoding results.

Measurement(s)	DNA
Technology Type(s)	bioinformatics analysis
Factor Type(s)	DNA sequence
Sample Characteristic - Organism	Eukaryota
Sample Characteristic - Environment	marine environment

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12324122

Comparing diversity levels in environmental samples: DNA sequence capture and metabarcoding approaches using 18S and COI genes

Environmental DNA studies targeting multiple taxa using metabarcoding provide remarkable insights into levels of species diversity in any habitat. The main drawbacks are the presence of primer bias and difficulty in identifying rare species. We tested a DNA sequence-capture method in parallel with the metabarcoding approach to reveal possible advantages of one method over the other. Both approaches were performed using the same eDNA samples and the same 18S and COI regions, followed by high throughput sequencing. Metabarcoded eDNA libraries were PCR amplified with one primer pair from 18S and COI genes. DNA sequence-capture libraries were enriched with 3,639 baits targeting the same gene regions. We tested amplicon sequence variants (ASVs) and operational taxonomic units (OTUs) in silico approaches for both markers and methods, using for this purpose the metabarcoding data set. ASVs methods uncovered more species for the COI gene, whereas the opposite occurred for the 18S gene, suggesting that clustering reads into OTUs could bias diversity richness especially using 18S with relaxed thresholds. Additionally, metabarcoding and DNA sequence-capture recovered 80%-90% of the control sample species. DNA sequence-capture was 8x more expensive, nonetheless it identified 1.5x more species for COI and 13x more genera for 18S than metabarcoding. Both approaches offer reliable results, sharing ca. 40% species and 72% families and retrieve more taxa when nuclear and mitochondrial markers are combined. eDNA metabarcoding is quite well established and low-cost, whereas DNA-sequence capture for biodiversity assessment is still in its infancy, is more time-consuming but provides more taxonomic assignments.

Nourished, Exposed Beaches Exhibit Altered Sediment Structure and Meiofaunal Communities

To retain recreational uses and shoreline protection, a large proportion of ocean beaches have been, and continue to be, nourished. Adding sand from subtidal and terrestrial sources to nourish beaches rarely re-creates the original sediment structure of the beach. Numerous studies have demonstrated that meiofaunal communities are altered by changes in sediment composition in low-energy substrates, therefore, we have explored whether beach nourishment has affected exposed, ocean beach meiofaunal communities. Since the early 2000s, we have conducted a series of sampling and experimental studies on meiofauna and sediments on nourished beaches in coastal North Carolina USA that had been sampled previously in the early 1970s, prior to any beach nourishment. Most of our studies consider meiofauna at the level of major taxa only. However, a few studies examine free-living flatworm (turbellarian) species in detail because of the existence of historical studies examining this group. Comparison of contemporary results to historical data and of heavily nourished versus lightly nourished beaches reveals extensive changes to beach sediment structure and meiofaunal community composition, indicating that the beaches are a more heterogeneous habitat than in the past. The effects of these substantial physical and biological changes to the production of beach ecosystem services are unlikely to be inconsequential.

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.

Environmental DNA survey captures patterns of fish and invertebrate diversity across a tropical seascape

Accurate, rapid, and comprehensive biodiversity assessments are critical for investigating ecological processes and supporting conservation efforts. Environmental DNA (eDNA) surveys show promise as a way to effectively characterize fine-scale patterns of community composition. We tested whether a single PCR survey of eDNA in seawater using a broad metazoan primer could identify differences in community composition between five adjacent habitats at 19 sites across a tropical Caribbean bay in Panama. We paired this effort with visual fish surveys to compare methods for a conspicuous taxonomic group. eDNA revealed a tremendous diversity of animals (8,586 operational taxonomic units), including many small taxa that would be undetected in traditional in situ surveys. Fish comprised only 0.07% of the taxa detected by a broad COI primer, yet included 43 species not observed in the visual survey. eDNA revealed significant differences in fish and invertebrate community composition across adjacent habitats and areas of the bay driven in part by taxa known to be habitat-specialists or tolerant to wave action. Our results demonstrate the ability of broad eDNA surveys to identify biodiversity patterns in the ocean.

DNA Metabarcoding of Deep-Sea Sediment Communities Using COI: Community Assessment, Spatio-Temporal Patterns and Comparison with 18S rDNA

Among the complex ecosystems and habitats that form the deep sea, submarine canyons and open slope systems are regarded as potential hotspots of biodiversity. We assessed the spatial and temporal patterns of biodiversity in sediment communities of a NW Mediterranean Canyon and its adjacent open slope (Blanes Canyon) with DNA metabarcoding. We sampled three layers of sediment and four different depths (900–1750 m) at two seasons, and used a fragment of the mitochondrial gene cytochrome c oxidase subunit I (COI) as a metabarcoding marker. The final dataset contained a total of 15,318 molecular operational taxonomic units (MOTUs). Metazoa, Stramenopiles and Archaeplastida were the dominant taxa and, within metazoans, Arthropoda, Nematoda and Cnidaria were the most diverse. There was a trend towards decreasing MOTU richness and diversity in the first few cm (1 to 5) of the sediment, with only 26.3% of the MOTUs shared across sediment layers. Our results show the presence of heterogeneous communities in the studied area, which was significantly different between zones, depths and seasons. We compared our results with the ones presented in a previous study, obtained using the v7 region of the 18S rRNA gene in the same samples. There were remarkable differences in the total number of MOTUs and in the most diverse taxa. COI recovered a higher number of MOTUs, but more remained unassigned taxonomically. However, the broad spatio-temporal patterns elucidated from both datasets coincided, with both markers retrieving the same ecological information. Our results showed that COI can be used to accurately characterize the studied communities and constitute a high-resolution method to detect ecological shifts. We also highlight that COI reference databases for deep-sea organisms have important gaps, and their completeness is essential in order to successfully apply metabarcoding techniques.

Simultaneous Metabarcoding of Eukaryotes and Prokaryotes to Elucidate the Community Structures within Tardigrade Microhabitats

Tardigrades are microscopic invertebrates that can withstand complete desiccation, but their interspecies interactions with prokaryotes and eukaryotes within their microhabitat remain relatively unexplored. Here, I utilized combined metabarcoding of eukaryotes and prokaryotes to simultaneously identify entire community structures within xeric and mesic mosses that harbor tardigrades. The populations of organisms within the microecosystems were successfully determined in 45 xeric moss samples and 47 mesic moss samples. Organismal composition was largely consistent regardless of the moss/lichen substrate, but significantly varied in the two tested locations, possibly because of the differences in environmental humidity. Xeric mosses containing xerophilic tardigrades and other anhydrobiotic invertebrates tended to have significantly limited biological diversity and prokaryotic population dominated by cyanobacteria, suggesting a selection due to extreme desiccation. A combined metabarcoding approach to identify both eukaryotes and prokaryotes can successfully elucidate community structures within microscopic ecosystems, and this can be a potential approach to study the microecology of meiofauna, including tardigrades.

Validating metabarcoding-based biodiversity assessments with multi-species occupancy models: A case study using coastal marine eDNA

Environmental DNA (eDNA) metabarcoding is an increasingly popular method for rapid biodiversity assessment. As with any ecological survey, false negatives can arise during sampling and, if unaccounted for, lead to biased results and potentially misdiagnosed environmental assessments. We developed a multi-scale, multi-species occupancy model for the analysis of community biodiversity data resulting from eDNA metabarcoding; this model accounts for imperfect detection and additional sources of environmental and experimental variation. We present methods for model assessment and model comparison and demonstrate how these tools improve the inferential power of eDNA metabarcoding data using a case study in a coastal, marine environment. Using occupancy models to account for factors often overlooked in the analysis of eDNA metabarcoding data will dramatically improve ecological inference, sampling design, and methodologies, empowering practitioners with an approach to wield the high-resolution biodiversity data of next-generation sequencing platforms.

From metabarcoding to metaphylogeography: separating the wheat from the chaff

Metabarcoding is by now a well-established method for biodiversity assessment in terrestrial, freshwater, and marine environments. Metabarcoding data sets are usually used for α- and β-diversity estimates, that is, interspecies (or inter-MOTU [molecular operational taxonomic unit]) patterns. However, the use of hypervariable metabarcoding markers may provide an enormous amount of intraspecies (intra-MOTU) information-mostly untapped so far. The use of cytochrome oxidase (COI) amplicons is gaining momentum in metabarcoding studies targeting eukaryote richness. COI has been for a long time the marker of choice in population genetics and phylogeographic studies. Therefore, COI metabarcoding data sets may be used to study intraspecies patterns and phylogeographic features for hundreds of species simultaneously, opening a new field that we suggest to name metaphylogeography. The main challenge for the implementation of this approach is the separation of erroneous sequences from true intra-MOTU variation. Here, we develop a cleaning protocol based on changes in entropy of the different codon positions of the COI sequence, together with co-occurrence patterns of sequences. Using a data set of community DNA from several benthic littoral communities in the Mediterranean and Atlantic seas, we first tested by simulation on a subset of sequences a two-step cleaning approach consisting of a denoising step followed by a minimal abundance filtering. The procedure was then applied to the whole data set. We obtained a total of 563 MOTUs that were usable for phylogeographic inference. We used semiquantitative rank data instead of read abundances to perform AMOVAs and haplotype networks. Genetic variability was mainly concentrated within samples, but with an important between seas component as well. There were intergroup differences in the amount of variability between and within communities in each sea. For two species, the results could be compared with traditional Sanger sequence data available for the same zones, giving similar patterns. Our study shows that metabarcoding data can be used to infer intra- and interpopulation genetic variability of many species at a time, providing a new method with great potential for basic biogeography, connectivity and dispersal studies, and for the more applied fields of conservation genetics, invasion genetics, and design of protected areas.

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.

Habitat Heterogeneity and Connectivity: Effects on the Planktonic Protist Community Structure at Two Adjacent Coastal Sites (the Lagoon and the Gulf of Venice, Northern Adriatic Sea, Italy) Revealed by Metabarcoding

The Lagoon of Venice (LoV) and the Gulf of Venice (GoV), two adjacent coastal Long Term Ecological Research (LTER) sites in the northern Adriatic Sea, represent a transitional/marine coupled ecosystem under the influence of regional and local factors. In this study, these sites were sampled on four dates from April 2016 to February 2017 for environmental DNA and relevant abiotic variables, aiming to assess the relative importance of habitat heterogeneity and connectivity in structuring the protist community. High Throughput Sequencing of V4-18S rRNA gene from 56 samples collected at seven stations produced ca 6 million reads, grouped into 7,336 Operational Taxonomic Units (OTUs) at 97% similarity, which were affiliated to protists belonging to 34 taxonomic groups. The whole community was dominated by Bacillariophyta, especially in spring-summer in the LoV, and by Dinophyta, mainly in the GoV. Ciliophora, Syndiniales, and Cryptophyceae were the next more abundant groups. The community structure varied across the seasons and was different in the two ecosystems, which shared 96% of the reads but showed a high proportion of OTUs distributed preferentially in one of the two sites (specialists) and a different partitioning of trophic categories. GoV specialists were mainly Dinophyceae (>56%), followed by Syndiniales and Bacillariophyta, while the LoV specialists were distributed among several groups, including Bacillariophyta, Syndiniales, Ciliophora, Cryptophyceae, and Trebouxiophyceae. The main abiotic drivers of the differences between protist communities were salinity and temperature, which however explained a minor part of the variance (17%), pointing at a higher relevance of biotic factors and inter-taxa relationships. This was more evident in the LoV, where the network analysis highlighted a higher number of OTUs’ connections than in the GoV. Overall, the metabarcoding approach allowed to depict the composition of the whole protist community in the lagoon and adjacent coastal waters with high resolution, revealing many taxa so far not reported in the area. In addition, despite no clear barrier to dispersal processes, differences in the relative abundance and temporal variability of local protist communities indicate that environmental heterogeneity, in these adjacent and connected ecosystems, can be strong enough to allow for ecological segregation.

GenBank is a reliable resource for 21st century biodiversity research

Traditional methods of characterizing biodiversity are increasingly being supplemented and replaced by approaches based on DNA sequencing alone. These approaches commonly involve extraction and high-throughput sequencing of bulk samples from biologically complex communities or samples of environmental DNA (eDNA). In such cases, vouchers for individual organisms are rarely obtained, often unidentifiable, or unavailable. Thus, identifying these sequences typically relies on comparisons with sequences from genetic databases, particularly GenBank. While concerns have been raised about biases and inaccuracies in laboratory and analytical methods, comparatively little attention has been paid to the taxonomic reliability of GenBank itself. Here we analyze the metazoan mitochondrial sequences of GenBank using a combination of distance-based clustering and phylogenetic analysis. Because of their comparatively rapid evolutionary rates and consequent high taxonomic resolution, mitochondrial sequences represent an invaluable resource for the detection of the many small and often undescribed organisms that represent the bulk of animal diversity. We show that metazoan identifications in GenBank are surprisingly accurate, even at low taxonomic levels (likely <1% error rate at the genus level). This stands in contrast to previously voiced concerns based on limited analyses of particular groups and the fact that individual researchers currently submit annotated sequences to GenBank without significant external taxonomic validation. Our encouraging results suggest that the rapid uptake of DNA-based approaches is supported by a bioinformatic infrastructure capable of assessing both the losses to biodiversity caused by global change and the effectiveness of conservation efforts aimed at slowing or reversing these losses.