Environmental DNA analysis shows high potential as a tool for estimating intraspecific genetic diversity in a wild fish population

New insights into biologic interpretation of bioinformatic pipelines for fish eDNA metabarcoding: A case study in Pearl River estuary

Environmental DNA (eDNA) metabarcoding is an emerging tool for monitoring biological communities in aquatic ecosystems. The selection of bioinformatic pipelines significantly impacts the results of biodiversity assessments. However, there is currently no consensus on the appropriate bioinformatic pipelines for fish community analysis in eDNA metabarcoding. In this study, we compared three bioinformatic pipelines (Uparse, DADA2, and UNOISE3) using real and mock (constructed with 15/30 known fish) communities to investigate the differences in biological interpretation during the data analysis process in eDNA metabarcoding. Performance evaluation and diversity analyses revealed that the choice of bioinformatic pipeline could impact the biological results of metabarcoding experiments. Among the three pipelines, the operational taxonomic units (OTU)-based pipeline (Uparse) showed the best performance (sensitivity: 0.6250 ± 0.0166; compositional similarity: 0.4000 ± 0.0571), the highest richness (25-102) and minimal inter-group differences in alpha diversity. It suggested the OTU-based pipeline possessed superior capability in fish diversity monitoring compared to ASV/ZOTU-based pipeline. Additionally, the Bray-Curtis distance matrix achieved the highest discriminative effect in the PCoA (43.3%-53.89%) and inter-group analysis (P < 0.01), indicating it was better at distinguishing compositional differences or specific genera of fish community at different sampling sites than other distance matrices. These findings provide new insights into fish community monitoring through eDNA metabarcoding in estuarine environments.

Haplotype-level metabarcoding of freshwater macroinvertebrate species: A prospective tool for population genetic analysis

Metabarcoding is a molecular-based tool capable of large quantity high-throughput species identification from bulk samples that is a faster and more cost-effective alternative to conventional DNA-sequencing approaches. Still, further exploration and assessment of the laboratory and bioinformatics strategies are required to unlock the potential of metabarcoding-based inference of haplotype information. In this study, we assessed the inference of freshwater macroinvertebrate haplotypes from metabarcoding data in a mock sample. We also examined the influence of DNA template concentration and PCR cycle on detecting true and spurious haplotypes. We tested this strategy on a mock sample containing twenty individuals from four species with known haplotypes based on the 658-bp Folmer region of the mitochondrial cytochrome c oxidase gene. We recovered fourteen zero-radius operational taxonomic units (zOTUs) of 421-bp length, with twelve zOTUs having a 100% match with the Sanger haplotype sequences. High-quality reads relatively increased with increasing PCR cycles, and the relative abundance of each zOTU was consistent for each cycle. This suggests that increasing the PCR cycles from 24 to 64 did not affect the relative abundance of each zOTU. As metabarcoding becomes more established and laboratory protocols and bioinformatic pipelines are continuously being developed, our study demonstrated the method's ability to infer intraspecific variability while highlighting the challenges that must be addressed before its eventual application for population genetic studies.

Unsupervised biological integrity assessment by eDNA biomonitoring of multi-trophic aquatic taxa

The biological integrity of global freshwater ecosystems is threatened by ever-increasing environmental stressors due to increased human activities, such as land-use change, eutrophication, toxic pollutants, overfishing, and exploitation. Traditional ecological assessments of lake or riverine ecosystems often require human supervision of a pre-selected reference area, using the current state of the reference area as the expected state. However, selecting an appropriate reference area has become increasingly difficult with the expansion of human activities. Here, an unsupervised biological integrity assessment framework based on environmental DNA metabarcoding without a prior reference area is proposed. Taxon richness, species dominance, co-occurrence network density, and phylogenetic distance were used to assess the aquatic communities in the Taihu Lake basin. Multi-gene metabarcoding revealed comprehensive biodiversity at multiple trophic levels including algae, protists, zooplankton, and fish. Fish sequences were mainly derived from 12S, zooplankton mainly from mitochondrial cytochrome C oxidase subunit I, and algae and protists mainly from 18S. There were significant differences in community composition among lakes, rivers, and reservoirs but no significant differences in the four fundamental biological indicators. The algal and zooplankton integrities were positively correlated with protist and fish integrities, respectively. Additionally, the algal integrity of lakes was found to be significantly lower than that of rivers. The unsupervised assessment framework proposed in this study allows different ecosystems, including the same ecosystem in different seasons, to adopt the same indicators and assessment methods, which is more convenient for environmental management and decision-making.

Estimating number of European eel (Anguilla anguilla) individuals using environmental DNA and haplotype count in small rivers

Knowledge about population genetic data is important for effective conservation management. Genetic research traditionally requires sampling directly from the organism, for example tissue, which can be challenging, time-consuming, and harmful to the animal. Environmental DNA (eDNA) approaches offer a way to sample genetic material noninvasively. In attempts to estimate population size of aquatic species using eDNA, researchers have found positive correlations between biomass and eDNA concentrations, but the approach is debated because of variations in the production and degrading of DNA in water. Recently, a more accurate eDNA-approach has emerged, focusing on the genomic differences between individuals. In this study, we used eDNA from water samples to estimate the number of European eel (Anguilla anguilla) individuals by examining haplotypes in the mitochondrial D-loop region, both in a closed aquatic environment with 10 eels of known haplotypes and in three rivers. The results revealed that it was possible to find every eel haplotype in the eDNA sample collected from the closed environment. We also found 13 unique haplotypes in the eDNA samples from the three rivers, which probably represent 13 eel individuals. This means that it is possible to obtain genomic information from European eel eDNA in water; however, more research is needed to develop the approach into a possible future tool for population quantification.

Environmental DNA study on aquatic ecosystem monitoring and management: Recent advances and prospects

Environmental DNA (eDNA) is organismal DNA that can be detected in the environment and is derived from cellular material of organisms shed into aquatic or terrestrial environments. It can be sampled and monitored using molecular methods, which is important for the early detection of invasive and native species as well as the discovery of rare and cryptic species. While few reviews have summarized the latest findings on eDNA for most aquatic animal categories in the aquatic ecosystem, especially for aquatic eDNA processing and application. In the present review, we first performed a bibliometric network analysis of eDNA studies on aquatic animals. Subsequently, we summarized the abiotic and biotic factors affecting aquatic eDNA occurrence. We also systematically discussed the relevant experiments and analyses of aquatic eDNA from various aquatic organisms, including fish, molluscans, crustaceans, amphibians, and reptiles. Subsequently, we discussed the major achievements of eDNA application in studies on the aquatic ecosystem and environment. The application of eDNA will provide an entirely new paradigm for biodiversity conservation, environment monitoring, and aquatic species management at a global scale.

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.

Environmental DNA analysis for macro-organisms: species distribution and more

In an era of severe biodiversity loss, biological monitoring is becoming increasingly essential. The analysis of environmental DNA (eDNA) has emerged as a new approach that could revolutionize the biological monitoring of aquatic ecosystems. Over the past decade, macro-organismal eDNA analysis has undergone significant developments and is rapidly becoming established as the golden standard for non-destructive and non-invasive biological monitoring. In this review, I summarize the development of macro-organismal eDNA analysis to date and the techniques used in this field. I also discuss the future perspective of these analytical methods in combination with sophisticated analytical techniques for DNA research developed in the fields of molecular biology and molecular genetics, including genomics, epigenomics, and single-cell technologies. eDNA analysis, which to date has been used primarily for determining the distribution of organisms, is expected to develop into a tool for elucidating the physiological state and behaviour of organisms. The fusion of microbiology and macrobiology through an amalgamation of these technologies is anticipated to lead to the future development of an integrated biology.

eDNA biomonitoring revealed the ecological effects of water diversion projects between Yangtze River and Tai Lake

Water diversion has been widely used to address water shortages and security issues. However, its long-term ecological impacts, particularly on the biodiversity and structure of the local community, have often been neglected due to limitations of conventional biomonitoring. Taking the water diversion projects from Yangtze River to Tai Lake (WDYT) as examples, environmental DNA (eDNA) metabarcoding was used to investigate the potential ecological impact of water diversion on the connected basins. Firstly, 136 phytoplankton genera/species, including 31 cyanobacteria and 105 eukaryotic phytoplankton (Euk-phytoplankton), were identified from 26 sites by metabarcoding of 16S rDNA V3 and 18S rDNA V9 regions. eDNA metabarcoding showed an obvious advantage in detecting nano/pico-plankton (< 20 μm in size) compared with the morphological approach. Secondly, more shared taxa and higher similarity of community composition were observed in Gonghu Bay/Zhushan Bay with its connected river than with the center of Tai Lake, indicating that water diversions were accelerating the biotic homogenization between different waterbodies. Skeletonema potamos, the native species of Yangtze River (4.04% of the total Euk-phytoplankton reads) was detected in different connecting regions of Tai Lake (0.03%-0.54% of the total Euk-phytoplankton reads), where its relative abundance was consistent with the influence of water diversion from Yangtze River. Furthermore, the introduction of S. potamos significantly affected the local community compositions of phytoplankton in Tai Lake. Finally, the ecological effect (e.g., taxa richness, community composition and species invasion) of the WDYT on phytoplankton in the west of Tai Lake was more significant than that in the east, which was consistent with the scale (volume and duration) of the water diversion projects. Overall, this study highlights the value of eDNA biomonitoring in the ecological impact assessment of water transfer projects.

Estimating the genetic diversity of Pacific salmon and trout using multigene eDNA metabarcoding

Genetic diversity underpins species conservation and management goals, and ultimately determines a species' ability to adapt. Using freshwater environmental DNA (eDNA) samples, we examined mitochondrial genetic diversity using multigene metabarcode sequence data from four Oncorhynchus species across 16 sites in Oregon and northern California. Our multigene metabarcode panel included targets commonly used in population genetic NADH dehydrogenase 2 (ND2), phylogenetic cytochrome c oxidase subunit 1 (COI) and eDNA (12S ribosomal DNA) screening. The ND2 locus showed the greatest within-species haplotype diversity for all species, followed by COI and then 12S rDNA for all species except Oncorhynchus kisutch. Sequences recovered for O. clarkii clarkii were either identical to, or one mutation different from, previously characterized haplotypes (95.3% and 4.5% of reads, respectively). The greatest diversity in O. c. clarkii was among coastal watersheds, and subsets of this diversity were shared with fish in inland watersheds. However, coastal streams and the Umpqua River watershed appear to harbour unique haplotypes. Sequences from O. mykiss revealed a disjunction between the Willamette watershed and southern watersheds suggesting divergent histories. We also identified similarities between populations in the northern Deschutes and southern Klamath watersheds, consistent with previously hypothesized connections between the two via inland basins. Oncorhynchus kisutch was only identified in coastal streams and the Klamath River watershed, with most diversity concentrated in the coastal Coquille watershed. Oncorhynchus tshawytscha was only observed at one site, but contained multiple haplotypes at each locus. The characterization of genetic diversity at multiple loci expands the knowledge gained from eDNA sampling and provides crucial information for conservation actions and genetic management.

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.

Water pre-filtration methods to improve environmental DNA detection by real-time PCR and metabarcoding

Environmental DNA (eDNA) analysis is a novel approach for biomonitoring and has been mostly used in clear water. It is difficult to detect eDNA in turbid water as filter clogging occurs, and environmental samples contain various substances that inhibit the polymerase chain reaction (PCR) and affect the accuracy of eDNA analysis. Therefore, we applied a pre-filtration method to better detect the fish species (particularly pale chub, Opsariichthys platypus) present in a water body by measuring eDNA in environmental samples containing PCR inhibitors. Upon conducting 12S rRNA metabarcoding analysis (MiFish), we found that pre-filtration did not affect the number or identities of fish species detected in our samples, but pre-filtration through pore sizes resulted in significantly reduced variance among replicate samples. Additionally, PCR amplification was improved by the pre-filtration of environmental samples containing PCR inhibitors such as humic substances. Although this study may appear to be a conservative and ancillary experiment, pre-filtration is a simple technique that can not only improve the physical properties of water, such as turbidity, but also the quality of eDNA biomonitoring.

High-throughput sequencing on preservative ethanol is effective at jointly examining infraspecific and taxonomic diversity, although bioinformatics pipelines do not perform equally

High-throughput sequencing of amplicons (HTSA) has been proposed as an effective approach to evaluate taxonomic and genetic diversity at the same time. However, there are still uncertainties as to how the results produced by different bioinformatics treatments impact the conclusions drawn on biodiversity and population genetics indices.We evaluated the ability of six bioinformatics pipelines to recover taxonomic and genetic diversity from HTSA data obtained from controlled assemblages. To that end, 20 assemblages were produced using 354 colonies of Botrylloides spp., sampled in the wild in ten marinas around Brittany (France). We used DNA extracted from preservative ethanol (ebDNA) after various time of storage (3, 6, and 12 months), and from a bulk of preserved specimens (bulkDNA). DNA was amplified with primers designed for targeting this ascidian genus. Results obtained from HTSA data were compared with Sanger sequencing on individual zooids (i.e., individual barcoding).Species identification and relative abundance determined with HTSA data from either ebDNA or bulkDNA were similar to those obtained with traditional individual barcoding. However, after 12 months of storage, the correlation between HTSA and individual-based data was lower than after shorter durations. The six bioinformatics pipelines were able to depict accurately the genetic diversity using standard population genetics indices (HS and FST), despite producing false positives and missing rare haplotypes. However, they did not perform equally and dada2 was the only pipeline able to retrieve all expected haplotypes.This study showed that ebDNA is a nondestructive alternative for both species identification and haplotype recovery, providing storage does not last more than 6 months before DNA extraction. Choosing the bioinformatics pipeline is a matter of compromise, aiming to retrieve all true haplotypes while avoiding false positives. We here recommend to process HTSA data using dada2, including a chimera-removal step. Even if the possibility to use multiplexed primer sets deserves further investigation to expand the taxonomic coverage in future similar studies, we showed that primers targeting a particular genus allowed to reliably analyze this genus within a complex community.

Estimation of tuna population by the improved analytical pipeline of unique molecular identifier-assisted HaCeD-Seq (haplotype count from eDNA)

Many studies have investigated the ability to identify species from environmental DNA (eDNA). However, even when individual species are identified, the accurate estimation of their abundances by traditional eDNA analyses has been still difficult. We previously developed a novel analytical method called HaCeD-Seq (Haplotype Count from eDNA), which focuses on the mitochondrial D-loop sequence. The D-loop is a rapidly evolving sequence and has been used to estimate the abundance of eel species in breeding water. In the current study, we have further improved this method by applying unique molecular identifier (UMI) tags, which eliminate the PCR and sequencing errors and extend the detection range by an order of magnitude. Based on this improved HaCeD-Seq pipeline, we computed the abundance of Pacific bluefin tuna (Thunnus orientalis) in aquarium tanks at the Tokyo Sea Life Park (Kasai, Tokyo, Japan). This tuna species is commercially important but is at high risk of resource depletion. With the developed UMI tag method, 90 out of 96 haplotypes (94%) were successfully detected from Pacific bluefin tuna eDNA. By contrast, only 29 out of 96 haplotypes (30%) were detected when UMI tags were not used. Our findings indicate the potential for conducting non-invasive fish stock surveys by sampling eDNA.

Fingerprinting of hatchery haplotypes and acquisition of genetic information by whole-mitogenome sequencing of masu salmon, Oncorhynchus masou masou, in the Kase River system, Japan

Stocking hatchery fish can lead to disturbance and extinction of the local indigenous population. Masu salmon Oncorhynchus masou masou, which is endemic across Japan, is a commonly stocked fish for recreational fishing in Japan. To conserve the indigenous resource, their genetic information is required, however, especially on Kyushu Island, the paucity of genetic information for this species has hindered proper resource management. Here, to identify hatchery mitogenome haplotypes of this species, stocked in the Kase River system, Kyushu Island, Japan, and to provide mitogenomic information for the resource management of this species, we analyzed the whole-mitogenome of masu salmon in this river system and several hatcheries potentially used for stocking. Whole-mitogenome sequencing clearly identified hatchery haplotypes, like fingerprints: among the 21 whole-mitogenome haplotypes obtained, six were determined to be hatchery haplotypes. These hatchery haplotypes were distributed in 13 out of 17 sites, suggesting that informal stocking of O. m. masou has been performed widely across this river system. The population of no hatchery haplotypes mainly belonged to clade I, a clade not found in Hokkaido Island in previous studies. Sites without hatchery haplotypes, and the non-hatchery haplotypes in clade I might be candidates for conservation as putative indigenous resources. The whole-mitogenome haplotype analysis also clarified that the same reared strain was used in multiple hatcheries. Analysis of molecular variance suggested that stocked hatchery haplotypes reduce the genetic variation among populations in this river system. It will be necessary to pay attention to genetic fluctuations so that the resources of this river system will not deteriorate further. The single nucleotide polymorphism data obtained here could be used for resource management in this and other rivers: e.g., for monitoring of informal stocking and stocked hatchery fishes, and/or putative indigenous resources.