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

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.

ExTaxsI: an exploration tool of biodiversity molecular data

BACKGROUND

The increasing availability of multi-omics data is leading to regularly revised estimates of existing biodiversity data. In particular, the molecular data enable novel species to be characterized and the information linked to those already observed to be increased with new genomics data. For this reason, the management and visualization of existing molecular data, and their related metadata, through the implementation of easy-to-use IT tools have become a key point to design future research. The more users are able to access biodiversity-related information, the greater the ability of the scientific community to expand its knowledge in this area.

RESULTS

In this article we focus on the development of ExTaxsI (Exploring Taxonomy Information), an IT tool that can retrieve biodiversity data stored in NCBI databases and provide a simple and explorable visualization. We use 3 case studies to show how an efficient organization of the available data can lead to obtaining new information that is fundamental as a starting point for new research. Using this approach highlights the limits in the distribution of data availability, a key factor to consider in the experimental design phase of broad-spectrum studies such as metagenomics.

CONCLUSIONS

ExTaxsI can easily retrieve molecular data and its metadata with an explorable visualization, with the aim of helping researchers to improve experimental designs and highlight the main gaps in the coverage of available data.

Meta-analysis shows both congruence and complementarity of DNA and eDNA metabarcoding to traditional methods for biological community assessment

DNA metabarcoding is increasingly used for the assessment of aquatic communities, and numerous studies have investigated the consistency of this technique with traditional morpho‐taxonomic approaches. These individual studies have used DNA metabarcoding to assess diversity and community structure of aquatic organisms both in marine and freshwater systems globally over the last decade. However, a systematic analysis of the comparability and effectiveness of DNA‐based community assessment across all of these studies has hitherto been lacking. Here, we performed the first meta‐analysis of available studies comparing traditional methods and DNA metabarcoding to measure and assess biological diversity of key aquatic groups, including plankton, microphytobentos, macroinvertebrates, and fish. Across 215 data sets, we found that DNA metabarcoding provides richness estimates that are globally consistent to those obtained using traditional methods, both at local and regional scale. DNA metabarcoding also generates species inventories that are highly congruent with traditional methods for fish. Contrastingly, species inventories of plankton, microphytobenthos and macroinvertebrates obtained by DNA metabarcoding showed pronounced differences to traditional methods, missing some taxa but at the same time detecting otherwise overseen diversity. The method is generally sufficiently advanced to study the composition of fish communities and replace more invasive traditional methods. For smaller organisms, like macroinvertebrates, plankton and microphytobenthos, DNA metabarcoding may continue to give complementary rather than identical estimates compared to traditional approaches. Systematic and comparable data collection will increase the understanding of different aspects of this complementarity, and increase the effectiveness of the method and adequate interpretation of the results.

DnoisE: distance denoising by entropy. An open-source parallelizable alternative for denoising sequence datasets

DNA metabarcoding is broadly used in biodiversity studies encompassing a wide range of organisms. Erroneous amplicons, generated during amplification and sequencing procedures, constitute one of the major sources of concern for the interpretation of metabarcoding results. Several denoising programs have been implemented to detect and eliminate these errors. However, almost all denoising software currently available has been designed to process non-coding ribosomal sequences, most notably prokaryotic 16S rDNA. The growing number of metabarcoding studies using coding markers such as COI or RuBisCO demands a re-assessment and calibration of denoising algorithms. Here we present DnoisE, the first denoising program designed to detect erroneous reads and merge them with the correct ones using information from the natural variability (entropy) associated to each codon position in coding barcodes. We have developed an open-source software using a modified version of the UNOISE algorithm. DnoisE implements different merging procedures as options, and can incorporate codon entropy information either retrieved from the data or supplied by the user. In addition, the algorithm of DnoisE is parallelizable, greatly reducing runtimes on computer clusters. Our program also allows different input file formats, so it can be readily incorporated into existing metabarcoding pipelines.

Land-use type temporarily affects active pond community structure but not gene expression patterns

Changes in land use and agricultural intensification threaten biodiversity and ecosystem functioning of small water bodies. We studied 67 kettle holes (KH) in an agricultural landscape in northeastern Germany using landscape-scale metatranscriptomics, to understand the responses of active bacterial, archaeal, and eukaryotic communities, to land-use type. These KH are proxies of the millions of small standing water bodies of glacial origin spread across the northern hemisphere. Like other landscapes in Europe, the study area has been used for intensive agriculture since the 1950s. In contrast to a parallel eDNA study which revealed the homogenization of biodiversity across KH conceivably resulting from long-lasting intensive agriculture, land-use type affected the structure of the active KH communities during spring crop fertilization, but not a month later. This effect was more pronounced in eukaryotes than in bacteria. In contrast, gene expression patterns did not differ between months or across land-use type, suggesting a high degree of functional redundancy across the KH communities. Variability in gene expression was best explained by active bacterial and eukaryotic community structures, suggesting that these changes in functioning are primarily driven by interactions between organisms. Our results show that influences of the surrounding landscape result in temporary changes in the activity of different community members. Thus, even in KH where biodiversity has been homogenized, communities continue to respond to land management. This needs to be considered when developing sustainable management options for restoration purposes and for successful mitigation of further biodiversity loss in agricultural landscapes.

Landscape Characteristics Affecting Small Mammal Occurrence in Heterogeneous Olive Grove Agro-Ecosystems

Understanding how small mammals (SM) are associated with environmental characteristics in olive groves is important to identify potential threats to agriculture and assess the overall conservation value and functioning of agro-ecosystems. Here, we provide first insights on this topic applied to traditional olive groves in northeast (NE) Portugal by assessing the landscape attributes that determine SM occurrence, focusing on one species of conservation concern (Microtus cabrerae Thomas 1906) and one species often perceived as a potential pest of olives (Microtus lusitanicus Gerbe 1879). Based on SM genetic non-invasive sampling in 51 olive groves and surrounding habitats, we identified seven rodent species and one insectivore. Occupancy modelling indicated that SM were generally less detected within olive groves than in surrounding habitats. The vulnerable M. cabrerae reached a mean occupancy (95% CI) of 0.77 (0.61–0.87), while M. lusitanicus stood at 0.37 (0.24–0.52). M. cabrerae was more likely to occur in land mosaics with high density of agricultural field edges, while M. lusitanicus was more associated with high density of pastureland patches. Overall, our study suggests that the complex structure and spatial heterogeneity of traditionally managed olive grove agro-ecosystems may favor the occurrence of species-rich SM communities, possibly including well-established populations of species of conservation importance, while keeping potential pest species at relatively low occupancy rates.

eDNA profiling of mammals, birds, and fish of surface waters by mitochondrial metagenomics: application for source tracking of fecal contamination in surface waters

Knowing the composition of animals present in aquatic ecosystems can tell us about the anthropic pressures on these environments. One of these pressures is the occurrence of fecal contamination. However, this contamination can originate from more than one animal species in areas where urban and agricultural activities overlap. Mitochondrial DNA (mtDNA) has become the standard barcoding tool to identify the presence of animal species in environment. Amplicon-sequencing metagenomics is a powerful approach to derive the animal profile in an environment. However, PCR primers targeting mtDNA of a broad range of animals are highly degenerate or generate short DNA fragments that could cause ambiguous affiliation. Here we report the development of a new set of primers targeting the mitochondrial 16S ribosomal RNA genes of a broad range of terrestrial and aquatic animals, which include mammals, birds, and fishes. These primers successfully amplified mtDNA from environmental DNA (eDNA) extracted from surface waters. Sequencing the resulting amplicons revealed the presence of mammals and birds that may contribute in fecal contamination of surface water. In one of the river samples high in fecal indicator bacteria, human and bovine mtDNA accounted for 40.5% and 4.1% of the sequences, respectively, suggesting fecal contamination by these two animals. These findings indicate that our PCR primers coupled with amplicon-sequencing metagenomics contribute in profiling the animal diversity in the surface waters and its surrounding. This approach could be a valuable tool to identify simultaneously the potential contribution of various animals as sources of fecal contamination in surface waters.

Unveiling biogeographical patterns of the ichthyofauna in the Tuichi basin, a biodiversity hotspot in the Bolivian Amazon, using environmental DNA

To date, more than 2400 valid fish species have been recorded in the Amazon basin. However, some regions remain poorly documented. This is the case in the Beni basin and in particular in one of its main sub-basins, the Tuichi, an Andean foothills rivers flowing through the Madidi National Park in the Bolivian Amazonia. The knowledge of its ichthyological diversity is, however, essential for the management and protection of aquatic ecosystems, which are threatened by the development of infrastructures (dams, factories and cities), mining and deforestation. Environmental DNA (eDNA) has been relatively little used so far in the Amazon basin. We sampled eDNA from water in 34 sites in lakes and rivers in the Beni basin including 22 sites in the Tuichi sub-basin, during the dry season. To assess the biogeographical patterns of the amazonian ichthyofauna, we implemented a metabarcoding approach using two pairs of specific primers designed and developed in our laboratory to amplify two partially overlapping CO1 fragments, one of 185bp and another of 285bp. We detected 252 fish taxa (207 at species level) among which 57 are newly identified for the Beni watershed. Species compositions are significantly different between lakes and rivers but also between rivers according to their hydrographic rank and altitude. Furthermore, the diversity patterns are related to the different hydro-ecoregions through which the Tuichi flows. The eDNA approach makes it possible to identify and complete the inventory of the ichthyofauna in this still poorly documented Amazon basin. However, taxonomic identification remains constrained by the lack of reference barcodes in public databases and does not allow the assignment of all OTUs. Our results can be taken into account in conservation and management strategies and could serve as a baseline for future studies, including on other Andean tributaries.

Environmental DNA Metabarcoding: A Novel Method for Biodiversity Monitoring of Marine Fish Communities

Environmental DNA (eDNA) is genetic material that has been shed from macroorganisms. It has received increased attention as an indirect marker for biodiversity monitoring. This article reviews the current status of eDNA metabarcoding (simultaneous detection of multiple species) as a noninvasive and cost-effective approach for monitoring marine fish communities and discusses the prospects for this growing field. eDNA metabarcoding coamplifies short fragments of fish eDNA across a wide variety of taxa and, coupled with high-throughput sequencing technologies, allows massively parallel sequencing to be performed simultaneously for dozens to hundreds of samples. It can predict species richness in a given area, detect habitat segregation and biogeographic patterns from small to large spatial scales, and monitor the spatiotemporal dynamics of fish communities. In addition, it can detect an anthropogenic impact on fish communities through evaluation of their functional diversity. Recognizing the strengths and limitations of eDNA metabarcoding will help ensure that continuous biodiversity monitoring at multiple sites will be useful for ecosystem conservation and sustainable use of fishery resources, possibly contributing to achieving the targets of the United Nations' Sustainable Development Goal 14 for 2030.

Phytobiome Metabarcoding: A Tool to Help Identify Prokaryotic and Eukaryotic Causal Agents of Undiagnosed Tree Diseases

Recent advancements in high-throughput sequencing have provided scientists with vastly enhanced tools to diagnose unknown tree diseases. One of these techniques is referred to as metabarcoding, which uses phylogenetically informative reference genes to taxonomically classify short DNA sequences amplified from environmental samples. Using metabarcoding, we are able to compare the microbiota of symptomatic and asymptomatic (including presumably naïve) samples and identify microbe(s) that are only present in symptomatic samples and could therefore be responsible for the undiagnosed disease. Metabarcoding involves two main steps: library preparation and bioinformatic processing. For library preparation, the appropriate reference gene for the organism of interest (i.e., bacteria, phytoplasma, fungi, or other eukaryotes, such as nematodes) is amplified from the DNA extracted from the environmental samples using PCR and prepared for sequencing. The bioinformatic processing includes four major steps: (1) quality check and cleanup on raw reads; (2) classification of the sequences into taxonomically informative groups (ASVs or OTUs); (3) taxonomy assignments based on the reference database; and (4) differential abundance and diversity analyses to identify microbes that are significantly associated with just symptomatic samples and that point toward the putative causal agent of the disease.

Experimental validation of taxon-specific mini-barcode primers for metabarcoding of zooplankton

Metabarcoding to determine the species composition and diversity of marine zooplankton communities is a fast-developing field in which the standardization of methods is yet to be fully achieved. The selection of genetic markers and primer choice are particularly important because they substantially influence species detection rates and accuracy. Validation is therefore an important step in the design of metabarcoding protocols. We developed taxon-specific mini-barcode primers for the cytochrome c oxidase subunit I (COI) gene region and used an experimental approach to test species detection rates and primer accuracy of the newly designed primers for prawns, shrimps and crabs and published primers for marine lobsters and fish. Artificially assembled mock communities (with known species ratios) and unsorted coastal tow-net zooplankton samples were sequenced and the detected species were compared with those seeded in mock communities to test detection rates. Taxon-specific primers increased detection rates of target taxa compared with a universal primer set. Primer cocktails (multiple primer sets) significantly increased species detection rates compared with single primer pairs and could detect up to 100% of underrepresented target taxa in mock communities. Taxon-specific primers recovered fewer false-positive or false-negative results than the universal primer. The methods used to design taxon-specific mini-barcodes and the experimental mock community validation protocols shown here can easily be applied to studies on other groups and will allow for a level of standardization among studies undertaken in different ecosystems or geographic locations.

Assessment of marine benthic diatom communities: insights from a combined morphological-metabarcoding approach in Mediterranean shallow coastal waters

We investigated the advantages and disadvantages of light microscope (LM)-based identifications and DNA metabarcoding, based on a 312-bp rbcL marker, for examining benthic diatom communities from Mediterranean shallow coastal environments. For this, we used biofilm samples collected from different substrata in the Ebro delta bays. We show that 1) Ebro delta bays harbour high-diversity diatom communities [LM identified 249 taxa] and 2) DNA metabarcoding effectively reflects this diversity at genus- but not species level, because of the incompleteness of the DNA reference library. Nevertheless, DNA metabarcoding offers new opportunities for detecting small, delicate and rare diatom species missed by LM and diatoms that lack silica frustules. The primers used, though designed for diatoms, successfully amplified rarely reported members of other stramenopile groups. Combining LM and DNA approaches offers stronger support for ecological studies of benthic microalgal communities in shallow coastal environments than using either approach on its own.

Diversity of Fungi Present in Permafrost in the South Shetland Islands, Maritime Antarctic

We assess the fungal diversity present in permafrost from different islands in the South Shetland Islands archipelago, maritime Antarctic, using next-generation sequencing (NGS). We detected 1,003,637 fungal DNA reads representing, in rank abundance order, the phyla Ascomycota, Mortierellomycota, Basidiomycota, Chytridiomycota, Rozellomycota, Mucoromycota, Calcarisporiellomycota and Zoopagomycota. Ten taxa were dominant these being, in order of abundance, Pseudogymnoascus appendiculatus, Penicillium sp., Pseudogymnoascus roseus, Penicillium herquei, Curvularia lunata, Leotiomycetes sp., Mortierella sp. 1, Mortierella fimbricystis, Fungal sp. 1 and Fungal sp. 2. A further 38 taxa had intermediate abundance and 345 were classified as rare. The total fungal community detected in the permafrost showed high indices of diversity, richness and dominance, although these varied between the sampling locations. The use of a metabarcoding approach revealed the presence of DNA of a complex fungal assemblage in the permafrost of the South Shetland Islands including taxa with a range of ecological functions among which were multiple animal, human and plant pathogenic fungi. Further studies are required to determine whether the taxa identified are present in the form of viable cells or propagules and which might be released from melting permafrost to other Antarctic habitats and potentially dispersed more widely.

Mapping differences in mammalian distributions and diversity using environmental DNA from rivers

Finding more efficient ways to monitor and estimate the diversity of mammalian communities is a major step towards their management and conservation. Environmental DNA (eDNA) from river water has recently been shown to be a viable method for biomonitoring mammalian communities. Most of the studies to date have focused on the potential for eDNA to detect individual species, with little focus on describing patterns of community diversity and structure. Here, we first focus on the sampling effort required to reliably map the diversity and distribution of semi-aquatic and terrestrial mammals and allow inferences of community structure surrounding two rivers in southeastern England. Community diversity and composition was then assessed based on species richness and β-diversity, with differences between communities partitioned into nestedness and turnover, and the sampling effort required to rapidly detect semi-aquatic and terrestrial species was evaluated based on species accumulation curves and occupancy modelling. eDNA metabarcoding detected 25 wild mammal species from five orders, representing the vast majority (82%) of the species expected in the area. The required sampling effort varied between orders, with common species (generally rodents, deer and lagomorphs) more readily detected, with carnivores detected less frequently. Measures of species richness differed between rivers (both overall and within each mammalian order) and patterns of β-diversity revealed the importance of species replacement in sites within each river, against a pattern of species loss between the two rivers. eDNA metabarcoding demonstrated its capability to rapidly detect mammal species, allowing inferences of community composition that will better inform future sampling strategies for this Class. Importantly, this study highlights the potential use of eDNA data for investigating mammalian community dynamics over different spatial scales.

Phytool, a ShinyApp to homogenise taxonomy of freshwater microalgae from DNA barcodes and microscopic observations

Methods for biomonitoring of freshwater phytoplankton are evolving rapidly with eDNA-based methods, offering great complementarity with microscopy. Metabarcoding approaches have been more commonly used over the last years, with a continuous increase in the amount of data generated. Depending on the researchers and the way they assigned barcodes to species (bioinformatic pipelines and molecular reference databases), the taxonomic assignment obtained for HTS DNA reads might vary. This is also true for traditional taxonomic studies by microscopy with regular adjustments of the classification and taxonomy.

For those reasons (leading to non-homogeneous taxonomies), gap-analyses and comparisons between studies become even more challenging and the curation processes to find potential consensus names are time-consuming. Here, we present a web-based application (Phytool), developed with ShinyApp (Rstudio), that aims to make the harmonisation of taxonomy easier and in a more efficient way, using a complete and up-to-date taxonomy reference database for freshwater microalgae. Phytool allows users to homogenise and update freshwater phytoplankton taxonomical names from sequence files and data tables directly uploaded in the application. It also gathers barcodes from curated references in a user-friendly way in which it is possible to search for specific organisms. All the data provided are downloadable with the possibility to apply filters in order to select only the required taxa and fields (e.g. specific taxonomic ranks). The main goal is to make accessible to a broad range of users the connection between microscopy and molecular biology and taxonomy through different ready-to-use functions. This study estimates that only 25% of species of freshwater phytoplankton in Phytobs are associated with a barcode. We plead for an increased effort to enrich reference databases by coupling taxonomy and molecular methods. Phytool should make this crucial work more efficient.

The application is available at https://caninuzzo.shinyapps.io/phytool_v1/

Current limitations and future prospects of detection and biomonitoring of NIS in the Mediterranean Sea through environmental DNA

The biodiversity of the Mediterranean Sea is currently threatened by the introduction of Non-Indigenous Species (NIS). Therefore, monitoring the distribution of NIS is of utmost importance to preserve the ecosystems. A promising approach for the identification of species and the assessment of biodiversity is the use of DNA barcoding, as well as DNA and eDNA metabarcoding. Currently, the main limitation in the use of genomic data for species identification is the incompleteness of the DNA barcode databases. In this research, we assessed the availability of DNA barcodes in the main reference libraries for the most updated inventory of 665 confirmed NIS in the Mediterranean Sea, with a special focus on the cytochrome oxidase I (COI) barcode and primers. The results of this study show that there are no barcodes for 33.18% of the species in question, and that 45.30% of the 382 species with COI barcode, have no primers publicly available. This highlights the importance of directing scientific efforts to fill the barcode gap of specific taxonomic groups in order to help in the effective application of the eDNA technique for investigating the occurrence and the distribution of NIS in the Mediterranean Sea.

Deep Learning Classification of Lake Zooplankton

Plankton are effective indicators of environmental change and ecosystem health in freshwater habitats, but collection of plankton data using manual microscopic methods is extremely labor-intensive and expensive. Automated plankton imaging offers a promising way forward to monitor plankton communities with high frequency and accuracy in real-time. Yet, manual annotation of millions of images proposes a serious challenge to taxonomists. Deep learning classifiers have been successfully applied in various fields and provided encouraging results when used to categorize marine plankton images. Here, we present a set of deep learning models developed for the identification of lake plankton, and study several strategies to obtain optimal performances, which lead to operational prescriptions for users. To this aim, we annotated into 35 classes over 17900 images of zooplankton and large phytoplankton colonies, detected in Lake Greifensee (Switzerland) with the Dual Scripps Plankton Camera. Our best models were based on transfer learning and ensembling, which classified plankton images with 98% accuracy and 93% F1 score. When tested on freely available plankton datasets produced by other automated imaging tools (ZooScan, Imaging FlowCytobot, and ISIIS), our models performed better than previously used models. Our annotated data, code and classification models are freely available online.

Current status of Phytophthora in Australia

Among the most economically relevant and environmentally devastating diseases globally are those caused by Phytophthora species. In Australia, production losses in agriculture and forestry result from several well-known cosmopolitan Phytophthora species and infestation of natural ecosystems by Phytophthora cinnamomi have caused irretrievable loss to biodiversity especially in proteaceous dominated heathlands. For this review, all available records of Phytophthora in Australia were collated and curated, resulting in a database of 7 869 records, of which 2 957 have associated molecular data. Australian databases hold records for 99 species, of which 20 are undescribed. Eight species have no records linked to molecular data, and their presence in Australia is considered doubtful. The 99 species reside in 10 of the 12 clades recognised within the complete phylogeny of Phytophthora. The review includes discussion on each of these species' status and additional information provided for another 29 species of concern. The first species reported in Australia in 1900 was Phytophthora infestans. By 2000, 27 species were known, predominantly from agriculture. The significant increase in species reported in the subsequent 20 years has coincided with extensive surveys in natural ecosystems coupled with molecular taxonomy and the recognition of numerous new phylogenetically distinct but morphologically similar species. Routine and targeted surveys within Australian natural ecosystems have resulted in the description of 27 species since 2009. Due to the new species descriptions over the last 20 years, many older records have been reclassified based on molecular identification. The distribution of records is skewed toward regions with considerable activity in high productivity agriculture, horticulture and forestry, and native vegetation at risk from P. cinnamomi. Native and exotic hosts of different Phytophthora species are found throughout the phylogeny; however, species from clades 1, 7 and 8 are more likely to be associated with exotic hosts. One of the most difficult challenges to overcome when establishing a pest status is a lack of reliable data on the current state of a species in any given country or location. The database compiled here for Australia and the information provided for each species overcomes this challenge. This review will aid federal and state governments in risk assessments and trade negotiations by providing a comprehensive resource on the current status of Phytophthora species in Australia. Citation: Burgess TI, Edwards J, Drenth A, et al. 2021. Current status of Phytophthora in Australia. Persoonia 47: 151-177. https://doi.org/10.3767/persoonia.2021.47.05.

Using environmental DNA methods to improve detectability in an endangered sturgeon (Acipenser sinensis) monitoring program

Background

To determine the presence and abundance of an aquatic species in large waterbodies, especially when populations are at low densities, is highly challenging for conservation biologists. Environmental DNA (eDNA) has the potential to offer a noninvasive and cost-effective method to complement traditional population monitoring, however, eDNA has not been extensively applied to study large migratory species. Chinese sturgeon (Acipenser sinensis), is the largest anadromous migratory fish in the Yangtze River, China, and in recent years its population has dramatically declined and spawning has failed, bringing this species to the brink of extinction. In this study, we aim to test the detectability of eDNA methods to determine the presence and relative abundance of reproductive stock of the species and whether eDNA can be used as a tool to reflect behavioral patterns. Chinese sturgeon eDNA was collected from four sites along the spawning ground across an eight month period, to investigate the temporal and spatial distribution using droplet digital PCR (ddPCR).

Results

We designed a pair of specific primers for Chinese sturgeon and demonstrated the high sensitivity of ddPCR to detect and quantify the Chinese sturgeon eDNA concentration with the limit of detection 0.17 copies/μl, with Chinese sturgeon eDNA been intermittently detected at all sampling sites. There was a consistent temporal pattern among four of the sampling sites that could reflect the movement characteristics of the Chinese sturgeon in the spawning ground, but without a spatial pattern. The eDNA concentration declined by approximately 2–3 × between December 2018 and December 2019.

Conclusions

The results prove the efficacy of eDNA for monitoring reproductive stock of the Chinese sturgeon and the e decreased eDNA concentration reflect that Chinese sturgeon may survive with an extremely small number of reproductive stock in the Yangtze River. Accordingly, we suggest future conservation measures should focus on both habitat restoration and matured fish restocking to ensure successful spawning. Overall, this study provides encouraging support for the application of eDNA methods to monitor endangered aquatic species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01948-w.

DNA metabarcoding reveals differences in distribution patterns and ecological preferences among genetic variants within some key freshwater diatom species

Our study evaluates differences in the distribution and ecology of genetic variants within several ecologically important diatom species that are also key for Water Framework Directive monitoring of European rivers: Fistulifera saprophila (FSAP), Achnanthidium minutissimum (ADMI), Nitzschia inconspicua (NINC) and Nitzschia soratensis (NSTS). We used DADA2 to infer amplicon sequence variants (ASVs) of a short rbcL barcode in 531 environmental samples from biomonitoring campaigns in Catalonia and France. ASVs within each species showed different distribution patterns. Threshold Indicator Taxa ANalysis revealed three ecological groupings of ASVs in both ADMI and FSAP. Two of these in each species were separated by opposite responses to calcium and conductivity. Boosted regression trees additionally showed that both variables greatly influenced the occurrence of these groupings. A third grouping in FSAP was characterized by a negative response to total organic carbon and hence was better represented in waters with higher ecological status than the other FSAP ASVs, contrasting with what is generally assumed for the species. In the two Nitzschia species, our analyses confirmed earlier studies: NINC preferred higher levels of calcium and conductivity. Our findings suggest that the broad ecological tolerance of some diatom species results from overlapping preferences among genetic variants, which individually show much more restricted preferences and distributions. This work shows the importance of studying the ecological preferences of genetic variants within species complexes, now possible with DNA metabarcoding. The results will help reveal and understand biogeographical distributions and facilitate the development of more accurate biological indexes for biomonitoring programmes.

Use of environmental DNA in assessment of fish functional and phylogenetic diversity

Assessing the impact of global changes and protection effectiveness is a key step in monitoring marine fishes. Most traditional census methods are demanding or destructive. Nondisturbing and nonlethal approaches based on video and environmental DNA are alternatives to underwater visual census or fishing. However, their ability to detect multiple biodiversity factors beyond traditional taxonomic diversity is still unknown. For bony fishes and elasmobranchs, we compared the performance of eDNA metabarcoding and long-term remote video to assess species' phylogenetic and functional diversity. We used 10 eDNA samples from 30 L of water each and 25 hr of underwater videos over 4 days on Malpelo Island (pacific coast of Colombia), a remote marine protected area. Metabarcoding of eDNA detected 66% more molecular operational taxonomic units (MOTUs) than species on video. We found 66 and 43 functional entities with a single eDNA marker and videos, respectively, and higher functional richness for eDNA than videos. Despite gaps in genetic reference databases, eDNA also detected a higher fish phylogenetic diversity than videos; accumulation curves showed how 1 eDNA transect detected as much phylogenetic diversity as 25 hr of video. Environmental DNA metabarcoding can be used to affordably, efficiently, and accurately census biodiversity factors in marine systems. Although taxonomic assignments are still limited by species coverage in genetic reference databases, use of MOTUs highlights the potential of eDNA metabarcoding once reference databases have expanded.

Assessing the effects of cascade dams on river ecological status using multi-species interaction-based index of biotic integrity (Mt-IBI)

Cascade dams have exerted significant effects on river ecosystems. To quantitatively assess dam-induced effects on river ecological status, a novel multi-species interaction-based index of biotic integrity (Mt-IBI) was developed. Benthic microbiota was selected as a bio-indicator for its sensitivity to the environmental disturbance. An environmental DNA metabarcoding tool was used to identify microbiota (bacteria, protozoan, and metazoan). The Mt-IBI was applied to assess the ecological status of the Hanjiang River, a representative dam-affected river in China. Fifteen sampling sites along the Hanjiang River were sampled in June 2018. Seven core metrics were screened from a total of 364 candidate metrics to calculate the value of the Mt-IBI. The Mt-IBI of the Hanjiang River ranged from 1.90 to 6.39, with a mean value of 4.02. The mean values of Mt-IBI at the reservoir and riverine side of dams were 2.11 and 3.81, respectively. The downstream reach without dam constructions had the highest mean Mt-IBI (5.79). Thus, the continuity of the river was strongly related to the Mt-IBI. Structural equation models (SEMs) were further established to identify the dominant environmental variables in the dam-affected river. The SEMs indicated that flow velocity (coefficient 0.749) was the most important determinant of ecological status in the Hanjiang River. Water organic matter also played a vital role in determining the ecological status of the Hanjiang River, and exerted the strongest direct effect (P < 0.001, r = 0.712). The reliability of SEMs was verified by building a support vector regression model (R² = 0.8141). This study can provide new tools for ecological assessment and diagnosis, and provide a new perspective for the management of cascade dams.

Mitochondrial Genetic Diversity, Population Structure and Detection of Antillean and Amazonian Manatees in Colombia: New Areas and New Techniques

The Antillean manatee (Trichechus manatus) and the Amazonian manatee (Trichechus inunguis) are distributed in rivers in the Caribbean and Amazonian region of Colombia respectively. For 30 years, genetic information has been obtained from these populations in order to inform conservation programs for these endangered species and decide on the location to release them back to the wild. However, in previous studies, samples from rivers in some areas of the country were not included, given the difficulties to access these regions due to either logistic or safety issues. In this study, we analyzed mitochondrial DNA (mtDNA) control region (CR) sequences of from samples of T. manatus (n = 37) and T. inunguis (n = 4) (410 and 361 bp, respectively), obtained in new and previously unexplored rivers and bays in the country, including Santa Marta, Urabá Gulf, Ayapel Marsh (San Jorge River Basin), Meta River and Magdalena Medio and the low Magdalena River (Cesar Province and Canal del Dique) as well as additional samples from Puerto Nariño in the Colombian Amazon. Our results included the discovery of two newly described mtDNA CR haplotypes for T. manatus. In addition, we confirmed significant population differentiation at the mitochondrial level between the Magdalena and Sinú rivers and differentiation among areas of the same river, including the middle and low Magdalena River. This differentiation may be related to anthropic changes in the river since construction of the Canal del Dique in the XVI century. We also tested environmental DNA sampling and analyses techniques to evaluate its potential use for manatee detection and monitoring in bodies of water in Colombia, in order to evaluate new areas for future manatee conservation initiatives. We emphasize the need to continue using genetic information to provide evidence on the potential best locations to undertake animal release to prevent outbreeding depression.

Monitoring newt communities in urban area using eDNA metabarcoding

Newts are amphibians commonly present in small ponds or garden pools in urban areas. They are protected in many countries and their presence is monitored through visual observation and/or trapping. However, newts are not easy to spot as they are small, elusive and often hidden at the bottom of water bodies. In recent years, environmental DNA (eDNA) has become a popular tool for detecting newts, with a focus on individual species using qPCR assays. Here, we assess the effectiveness of eDNA metabarcoding compared to conventional visual surveys of newt diversity in 45 ponds within urban areas of Geneva canton, Switzerland. We designed newt-specific mitochondrial 16S rRNA primers, which assign the majority of amplicons to newts, and were able to detect four species known to be present in the region, including the invasive subspecies Lissotriton vulgaris meridionalis, native to the Italian peninsula, that has been introduced in the Geneva area recently. The obtained eDNA results were congruent overall with conventional surveys, confirming the morphological observations in the majority of cases (67%). In 25% of cases, a species was only detected genetically, while in 8% of cases, the observations were not supported by eDNA metabarcoding. Our study confirms the usefulness of eDNA metabarcoding as a tool for the effective and non-invasive monitoring of newt community and suggests its broader use for the survey of newt diversity in urban area at larger scales.

Reinforcement of Environmental DNA Based Methods (Sensu Stricto) in Biodiversity Monitoring and Conservation: A Review

Simple Summary

Worldwide biodiversity loss points to a necessity of upgrading to a fast and effective monitoring method that can provide quick conservation action. Newly developed environmental DNA (eDNA) based method found to be more cost-effective, non-invasive, quick, and accurate than traditional monitoring (spot identification, camera trapping). Although the eDNA based methods are proliferating rapidly, as a newly developed branch, it needs more standardization and practitioner adaptation. The present study aims to evaluate the eDNA based methods, and their potential achievements in biodiversity monitoring, and conservation for quick practitioners’ adaption. The investigation shows that the eDNA technique is applicable largely in (i) early detection of invasive species, (ii) species detection for conservation, (iii) community-level biodiversity monitoring, (iv) ecosystem health monitoring, (v) study on trophic interactions, etc. Thus, the eDNA technique shows a great promise with its high accuracy and authenticity, and will be applicable alone or alongside other methods in the near future.

Abstract

Recently developed non-invasive environmental DNA-based (eDNA) techniques have enlightened modern conservation biology, propelling the monitoring/management of natural populations to a more effective and efficient approach, compared to traditional surveys. However, due to rapid-expansion of eDNA, confusion in terminology and collection/analytical pipelines can potentially jeopardize research progression, methodological standardization, and practitioner adoption in several ways. Present investigation reflects the developmental progress of eDNA (sensu stricto) including highlighting the successful case studies in conservation management. The eDNA technique is successfully relevant in several areas of conservation research (invasive/conserve species detection) with a high accuracy and authentication, which gradually upgrading modern conservation approaches. The eDNA technique related bioinformatics (e.g., taxon-specific-primers MiFish, MiBird, etc.), sample-dependent methodology, and advancement of sequencing technology (e.g., oxford-nanopore-sequencing) are helping in research progress. The investigation shows that the eDNA technique is applicable largely in (i) early detection of invasive species, (ii) species detection for conservation, (iii) community level biodiversity monitoring, (iv) ecosystem health monitoring, (v) study on trophic interactions, etc. Thus, the eDNA technique with a high accuracy and authentication can be applicable alone or coupled with traditional surveys in conservation biology. However, a comprehensive eDNA-based monitoring program (ecosystem modeling and function) is essential on a global scale for future management decisions.

RESCRIPt: Reproducible sequence taxonomy reference database management

Nucleotide sequence and taxonomy reference databases are critical resources for widespread applications including marker-gene and metagenome sequencing for microbiome analysis, diet metabarcoding, and environmental DNA (eDNA) surveys. Reproducibly generating, managing, using, and evaluating nucleotide sequence and taxonomy reference databases creates a significant bottleneck for researchers aiming to generate custom sequence databases. Furthermore, database composition drastically influences results, and lack of standardization limits cross-study comparisons. To address these challenges, we developed RESCRIPt, a Python 3 software package and QIIME 2 plugin for reproducible generation and management of reference sequence taxonomy databases, including dedicated functions that streamline creating databases from popular sources, and functions for evaluating, comparing, and interactively exploring qualitative and quantitative characteristics across reference databases. To highlight the breadth and capabilities of RESCRIPt, we provide several examples for working with popular databases for microbiome profiling (SILVA, Greengenes, NCBI-RefSeq, GTDB), eDNA and diet metabarcoding surveys (BOLD, GenBank), as well as for genome comparison. We show that bigger is not always better, and reference databases with standardized taxonomies and those that focus on type strains have quantitative advantages, though may not be appropriate for all use cases. Most databases appear to benefit from some curation (quality filtering), though sequence clustering appears detrimental to database quality. Finally, we demonstrate the breadth and extensibility of RESCRIPt for reproducible workflows with a comparison of global hepatitis genomes. RESCRIPt provides tools to democratize the process of reference database acquisition and management, enabling researchers to reproducibly and transparently create reference materials for diverse research applications. RESCRIPt is released under a permissive BSD-3 license at https://github.com/bokulich-lab/RESCRIPt.

Generating and managing sequence and taxonomy reference data presents a bottleneck to many researchers, whether they are generating custom databases or attempting to format existing, curated reference databases for use with standard sequence analysis tools. Evaluating database quality and choosing the “best” database can be an equally formidable challenge. We developed RESCRIPt to alleviate this bottleneck, supporting reproducible, streamlined generation, curation, and evaluation of reference sequence databases. RESCRIPt uses QIIME 2 artifact file formats, which store all processing steps as data provenance within each file, allowing researchers to retrace the computational steps used to generate any given file. We used RESCRIPt to benchmark several commonly used marker-gene sequence databases for 16S rRNA genes, ITS, and COI sequences, demonstrating both the utility of RESCRIPt to streamline use of these databases, but also to evaluate several qualitative and quantitative characteristics of each database. We show that larger databases are not always best, and curation steps to reduce redundancy and filter out noisy sequences may be beneficial for some applications. We anticipate that RESCRIPt will streamline the use, management, and evaluation/selection of reference database materials for microbiomics, diet metabarcoding, eDNA, and other diverse applications.

Environmental DNA (eDNA) Metabarcoding in the Fish Market and Nearby Seafood Restaurants in Taiwan Reveals the Underestimation of Fish Species Diversity in Seafood

Seafood, especially the traditional one in Taiwan, is rarely sourced from a fixed species and routinely from similar species depending on their availability. Hence, the species composition of seafood can be complicated. While a DNA-based approach has been routinely utilized for species identification, a large scale of seafood identification in fish markets and restaurants could be challenging (e.g., elevated cost and time-consuming only for a limited number of species identification). In the present study, we aimed to identify the majority of fish species potentially consumed in fish markets and nearby seafood restaurants using environmental DNA (eDNA) metabarcoding. Four eDNA samplings from a local fish market and nearby seafood restaurants were conducted using Sterivex cartridges. Nineteen universal primers previously validated for fish species identification were utilized to amplify the fragments of mitochondrial DNA (12S, COI, ND5) of species in eDNA samples and sequenced with NovaSeq 6000 sequencing. A total of 153 fish species have been identified based on 417 fish related operational taxonomic units (OTUs) generated from 50,534,995 reads. Principal Coordinate Analysis (PCoA) further showed the differences in fish species between the sampling times and sampling sites. Of these fish species, 22 chondrichthyan fish, 14 Anguilliformes species, and 15 Serranidae species were respectively associated with smoked sharks, braised moray eels, and grouper fish soups. To our best knowledge, this work represents the first study to demonstrate the feasibility of a large scale of seafood identification using eDNA metabarcoding approach. Our findings also imply the species diversity in traditional seafood might be seriously underestimated and crucial for the conservation and management of marine resources.

High-Throughput Sequencing of Environmental DNA as a Tool for Monitoring Eukaryotic Communities and Potential Pathogens in a Coastal Upwelling Ecosystem

The marine environment includes diverse microeukaryotic organisms that play important functional roles in the ecosystem. With molecular approaches, eukaryotic taxonomy has been improved, complementing classical analysis. In this study, DNA metabarcoding was performed to describe putative pathogenic eukaryotic microorganisms in sediment and marine water fractions collected in Galicia (NW Spain) from 2016 to 2018. The composition of eukaryotic communities was distinct between sediment and water fractions. Protists were the most diverse group, with the clade TSAR (Stramenopiles, Alveolata, Rhizaria, and Telonemida) as the primary representative organisms in the environment. Harmful algae and invasive species were frequently detected. Potential pathogens, invasive pathogenic organisms as well as the causative agents of harmful phytoplanktonic blooms were identified in this marine ecosystem. Most of the identified pathogens have a crucial impact on the aquacultural sector or affect to relevant species in the marine ecosystem, such as diatoms. Moreover, pathogens with medical and veterinary importance worldwide were also found, as well as pathogens that affect diatoms. The evaluation of the health of a marine ecosystem that directly affects the aquacultural sector with a zoonotic concern was performed with the metabarcoding assay.

The Dark mAtteR iNvestigator (DARN) tool: getting to know the known unknowns in COI amplicon data

The mitochondrial cytochrome C oxidase subunit I gene (COI) is commonly used in environmental DNA (eDNA) metabarcoding studies, especially for assessing metazoan diversity. Yet, a great number of COI operational taxonomic units (OTUs) or/and amplicon sequence variants (ASVs) retrieved from such studies do not get a taxonomic assignment with a reference sequence. To assess and investigate such sequences, we have developed the Dark mAtteR iNvestigator (DARN) software tool. For this purpose, a reference COI-oriented phylogenetic tree was built from 1,593 consensus sequences covering all the three domains of life. With respect to eukaryotes, consensus sequences at the family level were constructed from 183,330 sequences retrieved from the Midori reference 2 database, which represented 70% of the initial number of reference sequences. Similarly, sequences from 431 bacterial and 15 archaeal taxa at the family level (29% and 1% of the initial number of reference sequences respectively) were retrieved from the BOLD and the PFam databases. DARN makes use of this phylogenetic tree to investigate COI pre-processed sequences of amplicon samples to provide both a tabular and a graphical overview of their phylogenetic assignments. To evaluate DARN, both environmental and bulk metabarcoding samples from different aquatic environments using various primer sets were analysed. We demonstrate that a large proportion of non-target prokaryotic organisms, such as bacteria and archaea, are also amplified in eDNA samples and we suggest prokaryotic COI sequences to be included in the reference databases used for the taxonomy assignment to allow for further analyses of dark matter. DARN source code is available on GitHub at https://github.com/hariszaf/darn and as a Docker image at https://hub.docker.com/r/hariszaf/darn.

A unified model of species abundance, genetic diversity, and functional diversity reveals the mechanisms structuring ecological communities

Biodiversity accumulates hierarchically by means of ecological and evolutionary processes and feedbacks. Within ecological communities drift, dispersal, speciation, and selection operate simultaneously to shape patterns of biodiversity. Reconciling the relative importance of these is hindered by current models and inference methods, which tend to focus on a subset of processes and their resulting predictions. Here we introduce massive ecoevolutionary synthesis simulations (MESS), a unified mechanistic model of community assembly, rooted in classic island biogeography theory, which makes temporally explicit joint predictions across three biodiversity data axes: (i) species richness and abundances, (ii) population genetic diversities, and (iii) trait variation in a phylogenetic context. Using simulations we demonstrate that each data axis captures information at different timescales, and that integrating these axes enables discriminating among previously unidentifiable community assembly models. MESS is unique in generating predictions of community-scale genetic diversity, and in characterizing joint patterns of genetic diversity, abundance, and trait values. MESS unlocks the full potential for investigation of biodiversity processes using multidimensional community data including a genetic component, such as might be produced by contemporary eDNA or metabarcoding studies. We combine MESS with supervised machine learning to fit the parameters of the model to real data and infer processes underlying how biodiversity accumulates, using communities of tropical trees, arthropods, and gastropods as case studies that span a range of data availability scenarios, and spatial and taxonomic scales.

Environmental DNA sampling protocols for the surveillance of marine non-indigenous species in Irish coastal waters

Understanding the spread and distribution of Non-Indigenous Species (NIS) is key when implementing legislation to maintain good ecosystem health. Environmental DNA (eDNA) has shown great potential to detect aquatic organisms in a rapid and cost-effective way, however their applicability to new environments must be validated prior to their implementation. Here, we tested different field sampling methods in combination with eDNA metabarcoding to develop a tool to detect NIS. Large and small volumes of seawater were filtered, in addition to the collection of sediment and horizontal tow net samples at 12 locations across four distinct geographic areas in Ireland. The biggest dissimilarity in the species recovered was found between sediment and town net samples. Tow nets showed to be the most efficient. A total of 357 taxa were identified, including 16 NIS. Fine mesh tow nets were identified as the most cost-efficient for large-scale monitoring and surveillance of NIS.

Defining gut mycobiota for wild animals: a need for caution in assigning authentic resident fungal taxa

Animal gut mycobiota, the community of fungi that reside within the gastrointestinal tract, make an important contribution to host health. Accordingly, there is an emerging interest to quantify the gut mycobiota of wild animals. However, many studies of wild animal gut mycobiota do not distinguish between the fungi that likely can reside within animal gastrointestinal tracts from the fungal taxa that are non-residents, such as macrofungi, lichens or plant symbionts/pathogens that can be ingested as part of the host's diet. Confounding the non-resident and resident gut fungi may obscure attempts to identify processes associated with the authentic, resident gut mycobiota per se. To redress this problem, we propose some strategies to filter the taxa identified within an apparent gut mycobiota based on an assessment of host ecology and fungal traits. Consideration of the different sources and roles of fungi present within the gastrointestinal tract should facilitate a more precise understanding of the causes and consequences of variation in wild animal gut mycobiota composition.

Going out for dinner-The consumption of agriculture pests by bats in urban areas

Insectivorous bats provide ecosystem services in agricultural and urban landscapes by consuming arthropods that are considered pests. Bat species inhabiting cities are expected to consume insects associated with urban areas, such as mosquitoes, flying termites, moths, and beetles. We captured insectivorous bats in the Federal District of Brazil and used fecal DNA metabarcoding to investigate the arthropod consumed by five bat species living in colonies in city buildings, and ascertained whether their predation was related to ecosystem services. These insectivorous bat species were found to consume 83 morphospecies of arthropods and among these 41 were identified to species, most of which were agricultural pests. We propose that bats may roost in the city areas and forage in the nearby agricultural fields using their ability to fly over long distances. We also calculated the value of the pest suppression ecosystem service by the bats. By a conservative estimation, bats save US$ 94 per hectare of cornfields, accounting for an annual savings of US$ 390.6 million per harvest in Brazil. Our study confirms that, regardless of their roosting location, bats are essential for providing ecosystem services in the cities, with extensive impacts on crops and elsewhere, in addition to significant savings in the use of pesticides.