Environmental DNA illuminates the dark diversity of sharks

Use of environmental DNA in assessment of fish functional and phylogenetic diversity

Assessing the impact of global changes and protection effectiveness is a key step in monitoring marine fishes. Most traditional census methods are demanding or destructive. Nondisturbing and nonlethal approaches based on video and environmental DNA are alternatives to underwater visual census or fishing. However, their ability to detect multiple biodiversity factors beyond traditional taxonomic diversity is still unknown. For bony fishes and elasmobranchs, we compared the performance of eDNA metabarcoding and long-term remote video to assess species' phylogenetic and functional diversity. We used 10 eDNA samples from 30 L of water each and 25 hr of underwater videos over 4 days on Malpelo Island (pacific coast of Colombia), a remote marine protected area. Metabarcoding of eDNA detected 66% more molecular operational taxonomic units (MOTUs) than species on video. We found 66 and 43 functional entities with a single eDNA marker and videos, respectively, and higher functional richness for eDNA than videos. Despite gaps in genetic reference databases, eDNA also detected a higher fish phylogenetic diversity than videos; accumulation curves showed how 1 eDNA transect detected as much phylogenetic diversity as 25 hr of video. Environmental DNA metabarcoding can be used to affordably, efficiently, and accurately census biodiversity factors in marine systems. Although taxonomic assignments are still limited by species coverage in genetic reference databases, use of MOTUs highlights the potential of eDNA metabarcoding once reference databases have expanded.

Development of an eDNA metabarcoding tool for surveying the world’s largest amphibian

Due to the overexploitation of farming, as well as habitat destruction, the wild population of Chinese giant salamander (CGS) Andrias davidianus, a species with seven genetically distinct lineages, has decreased by over 80% in the past 70 years. Traditional survey methods have proven to be unsuitable for finding this rare and elusive species. We evaluated the efficacy of environmental DNA (eDNA) sampling to detect CGS indirectly from its aquatic environment. We developed several species-specific primer sets; validated their specificity and sensitivity; and assessed their utility in silico, in the laboratory, and at two field sites harboring released farm-bred CGS. We detected the presence of CGS DNA by using polymerase chain reaction and Sanger sequencing. We also sequenced an amplicon mixture of seven haplotype-represented samples using high-throughput sequencing. Our eDNA methods could detect the presence of CGS at moderate densities reported across its range, proving them as a cost-effective way to establish broad-scale patterns of occupancy for CGS. In addition, our primers enabled the detection of mitochondrial lineage mixture or introduced individuals from geographically isolated populations of CGS.

How many replicates to accurately estimate fish biodiversity using environmental DNA on coral reefs?

Quantifying fish species diversity in rich tropical marine environments remains challenging. Environmental DNA (eDNA) metabarcoding is a promising tool to face this challenge through the filtering, amplification, and sequencing of DNA traces from water samples. However, because eDNA concentration is low in marine environments, the reliability of eDNA to detect species diversity can be limited. Using an eDNA metabarcoding approach to identify fish Molecular Taxonomic Units (MOTUs) with a single 12S marker, we aimed to assess how the number of sampling replicates and filtered water volume affect biodiversity estimates. We used a paired sampling design of 30 L per replicate on 68 reef transects from 8 sites in 3 tropical regions. We quantified local and regional sampling variability by comparing MOTU richness, compositional turnover, and compositional nestedness. We found strong turnover of MOTUs between replicated pairs of samples undertaken in the same location, time, and conditions. Paired samples contained non-overlapping assemblages rather than subsets of one another. As a result, non-saturated localized diversity accumulation curves suggest that even 6 replicates (180 L) in the same location can underestimate local diversity (for an area <1 km). However, sampling regional diversity using ~25 replicates in variable locations (often covering 10 s of km) often saturated biodiversity accumulation curves. Our results demonstrate variability of diversity estimates possibly arising from heterogeneous distribution of eDNA in seawater, highly skewed frequencies of eDNA traces per MOTU, in addition to variability in eDNA processing. This high compositional variability has consequences for using eDNA to monitor temporal and spatial biodiversity changes in local assemblages. Avoiding false-negative detections in future biomonitoring efforts requires increasing replicates or sampled water volume to better inform management of marine biodiversity using eDNA.

Observations of sharks (Elasmobranchii) at Europa Island, a remote marine protected area important for shark conservation in the southern Mozambique Channel

Sharks have declined worldwide and remote sanctuaries are becoming crucial for shark conservation. The southwest Indian Ocean is a hotspot of both terrestrial and marine biodiversity mostly impacted by anthropogenic damage. Sharks were observed during surveys performed from April to June 2013 in the virtually pristine coral reefs around Europa Island, a remote Marine Protected Area located in the southern Mozambique Channel. Observation events comprised 67 1-hour scientific dives between 5 – 35m depth and 7 snorkeling inspections, as well as 4 dinghy-based observations in the shallow lagoon. In a period of 24 days, 475 sharks were tallied. Carcharhinus galapagensis was most encountered and contributed 20% of the abundance during diving, followed by C. albimarginatus (10%). Both species were more abundant between 11-14h, and on the exposed sides of the island. Numbers of Sphyrna lewini were highest with 370 individuals windward and leeward, mostly schooling. S. lewini aggregations in the area are hypothesized to be attracted to the seamount archipelago offering favorable conditions for deep incursions and of which Europa Island forms part. C. amblyrhynchos, Galeocerdo cuvier and S. mokarran were uncommon, while there was an additional observation of Rhincodon typus. The lagoon of Europa was a nursery ground for C. melanopterus where it was the only species present. A total of 8 species was recorded, contributing to the shark diversity of 15 species reported from Europa since 1952 in the scientific and gray literature. Overall, with the occurrence of several species of apex predators in addition to that of R. typus, large schools of S. lewini, fair numbers of reef sharks and a nursery of C. melanopterus, Europa’s sharks constitute a significant reservoir of biodiversity, which contributes to preserve the functioning of the ecosystem. Our observations highlight the relevance of Europa Island for shark conservation and the need for shark-targeted management in the EEZ of both Europa and Bassas da India.

Dark diversity of flea assemblages of small mammalian hosts: effects of environment, host traits and host phylogeny

An assemblage of species in a locality comprises two components, namely (i) species that are present (realised diversity) and (ii) species from the regional pool that may potentially inhabit this locality due to suitable ecological conditions, but that are absent (dark diversity). We investigated factors affecting the dark diversity of component communities of fleas parasitic on small mammals in the northern Palearctic at two scales. First, we considered the dark diversity of flea assemblages of the same host (for 13 host species) across regions and tested for the effects of environmental factors and the number of available host species on the dark diversity of within-region flea assemblages. Second, we considered the dark diversity of fleas across host species within a region (for 20 regions) and asked whether within-host dark diversity is associated with host phylogeny and/or traits. We found that the dark diversity of flea assemblages harboured by small mammals varied substantially (i) within the same host species across space (in 12 of 13 host species) and (ii) between host species within a region (in eight of 20 regions). The size of the dark diversity of flea assemblages of the same host across regions was generally affected by environmental factors (mainly by the amount of green vegetation), whereas the size of the dark diversity of flea assemblages of a host species within a region was affected by host traits (mainly by the degree of host sociality and the structure of its shelter and, to a lesser degree, by its geographic range size) but was not associated with host phylogenetic affinities. We conclude that application of the dark diversity concept to parasite communities across space or hosts allows a better understanding of the factors affecting the species richness and composition of these communities.

Seeing through sedimented waters: environmental DNA reduces the phantom diversity of sharks and rays in turbid marine habitats

BACKGROUND

Sharks and rays are some of the most threatened marine taxa due to the high levels of bycatch and significant demand for meat and fin-related products in many Asian communities. At least 25% of shark and ray species are considered to be threatened with extinction. In particular, the density of reef sharks in the Pacific has declined to 3-10% of pre-human levels. Elasmobranchs are thought to be sparse in highly urbanised and turbid environments. Low visibility coupled with the highly elusive behaviour of sharks and rays pose a challenge to diversity estimation and biomonitoring efforts as sightings are limited to chance encounters or from carcasses ensnared in nets. Here we utilised an eDNA metabarcoding approach to enhance the precision of elasmobranch diversity estimates in urbanised marine environments.

RESULTS

We applied eDNA metabarcoding on seawater samples to detect elasmobranch species in the hyper-urbanised waters off Singapore. Two genes-vertebrate 12S and elasmobranch COI-were targeted and amplicons subjected to Illumina high-throughput sequencing. With a total of 84 water samples collected from nine localities, we found 47 shark and ray molecular operational taxonomic units, of which 16 had species-level identities. When data were compared against historical collections and contemporary sightings, eDNA detected 14 locally known species as well as two potential new records.

CONCLUSIONS

Local elasmobranch richness uncovered by eDNA is greater than the seven species sighted over the last two decades, thereby reducing phantom diversity. Our findings demonstrate that eDNA metabarcoding is effective in detecting shark and ray species despite the challenges posed by the physical environment, granting a more consistent approach to monitor these highly elusive and threatened species.

Modelling Marine Predator Habitat Using the Abundance of Its Pelagic Prey in the Tropical South-Western Pacific

Understanding the ecological mechanisms underpinning distribution patterns is vital in managing populations of mobile marine species. This study is a first step towards an integrated description of the habitats and spatial distributions of marine predators in the Natural Park of the Coral Sea, one of the world’s largest marine-protected areas at about 1.3 million km2, covering the entirety of New Caledonia’s pelagic waters. The study aims to quantify the benefit of including a proxy for prey abundance in predator niche modelling, relative to other marine physical variables. Spatial distributions and relationships with environmental data were analysed using catch per unit of effort data for three fish species (albacore tuna, yellowfin tuna and dolphinfish), sightings collected from aerial surveys for three cetacean guilds (Delphininae, Globicephalinae and Ziphiidae) and foraging locations identified from bio-tracking for three seabird species (wedge-tailed shearwater, Tahiti petrel and red-footed booby). Predator distributions were modelled as a function of a static covariate (bathymetry), oceanographic covariates (sea surface temperature, chlorophyll-a concentration and 20 °C-isotherm depth) and an acoustically derived micronekton preyscape covariate. While distributions were mostly linked to bathymetry for seabirds, and chlorophyll and temperature for fish and cetaceans, acoustically derived prey abundance proxies slightly improved distribution models for all fishes and seabirds except the Tahiti petrel, but not for the cetaceans. Predicted spatial distributions showed that pelagic habitats occupied by predator fishes did not spatially overlap. Finally, predicted habitats and the use of the preyscapes in predator habitat modelling were discussed.

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

DNA metabarcoding is an important tool for molecular ecology. However, its effectiveness hinges on the quality of reference sequence databases and classification parameters employed. Here we evaluate the performance of MiFish 12S taxonomic assignments using a case study of California Current Large Marine Ecosystem fishes to determine best practices for metabarcoding. Specifically, we use a taxonomy cross-validation by identity framework to compare classification performance between a global database comprised of all available sequences and a curated database that only includes sequences of fishes from the California Current Large Marine Ecosystem. We demonstrate that the regional database provides higher assignment accuracy than the comprehensive global database. We also document a tradeoff between accuracy and misclassification across a range of taxonomic cutoff scores, highlighting the importance of parameter selection for taxonomic classification. Furthermore, we compared assignment accuracy with and without the inclusion of additionally generated reference sequences. To this end, we sequenced tissue from 597 species using the MiFish 12S primers, adding 252 species to GenBank's existing 550 California Current Large Marine Ecosystem fish sequences. We then compared species and reads identified from seawater environmental DNA samples using global databases with and without our generated references, and the regional database. The addition of new references allowed for the identification of 16 additional native taxa representing 17.0% of total reads from eDNA samples, including species with vast ecological and economic value. Together these results demonstrate the importance of comprehensive and curated reference databases for effective metabarcoding and the need for locus-specific validation efforts.

Species on the move around the Australian coastline: A continental-scale review of climate-driven species redistribution in marine systems

Climate-driven changes in the distribution of species are a pervasive and accelerating impact of climate change, and despite increasing research effort in this rapidly emerging field, much remains unknown or poorly understood. We lack a holistic understanding of patterns and processes at local, regional and global scales, with detailed explorations of range shifts in the southern hemisphere particularly under-represented. Australian waters encompass the world's third largest marine jurisdiction, extending from tropical to sub-Antarctic climate zones, and have waters warming at rates twice the global average in the north and two to four times in the south. Here, we report the results of a multi-taxon continent-wide review describing observed and predicted species redistribution around the Australian coastline, and highlight critical gaps in knowledge impeding our understanding of, and response to, these considerable changes. Since range shifts were first reported in the region in 2003, 198 species from nine Phyla have been documented shifting their distribution, 87.3% of which are shifting poleward. However, there is little standardization of methods or metrics reported in observed or predicted shifts, and both are hindered by a lack of baseline data. Our results demonstrate the importance of historical data sets and underwater visual surveys, and also highlight that approximately one-fifth of studies incorporated citizen science. These findings emphasize the important role the public has had, and can continue to play, in understanding the impact of climate change. Most documented shifts are of coastal fish species in sub-tropical and temperate systems, while tropical systems in general were poorly explored. Moreover, most distributional changes are only described at the poleward boundary, with few studies considering changes at the warmer, equatorward range limit. Through identifying knowledge gaps and research limitations, this review highlights future opportunities for strategic research effort to improve the representation of Australian marine species and systems in climate-impact research.

Maximizing regional biodiversity requires a mosaic of protection levels

Protected areas are the flagship management tools to secure biodiversity from anthropogenic impacts. However, the extent to which adjacent areas with distinct protection levels host different species numbers and compositions remains uncertain. Here, using reef fishes, European alpine plants, and North American birds, we show that the composition of species in adjacent Strictly Protected, Restricted, and Non-Protected areas is highly dissimilar, whereas the number of species is similar, after controlling for environmental conditions, sample size, and rarity. We find that between 12% and 15% of species are only recorded in Non-Protected areas, suggesting that a non-negligible part of regional biodiversity occurs where human activities are less regulated. For imperiled species, the proportion only recorded in Strictly Protected areas reaches 58% for fishes, 11% for birds, and 7% for plants, highlighting the fundamental and unique role of protected areas and their environmental conditions in biodiversity conservation.

This study shows that the dissimilarity in species composition between sites with different levels of protection is consistently high, suggesting that adjacent and connected areas with different protection levels host very dissimilar species assemblages.

Environmental DNA monitoring of oncogenic viral shedding and genomic profiling of sea turtle fibropapillomatosis reveals unusual viral dynamics

Pathogen-induced cancers account for 15% of human tumors and are a growing concern for endangered wildlife. Fibropapillomatosis is an expanding virally and environmentally co-induced sea turtle tumor epizootic. Chelonid herpesvirus 5 (ChHV5) is implicated as a causative virus, but its transmission method and specific role in oncogenesis and progression is unclear. We applied environmental (e)DNA-based viral monitoring to assess viral shedding as a direct means of transmission, and the relationship between tumor burden, surgical resection and ChHV5 shedding. To elucidate the abundance and transcriptional status of ChHV5 across early, established, regrowth and internal tumors we conducted genomics and transcriptomics. We determined that ChHV5 is shed into the water column, representing a likely transmission route, and revealed novel temporal shedding dynamics and tumor burden correlations. ChHV5 was more abundant in the water column than in marine leeches. We also revealed that ChHV5 is latent in fibropapillomatosis, including early stage, regrowth and internal tumors; higher viral transcription is not indicative of poor patient outcome, and high ChHV5 loads predominantly arise from latent virus. These results expand our knowledge of the cellular and shedding dynamics of ChHV5 and can provide insights into temporal transmission dynamics and viral oncogenesis not readily investigable in tumors of terrestrial species.

Environmental DNA metabarcoding reveals and unpacks a biodiversity conservation paradox in Mediterranean marine reserves

Although we are currently experiencing worldwide biodiversity loss, local species richness does not always decline under anthropogenic pressure. This conservation paradox may also apply in protected areas but has not yet received conclusive evidence in marine ecosystems. Here, we survey fish assemblages in six Mediterranean no-take reserves and their adjacent fishing grounds using environmental DNA (eDNA) while controlling for environmental conditions. We detect less fish species in marine reserves than in nearby fished areas. The paradoxical gradient in species richness is accompanied by a marked change in fish species composition under different managements. This dissimilarity is mainly driven by species that are often overlooked by classical visual surveys but detected with eDNA: cryptobenthic, pelagic, and rare fishes. These results do not negate the importance of reserves in protecting biodiversity but shed new light on how under-represented species groups can positively react to fishing pressure and how conservation efforts can shape regional biodiversity patterns.

Bioinformatic pipelines combining denoising and clustering tools allow for more comprehensive prokaryotic and eukaryotic metabarcoding

Environmental DNA metabarcoding is a powerful tool for studying biodiversity. However, bioinformatic approaches need to adjust to the diversity of taxonomic compartments targeted as well as to each barcode gene specificities. We built and tested a pipeline based on read correction with DADA2 allowing analysing metabarcoding data from prokaryotic (16S) and eukaryotic (18S, COI) life compartments. We implemented the option to cluster amplicon sequence variants (ASVs) into operational taxonomic units (OTUs) with swarm, a network-based clustering algorithm, and the option to curate ASVs/OTUs using LULU. Finally, taxonomic assignment was implemented via the Ribosomal Database Project Bayesian classifier (RDP) and BLAST. We validated this pipeline with ribosomal and mitochondrial markers using metazoan mock communities and 42 deep-sea sediment samples. The results show that ASVs and OTUs describe different levels of biotic diversity, the choice of which depends on the research questions. They underline the advantages and complementarity of clustering and LULU-curation for producing metazoan biodiversity inventories at a level approaching the one obtained using morphological criteria. While clustering removes intraspecific variation, LULU effectively removes spurious clusters, originating from errors or intragenomic variability. Swarm clustering affected alpha and beta diversity differently depending on genetic marker. Specifically, d-values > 1 appeared to be less appropriate with 18S for metazoans. Similarly, increasing LULU's minimum ratio level proved essential to avoid losing species in sample-poor data sets. Comparing BLAST and RDP underlined that accurate assignments of deep-sea species can be obtained with RDP, but highlighted the need for a concerted effort to build comprehensive, ecosystem-specific databases.

The Promise and Pitfalls of Environmental DNA and RNA Approaches for the Monitoring of Human and Animal Pathogens from Aquatic Sources

Novel forensics-inspired molecular approaches have revolutionized species detection in the wild and are particularly useful for tracing endangered or invasive species. These new environmental DNA or RNA (eDNA or eRNA)–based techniques are now being applied to human and animal pathogen surveillance, particularly in aquatic environments. They allow better disease monitoring (presence or absence and geographical spread) and understanding of pathogen occurrence and transmission, benefitting species conservation and, more recently, our understanding of the COVID-19 global human pandemic. In the present article, we summarize the benefits of eDNA-based monitoring, highlighted by two case studies: The first is a fibropapillomatosis tumor-associated herpesvirus (chelonid herpesvirus 5) driving a sea turtle panzootic, and the second relates to eRNA-based detection of the SARS-CoV-2 coronavirus driving the COVID-19 human pandemic. The limitations of eDNA- or eRNA-based approaches are also summarized, and future directions and recommendations of the field are discussed. Continuous eDNA- or eRNA-based monitoring programs can potentially improve human and animal health by predicting disease outbreaks in advance, facilitating proactive rather than reactive responses.

Detection of the elusive Dwarf sperm whale (Kogia sima) using environmental DNA at Malpelo island (Eastern Pacific, Colombia)

Monitoring large marine mammals is challenging due to their low abundances in general, an ability to move over large distances and wide geographical range sizes.The distribution of the pygmy (Kogia breviceps) and dwarf (Kogia sima) sperm whales is informed by relatively rare sightings, which does not permit accurate estimates of their distribution ranges. Hence, their conservation status has long remained Data Deficient (DD) in the Red list of the International Union for Conservation of Nature (IUCN), which prevent appropriate conservation measures.Environmental DNA (eDNA) metabarcoding uses DNA traces left by organisms in their environments to detect the presence of targeted taxon, and is here proved to be useful to increase our knowledge on the distribution of rare but emblematic megafauna.Retrieving eDNA from filtered surface water provides the first detection of the Dwarf sperm whale (Kogia sima) around the remote Malpelo island (Colombia).Environmental DNA collected during oceanic missions can generate better knowledge on rare but emblematic animals even in regions that are generally well sampled for other taxa.

Integrating eDNA metabarcoding and simultaneous underwater visual surveys to describe complex fish communities in a marine biodiversity hotspot

Marine biodiversity can be surveyed using underwater visual censuses and recently with eDNA metabarcoding. Although a promising tool, eDNA studies have shown contrasting results related to its detection scale and the number of species identified compared to other survey methods. Also, its accuracy relies on complete reference databases used for taxonomic assignment and, as other survey methods, species detection may show false-negative and false-positive errors. Here, we compared results from underwater visual censuses and simultaneous eDNA metabarcoding fish surveys in terms of observed species and community composition. We also assess the effect of a custom reference database in the taxonomic assignment, and evaluate occupancy, capture and detection probabilities, as well as error rates of eDNA survey data. We amplified a 12S rRNA fish barcode from 24 sampling sites in the gulf of California. More species were detected with eDNA metabarcoding than with UVC. Because each survey method largely detected different sets of species, the combined approach doubled the number of species registered. Both survey methods recovered a known biodiversity gradient and a biogeographic break, but eDNA captured diversity over a broader geographic and bathymetric scale. Furthermore, the use of a modest-sized custom reference database significantly increased taxonomic assignment. In a subset of species, occupancy models revealed eDNA surveys provided similar or higher detection probabilities compared to UVC. The occupancy value of each species had a large influence on eDNA detectability, and in the false positive and negative error. Overall, these results highlight the potential of eDNA metabarcoding in complementing other established ecological methods for studies of marine fishes.

Development and validation of a quantitative qPCR assay for detecting Natterjack toad (Epidalea calamita) eDNA samples

The Natterjack toad (Epidalea calamita) is the rarest amphibian species in Ireland, regionally Red-Listed as Endangered. We applied an eDNA approach to detect species presence in breeding pond water samples. We developed a species-specific qPCR assay targeting the cytochrome c oxidase subunit I (COI). The assay was tested in silico, in vitro (DNA extracted from tissue) and in vivo (DNA extracted from water samples). Water samples were collected from five ponds with known Natterjack toad presence or absence to validate the sensitivity and specificity of the assay. The assay was shown to be highly specific to the Natterjack toad and tested positive only against toad tissue samples and eDNA samples from ponds with known species presence. We believe this method can be used for rapid assessment of species occurrence.

Methodology of fish eDNA and its applications in ecology and environment

Fish environmental DNA (eDNA) studies have made substantial progress during the past decade, and significant advances in monitoring fishes have been gained by taking advantage of this technology. Although a number of reviews concerning eDNA are available and some recent fish eDNA reviews focused on fisheries or standard method have been published, a systematic review of methodology of fish eDNA and its applications in ecology and environment has not yet been published. To our knowledge, this is the first review of fish eDNA for solving ecological and environmental issues. First, the most comprehensive literature analysis of fish eDNA was presented and analyzed. Then, we systematically discuss the relevant experiments and analyses of fish eDNA, and infers that standard workflow is on the way to consensus. We additionally provide reference sequence databases and the primers used to amplify the reference sequences or detecting fish eDNA. The abiotic and biotic conditions affecting fish eDNA persistence are also summarized in a schematic diagram. Subsequently, we focus on the major achievements of fish eDNA in ecology and environment. We additionally highlight the exciting new tools, including in situ autonomous monitoring devices, CRISPR nucleic acid detection technology, and meta-omics technology for fish eDNA detection in future. Ultimately, methodology of fish eDNA will provide a wholly new paradigm for conservation actions of fishes, ecological and environmental management at a global scale.

Genome-scale target capture of mitochondrial and nuclear environmental DNA from water samples

Environmental DNA (eDNA) provides a promising supplement to traditional sampling methods for population genetic inferences, but current studies have almost entirely focused on short mitochondrial markers. Here, we develop one mitochondrial and one nuclear set of target capture probes for the whale shark (Rhincodon typus) and test them on seawater samples collected in Qatar to investigate the potential of target capture for eDNA‐based population studies. The mitochondrial target capture successfully retrieved ~235× (90× − 352× per base position) coverage of the whale shark mitogenome. Using a minor allele frequency of 5%, we find 29 variable sites throughout the mitogenome, indicative of at least five contributing individuals. We also retrieved numerous mitochondrial reads from an abundant nontarget species, mackerel tuna (Euthynnus affinis), showing a clear relationship between sequence similarity to the capture probes and the number of captured reads. The nuclear target capture probes retrieved only a few reads and polymorphic variants from the whale shark, but we successfully obtained millions of reads and thousands of polymorphic variants with different allele frequencies from E. affinis. We demonstrate that target capture of complete mitochondrial genomes and thousands of nuclear loci is possible from aquatic eDNA samples. Our results highlight that careful probe design, taking into account the range of divergence between target and nontarget sequences as well as presence of nontarget species at the sampling site, is crucial to consider. eDNA sampling coupled with target capture approaches provide an efficient means with which to retrieve population genomic data from aggregating and spawning aquatic species.

Using environmental DNA and occupancy modelling to estimate rangewide metapopulation dynamics

We demonstrate the power of combining two emergent tools for resolving rangewide metapopulation dynamics. First, we employed environmental DNA (eDNA) surveys to efficiently generate multiseason rangewide site occupancy histories. Second, we developed a novel dynamic, spatial multiscale occupancy model to estimate metapopulation dynamics. The model incorporates spatial relationships, explicitly accounts for non-detection bias and allows direct evaluation of the drivers of extinction and colonization. We applied these tools to examine metapopulation dynamics of endangered tidewater goby, a species endemic to California estuarine habitats. We analysed rangewide eDNA data from 190 geographically isolated sites (813 total water samples) surveyed from 2 years (2016 and 2017). Rangewide estimates of the proportion of sites that were occupied varied little between 2016 (0.52) and 2017 (0.51). However, there was evidence of extinction and colonization dynamics. The probability of extinction of an occupied site (0.106) and probability of colonization of an unoccupied site (0.085) were nearly equal. Stability in site occupancy proportions combined with nearly equal rates of extinction and colonization suggests a dynamic equilibrium between the 2 years surveyed. Assessment of covariate effects revealed that colonization probability increased as the number of occupied neighbouring sites increased and as distance between occupied sites decreased. We show that eDNA surveys can rapidly provide a snapshot of a species distribution over a broad geographic range and, when these surveys are paired with occupancy modelling, can uncover metapopulation dynamics and their drivers.

Environmental DNA effectively captures functional diversity of coastal fish communities

Robust assessments of taxonomic and functional diversity are essential components of research programmes aimed at understanding current biodiversity patterns and forecasting trajectories of ecological changes. Yet, evaluating marine biodiversity along its dimensions is challenging and dependent on the power and accuracy of the available data collection methods. Here we combine three traditional survey methodologies (underwater visual census strip transects [UVCt], baited underwater videos [BUV] and small-scale fishery catches [SSFc]), and one novel molecular technique (environmental DNA metabarcoding [eDNA]-12S rRNA and cytochrome oxidase subunit 1 [COI]) to investigate their efficiency and complementarity in assessing fish diversity. We analysed 1,716 multimethod replicates at a basin scale to measure the taxonomic and functional diversity of Mediterranean fish assemblages. Taxonomic identities were investigated at species, genus and family levels. Functional identities were assessed using combinations of morphological, behavioural and trophic traits. We show that: (a) SSFc provided the higher taxonomic diversity estimates followed by eDNA, and then UVCt and BUV; (b) eDNA was the only method able to gather the whole spectrum of considered functional traits, showing the most functionally diversified and least redundant fish assemblages; and (c) the effectiveness of eDNA in describing functional structure reflected its lack of selectivity towards any considered functional trait. Our findings suggest that the reach of eDNA analysis stretches beyond taxon detection efficiency and provides new insights into the potential of metabarcoding in ecological studies.

Lost and found: Frogs in a biodiversity hotspot rediscovered with environmental DNA

Declines and extinctions are increasing globally and challenge conservationists to keep pace with biodiversity monitoring. Organisms leave DNA traces in the environment, e.g., in soil, water, and air. These DNA traces are referred to as environmental DNA (eDNA). The analysis of eDNA is a highly sensitive method with the potential to rapidly assess local diversity and the status of threatened species. We searched for DNA traces of 30 target amphibian species of conservation concern, at different levels of threat, using an environmental DNA metabarcoding approach, together with an extensive sequence reference database to analyse water samples from six montane sites in the Atlantic Coastal Forest and adjacent Cerrado grasslands of Brazil. We successfully detected DNA traces of four declined species (Hylodes ornatus, Hylodes regius, Crossodactylus timbuhy, and Vitreorana eurygnatha); two locally disappeared (Phasmahyla exilis and Phasmahyla guttata); and one species that has not been seen since 1968 (putatively assigned to Megaelosia bocainensis). We confirm the presence of species undetected by traditional methods, underscoring the efficacy of eDNA metabarcoding for biodiversity monitoring at low population densities, especially in megadiverse tropical sites. Our results support the potential application of eDNA in conservation biology, to evaluate persistence and distribution of threatened species in surveyed habitats or sites, and improve accuracy of red lists, especially for species undetected over long periods.

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.

Environmental DNA detection tracks established seasonal occurrence of blacktip sharks (Carcharhinus limbatus) in a semi-enclosed subtropical bay

The integration of eDNA analysis into the population assessment and monitoring of sharks could greatly improve temporal and spatial data used for management purposes. This study aimed to compare eDNA detection against well-established seasonal changes in blacktip shark (Carcharhinus limbatus) abundance in Terra Ceia Bay (FL, USA). We used a species-specific real-time PCR approach to detect C. limbatus eDNA in the bay on a near monthly basis from spring through mid-fall in 2018 and 2019. Previous studies have shown that C. limbatus give birth in the bay in early summer and immature sharks occur there until late fall, when decreasing water temperatures cause them to move offshore and southwards. Water samples (2 L) were collected (4–6 per month) and filtered in the field, with each then being subjected to real-time PCR. Carcharhinus limbatus ‘positive’ filters were significantly more commonly collected during the April-July sampling period than during the August-October sampling period. While following the predicted pattern, eDNA concentration was generally too low for accurate quantification. Our results show that C. limbatus eDNA detection follows known seasonal residency patterns consistently over 2 years of monitoring. Species-specific eDNA analysis using real-time PCR could therefore represent a cost-effective, scalable sampling tool to facilitate improved shark population monitoring in semi-enclosed marine habitats.

Accumulation curves of environmental DNA sequences predict coastal fish diversity in the coral triangle

Environmental DNA (eDNA) has the potential to provide more comprehensive biodiversity assessments, particularly for vertebrates in species-rich regions. However, this method requires the completeness of a reference database (i.e. a list of DNA sequences attached to each species), which is not currently achieved for many taxa and ecosystems. As an alternative, a range of operational taxonomic units (OTUs) can be extracted from eDNA metabarcoding. However, the extent to which the diversity of OTUs provided by a limited eDNA sampling effort can predict regional species diversity is unknown. Here, by modelling OTU accumulation curves of eDNA seawater samples across the Coral Triangle, we obtained an asymptote reaching 1531 fish OTUs, while 1611 fish species are recorded in the region. We also accurately predict (R² = 0.92) the distribution of species richness among fish families from OTU-based asymptotes. Thus, the multi-model framework of OTU accumulation curves extends the use of eDNA metabarcoding in ecology, biogeography and conservation.

Temporal lags in observed and dark diversity in the Anthropocene

Understanding biodiversity changes in the Anthropocene (e.g. due to climate and land-use change) is an urgent ecological issue. This important task is challenging because global change effects and species responses are dependent on the spatial scales considered. Furthermore, responses are often not immediate. However, both scale and time delay issues can be tackled when, at each study site, we consider dynamics in both observed and dark diversity. Dark diversity includes those species in the region that can potentially establish and thrive in the local sites' conditions but are currently locally absent. Effectively, dark diversity connects biodiversity at the study site to the regional scales and defines the site-specific species pool (observed and dark diversity together). With dark diversity, it is possible to decompose species gains and losses into two space-related components: one associated with local dynamics (species moving from observed to dark diversity and vice versa) and another related to gains and losses of site-specific species pool (species moving to and from the pool after regional immigration, regional extinction or change in local ecological conditions). Extinction debt and immigration credit are useful to understand dynamics in observed diversity, but delays might happen in species pool changes as well. In this opinion piece we suggest that considering both observed and dark diversity and their temporal dynamics provides a deeper understanding of biodiversity changes. Considering both observed and dark diversity creates opportunities to improve conservation by allowing to identify species that are likely to go regionally extinct as well as foreseeing which of the species that newly arrive to the region are more likely to colonize local sites. Finally, by considering temporal lags and species gains and losses in observed and dark diversity, we combine phenomena at both spatial and temporal scales, providing a novel tool to examine biodiversity change in the Anthropocene.

Using vertebrate environmental DNA from seawater in biomonitoring of marine habitats

Conservation and management of marine biodiversity depends on biomonitoring of marine habitats, but current approaches are resource-intensive and require different approaches for different organisms. Environmental DNA (eDNA) extracted from water samples is an efficient and versatile approach to detecting aquatic animals. In the ocean, eDNA composition reflects local fauna at fine spatial scales, but little is known about the effectiveness of eDNA-based monitoring of marine communities at larger scales. We investigated the potential of eDNA to characterize and distinguish marine communities at large spatial scales by comparing vertebrate species composition among marine habitats in Qatar, the Arabian Gulf (also known as the Persian Gulf), based on eDNA metabarcoding of seawater samples. We conducted species accumulation analyses to estimate how much of the vertebrate diversity we detected. We obtained eDNA sequences from a diverse assemblage of marine vertebrates, spanning 191 taxa in 73 families. These included rare and endangered species and covered 36% of the bony fish genera previously recorded in the Gulf. Sites of similar habitat type were also similar in eDNA composition. The species accumulation analyses showed that the number of sample replicates was insufficient for some sampling sites but suggested that a few hundred eDNA samples could potentially capture >90% of the marine vertebrate diversity in the study area. Our results confirm that seawater samples contain habitat-characteristic molecular signatures and that eDNA monitoring can efficiently cover vertebrate diversity at scales relevant to national and regional conservation and management.

Dark diversity reveals importance of biotic resources and competition for plant diversity across habitats

Species richness is the most commonly used metric to quantify biodiversity. However, examining dark diversity, the group of missing species which can potentially inhabit a site, can provide a more thorough understanding of the processes influencing observed biodiversity and help evaluate the restoration potential of local habitats. So far, dark diversity has mainly been studied for specific habitats or large-scale landscapes, while less attention has been given to variation across broad environmental gradients or as a result of local conditions and biotic interactions. In this study, we investigate the importance of local environmental conditions in determining dark diversity and observed richness in plant communities across broad environmental gradients. Using the ecospace concept, we investigate how these biodiversity measures relate to abiotic gradients (defined as position), availability of biotic resources (defined as expansion), spatiotemporal extent of habitats (defined as continuity), and species interactions through competition. Position variables were important for both observed diversity and dark diversity, some with quadratic relationships, for example, plant richness showing a unimodal response to soil fertility corresponding to the intermediate productivity hypothesis. Interspecific competition represented by community mean Grime C had a negative effect on plant species richness. Besides position-related variables, organic carbon was the most important variable for dark diversity, indicating that in late-succession habitats such as forests and shrubs, dark diversity is generally low. The importance of highly competitive species indicates that intermediate disturbance, such as grazing, may facilitate higher species richness and lower dark diversity.

Environmental DNA can act as a biodiversity barometer of anthropogenic pressures in coastal ecosystems

Loss of biodiversity from lower to upper trophic levels reduces overall productivity and stability of coastal ecosystems in our oceans, but rarely are these changes documented across both time and space. The characterisation of environmental DNA (eDNA) from sediment and seawater using metabarcoding offers a powerful molecular lens to observe marine biota and provides a series of ‘snapshots’ across a broad spectrum of eukaryotic organisms. Using these next-generation tools and downstream analytical innovations including machine learning sequence assignment algorithms and co-occurrence network analyses, we examined how anthropogenic pressures may have impacted marine biodiversity on subtropical coral reefs in Okinawa, Japan. Based on 18 S ribosomal RNA, but not ITS2 sequence data due to inconsistent amplification for this marker, as well as proxies for anthropogenic disturbance, we show that eukaryotic richness at the family level significantly increases with medium and high levels of disturbance. This change in richness coincides with compositional changes, a decrease in connectedness among taxa, an increase in fragmentation of taxon co-occurrence networks, and a shift in indicator taxa. Taken together, these findings demonstrate the ability of eDNA to act as a barometer of disturbance and provide an exemplar of how biotic networks and coral reefs may be impacted by anthropogenic activities.

Assessing species diversity of Coral Triangle artisanal fisheries: A DNA barcode reference library for the shore fishes retailed at Ambon harbor (Indonesia)

The Coral Triangle (CT), a region spanning across Indonesia and Philippines, is home to about 4,350 marine fish species and is among the world's most emblematic regions in terms of conservation. Threatened by overfishing and oceans warming, the CT fisheries have faced drastic declines over the last decades. Usually monitored through a biomass-based approach, fisheries trends have rarely been characterized at the species level due to the high number of taxa involved and the difficulty to accurately and routinely identify individuals to the species level. Biomass, however, is a poor proxy of species richness, and automated methods of species identification are required to move beyond biomass-based approaches. Recent meta-analyses have demonstrated that species richness peaks at intermediary levels of biomass. Consequently, preserving biomass is not equal to preserving biodiversity. We present the results of a survey to estimate the shore fish diversity retailed at the harbor of Ambon Island, an island located at the center of the CT that display exceptionally high biomass despite high levels of threat, while building a DNA barcode reference library of CT shore fishes targeted by artisanal fisheries. We sampled 1,187 specimens and successfully barcoded 696 of the 760 selected specimens that represent 202 species. Our results show that DNA barcodes were effective in capturing species boundaries for 96% of the species examined, which opens new perspectives for the routine monitoring of the CT fisheries.

Effects of exposure to large sharks on the abundance and behavior of mobile prey fishes along a temperate coastal gradient

Top predators can exert strong influences on community structure and function, both via direct, consumptive effects, as well as through non-consumptive, fear-based effects (i.e. predation risk). However, these effects are challenging to quantify, particularly for mobile predators in marine ecosystems. To advance this field of research, here we used baited remote underwater video stations (BRUVs) to assess how the behavior of mobile fish species off Cape Cod, Massachusetts, was affected by exposure to large sharks. We categorized sites into three levels of differential shark predation exposure (white sharks, Carcharodon carcharias) and quantified the relative abundance and arrival times (elapsed time before appearing on screen) for six mobile fish prey groups to the BRUV stations. Increased large shark exposure was associated with a decrease in overall prey abundance, but the overall response was prey group-specific. Foraging of smooth dogfish, a likely important prey item for large sharks in the system, was significantly reduced in areas frequented by white sharks. Specifically, the predicted probabilities of smooth dogfish bait contacts or bite attempts occurring were reduced by factors of 5.7 and 8.4, respectively, in areas of high exposure as compared to low exposure. These modifications were underscored by a decrease in smooth dogfish abundance in areas of high exposure as well. Our results suggest that populations of large, roving sharks may induce food-related costs in prey. We discuss the implications of this work within the context of the control of risk (COR) hypothesis, for the purposes of advancing our understanding of the ecological role and effects of large sharks on coastal marine ecosystems.

eDNA metabarcoding survey reveals fine-scale coral reef community variation across a remote, tropical island ecosystem

Environmental DNA (eDNA) metabarcoding, a technique for retrieving multispecies DNA from environmental samples, can detect a diverse array of marine species from filtered seawater samples. There is a growing potential to integrate eDNA alongside existing monitoring methods in order to establish or improve the assessment of species diversity. Remote island reefs are increasingly vulnerable to climate-related threats and as such there is a pressing need for cost-effective whole-ecosystem surveying to baseline biodiversity, study assemblage changes and ultimately develop sustainable management plans. We investigated the utility of eDNA metabarcoding as a high-resolution, multitrophic biomonitoring tool at the Cocos (Keeling) Islands, Australia (CKI)-a remote tropical coral reef atoll situated within the eastern Indian Ocean. Metabarcoding assays targeting the mitochondrial 16S rRNA and CO1 genes, as well as the 18S rRNA nuclear gene, were applied to 252 surface seawater samples collected from 42 sites within a 140 km² area. Our assays successfully detected a wide range of bony fish and elasmobranchs (244 taxa), crustaceans (88), molluscs (37) and echinoderms (7). Assemblage composition varied significantly between sites, reflecting habitat partitioning across the island ecosystem and demonstrating the localisation of eDNA signals, despite extensive tidal and oceanic movements. In addition, we document putative new occurrence records for 46 taxa and compare the efficiency of our eDNA approach to visual survey techniques at CKI. Our study demonstrates the utility of a multimarker metabarcoding approach in capturing multitrophic biodiversity across an entire coral reef atoll and sets an important baseline for ongoing monitoring and management.

Marine environmental DNA: Approaches, applications, and opportunities

Environmental DNA (eDNA) is increasingly being used to document species distributions and habitat use in marine systems, with much of the recent effort focused on leveraging advances in next-generation DNA sequencing to assess and track biodiversity across taxonomic groups. Environmental DNA offers a number of important advantages over traditional survey techniques, including non-invasive sampling, sampling where traditional approaches are impractical or inefficient (e.g. deep oceans), reduced cost, and increased detection sensitivity. However, eDNA applications are currently limited because of an insufficient understanding of the influence of sample source, analytical approach, and marker type on eDNA detections. Because approaches vary considerably among eDNA studies, we present a summary of the current state of the field and emerging best practices. The impact of observed variation in rates of eDNA production, persistence, and transport are also discussed and future research needs are highlighted with the goal of expanding eDNA applications, including the development of statistical models to improve the predictability of eDNA detection and quantification.

Zooplankton biodiversity monitoring in polluted freshwater ecosystems: A technical review

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

•

Freshwater ecosystems and associated biodiversity have been highly degraded.

•

Biodiversity monitoring is crucial for diagnosing degradation degrees.

•

Here we review available methods for monitoring zooplankton biodiversity.

•

We propose possible solutions for existing technical issues.

Biodiversity assessment of tropical shelf eukaryotic communities via pelagic eDNA metabarcoding

Our understanding of marine communities and their functions in an ecosystem relies on the ability to detect and monitor species distributions and abundances. Currently, the use of environmental DNA (eDNA) metabarcoding is increasingly being applied for the rapid assessment and monitoring of aquatic species. Most eDNA metabarcoding studies have either focussed on the simultaneous identification of a few specific taxa/groups or have been limited in geographical scope. Here, we employed eDNA metabarcoding to compare beta diversity patterns of complex pelagic marine communities in tropical coastal shelf habitats spanning the whole Caribbean Sea. We screened 68 water samples using a universal eukaryotic COI barcode region and detected highly diverse communities, which varied significantly among locations, and proved good descriptors of habitat type and environmental conditions. Less than 15% of eukaryotic taxa were assigned to metazoans, most DNA sequences belonged to a variety of planktonic "protists," with over 50% of taxa unassigned at the phylum level, suggesting that the sampled communities host an astonishing amount of micro-eukaryotic diversity yet undescribed or absent from COI reference databases. Although such a predominance of micro-eukaryotes severely reduces the efficiency of universal COI markers to investigate vertebrate and other metazoans from aqueous eDNA, the study contributes to the advancement of rapid biomonitoring methods and brings us closer to a full inventory of extant marine biodiversity.

A practical guide to sample preservation and pre-PCR processing of aquatic environmental DNA

Environmental DNA (eDNA) is rapidly growing in popularity as a tool for community assessments and species detection. While eDNA approaches are now widely applied, there is not yet agreement on best practices for sample collection and processing. Investigators looking to integrate eDNA approaches into their research programme are required to examine a growing collection of disparate studies to make an often uncertain decision about which protocols best fit their needs. To promote the application of eDNA approaches and to encourage the generation of high-quality data, here we review the most common techniques for the collection, preservation and extraction of metazoan eDNA from water samples. Specifically, we focus on experimental studies that compare various methods and outline the numerous challenges associated with eDNA. While the diverse applications of eDNA do not lend themselves to a one-size-fits-all recommendation, in most cases, capture/concentration of eDNA on cellulose nitrate filters (with pore size determined by water turbidity), followed by storage of filters in Longmire's buffer and extraction with a DNeasy Blood & Tissue Kit (or similar) has been shown to provide sufficient, high-quality DNA. However, we also emphasize the importance of testing and optimizing protocols for the system of interest.

Particle Size Distribution of Environmental DNA from the Nuclei of Marine Fish

Environmental DNA (eDNA) analyses have enabled a more efficient surveillance of species distribution and composition than conventional methods. However, the characteristics and dynamics of eDNA (e.g., origin, state, transport, and fate) remain unknown. This is especially limited for the eDNA derived from nuclei (nu-eDNA), which has recently been used in eDNA analyses. Here, we compared the particle size distribution (PSD) of nu-eDNA from Japanese Jack Mackerel (Trachurus japonicus) with that of mt-eDNA (eDNA derived from mitochondria) reported in previous studies. We repeatedly sampled rearing water from the tanks under multiple temperatures and fish biomass levels, and quantified the copy numbers of size-fractioned nu-eDNA. We found that the concentration of nu-eDNA was higher than that of mt-eDNA at 3-10 μm size fraction. Moreover, at the 0.8-3 μm and 0.4-0.8 μm size fractions, eDNA concentrations of both types increased with higher temperature and their degradation tended to be suppressed. These results imply that the production of eDNA from large to small size fractions could buffer the degradation of small-sized eDNA, which could improve its persistence in water. Our findings will contribute to refine the difference between nu- and mt-eDNA properties, and assist eDNA analyses as an efficient tool for the conservation of aquatic species.