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

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

Background

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

Findings

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

Conclusion

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

Demystifying eDNA validation

As environmental DNA (eDNA) approaches gain momentum for biodiversity analysis, validation becomes a key consideration. I focus on four facets of eDNA validation. Validation through technical processes, legal use, official statements, and 'good enough' scenarios can advance the field to aid societal issues such as climate emergency and biodiversity crisis.

The potential of metabarcoding plant components of Malaise trap samples to enhance knowledge of plant-insect interactions

The worldwide rapid declines in insect and plant abundance and diversity that have occurred in the past decades have gained public attention and demand for political actions to counteract these declines are growing. Rapid large-scale biomonitoring can aid in observing these changes and provide information for decisions for land management and species protection. Malaise traps have long been used for insect sampling and when insects are captured in these traps, they carry traces of plants they have visited on the body surface or as digested food material in the gut contents. Metabarcoding offers a promising method for identifying these plant traces, providing insight into the plants with which insects are directly interacting at a given time. To test the efficacy of DNA metabarcoding with these sample types, 79 samples from 21 sites across Germany were analysed with the ITS2 barcode. This study, to our knowledge, is the first examination of metabarcoding plant DNA traces from Malaise trap samples. Here, we report on the feasibility of sequencing these sample types, analysis of the resulting taxa, the usage of cultivated plants by insects near nature conservancy areas and the detection of rare and neophyte species. Due to the frequency of contamination and false positive reads, isolation and PCR negative controls should be used in every reaction. Metabarcoding has advantages in efficiency and resolution over microscopic identification of pollen and is the only possible identification method for the other plant traces from Malaise traps and could provide a broad utility for future studies of plant-insect interactions.

Constructing ecological indices for urban environments using species distribution models

In an increasingly urbanized world, there is a need to study urban areas as their own class of ecosystems as well as assess the impacts of anthropogenic impacts on biodiversity. However, collecting a sufficient number of species observations to estimate patterns of biodiversity in a city can be costly. Here we investigated the use of community science-based data on species occurrences, combined with species distribution models (SDMs), built using MaxEnt and remotely-sensed measures of the environment, to predict the distribution of a number of species across the urban environment of Los Angeles. By selecting species with the most accurate SDMs, and then summarizing these by class, we were able to produce two species richness models (SRMs) to predict biodiversity patterns for species in the class Aves and Magnoliopsida and how they respond to a variety of natural and anthropogenic environmental gradients.

We found that species considered native to Los Angeles tend to have significantly more accurate SDMs than their non-native counterparts. For all species considered in this study we found environmental variables describing anthropogenic activities, such as housing density and alterations to land cover, tend to be more influential than natural factors, such as terrain and proximity to freshwater, in shaping SDMs. Using a random forest model we found our SRMs could account for approximately 54% and 62% of the predicted variation in species richness for species in the classes Aves and Magnoliopsida respectively. Using community science-based species occurrences, SRMs can be used to model patterns of urban biodiversity and assess the roles of environmental factors in shaping them.

Curation of a reference database of COI sequences for insect identification through DNA metabarcoding: COins

DNA metabarcoding is a widespread approach for the molecular identification of organisms. While the associated wet-lab and data processing procedures are well established and highly efficient, the reference databases for taxonomic assignment can be implemented to improve the accuracy of identifications. Insects are among the organisms for which DNA-based identification is most commonly used; yet, a DNA-metabarcoding reference database specifically curated for their species identification using software requiring local databases is lacking. Here, we present COins, a database of 5' region cytochrome c oxidase subunit I sequences (COI-5P) of insects that includes over 532 000 representative sequences of >106 000 species specifically formatted for the QIIME2 software platform. Through a combination of automated and manually curated steps, we developed this database starting from all COI sequences available in the Barcode of Life Data System for insects, focusing on sequences that comply with several standards, including a species-level identification. COins was validated on previously published DNA-metabarcoding sequences data (bulk samples from Malaise traps) and its efficiency compared with other publicly available reference databases (not specific for insects). COins can allow an increase of up to 30% of species-level identifications and thus can represent a valuable resource for the taxonomic assignment of insects' DNA-metabarcoding data, especially when species-level identification is needed https://doi.org/10.6084/m9.figshare.19130465.v1.

Similarities and Differences in Fish Community Composition Accessed by Electrofishing, Gill Netting, Seining, Trawling, and Water eDNA Metabarcoding in Temperate Reservoirs

It is difficult to understand the composition and diversity of biological communities in complex and heterogeneous environments using traditional sampling methods. Recently, developments in environmental DNA metabarcoding have emerged as a powerful, non-invasive method for comprehensive community characterization and biodiversity monitoring in different types of aquatic ecosystems. In this study, water eDNA targeting fish (wf-eDNA) and four traditional fish sampling methods (electrofishing, gill netting, seining, trawling) were compared to evaluate the reliability and efficiency of wf-eDNA (vertebrate mitochondrial 12S ribosomal RNA (rRNA) as an alternative approach to assess the diversity and composition of freshwater fish communities. The results of wf-eDNA showed a consistency between the traditional sampling methods regarding species detection. However, some fish species detected using wf-eDNA assay were not detected using traditional sampling methods and vice versa. Comparison of wf-eDNA and traditional sampling methods revealed spatial homogeneity in fish community composition in all reservoirs. Ordination analysis showed that the wf-eDNA approach covers all traditional sampling methods and occupies an intermediate position. In addition, based on the Shannon diversity index, we found that in one reservoir the wf-eDNA method yielded similar fish community diversity to traditional sampling methods. However, in other reservoirs, the calculated Shannon diversity index of the wf-eDNA method was significantly higher than traditional sampling methods. In general, significant positive correlations were found between the wf-eDNA method and almost all traditional sampling methods. We conclude that wf-eDNA seems to be a reliable and complementary approach for biomonitoring and ecosystem management of freshwater ichthyofauna.

eDNAssay: a machine learning tool that accurately predicts qPCR cross-amplification

Environmental DNA (eDNA) sampling is a highly sensitive and cost-effective technique for wildlife monitoring, notably through the use of qPCR assays. However, it can be difficult to ensure assay specificity when many closely related species cooccur. In theory, specificity may be assessed in silico by determining whether assay oligonucleotides have enough base-pair mismatches with nontarget sequences to preclude amplification. However, the mismatch qualities required are poorly understood, making in silico assessments difficult and often necessitating extensive in vitro testing-typically the greatest bottleneck in assay development. Increasing the accuracy of in silico assessments would therefore streamline the assay development process. In this study, we paired 10 qPCR assays with 82 synthetic gene fragments for 530 specificity tests using SYBR Green intercalating dye (n = 262) and TaqMan hydrolysis probes (n = 268). Test results were used to train random forest classifiers to predict amplification. The primer-only model (SYBR Green-based) and full-assay model (TaqMan probe-based) were 99.6% and 100% accurate, respectively, in cross-validation. We further assessed model performance using six independent assays not used in model training. In these tests the primer-only model was 92.4% accurate (n = 119) and the full-assay model was 96.5% accurate (n = 144). The high performance achieved by these models makes it possible for eDNA practitioners to more quickly and confidently develop assays specific to the intended target. Practitioners can access the full-assay model via eDNAssay (https://NationalGenomicsCenter.shinyapps.io/eDNAssay), a user-friendly online tool for predicting qPCR cross-amplification.

Applying convolutional neural networks to speed up environmental DNA annotation in a highly diverse ecosystem

High-throughput DNA sequencing is becoming an increasingly important tool to monitor and better understand biodiversity responses to environmental changes in a standardized and reproducible way. Environmental DNA (eDNA) from organisms can be captured in ecosystem samples and sequenced using metabarcoding, but processing large volumes of eDNA data and annotating sequences to recognized taxa remains computationally expensive. Speed and accuracy are two major bottlenecks in this critical step. Here, we evaluated the ability of convolutional neural networks (CNNs) to process short eDNA sequences and associate them with taxonomic labels. Using a unique eDNA data set collected in highly diverse Tropical South America, we compared the speed and accuracy of CNNs with that of a well-known bioinformatic pipeline (OBITools) in processing a small region (60 bp) of the 12S ribosomal DNA targeting freshwater fishes. We found that the taxonomic labels from the CNNs were comparable to those from OBITools, with high correlation levels for the composition of the regional fish fauna. The CNNs enabled the processing of raw fastq files at a rate of approximately 1 million sequences per minute, which was about 150 times faster than with OBITools. Given the good performance of CNNs in the highly diverse ecosystem considered here, the development of more elaborate CNNs promises fast deployment for future biodiversity inventories using eDNA.

Sensitive and accurate DNA metabarcoding of parasitic helminth mock communities using the mitochondrial rRNA genes

Next-generation sequencing technologies have accelerated the pace of helminth DNA metabarcoding research, enabling species detection in bulk community samples. However, finding suitable genetic markers with robust species-level resolution and primers targeting a broad species range among parasitic helminths are some of the challenges faced. This study aimed to demonstrate the potential use of the mitochondrial 12S and 16S rRNA genes for parasitic helminth (nematodes, trematodes, cestodes) DNA metabarcoding. To demonstrate the robustness of the 12S and 16S rRNA genes for DNA metabarcoding, we determined the proportion of species successfully recovered using mock helminth communities without environment matrix and mock helminth communities artificially spiked with environmental matrices. The environmental matrices are human fecal material, garden soil, tissue, and pond water. Our results revealed the robustness of the mitochondrial rRNA genes, through the high sensitivity of the 12S rRNA gene, and the effectiveness of the 12S and 16S primers targeting platyhelminths. With the mitochondrial rRNA genes, a broad range of parasitc helminths were successfully detected to the species level. The potential of the mitochondrial rRNA genes for helminth DNA metabarcoding was demonstrated, providing a valuable gateway for future helminth DNA metabarcoding applications like helminth detection and biodiversity studies.

DNA metabarcoding across disciplines: sequencing our way to greater understanding across scales of biological organization

DNA metabarcoding describes the use of targeted DNA (i.e., amplicon) sequencing to identify community constituents from a complex sample containing genetic material from multiple organisms, such as water, soil, gut contents, microbiomes, or biofilms. This molecular approach for characterizing mixed DNA samples relies on the development of "universal primers" that allow for effective amplification of target sequences across a broad range of taxa. Armed with optimized lab protocols and rigorous bioinformatics tools, DNA metabarcoding can produce a wealth of information about the hidden biodiversity of various sample types by probing for organisms' molecular footprints. DNA metabarcoding has received considerable popular press over the last few years because of gut microbiome studies in humans and beyond. However, there are many other applications that are continually integrating molecular biology with other fields of study to address questions that have previously been unanswerable, for both prokaryotic and eukaryotic targets. For example, we can now sample mostly-digested gut contents from virtually any organism to learn about ontogeny and foraging ecology. Water samples collected from different locations can be filtered to extract eDNA (i.e., environmental DNA), revealing the biodiversity of fishes and other taxa targeted by carefully selected primer sets. This universal primer metabarcoding approach has even been extended to looking at diverse gene families within single species, which is particularly useful for complex immune system genetics. The purpose of this SICB symposium was to bring together researchers using DNA metabarcoding approaches to (a) showcase the diversity of applications of this technique for addressing questions spanning ecology, evolution, and physiology, and (b) to spark connections among investigators from different fields that are utilizing similar approaches to facilitate optimization and standardization of metabarcoding methods and analyses. The resulting manuscripts from this symposium represent a great diversity of metabarcoding applications and taxonomic groups of interest.

Tracking an invasion front with environmental DNA

Data from environmental DNA (eDNA) may revolutionize environmental monitoring and management, providing increased detection sensitivity at reduced cost and survey effort. However, eDNA data are rarely used in decision-making contexts, mainly due to uncertainty around (1) data interpretation and (2) whether and how molecular tools dovetail with existing management efforts. We address these challenges by jointly modeling eDNA detection via qPCR and traditional trap data to estimate the density of invasive European green crab (Carcinus maenas), a species for which, historically, baited traps have been used for both detection and control. Our analytical framework simultaneously quantifies uncertainty in both detection methods and provides a robust way of integrating different data streams into management processes. Moreover, the joint model makes clear the marginal information benefit of adding eDNA (or any other) additional data type to an existing monitoring program, offering a path to optimizing sampling efforts for species of management interest. Here, we document green crab eDNA beyond the previously known invasion front and find that the value of eDNA data dramatically increases with low population densities and low traditional sampling effort, as is often the case at leading-edge locations. We also highlight the detection limits of the molecular assay used in this study, as well as scenarios under which eDNA sampling is unlikely to improve existing management efforts.

Metagenomic Analysis Using Phylogenetic Placement-A Review of the First Decade

Phylogenetic placement refers to a family of tools and methods to analyze, visualize, and interpret the tsunami of metagenomic sequencing data generated by high-throughput sequencing. Compared to alternative (e. g., similarity-based) methods, it puts metabarcoding sequences into a phylogenetic context using a set of known reference sequences and taking evolutionary history into account. Thereby, one can increase the accuracy of metagenomic surveys and eliminate the requirement for having exact or close matches with existing sequence databases. Phylogenetic placement constitutes a valuable analysis tool per se, but also entails a plethora of downstream tools to interpret its results. A common use case is to analyze species communities obtained from metagenomic sequencing, for example via taxonomic assignment, diversity quantification, sample comparison, and identification of correlations with environmental variables. In this review, we provide an overview over the methods developed during the first 10 years. In particular, the goals of this review are 1) to motivate the usage of phylogenetic placement and illustrate some of its use cases, 2) to outline the full workflow, from raw sequences to publishable figures, including best practices, 3) to introduce the most common tools and methods and their capabilities, 4) to point out common placement pitfalls and misconceptions, 5) to showcase typical placement-based analyses, and how they can help to analyze, visualize, and interpret phylogenetic placement data.

The use of citizen science in fish eDNA metabarcoding for evaluating regional biodiversity in a coastal marine region: A pilot study

To test the feasibility of a citizen science program for fish eDNA metabarcoding in coastal marine environments, we recruited six groups of voluntary citizens for a science education course at a natural history museum. We held a seminar on eDNA and a workshop for seawater sampling and on-site filtration using syringes and filter cartridges for the participants. After that, they selected single survey sites following the guidelines for conducting a safe field trip. They performed seawater sampling and on-site filtration at these sites during their summer holidays. The six selected sites unexpectedly included diverse coastal habitats within a 40 km radius, located at temperate latitudes in central Japan (~35°N). After the field trips, they returned filtered cartridges to the museum, and we extracted eDNA from the filters. We performed fish eDNA metabarcoding, along with data analysis. Consequently, we identified 140 fish species across 66 families and 118 genera from the six samples, with species richness ranging from 14 to 66. Despite its limited sample size, such a diverse taxonomic range of fish species exhibited spatial biodiversity patterns within the region, which are consistent with species distribution. These include north-south and urbanization gradients of species richness, geographic structure of the fish communities, and varying salinity preferences of the component species. This case study demonstrates the potential of fish eDNA metabarcoding as an educational and scientific tool to raise public awareness and perform large-scale citizen science initiatives encompassing regional, national, or global fauna.

Bee foraging preferences, microbiota and pathogens revealed by direct shotgun metagenomics of honey

Honeybees (Apis mellifera) continue to succumb to human and environmental pressures despite their crucial role in providing essential ecosystem services. Owing to their foraging and honey production activities, honeybees form complex relationships with species across all domains, such as plants, viruses, bacteria and other hive pests, making honey a valuable biomonitoring tool for assessing their ecological niche. Thus, the application of honey shotgun metagenomics (SM) has paved the way for a detailed description of the species honeybees interact with. Nevertheless, SM bioinformatics tools and DNA extraction methods rely on resources not necessarily optimized for honey. In this study, we compared five widely used taxonomic classifiers using simulated species communities commonly found in honey. We found that Kraken 2 with a threshold of 0.5 performs best in assessing species distribution. We also optimized a simple NaOH-based honey DNA extraction methodology (Direct-SM), which profiled species seasonal variability similarly to an established column-based DNA extraction approach (SM). Both approaches produce results consistent with melissopalinology analysis describing the botanical landscape surrounding the apiary. Interestingly, we detected a strong stability of the bacteria constituting the core and noncore gut microbiome across seasons, pointing to the potential utility of honey for noninvasive assessment of bee microbiota. Finally, the Direct-SM approach to detect Varroa correlates well with the biomonitoring of mite infestation observed in hives. These observations suggest that Direct-SM methodology has the potential to comprehensively describe honeybee ecological niches and can be tested as a building block for large-scale studies to assess bee health in the field.

Diversity, distribution and ecology of fungal communities present in Antarctic lake sediments uncovered by DNA metabarcoding

We assessed fungal diversity in sediments obtained from four lakes in the South Shetland Islands and James Ross Island, Antarctica, using DNA metabarcoding. We detected 218 amplicon sequence variants (ASVs) dominated by the phyla Ascomycota, Basidiomycota, Mortierellomycota, Mucoromycota and Chytridiomycota. In addition, the rare phyla Aphelidiomycota, Basidiobolomycota, Blastocladiomycota, Monoblepharomycota, Rozellomycota and Zoopagomycota as well as fungal-like Straminopila belonging to the phyla Bacillariophyta and Oomycota were detected. The fungal assemblages were dominated by unknown fungal taxa (Fungal sp. 1 and Fungal sp. 2), followed by Talaromyces rubicundus and Dactylonectria anthuriicola. In general, they displayed high diversity, richness and moderate dominance. Sequences representing saprophytic, pathogenic and symbiotic fungi were detected, including the phytopathogenic fungus D. anthuriicola that was abundant, in the relatively young Soto Lake on Deception Island. The lake sediments studied contained the DNA of rich, diverse and complex fungal communities, including both fungi commonly reported in Antarctica and other taxa considered to be rare. However, as the study was based on the use of environmental DNA, which does not unequivocally confirm the presence of active or viable organisms, further studies using other approaches such as shotgun sequencing are required to elucidate the ecology of fungi in these Antarctic lake sediments.

Comparison of materials for rapid passive collection of environmental DNA

Passive collection is an emerging sampling method for environmental DNA (eDNA) in aquatic systems. Passive eDNA collection is inexpensive and efficient, and requires minimal equipment, making it suited to high‐density sampling and remote deployment. Here, we compare the effectiveness of nine membrane materials for passively collecting fish eDNA from a 3‐million‐litre marine mesocosm. We submerged materials (cellulose, cellulose with 1% and 3% chitosan, cellulose overlayed with electrospun nanofibres and 1% chitosan, cotton fibres, hemp fibres, and sponge with either zeolite or active carbon) for intervals between 5 and 1080 min. We show that for most materials, with as little as 5 min of submersion, mitochondrial fish eDNA measured with qPCR, and fish species richness measured with metabarcoding, was comparable to that collected by conventional filtering. Furthermore, PCR template DNA concentrations and species richness were generally not improved significantly by longer submersion. Species richness detected for all materials ranged between 11 and 37 species, with a median of 27, which was comparable to the range for filtered eDNA (19–32). Using scanning electron microscopy, we visualized biological matter adhering to the surface of materials, rather than entrapped, with images also revealing a diversity in size and structure of putative eDNA particles. eDNA can be collected rapidly from seawater with a passive approach and using a variety of materials. This will suit cost‐ and time‐sensitive biological surveys, and where access to equipment is limited.

Similar fish species composition despite larger environmental heterogeneity during severe hypoxia in a coastal ecosystem

Environmental heterogeneity is one of the most influential factors that create compositional variation among local communities. Greater compositional variation is expected when an environmental gradient encompasses the most severe conditions where species sorting is more likely to operate. However, evidence for stronger species sorting at severer environment has typically been obtained for less mobile organisms and tests are scarce for those with higher dispersal ability that allows individuals to sensitively respond to environmental stress. Here, with the dynamics of fish communities in a Japanese bay revealed by environmental DNA metabarcoding analyses as a model case, we tested the hypothesis that larger environmental heterogeneity caused by severe seasonal hypoxia (lower concentration of oxygen in bottom waters in summer) leads to larger variation of species composition among communities. During summer, fish species richness was lower in the bottom layer, suggesting the severity of the hypoxic bottom water. In contrast to the prediction, we found that although the environmental parameters of bottom and surface water was clearly distinct in summer, fish species composition was more similar between the two layers. Our null model analysis suggested that the higher compositional similarity during hypoxia season was not a result of the sampling effect reflecting differences in the alpha or gamma diversity. Furthermore, a shift in the species occurrence from bottom to surface layers was observed during hypoxia season, which was consistent across species, suggesting that the severe condition in the bottom adversely affected fish species irrespective of their identity. These results suggest that larger environmental heterogeneity does not necessarily lead to higher compositional variation once the environmental gradient encompasses extremely severe conditions. This is most likely because individual organisms actively avoided the severity quasi‐neutrally, which induced mass effect‐like dispersal and lead to the mixing of species composition across habitats. By showing counter evidence against the prevailing view, we provide novel insights into how species sorting by environment acts in heterogeneous and severe conditions.

met v1: expanding on old estimations of biodiversity from eDNA with a new database framework

A long-standing problem in environmental DNA has been the inability to compute across large number of datasets. Here we introduce an open-source software framework that can store a large number of environmental DNA datasets, as well as provide a platform for analysis, in an easily customizable way. We show the utility of such an approach by analyzing over 1400 arthropod metabarcode datasets. This article introduces a new software framework, met, which utilizes large numbers of metabarcode datasets to draw conclusions about patterns of diversity at large spatial scales. Given more accurate estimations on the distribution of variance in metabarcode datasets, this software framework could facilitate novel analyses that are outside the scope of currently available similar platforms.

Database URL https://osf.io/spb8v/

eDNA in subterranean ecosystems: Applications, technical aspects, and future prospects

Monitoring of biota is pivotal for the assessment and conservation of ecosystems. Environments worldwide are being continuously and increasingly exposed to multiple adverse impacts, and the accuracy and reliability of the biomonitoring tools that can be employed shape not only the present, but more importantly, the future of entire habitats. The analysis of environmental DNA (eDNA) metabarcoding data provides a quick, affordable, and reliable molecular approach for biodiversity assessments. However, while extensively employed in aquatic and terrestrial surface environments, eDNA-based studies targeting subterranean ecosystems are still uncommon due to the lack of accessibility and the cryptic nature of these environments and their species. Recent advances in genetic and genomic analyses have established a promising framework for shedding new light on subterranean biodiversity and ecology. To address current knowledge and the future use of eDNA methods in groundwaters and caves, this review explores conceptual and technical aspects of the application and its potential in subterranean systems. We briefly introduce subterranean biota and describe the most used traditional sampling techniques. Next, eDNA characteristics, application, and limitations in the subsurface environment are outlined. Last, we provide suggestions on how to overcome caveats and delineate some of the research avenues that will likely shape this field in the near future. We advocate that eDNA analyses, when carefully conducted and ideally combined with conventional sampling techniques, will substantially increase understanding and enable crucial expansion of subterranean community characterisation. Given the importance of groundwater and cave ecosystems for nature and humans, eDNA can bring to the surface essential insights, such as study of ecosystem assemblages and rare species detection, which are critical for the preservation of life below, as well as above, the ground.

Temporal Patterns of Honeybee Foraging in a Diverse Floral Landscape Revealed Using Pollen DNA Metabarcoding of Honey

Understanding the plants pollinators use through the year is vital to support pollinator populations and mitigate for declines in floral resources due to habitat loss. DNA metabarcoding allows the temporal picture of nectar and pollen foraging to be examined in detail. Here, we use DNA metabarcoding to examine the forage use of honeybees (Apis mellifera L.) within a florally diverse landscape within the UK, documenting the key forage plants used and seasonal progression over two years. The total number of plant taxa detected in the honey was 120, but only 16 of these were found with a high relative read abundance of DNA, across the main foraging months (April-September). Only a small proportion of the available flowering genera in the landscape were used by the honeybees. The greatest relative read abundance came from native or near-native plants, including Rubus spp., Trifolium repens, the Maleae tribe including Crataegus, Malus and Cotoneaster, and Hedera helix. Tree species were important forage in the spring months, followed by increased use of herbs and shrubs later in the foraging season. Garden habitat increased the taxon richness of native, near-native and horticultural plants found in the honey. Although horticultural plants were rarely found abundantly within the honey samples, they may be important for increasing nutritional diversity of the pollen forage.

A triad of kicknet sampling, eDNA metabarcoding, and predictive modeling to assess richness of mayflies, stoneflies and caddisflies in rivers

Monitoring biodiversity is essential to understand the impacts of human activities and for effective management of ecosystems. Thereby, biodiversity can be assessed through direct collection of targeted organisms, through indirect evidence of their presence (e.g. signs, environmental DNA, camera trap, etc.), or through extrapolations from species distribution and species richness models. Differences in approaches used in biodiversity assessment, however, may come with individual challenges and hinder cross-study comparability. In the context of rapidly developing techniques, we compared three different approaches in order to better understand assessments of aquatic macroinvertebrate diversity. Specifically, we compared the community composition and species richness of three orders of aquatic macroinvertebrates (mayflies, stoneflies, and caddisflies, hereafter EPT) obtained via eDNA metabarcoding and via traditional in situ kicknet sampling to catchment-level based predictions of a species richness model. We used kicknet data from 24 sites in Switzerland and compared taxonomic lists to those obtained using eDNA amplified with two different primer sets. Richness detected by these methods was compared to the independent predictions made by a statistical species richness model, that is, a generalized linear model using landscape-level features to estimate EPT diversity. Despite the ability of eDNA to consistently detect some EPT species found by traditional sampling, we found important discrepancies in community composition between the kicknet and eDNA approaches, particularly at a local scale. We found the EPT-specific primer set fwhF2/EPTDr2n, detected a greater number of targeted EPT species compared to the more general primer set mlCOIintF/HCO2198. Moreover, we found that the species richness measured by eDNA from either primer set was poorly correlated to the richness measured by kicknet sampling (Pearson correlation = 0.27) and that the richness estimated by eDNA and kicknet were poorly correlated with the prediction of the species richness model (Pearson correlation = 0.30 and 0.44, respectively). The weak relationships between the traditional kicknet sampling and eDNA with this model indicates inherent limitations in upscaling species richness estimates, and possibly a limited ability of the model to meet real world expectations. It is also possible that the number of replicates was not sufficient to detect ambiguous correlations. Future challenges include improving the accuracy and sensitivity of each approach individually, yet also acknowledging their respective limitations, in order to best meet stakeholder demands and address the biodiversity crisis we are facing.

Spider webs as eDNA samplers: biodiversity assessment across the tree of life

The concept of environmental DNA (eDNA) utilizes nucleic acids of organisms directly from the environment. Recent breakthrough studies have successfully detected a wide spectrum of prokaryotic and eukaryotic eDNA from a variety of environments, ranging from ancient to modern, and from terrestrial to aquatic. With their diversity and ubiquity in nature, spider webs might act as powerful biofilters and could thus represent a promising new source of eDNA, but their utility under natural field conditions is severely understudied. Here, we bridge this knowledge gap to establish spider webs as a source of eDNA with far reaching implications. First, we conducted a field study to track specific arthropod targets from different spider webs. We then employed high-throughput amplicon sequencing of taxonomic barcodes to investigate the utility of spider web eDNA for biodiversity monitoring of animals, fungi, and bacteria. Our results show that genetic remains on spider webs allow the detection of even the smallest target organisms. We also demonstrate that eDNA from spider webs is useful in research of community compositions across the different domains of life, with potentially highly detailed temporal and spatial information.

The Multiple States of Environmental DNA and What Is Known about Their Persistence in Aquatic Environments

Increased use of environmental DNA (eDNA) analysis for indirect species detection has spurred the need to understand eDNA persistence in the environment. Understanding the persistence of eDNA is complex because it exists in a mixture of different states (e.g., dissolved, particle adsorbed, intracellular, and intraorganellar), and each state is expected to have a specific decay rate that depends on environmental parameters. Thus, improving knowledge about eDNA conversion rates between states and the reactions that degrade eDNA in different states is needed. Here, we focus on eukaryotic extraorganismal eDNA, outline how water chemistry and suspended mineral particles likely affect conversion among each eDNA state, and indicate how environmental parameters affect persistence of states in the water column. On the basis of deducing these controlling parameters, we synthesized the eDNA literature to assess whether we could already derive a general understanding of eDNA states persisting in the environment. However, we found that these parameters are often not being measured or reported when measured, and in many cases very few experimental data exist from which to draw conclusions. Therefore, further study of how environmental parameters affect eDNA state conversion and eDNA decay in aquatic environments is needed. We recommend analytic controls that can be used during the processing of water to assess potential losses of different eDNA states if all were present in a water sample, and we outline future experimental work that would help determine the dominant eDNA states in water.

Fungal diversity in a sediment core from climate change impacted Boeckella Lake, Hope Bay, north-eastern Antarctic Peninsula assessed using metabarcoding

We studied the fungal DNA present in a lake sediment core obtained from Trinity Peninsula, Hope Bay, north-eastern Antarctic Peninsula, using metabarcoding through high-throughput sequencing (HTS). Sequences obtained were assigned to 146 amplicon sequence variants (ASVs) primarily representing unknown fungi, followed by the phyla Ascomycota, Rozellomycota, Basidiomycota, Chytridiomycota and Mortierellomycota. The most abundant taxa were assigned to Fungal sp., Pseudeurotium hygrophilum, Rozellomycota sp. 1, Pseudeurotiaceae sp. 1 and Chytridiomycota sp. 1. The majority of the DNA reads, representing 40 ASVs, could only be assigned at higher taxonomic levels and may represent taxa not currently included in the sequence databases consulted and/or be previously undescribed fungi. Different sections of the core were characterized by high sequence diversity, richness and moderate ecological dominance indices. The assigned diversity was dominated by cosmopolitan cold-adapted fungi, including known saprotrophic, plant and animal pathogenic and symbiotic taxa. Despite the overall dominance of Ascomycota and Basidiomycota and psychrophilic Mortierellomycota, members of the cryptic phyla Rozellomycota and Chytridiomycota were also detected in abundance. As Boeckella Lake may cease to exist in approaching decades due the effects of local climatic changes, it also an important location for the study of the impacts of these changes on Antarctic microbial diversity.

Sampling multiple life stages significantly increases estimates of marine biodiversity

Biodiversity assessments are critical for setting conservation priorities, understanding ecosystem function and establishing a baseline to monitor change. Surveys of marine biodiversity that rely almost entirely on sampling adult organisms underestimate diversity because they tend to be limited to habitat types and individuals that can be easily surveyed. Many marine animals have planktonic larvae that can be sampled from the water column at shallow depths. This life stage often is overlooked in surveys but can be used to relatively rapidly document diversity, especially for the many species that are rare or live cryptically as adults. Using DNA barcode data from samples of nemertean worms collected in three biogeographical regions—Northeastern Pacific, the Caribbean Sea and Eastern Tropical Pacific—we found that most species were collected as either benthic adults or planktonic larvae but seldom in both stages. Randomization tests show that this deficit of operational taxonomic units collected as both adults and larvae is extremely unlikely if larvae and adults were drawn from the same pool of species. This effect persists even in well-studied faunas. These results suggest that sampling planktonic larvae offers access to a different subset of species and thus significantly increases estimates of biodiversity compared to sampling adults alone. Spanish abstract is available in the electronic supplementary material.

Sarracenia pitcher plant-associated microbial communities differ primarily by host species across a longitudinal gradient

Plant-associated microbial communities can profoundly affect plant health and success, and research is still uncovering factors driving the assembly of these communities. Here, we examine how geography versus host species affect microbial community structure and differential abundances of individual taxa. We use metabarcoding to characterize the bacteria and eukaryotes associated with five, often co-occurring species of Sarracenia pitcher plants (Sarraceniaceae) and three natural hybrids along the longitudinal gradient of the U.S. Gulf Coast, as well as samples from S. purpurea in Massachusetts. To tease apart the effects of geography versus host species, we focus first on sites with co-occurring species and then on species located across different sites. Our analyses show that bacterial and eukaryotic community structures are clearly and consistently influenced by host species identity, with geographic factors also playing a role. Naturally-occurring hybrids appear to also host unique communities, that are in some ways intermediate between their parent species. We see significant effects of geography (site and longitude), but these generally explain less of the variation among pitcher communities. Overall, in Sarracenia pitchers, host plant phenotype significantly affects the pitcher microbiomes and other associated organisms. This article is protected by copyright. All rights reserved.

Development and validation of an eDNA protocol for monitoring endemic Asian spiny frogs in the Himalayan region of Pakistan

Wildlife monitoring programs are instrumental for the assessment of species, habitat status, and for the management of factors affecting them. This is particularly important for species found in freshwater ecosystems, such as amphibians, as they have higher estimated extinction rates than terrestrial species. We developed and validated two species-specific environmental DNA (eDNA) protocols and applied them in the field to detect the Hazara Torrent Frog (Allopaa hazarensis) and Murree Hills Frog (Nanorana vicina). Additionally, we compared eDNA surveys with visual encounter surveys and estimated site occupancy. eDNA surveys resulted in higher occurrence probabilities for both A. hazarensis and N. vicina than for visual encounter surveys. Detection probability using eDNA was greater for both species, particularly for A. hazarensis. The top-ranked detection model for visual encounter surveys included effects of both year and temperature on both species, and the top-ranked occupancy model included effects of elevation and year. The top-ranked detection model for eDNA data was the null model, and the top-ranked occupancy model included effects of elevation, year, and wetland type. To our knowledge, this is the first time an eDNA survey has been used to monitor amphibian species in the Himalayan region.

Fungi in the Antarctic Cryosphere: Using DNA Metabarcoding to Reveal Fungal Diversity in Glacial Ice from the Antarctic Peninsula Region

We assessed fungal diversity present in glacial from the Antarctic Peninsula using DNA metabarcoding through high-throughput sequencing (HTS). We detected a total of 353,879 fungal DNA reads, representing 94 genera and 184 taxa, in glacial ice fragments obtained from seven sites in the north-west Antarctic Peninsula and South Shetland Islands. The phylum Ascomycota dominated the sequence diversity, followed by Basidiomycota and Mortierellomycota. Penicillium sp., Cladosporium sp., Penicillium atrovenetum, Epicoccum nigrum, Pseudogymnoascus sp. 1, Pseudogymnoascus sp. 2, Phaeosphaeriaceae sp. and Xylaria grammica were the most dominant taxa, respectively. However, the majority of the fungal diversity comprised taxa of rare and intermediate relative abundance, predominately known mesophilic fungi. High indices of diversity and richness were calculated, along with moderate index of dominance, which varied among the different sampling sites. Only 26 (14%) of the total fungal taxa detected were present at all sampling sites. The identified diversity was dominated by saprophytic taxa, followed by known plant and animal pathogens and a low number of symbiotic fungi. Our data suggest that Antarctic glacial ice may represent a hotspot of previously unreported fungal diversity; however, further studies are required to integrate HTS and culture approaches to confirm viability of the taxa detected.

Worms and bugs of the gut: the search for diagnostic signatures using barcoding, and metagenomics-metabolomics

Gastrointestinal (GI) helminth infections cause significant morbidity in both humans and animals worldwide. Specific and sensitive diagnosis is central to the surveillance of such infections and to determine the effectiveness of treatment strategies used to control them. In this article, we: (i) assess the strengths and limitations of existing methods applied to the diagnosis of GI helminth infections of humans and livestock; (ii) examine high-throughput sequencing approaches, such as targeted molecular barcoding and shotgun sequencing, as tools to define the taxonomic composition of helminth infections; and (iii) discuss the current understanding of the interactions between helminths and microbiota in the host gut. Stool-based diagnostics are likely to serve as an important tool well into the future; improved diagnostics of helminths and their environment in the gut may assist the identification of biomarkers with the potential to define the health/disease status of individuals and populations, and to identify existing or emerging anthelmintic resistance.

DNA in extracellular vesicles: from evolution to its current application in health and disease

Extracellular vesicle (EV) secretion is a highly conserved evolutionary trait in all organisms in the three domains of life. The packaging and release of EVs appears to be a bulk-flow process which takes place mainly under extreme conditions. EVs participate in horizontal gene transfer, which supports the survival of prokaryotic and eukaryotic microbes. In higher eukaryotes, almost all cells secrete a heterogeneous population of EVs loaded with various biomolecules. EV secretion is typically higher in cancer microenvironments, promoting tumor progression and metastasis. EVs are now recognized as additional mediators of autocrine and paracrine communication in health and disease. In this context, proteins and RNAs have been studied the most, but extracellular vesicle DNA (EV-DNA) has started to gain in importance in the last few years. In this review, we summarize new findings related to the loading mechanism(s), localization, and post-shedding function of EV-DNA. We also discuss the feasibility of using EV-DNA as a biomarker when performing a liquid biopsy, at the same time emphasizing the lack of data from clinical trials in this regard. Finally, we outline the potential of EV-DNA uptake and its interaction with the host genome as a promising tool for understanding the mechanisms of cancer evolution.

Protecting DNA in membrane vesicles seems to be a conserved phenomenon for the horizontal genetic flux between prokaryotes and lower eukaryotes.

Capturing and analyzing this vesicular DNA enables quick and non-invasive monitoring of natural ecosystems.

Cancer-derived extracellular vesicles containing DNA open up novel directions in cell-to-cell communication and therefore disease monitoring.

Complex and fluctuating conditions of the tumor microenvironment, mimicking natural ecosystems, could favor EV-DNA release, mediating tumor multi-clonal evolution and providing survival benefits.

Using DNA Metabarcoding to Identify Floral Visitation by Pollinators

The identification of floral visitation by pollinators provides an opportunity to improve our understanding of the fine-scale ecological interactions between plants and pollinators, contributing to biodiversity conservation and promoting ecosystem health. In this review, we outline the various methods which can be used to identify floral visitation, including plant-focused and insect-focused methods. We reviewed the literature covering the ways in which DNA metabarcoding has been used to answer ecological questions relating to plant use by pollinators and discuss the findings of this research. We present detailed methodological considerations for each step of the metabarcoding workflow, from sampling through to amplification, and finally bioinformatic analysis. Detailed guidance is provided to researchers for utilisation of these techniques, emphasising the importance of standardisation of methods and improving the reliability of results. Future opportunities and directions of using molecular methods to analyse plant–pollinator interactions are then discussed.

Highly sensitive environmental DNA detection of topmouth gudgeon, Pseudorasbora parva: a comparison of qPCR and microfluidic qdPCR

Topmouth gudgeon is a freshwater fish species native to East Asia. Nowadays, P. parva is spread throughout Europe which is of concern because besides being considered one of the worst aquatic Invasive Alien Species (IAS) in Europe it is also a known vector of Spherotecum destruens, the rosette-like parasite lethal to other fish species. The present study describes the development and validation of a new species-specific assay based on hydrolysis probe chemistry to detect P. parva environmental DNA (eDNA) in water samples collected in a northern region of Italy (Friuli Venezia Giulia). Water samples were collected from 55 sites in an area where partial information on the occurrence of the species is available. eDNA was isolated from all samples and the presence of P. parva eDNA was tested by means of qPCR (quantitative PCR) and microfluidic qdPCR (quantitative digital PCR) techniques. Field results for both qPCR and qdPCR were largely in agreement in terms of detection (presence/absence). Thus, we judged the presence/absence by combining the results from the two methods and found that nine sites showed “strong positive” signal of P. parva eDNA (at least 2 positive replicates), 3 showed “suspected” (only 1 positive replicate), and 42 showed “absent”. The current study shows the strong potential of the newly developed eDNA approach to be a valuable addition to the monitoring of the highly invasive topmouth gudgeon in freshwater ecosystems.

Managing human-mediated range shifts: understanding spatial, temporal and genetic variation in marine non-native species

The use of molecular tools to manage natural resources is increasingly common. However, DNA-based methods are seldom used to understand the spatial and temporal dynamics of species' range shifts. This is important when managing range shifting species such as non-native species (NNS), which can have negative impacts on biotic communities. Here, we investigated the ascidian NNS Ciona robusta, Clavelina lepadiformis, Microcosmus squamiger and Styela plicata using a combined methodological approach. We first conducted non-molecular biodiversity surveys for these NNS along the South African coastline, and compared the results with historical surveys. We detected no consistent change in range size across species, with some displaying range stability and others showing range shifts. We then sequenced a section of cytochrome c oxidase subunit I (COI) from tissue samples and found genetic differences along the coastline but no change over recent times. Finally, we found that environmental DNA metabarcoding data showed broad congruence with both the biodiversity survey and the COI datasets, but failed to capture the complete incidence of all NNS. Overall, we demonstrated how a combined methodological approach can effectively detect spatial and temporal variation in genetic composition and range size, which is key for managing both thriving NNS and threatened species.

This article is part of the theme issue ‘Species’ ranges in the face of changing environments (part I)’.

Dynamics of microbiotic patterns reveal surface water groundwater interactions in intermittent and perennial streams

Exchange between groundwater (GW), hyporheic zone waters (HZ) and surface waters (SW) is critical for water quality, quantity, and the ecological health and functioning of all three ecosystems. Hydrological exchange is particularly important in intermittent creeks, such as in the Murray Darling Basin, Australia, where stream reaches shift from losing to gaining depending on the volume of surface flows. In this study we used hydrochemistry to identify SW-GW exchange and combined this with eDNA data to analyse the response of eukaryote and prokaryote communities to differing flow conditions within intermittent and perennial stream reaches. Our study suggested that SW and GW microbial communities were only around 30% similar. Differences in microbiota between SW, HZ and GW habitats were driven by changes in relative abundances of surface water dominant organisms (such as those capable of photosynthesis) as well as anaerobic taxa typical of GW environments (e.g., methanogens), with GW and HZ microbial communities becoming increasingly different to those in SW as flow ceased in intermittent creeks. Fine-scale hydrologic changes were identified through microbial communities in the perennial Maules Creek, indicating the importance of GW-SW exchange to biotic communities. This study highlights the importance of flow in shaping microbial communities and biogeochemical cycling within intermittent creeks and their connected alluvial aquifers. Our results suggest that microbiota may prove a useful indicator of SW-GW exchange, and in some circumstances, may be more sensitive in demonstrating fine-scale changes in SW-GW interactions than water chemistry. This knowledge furthers our understanding of GW-SW exchange and its impacts on ecological health.

Evaluating eDNA for Use within Marine Environmental Impact Assessments

In this review, the use of environmental DNA (eDNA) within Environmental Impact Assessment (EIA) is evaluated. EIA documents provide information required by regulators to evaluate the potential impact of a development project. Currently eDNA is being incorporated into biodiversity assessments as a complementary method for detecting rare, endangered or invasive species. However, questions have been raised regarding the maturity of the field and the suitability of eDNA information as evidence for EIA. Several key issues are identified for eDNA information within a generic EIA framework for marine environments. First, it is challenging to define the sampling unit and optimal sampling strategy for eDNA with respect to the project area and potential impact receptor. Second, eDNA assay validation protocols are preliminary at this time. Third, there are statistical issues around the probability of obtaining both false positives (identification of taxa that are not present) and false negatives (non-detection of taxa that are present) in results. At a minimum, an EIA must quantify the uncertainty in presence/absence estimates by combining series of Bernoulli trials with ad hoc occupancy models. Finally, the fate and transport of DNA fragments is largely unknown in environmental systems. Shedding dynamics, biogeochemical and physical processes that influence DNA fragments must be better understood to be able to link an eDNA signal with the receptor’s state. The biggest challenge is that eDNA is a proxy for the receptor and not a direct measure of presence. Nonetheless, as more actors enter the field, technological solutions are likely to emerge for these issues. Environmental DNA already shows great promise for baseline descriptions of the presence of species surrounding a project and can aid in the identification of potential receptors for EIA monitoring using other methods.

Using metabarcoding to assess Viridiplantae sequence diversity present in Antarctic glacial ice

Antarctica contains most of the glacial ice on the planet, a habitat that is largely unexplored by biologists. Recent warming in parts of Antarctica, particularly the Antarctic Peninsula region, is leading to widespread glacial retreat, releasing melt water and, potentially, contained biological material and propagules. In this study, we used a DNA metabarcoding approach to characterize Viridiplantae DNA present in Antarctic glacial ice. Ice samples from six glaciers in the South Shetland Islands and Antarctic Peninsula were analysed, detecting the presence of DNA representing a total of 16 taxa including 11 Chlorophyta (green algae) and five Magnoliophyta (flowering plants). The green algae may indicate the presence of a viable algal community in the ice or simply of preserved DNA, and the sequence diversity assigned included representatives of Chlorophyta not previously recorded in Antarctica. The presence of flowering plant DNA is most likely to be associated with pollen or tissue fragments introduced by humans.

Passive sampling of environmental DNA in aquatic environments using 3D-printed hydroxyapatite samplers

The study of environmental DNA released by aquatic organisms in their habitat offers a fast, non-invasive and sensitive approach to monitor their presence. Common eDNA sampling methods such as water filtration and DNA precipitation are time consuming, require difficult-to-handle equipment and partially integrate eDNA signals. To overcome these limitations, we created the first proof of concept of a passive, 3D-printed and easy-to-use eDNA sampler. We designed the samplers from hydroxyapatite (HAp samplers), a natural mineral with a high DNA adsorption capacity. The porous structure and shape of the samplers were designed to optimise DNA adsorption and facilitate their handling in the laboratory and in the field. Here we show that HAp samplers can efficiently collect genomic DNA in controlled set-ups, but can also collect animal eDNA under controlled and natural conditions with yields similar to conventional methods. However, we also observed large variations in the amount of DNA collected even under controlled conditions. A better understanding of the DNA-hydroxyapatite interactions on the surface of the samplers is now necessary to optimise the eDNA adsorption and to allow the development of a reliable, easy-to-use and reusable eDNA sampling tool.

Development and evaluation of PCR primers for environmental DNA (eDNA) metabarcoding of Amphibia

Biodiversity monitoring is important for the conservation of natural ecosystems in general, but particularly for amphibians, whose populations are pronouncedly declining. However, amphibians’ ecological traits (e.g. nocturnal or aquatic) often prevent their precise monitoring. Environmental DNA (eDNA) metabarcoding – analysis of extra-organismal DNA released into the environment – allows the easy and effective monitoring of the biodiversity of aquatic organisms. Here, we developed and tested the utility of original PCR primer sets. First, we conducted in vitro PCR amplification tests with universal primer candidates using total DNA extracted from amphibian tissues. Five primer sets successfully amplified the target DNA fragments (partial 16S rRNA gene fragments of 160–311 bp) from all 16 taxa tested (from the three living amphibian orders Anura, Caudata and Gymnophiona). Next, we investigated the taxonomic resolution retrieved using each primer set. The results revealed that the universal primer set “Amph16S” had the highest resolution amongst the tested sets. Finally, we applied Amph16S to the water samples collected in the field and evaluated its detection capability by comparing the species detected using eDNA and physical survey (capture-based sampling and visual survey) in multiple agricultural ecosystems across Japan (160 sites in 10 areas). The eDNA metabarcoding with Amph16S detected twice as many species as the physical surveys (16 vs. 8 species, respectively), indicating the effectiveness of Amph16S in biodiversity monitoring and ecological research for amphibian communities.

Environmental DNA captures native and non-native fish community variations across the lentic and lotic systems of a megacity

Globally, urbanization poses a major threat to terrestrial biodiversity, yet its impact on fish diversity is poorly understood, mainly because of surveying difficulties. In this study, environmental DNA metabarcoding was used to survey fish communities at 109 lentic and lotic sites across Beijing, and how environmental variables affect fish biodiversity at fine urban spatial scales was investigated. We identified 52 native and 23 non-native taxa, with lentic and lotic waters harboring both common and habitat-specific species. Water quality strongly affected native fish diversity, especially in lentic systems, but had little influence on non-native diversity. Fish diversity showed little response to urban land cover variation, but the relative sequence abundance of non-natives in lotic waters increased linearly with distance from the city center. Our findings illustrate the complex effects of urbanization on native versus non-native fishes in different aquatic habitats and highlight the distinctive considerations needed to conserve urban aquatic biodiversity.

Agricultural practices linked to shifts in groundwater microbial structure and denitrifying bacteria

Irrigation enhances the connectivity between the surface and groundwater by facilitating the transport of energy sources and oxygen. When combined with fertilisers, the impact on groundwater microbial communities and their interactions with nitrogen cycling in aquifers is poorly understood. This study examines the impact of different landuses (irrigated and non-irrigated) on groundwater microbial communities. A total of 38 wells accessing shallow aquifers in three sub-catchments of the Murray Darling Basin, Australia, were sampled for water chemistry and microbial community structure using environmental DNA (eDNA) techniques. All sub-catchments showed evidence of intense irrigation and groundwater contamination with total nitrogen, nitrates and phosphorus concentrations often well above background, with total nitrogen concentrations up to 70 mg/L and nitrate concentration up to 18 mg/L. Across sub-catchments there was high microbial diversity, with differences in community structure and function between catchments and landuses. Of the 1100 operational taxonomic units (OTUs) recorded, 47 OTUs were common across catchments with species from Woesearchaeota, Nitrospirales, Nitrosopumilales and Acidobacter taxonomic groups contributing greatly to groundwater microbial communities. Within non-irrigated sites, groundwaters contained similar proportions of nitrifying and denitrifying capable taxa, whereas irrigated sites had significantly higher abundances of microbes with nitrifying rather than denitrifying capabilities. Microbial diversity was lower in irrigated sites in the Macquarie catchment. These results indicate that irrigated landuses impact microbial community structure and diversity within groundwaters and suggest that the ratios of denitrifying to nitrifying capable microbes as well as specific orders (e.g., Nitrososphaerales) may be useful to indicate long-term nitrogen contamination of groundwaters. Such research is important for understanding the biogeochemical processes that are key predictors of redox state and contamination of groundwater by N species and other compounds. This will help to predict human impacts on groundwater microbial structure, diversity, and ecosystem functions, aiding the long-term management groundwater resources.

MetaCOXI: an integrated collection of metazoan mitochondrial cytochrome oxidase subunit-I DNA sequences

Nucleotide sequences reference collections or databases are fundamental components in DNA barcoding and metabarcoding data analyses pipelines. In such analyses, the accurate taxonomic assignment is a crucial aspect, relying directly on the availability of comprehensive and curated reference sequence collection and its taxonomy information. The currently wide use of the mitochondrial cytochrome oxidase subunit-I (COXI) as a standard DNA barcode marker in metazoan biodiversity studies highlights the need to shed light on the availability of the related relevant information from different data sources and their eventual integration. To adequately address data integration process, many aspects should be markedly considered starting from DNA sequence curation followed by taxonomy alignment with solid reference backbone and metadata harmonization according to universal standards. Here, we present MetaCOXI, an integrated collection of curated metazoan COXI DNA sequences with their associated harmonized taxonomy and metadata. This collection was built on the two most extensive available data resources, namely the European Nucleotide Archive (ENA) and the Barcode of Life Data System (BOLD). The current release contains more than 5.6 million entries (39.1% unique to BOLD, 3.6% unique to ENA, and 57.2% shared between both), their related taxonomic classification based on NCBI reference taxonomy, and their available main metadata relevant to environmental DNA studies, such as geographical coordinates, sampling country and host species. MetaCOXI is available in standard universal formats (‘fasta’ for sequences & ‘tsv’ for taxonomy and metadata), which can be easily incorporated in standard or specific DNA barcoding and/or metabarcoding data analysis pipelines.

Database URL: https://github.com/bachob5/MetaCOXI

Unearthing soil arthropod diversity through DNA metabarcoding

DNA metabarcoding has the potential to greatly advance understanding of soil biodiversity, but this approach has seen limited application for the most abundant and species-rich group of soil fauna-the arthropods. This study begins to address this gap by comparing information on species composition recovered from metabarcoding two types of bulk samples (specimens, soil) from a temperate zone site and from bulk soil samples collected at eight sites in the Arctic. Analysis of 22 samples (3 specimen, 19 soil) revealed 410 arthropod OTUs belonging to 112 families, 25 orders, and nine classes. Studies at the temperate zone site revealed little overlap in species composition between soil and specimen samples, but more overlap at higher taxonomic levels (families, orders) and congruent patterns of α- and β-diversity. Expansion of soil analyses to the Arctic revealed locally rich, highly dissimilar, and spatially structured assemblages compatible with dispersal limited and environmentally driven assembly. The current study demonstrates that DNA metabarcoding of bulk soil enables rapid, large-scale assessments of soil arthropod diversity. However, deep sequence coverage is required to adequately capture the species present in these samples, and expansion of the DNA barcode reference library is necessary to improve taxonomic resolution of the sequences recovered through this approach.