Environmental DNA metabarcoding: Transforming how we survey animal and plant communities

Use of passive sampling in environmental DNA metabarcoding technology: Monitoring of fish diversity in the Jiangmen coastal waters

Environmental DNA (eDNA) metabarcoding technology is proving to be the most up-to-date and promising method for monitoring marine fish diversity. Fish eDNA is usually collected on a filter membrane after the filtration of water. Not only does this require the use of specialized equipment, but the amount of filtered water needed is also difficult to meet. The recently proposed passive eDNA collection method can expand the sampling scale, providing new perspectives for monitoring marine biodiversity. The role of collection methods in eDNA surveys, however, remains unclear. In this study, a low-cost custom framework with two types of filter membrane materials was used to conduct passive submersion samplings at the north and south ends of Shangchuan Island, Jiangmen, China. After defined periods of submersion, the filter membranes were recovered and eDNA extracted. Metabarcoding techniques were applied to detect fish species information in the eDNA samples. A total of 106 marine fish species from 27 orders, 53 families, and 92 genera, including one cartilaginous fish, were identified in the samples. The majority of fish detected by active filtration were also found in the passively collected samples, within the same location. Both sampling methods, therefore, showed similar species richness. Passive sampling was effective in identifying fish species diversity and provided a higher spatial resolution owing to the sample replicates. Passive sampling was also more sensitive in detecting species that differ significantly in abundance (biomarkers) between different sampling depths. When active filtration is not possible, or when large-scale sampling is the purpose of the study, passive sampling methods certainly provide a promising alternative. The findings of our study provide guidance for fish surveys and continuous bio-stereoscopic monitoring in coastal waters.

Freshwater eDNA reveals dramatic biological shifts linked to deforestation of New Zealand

Deforestation is considered a major threat to biodiversity across many parts of the globe, but the biological impacts of this dramatic ecosystem disturbance often remain incompletely understood. In New Zealand - the world's last major landmass to be colonised by humans - widespread deforestation over recent centuries has left a highly fragmented suite of relict forest stands, ideal for assessing anthropogenic biological change. We hypothesise that this widespread environmental disturbance has underpinned repeated and predictable ecological shifts across distinct rivers and regions. Here we use freshwater environmental DNA (eDNA) data (113 samples across 38 locations; 89 insect taxa) to test for concordant biological shifts linked to this deforestation. eDNA analyses highlight consistent compositional and functional differentiation between forested versus deforested assemblages, including turnover of 'cryptic' congeneric taxa that are morphologically similar yet ecologically and genetically distinct. These dramatic biological shifts are evident even over fine spatial scales within streams, emphasising the widespread emergence of a novel 'deforested' assemblage. Our results illustrate that environmental change can drive predictable biological shifts across broad geographic regions, and highlight the power of eDNA for assessing anthropogenic ecosystem change over large geographic scales.

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

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

Exploration and Retrieval of Virus-Related Molecular Data Using ExTaxsI: The Monkeypox Use Case

Retrieval and visualization of biological data are essential for understanding complex systems. With the increasing volume of data generated from high-throughput sequencing technologies, effective and optimized data visualization tools have become indispensable. This is particularly relevant in the COVID-19 postpandemic period, where understanding the diversity and interactions of microbial communities (i.e., viral and bacterial) constitutes an important asset to develop and plan suitable interventions.In this chapter, we show the usage and the potentials of ExTaxsI (Exploring Taxonomy Information) tool to retrieve viral biodiversity data stored in National Center for Biotechnology Information (NCBI) databases and create the related visualization. In addition, by integrating different functions and modules, the tool generates relevant types of visualization plots to facilitate the exploration of microbial biodiversity communities useful to deep dive into ecological and taxonomic relationships among different species and identify potential significant targets.Using the Monkeypox virus as a case study, this work points out significant perspectives on biological data visualization, which can be used to gain insights into the ecology, evolution, and pathogenesis of viruses. Accordingly, we show the potentiality of ExTaxsI to organize and describe the available/downloaded data in an easy, simple, and interpretable way allowing the user to interact dynamically with the visualization plots through specific filters, zoom, and explore functions.

Optimizing detectability of the endangered fan mussel using eDNA and ddPCR

Spatial and temporal monitoring of species threatened with extinction is of critical importance for conservation and ecosystem management. In the Mediterranean coast, the fan mussel (Pinna nobilis) is listed as critically endangered after suffering from a mass mortality event since 2016, leading to 100% mortality in most marine populations. Conventional monitoring for this macroinvertebrate is done using scuba, which is challenging in dense meadows or with low visibility. Here we developed an environmental DNA assay targeting the fan mussel and assessed the influence of several environmental parameters on the species detectability in situ. We developed and tested an eDNA molecular marker and collected 48 water samples in two sites at the Thau lagoon (France) with distinct fan mussel density, depths and during two seasons (summer and autumn). Our marker can amplify fan mussel DNA but lacks specificity since it also amplifies a conspecific species (Pinna rudis). We successfully amplified fan mussel DNA from in situ samples with 46 positive samples (out of 48) using ddPCR, although the DNA concentrations measured were low over almost all samples. Deeper sampling depth slightly increased DNA concentrations, but no seasonal effect was found. We highlight a putative spawning event on a single summer day with much higher DNA concentration compared to all other samples. We present an eDNA molecular assay able to detect the endangered fan mussel and provide guidelines to optimize the sampling protocol to maximize detectability. Effective and non-invasive monitoring tools for endangered species are promising to monitor remaining populations and have the potential of ecological restoration or habitat recolonization following a mass mortality event.

Environmental DNA biomonitoring reveals the human impacts on native and non-native fish communities in subtropical river systems

Subtropical rivers are one of the hotspots of global biodiversity, facing increased risks of fish diversity changes and species extinction. However, until now, human impacts on native and non-native fish communities in subtropical rivers still lack sufficient effort. Here, we used the environmental DNA (eDNA) approach to investigate fish communities in the Dongjiang River of southeast China, a typical subtropical river, and explored the effects of regional land use and local water pollution on fish taxonomic and functional diversity. Our data showed that 90 species or genera of native fish and 15 species or genera of non-native fish were detected by the eDNA approach, and there was over 85% overlap between eDNA datasets and historical records. The taxonomic and functional diversity of all, native and non-native fish communities showed consistent spatial patterns, that is, the upstream of the tributary was significantly higher than that of the mainstream and downstream. Land use and water pollution such as COD and TP were the determinants in shaping the spatial structure of fish communities, and water pollution explained 31.56%, 29.88%, and 27.80% of the structural variation in all, native and non-native fish communities, respectively. The Shannon diversity and functional richness of native fish showed a significant downward trend driven by COD (pShannon = 0.0374; pfunctional = 0.0215) and land use (pShannon = 0.0159; pfunctional = 0.0441), but they did not have significant impacts on non-native fish communities. Overall, this study emphasizes the inconsistent response of native and non-native fish communities to human impacts in subtropical rivers, and managers need to develop strategies tailored to specific fish species to effectively protect water security and rivers.

Assessment of the impact of dams on aquatic food webs using stable isotopes: Current progress and future challenges

Dams have disrupted natural river systems worldwide and although population and community level effects on aquatic biota have been well documented, food web responses remain poorly understood and difficult to characterize. The application of stable isotope analysis (SIA) provides a means to assess the effect of dams on food webs. Here we review the effect of dams on aquatic food webs using SIA, aiming to detect knowledge gaps in the field of dam impacts on aquatic food webs and propose a conceptual framework to help formulate hypotheses about dam impacts on food webs guided by food web theory. Dams can affect aquatic food webs via two pathways: a bottom-up pathway with altered basal food sources and their transfer to consumers through changes in flow, nutrients, temperature and sediment, and a top-down pathway with consumer species composition altered mainly through habitat fragmentation and related physiochemical changes. Taking these mechanisms into consideration, the impact of dams on food web attributes derived from SIA was evaluated. These studies generally apply mixing models to determine how dams alter the dominant carbon sources supporting food webs, use δ15N to examine how dams alter food-chain length, or use Layman metrics of isotope variability to assess niche changes for invertebrate and fish assemblages. Most studies compare the patterns of SIA metrics spatially (e.g. upstream vs reservoir vs downstream of dams; regulated vs unregulated rivers) and temporally (before vs after dam construction), without explicit hypotheses and/or links to theoretical concepts of food webs. We propose several steps to make SIA studies of dam impacts more rigorous and enhance their potential for producing novel insights. Future studies should quantify the shape and strength of the effect of dams on SIA-measured food web response, be conducted at larger temporal and spatial scales (particularly along the river longitudinal continuum and the lateral connected ecosystems (e.g., floodplains)), and consider effects of dams on food web resilience and tipping points.

Current status and topical issues on the use of eDNA-based targeted detection of rare animal species

Animal detection through DNA present in environmental samples (eDNA) is a valuable tool for detecting rare species, that are difficult to observe and monitor. eDNA-based tools are underpinned by molecular evolutionary principles, key to devising tools to efficiently single out a targeted species from an environmental sample. Here, we present a comprehensive review of the use of eDNA-based methods for the detection of targeted animal species, such as rare, endangered, or invasive species, through the analysis of 549 publications (2008-2022). Aquatic ecosystems have been the most surveyed, in particular, freshwaters (74 %), and to a less extent marine (14 %) and terrestrial systems (10 %). Vertebrates, in particular, fish (38 %), and endangered species, have been the focus of most of these studies, and Cytb and COI are the most employed markers. Among invertebrates, assays have been mainly designed for Mollusca and Crustacea species (21 %), in particular, to target invasive species, and COI the most employed marker. Targeted molecular approaches, in particular qPCR, have been the most adopted (75 %), while eDNA metabarcoding has been rarely used to target single or few species (approx. 6 %). However, less attention has been given in these studies to the effects of environmental factors on the amount of shed DNA, the differential amount of shed DNA among species, or the sensitivity of the markers developed, which may impact the design of the assays, particularly to warrant the required detection level and avoid false negatives and positives. The accuracy of the assays will also depend on the availability of genetic data and vouchered tissue or DNA samples from closely related species to assess both marker and primers' specificity. In addition, eDNA-based assays developed for a particular species may have to be refined for use in a new geographic area taking into account site-specific populations, as well as any intraspecific variation.

Inferring the ecological and evolutionary determinants of community genetic diversity

Understanding the relative contributions of ecological and evolutionary processes to the structuring of ecological communities is needed to improve our ability to predict how communities may respond to future changes in an increasingly human-modified world. Metabarcoding methods make it possible to gather population genetic data for all species within a community, unlocking a new axis of data to potentially unveil the origins and maintenance of biodiversity at local scales. Here, we present a new eco-evolutionary simulation model for investigating community assembly dynamics using metabarcoding data. The model makes joint predictions of species abundance, genetic variation, trait distributions and phylogenetic relationships under a wide range of parameter settings (e.g. high speciation/low dispersal or vice versa) and across a range of community states, from pristine and unmodified to heavily disturbed. We first demonstrate that parameters governing metacommunity and local community processes leave detectable signatures in simulated biodiversity data axes. Next, using a simulation-based machine learning approach we show that neutral and non-neutral models are distinguishable and that reasonable estimates of several model parameters within the local community can be obtained using only community-scale genetic data, while phylogenetic information is required to estimate those describing metacommunity dynamics. Finally, we apply the model to soil microarthropod metabarcoding data from the Troodos mountains of Cyprus, where we find that communities in widespread forest habitats are structured by neutral processes, while high-elevation and isolated habitats act as an abiotic filter generating non-neutral community structure. We implement our model within the ibiogen R package, a package dedicated to the investigation of island, and more generally community-scale, biodiversity using community-scale genetic data.

Environment and shipping drive environmental DNA beta-diversity among commercial ports

The spread of nonindigenous species by shipping is a large and growing global problem that harms coastal ecosystems and economies and may blur coastal biogeographical patterns. This study coupled eukaryotic environmental DNA (eDNA) metabarcoding with dissimilarity regression to test the hypothesis that ship-borne species spread homogenizes port communities. We first collected and metabarcoded water samples from ports in Europe, Asia, Australia and the Americas. We then calculated community dissimilarities between port pairs and tested for effects of environmental dissimilarity, biogeographical region and four alternative measures of ship-borne species transport risk. We predicted that higher shipping between ports would decrease community dissimilarity, that the effect of shipping would be small compared to that of environment dissimilarity and shared biogeography, and that more complex shipping risk metrics (which account for ballast water and stepping-stone spread) would perform better. Consistent with our hypotheses, community dissimilarities increased significantly with environmental dissimilarity and, to a lesser extent, decreased with ship-borne species transport risks, particularly if the ports had similar environments and stepping-stone risks were considered. Unexpectedly, we found no clear effect of shared biogeography, and that risk metrics incorporating estimates of ballast discharge did not offer more explanatory power than simpler traffic-based risks. Overall, we found that shipping homogenizes eukaryotic communities between ports in predictable ways, which could inform improvements in invasive species policy and management. We demonstrated the usefulness of eDNA metabarcoding and dissimilarity regression for disentangling the drivers of large-scale biodiversity patterns. We conclude by outlining logistical considerations and recommendations for future studies using this approach.

Metabarcoding and applied ecology with hyperdiverse organisms: Recommendations for biological control research

Metabarcoding is revolutionizing fundamental research in ecology by enabling large-scale detection of species and producing data that are rich with community context. However, the benefits of metabarcoding have yet to be fully realized in fields of applied ecology, especially those such as classical biological control (CBC) research that involve hyperdiverse taxa. Here, we discuss some of the opportunities that metabarcoding provides CBC and solutions to the main methodological challenges that have limited the integration of metabarcoding in existing CBC workflows. We focus on insect parasitoids, which are popular and effective biological control agents (BCAs) of invasive species and agricultural pests. Accurately identifying native, invasive and BCA species is paramount, since misidentification can undermine control efforts and lead to large negative socio-economic impacts. Unfortunately, most existing publicly accessible genetic databases cannot be used to reliably identify parasitoid species, thereby limiting the accuracy of metabarcoding in CBC research. To address this issue, we argue for the establishment of authoritative genetic databases that link metabarcoding data to taxonomically identified specimens. We further suggest using multiple genetic markers to reduce primer bias and increase taxonomic resolution. We also provide suggestions for biological control-specific metabarcoding workflows intended to track the long-term effectiveness of introduced BCAs. Finally, we use the example of an invasive pest, Drosophila suzukii, in a reflective "what if" thought experiment to explore the potential power of community metabarcoding in CBC.

Genomic insights into the conservation of wild and domestic animal diversity: A review

Due to environmental change and anthropogenic activities, global biodiversity has suffered an unprecedented loss, and the world is now heading toward the sixth mass extinction event. This urges the need to step up our efforts to promote the sustainable use of animal genetic resources and plan effective strategies for their conservation. Although habitat preservation and restoration are the primary means of conserving biodiversity, genomic technologies offer a variety of novel tools for identifying biodiversity hotspots and thus, support conservation efforts. Conservation genomics is a broad area of science that encompasses the application of genomic data from thousands or tens of thousands of genome-wide markers to address important conservation biology concerns. Genomic approaches have revolutionized the way we understand and manage animal populations, providing tools to identify and preserve unique genetic variants and alleles responsible for adaptive genetic variation, reducing the deleterious consequences of inbreeding, and increasing the adaptive potential of threatened species. The advancement of genomic technologies, particularly comparative genomic approaches, and the increased accessibility of genomic resources in the form of genome-enabled taxa for non-model organisms, provides a distinct advantage in defining conservation units over traditional genetics approaches. The objective of this review is to provide an exhaustive overview of the concept of conservation genomics, discuss the rationale behind the transition from conservation genetics to genomic approaches, and emphasize the potential applications of genomic techniques for conservation purposes. We also highlight interesting case studies in both livestock and wildlife species where genomic techniques have been used to accomplish conservation goals. Finally, we address some challenges and future perspectives in this field.

Environmental DNA-based biomonitoring of Cuban Crocodylus and their accompanying vertebrate fauna from Zapata Swamp, Cuba

Crocodylians globally face considerable challenges, including population decline and extensive habitat modification. Close monitoring of crocodylian populations and their habitats is imperative for the timely detection of population trends, especially in response to management interventions. Here we use eDNA metabarcoding to identify the Critically Endangered Crocodylus rhombifer and the Vulnerable C. acutus, as well as vertebrate community diversity, in Cuba's Zapata Swamp. We tested four different primer sets, including those used previously in Crocodylus population genetic and phylogenetic research, for their efficiency at detecting crocodylian eDNA. We detected C. rhombifer eDNA in 11 out of 15 sampled locations within its historical geographic distribution. We found that data analyses using the VertCOI primers and the mBRAVE bioinformatics pipeline were the most effective molecular marker and pipeline combination for identifying this species from environmental samples. We also identified 55 vertebrate species in environmental samples across the four bioinformatics pipelines- ~ 85% known to be present in the Zapata ecosystem. Among them were eight species previously undetected in the area and eight alien species, including known predators of hatchling crocodiles (e.g., Clarias sp.) and egg predators (e.g., Mus musculus). This study highlights eDNA metabarcoding as a powerful tool for crocodylian biomonitoring within fragile and diverse ecosystems, particularly where fast, non-invasive methods permit detection in economically important areas and will lead to a better understanding of complex human-crocodile interactions and evaluate habitat suitability for potential reintroductions or recovery programs for threatened crocodylian species.

Investigating non-fungal eukaryotic diversity in snow in the Antarctic Peninsula region using DNA metabarcoding

Snow is a unique microhabitat, despite being a harsh environment, multiple life forms have adapted to survive in it. While algae, bacteria and fungi are dominant microorganisms in Antarctic snow, little is known about other organisms that may be present in this habitat. We used metabarcoding to investigate DNA sequence diversity of non-fungal eukaryotes present in snow obtained from six different sites across the Maritime Antarctica. A total of 20 taxa were assigned to obtained sequences, representing five Kingdoms (Chromista, Protozoa, Viridiplantae and Metazoa) and four phyla (Ciliophora, Cercozoa, Chlorophyta and Cnidaria). The highest diversity indices were detected in Trinity Peninsula followed by Robert Island, Arctowski Peninsula, Deception Island, King George Island and Snow Island. The most abundant assignments were to Trebouxiophyceae, followed by Chlamydomonas nivalis and Chlamidomonadales. No taxa were detected at all sites. Three potentially new records for Antarctica were detected: two Ciliophora (Aspidisca magna and Stokesia sp.) and the green algae Trebouxia potteri. Our data suggested that similarities found between the sites may be more related with snow physicochemical properties rather than geographic proximity or latitude. This study provides new insights into the diversity and distribution of eukaryotic organisms in Antarctic snow.

Unraveling the potential of environmental DNA for deciphering recent advances in plant-animal interactions: a systematic review

MAIN CONCLUSION

Environmental DNA-based monitoring provides critical insights for enhancing our understanding of plant-animal interactions in the context of worldwide biodiversity decrease for developing a global framework for effective plant biodiversity conservation. To understand the ecology and evolutionary patterns of plant-animal interactions (PAI) and their pivotal roles in ecosystem functioning advances in molecular ecology tools such as Environmental DNA (eDNA) provide unprecedented research avenues. These methods being non-destructive in comparison to traditional biodiversity monitoring methods, enhance the discernment of ecosystem health, integrity, and complex interactions. This review intends to offer a systematic and critical appraisal of the prospective of eDNA for investigating PAI. The review thoroughly discusses and analyzes the recent reports (2015-2022) employing preferred reporting items for systematic reviews and meta-analyses (PRISMA) to outline the recent progression in eDNA approaches for elucidating PAI. The current review envisages that eDNA has a significant potential to monitor both plants and associated cohort of prospective pollinators (avian species and flowering plants, bees and plants, arthropods and plants, bats and plants, etc.). Furthermore, a brief description of the factors that influence the utility and interpretation of PAI eDNA is also presented. The review establishes that factors such as biotic and abiotic, primer selection and taxonomic resolution, and indeterminate spatio-temporal scales impact the availability and longevity of eDNA. The study also identified the limitations that influence PAI detection and suggested possible solutions for better execution of these molecular approaches. Overcoming these research caveats will augment the assortment of PAI analysis through eDNA that could be vital for ecosystem health and integrity. This review forms a critical guide and offers prominent insights for ecologists, environmental managers and researchers to assess and evaluate plant-animal interaction through environmental DNA.

Assessment of multiple fecal contamination sources in surface waters using environmental mitochondrial DNA metabarcoding

Waterborne diseases are transmitted to humans through the fecal contamination of water, where homeothermic species are the main reservoir. Fecal indicator bacteria (FIB) are often used to determine the occurrence of fecal contamination. However, FIB cannot provide the source of fecal contamination. Furthermore, as fecal inputs and contamination could originate from multiple sources (e.g., human, livestock, wildlife), multiple source tracking markers are required to identify fecal sources. From a previous study, we developed a mitochondrial DNA (mtDNA) metabarcoding approach to assess the presence of multiple homeotherms in four surface waters. Here, we have broadened our approach by sampling 86 surface water samples from the L'Assomption River and Ville-Marie watersheds (Province of Quebec, Canada). Fecal coliform levels were higher than the expected sanitary recommendations for recreational water (> 200 CFU/100 mL) in 73 % samples. The occurrence of mtDNA from human, livestock, domestic animals, wild mammals and wild birds was found in 40-88 % of the samples. Multivariate analyses showed significant covariations between homeothermic taxa and fecal coliforms, enterococci, β-D-glucuronidase, conductivity, the human-specific Bacteroidales Hf183 genetic marker, and the human population, in the watersheds of L'Assomption River (p = 0.001) and Ville-Marie (p = 0.015) (Province of Quebec, Canada). Through the application of Bayes Theorem, it was determined that fecal coliforms co-occurred with the detection of bovine, beaver, robin and chicken mtDNA in 100 % of cases in the L'Assomption River watershed, and human mtDNA co-occurred with fecal coliforms in 93 % and 76 % of cases in L'Assomption River watershed and Ville-Marie sub-catchment, respectively. This study suggests that fecal contamination could be the result of multiple species, among which some wild animals may contribute to fecal inputs in surface waters, resulting in potential risk to human health. This reinforces the necessity of using the mtDNA metabarcoding method to monitor multi-animal species.

Insights in Pharmaceutical Pollution: The Prospective Role of eDNA Metabarcoding

Environmental pollution is a growing threat to natural ecosystems and one of the world's most pressing concerns. The increasing worldwide use of pharmaceuticals has elevated their status as significant emerging contaminants. Pharmaceuticals enter aquatic environments through multiple pathways related to anthropogenic activity. Their high consumption, insufficient waste treatment, and the incapacity of organisms to completely metabolize them contribute to their accumulation in aquatic environments, posing a threat to all life forms. Various analytical methods have been used to quantify pharmaceuticals. Biotechnology advancements based on next-generation sequencing (NGS) techniques, like eDNA metabarcoding, have enabled the development of new methods for assessing and monitoring the ecotoxicological effects of pharmaceuticals. eDNA metabarcoding is a valuable biomonitoring tool for pharmaceutical pollution because it (a) provides an efficient method to assess and predict pollution status, (b) identifies pollution sources, (c) tracks changes in pharmaceutical pollution levels over time, (d) assesses the ecological impact of pharmaceutical pollution, (e) helps prioritize cleanup and mitigation efforts, and (f) offers insights into the diversity and composition of microbial and other bioindicator communities. This review highlights the issue of aquatic pharmaceutical pollution while emphasizing the importance of using modern NGS-based biomonitoring actions to assess its environmental effects more consistently and effectively.

New deep learning-based methods for visualizing ecosystem properties using environmental DNA metabarcoding data

Environmental DNA (eDNA) metabarcoding provides an efficient approach for documenting biodiversity patterns in marine and terrestrial ecosystems. The complexity of these data prevents current methods from extracting and analyzing all the relevant ecological information they contain, and new methods may provide better dimensionality reduction and clustering. Here we present two new deep learning-based methods that combine different types of neural networks (NNs) to ordinate eDNA samples and visualize ecosystem properties in a two-dimensional space: the first is based on variational autoencoders and the second on deep metric learning. The strength of our new methods lies in the combination of two inputs: the number of sequences found for each molecular operational taxonomic unit (MOTU) detected and their corresponding nucleotide sequence. Using three different datasets, we show that our methods accurately represent several biodiversity indicators in a two-dimensional latent space: MOTU richness per sample, sequence α-diversity per sample, Jaccard's and sequence β-diversity between samples. We show that our nonlinear methods are better at extracting features from eDNA datasets while avoiding the major biases associated with eDNA. Our methods outperform traditional dimension reduction methods such as Principal Component Analysis, t-distributed Stochastic Neighbour Embedding, Nonmetric Multidimensional Scaling and Uniform Manifold Approximation and Projection for dimension reduction. Our results suggest that NNs provide a more efficient way of extracting structure from eDNA metabarcoding data, thereby improving their ecological interpretation and thus biodiversity monitoring.

Creating, curating and evaluating a mitogenomic reference database to improve regional species identification using environmental DNA

Species detection using eDNA is revolutionizing global capacity to monitor biodiversity. However, the lack of regional, vouchered, genomic sequence information-especially sequence information that includes intraspecific variation-creates a bottleneck for management agencies wanting to harness the complete power of eDNA to monitor taxa and implement eDNA analyses. eDNA studies depend upon regional databases of mitogenomic sequence information to evaluate the effectiveness of such data to detect and identify taxa. We created the Oregon Biodiversity Genome Project to create a database of complete, nearly error-free mitogenomic sequences for all of Oregon's fishes. We have successfully assembled the complete mitogenomes of 313 specimens of freshwater, anadromous and estuarine fishes representing 24 families, 55 genera and 129 species and lineages. Comparative analyses of these sequences illustrate that many regions of the mitogenome are taxonomically informative, that the short (~150 bp) mitochondrial 'barcode' regions typically used for eDNA assays do not consistently diagnose for species and that complete single or multiple genes of the mitogenome are preferable for identifying Oregon's fishes. This project provides a blueprint for other researchers to follow as they build regional databases, illustrates the taxonomic value and limits of complete mitogenomic sequences and offers clues as to how current eDNA assays and environmental genomics methods of the future can best leverage this information.

The secret hidden in dust: Assessing the potential to use biological and chemical properties of the airborne fraction of soil for provenance assignment and forensic casework

The airborne fraction of soil (dust) is both ubiquitous in nature and contains localised biological and chemical signatures, making it a potential medium for forensic intelligence. Metabarcoding of dust can yield biological communities unique to the site of interest, similarly, geochemical analyses can uncover elements and minerals within dust that can be matched to a geographic location. Combining these analyses presents multiple lines of evidence as to the origin of dust collected from items of interest. In this work, we investigated whether bacterial and fungal communities in dust change through time and whether they are comparable to soil samples of the same site. We integrated dust metabarcoding into a framework amenable to forensic casework, (i.e., using calibrated log-likelihood ratios) to predict the origin of dust samples using models constructed from both dust samples and soil samples from the same site. Furthermore, we tested whether both metabarcoding and geochemical/mineralogical analyses could be conducted on a single swabbed sample, for situations where sampling is limited. We found both analyses could generate results from a single swabbed sample and found biological and chemical signatures unique to sites. However, we did find significant variation within sites, where this did not always correlate with time but was a random effect of sampling. This variation within sites was not greater than between sites and so did not influence site discrimination. When modelling bacterial and fungal diversity using calibrated log-likelihood ratios, we found samples were correctly predicted using dust 67% and 56% of the time and using soil 56% and 22% of the time for bacteria and fungi communities respectively. Incorrect predictions were related to within site variability, highlighting limitations to assigning dust provenance using metabarcoding of soil.

Integrating citizen science and environmental DNA metabarcoding to study biodiversity of groundwater amphipods in Switzerland

Groundwater is the physically largest freshwater ecosystem, yet one of the least explored habitats on earth, both because of accessing difficulties and the scarcity of the organisms inhabiting it. Here, we demonstrate how a two-fold approach provides complementary information on the occurrence and diversity of groundwater amphipods. Firstly, we used a citizen science approach in collaboration with municipal water providers who sampled groundwater organisms in their spring catchment boxes over multiple weeks, followed by DNA barcoding. Secondly, we collected four 10 L water samples at each site, in one sampling event, for environmental DNA (eDNA) metabarcoding. We found that citizen science was very effective in describing the distribution and abundance of groundwater amphipods. Although the single time-point of eDNA sampling did not detect as many amphipods, it allowed the assessment of the entire groundwater community, including microorganisms. By combining both methods, we found different amphipod species co-occurring with distinct sequences from the eDNA-metabarcoding dataset, representing mainly micro-eukaryotic species. We also found a distinct correlation between the diversity of amphipods and the overall biodiversity of groundwater organisms detected by eDNA at each site. We thus suggest that these approaches can be used to get a better understanding of subterranean biodiversity.

Environmental DNA metabarcoding reveals distinct fish assemblages supported by seagrass (Zostera marina and Zostera pacifica) beds in different geographic settings in Southern California

Seagrass beds are disappearing at a record pace despite their known value to our oceans and coastal communities. Simultaneously, our coastlines are under the constant pressure of climate change which is impacting their chemical, physical and biological characteristics. It is thus pertinent to evaluate and record habitat use so we can understand how these different environments contribute to local biodiversity. This study evaluates the assemblages of fish found at five Zostera beds in Southern California using environmental DNA (eDNA) metabarcoding. eDNA is a powerful biodiversity monitoring tool that offers key advantages to conventional monitoring. Results from our eDNA study found 78 species of fish that inhabit these five beds around Southern California representing embayment, open coastal mainland and open coastal island settings. While each bed had the same average number of species found throughout the year, the composition of these fish assemblages was strongly site dependent. There were 35 fish that were found at both open coast and embayment seagrass beds, while embayment seagrass sites had 20 unique fish and open coast sites had 23 unique fish. These results demonstrate that seagrass fish assemblages are heterogenous based on their geographic positioning and that marine managers must take this into account for holistic conservation and restoration efforts.

Microbial Community Colonization Process Unveiled through eDNA-PFU Technology in Mesocosm Ecosystems

Microbial communities are essential components of aquatic ecosystems and are widely employed for the detection, protection, and restoration of water ecosystems. The polyurethane foam unit (PFU) method, an effective and widely used environmental monitoring technique, has been improved with the eDNA-PFU method, offering efficiency, rapidity, and standardization advantages. This research aimed to explore the colonization process of microbial communities within PFUs using eDNA-PFU technology. To achieve this, we conducted ten-day monitoring and sequencing of microbial communities within PFUs in a stable and controlled artificial aquatic ecosystem, comparing them with water environmental samples (eDNA samples). Results showed 1065 genera in eDNA-PFU and 1059 in eDNA, with eDNA-PFU detecting 99.95% of eDNA-identified species. Additionally, the diversity indices of bacteria and eukaryotes in both methods showed similar trends over time in the colonization process; however, relative abundance differed. We further analyzed the colonization dynamics of microbes in eDNA-PFU and identified four clusters with varying colonization speeds. Notably, we found differences in colonization rates between bacteria and eukaryotes. Furthermore, the Molecular Ecological Networks (MEN) showed that the network in eDNA-PFU was more modular, forming a unique microbial community differentiated from the aquatic environment. In conclusion, this study, using eDNA-PFU, comprehensively explored microbial colonization and interrelationships in a controlled mesocosm system, providing foundational data and reference standards for its application in aquatic ecosystem monitoring and beyond.

Environmental transcriptomics under heat stress: Can environmental RNA reveal changes in gene expression of aquatic organisms?

To safeguard biodiversity in a changing climate, taxonomic information about species turnover and insights into the health of organisms are required. Environmental DNA approaches are increasingly used for species identification, but cannot provide functional insights. Transcriptomic methods reveal the physiological states of macroorganisms, but are currently species-specific and require tissue sampling or animal sacrifice, making community-wide assessments challenging. Here, we test whether broad functional information (expression level of the transcribed genes) can be harnessed from environmental RNA (eRNA), which includes extra-organismal RNA from macroorganisms along with whole microorganisms. We exposed Daphnia pulex as well as phytoplankton prey and microorganism colonizers to control (20°C) and heat stress (28°C) conditions for 7 days. We sequenced eRNA from tank water (after complete removal of Daphnia) as well as RNA from Daphnia tissue, enabling comparisons of extra-organismal and organismal RNA-based gene expression profiles. Both RNA types detected similar heat stress responses of Daphnia. Using eRNA, we identified 32 Daphnia genes to be differentially expressed following heat stress. Of these, 17 were also differentially expressed and exhibited similar levels of relative expression in organismal RNA. In addition to the extra-organismal Daphnia response, eRNA detected community-wide heat stress responses consisting of distinct functional profiles and 121 differentially expressed genes across eight taxa. Our study demonstrates that environmental transcriptomics based on extra-organismal eRNA can noninvasively reveal gene expression responses of macroorganisms following environmental changes, with broad potential implications for the biomonitoring of health across the trophic chain.

SPECIES IN THE FECES: DNA METABARCODING TO DETECT POTENTIAL GASTROPOD HOSTS OF PARELAPHOSTRONGYLUS TENUIS Consumed By Moose (Alces Alces)

Our understanding of wildlife multihost pathogen transmission systems is often incomplete due to the difficulty of observing contact between hosts. Understanding these interactions can be critical for preventing disease-induced wildlife declines. The proliferation of high-throughput sequencing technologies provides new opportunities to better explore these cryptic interactions. Parelaphostrongylus tenuis, a multihost parasite, is a leading cause of death in some moose (Alces alces) populations threatened by local extinction in the midwestern and northeastern US and southeastern Canada. Moose contract P. tenuis by consuming infected gastropod intermediate hosts, but little is known about which gastropod species moose consume. To gain more insight, we used a genetic metabarcoding approach on 258 georeferenced and temporally stratified moose fecal samples collected May-October 2017-20 from a declining population in the north-central US. We detected moose consumption of three species of gastropods across five positive samples. Two of these (Punctum minutissimum and Helisoma sp.) have been minimally investigated for the ability to host P. tenuis, while one (Zonitoides arboreus) is a well-documented host. Moose consumption of gastropods documented herein occurred in June and September. Our findings prove that moose consume gastropod species known to become infected by P. tenuis and demonstrate that fecal metabarcoding can provide novel insight on interactions between hosts of a multispecies pathogen transmission system. After determining and improving test sensitivity, these methods may also be extended to document important interactions in other multihost disease systems.

eDNA-based diversity and multitrophic network reveal the effects of land use and pollutants on the subtropical Dongjiang River systems

Biodiversity and its constituted multitrophic network in rivers are accelerating change under human land use and pollutants. However, due to the lack of complete datasets across taxa limited by traditional morphological biomonitoring, the change patterns of biodiversity and multitrophic networks are still unclear. Here, we used the eDNA approach to capture multitrophic communities (including fish, aquatic insects, protozoa, diatom and bacteria) in the Dongjiang River, a typical subtropical river in southeast China, and analyzed the changing patterns of biodiversity and multitrophic networks in relation to land use and water pollution. First, our data showed that the eDNA approach provided a snapshot of the multitrophic communities in the Dongjiang River, and the monitored 5833 OTUs were annotated to 55 phyla, 144 classes, 329 orders, 521 families, 945 genera and 406 species. Second, the multitrophic diversity index had similar patterns on the longitudinal scale of rivers, with significant decreases from the upstream to the downstream, while individual taxonomic groups exhibited variable spatial patterns. While there were similar spatial patterns between network metrics and diversity index, the former had stronger relationships with the spatial distance. Third, the multitrophic diversity and networks were significantly negatively correlated with land use and water pollution (e.g., CODMn), and network structures often had stronger and non-linear responses. Overall, this study highlights that eDNA biomonitoring of multitrophic communities and networks can provide deeper insights into ecosystem changes and help develop more targeted management strategies.

Composition and diversity of bacterial communities associated with honey bee foragers from two contrasting environments

The honey bee is associated with a diverse community of microbes (viruses, bacteria, fungi, and protists), commonly known as the microbiome. Here, we present data on honey bee microbiota from two localities having different surrounding landscapes - mountain (the Rhodope Mountains) and lowland (the Danube plain). The bacterial communities of abdomen of adult bees were studied using amplicon sequencing of the 16S rRNA gene. The composition and dominance structure and their variability within and between localities, alpha and beta diversity, and core and differential taxa were compared at different hierarchical levels (operational taxonomic units to phylum). Seven genera (Lactobacillus, Gilliamella, Bifidobacterium, Commensalibacter, Bartonella, Snodgrassella, and Frischella), known to include core gut-associated phylotypes or species clusters, dominated (92-100%) the bacterial assemblages. Significant variations were found in taxa distribution across both geographical regions and within each apiary. Lactobacillus (Firmicutes) prevailed significantly in the mountain locality followed by Gilliamella and Bartonella (Proteobacteria). Bacteria of four genera, core (Bartonella and Lactobacillus) and non-core (Pseudomonas and Morganella), dominated the bee-associated assemblages of the Danube plain locality. Several ubiquitous bacterial genera (e.g., Klebsiella, Serratia, and Providencia), some species known also as potential and opportunistic bee pathogens, had been found in the lowland locality. Beta diversity analyses confirmed the observed differences in the bacterial communities from both localities. The occurrence of non-core taxa contributes substantially to higher microbial richness and diversity in bees from the Danube plain locality. We assume that the observed differences in the microbiota of honey bees from both apiaries are due to a combination of factors specific for each region. The surrounding landscape features of both localities and related vegetation, anthropogenic impact and land use intensity, the beekeeping management practices, and bee health status might all contribute to observed differences in bee microbiota traits.

Field and laboratory guidelines for reliable bioinformatic and statistical analysis of bacterial shotgun metagenomic data

Shotgun metagenomics is an increasingly cost-effective approach for profiling environmental and host-associated microbial communities. However, due to the complexity of both microbiomes and the molecular techniques required to analyze them, the reliability and representativeness of the results are contingent upon the field, laboratory, and bioinformatic procedures employed. Here, we consider 15 field and laboratory issues that critically impact downstream bioinformatic and statistical data processing, as well as result interpretation, in bacterial shotgun metagenomic studies. The issues we consider encompass intrinsic properties of samples, study design, and laboratory-processing strategies. We identify the links of field and laboratory steps with downstream analytical procedures, explain the means for detecting potential pitfalls, and propose mitigation measures to overcome or minimize their impact in metagenomic studies. We anticipate that our guidelines will assist data scientists in appropriately processing and interpreting their data, while aiding field and laboratory researchers to implement strategies for improving the quality of the generated results.

Quantification by droplet digital PCR and species identification by metabarcoding of environmental (e)DNA from Blainville's beaked whales, with assisted localization from an acoustic array

Detection and identification of species, subspecies or stocks of whales, dolphins and porpoises at sea remain challenging, particularly for cryptic or elusive species like beaked whales (Family: Ziphiidae). Here we investigated the potential for using an acoustically assisted sampling design to collect environmental (e)DNA from beaked whales on the U.S. Navy's Atlantic Undersea Test and Evaluation Center (AUTEC) in The Bahamas. During 12 days of August 2019, we conducted 9 small-boat surveys and collected 56 samples of seawater (paired subsamples of 1L each, including controls) using both a spatial collection design in the absence of visual confirmation of whales, and a serial collection design in the proximity of whales at the surface. There were 7 sightings of whales, including 11 Blainville's beaked whales (Mesoplodon densirostris). All whales were located initially with the assistance of information from a bottom-mounted acoustic array available on the AUTEC range. Quantification by droplet digital (dd)PCR from the four spatial design collections showed no samples of eDNA above the threshold of detection and none of these 20 samples yielded amplicons for conventional or next-generation sequencing. Quantification of the 31 samples from four serial collections identified 11 likely positive detections. eDNA barcoding by conventional sequencing and eDNA metabarcoding by next-generation sequencing confirmed species identification for 9 samples from three of the four serial collections. We further resolved five intra-specific variants (i.e., haplotypes), two of which showed an exact match to previously published haplotypes and three that have not been reported previously to the international repository, GenBank. A minimum spanning network of the five eDNA haplotypes, with all other published haplotypes of Blainville's beaked whales, suggested the potential for further resolution of differences between oceanic populations.

Interpreting phylogenetic placements for taxonomic assignment of environmental DNA

Taxonomic assignment of operational taxonomic units (OTUs) is an important bioinformatics step in analyzing environmental sequencing data. Pairwise alignment and phylogenetic-placement methods represent two alternative approaches to taxonomic assignments, but their results can differ. Here we used available colpodean ciliate OTUs from forest soils to compare the taxonomic assignments of VSEARCH (which performs pairwise alignments) and EPA-ng (which performs phylogenetic placements). We showed that when there are differences in taxonomic assignments between pairwise alignments and phylogenetic placements at the subtaxon level, there is a low pairwise similarity of the OTUs to the reference database. We then showcase how the output of EPA-ng can be further evaluated using GAPPA to assess the taxonomic assignments when there exist multiple equally likely placements of an OTU, by taking into account the sum over the likelihood weights of the OTU placements within a subtaxon, and the branch distances between equally likely placement locations. We also inferred the evolutionary and ecological characteristics of the colpodean OTUs using their placements within subtaxa. This study demonstrates how to fully analyze the output of EPA-ng, by using GAPPA in conjunction with knowledge of the taxonomic diversity of the clade of interest.

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.

A combination of machine-learning and eDNA reveals the genetic signature of environmental change at the landscape levels

The current advances of environmental DNA (eDNA) bring profound changes to ecological monitoring and provide unique insights on the biological diversity of ecosystems. The very nature of eDNA data is challenging yet also revolutionizing how biological monitoring information is analysed. In particular, new metrics and approaches should take full advantage of the extent and detail of molecular data produced by genetic methods. In this perspective, machine learning algorithms are particularly promising as they can capture complex relationships between the multiple environmental pressures and the diversity of biological communities. We investigated the potential of a new generation of biomonitoring tools that implement machine-learning techniques to fully exploit eDNA datasets. We trained a machine learning model to discriminate between reference and impacted communities of freshwater macroinvertebrates and assessed its performances using a large eDNA dataset collected at 64 standard federal monitoring sites across Switzerland. We show that a model trained on eDNA is significantly better than a naive model and performs similarly to a model trained on traditional data. Our proof-of-concept shows that such a combination of eDNA and machine learning approaches has the potential to complement or even replace traditional environmental monitoring, and could be scaled along temporal or spatial dimensions.

Trapped DNA fragments in marine sponge specimens unveil North Atlantic deep-sea fish diversity

Sponges pump water to filter feed and for diffusive oxygen uptake. In doing so, trace DNA fragments from a multitude of organisms living around them are trapped in their tissues. Here we show that the environmental DNA retrieved from archived marine sponge specimens can reconstruct the fish communities at the place of sampling and discriminate North Atlantic assemblages according to biogeographic region (from Western Greenland to Svalbard), depth habitat (80-1600 m), and even the level of protection in place. Given the cost associated with ocean biodiversity surveys, we argue that targeted and opportunistic sponge samples - as well as the specimens already stored in museums and other research collections - represent an invaluable trove of biodiversity information that can significantly extend the reach of ocean monitoring.

THAPBI PICT-a fast, cautious, and accurate metabarcoding analysis pipeline

THAPBI PICT is an open source software pipeline for metabarcoding analysis of Illumina paired-end reads, including cases of multiplexing where more than one amplicon is amplified per DNA sample. Initially a Phytophthora ITS1 Classification Tool (PICT), we demonstrate using worked examples with our own and public data sets how, with appropriate primer settings and a custom database, it can be applied to other amplicons and organisms, and used for reanalysis of existing datasets. The core dataflow of the implementation is (i) data reduction to unique marker sequences, often called amplicon sequence variants (ASVs), (ii) dynamic thresholds for discarding low abundance sequences to remove noise and artifacts (rather than error correction by default), before (iii) classification using a curated reference database. The default classifier assigns a label to each query sequence based on a database match that is either perfect, or a single base pair edit away (substitution, deletion or insertion). Abundance thresholds for inclusion can be set by the user or automatically using per-batch negative or synthetic control samples. Output is designed for practical interpretation by non-specialists and includes a read report (ASVs with classification and counts per sample), sample report (samples with counts per species classification), and a topological graph of ASVs as nodes with short edit distances as edges. Source code available from https://github.com/peterjc/thapbi-pict/ with documentation including installation instructions.

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.

Catchment-based sampling of river eDNA integrates terrestrial and aquatic biodiversity of alpine landscapes

Monitoring of terrestrial and aquatic species assemblages at large spatial scales based on environmental DNA (eDNA) has the potential to enable evidence-based environmental policymaking. The spatial coverage of eDNA-based studies varies substantially, and the ability of eDNA metabarcoding to capture regional biodiversity remains to be assessed; thus, questions about best practices in the sampling design of entire landscapes remain open. We tested the extent to which eDNA sampling can capture the diversity of a region with highly heterogeneous habitat patches across a wide elevation gradient for five days through multiple hydrological catchments of the Swiss Alps. Using peristaltic pumps, we filtered 60 L of water at five sites per catchment for a total volume of 1800 L. Using an eDNA metabarcoding approach focusing on vertebrates and plants, we detected 86 vertebrate taxa spanning 41 families and 263 plant taxa spanning 79 families across ten catchments. For mammals, fishes, amphibians and plants, the detected taxa covered some of the most common species in the region according to long-term records while including a few more rare taxa. We found marked turnover among samples from distinct elevational classes indicating that the biological signal in alpine rivers remains relatively localised and is not aggregated downstream. Accordingly, species compositions differed between catchments and correlated with catchment-level forest and grassland cover. Biomonitoring schemes based on capturing eDNA across rivers within biologically integrated catchments may pave the way toward a spatially comprehensive estimation of biodiversity.

Monitoring mosquito richness in an understudied area: can environmental DNA metabarcoding be a complementary approach to adult trapping?

Mosquito surveillance programmes are essential to assess the risks of local vector-borne disease outbreaks as well as for early detection of mosquito invasion events. Surveys are usually performed with traditional sampling tools (i.e., ovitraps and dipping method for immature stages or light or decoy traps for adults). Over the past decade, numerous studies have highlighted that environmental DNA (eDNA) sampling can enhance invertebrate species detection and provide community composition metrics. However, the usefulness of eDNA for detection of mosquito species has, to date, been largely neglected. Here, we sampled water from potential larval breeding sites along a gradient of anthropogenic perturbations, from the core of an oil palm plantation to the rainforest on São Tomé Island (Gulf of Guinea, Africa). We showed that (i) species of mosquitoes could be detected via metabarcoding mostly when larvae were visible, (ii) larvae species richness was greater using eDNA than visual identification and (iii) new mosquito species were also detected by the eDNA approach. We provide a critical discussion of the pros and cons of eDNA metabarcoding for monitoring mosquito species diversity and recommendations for future research directions that could facilitate the adoption of eDNA as a tool for assessing insect vector communities.

Statistical methods to identify mechanisms in studies of eco-evolutionary dynamics

While the reciprocal effects of ecological and evolutionary dynamics are increasingly recognized as an important driver for biodiversity, detection of such eco-evolutionary feedbacks, their underlying mechanisms, and their consequences remains challenging. Eco-evolutionary dynamics occur at different spatial and temporal scales and can leave signatures at different levels of organization (e.g., gene, protein, trait, community) that are often difficult to detect. Recent advances in statistical methods combined with alternative hypothesis testing provides a promising approach to identify potential eco-evolutionary drivers for observed data even in non-model systems that are not amenable to experimental manipulation. We discuss recent advances in eco-evolutionary modeling and statistical methods and discuss challenges for fitting mechanistic models to eco-evolutionary data.

Pooling of intra-site measurements inflates variability of the correlation between environmental DNA concentration and organism abundance

Environmental DNA (eDNA) analysis can promote efficient ecosystem monitoring and resource management. However, limited knowledge of the factors affecting the relationship between eDNA concentration and organism abundance causes uncertainty in relative abundance estimates based on eDNA concentration. Pooling of data points obtained from multiple locations within a site has been used to mitigate intra-site variation in eDNA and abundance estimates, but decreases the sample size used for estimating the relationship. I here assessed how the pooling of intra-site measurements of eDNA concentration and organism abundance impacted the reliability of the correlative relationship between eDNA concentration and organism abundance. Mathematical models were developed to simulate measurements of eDNA concentrations and organism abundances from multiple locations in a given survey site, and the CVs (coefficient of variability) of the correlations were compared depending on whether data points from different locations were individually treated or pooled. Although the mean and median values of the correlation coefficients were similar between the scenarios, the CVs of the simulated correlations were substantially higher under the pooled scenario than the individual scenario. Additionally, I re-analyzed two empirical studies conducted in lakes, both showing higher CVs of the correlations by pooling intra-site measurements. This study suggests that it would make eDNA-based abundance estimation more reliable and reproducible to individually analyze target eDNA concentrations and organism abundance estimates.

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.

Current Status of Omics in Biological Quality Elements for Freshwater Biomonitoring

Freshwater ecosystems have been experiencing various forms of threats, mainly since the last century. The severity of this adverse scenario presents unprecedented challenges to human health, water supply, agriculture, forestry, ecological systems, and biodiversity, among other areas. Despite the progress made in various biomonitoring techniques tailored to specific countries and biotic communities, significant constraints exist, particularly in assessing and quantifying biodiversity and its interplay with detrimental factors. Incorporating modern techniques into biomonitoring methodologies presents a challenging topic with multiple perspectives and assertions. This review aims to present a comprehensive overview of the contemporary advancements in freshwater biomonitoring, specifically by utilizing omics methodologies such as genomics, metagenomics, transcriptomics, proteomics, metabolomics, and multi-omics. The present study aims to elucidate the rationale behind the imperative need for modernization in this field. This will be achieved by presenting case studies, examining the diverse range of organisms that have been studied, and evaluating the potential benefits and drawbacks associated with the utilization of these methodologies. The utilization of advanced high-throughput bioinformatics techniques represents a sophisticated approach that necessitates a significant departure from the conventional practices of contemporary freshwater biomonitoring. The significant contributions of omics techniques in the context of biological quality elements (BQEs) and their interpretations in ecological problems are crucial for biomonitoring programs. Such contributions are primarily attributed to the previously overlooked identification of interactions between different levels of biological organization and their responses, isolated and combined, to specific critical conditions.

A comparison of biomonitoring methodologies for surf zone fish communities

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

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

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

Metabarcoding of fecal pellets in wild muskox populations reveals negative relationships between microbiome and diet alpha diversity

Microbiome diversity and diet composition concomitantly influence species health, fitness, immunity, and digestion. In environments where diet varies spatially and temporally, microbiome plasticity may promote rapid host adaptation to available resources. For northern ungulates in particular, metabarcoding of noninvasively collected fecal pellets presents unprecedented insights into their diverse ecological requirements and niches by clarifying the interrelationships of microbiomes, key to deriving nutrients, in context of altered forage availability in changing climates. Muskoxen (Ovibos moschatus) are Arctic-adapted species that experience fluctuating qualities and quantities of vegetation. Geography and seasonality have been noted to influence microbiome composition and diversity in muskoxen, yet it is unclear how their microbiomes intersect with diet. Following observations from other species, we hypothesized increasing diet diversity would result in higher microbiome diversity in muskoxen. We assessed diet composition in muskoxen using three common plant metabarcoding markers and explored correlations with microbiome data. Patterns of dietary diversity and composition were not fully concordant among the markers used, yet all reflected the primary consumption of willows and sedges. Individuals with similar diets had more similar microbiomes, yet in contrast to most literature, yielded negative relationships between microbiome and diet alpha diversity. This negative correlation may reflect the unique capacities of muskoxen to survive solely on high-fiber Arctic forage and provide insight into their resiliency to exploit changing dietary resources in a rapidly warming Arctic altering vegetation diversity.

Utilizing the state of environmental DNA (eDNA) to incorporate time-scale information into eDNA analysis

Environmental DNA (eDNA) analysis allows cost-effective and non-destructive biomonitoring with a high detection sensitivity in terrestrial and aquatic environments. However, the eDNA results can sometimes include false-positive inferences of target organisms owing to the detection of aged eDNA that has long since been released from the individual and is more likely to be detected at a site further away from its source. In order to address the issue, this manuscript focuses on the state of eDNA, proposing new methodologies to estimate the age of eDNA: (1) DNA damage rate, (2) eDNA particle size distribution, and (3) viable cell-derived eDNA. In addition, the manuscript also focuses on the shorter persistence of environmental RNA (eRNA) compared with eDNA, highlighting the application of eRNA and environmental nucleic acid ratio for assessing the age of the genetic materials in water. Although substantial further research is essential to support the feasibility of these methodologies, incorporating time-scale information into eDNA analysis would update current eDNA analysis, improve the accuracy and reliability of eDNA-based monitoring, and further refine eDNA analysis as a useful monitoring tool in ecology, fisheries and various environmental sciences.

Modelling environmental DNA transport in rivers reveals highly resolved spatio-temporal biodiversity patterns

The ever-increasing threats to riverine ecosystems call for novel approaches for highly resolved biodiversity assessments across taxonomic groups and spatio-temporal scales. Recent advances in the joint use of environmental DNA (eDNA) data and eDNA transport models in rivers (e.g., eDITH) allow uncovering the full structure of riverine biodiversity, hence elucidating ecosystem processes and supporting conservation measures. We applied eDITH to a metabarcoding dataset covering three taxonomic groups (fish, invertebrates, bacteria) and three seasons for a catchment sampled for eDNA at 73 sites. We upscaled eDNA-based biodiversity predictions to approximately 1900 reaches, and assessed α- and β-diversity patterns across seasons and taxonomic groups over the whole network. Genus richness predicted by eDITH was generally higher than values from direct eDNA analysis. Both predicted α- and β-diversity varied depending on season and taxonomic group. Predicted fish α-diversity increased downstream in all seasons, while invertebrate and bacteria α-diversity either decreased downstream or were unrelated to network position. Spatial β-diversity mostly decreased downstream, especially for bacteria. The eDITH model yielded a more refined assessment of freshwater biodiversity as compared to raw eDNA data, both in terms of spatial coverage, diversity patterns and effect of covariates, thus providing a more complete picture of freshwater biodiversity.

Serratia marcescens Colonization in a Neonatal Intensive Care Unit Has Multiple Sources, with Sink Drains as a Major Reservoir

Compelling evidence suggests a contribution of the sink environment to the transmission of opportunistic pathogens from the hospital environment to patients in neonatal intensive care units (NICU). In this study, the distribution of the opportunistic pathogen Serratia marcescens in the sink environment and newborns in a NICU was investigated. More than 500 sink drain and faucet samples were collected over the course of five sampling campaigns undertaken over 3 years. Distribution and diversity of S. marcescens were examined with a modified MacConkey medium and a high-throughput short-sequence typing (HiSST) method. Sink drains were an important reservoir of S. marcescens, with an average of 44% positive samples, whereas no faucet sample was positive. The genotypic diversity of S. marcescens was moderate, with an average of two genotypes per drain, while the spatial distribution of S. marcescens was heterogeneous. The genotypic profiles of 52 clinical isolates were highly heterogeneous, with 27 unique genotypes, of which 71% of isolates were found in more than one patient. S. marcescens acquisition during the first outbreaks was mainly caused by horizontal transmissions. HiSST analyses revealed 10 potential cases of patient-to-patient transmission of S. marcescens, five cases of patient-to-sink transmission, and one bidirectional transfer between sink and patient. Environmental and clinical isolates were found in sink drains up to 1 year after the first detection, supporting persisting drain colonization. This extensive survey suggests multiple reservoirs of S. marcescens within the NICU, including patients and sink drains, but other external sources should also be considered. IMPORTANCE The bacterium Serratia marcescens is an important opportunistic human pathogen that thrives in many environments, can become multidrug resistant, and is often involved in nosocomial outbreaks in neonatal intensive care units (NICU). We evaluated the role of sinks during five suspected S. marcescens outbreaks in a NICU. An innovative approach combining molecular and culture methods was used to maximize the detection and typing of S. marcescens in the sink environment. Our results indicate multiple reservoirs of S. marcescens within the NICU, including patients, sink drains, and external sources. These results highlight the importance of sinks as a major reservoir of S. marcescens and potential sources of future outbreaks.

Heterogeneity of Fish Taxonomic and Functional Diversity Evaluated by eDNA and Gillnet along a Mangrove-Seagrass-Coral Reef Continuum

The effective and reliable monitoring of fish communities is important for the management and protection of marine ecosystems. Environmental DNA (eDNA) metabarcoding is a relatively new method that has been widely used in recent years, while traditional sampling via fish catching (i.e., gillnets) is one of the most common and reliable fish monitoring methods used to date. We compared the taxonomic and functional diversity of fish detected within a mangrove-seagrass-coral reef continuum using both survey methods. One liter seawater and gillnet samples were collected in August 2021 from mangrove forests, seagrass meadows and coral reef habitats (n = 3 each) in Hainan, China. Surveys using eDNA and gillnets identified 139 genera belonging to 66 families and 58 genera belonging to 42 families, respectively. Regardless of the survey method, fish detected in mangrove, seagrass and coral reef habitats were heterogeneous in their communities; however, the shared species between habitats suggest some degree of connectivity. There were no significant differences between habitats in terms of taxonomic and functional diversity, but a higher taxonomic diversity was detected using eDNA. Both methods were able to distinguish fish assemblages between different habitats; however, gillnet surveys performed better than eDNA surveys for distinguishing mangrove from seagrass assemblages. Therefore, the concurrent use of eDNA and gillnet survey methods provides a more comprehensive approach to understanding the heterogeneity of fish taxonomic and functional diversity along mangrove-seagrass-coral reef continuums.