MiFish, a set of universal PCR primers for metabarcoding environmental DNA from fishes: detection of more than 230 subtropical marine species

Aquatic environmental DNA: A review of the macro-organismal biomonitoring revolution

Environmental DNA (eDNA) is the fastest growing biomonitoring tool fuelled by two key features: time efficiency and sensitivity. Technological advancements allow rapid biodiversity detection at both species and community levels with increasing accuracy. Concurrently, there has been a global demand to standardise eDNA methods, but this is only possible with an in-depth overview of the technological advancements and a discussion of the pros and cons of available methods. We therefore conducted a systematic literature review of 407 peer-reviewed papers on aquatic eDNA published between 2012 and 2021. We observed a gradual increase in the annual number of publications from four (2012) to 28 (2018), followed by a rapid growth to 124 publications in 2021. This was mirrored by a tremendous diversification of methods in all aspects of the eDNA workflow. For example, in 2012 only freezing was applied to preserve filter samples, whereas we recorded 12 different preservation methods in the 2021 literature. Despite an ongoing standardisation debate in the eDNA community, the field is seemingly moving fast in the opposite direction and we discuss the reasons and implications. Moreover, by compiling the largest PCR-primer database to date, we provide information on 522 and 141 published species-specific and metabarcoding primers targeting a wide range of aquatic organisms. This works as a user-friendly 'distillation' of primer information that was hitherto scattered across hundreds of papers, but the list also reflects which taxa are commonly studied with eDNA technology in aquatic environments such as fish and amphibians, and reveals that groups such as corals, plankton and algae are under-studied. Efforts to improve sampling and extraction methods, primer specificity and reference databases are crucial to capture these ecologically important taxa in future eDNA biomonitoring surveys. In a rapidly diversifying field, this review synthetises aquatic eDNA procedures and can guide eDNA users towards best practice.

Investigations of Fish Assemblages Using Two Methods in Three Terminal Reservoirs of the East Route of South-to-North Water Transfer Project, China

The terminal reservoirs of water transfer projects directly supply water for domestic, agricultural, and industrial applications, and the water quality of these reservoirs produce crucial effects on the achievement of project targets. Typically, fish assemblages are monitored as indicators of reservoir water quality, and can also be regulated for its improvement. In the present study, we compared traditional fish landing (TFL) and environmental DNA (eDNA) metabarcoding methods for monitoring fish assemblages in three terminal reservoirs of the East Route of the South-to-North Water Transfer Project, China. Results of TFL and eDNA showed similar assemblage structures and patterns of diversity and spatial distribution with obvious differences in fish composition across three examined reservoirs. Demersal and small fish were dominant in all reservoirs. In addition, a strong association between water transfer distance and assemblages and distribution of non-native fish was found. Our findings highlight the necessity of the fish assemblage monitoring and managing for water quality and revealed the impact of water diversion distance on the structure of fish assemblages and dispersal of alien species along the water transfer project.

DNA metabarcoding to assess prey overlap between tuna and seabirds in the Eastern tropical Atlantic: Implications for an ecosystem-based management

Overfishing has been drastically changing food webs in marine ecosystems, and it is pivotal to quantify these changes at the ecosystem level. This is especially important for ecosystems with a high diversity of top predators such as the Eastern Atlantic marine region. In this work we used high-throughput sequencing methods to describe the diet of the two most abundant tuna species, the Skipjack tuna (Katsuwonus pelamis) and the Yellowfin tuna (Thunnus albacares), highly targeted by fisheries off west Africa. We also explored prey diversity overlap between these tuna species and the seabird species breeding in Cabo Verde that are most likely to share prey preferences and suffer from bycatch, the Brown booby (Sula leucogaster) and Cape Verde shearwater (Calonectris edwardsii). Overall, the diet of both tuna species was more diverse than that of seabirds. Skipjack tuna diet was dominated by prey from lower trophic levels, such as krill, anchovies, and siphonophores, while the Yellowfin tuna diet was mainly based on epipelagic fish such as flying and halfbeak fishes. Some of the most abundant prey families detected in the Yellowfin tuna diet were shared with both seabird species, resulting in a high prey diversity overlap between this tuna species and seabirds These results have implications for the management of tuna fisheries in the Eastern Tropical Atlantic, because a large decrease of both tuna species might have cascading effects on both primary and secondary consumer levels, and the decrease of these underwater predators may have implications on the viability of tropical seabird populations.

The challenge of implementing environmental DNA metabarcoding to detect elasmobranchs in a resource-limited marine protected area

Elasmobranchs are threatened and eDNA metabarcoding is a powerful tool that can help efforts to better understand and conserve them. Nevertheless, the inter-calibration between optimal methodological practices and its implementation in resource-limited situations is still an issue. Based on promising results from recent studies, the authors applied a cost-effective protocol with parameters that could be easily replicated by any conservationist. Nonetheless, the results with fewer elasmobranchs detected than expected reveal that endorsed primers and sampling strategies still require further optimization, especially for applications in resource-limited conservation programmes.

Decoding dissolved information: environmental DNA sequencing at global scale to monitor a changing ocean

The use of environmental DNA (eDNA) technology for environmental monitoring is rapidly expanding, with applications for fisheries, coral reefs, harmful algal blooms, invasive and endangered species, and biodiversity monitoring. By enabling detection of species over space and time, eDNA fulfills a fundamental need of environmental surveys. Traditional surveys are expensive, require significant capital expenditure, and can be destructive; eDNA offers promise for cheaper, less invasive, and higher-resolution (i.e. genetic) assessments of environments and stocks. However, challenges in quantification, detection limits, biobanking capacity, reference databases, and data management and integration remain significant hurdles to efficient eDNA monitoring at global and decadal scale. Here, we consider the current state of eDNA technology and its suitability for the problems for which it is being used. We explore the current best practices, the logistical and social challenges that prevent eDNA from widespread adoption and benefit, and the emerging technologies that may address those challenges.

Fish eDNA detections in ports mirror fishing fleet activities and highlight the spread of non-indigenous species in the Mediterranean Sea

Gathering comprehensive marine biodiversity data can be difficult, costly and time consuming, preventing adequate knowledge of diversity patterns in many areas worldwide. We propose fishing ports as "natural" sinks of biodiversity information collected by fishing vessels probing disparate habitats, depths, and environments. By combining rapid environmental DNA metabarcoding (eDNA) surveys and data from public registers and Automatic Identification Systems, we show significant positive relationships between fishing fleet activities (i.e. fishing effort and characteristics of the fishing grounds) and the taxonomic fish assemblage composition in eleven Mediterranean fishing ports. Overall, we identified 160 fish and 123 invertebrate OTUs, including at least seven non-indigenous species, in some instances well beyond their known distribution areas. Our findings suggest that eDNA assessments of fishing harbours' waters might offer a rapid way to monitor marine biodiversity in unknown or under-sampled areas, as well as to reconstruct fishing catches, often underreported in several regions.

crabs-A software program to generate curated reference databases for metabarcoding sequencing data

The measurement of biodiversity is an integral aspect of life science research. With the establishment of second- and third-generation sequencing technologies, an increasing amount of metabarcoding data is being generated as we seek to describe the extent and patterns of biodiversity in multiple contexts. The reliability and accuracy of taxonomically assigning metabarcoding sequencing data have been shown to be critically influenced by the quality and completeness of reference databases. Custom, curated, eukaryotic reference databases, however, are scarce, as are the software programs for generating them. Here, we present crabs (Creating Reference databases for Amplicon-Based Sequencing), a software package to create custom reference databases for metabarcoding studies. crabs includes tools to download sequences from multiple online repositories (i.e., NCBI, BOLD, EMBL, MitoFish), retrieve amplicon regions through in silico PCR analysis and pairwise global alignments, curate the database through multiple filtering parameters (e.g., dereplication, sequence length, sequence quality, unresolved taxonomy, inclusion/exclusion filter), export the reference database in multiple formats for immediate use in taxonomy assignment software, and investigate the reference database through implemented visualizations for diversity, primer efficiency, reference sequence length, database completeness and taxonomic resolution. crabs is a versatile tool for generating curated reference databases of user-specified genetic markers to aid taxonomy assignment from metabarcoding sequencing data. crabs can be installed via docker and is available for download as a conda package and via GitHub (https://github.com/gjeunen/reference_database_creator).

Evaluating the sampling effort for the metabarcoding‐based detection of fish environmental DNA in the open ocean

Clarifying the effect of the sampling protocol on the detection of environmental DNA (eDNA) is essential for appropriately designing biodiversity research. However, technical issues influencing eDNA detection in the open ocean, which consists of water masses with varying environmental conditions, have not been thoroughly investigated. This study evaluated the sampling effort for the metabarcoding‐based detection of fish eDNA using replicate sampling with filters of different pore sizes (0.22 and 0.45 μm) in the subtropical and subarctic northwestern Pacific Ocean and Arctic Chukchi Sea. The asymptotic analysis predicted that the accumulation curves for detected taxa did not saturate in most cases, indicating that our sampling effort (7 or 8 replicates, corresponding to 10.5–40 L of filtration in total) was insufficient to fully assess the species diversity in the open ocean and that tens of replicates or a substantial filtration volume were required. The Jaccard dissimilarities between filtration replicates were comparable with those between the filter types at any site. In subtropical and subarctic sites, turnover dominated the dissimilarity, suggesting that the filter pore size had a negligible effect. In contrast, nestedness dominated the dissimilarity in the Chukchi Sea, implying that the 0.22 μm filter could collect a broader range of eDNA than the 0.45 μm filter. Therefore, the effect of filter selection on the collection of fish eDNA likely varies depending on the region. These findings highlight the highly stochastic nature of fish eDNA collection in the open ocean and the difficulty of standardizing the sampling protocol across various water masses.

Pilot study of a comprehensive resource estimation method from environmental DNA using universal D-loop amplification primers

Many studies have investigated the ability of environmental DNA (eDNA) to identify the species. However, when individual species are to be identified, accurate estimation of their abundance using traditional eDNA analyses is still difficult. We previously developed a novel analytical method called HaCeD-Seq (haplotype count from eDNA by sequencing), which focuses on the mitochondrial D-loop sequence for eels and tuna. In this study, universal D-loop primers were designed to enable the comprehensive detection of multiple fish species by a single sequence. To sequence the full-length D-loop with high accuracy, we performed nanopore sequencing with unique molecular identifiers (UMI). In addition, to determine the D-loop reference sequence, whole genome sequencing was performed with thin coverage, and complete mitochondrial genomes were determined. We developed a UMI-based Nanopore D-loop sequencing analysis pipeline and released it as open-source software. We detected 5 out of 15 species (33%) and 10 haplotypes out of 35 individuals (29%) among the detected species. This study demonstrates the possibility of comprehensively obtaining information related to population size from eDNA. In the future, this method can be used to improve the accuracy of fish resource estimation, which is currently highly dependent on fishing catches.

Mare-MAGE curated reference database of fish mitochondrial genes

Biodiversity assessment approaches based on molecular biology techniques such as metabarcoding, RAD-seq, or SnaPshot sequencing are increasingly applied in assessing marine and aquatic ecosystems. Here we present a new reference database for fish meta-barcoding based on mitochondrial genes. The Mare-MAGE database contains quality-checked sequences of the mitochondrial 12S ribosomal RNA and Cytochrome c Oxidase I gene. All sequences were obtained from the National Center for Biotechnology Information- GenBank (NBCI-GenBank), the European Nucleotide Archive (ENA), AquaGene Database and BOLD database, and have undergone intensive processing. They were checked for false annotations and non-target anomalies, according to the Integrated Taxonomic Information System (ITIS) and FishBase. The dataset is compiled in ARB-Home, FASTA and Qiime2 formats, and is publicly available from the Mare-MAGE database website ( http://mare-mage.weebly.com/ ). It includes altogether 231,333 COI and 12S rRNA gene sequences of fish, covering 19,506 species of 4,058 genera and 586 families.

MitoFish, MitoAnnotator, and MiFish Pipeline: Updates in 10 Years

MitoFish, MitoAnnotator, and MiFish Pipeline are comprehensive databases of fish mitochondrial genomes (mitogenomes), accurate annotation software of fish mitogenomes, and a web platform for metabarcoding analysis of fish mitochondrial environmental DNA (eDNA), respectively. The MitoFish Suite currently receives over 48,000 visits worldwide every year; however, the performance and usefulness of the online platforms can still be improved. Here, we present essential updates on these platforms, including an enrichment of the reference data sets, an enhanced searching function, substantially faster genome annotation and eDNA analysis with the denoising of sequencing errors, and a multisample comparative analysis function. These updates have made our platform more intuitive, effective, and reliable. These updated platforms are freely available at http://mitofish.aori.u-tokyo.ac.jp/.

New insights into the diversity of cryptobenthic Cirripectes blennies in the Mascarene Archipelago sampled using Autonomous Reef Monitoring Structures (ARMS)

Autonomous Reef Monitoring Structures (ARMS) are artificial mini-reefs designed for standardized sampling of sessile and small motile cryptobenthic organisms. ARMS are also effective for collecting small cryptobenthic fishes, such as the combtooth blennies of the genus Cirripectes. Recent studies discovered several Cirripectes species endemic to islands or archipelagos, in spite of the generally broad distributions of tropical and subtropical blennies. Thus, to evaluate the diversity and distribution of Cirripectes species in the Mascarene Archipelago, a little-studied region but an important biodiversity hotspot, complete mitochondrial genomes, and nuclear rhodopsin genes were sequenced for 39 specimens collected with ARMS deployed on outer reef slopes at Reunion and Rodrigues islands. Mitochondrial COI sequences were analyzed to integrate these specimens within the largest dataset of publicly available sequences. Three species were found in the Mascarene Archipelago, Cirripectes castaneus, Cirripectes randalli, and Cirripectes stigmaticus. C. castaneus and C. stigmaticus both have an Indo-Pacific distribution with several haplotypes shared among distant localities. In agreement with the literature, C. randalli shows a small-range endemism restricted to the Mascarenes. We confirmed the presence of C. castaneus, C. randalli, and C. stigmaticus in Rodrigues, and the presence of C. stigmaticus in Reunion. This study contributes to filling the gaps in taxonomic and molecular knowledge of the reef cryptobiome in the South-West Indian Ocean, and provides the first complete mitogenomes for the genus, a crucial step for future molecular-based inventories (e.g., eDNA).

Biodiversity and spatial distribution of ascidian using environmental DNA metabarcoding

Monitoring studies are necessary to understand the biodiversity of marine ecosystems and are useful for identifying and managing rare or invasive species. Because monitoring has traditionally relied only on visual surveys (e.g., trapping, netting, electrofishing, and SCUBA diving) with limited time and physical resources, environmental DNA (eDNA) analysis is being applied as an efficient monitoring method. This study compared whether the eDNA metabarcoding technique can replace the traditional visual survey in an ascidian fauna study. We designed ascidian-specific primers and identified a clear gap (3.75%) by barcoding gap analysis. Then, we collected seawater samples for eDNA analysis during the summer (August-September) of 2021 at three sites (Mokpo, Yeosu, and Uljin) in South Korea. In the survey sites of this study, 25 species were observed through literature and visual survey, among which 9 species were detected by metabarcoding and 16 species were not detected. On the other hand, 10 species were detected only by metabarcoding, and one of them was identified as Pyura mirabilis, an unrecorded species in South Korea. This study succeeded in detecting cryptic or rare species with one seawater collection, which can be used to determine their unexplored habitat. Therefore, we conclude that monitoring using eDNA is more efficient than visual surveys for detecting rare or cryptic ascidian species. We also suggest that, when combined with traditional monitoring methods, it could be a tool to complement ascidian fauna studies.

Development of two new sets of PCR primers for eDNA metabarcoding of brittle stars (Echinodermata, Ophiuroidea)

Brittle stars (class Ophiuroidea) are marine invertebrates comprising approximately 2,100 extant species, and are considered to constitute the most diverse taxon of the phylum Echinodermata. As a non-invasive method for monitoring biodiversity, we developed two new sets of PCR primers for metabarcoding environmental DNA (eDNA) from brittle stars. The new primer sets were designed to amplify 2 short regions of the mitochondrial 16S rRNA gene, comprising a conserved region (111–115 bp, 112 bp on average; named “16SOph1”) and a hyper-variable region (180–195 bp, 185 bp on average; named “16SOph2”) displaying interspecific variation. The performance of the primers was tested using eDNA obtained from two sources: a) rearing water of an 2.5 or 170 L aquarium tanks containing 15 brittle star species and b) from natural seawater collected around Misaki, the Pacific coast of central Japan, at depths ranging from shallow (2 m) to deep (> 200 m) sea. To build a reference library, we obtained 16S rRNA sequences of brittle star specimens collected from around Misaki and from similar depths in Japan, and sequences registered in International Nucleotide Sequence Database Collaboration. As a result of comparison of the obtained eDNA sequences with the reference library 37 (including cryptic species) and 26 brittle star species were detected with certain identities by 16SOph1 and 16SOph2 analyses, respectively. In shallow water, the number of species and reads other than the brittle stars detected with 16SOph1 was less than 10% of the total number. On the other hand, the number of brittle star species and reads detected with 16SOph2 was less than half of the total number, and the number of detected non-brittle star metazoan species ranged from 20 to 46 species across 6 to 8 phyla (only the reads at the “Tank” were less than 0.001%). The number of non-brittle star species and reads at 80 m was less than 10% with both of the primer sets. These findings suggest that 16SOph1 is specific to the brittle star and 16SOph2 is suitable for a variety of marine metazoans. It appears, however, that further optimization of primer sequences would still be necessary to avoid possible PCR dropouts from eDNA extracts. Moreover, a detailed elucidation of the brittle star fauna in the examined area, and the accurate identification of brittle star species in the current DNA databank is required.

Can the Eurasian otter (Lutra lutra) be used as an effective sampler of fish diversity? Using molecular assessment of otter diet to survey fish communities

The Eurasian otter Lutra lutra is a generalist carnivore that is widely distributed in many aquatic ecosystems. Based on its inherent attributes of opportunistic foraging behaviour and broad dietary range, it is naturally considered a potential sampler of the diversity of aquatic vertebrates. To test the ability and efficiency of otters as a diversity sampler, we used DNA metabarcoding to investigate the composition in vertebrates of the diet of otters that inhabit a forest stream area in northeast China. Twenty vertebrate prey taxa were detected in 98 otter spraints. Otter diet mainly comprised aquatic fishes (59.4%) and amphibians (39.0%). We also used traditional approaches to investigate fish communities at 60 sampling sites in the same area to determine the relationship between fish population composition in the environment and otter diet. The comparison revealed that 28 species of fish were distributed in this area, of which five are simultaneously detected in otter spraints. This indicates that molecular analysis of the diet of otters is not an ideal approach for investigating fish diversity, at least when using the 12SV5 primer pair. Based on a review of the available molecular research on otter diet, we conclude that the low species resolution may be due to the presence of many closely-related prey species in native habitats and lack of suitable barcodes. Considering the remarkable power of diet metabarcoding analysis in capturing elusive and rare species, it represents an approach that can compensate for the defects associated with fishing methods and we suggest that it can be used as an auxiliary means of measuring traditional fish diversity.

Pragmatic applications of DNA barcoding markers in identification of fish species – a review

DNA barcoding and mini barcoding involve Cytochrome Oxidase Subunit I (COI) gene in mitochondrial genome and is used for accurate identification of species and biodiversity. The basic goal of the current study is to develop a complete reference database of fishes. It also evaluates the applicability of COI gene to identify fish at the species level with other aspects i.e., as Kimura 2 parameter (K2P) distance. The mean observed length of the sequence was ranging between 500 to 700 base pairs for fish species in DNA barcoding and 80 to 650 base pairs for DNA mini barcoding. This method describes the status of known to unknown samples but it also facilitates the detection of previously un-sampled species at distinct level. So, mini-barcoding is a method focuses on the analysis of short-length DNA markers has been demonstrated to be effective for species identification of processed food containing degraded DNA. While DNA meta-barcoding refers to the automated identification of multiple species from a single bulk sample. The may contain entire organisms or a single environmental sample containing degraded DNA. Despite DNA barcoding, mini barcoding and meta-barcoding are efficient methods for species identification which are helpful in conservation and proper management of biodiversity. It aids researchers to take an account of genetic as well as evolutionary relationships by collecting their morphological, distributional and molecular data. Overall, this paper discusses DNA barcoding technology and how it has been used to various fish species, as well as its universality, adaptability, and novel approach to DNA-based species identification.

Detection of the Endangered Siamese Bat Catfish (Oreoglanis siamensis Smith, 1933) in Doi Inthanon National Park Using Environmental DNA

Siamese bat catfish (Oreoglanis siamensis Smith, 1993) has been listed as an endangered species, and its abundance has been severely declining due to habitat degradation and overfishing. To establish an appropriate management strategy, it is crucial to gain information about the distribution of this endangered species. As O. siamensis live under rocks in streams, detecting their presence is difficult. Recently, environmental DNA (eDNA)-based detection has been demonstrated to be a valid tool for monitoring rare species, such as O. siamensis. Therefore, this study developed an eDNA assay targeting a 160 bp fragment of the COI region to detect the presence of this species in its natural habitat. An amount of 300 mL of water samples (0.7 μm filtered) were collected from 15 sites in the Mae Klang sub-basin, where this fish species was visually detected at two locations. O. siamensis eDNA was detected at 12 of the 15 sites sampled with varying concentrations (0.71-20.27 copies/mL), including at the sites where this species was visually detected previously. The developed O. siamensis eDNA assay was shown to be effective for detecting the presence of this endangered species in the Klang Phat and Klang Rivers within the Doi Inthanon National Park.

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.

Toward quantitative metabarcoding

Amplicon-sequence data from environmental DNA (eDNA) and microbiome studies provide important information for ecology, conservation, management, and health. At present, amplicon-sequencing studies-known also as metabarcoding studies, in which the primary data consist of targeted, amplified fragments of DNA sequenced from many taxa in a mixture-struggle to link genetic observations to the underlying biology in a quantitative way, but many applications require quantitative information about the taxa or systems under scrutiny. As metabarcoding studies proliferate in ecology, it becomes more important to develop ways to make them quantitative to ensure that their conclusions are adequately supported. Here we link previously disparate sets of techniques for making such data quantitative, showing that the underlying polymerase chain reaction mechanism explains the observed patterns of amplicon data in a general way. By modeling the process through which amplicon-sequence data arise, rather than transforming the data post hoc, we show how to estimate the starting DNA proportions from a mixture of many taxa. We illustrate how to calibrate the model using mock communities and apply the approach to simulated data and a series of empirical examples. Our approach opens the door to improve the use of metabarcoding data in a wide range of applications in ecology, public health, and related fields.

Quantitative monitoring of diverse fish communities on a large scale combining eDNA metabarcoding and qPCR

Environmental DNA (eDNA) metabarcoding is an effective method for studying fish communities but allows only an estimation of relative species abundance (density/biomass). Here, we combine metabarcoding with an estimation of the total abundance of eDNA amplified by our universal marker (teleo) using a quantitative (q)PCR approach to infer the absolute abundance of fish species. We carried out a 2850-km eDNA survey within the Danube catchment using a spatial integrative sampling protocol coupled with traditional electrofishing for fish biomass and density estimation. Total fish eDNA concentrations and total fish abundance were highly correlated. The correlation between eDNA concentrations per taxon and absolute specific abundance was of comparable strength when all sites were pooled and remained significant when the sites were considered separately. Furthermore, a nonlinear mixed model showed that species richness was underestimated when the amount of teleo-DNA extracted from a sample was below a threshold of 0.65 × 10⁶ copies of eDNA. This result, combined with the decrease in teleo-DNA concentration by several orders of magnitude with river size, highlights the need to increase sampling effort in large rivers. Our results provide a comprehensive description of longitudinal changes in fish communities and underline our combined metabarcoding/qPCR approach for biomonitoring and bioassessment surveys when a rough estimate of absolute species abundance is sufficient.

The global depth range of marine fishes and their genetic coverage for environmental DNA metabarcoding

The bathymetric and geographical distribution of marine species represent a key information in biodiversity conservation. Yet, deep-sea ecosystems are among the least explored on Earth and are increasingly impacted by human activities. Environmental DNA (eDNA) metabarcoding has emerged as a promising method to study fish biodiversity but applications to the deep-sea are still scarce. A major limitation in the application of eDNA metabarcoding is the incompleteness of species sequences available in public genetic databases which reduces the extent of detected species. This incompleteness by depth is still unknown. Here, we built the global bathymetric and geographical distribution of 10,826 actinopterygian and 960 chondrichthyan fish species. We assessed their genetic coverage by depth and by ocean for three main metabarcoding markers used in the literature: teleo and MiFish-U/E. We also estimated the number of primer mismatches per species amplified by in silico polymerase chain reaction which influence the probability of species detection. Actinopterygians show a stronger decrease in species richness with depth than Chondrichthyans. These richness gradients are accompanied by a continuous species turnover between depths. Fish species coverage with the MiFish-U/E markers is higher than with teleo while threatened species are more sequenced than the others. "Deep-endemic" species, those not ascending to the shallow depth layer, are less sequenced than not threatened species. The number of primer mismatches is not higher for deep-sea species than for shallower ones. eDNA metabarcoding is promising for species detection in the deep-sea to better account for the 3-dimensional structure of the ocean in marine biodiversity monitoring and conservation. However, we argue that sequencing efforts on "deep-endemic" species are needed.

Assessing the utility of marine filter feeders for environmental DNA (eDNA) biodiversity monitoring

Aquatic environmental DNA (eDNA) surveys are transforming how marine ecosystems are monitored. The time-consuming preprocessing step of active filtration, however, remains a bottleneck. Hence, new approaches that eliminate the need for active filtration are required. Filter-feeding invertebrates have been proven to collect eDNA, but side-by-side comparative studies to investigate the similarity between aquatic and filter-feeder eDNA signals are essential. Here, we investigated the differences among four eDNA sources (water; bivalve gill-tissue; sponges; and ethanol in which filter-feeding organisms were stored) along a vertically stratified transect in Doubtful Sound, New Zealand using three metabarcoding primer sets targeting fish and vertebrates. Combined, eDNA sources detected 59 vertebrates, while concurrent diver surveys observed eight fish species. There were no significant differences in alpha and beta diversity between water and sponge eDNA and both sources were highly correlated. Vertebrate eDNA was successfully extracted from the ethanol in which sponges were stored, although a reduced number of species were detected. Bivalve gill-tissue dissections, on the other hand, failed to reliably detect eDNA. Overall, our results show that vertebrate eDNA signals obtained from water samples and marine sponges are highly concordant. The strong similarity in eDNA signals demonstrates the potential of marine sponges as an additional tool for eDNA-based marine biodiversity surveys, by enabling the incorporation of larger sample numbers in eDNA surveys, reducing plastic waste, simplifying sample collection, and as a cost-efficient alternative. However, we note the importance to not detrimentally impact marine communities by, for example, nonlethal subsampling, specimen cloning, or using bycatch specimens.

Integrating technologies provides insight into the subsurface foraging behaviour of a humpback whale (Megaptera novaeangliae) feeding on walleye pollock (Gadus chalcogrammus) in Juan de Fuca Strait, Canada

Subsurface foraging is an important proportion of the activity budget of rorqual whales, yet information on their behaviour underwater remains challenging to obtain. Rorquals are assumed to feed throughout the water column and to select prey as a function of depth, availability and density, but there remain limitations in the precise identification of targeted prey. Current data on rorqual foraging in western Canadian waters have thus been limited to observations of prey species amenable to surface feeding, such as euphausiids and Pacific herring (Clupea pallasii), with no information on deeper alternative prey sources. We measured the foraging behaviour of a humpback whale (Megaptera novaeangliae) in Juan de Fuca Strait, British Columbia, using three complimentary methods: whale-borne tag data, acoustic prey mapping, and fecal sub-sampling. Acoustically detected prey layers were near the seafloor and consistent with dense schools of walleye pollock (Gadus chalcogrammus) distributed above more diffuse aggregations of pollock. Analysis of a fecal sample from the tagged whale confirmed that it had been feeding on pollock. Integrating the dive profile with the prey data revealed that the whale's foraging effort followed the general pattern of areal prey density, wherein the whale had a higher lunge-feeding rate at the highest prey abundance and stopped feeding when prey became limited. Our findings of a humpback whale feeding on seasonally energy-dense fish like walleye pollock, which are potentially abundant in British Columbia, suggests that pollock may be an important prey source for this rapidly growing whale population. This result is informative when assessing regional fishing activities for semi-pelagic species as well as the whales' vulnerability to fishing gear entanglements and feeding disturbances during a narrow window of prey acquisition.

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.

High-Throughput DNA Metabarcoding as an Approach for Ichthyoplankton Survey in Oujiang River Estuary, China

High-throughput DNA metabarcoding of mitochondrial 12S rRNA and Cyt b gene sequences was coupled with a morphology-based identification tool to assess ichthyoplankton community structure in Oujiang River Estuary, China. The performances of 12S and Cyt b barcoding markers were compared in terms of taxonomic resolution, detection and coverage, and their suitability was established for use as a quick and powerful ichthyoplankton assessment tool. A total of 30,138 ichthyoplankton (2462 eggs and 27,676 larvae) samples were collected from April to August 2015 and identified to 145 taxa belonging to 57 families and 105 genera. June and July were the main spawning months. Ichthyoplankton were more abundant around Lingkun and Qidu Islands and the upper parts of Oujiang River Estuary. The 12S gene marker presented higher species coverage and detection rate than Cyt b. DNA metabarcoding exhibited more representative species identification power than morphology. The findings reported in this study provided a key attempt towards the development of time-efficient and cost-effective ichthyoplankton identification and assessment tool.

Quantitative environmental DNA metabarcoding shows high potential as a novel approach to quantitatively assess fish community

The simultaneous conservation of species richness and evenness is important to effectively reduce biodiversity loss and keep ecosystem health. Environmental DNA (eDNA) metabarcoding has been used as a powerful tool for identifying community composition, but it does not necessarily provide quantitative information due to several methodological limitations. Thus, the quantification of eDNA through metabarcoding is an important frontier of eDNA-based biomonitoring. Particularly, the qMiSeq approach has recently been developed as a quantitative metabarcoding method and has attracted much attention due to its usefulness. The aim here was to evaluate the performance of the qMiSeq approach as a quantitative monitoring tool for fish communities by comparing the quantified eDNA concentrations with the results of fish capture surveys. The eDNA water sampling and the capture surveys using the electrical shocker were conducted at a total of 21 sites in four rivers in Japan. As a result, we found significant positive relationships between the eDNA concentrations of each species quantified by qMiSeq and both the abundance and biomass of each captured taxon at each site. Furthermore, for seven out of eleven taxa, a significant positive relationship was observed between quantified DNA concentrations by sample and the abundance and/or biomass. In total, our results demonstrated that eDNA metabarcoding with the qMiSeq approach is a suitable and useful tool for quantitative monitoring of fish communities. Due to the simplicity of the eDNA analysis, the eDNA metabarcoding with qMiSeq approach would promote further growth of quantitative monitoring of biodiversity.

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.

Fine-scale diversity of prey detected in humpback whale feces

Predator diets are largely influenced by prey availability and abundance. Yet, in heterogenous marine environments, identifying the prey species consumed by diving mammals remains a fundamental challenge. For rorqual whales, the energetic constraints of prey engulfment require that the whales seek areas of high prey abundance and execute discrete lunge feeding events on patches of high-density prey. Prey occurrences in feces should therefore provide meaningful insight into the dominant taxa in food patches selected by the animal. We investigated the prey consumed by humpback whales in three regions in southern British Columbia (BC), Canada, using opportunistic fecal sampling, microscopy, and DNA metabarcoding of 14 fecal samples. Fish including Pacific herring (Clupea pallasii), hake (Merluccius productus), and eulachon (Thaleichthys pacificus) were the most common fish species potentially targeted by humpback whales in two regions. The krill Euphausia pacifica was the most prevalent invertebrate DNA detected in all three regions, while sergestid and mysid shrimp may also be important. High DNA read abundances from walleye pollock (Gadus chalcogrammus) and sablefish (Anoplopoma fimbria) were also recovered in one sample each, suggesting that juveniles of these semi-pelagic species may occasionally be targeted. In general, we observed heavily digested fecal material that drove substantial dissimilarities in taxonomic resolution between polymerase chain reaction-based and morphological analyses of the feces. Pacific herring and walleye pollock were the only prey species confirmed by both methods. Our results highlight that molecular and visual analyses of fecal samples provide a complementary approach to diet analysis, with each method providing unique insight into prey diversity.

Dynamic analysis of the microbial communities and metabolome of healthy banana rhizosphere soil during one growth cycle

Background

The banana-growing rhizosphere soil ecosystem is very complex and consists of an entangled network of interactions between banana plants, microbes and soil, so identifying key components in banana production is difficult. Most of the previous studies on these interactions ignore the role of the banana plant. At present, there is no research on the the micro-ecological environment of the banana planting growth cycle.

Methods

Based on high-throughput sequencing technology and metabolomics technology, this study analyzed the rhizosphere soil microbial community and metabolic dynamics of healthy banana plants during one growth cycle.

Results

Assessing the microbial community composition of healthy banana rhizosphere soil, we found that the bacteria with the highest levels were Proteobacteria, Chloroflexi, and Acidobacteria, and the dominant fungi were Ascomycota, Basidiomycota, and Mortierellomycota. The metabolite profile of healthy banana rhizosphere soil showed that sugars, lipids and organic acids were the most abundant, accounting for about 50% of the total metabolites. The correlation network between fungi and metabolites was more complex than that of bacteria and metabolites. In a soil environment with acidic pH, bacterial genera showed a significant negative correlation with pH value, while fungal genera showed no significant negative correlation with pH value. The network interactions between bacteria, between fungi, and between bacteria and fungi were all positively correlated.

Conclusions

Healthy banana rhizosphere soil not only has a stable micro-ecology, but also has stable metabolic characteristics. The microorganisms in healthy banana rhizosphere soil have mutually beneficial rather than competitive relationships.

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.

Environmental nucleic acids: A field-based comparison for monitoring freshwater habitats using eDNA and eRNA

Nucleic acids released by organisms and isolated from environmental substrates are increasingly being used for molecular biomonitoring. While environmental DNA (eDNA) has received much attention, the potential of environmental RNA as a biomonitoring tool remains under-explored. Several recent studies using paired DNA and RNA metabarcoding of bulk samples suggest that RNA might better reflect "metabolically active" parts of the community. However, such studies mainly capture organismal eDNA and eRNA. For larger eukaryotes, isolation of extra-organismal RNA will be important, but viability needs to be examined in a field-based setting. In this study we evaluate (a) whether extra-organismal eRNA release from macroeukaryotes can be detected given its supposedly rapid degradation, and (b) if the same field collection methods for eDNA can be applied to eRNA. We collected eDNA and eRNA from water in lakes where fish community composition is well documented, enabling a comparison between the two nucleic acids in two different seasons with monitoring using conventional methods. We found that eRNA is released from macroeukaryotes and can be filtered from water and metabarcoded in a similar manner as eDNA to reliably provide species composition information. eRNA had a small but significantly greater true positive rate than eDNA, indicating that it correctly detects more species known to exist in the lakes. Given relatively small differences between the two molecules in describing fish community composition, we conclude that if eRNA provides significant advantages in terms of lability, it is a strong candidate to add to the suite of molecular monitoring tools.

Environmental DNA metabarcoding is a promising method for assaying fish diversity in cenotes of the Yucatán Peninsula, Mexico

The karst aquifer of the Yucatán Peninsula (YP) in southeastern Mexico is a unique ecosystem in which water-filled sinkholes, locally known as cenotes, connect subterranean waters with the surface. This system is home to around 20 species of freshwater fishes, including several that are endemic and/or threatened. Studies on this unique ichthyofauna have been partially hampered by the technical difficulties associated with sampling these habitats, particularly submerged caves. In this proof-of-concept study, we use environmental DNA (eDNA) metabarcoding to survey the diversity of freshwater fishes associated with the YP karst aquifer by sampling six cenotes from across the Ring of Cenotes region in northwestern Yucatán, a 180-km-diameter semicircular band of abundant karst sinkholes. Through a combination of conventional sampling (direct observation, fishing) and eDNA metabarcoding, we detected eight species of freshwater fishes across the six sampled cenotes. Overall, our eDNA metabarcoding approach was effective at detecting the presence of fishes from cenote water samples, including one of the two endemic cave-dwelling fish species restricted to the subterranean section of the aquifer. Although our study was focused on detecting fishes via eDNA, we also recovered DNA from several other vertebrate groups, particularly bats. These results suggest that the eDNA metabarcoding approach represents a promising and largely noninvasive method to assay aquatic biodiversity in these vulnerable habitats, allowing more effective, frequent, and wide-ranging surveys. Our detection of DNA from aerial and terrestrial vertebrate fauna implies that eDNA from cenotes, besides being a means to survey aquatic fauna, may also offer an effective way to quickly survey non-aquatic biodiversity associated with these persistent water bodies.

Saint Peter and Saint Paul Archipelago barcoded: Fish diversity in the remoteness and DNA barcodes reference library for metabarcoding monitoring

In order to monitor the effects of anthropogenic pressures in ecosystems, molecular techniques can be used to characterize species composition. Among molecular markers capable of identifying species, the cytochrome c oxidase I (COI) is the most used. However, new possibilities of biodiversity profiling have become possible, in which molecular fragments of medium and short-length can now be analyzed in metabarcoding studies. Here, a survey of fishes from the Saint Peter and Saint Paul Archipelago was barcoded using the COI marker, which allowed the identification of 21 species. This paved the way to further investigate the fish biodiversity of the archipelago, transitioning from barcoding to metabarcoding analysis. As preparatory steps for future metabarcoding studies, the first extensive COI library of fishes listed for these islands was constructed and includes new data generated in this survey as well as previously available data, resulting in a final database with 9,183 sequences from 169 species and 63 families of fish. A new primer specifically designed for those fishes was tested in silico to amplify a region of 262 bp. The new approach should guarantee a reliable surveillance of the archipelago and can be used to generate policies that will enhance the archipelago’s protection.

Assessing a megadiverse but poorly known community of fishes in a tropical mangrove estuary through environmental DNA (eDNA) metabarcoding

Biodiversity surveys are crucial for monitoring the status of threatened aquatic ecosystems, such as tropical estuaries and mangroves. Conventional monitoring methods are intrusive, time-consuming, substantially expensive, and often provide only rough estimates in complex habitats. An advanced monitoring approach, environmental DNA (eDNA) metabarcoding, is promising, although only few applications in tropical mangrove estuaries have been reported. In this study, we explore the advantages and limitations of an eDNA metabarcoding survey on the fish community of the Merbok Estuary (Peninsular Malaysia). COI and 12S eDNA metabarcoding assays collectively detected 178 species from 127 genera, 68 families, and 25 orders. Using this approach, significantly more species have been detected in the Merbok Estuary over the past decade (2010–2019) than in conventional surveys, including several species of conservation importance. However, we highlight three limitations: (1) in the absence of a comprehensive reference database the identities of several species are unresolved; (2) some of the previously documented specimen-based diversity was not captured by the current method, perhaps as a consequence of PCR primer specificity, and (3) the detection of non-resident species—stenohaline freshwater taxa (e.g., cyprinids, channids, osphronemids) and marine coral reef taxa (e.g., holocentrids, some syngnathids and sharks), not known to frequent estuaries, leading to the supposition that their DNA have drifted into the estuary through water movements. The community analysis revealed that fish diversity along the Merbok Estuary is not homogenous, with the upstream more diverse than further downstream. This could be due to the different landscapes or degree of anthropogenic influences along the estuary. In summary, we demonstrated the practicality of eDNA metabarcoding in assessing fish community and structure within a complex and rich tropical environment within a short sampling period. However, some limitations need to be considered and addressed to fully exploit the efficacy of this approach.

The relationship between eDNA density distribution and current fields around an artificial reef in the waters of Tateyama Bay, Japan

Monitoring of artificial reefs (ARs) has been conducted through such methods as visual censuses, surveys using fishing gear, and echo sounder. These methods have disadvantages: visual census is not possible at ARs in deeper waters, fishing gear surveys are invasive to fish individuals, and echo sounders have difficulty in species identification. A new AR monitoring method is required to compensate for these disadvantages. While eDNA has become a valid monitoring tool for marine biodiversities, it is influenced by degradation and transport of the molecules that affect information about the spatio-temporal distribution of fish. An understanding of the relationship between current fields and eDNA distribution, particularly in open waters, is critical when using eDNA as an index for fish aggregation at ARs. We investigated the relationship between eDNA distribution and current fields around an AR for four dominant species (Engraulis japonicus, Parapristipoma trilineatum, Scomber spp and Trachurus japonicus) in Tateyama Bay, Japan. The highest density of fish schools is formed directly above or at the upstream side of ARs. If we assume that the center of eDNA originates at these locations at an AR and eDNA is simply transported by currents, a higher density of eDNA would distribute downstream from the AR. However, our results indicate that eDNA distribution is in accord with actual fish distribution, namely eDNA densities are more abundant in the upstream side of ARs. We thus consider that eDNA distribution is more influenced by actual distribution patterns than by the transport processes.

Coevolution of Rumen Epithelial circRNAs with Their Microbiota and Metabolites in Response to Cold-Season Nutritional Stress in Tibetan Sheep

This study explores the effects of the coevolution of the host genome (the first genome) and gut microbiome (the second genome) on nutrition stress in Tibetan sheep during the cold season. The rumen epithelial tissue of six Tibetan sheep (Oula-type) was collected as experimental samples during the cold and warm seasons and the study lasted for half a year. The cDNA library was constructed and subjected to high-throughput sequencing. The circRNAs with significant differential expression were identified through bioinformatics analysis and functional prediction, and verified by real-time quantitative PCR (qRT-PCR). The results showed that a total of 56 differentially expressed (DE) circRNAs of rumen epithelial tissue were identified using RNA-seq technology, among which 29 were significantly upregulated in the cold season. The circRNA-miRNA regulatory network showed that DE circRNAs promoted the adaptation of Tibetan sheep in the cold season by targeting miR-150 and oar-miR-370-3p. The results of correlation analysis among circRNAs, microbiota, and metabolites showed that the circRNA NC_040275.1:28680890|28683112 had a very significant positive correlation with acetate, propionate, butyrate, and total volatile fatty acid (VFA) (p < 0.01), and had a significant positive correlation with Ruminococcus-1 (p < 0.05). In addition, circRNA NC_040256.1:78451819|78454934 and metabolites were enriched in the same KEGG pathway biosynthesis of amino acids (ko01230). In conclusion, the host genome and rumen microbiome of Tibetan sheep co-encoded a certain glycoside hydrolase (β-glucosidase) and coevolved efficient VFA transport functions and amino acid anabolic processes; thus, helping Tibetan sheep adapt to nutrient stress in the cold season in high-altitude areas.

Detection of fish species from marine protected areas of the North Sea using environmental DNA

This report describe the first application of environmental DNA-metabarcoding approach for the assessment of fish species diversity in two marine protected areas of the North Sea: the Doggerbank and the Sylt Outer Reef. We collected 64 water samples and detected 24 fish species. We discuss qualitative differences between MPAs and compare the results with those obtained from bottom-trawl surveys in the same areas. We found three additional species to those documented in the same year with trawls, including the critically endangered European eel.

Skimming for barcodes: rapid production of mitochondrial genome and nuclear ribosomal repeat reference markers through shallow shotgun sequencing

DNA barcoding is critical to conservation and biodiversity research, yet public reference databases are incomplete. Existing barcode databases are biased toward cytochrome oxidase subunit I (COI) and frequently lack associated voucher specimens or geospatial metadata, which can hinder reliable species assignments. The emergence of metabarcoding approaches such as environmental DNA (eDNA) has necessitated multiple marker techniques combined with barcode reference databases backed by voucher specimens. Reference barcodes have traditionally been generated by Sanger sequencing, however sequencing multiple markers is costly for large numbers of specimens, requires multiple separate PCR reactions, and limits resulting sequences to targeted regions. High-throughput sequencing techniques such as genome skimming enable assembly of complete mitogenomes, which contain the most commonly used barcoding loci (e.g., COI, 12S, 16S), as well as nuclear ribosomal repeat regions (e.g., ITS1&2, 18S). We evaluated the feasibility of genome skimming to generate barcode references databases for marine fishes by assembling complete mitogenomes and nuclear ribosomal repeats. We tested genome skimming across a taxonomically diverse selection of 12 marine fish species from the collections of the National Museum of Natural History, Smithsonian Institution. We generated two sequencing libraries per species to test the impact of shearing method (enzymatic or mechanical), extraction method (kit-based or automated), and input DNA concentration. We produced complete mitogenomes for all non-chondrichthyans (11/12 species) and assembled nuclear ribosomal repeats (18S-ITS1-5.8S-ITS2-28S) for all taxa. The quality and completeness of mitogenome assemblies was not impacted by shearing method, extraction method or input DNA concentration. Our results reaffirm that genome skimming is an efficient and (at scale) cost-effective method to generate all mitochondrial and common nuclear DNA barcoding loci for multiple species simultaneously, which has great potential to scale for future projects and facilitate completing barcode reference databases for marine fishes.

Analysis of the Persistence and Particle Size Distributional Shift of Sperm-Derived Environmental DNA to Monitor Jack Mackerel Spawning Activity

Environmental DNA (eDNA) analysis holds great promise as an efficient and noninvasive method to monitor not only the distribution of organisms but also their spawning activity. In eDNA analysis-based monitoring of spawning activity, the detection of sperm-derived eDNA is a key point; however, its characteristics and dynamics are completely unknown. The present study focuses on the persistence and particle size distribution (PSD) of eDNA derived from the sperm of Japanese jack mackerel. First, we investigated the time-dependent degradation and the PSD of sperm-derived eDNA by artificially adding sperm to seawater. Next, we kept fish in tanks and examined the changes in eDNA concentration and PSD before and after spawning. The results of two experiments showed that the degradation of sperm-derived eDNA proceeded rapidly, with PSD shifting to a smaller size regardless of the DNA region (Cyt b or ITS1). Additionally, it was shown that the nuclei and mitochondria released from sperm through degradation had a size distribution that was not simply dependent on each organelle size. These results will contribute to elucidating the characteristics and dynamics of eDNA specifically during the spawning season and to further developing eDNA analysis as a powerful tool for the monitoring of spawning activity.

Environmental DNA analysis as an emerging non-destructive method for plant biodiversity monitoring: a review

Environmental DNA (eDNA) analysis has recently transformed and modernized biodiversity monitoring. The accurate detection, and to some extent quantification, of organisms (individuals/populations/communities) in environmental samples is galvanizing eDNA as a successful cost and time-efficient biomonitoring technique. Currently, eDNA's application to plants remains more limited in implementation and scope compared to animals and microorganisms. This review evaluates the development of eDNA-based methods for (vascular) plants, comparing its performance and power of detection with that of traditional methods, to critically evaluate and advise best-practices needed to innovate plant biomonitoring. Recent advancements, standardization and field applications of eDNA-based methods have provided enough scope to utilize it in conservation biology for numerous organisms. Despite our review demonstrating only 13% of all eDNA studies focus on plant taxa to date, eDNA has considerable environmental DNA has considerable potential for plants, where successful detection of invasive, endangered and rare species, and community-level interpretations have provided proof-of-concept. Monitoring methods using eDNA were found to be equal or more effective than traditional methods; however, species detection increased when both methods were coupled. Additionally, eDNA methods were found to be effective in studying species interactions, community dynamics and even effects of anthropogenic pressure. Currently, elimination of potential obstacles (e.g. lack of relevant DNA reference libraries for plants) and the development of user-friendly protocols would greatly contribute to comprehensive eDNA-based plant monitoring programs. This is particularly needed in the data-depauperate tropics and for some plant groups (e.g., Bryophytes and Pteridophytes). We further advocate to coupling traditional methods with eDNA approaches, as the former is often cheaper and methodologically more straightforward, while the latter offers non-destructive approaches with increased discrimination ability. Furthermore, to make a global platform for eDNA, governmental and academic-industrial collaborations are essential to make eDNA surveys a broadly adopted and implemented, rapid, cost-effective and non-invasive plant monitoring approach.

Tourmaline: A containerized workflow for rapid and iterable amplicon sequence analysis using QIIME 2 and Snakemake

Background

Amplicon sequencing (metabarcoding) is a common method to survey diversity of environmental communities whereby a single genetic locus is amplified and sequenced from the DNA of whole or partial organisms, organismal traces (e.g., skin, mucus, feces), or microbes in an environmental sample. Several software packages exist for analyzing amplicon data, among which QIIME 2 has emerged as a popular option because of its broad functionality, plugin architecture, provenance tracking, and interactive visualizations. However, each new analysis requires the user to keep track of input and output file names, parameters, and commands; this lack of automation and standardization is inefficient and creates barriers to meta-analysis and sharing of results.

Findings

We developed Tourmaline, a Python-based workflow that implements QIIME 2 and is built using the Snakemake workflow management system. Starting from a configuration file that defines parameters and input files—a reference database, a sample metadata file, and a manifest or archive of FASTQ sequences—it uses QIIME 2 to run either the DADA2 or Deblur denoising algorithm; assigns taxonomy to the resulting representative sequences; performs analyses of taxonomic, alpha, and beta diversity; and generates an HTML report summarizing and linking to the output files. Features include support for multiple cores, automatic determination of trimming parameters using quality scores, representative sequence filtering (taxonomy, length, abundance, prevalence, or ID), support for multiple taxonomic classification and sequence alignment methods, outlier detection, and automated initialization of a new analysis using previous settings. The workflow runs natively on Linux and macOS or via a Docker container. We ran Tourmaline on a 16S ribosomal RNA amplicon data set from Lake Erie surface water, showing its utility for parameter optimization and the ability to easily view interactive visualizations through the HTML report, QIIME 2 viewer, and R- and Python-based Jupyter notebooks.

Conclusion

Automated workflows like Tourmaline enable rapid analysis of environmental amplicon data, decreasing the time from data generation to actionable results. Tourmaline is available for download at github.com/aomlomics/tourmaline.