Replication levels, false presences and the estimation of the presence/absence from eDNA metabarcoding data

The development of terrestrial ecosystems emerging after glacier retreat

The global retreat of glaciers is dramatically altering mountain and high-latitude landscapes, with new ecosystems developing from apparently barren substrates1-4. The study of these emerging ecosystems is critical to understanding how climate change interacts with microhabitat and biotic communities and determines the future of ice-free terrains1,5. Here, using a comprehensive characterization of ecosystems (soil properties, microclimate, productivity and biodiversity by environmental DNA metabarcoding6) across 46 proglacial landscapes worldwide, we found that all the environmental properties change with time since glaciers retreated, and that temperature modulates the accumulation of soil nutrients. The richness of bacteria, fungi, plants and animals increases with time since deglaciation, but their temporal patterns differ. Microorganisms colonized most rapidly in the first decades after glacier retreat, whereas most macroorganisms took longer. Increased habitat suitability, growing complexity of biotic interactions and temporal colonization all contribute to the increase in biodiversity over time. These processes also modify community composition for all the groups of organisms. Plant communities show positive links with all other biodiversity components and have a key role in ecosystem development. These unifying patterns provide new insights into the early dynamics of deglaciated terrains and highlight the need for integrated surveillance of their multiple environmental properties5.

Detection, movement and persistence of invertebrate eDNA in groundwater

Sampling groundwater biodiversity is difficult because of limited access and issues with species identification. Environmental DNA (eDNA) provides a viable alternative to traditional sampling approaches, however limited knowledge of the abundance and fate of DNA in groundwater hinders the interpretation of data from these environments. Groundwater environments are dark and have lower oxygen concentrations and microbial activity than surface waters. Consequently, assumptions about DNA fate in surface ecosystems may not apply to groundwaters. Here, we test the longevity and transport of eDNA in groundwater within a static microcosm and a flow-through mesocosm. A variety of invertebrates were placed within a mesocosm and microcosm to enable DNA shedding, and then removed. DNA persisted for up to 5 weeks after their removal in the static experiment and was detected between 9 and 33 days in the flow-through experiment. Sediments and water both proved important for eDNA detection. Crustacean DNA was detected sporadically and unpredictably, whereas non-crustacean DNA was detected more frequently despite their lower densities. We suggest that detecting crustaceans poses a challenge to utilising eDNA approaches for stygofauna monitoring. This is confounded by the scarcity of sequences for stygofauna in reference databases. Further research is needed before eDNA alone can be routinely employed for stygofauna detection.

Using ancient sedimentary DNA to forecast ecosystem trajectories under climate change

Ecosystem response to climate change is complex. In order to forecast ecosystem dynamics, we need high-quality data on changes in past species abundance that can inform process-based models. Sedimentary ancient DNA (sedaDNA) has revolutionised our ability to document past ecosystems' dynamics. It provides time series of increased taxonomic resolution compared to microfossils (pollen, spores), and can often give species-level information, especially for past vascular plant and mammal abundances. Time series are much richer in information than contemporary spatial distribution information, which have been traditionally used to train models for predicting biodiversity and ecosystem responses to climate change. Here, we outline the potential contribution of sedaDNA to forecast ecosystem changes. We showcase how species-level time series may allow quantification of the effect of biotic interactions in ecosystem dynamics, and be used to estimate dispersal rates when a dense network of sites is available. By combining palaeo-time series, process-based models, and inverse modelling, we can recover the biotic and abiotic processes underlying ecosystem dynamics, which are traditionally very challenging to characterise. Dynamic models informed by sedaDNA can further be used to extrapolate beyond current dynamics and provide robust forecasts of ecosystem responses to future climate change. This article is part of the theme issue 'Ecological novelty and planetary stewardship: biodiversity dynamics in a transforming biosphere'.

eDNA Metabarcoding Reveals the Species-Area Relationship of Amphibians on the Zhoushan Archipelago

The species-area relationship is important for understanding species diversity patterns at spatial scales, but few studies have examined the relationship using environmental DNA (eDNA) techniques. We investigated amphibian diversity on 21 islands of the Zhoushan Archipelago and nearby mainland areas in China using the combination of eDNA metabarcoding and the traditional line transect method (TLTM) and identified the species-area relationship for amphibians on the islands. The mean detection probability of eDNA is 0.54, while the mean detection probability of TLTM is 0.24. The eDNA metabarcoding detected eight amphibian species on the islands and nine species in the mainland areas, compared with seven species on the islands and nine species in the mainland areas that were identified by TLTM. Amphibian richness on the islands increased with island area and habitat diversity. The species-area relationship for amphibians in the archipelago was formulated as the power function (S = 0.47A0.21) or exponential function (S = 2.59 + 2.41 (logA)). Our results suggested that eDNA metabarcoding is more sensitive for the detection of amphibian species. The combined use of eDNA metabarcoding and the traditional line transect method may optimize the survey results for amphibians.

Dynamics and drivers of mycorrhizal fungi after glacier retreat

The development of terrestrial ecosystems depends greatly on plant mutualists such as mycorrhizal fungi. The global retreat of glaciers exposes nutrient-poor substrates in extreme environments and provides a unique opportunity to study early successions of mycorrhizal fungi by assessing their dynamics and drivers. We combined environmental DNA metabarcoding and measurements of local conditions to assess the succession of mycorrhizal communities during soil development in 46 glacier forelands around the globe, testing whether dynamics and drivers differ between mycorrhizal types. Mycorrhizal fungi colonized deglaciated areas very quickly (< 10 yr), with arbuscular mycorrhizal fungi tending to become more diverse through time compared to ectomycorrhizal fungi. Both alpha- and beta-diversity of arbuscular mycorrhizal fungi were significantly related to time since glacier retreat and plant communities, while microclimate and primary productivity were more important for ectomycorrhizal fungi. The richness and composition of mycorrhizal communities were also significantly explained by soil chemistry, highlighting the importance of microhabitat for community dynamics. The acceleration of ice melt and the modifications of microclimate forecasted by climate change scenarios are expected to impact the diversity of mycorrhizal partners. These changes could alter the interactions underlying biotic colonization and belowground-aboveground linkages, with multifaceted impacts on soil development and associated ecological processes.

Cost-effort analysis of Baited Remote Underwater Video (BRUV) and environmental DNA (eDNA) in monitoring marine ecological communities

Monitoring the diversity and distribution of species in an ecosystem is essential to assess the success of restoration strategies. Implementing biomonitoring methods, which provide a comprehensive assessment of species diversity and mitigate biases in data collection, holds significant importance in biodiversity research. Additionally, ensuring that these methods are cost-efficient and require minimal effort is crucial for effective environmental monitoring. In this study we compare the efficiency of species detection, the cost and the effort of two non-destructive sampling techniques: Baited Remote Underwater Video (BRUV) and environmental DNA (eDNA) metabarcoding to survey marine vertebrate species. Comparisons were conducted along the Sussex coast upon the introduction of the Nearshore Trawling Byelaw. This Byelaw aims to boost the recovery of the dense kelp beds and the associated biodiversity that existed in the 1980s. We show that overall BRUV surveys are more affordable than eDNA, however, eDNA detects almost three times as many species as BRUV. eDNA and BRUV surveys are comparable in terms of effort required for each method, unless eDNA analysis is carried out externally, in which case eDNA requires less effort for the lead researchers. Furthermore, we show that increased eDNA replication yields more informative results on community structure. We found that using both methods in conjunction provides a more complete view of biodiversity, with BRUV data supplementing eDNA monitoring by recording species missed by eDNA and by providing additional environmental and life history metrics. The results from this study will serve as a baseline of the marine vertebrate community in Sussex Bay allowing future biodiversity monitoring research projects to understand community structure as the ecosystem recovers following the removal of trawling fishing pressure. Although this study was regional, the findings presented herein have relevance to marine biodiversity and conservation monitoring programs around the globe.

In silico assessment of 18S rDNA metabarcoding markers for the characterization of nematode communities

Nematodes are keystone actors of soil, freshwater and marine ecosystems, but the complexity of morphological identification has limited broad-scale monitoring of nematode biodiversity. DNA metabarcoding is increasingly used to assess nematode diversity but requires universal primers with high taxonomic coverage and high taxonomic resolution. Several primers have been proposed for the metabarcoding of nematode diversity, many of which target the 18S rRNA gene. In silico analyses have a great potential to assess key parameters of primers, including taxonomic coverage, resolution and specificity. Based on a recently-available reference database, we tested in silico the performance of fourteen commonly used and one newly optimized primer for nematode metabarcoding. Most primers showed very good coverage, amplifying most of the sequences in the reference database, while four markers showed limited coverage. All primers showed good taxonomic resolution. Resolution was particularly good if the aim was the identification of higher-level taxa, such as genera or families. Overall, species-level resolution was higher for primers amplifying long fragments. None of the primers was highly specific for nematodes as, despite some variation, they all amplified a large number of other eukaryotes. Differences in performance across primers highlight the complexity of the choice of markers appropriate for the metabarcoding of nematodes, which depends on a trade-off between taxonomic resolution and the length of amplified fragments. Our in silico analyses provide new insights for the identification of the most appropriate primers, depending on the study goals and the origin of DNA samples. This represents an essential step to design and optimize metabarcoding studies assessing nematode diversity.

Comparative evaluation for the performance of environmental DNA and RNA analyses targeting mitochondrial and nuclear genes from ayu (Plecoglossus altivelis)

Environmental DNA and RNA (eDNA and eRNA; collectively eNA) analyses have the potential for non-invasive and cost-efficient biomonitoring compared with traditional capture-based surveys. Although various types of eNA particles, including not only mitochondrial eDNA but also nuclear eDNA and their transcripts, are present in the water, performances of eNA detection and quantification have not yet been evaluated sufficiently across multiple mitochondrial and nuclear genes. We conducted a tank experiment with ayu (Plecoglossus altivelis) to compare the detection sensitivity, yields per water sample, and quantification variability between replicates of each type of eNAs. The assay targeting the multi-copy nuclear gene exhibited a higher sensitivity than the assay targeting the mitochondrial gene, and both the target eDNA and eRNA concentrations per water sample were higher for the nuclear gene. On the contrary, variation in eRNA quantifications per sample does not necessarily correspond to that in eDNA, and the intra-sample quantification variability (represented as the CVs between PCR replicates) tended to be larger for eRNA than eDNA. Our results suggested that, even if suitable to the sensitive detection of species occurrence, the use of eRNA particularly derived from multi-copy nuclear gene may not be necessarily appropriate for the reliable assessment of species abundance. The findings in this study would help optimize eNA analyses for making biomonitoring and stock assessment in aquatic environments more efficient and reliable.

Interplay of competition and facilitation in grazing succession by migrant Serengeti herbivores

Competition, facilitation, and predation offer alternative explanations for successional patterns of migratory herbivores. However, these interactions are difficult to measure, leaving uncertainty about the mechanisms underlying body-size-dependent grazing-and even whether succession occurs at all. We used data from an 8-year camera-trap survey, GPS-collared herbivores, and fecal DNA metabarcoding to analyze the timing, arrival order, and interactions among migratory grazers in Serengeti National Park. Temporal grazing succession is characterized by a "push-pull" dynamic: Competitive grazing nudges zebra ahead of co-migrating wildebeest, whereas grass consumption by these large-bodied migrants attracts trailing, small-bodied gazelle that benefit from facilitation. "Natural experiments" involving intense wildfires and rainfall respectively disrupted and strengthened these effects. Our results highlight a balance between facilitative and competitive forces in co-regulating large-scale ungulate migrations.

Using eDNA to survey amphibians: Methods, applications, and challenges

In recent years, environmental DNA (eDNA) has received attention from biologists due to its sensitivity, convenience, labor and material efficiency, and lack of damage to organisms. The extensive application of eDNA has opened avenues for the monitoring and biodiversity assessment of amphibians, which are frequently small and difficult to observe in the field, in areas such as biodiversity survey assessment and detection of specific, rare and threatened, or alien invasive species. However, the accuracy of eDNA can be influenced by factors such as ambient temperature, pH, and false positives or false negatives, which makes eDNA an adjunctive tool rather than a replacement for traditional surveys. This review provides a concise overview of the eDNA method and its workflow, summarizes the differences between applying eDNA for detecting amphibians and other organisms, reviews the research progress in eDNA technology for amphibian monitoring, identifies factors influencing detection efficiency, and discusses the challenges and prospects of eDNA. It aims to serve as a reference for future research on the application of eDNA in amphibian detection.

The importance of species addition 'versus' replacement varies over succession in plant communities after glacier retreat

The mechanisms underlying plant succession remain highly debated. Due to the local scope of most studies, we lack a global quantification of the relative importance of species addition 'versus' replacement. We assessed the role of these processes in the variation (β-diversity) of plant communities colonizing the forelands of 46 retreating glaciers worldwide, using both environmental DNA and traditional surveys. Our findings indicate that addition and replacement concur in determining community changes in deglaciated sites, but their relative importance varied over time. Taxa addition dominated immediately after glacier retreat, as expected in harsh environments, while replacement became more important for late-successional communities. These changes were aligned with total β-diversity changes, which were more pronounced between early-successional communities than between late-successional communities (>50 yr since glacier retreat). Despite the complexity of community assembly during plant succession, the observed global pattern suggests a generalized shift from the dominance of facilitation and/or stochastic processes in early-successional communities to a predominance of competition later on.

Environmental drivers of Arctic communities based on metabarcoding of marine sediment eDNA

Our ability to assess biodiversity at relevant spatial and temporal scales for informing management is of increasing importance given this is foundational to identify and mitigate the impacts of global change. Collecting baseline information and tracking ecological changes are particularly important for areas experiencing rapid changes and representing data gaps such as Arctic marine ecosystems. Environmental DNA has the potential to provide such data. We extracted environmental DNA from 90 surface sediment samples to assess eukaryote diversity around Greenland and Svalbard using two separate primer pairs amplifying different sections of the 18S rRNA gene. We detected 27 different phyla and 99 different orders and found that temperature and the change in temperature explained the most variation in the community in a single linear model, while latitude, sea ice cover and change in temperature explained the most variation in the community when assessed by individual non-linear models. We identified potential indicator taxa for Arctic climate change, including a terebellid annelid worm. In conclusion, our study demonstrates that environmental DNA offers a feasible method to assess biodiversity and identifies warming as a key driver of differences in biodiversity across these remote ecosystems.

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

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

First pelagic fish biodiversity assessment of Cosmonaut Sea based on environmental DNA

The Cosmonaut Sea is a typical marginal sea in East Antarctica that has not yet been greatly impacted by climate change. As one of the least explored areas in the Southern Ocean, our knowledge regarding its fish taxonomy and diversity has been sparse. eDNA metabarcoding, as an emerging and promising tool for marine biodiversity research and monitoring, has been widely used across taxa and habitats. During the 38th Chinese Antarctic Research Expedition (CHINARE-38), we collected seawater and surface sediment samples from 38 stations in the Cosmonaut Sea and performed the first, to our knowledge, eDNA analysis of fish biodiversity in the Southern Ocean based on the molecular markers of 12S rRNA and 16S rRNA. There were 48 fish species detected by the two markers in total, with 30 and 34 species detected by the 12S rRNA and 16S rRNA marker, respectively. This was more than the trawling results (19 species) and historical survey records (16 species, "BROKE-West" cruise). With some nonsignificant differences between the Gunnerus Ridge and the Oceanic Area of Enderby Land, the Cosmonaut Sea had a richer fish biodiversity in this research compared with previous studies, and its overall composition and distribution patterns were consistent with what we know in East Antarctica. We also found that the eDNA composition of fish in the Cosmonaut Sea might be related to some environmental factors. Our study demonstrated that the use of the eDNA technique for Antarctic fish biodiversity research is likely to yield more information with less sampling effort than traditional methods. In the context of climate change, the eDNA approach will provide a novel and powerful tool that is complementary to traditional methods for polar ecology research.

Environmental DNA metabarcoding of water samples as a tool for monitoring Iberian freshwater fish composition

Environmental DNA (eDNA) metabarcoding has been increasingly used to monitor the community assemblages of a wide variety of organisms. Here, we test the efficacy of eDNA metabarcoding to assess the composition of Iberian freshwater fishes, one of the most endangered groups of vertebrates in Spain. For this purpose, we sampled 12 sampling sites throughout one of Spain's largest basins, the Duero, which is home to approximately 70% of the genera and 30% of the primary freshwater fish in Spain. We sampled these sampling sites in the summer by using electrofishing, a traditional sampling method, and eDNA metabarcoding of river water samples using the mitochondrial 12S rRNA gene (12S) as a marker. We also resampled four of these sampling sites in autumn by eDNA. We compared the results obtained through eDNA metabarcoding with those of electrofishing surveys (ones conducted for the present study and past ones) and assessed the suitability of 12S as an eDNA metabarcoding marker for this group of freshwater fishes. We found that the 12S fragment, analysed for 25 Iberian species, showed sufficient taxonomic resolution to be useful for eDNA approaches, and even showed population-level differences in the studied populations across the tissue samples for Achondrostoma arcasii. In most cases, a greater number of species was detected through eDNA metabarcoding than through electrofishing. Based on our results, eDNA metabarcoding is a powerful tool to study the freshwater fish composition in the Iberian Peninsula and to unmask cryptic diversity. However, we highlight the need to generate a local genetic database for 12S gene for such studies and to interpret the results with caution when studying only mitochondrial DNA. Finally, our survey shows that the high detection sensitivity of eDNA metabarcoding and the non-invasiveness of this method allows it to act as a detection system for species of low abundance, such as early invasive species or species in population decline, two key aspects of conservation management of Spanish freshwater fishes.

eDNA Metabarcoding- and Microscopic Analysis for Diet Determination in Waterfowl, a Comparative Study in Vejlerne, Denmark

Understanding diets and structural food webs are keys to the apprehension of ecological communities, upon which conservation and management biology are based. The understanding of grazing and habitat choice for waterfowl is one of the most important topics for avian ecologists today and can, to some degree, be answered by dietary analysis. Droppings collected from four waterfowl, the Eurasian wigeon (Anas penelope), Greylag goose (Anser anser), pink-footed goose (Anser brachyrhynchus) and Barnacle goose (Branta leucopsis) in Vejlerne (Denmark), were analysed microscopically and through eDNA metabarcoding with the use of next generation sequencing (NGS) to accumulate knowledge about the diet of these waterfowl. In total, 120 dropping samples were microscopically analysed, of which the eDNA metabarcoding analysis was done on 79 samples. The prey items were identified according to the taxonomic level of species, and a qualitative method, frequency of occurrence (FO) and FO calculated as a percentage, was used in order to compare the results from the two methods. As neither of the methods was able to encompass all species discovered when combining the two methods, it was concluded in this study that the two methods can support each other in a dietary analysis of waterfowl, but not replace one another.

Environmental DNA metabarcoding serves as a promising method for aquatic species monitoring and management: A review focused on its workflow, applications, challenges and prospects

Marine and freshwater biodiversity is under threat from both natural and manmade causes. Biological monitoring is currently a top priority for biodiversity protection. Given present limitations, traditional biological monitoring methods may not achieve the proposed monitoring aims. Environmental DNA metabarcoding technology reflects species information by capturing and extracting DNA from environmental samples, using molecular biology techniques to sequence and analyze the DNA, and comparing the obtained information with existing reference libraries to obtain species identification. However, its practical application has highlighted several limitations. This paper summarizes the main steps in the environmental application of eDNA metabarcoding technology in aquatic ecosystems, including the discovery of unknown species, the detection of invasive species, and evaluations of biodiversity. At present, with the rapid development of big data and artificial intelligence, certain advanced technologies and devices can be combined with environmental DNA metabarcoding technology to promote further development of aquatic species monitoring and management.

Urban landscape-level biodiversity assessments of aquatic and terrestrial vertebrates by environmental DNA metabarcoding

Globally, expansive urbanization profoundly alters natural habitats and the associated biota. Monitoring biodiversity in cities can provide essential information for conservation management, but the complexity of urban landscapes poses serious challenges to conventional observational and capture-based surveys. Here we assessed pan-vertebrate biodiversity, including both aquatic and terrestrial taxa, using environmental DNA (eDNA) sampled from 109 water sites across Beijing, China. Using eDNA metabarcoding with a single primer set (Tele02), we detected 126 vertebrate species, including 73 fish, 39 birds, 11 mammals, and 3 reptiles belonging to 91 genera, 46 families, and 22 orders. The probability of detection from eDNA varied substantially among species and was related to their lifestyle, as shown by the greater detectability of fish compared to that of terrestrial and arboreal (birds and mammals) groups, as well as the greater detectability of water birds compared to that of forest birds (Wilcoxon rank-sum test p = 0.007). Furthermore, the eDNA detection probabilities across all vertebrates (Wilcoxon rank-sum test p = 0.009), as well as for birds (p < 0.001), were higher at lentic sites in comparison with lotic sites. Also, the detected biodiversity was positively correlated with lentic waterbody size for fish (Spearman p = 0.012), but not for other groups. Our results demonstrate the capacity of eDNA metabarcoding to efficiently surveil diverse vertebrate communities across an extensive spatial scale in heterogenous urban landscapes. With further methodological development and optimization, the eDNA approach has great potential for non-invasive, efficient, economic, and timely assessments of biodiversity responses to urbanization, thus guiding city ecosystem conservation management.

Authentication of milk thistle commercial products using UHPLC-QTOF-ESI + MS metabolomics and DNA metabarcoding

BACKGROUND

Milk thistle is one of the most popular hepatoprotectants, and is often sold in combination with other ingredients. Botanical supplements are known to be vulnerable to contamination and adulteration, and emerging technologies show promise to improve their quality control.

METHODS

Untargeted and semi-targeted metabolomics based on UHPLC-QTOF-ESI⁺MS techniques, UV spectrometry, and DNA metabarcoding using Illumina MiSeq were used to authenticate eighteen milk thistle botanical formulations (teas, capsules, tablets, emulsion).

RESULTS

Untargeted metabolomics separated 217 molecules and by multivariate analysis the discrimination between the different preparations was established. The semi-targeted metabolomics focused on 63 phytochemicals, mainly silymarin flavonolignans and flavonoids, that may be considered as putative biomarkers of authenticity. All formulations contained molecules from silymarin complexes at different levels. The quantitative evaluation of silybins was done using in parallel UV spectrometry and UHPLC-QTOF-ESI⁺MS and their correlations were compared. DNA metabarcoding detected milk thistle in eleven out of sixteen retained preparations, whereas two others had incomplete evidence of milk thistle despite metabolomics validating specific metabolites, e.g., silymarin complex, identified and quantified in all samples. Meanwhile, the DNA metabarcoding provided insights into the total species composition allowing the interpretation of the results in a broad context.

CONCLUSION

Our study emphasizes that combining spectroscopic, chromatographic, and genetic techniques bring complementary information to guarantee the quality of the botanical formulations.

Using environmental DNA metabarcoding to monitor fish communities in small rivers and large brooks: Insights on the spatial scale of information

Monitoring fish communities is central to the evaluation of ecological health of rivers. Both presence/absence of fish species and their relative quantity in local fish assemblages are crucial parameters to measure. Fish communities in lotic systems are traditionally monitored via electrofishing, characterized by a known limited efficiency and high survey costs. Analysis of environmental DNA could serve as a non-destructive alternative for detection and quantification of lotic fish communities, but this approach still requires further insights in practical sampling schemes incorporating transport and dilution of the eDNA particles; optimization of predictive power and quality assurance of the molecular detection method. Via a controlled cage experiment, we aim to extend the knowledge on streamreach of eDNA in small rivers and large brooks, as laid out in the European Water Framework Directive's water typology. Using a high and low source biomass in two river transects of a species-poor river characterized by contrasting river discharge rates, we found strong and significant correlations between the eDNA relative species abundances and the relative biomass per species in the cage community. Despite a decreasing correlation over distance, the underlying community composition remained stable from 25 to 300 m, or up to 1 km downstream of the eDNA source, depending on the river discharge rate. Such decrease in similarity between relative source biomass and the corresponding eDNA-based community profile with increasing distance downstream from the source, might be attributed to variation in species-specific eDNA persistence. Our findings offer crucial insights on eDNA behaviour and characterization of riverine fish communities. We conclude that water sampled from a relatively small river offers an adequate eDNA snapshot of the total fish community in the 300-1000 m upstream transect. The potential application for other river systems is further discussed.

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.

Horizon scan of DNA-based methods for quality control and monitoring of herbal preparations

Herbal medicines and preparations are widely used in healthcare systems globally, but concerns remain about their quality and safety. New herbal products are constantly being introduced to the market under varying regulatory frameworks, with no global consensus on their definition or characterization. These biologically active mixtures are sold through complex globalized value chains, which create concerns around contamination and profit-driven adulteration. Industry, academia, and regulatory bodies must collaborate to develop innovative strategies for the identification and authentication of botanicals and their preparations to ensure quality control. High-throughput sequencing (HTS) has significantly improved our understanding of the total species diversity within DNA mixtures. The standard concept of DNA barcoding has evolved over the last two decades to encompass genomic data more broadly. Recent research in DNA metabarcoding has focused on developing methods for quantifying herbal product ingredients, yielding meaningful results in a regulatory framework. Techniques, such as loop-mediated isothermal amplification (LAMP), DNA barcode-based Recombinase Polymerase Amplification (BAR-RPA), DNA barcoding coupled with High-Resolution Melting (Bar-HRM), and microfluidics-based methods, offer more affordable tests for the detection of target species. While target capture sequencing and genome skimming are considerably increasing the species identification resolution in challenging plant clades, ddPCR enables the quantification of DNA in samples and could be used to detect intended and unwanted ingredients in herbal medicines. Here, we explore the latest advances in emerging DNA-based technologies and the opportunities they provide as taxa detection tools for evaluating the safety and quality of dietary supplements and herbal medicines.

The Next Generation of Microbial Ecology and Its Importance in Environmental Sustainability

Collectively, we have been reviewers for microbial ecology, genetics and genomics studies that include environmental DNA (eDNA), microbiome studies, and whole bacterial genome biology for Microbial Ecology and other journals for about three decades. Here, we wish to point out trends and point to areas of study that readers, especially those moving into the next generation of microbial ecology research, might learn and consider. In this communication, we are not saying the work currently being accomplished in microbial ecology and restoration biology is inadequate. What we are saying is that a significant milestone in microbial ecology has been reached, and approaches that may have been overlooked or were unable to be completed before should be reconsidered in moving forward into a new more ecological era where restoration of the ecological trajectory of systems has become critical. It is our hope that this introduction, along with the papers that make up this special issue, will address the sense of immediacy and focus needed to move into the next generation of microbial ecology study.

DNA High-Throughput Sequencing for Arthropod Gut Content Analysis to Evaluate Effectiveness and Safety of Biological Control Agents

The search for effective biological control agents without harmful non-target effects has been constrained by the use of impractical (field direct observation) or imprecise (cage experiments) methods. While advances in the DNA sequencing methods, more specifically the development of high-throughput sequencing (HTS), have been quickly incorporated in biodiversity surveys, they have been slow to be adopted to determine arthropod prey range, predation rate and food web structure, and critical information to evaluate the effectiveness and safety of a biological control agent candidate. The lack of knowledge on how HTS methods could be applied by ecological entomologists constitutes part of the problem, although the lack of expertise and the high cost of the analysis also are important limiting factors. In this review, we describe how the latest HTS methods of metabarcoding and Lazaro, a method to identify prey by mapping unassembled shotgun reads, can serve biological control research, showing both their power and limitations. We explain how they work to determine prey range and also how their data can be used to estimate predation rates and subsequently be translated into food webs of natural enemy and prey populations helping to elucidate their role in the community. We present a brief history of prey detection through molecular gut content analysis and also the attempts to develop a more precise formula to estimate predation rates, a problem that still remains. We focused on arthropods in agricultural ecosystems, but most of what is covered here can be applied to natural systems and non-arthropod biological control candidates as well.

Seasonal variation of population and individual dietary niche in the avivorous bat, Ia io

The variation in niche breadth can affect how species respond to environmental and resource changes. However, there is still no clear understanding of how seasonal variability in food resources impacts the variation of individual dietary diversity, thereby affecting the dynamics of a population's dietary niche breadth. Optimal foraging theory (OFT) and the niche variation hypothesis (NVH) predict that when food resources are limited, the population niche breadth will widen or narrow due to increased within-individual dietary diversity and individual specialization or reduced within-individual dietary diversity, respectively. Here, we used DNA metabarcoding to examine the composition and seasonality of diets of the avivorous bat Ia io. Furthermore, we investigated how the dietary niches changed among seasons and how the population niche breadth changed when the availability of insect resources was reduced in autumn. We found that there was differentiation in dietary niches among seasons and a low degree of overlap, and the decrease of insect resource availability and the emergence of ecological opportunities of nocturnal migratory birds might drive dietary niche shifts toward birds in I. io. However, the population's dietary niche breadth did not broaden by increasing the within-individual dietary diversity or individual specialization, but rather became narrower by reducing dietary diversity via predation on bird resources that served as an ecological opportunity when insect resources were scarce in autumn. Our findings were consistent with the predictions of OFT, because birds as prey for bats provided extremely different resources from those of insects in size and nutritional value. Our work highlights the importance of size and quality of prey resources along with other factors (i.e., physiological, behavioral, and life-history traits) in dietary niche variation.

Environmental DNA Assay for the Detection of the American Bullfrog (Lithobates catesbeianus) in the Early Stages of the Invasion in the Ebre Delta

The American bullfrog (Lithobates catesbeianus) is considered to be one of the most harmful invasive species. In the Iberian Peninsula, this species had been cited occasionally until the year 2018, when L. catesbeianus appeared in the Ebre Delta, and, for the first time, it started breeding in a territory of the Peninsula. Using environmental DNA (eDNA) analysis and visual surveys, the American bullfrog invasion in the Ebre Delta was monitored across two consecutive years (2019-2020). No specimens were observed in 2019, and results for the eDNA survey also failed to detect this species in the Delta. In 2020, two individuals were captured and, under the most conservative criteria to constrain the number of positive detections, eDNA analyses detected the presence of the American bullfrog in at least five locations. Performing an eDNA assay yielded a higher sensitivity with a lower sampling effort than traditional methods. Although the American bullfrog does not appear to still be well-established in the Ebre Delta, only a few bullfrog individuals could be enough for their establishment in suitable habitats. In this context, eDNA assays are essential tools to facilitate the detection, control, and eradication of this species in the first stage of the invasion process.

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.

Maximising signal-to-noise ratios in environmental DNA-based monitoring

Man's impacts on global ecosystems are increasing and there is a growing demand that these activities be appropriately monitored. Monitoring requires measurement of a response metric ('signal') that changes maximally and consistently in response to the monitored activity irrespective of other factors ('noise'), thus maximising the signal-to-noise ratio. Indices derived from time-consuming morphology-based taxonomic identification of organisms are a core part of many monitoring programmes. Metabarcoding is an alternative to morphology-based identification and involves the sequencing of short fragments of DNA ('markers') from multiple taxa simultaneously. DNA suitable for metabarcoding includes that extracted from environmental samples (eDNA). Metabarcoding outputs DNA sequences that can be identified (annotated) by matching them against archived annotated sequences. However, sequences from most organisms are not archived - preventing annotation and potentially limiting metabarcoding in monitoring applications. Consequently, there is growing interest in using unannotated sequences as response metrics in monitoring programmes. We compared the sequences from three commonly used markers (16S (V3/V4 regions), 18S (V1/V2 regions) and COI) and, sampling along steep impact gradients, showed that the 16S and COI sequences were associated with the largest and smallest signal-to-noise ratio respectively. We trialled four separate, intuitive, noise-reduction approaches and demonstrated that removing less frequent sequences improved the signal-to-noise ratio, partitioning an additional 25 % from noise to explanatory factors in non-parametric ANOVA (NPA) and reducing dispersion in the data. For the 16S marker, retaining only the most frequently observed sequence, per sample, resulting in nine sequences across 150 samples, generated a near-maximal signal-to-noise ratio (95 % of the variance explained in NPA). We recommend that NPA, combined with rigorous elimination of less frequent sequences, be used to pre-filter sequences/taxa being used in monitoring applications. Our approach will simplify downstream analysis, for example the identification of key taxa and functional associations.

Detection of the threatened snail darter Percina tanasi in the Tennessee River system using environmental DNA

The distribution of many fishes that occupy large rivers is poorly known, in part due to the difficulties of sampling for them. This is especially true for small-bodied or rare species, such as the snail darter Percina tanasi, 44, 469-488; 1976). This federally listed (threatened) species has a limited distribution in the Tennessee River system in Alabama and Tennessee, where it is known from a few large tributaries or small rivers. In Alabama, P. tanasi was previously known from only one locality, but has recently been found in two additional, widely separated systems. These new records raise questions regarding the accuracy of our current understanding of the range for this species. Particularly, is P. tanasi present throughout the main stem Tennessee River, and is this species dispersing into new areas from source populations in the river? To clarify the distribution of P. tanasi in Alabama, 83 unique sites were surveyed using environmental DNA analysis. This cost-effective detection tool reduces the difficulty associated with empirically sampling large rivers for small fishes. Approximately 42% of sites sampled were positive for P. tanasi DNA. This study confirmed the known localities of P. tanasi in the Bear Creek, Elk River and Paint Rock River. Several new localities were also discovered throughout the main stem Tennessee River and in Shoal Creek, near Florence, Alabama. These findings can inform biologists about where to prioritize conservation efforts and further could lead to studies assessing movement and relatedness between populations in this system.

VTAM: A robust pipeline for validating metabarcoding data using controls

To obtain accurate estimates for biodiversity and ecological studies, metabarcoding studies should be carefully designed to minimize both false positive (FP) and false negative (FN) occurrences. Internal controls (mock samples and negative controls), replicates, and overlapping markers allow controlling metabarcoding errors but current metabarcoding software packages do not explicitly integrate these additional experimental data to optimize filtering. We have developed the metabarcoding analysis software VTAM, which uses explicitly these elements of the experimental design to find optimal parameter settings that minimize FP and FN occurrences. VTAM showed similar sensitivity, but a higher precision compared to two other pipelines using three datasets and two different markers (COI, 16S). The stringent filtering procedure implemented in VTAM aims to produce robust metabarcoding data to obtain accurate ecological estimates and represents an important step towards a non-arbitrary and standardized validation of metabarcoding data for conducting ecological studies. VTAM is implemented in Python and available from: https://github.com/aitgon/vtam. The VTAM benchmark code is available from: https://github.com/aitgon/vtam_benchmark.

Environmental DNA (eDNA) reveals potential for interoceanic fish invasions across the Panama Canal

Interoceanic canals can facilitate biological invasions as they connect the world's oceans and remove dispersal barriers between bioregions. As a consequence, multiple opportunities for biotic exchange arise and the resulting establishment of migrant species often causes adverse ecological and economic impacts. The Panama Canal is a key region for biotic exchange as it connects the Pacific and Atlantic Oceans in Central America. In this study, we used two complementary methods (environmental DNA (eDNA) metabarcoding and gillnetting) to survey fish communities in this unique waterway. Using COI (cytochrome oxidase subunit I) metabarcoding, we detected a total of 142 fish species, including evidence for the presence of sixteen Atlantic and eight Pacific marine fish in different freshwater sections of the Canal. Of these, nine are potentially new records. Molecular data did not capture all species caught with gillnets, but generally provided a more complete image of the known fish fauna as more small-bodied fish species were detected. Diversity indices based on eDNA surveys revealed significant differences across different sections of the Canal reflecting in part the prevailing environmental conditions. The observed increase in the presence of marine fish species in the Canal indicates a growing potential for interoceanic fish invasions. The potential ecological and evolutionary consequences of this increase in marine fishes are not only restricted to the fish fauna in the Canal as they could also impact adjacent ecosystems in the Pacific and Atlantic Oceans.

Optimal sample type and number vary in small shallow lakes when targeting non-native fish environmental DNA

Non-native fish have been shown to have deleterious impacts on freshwater ecosystems in New Zealand. Early detection is critical for their effective management. Traditional capture-based techniques may not detect newly introduced fish, especially if they are present in low abundance. Molecular techniques that target environmental DNA (eDNA) have been shown, in many instances, to be more sensitive, cost-effective and require lower sampling effort. However, appropriate sampling strategies are needed to ensure robust and interpretable data are obtained. In this study we used droplet digital PCR assays to investigate the presence of two non-native fish in New Zealand, the European perch (Perca fluviatilis) and rudd (Scardinius erythrophthalmus) in three small lakes. Samples were collected from water and surface sediment at near-shore and mid-lake sites. Probabilistic modelling was used to assess the occupancy of fish eDNA and develop guidance on sampling strategies. Based on the detection probability measures from the present study, at least six sites and five replicates per site are needed to reliably detect fish eDNA in sediment samples, and twelve sites with eight replicates per site for water samples. The results highlight the potential of developing monitoring and surveillance programs adapted to lakes, that include the use of assays targeting eDNA. This study focused on small shallow lakes, and it is likely that these recommendations may vary in larger, deeper, and more geomorphologically complex lakes, and this requires further research.

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.

The community ecology perspective of omics data

The measurement of uncharacterized pools of biological molecules through techniques such as metabarcoding, metagenomics, metatranscriptomics, metabolomics, and metaproteomics produces large, multivariate datasets. Analyses of these datasets have successfully been borrowed from community ecology to characterize the molecular diversity of samples (ɑ-diversity) and to assess how these profiles change in response to experimental treatments or across gradients (β-diversity). However, sample preparation and data collection methods generate biases and noise which confound molecular diversity estimates and require special attention. Here, we examine how technical biases and noise that are introduced into multivariate molecular data affect the estimation of the components of diversity (i.e., total number of different molecular species, or entities; total number of molecules; and the abundance distribution of molecular entities). We then explore under which conditions these biases affect the measurement of ɑ- and β-diversity and highlight how novel methods commonly used in community ecology can be adopted to improve the interpretation and integration of multivariate molecular data. Video Abstract.

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

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

Applications of environmental DNA (eDNA) in agricultural systems: Current uses, limitations and future prospects

Global food production, food supply chains and food security are increasingly stressed by human population growth and loss of arable land, becoming more vulnerable to anthropogenic and environmental perturbations. Numerous mutualistic and antagonistic species are interconnected with the cultivation of crops and livestock and these can be challenging to identify on the large scales of food production systems. Accurate identifications to capture this diversity and rapid scalable monitoring are necessary to identify emerging threats (i.e. pests and pathogens), inform on ecosystem health (i.e. soil and pollinator diversity), and provide evidence for new management practices (i.e. fertiliser and pesticide applications). Increasingly, environmental DNA (eDNA) is providing rapid and accurate classifications for specific organisms and entire species assemblages in substrates ranging from soil to air. Here, we aim to discuss how eDNA is being used for monitoring of agricultural ecosystems, what current limitations exist, and how these could be managed to expand applications into the future. In a systematic review we identify that eDNA-based monitoring in food production systems accounts for only 4 % of all eDNA studies. We found that the majority of these eDNA studies target soil and plant substrates (60 %), predominantly to identify microbes and insects (60 %) and are biased towards Europe (42 %). While eDNA-based monitoring studies are uncommon in many of the world's food production systems, the trend is most pronounced in emerging economies often where food security is most at risk. We suggest that the biggest limitations to eDNA for agriculture are false negatives resulting from DNA degradation and assay biases, as well as incomplete databases and the interpretation of abundance data. These require in silico, in vitro, and in vivo approaches to carefully design, test and apply eDNA monitoring for reliable and accurate taxonomic identifications. We explore future opportunities for eDNA research which could further develop this useful tool for food production system monitoring in both emerging and developed economies, hopefully improving monitoring, and ultimately food security.

Designing a surveillance program for early detection of alien plants and insects in Norway

Naturalized species of alien plants and animals comprise < 3% of biodiversity recorded in Norway but have had major impacts on natural ecosystems through displacement of native species. Encroachment of alien species has been especially problematic for coastal sites close to transport facilities and urban areas with high density housing. The goal of our field project was to design and test a surveillance program for early detection of alien species of vascular plants and terrestrial insects at the first phase of establishment in natural areas. In our 3-year project (2018–2020), we sampled 60 study plots in three counties in the Oslofjord region of southern Norway. Study plots (6.25 ha) were selected by two criteria: manual selection based on expert opinion (27 plots) or by random selection based on weights from a hotspot analysis of occurrence of alien species (33 plots). Vascular plants were surveyed by two experienced botanists who found a total of 239 alien species of vascular plants in 95 rounds of surveys. Insects and other invertebrates were captured with a single Malaise trap per site, with 3–4 rounds of repeated sampling. We used DNA-metabarcoding to identify invertebrates based on DNA extractions from crushed insects or from the preservative media. Over 3500 invertebrate taxa were detected in 255 rounds of sampling. We recorded 20 alien species of known risk, and 115 species that were new to Norway, including several ‘doorknocker’ species identified by previous risk assessments. We modeled the probabilities of occupancy (ψ) and detection (p) with occupancy models with repeated visits by multiple observers (vascular plants) or multiple rounds of sampling (insects). The two probabilities covaried with risk category for alien organisms and both were low for species categorized as no known or low risk (range = 0.052–0.326) but were higher for species categorized as severe risk (range = 0.318–0.651). Selecting sites at random or manually did not improve the probability of finding novel alien species, but occupancy had a weak positive relationship with housing density for some categories of alien plants and insects. We used our empirical estimates to test alternative sampling designs that would minimize the combined variance of occupancy and detection (A-optimality criterion). Sampling designs with 8–10 visits per site were best for surveillance of new alien species if the probabilities of occupancy and detection were both low, and provided low conditional probabilities of site occupancy ($$\hat{\psi }_{condl}$$ ψ ^ condl ≤ 0.032) and a high probabilities of cumulative detection ($$\hat{p}*$$ p ^ ∗  ≥ 0.943). Our field results demonstrate that early detection is feasible as a key component of a national surveillance program based on early detection and rapid response.

Environmental DNA-based profiling of benthic bacterial and eukaryote communities along a crude oil spill gradient in a coral reef in the Persian Gulf

Coral reef ecosystems in the Persian Gulf are frequently exposed to crude oil spills. We investigated benthic bacterial and eukaryote community structures at such coral reef sites subjected to different degrees of polycyclic aromatic hydrocarbon (PAH) pollution using environmental DNA (eDNA) metabarcoding. Both bacterial and eukaryote communities responded with pronounced shifts to crude oil pollution and distinguished control sites, moderately and heavily impacted sites with significant confidentiality. The observed community patterns were predominantly driven by Alphaproteobacteria and metazoans. Among these, we identified individual genera that were previously linked to oil spill stress, but also taxa, for which a link to hydrocarbon still remains to be established. Considering the lack of an early-warning system for the environmental status of coral reef ecosystems exposed to frequent crude-oil spills, our results encourage further research towards the development of an eDNA-based biomonitoring tool that exploits benthic bacterial and eukaryote communities as bioindicators.

Drastic reduction of false positive species in samples of insects by intersecting the default output of two popular metagenomic classifiers

The use of high-throughput sequencing to recover short DNA reads of many species has been widely applied on biodiversity studies, either as amplicon metabarcoding or shotgun metagenomics. These reads are assigned to taxa using classifiers. However, for different reasons, the results often contain many false positives. Here we focus on the reduction of false positive species attributable to the classifiers. We benchmarked two popular classifiers, BLASTn followed by MEGAN6 (BM) and Kraken2 (K2), to analyse shotgun sequenced artificial single-species samples of insects. To reduce the number of misclassified reads, we combined the output of the two classifiers in two different ways: (1) by keeping only the reads that were attributed to the same species by both classifiers (intersection approach); and (2) by keeping the reads assigned to some species by any classifier (union approach). In addition, we applied an analytical detection limit to further reduce the number of false positives species. As expected, both metagenomic classifiers used with default parameters generated an unacceptably high number of misidentified species (tens with BM, hundreds with K2). The false positive species were not necessarily phylogenetically close, as some of them belonged to different orders of insects. The union approach failed to reduce the number of false positives, but the intersection approach got rid of most of them. The addition of an analytic detection limit of 0.001 further reduced the number to ca. 0.5 false positive species per sample. The misidentification of species by most classifiers hampers the confidence of the DNA-based methods for assessing the biodiversity of biological samples. Our approach to alleviate the problem is straightforward and significantly reduced the number of reported false positive species.

Environmental DNA detects Spawning Habitat of an ephemeral migrant fish (Anadromous Rainbow Smelt: Osmerus mordax)

BACKGROUND

Anadromous rainbow smelt (Osmerus mordax) have experienced a large range reduction in recent decades and the status of remnant spawning populations is poorly known in Maine, where these fish have significant ecological, cultural, and commercial relevance. Defining the remnant range of anadromous smelt is more difficult than for many declining fish species because adults are only ephemerally present while spawning in small coastal streams at night during spring runoff periods when traditional assessments can be unreliable or even hazardous. We hypothesized that eDNA might facilitate improved survey efforts to define smelt spawning habitat, but that detection could also face challenges from adult eDNA quickly flushing out of these small stream systems. We combined daytime eDNA sampling with nighttime fyke netting to ascertain a potential window of eDNA detection before conducting eDNA surveys in four streams of varying abundance. Hierarchical occupancy modeling was in turn employed to estimate eDNA encounter probabilities relative to numbers of sampling events (date), samples within events, and qPCR replicates within samples.

RESULTS

Results from the combined eDNA and fyke net study indicated eDNA was detectable over an extended period, culminating approximately 8-13 days following peak spawning, suggesting developing smelt larvae might be the primary source of eDNA. Subsequently, smelt eDNA was readily detected in eDNA surveys of four streams, particularly following remediation of PCR inhibitors. Hierarchical occupancy modeling confirmed our surveys had high empirical detection for most sites, and that future surveys employing at least three sampling events, three samples per event, and six qPCR replicates can afford greater than 90% combined detection capability in low abundance systems.

CONCLUSIONS

These results demonstrate that relatively modest eDNA sampling effort has high capacity to detect this ephemerally present species of concern at low to moderate abundances. As such, smelt eDNA detection could improve range mapping by providing longer survey windows, safer sampling conditions, and lower field effort in low density systems, than afforded by existing visual and netting approaches.

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.

Postglacial species arrival and diversity buildup of northern ecosystems took millennia

What drives ecosystem buildup, diversity, and stability? We assess species arrival and ecosystem changes across 16 millennia by combining regional-scale plant sedimentary ancient DNA from Fennoscandia with near-complete DNA and trait databases. We show that postglacial arrival time varies within and between plant growth forms. Further, arrival times were mainly predicted by adaptation to temperature, disturbance, and light. Major break points in ecological trait diversity were seen between 13.9 and 10.8 calibrated thousand years before the present (cal ka BP), as well as break point in functional diversity at 12.0 cal ka BP, shifting from a state of ecosystem buildup to a state where most habitat types and biotic ecosystem components were in place. Trait and functional diversity stabilized around 8 cal ka BP, after which both remained stable, although changes in climate took place and species inflow continued. Our ecosystem reconstruction indicates a millennial-scale time phase of formation to reach stable and resilient levels of diversity and functioning.

Comparison of destructive and non-destructive DNA extraction methods for the metabarcoding of arthropod bulk samples

DNA metabarcoding is routinely used for biodiversity assessment, especially targeting highly diverse groups for which limited taxonomic expertise is available. Various protocols are currently in use, although standardization is key to its application in large-scale monitoring. DNA metabarcoding of arthropod bulk samples can be either conducted destructively from sample tissue, or non-destructively from sample fixative or lysis buffer. Non-destructive methods are highly desirable for the preservation of sample integrity but have yet to be experimentally evaluated in detail. Here, we compare diversity estimates from 14 size sorted Malaise trap samples processed consecutively with three non-destructive approaches (one using fixative ethanol and two using lysis buffers) and one destructive approach (using homogenized tissue). Extraction from commercial lysis buffer yielded comparable species richness and high overlap in species composition to the ground tissue extracts. A significantly divergent community was detected from preservative ethanol-based DNA extraction. No consistent trend in species richness was found with increasing incubation time in lysis buffer. These results indicate that non-destructive DNA extraction from incubation in lysis buffer could provide a comparable alternative to destructive approaches with the added advantage of preserving the specimens for post-metabarcoding taxonomic work but at a higher cost per sample.

Toward global integration of biodiversity big data: a harmonized metabarcode data generation module for terrestrial arthropods

Metazoan metabarcoding is emerging as an essential strategy for inventorying biodiversity, with diverse projects currently generating massive quantities of community-level data. The potential for integrating across such data sets offers new opportunities to better understand biodiversity and how it might respond to global change. However, large-scale syntheses may be compromised if metabarcoding workflows differ from each other. There are ongoing efforts to improve standardization for the reporting of inventory data. However, harmonization at the stage of generating metabarcode data has yet to be addressed. A modular framework for harmonized data generation offers a pathway to navigate the complex structure of terrestrial metazoan biodiversity. Here, through our collective expertise as practitioners, method developers, and researchers leading metabarcoding initiatives to inventory terrestrial biodiversity, we seek to initiate a harmonized framework for metabarcode data generation, with a terrestrial arthropod module. We develop an initial set of submodules covering the 5 main steps of metabarcode data generation: (i) sample acquisition; (ii) sample processing; (iii) DNA extraction; (iv) polymerase chain reaction amplification, library preparation, and sequencing; and (v) DNA sequence and metadata deposition, providing a backbone for a terrestrial arthropod module. To achieve this, we (i) identified key points for harmonization, (ii) reviewed the current state of the art, and (iii) distilled existing knowledge within submodules, thus promoting best practice by providing guidelines and recommendations to reduce the universe of methodological options. We advocate the adoption and further development of the terrestrial arthropod module. We further encourage the development of modules for other biodiversity fractions as an essential step toward large-scale biodiversity synthesis through harmonization.

Metabarcoding of Fish Larvae in the Merbok River Reveals Species Diversity and Distribution Along its Mangrove Environment

The Merbok River (north-west of Peninsular Malaysia) is a mangrove estuary that provides habitat for over 100 species of fish, which are economically and ecologically important. Threats such as habitat loss and overfishing are becoming a great concern for fisheries conservation and management. The identification of larval fish in this estuarine system is important to complement information on the adults. This is because the data could inform the spawning behaviour, reproductive biology, selection of nursery grounds and migration route of fish. Such information is invaluable for fisheries and aquatic environmental monitoring, and thus for their conservation and management. However, identifying fish larvae is a challenging task based only on morphology and even traditional DNA barcoding. To address this, DNA metabarcoding was utilised to detect the diversity of fish in the Merbok River. To complete the study, the fish larvae were collected at six sampling sites of the river. The extracted larval DNA was amplified for the Cytochrome Oxidase subunit 1 (COI) and 12S ribosomal RNA (12S rRNA) genes based on the metabarcoding approach using shotgun sequencing on the next-generation sequencing (NGS) Illumina MiSeq platform. Eighty-nine species from 65 genera and 41 families were detected, with Oryzias javanicus, Oryzias dancena, Lutjanus argentimaculatus and Lutjanus malabaricus among the most common species. The lower diversity observed from previous morphological studies is suggested to be mainly due to seasonal variation over the sampling period between the two methods and limited 12S rRNA sequences in current databases. The metabarcode data and a validation Sanger sequencing step using 15 species-specific primer pairs detected three species in common: Oryzias javanicus, Decapterus maruadsi and Pennahia macrocephalus. Several discrepancies observed between the two molecular approaches could be attributed to contaminants during sampling and DNA extraction, which could mask the presence of target species, especially when DNA from the contaminants is more abundant than the target organisms. In conclusion, this rapid and cost-effective identification method using DNA metabarcoding allowed the detection of numerous fish species from bulk larval samples in the Merbok River. This method can be applied to other sites and other organisms of interest.

Environmental RNA outperforms eDNA metabarcoding in assessing impact of marine pollution: A chromium-spiked mesocosm test

Environmental (e)DNA metabarcoding holds great promise for biomonitoring and ecotoxicological applications. However, few studies have compared the performance of eDNA versus eRNA metabarcoding in assessing organismal response to marine pollution, in experimental conditions. Here, we performed a chromium (Cr)-spiked mesocosm experimental test on benthic foraminiferal community to investigate the effects on species diversity by analysing both eDNA and eRNA metabarcoding data across different Cr concentrations in the sediment. Foraminiferal diversity in the eRNA data showed a significant negative correlation with the Cr concentration in the sediment, while a positive response was observed in the eDNA data. The foraminiferal OTUs exhibited a higher turnover rate in eRNA than in the eDNA-derived community. Furthermore, in the eRNA samples, OTUs abundance was significantly affected by the Cr gradient in the sediment (Pseudo-R² = 0.28, p = 0.05), while no significant trend was observed in the eDNA samples. The correlation between Cr concentration and foraminiferal diversity in eRNA datasets was stronger when the less abundant OTUs (<100 reads) were removed and the analyses were conducted exclusively on OTUs shared between eRNA and eDNA datasets. This indicates the importance of metabarcoding data filtering to capture ecological impacts, in addition to using the putatively active organisms in the eRNA dataset. The comparative analyses on foraminiferal diversity revealed that eRNA-based metabarcoding can better assess the response to heavy metal exposure in presence of subtle concentrations of the pollutant. Furthermore, our results suggest that to unlock the full potential for ecosystem assessment, eDNA and eRNA should be studied in parallel to control for potential sequence artifacts in routine ecosystem surveys.

Comparison of traditional and DNA metabarcoding samples for monitoring tropical soil arthropods (Formicidae, Collembola and Isoptera)

The soil fauna of the tropics remains one of the least known components of the biosphere. Long-term monitoring of this fauna is hampered by the lack of taxonomic expertise and funding. These obstacles may potentially be lifted with DNA metabarcoding. To validate this approach, we studied the ants, springtails and termites of 100 paired soil samples from Barro Colorado Island, Panama. The fauna was extracted with Berlese-Tullgren funnels and then either sorted with traditional taxonomy and known, individual DNA barcodes ("traditional samples") or processed with metabarcoding ("metabarcoding samples"). We detected 49 ant, 37 springtail and 34 termite species with 3.46 million reads of the COI gene, at a mean sequence length of 233 bp. Traditional identification yielded 80, 111 and 15 species of ants, springtails and termites, respectively; 98%, 37% and 100% of these species had a Barcode Index Number (BIN) allowing for direct comparison with metabarcoding. Ants were best surveyed through traditional methods, termites were better detected by metabarcoding, and springtails were equally well detected by both techniques. Species richness was underestimated, and faunal composition was different in metabarcoding samples, mostly because 37% of ant species were not detected. The prevalence of species in metabarcoding samples increased with their abundance in traditional samples, and seasonal shifts in species prevalence and faunal composition were similar between traditional and metabarcoding samples. Probable false positive and negative species records were reasonably low (13-18% of common species). We conclude that metabarcoding of samples extracted with Berlese-Tullgren funnels appear suitable for the long-term monitoring of termites and springtails in tropical rainforests. For ants, metabarcoding schemes should be complemented by additional samples of alates from Malaise or light traps.

Environmental DNA of insects and springtails from caves reveals complex processes of eDNA transfer in soils

Subterranean environments host a substantial amount of biodiversity, however assessing the distribution of species living underground is still extremely challenging. Environmental DNA (eDNA) metabarcoding is a powerful tool to estimate biodiversity in poorly known environments and has excellent performance for soil organisms. Here, we tested 1) whether eDNA metabarcoding from cave soils/sediments allows to successfully detect springtails (Hexapoda: Collembola) and insects (Hexapoda: Insecta); 2) whether eDNA mostly represents autochthonous (cave-dwelling) organisms or it also incorporates information from species living in surface environments; 3) whether eDNA detection probability changes across taxa with different ecology. Environmental DNA metabarcoding analyses detected a large number of Molecular Operational Taxonomic Units (MOTUs) for both insects and springtails. For springtails, detection probability was high, with a substantial proportion of hypogean species, suggesting that eDNA provides good information on the distribution of these organisms in caves. Conversely, for insects most of MOTUs represented taxa living outside caves, and the majority of them represented taxa/organisms living in freshwater environments (Ephemeroptera, Plecoptera and Trichoptera). The eDNA of freshwater insects was particularly abundant in deep sectors of caves, far from the entrance. Furthermore, average detection probability of insects was significantly lower than the one of springtails. This suggests that cave soils/sediments act as "conveyer belts of biodiversity information", possibly because percolating water lead to the accumulation of eDNA of organisms living in nearby areas. Cave soils hold a complex mix of autochthonous and allochthonous eDNA. eDNA provided unprecedented information on the understudied subterranean cave organisms; analyses of detection probability and occupancy can help teasing apart local eDNA from the eDNA representing spatially-integrated biodiversity for whole landscape.