Improving metabarcoding taxonomic assignment: A case study of fishes in a large marine ecosystem

Environmental DNA exploring the distribution of Indo-Pacific humpback dolphins and fish diversity from the Pearl River Estuary, China

Understanding the spatial distribution patterns of endangered species is crucial for their protection. However, gathering such information for cetaceans remains challenging due to their underwater life and elusive nature. Here, by conuding a three-year field ecological survey, we first employed environmental DNA (eDNA) technology to elucidate the distribution and its influencing factors of the Indo-Pacific humpback dolphins from their largest habitat, the Pearl River Estuary (PRE), China. The self-designed SCDloop primers exhibited high specificity for the dolphins and demonstrated elevated eDNA concentrations in proximity to areas with significant human activity in the PRE. Fish species detected in the PRE included most prey items consumed by the dolphins. A significant positive correlation between the occurrence of humpback dolphins and fish diversity was found, indicating the crucial role of fishery resources in the dolphins' habitat selection. Our findings support the use of eDNA technology as a supplementary tool for monitoring marine cetaceans.

Bridging the gaps through environmental DNA: A review of critical considerations for interpreting the biodiversity data in coral reef ecosystems

Coral reefs, the rainforests of the sea, are vital hotspots for marine biodiversity. However, the persistent challenge of climate change directly threatens the delicate balance of coral reef ecosystems, impacting myriad species and critical ecosystem services. Therefore, this comprehensive review critically discusses the associated challenges in assessing and preserving coral reef diversity, emphasizing the need for novel biomonitoring techniques due to the elusive and cryptic nature of many reef organisms. The review focuses on environmental DNA (eDNA) analysis as a non-invasive tool for coral species monitoring at various ecological levels. The review highlights that using eDNA in coral reef monitoring requires careful consideration of multiple factors, such as strategic assay development, optimization, and marker selection, substrate selection, and sample volume, which are critical for maximizing the probability of species detection. Moreover, integrating environmental RNA (eRNA) provides additional insights into temporal aspects advancing the coral reef biodiversity research and conservation efforts.

Enhancing metabarcoding efficiency and ecological insights through integrated taxonomy and DNA reference barcoding: A case study on beach meiofauna

Molecular techniques like metabarcoding, while promising for exploring diversity of communities, are often impeded by the lack of reference DNA sequences available for taxonomic annotation. Our study explores the benefits of combining targeted DNA barcoding and morphological taxonomy to improve metabarcoding efficiency, using beach meiofauna as a case study. Beaches are globally important ecosystems and are inhabited by meiofauna, microscopic animals living in the interstitial space between the sand grains, which play a key role in coastal biodiversity and ecosystem dynamics. However, research on meiofauna faces challenges due to limited taxonomic expertise and sparse sampling. We generated 775 new cytochrome c oxidase I DNA barcodes from meiofauna specimens collected along the Netherlands' west coast and combined them with the NCBI GenBank database. We analysed alpha and beta diversity in 561 metabarcoding samples from 24 North Sea beaches, a region extensively studied for meiofauna, using both the enriched reference database and the NCBI database without the additional reference barcodes. Our results show a 2.5-fold increase in sequence annotation and a doubling of species-level Operational Taxonomic Units (OTUs) identification when annotating the metabarcoding data with the enhanced database. Additionally, our analyses revealed a bell-shaped curve of OTU richness across the intertidal zone, aligning more closely with morphological analysis patterns, and more defined community dissimilarity patterns between supralittoral and intertidal sites. Our research highlights the importance of expanding molecular reference databases and combining morphological taxonomy with molecular techniques for biodiversity assessments, ultimately improving our understanding of coastal ecosystems.

Endemic Radiation of African Moonfish, Selene dorsalis (Gill 1863), in the Eastern Atlantic: Mitogenomic Characterization and Phylogenetic Implications of Carangids (Teleostei: Carangiformes)

This study offers an in-depth analysis of the mitochondrial genome of Selene dorsalis (Gill 1863), a species native to the Eastern Atlantic Ocean. The circular mitochondrial DNA molecule measures 16,541 base pairs and comprises 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA genes, and a control region (CR). The nucleotide composition exhibits a notable adenine-thymine (AT) bias, accounting for 53.13%, which aligns with other species in the Carangidae family. Most PCGs initiate with the ATG codon, with the exception of Cytochrome C oxidase subunit I, which starts with GTG. Analysis of relative synonymous codon usage reveals that leucine and serine are the most prevalent amino acids in the mitochondrial genome of S. dorsalis and its congeners (S. vomer and S. setapinnis). All tRNAs display the typical cloverleaf structure, though tRNA Serine (S1) lacks a dihydrouracil arm. Pairwise comparisons of synonymous and nonsynonymous substitutions for all PCGs yielded values below '1', indicating strong purifying selection. The CR spans 847 bp, representing 5.12% of the mitochondrial genome, and is characterized by high AT content (62.81%). It is situated between tRNA-Pro (TGG) and tRNA-Phe (GAA). The CR contains conserved sequence blocks, with CSB-1 being the longest at 22 bp and CSB-D the shortest at 18 bp. Phylogenetic analysis, using Bayesian and Maximum-likelihood trees constructed from concatenated PCGs across 72 species, successfully differentiates S. dorsalis from other carangids. This study also explores how ocean currents and gyres might influence lineage diversification and parapatric speciation of Selene species between the Atlantic and Pacific Oceans. These results highlight the importance of the mitochondrial genome in elucidating the structural organization and evolutionary dynamics of S. dorsalis and its relatives within marine ecosystems.

Examination of sequence variations in partial mitochondrial 12S gene amongst damselfish species as references for DNA barcoding

Accurate species identification, based on DNA barcoding, can be achieved when sufficient sequence variations are present amongst species in the sampled marker. In general, the ability to discriminate species decreases with shorter sequences; however, shorter regions have a merit in amplification success by the polymerase chain reaction. In either case, it is important to investigate sequence variations amongst species before barcoding to understand its reliability and limitations. In this study, we investigate how accurately short, but hypervariable portion of the mitochondrial 12S ribosomal RNA (12S) gene (MiFish region with approximately 180 bp) is used to identify each species in diversified pomacentrid fishes compared with the longer region of the same gene (approximately 750 bp). We prepared three datasets with 301 sequences of the MiFish region for 150 species, the same 301 of sequences of the longer 12S region and 476 sequences of the MiFish region for 183 species. Neighbour-joining (NJ) analyses and genetic distance analyses revealed several indistinguishable pairs of species in these DNA regions. Although the number of such pairs was larger in the MiFish region, 83.6% (153 of 183) of species possessed respective unique sequences even in the MiFish region (versus 96.0% [144 of 150 species] in the longer 12S region). A part of indistinguishable pairs of species might have caused by mitochondrial DNA introgressions and taxonomically unresolved problems. Our analysis clarified the effectiveness and limitations of species identification using DNA barcoding for Pomacentridae and the sequences we provided here contribute to the expansion of references for pomacentrid mitochondrial 12S sequences.

Using low volume eDNA methods to sample pelagic marine animal assemblages

Environmental DNA (eDNA) is an increasingly useful method for detecting pelagic animals in the ocean but typically requires large water volumes to sample diverse assemblages. Ship-based pelagic sampling programs that could implement eDNA methods generally have restrictive water budgets. Studies that quantify how eDNA methods perform on low water volumes in the ocean are limited, especially in deep-sea habitats with low animal biomass and poorly described species assemblages. Using 12S rRNA and COI gene primers, we quantified assemblages comprised of micronekton, coastal forage fishes, and zooplankton from low volume eDNA seawater samples (n = 436, 380-1800 mL) collected at depths of 0-2200 m in the southern California Current. We compared diversity in eDNA samples to concurrently collected pelagic trawl samples (n = 27), detecting a higher diversity of vertebrate and invertebrate groups in the eDNA samples. Differences in assemblage composition could be explained by variability in size-selectivity among methods and DNA primer suitability across taxonomic groups. The number of reads and amplicon sequences variants (ASVs) did not vary substantially among shallow (<200 m) and deep samples (>600 m), but the proportion of invertebrate ASVs that could be assigned a species-level identification decreased with sampling depth. Using hierarchical clustering, we resolved horizontal and vertical variability in marine animal assemblages from samples characterized by a relatively low diversity of ecologically important species. Low volume eDNA samples will quantify greater taxonomic diversity as reference libraries, especially for deep-dwelling invertebrate species, continue to expand.

Assigning the unassigned: A signature-based classification of rDNA metabarcodes reveals new deep-sea diversity

Environmental DNA metabarcoding reveals a vast genetic diversity of marine eukaryotes. Yet, most of the metabarcoding data remain unassigned due to the paucity of reference databases. This is particularly true for the deep-sea meiofauna and eukaryotic microbiota, whose hidden diversity is largely unexplored. Here, we tackle this issue by using unique DNA signatures to classify unknown metabarcodes assigned to deep-sea foraminifera. We analyzed metabarcoding data obtained from 311 deep-sea sediment samples collected in the Clarion-Clipperton Fracture Zone, an area of potential polymetallic nodule exploitation in the Eastern Pacific Ocean. Using the signatures designed in the 37F hypervariable region of the 18S rRNA gene, we were able to classify 802 unassigned metabarcodes into 61 novel lineages, which have been placed in 27 phylogenetic clades. The comparison of new lineages with other foraminiferal datasets shows that most novel lineages are widely distributed in the deep sea. Five lineages are also present in the shallow-water datasets; however, phylogenetic analysis of these lineages separates deep-sea and shallow-water metabarcodes except in one case. While the signature-based classification does not solve the problem of gaps in reference databases, this taxonomy-free approach provides insight into the distribution and ecology of deep-sea species represented by unassigned metabarcodes, which could be useful in future applications of metabarcoding for environmental monitoring.

Environmental DNA reveals temporal variation in mesophotic reefs of the Humboldt upwelling ecosystems of central Chile: Toward a baseline for biodiversity monitoring of unexplored marine habitats

Temperate mesophotic reef ecosystems (TMREs) are among the least known marine habitats. Information on their diversity and ecology is geographically and temporally scarce, especially in highly productive large upwelling ecosystems. Lack of information remains an obstacle to understanding the importance of TMREs as habitats, biodiversity reservoirs and their connections with better-studied shallow reefs. Here, we use environmental DNA (eDNA) from water samples to characterize the community composition of TMREs on the central Chilean coast, generating the first baseline for monitoring the biodiversity of these habitats. We analyzed samples from two depths (30 and 60 m) over four seasons (spring, summer, autumn, and winter) and at two locations approximately 16 km apart. We used a panel of three metabarcodes, two that target all eukaryotes (18S rRNA and mitochondrial COI) and one specifically targeting fishes (16S rRNA). All panels combined encompassed eDNA assigned to 42 phyla, 90 classes, 237 orders, and 402 families. The highest family richness was found for the phyla Arthropoda, Bacillariophyta, and Chordata. Overall, family richness was similar between depths but decreased during summer, a pattern consistent at both locations. Our results indicate that the structure (composition) of the mesophotic communities varied predominantly with seasons. We analyzed further the better-resolved fish assemblage and compared eDNA with other visual methods at the same locations and depths. We recovered eDNA from 19 genera of fish, six of these have also been observed on towed underwater videos, while 13 were unique to eDNA. We discuss the potential drivers of seasonal differences in community composition and richness. Our results suggest that eDNA can provide valuable insights for monitoring TMRE communities but highlight the necessity of completing reference DNA databases available for this region.

rCRUX: A Rapid and Versatile Tool for Generating Metabarcoding Reference libraries in R

The sequencing revolution requires accurate taxonomic classification of DNA sequences. Key to making accurate taxonomic assignments are curated, comprehensive reference barcode databases. However, the generation and curation of such databases has remained challenging given the large and continuously growing volumes of both DNA sequence data and novel reference barcode targets. Monitoring and research applications require a greater diversity of specialized gene regions and targeted taxa then are currently curated by professional staff. Thus there is a growing need for an easy to implement computational tool that can generate comprehensive metabarcoding reference libraries for any bespoke locus. We address this need by reimagining CRUX from the Anacapa Toolkit and present the rCRUX package in R which, like it's predecessor, relies on sequence homology and PCR primer compatibility instead of keyword-searches to avoid limitations of user-defined metadata. The typical workflow involves searching for plausible seed amplicons (get_seeds_local() or get_seeds_remote()) by simulating in silico PCR to acquire a set of sequences analogous to PCR products containing a user-defined set of primer sequences. Next, these seeds are used to iteratively blast search seed sequences against a local copy of the National Center for Biotechnology Information (NCBI) formatted nt database using a taxonomic-rank based stratified random sampling approach ( blast_seeds() ). This results in a comprehensive set of sequence matches. This database is dereplicated and cleaned (derep_and_clean_db()) by identifying identical reference sequences and collapsing the taxonomic path to the lowest taxonomic agreement across all matching reads. This results in a curated, comprehensive database of primer-specific reference barcode sequences from NCBI. Databases can then be compared (compare_db()) to determine read and taxonomic overlap. We demonstrate that rCRUX provides more comprehensive reference databases for the MiFish Universal Teleost 12S, Taberlet trnl, fungal ITS, and Leray CO1 loci than CRABS, MetaCurator, RESCRIPt, and ecoPCR reference databases. We then further demonstrate the utility of rCRUX by generating 24 reference databases for 20 metabarcoding loci, many of which lack dedicated reference database curation efforts. The rCRUX package provides a simple to use tool for the generation of curated, comprehensive reference databases for user-defined loci, facilitating accurate and effective taxonomic classification of metabarcoding and DNA sequence efforts broadly.

DNA metabarcoding reveals spatial and temporal variation of fish eye fluke communities in lake ecosystems

Eye flukes (Diplostomidae) are diverse and abundant trematode parasites that form multi-species communities in fish with negative effects on host fitness and survival. However, the environmental factors and host-related characteristics that determine species diversity, composition, and coexistence in such communities remain poorly understood. Here, we developed a cost-effective cox1 region-specific DNA metabarcoding approach to characterize parasitic diplostomid communities in two common fish species (Eurasian perch and common roach) collected from seven temperate lakes in Estonia. We found considerable inter- and intra-lake, as well as inter-host species, variation in diplostomid communities. Sympatric host species characterization revealed that parasite communities were typically more diverse in roach than perch. Additionally, we detected five positive and two negative diplostomid species associations in roach, whereas only a single negative association was observed in perch. These results indicate that diplostomid communities in temperate lakes are complex and dynamic systems exhibiting both spatial and temporal heterogeneity. They are influenced by various environmental factors and by host-parasite and inter-parasite interactions. We expect that the described methodology facilitates ecological and biodiversity research of diplostomid parasites. It is also adaptable to other parasite groups where it could serve to improve current understanding of diversity, distribution, and interspecies interactions of other understudied taxa.

A comparison of biomonitoring methodologies for surf zone fish communities

Surf zones are highly dynamic marine ecosystems that are subject to increasing anthropogenic and climatic pressures, posing multiple challenges for biomonitoring. Traditional methods such as seines and hook and line surveys are often labor intensive, taxonomically biased, and can be physically hazardous. Emerging techniques, such as baited remote underwater video (BRUV) and environmental DNA (eDNA) are promising nondestructive tools for assessing marine biodiversity in surf zones of sandy beaches. Here we compare the relative performance of beach seines, BRUV, and eDNA in characterizing community composition of bony (teleost) and cartilaginous (elasmobranch) fishes of surf zones at 18 open coast sandy beaches in southern California. Seine and BRUV surveys captured overlapping, but distinct fish communities with 50% (18/36) of detected species shared. BRUV surveys more frequently detected larger species (e.g. sharks and rays) while seines more frequently detected one of the most abundant species, barred surfperch (Amphistichus argenteus). In contrast, eDNA metabarcoding captured 88.9% (32/36) of all fishes observed in seine and BRUV surveys plus 57 additional species, including 15 that frequent surf zone habitats. On average, eDNA detected over 5 times more species than BRUVs and 8 times more species than seine surveys at a given site. eDNA approaches also showed significantly higher sensitivity than seine and BRUV methods and more consistently detected 31 of the 32 (96.9%) jointly observed species across beaches. The four species detected by BRUV/seines, but not eDNA were only resolved at higher taxonomic ranks (e.g. Embiotocidae surfperches and Sygnathidae pipefishes). In frequent co-detection of species between methods limited comparisons of richness and abundance estimates, highlighting the challenge of comparing biomonitoring approaches. Despite potential for improvement, results overall demonstrate that eDNA can provide a cost-effective tool for long-term surf zone monitoring that complements data from seine and BRUV surveys, allowing more comprehensive surveys of vertebrate diversity in surf zone habitats.

rCRUX: A Rapid and Versatile Tool for Generating Metabarcoding Reference libraries in R

Key to making accurate taxonomic assignments are curated, comprehensive reference barcode databases. However, the generation and curation of such databases has remained challenging given the large and continuously growing volumes of DNA sequence data and novel reference barcode targets. Monitoring and research applications require a greater diversity of specialized gene regions and targeted taxa to meet taxonomic classification goals then are currently curated by professional staff. Thus, there is a growing need for an easy to implement tool that can generate comprehensive metabarcoding reference libraries for any bespoke locus. We address this need by reimagining CRUX from the Anacapa Toolkit and present the rCRUX package in R. The typical workflow involves searching for plausible seed amplicons (get_seeds_local() or get_seeds_remote()) by simulating in silico PCR to acquire seed sequences containing a user-defined primer set. Next these seeds are used to iteratively blast search seed sequences against a local NCBI formatted database using a taxonomic rank based stratified random sampling approach (blast_seeds()) that results in a comprehensive set of sequence matches. This database is dereplicated and cleaned (derep_and_clean_db()) by identifying identical reference sequences and collapsing the taxonomic path to the lowest taxonomic agreement across all matching reads. This results in a curated, comprehensive database of primer specific reference barcode sequences from NCBI. We demonstrate that rCRUX provides more comprehensive reference databases for the MiFish Universal Teleost 12S, Taberlet trnl, and fungal ITS locus than CRABS, METACURATOR, RESCRIPt, and ECOPCR reference databases. We then further demonstrate the utility of rCRUX by generating 16 reference databases for metabarcoding loci that lack dedicated reference database curation efforts. The rCRUX package provides a simple to use tool for the generation of curated, comprehensive reference databases for user-defined loci, facilitating accurate and effective taxonomic classification of metabarcoding and DNA sequence efforts broadly.

Maximizing the reliability and the number of species assignments in metabarcoding studies using a curated regional library and a public repository

Biodiversity assessments relying on DNA have increased rapidly over the last decade. However, the reliability of taxonomic assignments in metabarcoding studies is variable and affected by the reference databases and the assignment methods used. Species level assignments are usually considered as reliable using regional libraries but unreliable using public repositories. In this study, we aimed to test this assumption for metazoan species detected in the Gulf of St. Lawrence in the Northwest Atlantic. We first created a regional library (GSL-rl) by data mining COI barcode sequences from BOLD, and included a reliability ranking system for species assignments. We then estimated 1) the accuracy and precision of the public repository NCBI-nt for species assignments using sequences from the regional library and 2) compared the detection and reliability of species assignments of a metabarcoding dataset using either NCBI-nt or the regional library and popular assignment methods. With NCBI-nt and sequences from the regional library, the BLAST-LCA (least common ancestor) method was the most precise method for species assignments, but the accuracy was higher with the BLAST-TopHit method (&gt;80% over all taxa, between 70% and 90% amongst taxonomic groups). With the metabarcoding dataset, the reliability of species assignments was greater using GSL-rl compared to NCBI-nt. However, we also observed that the total number of reliable species assignments could be maximized using both GSL-rl and NCBI-nt with different optimized assignment methods. The use of a two-step approach for species assignments, i.e., using a regional library and a public repository, could improve the reliability and the number of detected species in metabarcoding studies.

Using eDNA sampling for species-specific fish detection in tropical oceanic samples: limitations and recommendations for future use

Background

Over the past decade, environmental DNA (eDNA) has become a resourceful tool in conservation and biomonitoring. Environmental DNA has been applied in a variety of environments, but the application to studies of marine fish, particularly at tropical latitudes, are limited. Since many commercially important Caribbean fishes are overexploited, these species are optimal candidates to explore the use of this method as a biomonitoring tool. Specifically, for many of these species, the formation of fish spawning aggregations (FSAs) marks a critical life history event where fishes will gather in large numbers for reproduction. These FSAs are ephemeral in nature, lasting only a few days, but are predictable in time and space which makes them susceptible to overfishing.

Methods

In this study, we test the feasibility of using an eDNA sampling approach (water and sediment collection) to detect the presence of known FSAs off the west coast of Puerto Rico, with cytochrome c oxidase subunit 1 (CO1) and 12S rRNA (12S) primers designed to target specific species. A total of 290 eDNA samples were collected and, of those, 206 eDNA samples were processed. All eDNA samples varied in DNA concentration, both between replicates and collection methods. A total of 12 primer sets were developed and tested using traditional PCR and qPCR.

Results

Despite validation of primer accuracy and sample collection during known peak spawning times, the use of traditional PCR and qPCR with both molecular markers failed to produce species-specific amplification. Thus, a trial test was conducted using the CO1 primers in which target fish DNA was 'spiked' at various concentrations into the respective eDNA samples to determine the target species DNA concentration limit of detection. Upon successful amplification of the trial, results indicated that eDNA samples were below the detection threshold of our methods, suggesting that the number of fish present at the spawning aggregations was inadequate for single-species detection methods. In addition, elements such as the unavoidable presence of non-target DNA, oceanic environmental conditions, shedding rates of target fish, among other biotic and abiotic factors could have affected DNA persistence and degradation rates at the sites.

Conclusion

We provide recommendations for species-specific fish detection in lower latitudes, and suggestions for studies aiming to monitor or detect fish spawning aggregations using eDNA sampling.

Signal and noise in metabarcoding data

Metabarcoding is a powerful molecular tool for simultaneously surveying hundreds to thousands of species from a single sample, underpinning microbiome and environmental DNA (eDNA) methods. Deriving quantitative estimates of underlying biological communities from metabarcoding is critical for enhancing the utility of such approaches for health and conservation. Recent work has demonstrated that correcting for amplification biases in genetic metabarcoding data can yield quantitative estimates of template DNA concentrations. However, a major source of uncertainty in metabarcoding data stems from non-detections across technical PCR replicates where one replicate fails to detect a species observed in other replicates. Such non-detections are a special case of variability among technical replicates in metabarcoding data. While many sampling and amplification processes underlie observed variation in metabarcoding data, understanding the causes of non-detections is an important step in distinguishing signal from noise in metabarcoding studies. Here, we use both simulated and empirical data to 1) suggest how non-detections may arise in metabarcoding data, 2) outline steps to recognize uninformative data in practice, and 3) identify the conditions under which amplicon sequence data can reliably detect underlying biological signals. We show with both simulations and empirical data that, for a given species, the rate of non-detections among technical replicates is a function of both the template DNA concentration and species-specific amplification efficiency. Consequently, we conclude metabarcoding datasets are strongly affected by (1) deterministic amplification biases during PCR and (2) stochastic sampling of amplicons during sequencing-both of which we can model-but also by (3) stochastic sampling of rare molecules prior to PCR, which remains a frontier for quantitative metabarcoding. Our results highlight the importance of estimating species-specific amplification efficiencies and critically evaluating patterns of non-detection in metabarcoding datasets to better distinguish environmental signal from the noise inherent in molecular detections of rare targets.

A manager's guide to using eDNA metabarcoding in marine ecosystems

Environmental DNA (eDNA) metabarcoding is a powerful tool that can enhance marine ecosystem/biodiversity monitoring programs. Here we outline five important steps managers and researchers should consider when developing eDNA monitoring program: (1) select genes and primers to target taxa; (2) assemble or develop comprehensive barcode reference databases; (3) apply rigorous site occupancy based decontamination pipelines; (4) conduct pilot studies to define spatial and temporal variance of eDNA; and (5) archive samples, extracts, and raw sequence data. We demonstrate the importance of each of these considerations using a case study of eDNA metabarcoding in the Ports of Los Angeles and Long Beach. eDNA metabarcoding approaches detected 94.1% (16/17) of species observed in paired trawl surveys while identifying an additional 55 native fishes, providing more comprehensive biodiversity inventories. Rigorous benchmarking of eDNA metabarcoding results improved ecological interpretation and confidence in species detections while providing archived genetic resources for future analyses. Well designed and validated eDNA metabarcoding approaches are ideally suited for biomonitoring applications that rely on the detection of species, including mapping invasive species fronts and endangered species habitats as well as tracking range shifts in response to climate change. Incorporating these considerations will enhance the utility and efficacy of eDNA metabarcoding for routine biomonitoring applications.

Fishing for fish environmental DNA: Ecological applications, methodological considerations, surveying designs, and ways forward

Vast global declines of freshwater and marine fish diversity and population abundance pose serious threats to both ecosystem sustainability and human livelihoods. Environmental DNA (eDNA)-based biomonitoring provides robust, efficient, and cost-effective assessment of species occurrences and population trends in diverse aquatic environments. Thus, it holds great potential for improving conventional surveillance frameworks to facilitate fish conservation and fisheries management. However, the many technical considerations and rapid developments underway in the eDNA arena can overwhelm researchers and practitioners new to the field. Here, we systematically analysed 416 fish eDNA studies to summarize research trends in terms of investigated targets, research aims, and study systems, and reviewed the applications, rationales, methodological considerations, and limitations of eDNA methods with an emphasis on fish and fisheries research. We highlighted how eDNA technology may advance our knowledge of fish behaviour, species distributions, population genetics, community structures, and ecological interactions. We also synthesized the current knowledge of several important methodological concerns, including the qualitative and quantitative power eDNA has to recover fish biodiversity and abundance, and the spatial and temporal representations of eDNA with respect to its sources. To facilitate ecological applications implementing fish eDNA techniques, recent literature was summarized to generate guidelines for effective sampling in lentic, lotic, and marine habitats. Finally, we identified current gaps and limitations, and pointed out newly emerging research avenues for fish eDNA. As methodological optimization and standardization improve, eDNA technology should revolutionize fish monitoring and promote biodiversity conservation and fisheries management that transcends geographic and temporal boundaries.

Detection of human pathogenic bacteria in rectal DNA samples from Zalophus californianus in the Gulf of California, Mexico

Human intrusions into undisturbed wildlife areas greatly contribute to the emergence of infectious diseases. To minimize the impacts of novel emerging infectious diseases (EIDs) on human health, a comprehensive understanding of the microbial species that reside within wildlife species is required. The Gulf of California (GoC) is an example of an undisturbed ecosystem. However, in recent decades, anthropogenic activities within the GoC have increased. Zalophus californianus has been proposed as the main sentinel species in the GoC; hence, an assessment of sea lion bacterial microbiota may reveal hidden risks for human health. We evaluated the presence of potential human pathogenic bacterial species from the gastrointestinal (GI) tracts of wild sea lions through a metabarcoding approach. To comprehensively evaluate this bacterial consortium, we considered the genetic information of six hypervariable regions of 16S rRNA. Potential human pathogenic bacteria were identified down to the species level by integrating the RDP and Pplacer classifier outputs. The combined genetic information from all analyzed regions suggests the presence of at least 44 human pathogenic bacterial species, including Shigella dysenteriae and Bacillus anthracis. Therefore, the risks of EIDs from this area should be not underestimated.

Multi-locus evaluation of gastrointestinal bacterial communities from Zalophus californianus pups in the Gulf of California, México

Background

The gastrointestinal (GI) bacterial communities of sea lions described to date have occasionally revealed large intraspecific variability, which may originate from several factors including different methodological approaches. Indeed, GI bacterial community surveys commonly rely on the use of a single hypervariable region (HR) of 16S rRNA, which may result in misleading structural interpretations and limit comparisons among studies. Here, we considered a multi-locus analysis by targeting six HRs of 16S rRNA with the aims of (i) comprehensively assessing the GI bacterial consortium in rectal samples from Zalophus californianus pups and (ii) elucidating structural variations among the tested HRs. In addition, we evaluated which HRs may be most suitable for identifying intrinsic, structurally related microbiome characteristics, such as geographic variations or functional capabilities.

Methods

We employed a Short MUltiple Regions Framework (SMURF) approach using the Ion 16S™ Metagenomic Kit. This kit provides different proprietary primers designed to target six HRs of the 16S rRNA gene. To date, the only analytical pipeline available for this kit is the Ion Reporter™ Software of Thermo Fisher Scientific. Therefore, we propose an in-house pipeline to use with open-access tools, such as QIIME2 and PICRUSt 2, in downstream bioinformatic analyses.

Results

As hypothesized, distinctive bacterial community profiles were observed for each analyzed HR. A higher number of bacterial taxa were detected with the V3 and V6-V7 regions. Conversely, the V8 and V9 regions were less informative, as we detected a lower number of taxa. The synergistic information of these HRs suggests that the GI microbiota of Zalophus californianus pups is predominated by five bacterial phyla: Proteobacteria (~50%), Bacteroidetes (~20%), Firmicutes (~18%), Fusobacteria (~7%), and Epsilonbacteraeota (~4%). Notably, our results differ at times from previously reported abundance profiles, which may promote re-evaluations of the GI bacterial compositions in sea lions and other pinniped species that have been reported to date. Moreover, consistent geographic differences were observed only with the V3, V4, and V6-V7 regions. In addition, these HRs also presented higher numbers of predicted molecular pathways, although no significant functional changes were apparent. Together, our results suggests that multi-locus analysis should be encouraged in GI microbial surveys, as single-locus approaches may result in misleading structural results that hamper the identification of structurally related microbiome features.

Holistic Impact Evaluation of Human Activities on the Coastal Fish Biodiversity in the Chinese Coastal Environment

Ecological qualities and resources in coasts are threatened by various human activities, such as pollution and fishery. Impact evaluation of environmental stressors over a wide coastal stretch has been limited due to lack of efficient and standardizable biodiversity monitoring and assessment tools. Integrating environmental DNA (eDNA) and ecological traits, a holistic approach was developed to assess the impact of pollution and aquaculture on fish biodiversity in Chinese coastal areas. Taking the Yalujiang Estuary (YLJK) from the Yellow Sea and the Nan'ao Island Area (NAO) from the South China Sea as cases, the performance of the eDNA biomonitoring workflow was validated. First, the eDNA results of 22 sampling sites reached more than 85% of the asymptotes of species or ASVs in each area. A total of 115 fish species in both areas were detected and NAO was 1.8 times richer than YLJK using eDNA and the fish eDNA composition was consistent with the historical data. eDNA recovered distinct variations of fish sequence, taxonomic and functional diversity, and the corresponding trends following the offshore distance between the two areas. Fish sequence diversity was decreased primarily by estuarine pollution factors (chemical oxygen demand and zinc) in the YLJK. Compared with no breeding areas, lower fish sequence diversity was in breeding areas in the NAO. By integrating ecological traits, the eDNA approach offers promising opportunities for future fish biodiversity monitoring and assessment in national and global coastal environments.

Fish Diversity Monitored by Environmental DNA in the Yangtze River Mainstream

Surveys and assessments based on environmental DNA are not only efficient and time-saving, but also cause less harm to monitoring targets. Environmental DNA has become a common tool for the assessment and monitoring of aquatic organisms. In this study, we investigated fish resources in the Yangtze River mainstream using environmental DNA, and the variations in fish during two seasons (spring and autumn) were compared. The results showed that 13 species were identified in spring, and nine species of fish were identified in autumn. The fish with higher eDNA detection were Sinibotia superciliaris, Tachysurus fulvidraco, Cyprinus carpio, Ctenopharyngodon Idella, Monopterus albus, Acanthogobius hasta, Saurogobio dabryi, Oncorhynchus mykiss, Mugil cephalus, Odontamblyopus rubicundus. Seasonal variation between spring and autumn was not significant, and the environmental factors had different effects on fish assemblages during the two seasons. Our study used the eDNA technique to monitor the composition of fish in the spring and autumn in the Yangtze River mainstream, providing a new technology for the long-term management and protection of fishery resources in the region. Of course, problems such as pollution and insufficient databases are the current shortcomings of environmental DNA, which will be the focus of our future research and study.

Fish eDNA metabarcoding from aquatic biofilm samples: Methodological aspects

Fish eDNA metabarcoding is usually performed from filtered water samples. The volume of filtered water depends on the study scope and can rapidly become time consuming according to the number of samples that have to be processed. To avoid time allocated to filtration, passive DNA samplers have been used to recover fish eDNA from marine environments faster. In freshwater ecosystems, aquatic biofilms were used to catch eDNA from macroinvertebrates. Here, we test the capacity of aquatic biofilms to entrap fish eDNA in a large lake and, therefore, the possibility to perform fish eDNA metabarcoding from this matrix compared to the traditional fish eDNA approach from filtered water samples. Methodological aspects of the use of aquatic biofilms for fish eDNA metabarcoding (e.g. PCR replicates, biological replicates, bioinformatics pipeline, reference database and taxonomic assignment) were validated against a mock community. When using biofilms from habitats sheltered from wind and waves, biofilm and water approach provided similar inventories. Richness and diversity were comparable between both approaches. Approaches differed only for rare taxa. Our results illustrate the capacity of aquatic biofilms to act as passive eDNA samplers of fish eDNA and, therefore, the possibility to use biofilms to monitor fish communities efficiently from biofilms. Furthermore, our results open up avenues of research to study a diversity of biological groups (among which bioindicators as diatoms, macroinvertebrates and fish) from eDNA isolated from a single environmental matrix reducing sampling efforts, analysis time and costs.

eDNA captures depth partitioning in a kelp forest ecosystem

Environmental DNA (eDNA) metabarcoding is an increasingly important tool for surveying biodiversity in marine ecosystems. However, the scale of temporal and spatial variability in eDNA signatures, and how this variation may impact eDNA-based marine biodiversity assessments, remains uncertain. To address this question, we systematically examined variation in vertebrate eDNA signatures across depth (0 m to 10 m) and horizontal space (nearshore kelp forest and surf zone) over three successive days in Southern California. Across a broad range of teleost fish and elasmobranchs, results showed significant variation in species richness and community assemblages between surface and depth, reflecting microhabitat depth preferences of common Southern California nearshore rocky reef taxa. Community assemblages between nearshore and surf zone sampling stations at the same depth also differed significantly, consistent with known habitat preferences. Additionally, assemblages also varied across three sampling days, but 69% of habitat preferences remained consistent. Results highlight the sensitivity of eDNA in capturing fine-scale vertical, horizontal, and temporal variation in marine vertebrate communities, demonstrating the ability of eDNA to capture a highly localized snapshot of marine biodiversity in dynamic coastal environments.

Characterizing Industrial and Artisanal Fishing Vessel Catch Composition Using Environmental DNA and Satellite-Based Tracking Data

The decline in wild-caught fisheries paired with increasing global seafood demand is pushing the need for seafood sustainability to the forefront of national and regional priorities. Validation of species identity is a crucial early step, yet conventional monitoring and surveillance tools are limited in their effectiveness because they are extremely time-consuming and require expertise in fish identification. DNA barcoding methods are a versatile tool for the genetic monitoring of wildlife products; however, they are also limited by requiring individual tissue samples from target specimens which may not always be possible given the speed and scale of seafood operations. To circumvent the need to individually sample organisms, we pilot an approach that uses forensic environmental DNA (eDNA) metabarcoding to profile fish species composition from the meltwater in fish holds on industrial and artisanal fishing vessels in Ecuador. Fish identified genetically as present were compared to target species reported by each vessel’s crew. Additionally, we contrasted the geographic range of identified species against the satellite-based fishing route data of industrial vessels to determine if identified species could be reasonably expected in the catch.