Anacapa Toolkit : An environmental DNA toolkit for processing multilocus metabarcode datasets

Reference Sequence Browser: An R application with a user-friendly GUI to rapidly query sequence databases

Land managers, researchers, and regulators increasingly utilize environmental DNA (eDNA) techniques to monitor species richness, presence, and absence. In order to properly develop a biological assay for eDNA metabarcoding or quantitative PCR, scientists must be able to find not only reference sequences (previously identified sequences in a genomics database) that match their target taxa but also reference sequences that match non-target taxa. Determining which taxa have publicly available sequences in a time-efficient and accurate manner currently requires computational skills to search, manipulate, and parse multiple unconnected DNA sequence databases. Our team iteratively designed a Graphic User Interface (GUI) Shiny application called the Reference Sequence Browser (RSB) that provides users efficient and intuitive access to multiple genetic databases regardless of computer programming expertise. The application returns the number of publicly accessible barcode markers per organism in the NCBI Nucleotide, BOLD, or CALeDNA CRUX Metabarcoding Reference Databases. Depending on the database, we offer various search filters such as min and max sequence length or country of origin. Users can then download the FASTA/GenBank files from the RSB web tool, view statistics about the data, and explore results to determine details about the availability or absence of reference sequences.

A rapid phylogeny-based method for accurate community profiling of large-scale metabarcoding datasets

Environmental DNA (eDNA) is becoming an increasingly important tool in diverse scientific fields from ecological biomonitoring to wastewater surveillance of viruses. The fundamental challenge in eDNA analyses has been the bioinformatical assignment of reads to taxonomic groups. It has long been known that full probabilistic methods for phylogenetic assignment are preferable, but unfortunately, such methods are computationally intensive and are typically inapplicable to modern next-generation sequencing data. We present a fast approximate likelihood method for phylogenetic assignment of DNA sequences. Applying the new method to several mock communities and simulated datasets, we show that it identifies more reads at both high and low taxonomic levels more accurately than other leading methods. The advantage of the method is particularly apparent in the presence of polymorphisms and/or sequencing errors and when the true species is not represented in the reference database.

Environmental DNA reveals temporal variation in mesophotic reefs of the Humboldt upwelling ecosystems of central Chile: Toward a baseline for biodiversity monitoring of unexplored marine habitats

Temperate mesophotic reef ecosystems (TMREs) are among the least known marine habitats. Information on their diversity and ecology is geographically and temporally scarce, especially in highly productive large upwelling ecosystems. Lack of information remains an obstacle to understanding the importance of TMREs as habitats, biodiversity reservoirs and their connections with better-studied shallow reefs. Here, we use environmental DNA (eDNA) from water samples to characterize the community composition of TMREs on the central Chilean coast, generating the first baseline for monitoring the biodiversity of these habitats. We analyzed samples from two depths (30 and 60 m) over four seasons (spring, summer, autumn, and winter) and at two locations approximately 16 km apart. We used a panel of three metabarcodes, two that target all eukaryotes (18S rRNA and mitochondrial COI) and one specifically targeting fishes (16S rRNA). All panels combined encompassed eDNA assigned to 42 phyla, 90 classes, 237 orders, and 402 families. The highest family richness was found for the phyla Arthropoda, Bacillariophyta, and Chordata. Overall, family richness was similar between depths but decreased during summer, a pattern consistent at both locations. Our results indicate that the structure (composition) of the mesophotic communities varied predominantly with seasons. We analyzed further the better-resolved fish assemblage and compared eDNA with other visual methods at the same locations and depths. We recovered eDNA from 19 genera of fish, six of these have also been observed on towed underwater videos, while 13 were unique to eDNA. We discuss the potential drivers of seasonal differences in community composition and richness. Our results suggest that eDNA can provide valuable insights for monitoring TMRE communities but highlight the necessity of completing reference DNA databases available for this region.

Semi-automated sequence curation for reliable reference datasets in ITS2 vascular plant DNA (meta-)barcoding

One of the most critical steps for accurate taxonomic identification in DNA (meta)-barcoding is to have an accurate DNA reference sequence dataset for the marker of choice. Therefore, developing such a dataset has been a long-term ambition, especially in the Viridiplantae kingdom. Typically, reference datasets are constructed with sequences downloaded from general public databases, which can carry taxonomic and other relevant errors. Herein, we constructed a curated (i) global dataset, (ii) European crop dataset, and (iii) 27 datasets for the EU countries for the ITS2 barcoding marker of vascular plants. To that end, we first developed a pipeline script that entails (i) an automated curation stage comprising five filters, (ii) manual taxonomic correction for misclassified taxa, and (iii) manual addition of newly sequenced species. The pipeline allows easy updating of the curated datasets. With this approach, 13% of the sequences, corresponding to 7% of species originally imported from GenBank, were discarded. Further, 259 sequences were manually added to the curated global dataset, which now comprises 307,977 sequences of 111,382 plant species.

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

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

Principles for Constructing DNA Barcode Reference Libraries

All DNA barcode methods rely on reference sequences linked to well-curated voucher specimens. Definitions for and locations of DNA barcode reference libraries are not standardized, and vary throughout the literature. Standardizing, and centralizing reference specimens would provide an unambiguous source, analogous to reference genomes, to reproduce identifications and improve a library. This chapter proposes a working definition of a DNA barcode reference library, consistent with DNA barcode data standards, along with principles and methods to consider when producing or using such a library. These methods allow explicit traceback to sequence-sources which elevate the value of voucher specimens, and create a potential for community curation.

Environmental DNA reveals patterns of biological invasion in an inland sea

Non-native species have the potential to cause ecological and economic harm to coastal and estuarine ecosystems. Understanding which habitat types are most vulnerable to biological invasions, where invasions originate, and the vectors by which they arrive can help direct limited resources to prevent or mitigate ecological and socio-economic harm. Information about the occurrence of non-native species can help guide interventions at all stages of invasion, from first introduction, to naturalization and invasion. However, monitoring at relevant scales requires considerable investment of time, resources, and taxonomic expertise. Environmental DNA (eDNA) metabarcoding methods sample coastal ecosystems at broad spatial and temporal scales to augment established monitoring methods. We use COI mtDNA eDNA sampling to survey a diverse assemblage of species across distinct habitats in the Salish Sea in Washington State, USA, and classify each as non-native, native, or indeterminate in origin. The non-native species detected include both well-documented invaders and species not previously reported within the Salish Sea. We find a non-native assemblage dominated by shellfish and algae with native ranges in the temperate western Pacific, and find more-retentive estuarine habitats to be invaded at far higher levels than better-flushed rocky shores. Furthermore, we find an increase in invasion level with higher water temperatures in spring and summer across habitat types. This analysis contributes to a growing understanding of the biotic and abiotic factors that influence invasion level, and underscores the utility of eDNA surveys to monitor biological invasions and to better understand the factors that drive these invasions.

Delayed postglacial colonization of Betula in Iceland and the circum North Atlantic

As the Arctic continues to warm, woody shrubs are expected to expand northward. This process, known as 'shrubification,' has important implications for regional biodiversity, food web structure, and high-latitude temperature amplification. While the future rate of shrubification remains poorly constrained, past records of plant immigration to newly deglaciated landscapes in the Arctic may serve as useful analogs. We provide one new postglacial Holocene sedimentary ancient DNA (sedaDNA) record of vascular plants from Iceland and place a second Iceland postglacial sedaDNA record on an improved geochronology; both show Salicaceae present shortly after deglaciation, whereas Betulaceae first appears more than 1000 y later. We find a similar pattern of delayed Betulaceae colonization in eight previously published postglacial sedaDNA records from across the glaciated circum North Atlantic. In nearly all cases, we find that Salicaceae colonizes earlier than Betulaceae and that Betulaceae colonization is increasingly delayed for locations farther from glacial-age woody plant refugia. These trends in Salicaceae and Betulaceae colonization are consistent with the plant families' environmental tolerances, species diversity, reproductive strategies, seed sizes, and soil preferences. As these reconstructions capture the efficiency of postglacial vascular plant migration during a past period of high-latitude warming, a similarly slow response of some woody shrubs to current warming in glaciated regions, and possibly non-glaciated tundra, may delay Arctic shrubification and future changes in the structure of tundra ecosystems and temperature amplification.

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

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

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

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

Soil Microbial Community Composition and Tolerance to Contaminants in an Urban Brownfield Site

Brownfields are unused sites that contain hazardous substances due to previous commercial or industrial use. The sites are inhospitable for many organisms, but some fungi and microbes can tolerate and thrive in the nutrient-depleted and contaminated soils. However, few studies have characterized the impacts of long-term contamination on soil microbiome composition and diversity at brownfields. This study focuses on an urban brownfield-a former rail yard in Los Angeles that is contaminated with heavy metals, volatile organic compounds, and petroleum-derived pollutants. We anticipate that heavy metals and organic pollutants will shape soil microbiome diversity and that several candidate fungi and bacteria will be tolerant to the contaminants. We sequence three gene markers (16S ribosomal RNA, 18S ribosomal RNA, and the fungal internal transcribed spacer (FITS)) in 55 soil samples collected at five depths to (1) profile the composition of the soil microbiome across depths; (2) determine the extent to which hazardous chemicals predict microbiome variation; and (3) identify microbial taxonomic groups that may metabolize these contaminants. Detected contaminants in the samples included heavy metals, petroleum hydrocarbons, polycyclic aromatic hydrocarbons, and volatile organic compounds. Bacterial, eukaryotic, and fungal communities all varied with depth and with concentrations of arsenic, chromium, cobalt, and lead. 18S rRNA microbiome richness and fungal richness were positively correlated with lead and cobalt levels, respectively. Furthermore, bacterial Paenibacillus and Iamia, eukaryotic Actinochloris, and fungal Alternaria were enriched in contaminated soils compared to uncontaminated soils and represent taxa of interest for future bioremediation research. Based on our results, we recommend incorporating DNA-based multi-marker microbial community profiling at multiple sites and depths in brownfield site assessment standard methods and restoration.

Signal and noise in metabarcoding data

Metabarcoding is a powerful molecular tool for simultaneously surveying hundreds to thousands of species from a single sample, underpinning microbiome and environmental DNA (eDNA) methods. Deriving quantitative estimates of underlying biological communities from metabarcoding is critical for enhancing the utility of such approaches for health and conservation. Recent work has demonstrated that correcting for amplification biases in genetic metabarcoding data can yield quantitative estimates of template DNA concentrations. However, a major source of uncertainty in metabarcoding data stems from non-detections across technical PCR replicates where one replicate fails to detect a species observed in other replicates. Such non-detections are a special case of variability among technical replicates in metabarcoding data. While many sampling and amplification processes underlie observed variation in metabarcoding data, understanding the causes of non-detections is an important step in distinguishing signal from noise in metabarcoding studies. Here, we use both simulated and empirical data to 1) suggest how non-detections may arise in metabarcoding data, 2) outline steps to recognize uninformative data in practice, and 3) identify the conditions under which amplicon sequence data can reliably detect underlying biological signals. We show with both simulations and empirical data that, for a given species, the rate of non-detections among technical replicates is a function of both the template DNA concentration and species-specific amplification efficiency. Consequently, we conclude metabarcoding datasets are strongly affected by (1) deterministic amplification biases during PCR and (2) stochastic sampling of amplicons during sequencing-both of which we can model-but also by (3) stochastic sampling of rare molecules prior to PCR, which remains a frontier for quantitative metabarcoding. Our results highlight the importance of estimating species-specific amplification efficiencies and critically evaluating patterns of non-detection in metabarcoding datasets to better distinguish environmental signal from the noise inherent in molecular detections of rare targets.

A detailed workflow to develop QIIME2-formatted reference databases for taxonomic analysis of DNA metabarcoding data

Background

The DNA metabarcoding approach has become one of the most used techniques to study the taxa composition of various sample types. To deal with the high amount of data generated by the high-throughput sequencing process, a bioinformatics workflow is required and the QIIME2 platform has emerged as one of the most reliable and commonly used. However, only some pre-formatted reference databases dedicated to a few barcode sequences are available to assign taxonomy. If users want to develop a new custom reference database, several bottlenecks still need to be addressed and a detailed procedure explaining how to develop and format such a database is currently missing. In consequence, this work is aimed at presenting a detailed workflow explaining from start to finish how to develop such a curated reference database for any barcode sequence.

Results

We developed DB4Q2, a detailed workflow that allowed development of plant reference databases dedicated to ITS2 and rbcL, two commonly used barcode sequences in plant metabarcoding studies. This workflow addresses several of the main bottlenecks connected with the development of a curated reference database. The detailed and commented structure of DB4Q2 offers the possibility of developing reference databases even without extensive bioinformatics skills, and avoids ‘black box’ systems that are sometimes encountered. Some filtering steps have been included to discard presumably fungal and misidentified sequences. The flexible character of DB4Q2 allows several key sequence processing steps to be included or not, and downloading issues can be avoided. Benchmarking the databases developed using DB4Q2 revealed that they performed well compared to previously published reference datasets.

Conclusion

This study presents DB4Q2, a detailed procedure to develop custom reference databases in order to carry out taxonomic analyses with QIIME2, but also with other bioinformatics platforms if desired. This work also provides ready-to-use plant ITS2 and rbcL databases for which the prediction accuracy has been assessed and compared to that of other published databases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01067-5.

A manager's guide to using eDNA metabarcoding in marine ecosystems

Environmental DNA (eDNA) metabarcoding is a powerful tool that can enhance marine ecosystem/biodiversity monitoring programs. Here we outline five important steps managers and researchers should consider when developing eDNA monitoring program: (1) select genes and primers to target taxa; (2) assemble or develop comprehensive barcode reference databases; (3) apply rigorous site occupancy based decontamination pipelines; (4) conduct pilot studies to define spatial and temporal variance of eDNA; and (5) archive samples, extracts, and raw sequence data. We demonstrate the importance of each of these considerations using a case study of eDNA metabarcoding in the Ports of Los Angeles and Long Beach. eDNA metabarcoding approaches detected 94.1% (16/17) of species observed in paired trawl surveys while identifying an additional 55 native fishes, providing more comprehensive biodiversity inventories. Rigorous benchmarking of eDNA metabarcoding results improved ecological interpretation and confidence in species detections while providing archived genetic resources for future analyses. Well designed and validated eDNA metabarcoding approaches are ideally suited for biomonitoring applications that rely on the detection of species, including mapping invasive species fronts and endangered species habitats as well as tracking range shifts in response to climate change. Incorporating these considerations will enhance the utility and efficacy of eDNA metabarcoding for routine biomonitoring applications.

The Functional Biogeography of eDNA Metacommunities in the Post-Fire Landscape of the Angeles National Forest

Wildfires have continued to increase in frequency and severity in Southern California due in part to climate change. To gain a further understanding of microbial soil communities’ response to fire and functions that may enhance post-wildfire resilience, soil fungal and bacterial microbiomes were studied from different wildfire areas in the Gold Creek Preserve within the Angeles National Forest using 16S, FITS, 18S, 12S, PITS, and COI amplicon sequencing. Sequencing datasets from December 2020 and June 2021 samplings were analyzed using QIIME2, ranacapa, stats, vcd, EZBioCloud, and mixomics. Significant differences were found among bacterial and fungal taxa associated with different fire areas in the Gold Creek Preserve. There was evidence of seasonal shifts in the alpha diversity of the bacterial communities. In the sparse partial least squares analysis, there were strong associations (r > 0.8) between longitude, elevation, and a defined cluster of Amplicon Sequence Variants (ASVs). The Chi-square test revealed differences in fungi−bacteria (F:B) proportions between different trails (p = 2 × 10−16). sPLS results focused on a cluster of Green Trail samples with high elevation and longitude. Analysis revealed the cluster included the post-fire pioneer fungi Pyronema and Tremella. Chlorellales algae and possibly pathogenic Fusarium sequences were elevated. Bacterivorous Corallococcus, which secretes antimicrobials, and bacterivorous flagellate Spumella were associated with the cluster. There was functional redundancy in clusters that were differently composed but shared similar ecological functions. These results implied a set of traits for post-fire resiliency. These included photo-autotrophy, mineralization of pyrolyzed organic matter and aromatic/oily compounds, potential pathogenicity and parasitism, antimicrobials, and N-metabolism.

AncestralClust: clustering of divergent nucleotide sequences by ancestral sequence reconstruction using phylogenetic trees

MOTIVATION

Clustering is a fundamental task in the analysis of nucleotide sequences. Despite the exponential increase in the size of sequence databases of homologous genes, few methods exist to cluster divergent sequences. Traditional clustering methods have mostly focused on optimizing high speed clustering of highly similar sequences. We develop a phylogenetic clustering method which infers ancestral sequences for a set of initial clusters and then uses a greedy algorithm to cluster sequences.

RESULTS

We describe a clustering program AncestralClust, which is developed for clustering divergent sequences. We compare this method with other state-of-the-art clustering methods using datasets of homologous sequences from different species. We show that, in divergent datasets, AncestralClust has higher accuracy and more even cluster sizes than current popular methods.

AVAILABILITY AND IMPLEMENTATION

AncestralClust is an Open Source program available at https://github.com/lpipes/ancestralclust.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Revisiting the effect of PCR replication and sequencing depth on biodiversity metrics in environmental DNA metabarcoding

Environmental DNA (eDNA) metabarcoding is an increasingly popular tool for measuring and cataloguing biodiversity. Because the environments and substrates in which DNA is preserved differ considerably, eDNA research often requires bespoke approaches to generating eDNA data. Here, we explore how two experimental choices in eDNA study design-the number of PCR replicates and the depth of sequencing of PCR replicates-influence the composition and consistency of taxa recovered from eDNA extracts. We perform 24 PCR replicates from each of six soil samples using two of the most common metabarcodes for Fungi and Viridiplantae (ITS1 and ITS2), and sequence each replicate to an average depth of ~84,000 reads. We find that PCR replicates are broadly consistent in composition and relative abundance of dominant taxa, but that low abundance taxa are often unique to one or a few PCR replicates. Taxa observed in one out of 24 PCR replicates make up 21-29% of the total taxa detected. We also observe that sequencing depth or rarefaction influences alpha diversity and beta diversity estimates. Read sampling depth influences local contribution to beta diversity, placement in ordinations, and beta dispersion in ordinations. Our results suggest that, because common taxa drive some alpha diversity estimates, few PCR replicates and low read sampling depths may be sufficient for many biological applications of eDNA metabarcoding. However, because rare taxa are recovered stochastically, eDNA metabarcoding may never fully recover the true amplifiable alpha diversity in an eDNA extract. Rare taxa drive PCR replicate outliers of alpha and beta diversity and lead to dispersion differences at different read sampling depths. We conclude that researchers should consider the complexity and unevenness of a community when choosing analytical approaches, read sampling depths, and filtering thresholds to arrive at stable estimates.

Meta-Fish-Lib: A generalised, dynamic DNA reference library pipeline for metabarcoding of fishes

The accuracy and reliability of DNA metabarcoding analyses depend on the breadth and quality of the reference libraries that underpin them. However, there are limited options available to obtain and curate the huge volumes of sequence data that are available on public repositories such as NCBI and BOLD. Here, we provide a pipeline to download, clean and annotate mitochondrial DNA sequence data for a given list of fish species. Features of this pipeline include (a) support for multiple metabarcode markers; (b) searches on species synonyms and taxonomic name validation; (c) phylogeny assisted quality control for identification and removal of misannotated sequences; (d) automatically generated coverage reports for each new GenBank release update; and (e) citable, versioned DOIs. As an example we provide a ready-to-use curated reference library for the marine and freshwater fishes of the U.K. To augment this reference library for environmental DNA metabarcoding specifically, we generated 241 new MiFish-12S sequences for 88 U.K. marine species, and make available new primer sets useful for sequencing these. This brings the coverage of common U.K. species for the MiFish-12S fragment to 93%, opening new avenues for scaling up fish metabarcoding across wide spatial gradients. The Meta-Fish-Lib reference library and pipeline is hosted at https://github.com/genner-lab/meta-fish-lib.

Landscape analyses using eDNA metabarcoding and Earth observation predict community biodiversity in California

Ecosystems globally are under threat from ongoing anthropogenic environmental change. Effective conservation management requires more thorough biodiversity surveys that can reveal system-level patterns and that can be applied rapidly across space and time. Using modern ecological models and community science, we integrate environmental DNA and Earth observations to produce a time snapshot of regional biodiversity patterns and provide multi-scalar community-level characterization. We collected 278 samples in spring 2017 from coastal, shrub, and lowland forest sites in California, a complex ecosystem and biodiversity hotspot. We recovered 16,118 taxonomic entries from eDNA analyses and compiled associated traditional observations and environmental data to assess how well they predicted alpha, beta, and zeta diversity. We found that local habitat classification was diagnostic of community composition and distinct communities and organisms in different kingdoms are predicted by different environmental variables. Nonetheless, gradient forest models of 915 families recovered by eDNA analysis and using BIOCLIM variables, Sentinel-2 satellite data, human impact, and topographical features as predictors, explained 35% of the variance in community turnover. Elevation, sand percentage, and photosynthetic activities (NDVI32) were the top predictors. In addition to this signal of environmental filtering, we found a positive relationship between environmentally predicted families and their numbers of biotic interactions, suggesting environmental change could have a disproportionate effect on community networks. Together, these analyses show that coupling eDNA with environmental predictors including remote sensing data has capacity to test proposed Essential Biodiversity Variables and create new landscape biodiversity baselines that span the tree of life.

Improving metabarcoding taxonomic assignment: A case study of fishes in a large marine ecosystem

DNA metabarcoding is an important tool for molecular ecology. However, its effectiveness hinges on the quality of reference sequence databases and classification parameters employed. Here we evaluate the performance of MiFish 12S taxonomic assignments using a case study of California Current Large Marine Ecosystem fishes to determine best practices for metabarcoding. Specifically, we use a taxonomy cross-validation by identity framework to compare classification performance between a global database comprised of all available sequences and a curated database that only includes sequences of fishes from the California Current Large Marine Ecosystem. We demonstrate that the regional database provides higher assignment accuracy than the comprehensive global database. We also document a tradeoff between accuracy and misclassification across a range of taxonomic cutoff scores, highlighting the importance of parameter selection for taxonomic classification. Furthermore, we compared assignment accuracy with and without the inclusion of additionally generated reference sequences. To this end, we sequenced tissue from 597 species using the MiFish 12S primers, adding 252 species to GenBank's existing 550 California Current Large Marine Ecosystem fish sequences. We then compared species and reads identified from seawater environmental DNA samples using global databases with and without our generated references, and the regional database. The addition of new references allowed for the identification of 16 additional native taxa representing 17.0% of total reads from eDNA samples, including species with vast ecological and economic value. Together these results demonstrate the importance of comprehensive and curated reference databases for effective metabarcoding and the need for locus-specific validation efforts.

Characterizing Industrial and Artisanal Fishing Vessel Catch Composition Using Environmental DNA and Satellite-Based Tracking Data

The decline in wild-caught fisheries paired with increasing global seafood demand is pushing the need for seafood sustainability to the forefront of national and regional priorities. Validation of species identity is a crucial early step, yet conventional monitoring and surveillance tools are limited in their effectiveness because they are extremely time-consuming and require expertise in fish identification. DNA barcoding methods are a versatile tool for the genetic monitoring of wildlife products; however, they are also limited by requiring individual tissue samples from target specimens which may not always be possible given the speed and scale of seafood operations. To circumvent the need to individually sample organisms, we pilot an approach that uses forensic environmental DNA (eDNA) metabarcoding to profile fish species composition from the meltwater in fish holds on industrial and artisanal fishing vessels in Ecuador. Fish identified genetically as present were compared to target species reported by each vessel’s crew. Additionally, we contrasted the geographic range of identified species against the satellite-based fishing route data of industrial vessels to determine if identified species could be reasonably expected in the catch.

Short-lived detection of an introduced vertebrate eDNA signal in a nearshore rocky reef environment

Environmental DNA (eDNA) is increasingly used to measure biodiversity of marine ecosystems, yet key aspects of the temporal dynamics of eDNA remain unknown. Of particular interest is in situ persistence of eDNA signals in dynamic marine environments, as eDNA degradation rates have predominantly been quantified through mesocosm studies. To determine in situ eDNA residence times, we introduced an eDNA signal from a non-native fish into a protected bay of a Southern California rocky reef ecosystem, and then measured changes in both introduced and background eDNA signals across a fixed transect over 96 hours. Foreign eDNA signal was no longer detected only 7.5 hours after introduction, a time substantially shorter than the multi-day persistence times in laboratory studies. Moreover, the foreign eDNA signal spread along the entire 38 m transect within 1.5 hours after introduction, indicating that transport and diffusion play a role in eDNA detectability even in protected low energy marine environments. Similarly, native vertebrate eDNA signals varied greatly over the 96 hours of observation as well as within two additional nearby fixed transects sampled over 120 hours. While community structure did significantly change across time of day and tidal direction, neither accounted for the majority of observed variation. Combined, results show that both foreign and native eDNA signatures can exhibit substantial temporal heterogeneity, even on hourly time scales. Further work exploring eDNA decay from lagrangian perspective and quantifying effects of sample and technical replication are needed to better understand temporal variation of eDNA signatures in nearshore marine environments.

Benchmarking bioinformatic tools for fast and accurate eDNA metabarcoding species identification

Bioinformatic analysis of eDNA metabarcoding data is a crucial step toward rigorously assessing biodiversity. Many programs are now available for each step of the required analyses, but their relative abilities at providing fast and accurate species lists have seldom been evaluated. We used simulated mock communities and real fish eDNA metabarcoding data to evaluate the performance of 13 bioinformatic programs and pipelines to retrieve fish occurrence and read abundance using the 12S mt rRNA gene marker. We used four indices to compare the outputs of each program with the simulated samples: sensitivity, F-measure, root-mean-square error (RMSE) on read relative abundances, and execution time. We found marked differences among programs only for the taxonomic assignment step, both in terms of sensitivity, F-measure and RMSE. Running time was highly different between programs for each step. The fastest programs with best indices for each step were assembled into a pipeline. We compared this pipeline to pipelines constructed from existing toolboxes (OBITools, Barque, and QIIME 2). Our pipeline and Barque obtained the best performance for all indices and appear to be better alternatives to highly used pipelines for analysing fish eDNA metabarcoding data when a complete reference database is available. Analysis on real eDNA metabarcoding data also indicated differences for taxonomic assignment and execution time only. This study reveals major differences between programs during the taxonomic assignment step. The choice of algorithm for the taxonomic assignment can have a significant impact on diversity estimates and should be made according to the objectives of the study.

Environmental DNA Metabarcoding for Simultaneous Monitoring and Ecological Assessment of Many Harmful Algae

Harmful algae can have profound economic, environmental, and social consequences. As the timing, frequency, and severity of harmful algal blooms (HABs) change alongside global climate, efficient tools to monitor and understand the current ecological context of these taxa are increasingly important. Here we employ environmental DNA metabarcoding to identify patterns in a wide variety of potentially harmful algae and associated ecological communities in the Hood Canal of Puget Sound in Washington State, USA. Tracking trends of occurrence in a series of water samples over a period of 19 months, we find algal sequences from genera with harmful members in a majority of samples, suggesting that these groups are routinely present in local waters. We report patterns in variants of the economically important genus Pseudo-nitzschia (of which some members produce domoic acid; family Bacillariaceae), as well as multiple potentially harmful algal taxa previously unknown or poorly documented in the region, including a cold-water variant from the genus Alexandrium (of which some members produce saxitoxin; family Gonyaulacaceae), two variants from the genus Karlodinium (of which some members produce karlotoxins; family Kareniaceae), and one variant from the parasitic genus Hematodinium (family Syndiniaceae). We then use data on environmental variables and the biological community surrounding each algal taxon to illustrate the ecological context in which they are commonly found. Environmental DNA metabarcoding thus simultaneously (1) alerts us to potential new or cryptic occurrences of algae from harmful genera, (2) expands our knowledge of the co-occurring conditions and species associated with the growth of these organisms in changing marine environments, and (3) suggests a pathway for multispecies monitoring and management moving forward.

eDNAFlow, an automated, reproducible and scalable workflow for analysis of environmental DNA sequences exploiting Nextflow and Singularity

Metabarcoding of environmental DNA (eDNA) when coupled with high throughput sequencing is revolutionising the way biodiversity can be monitored across a wide range of applications. However, the large number of tools deployed in downstream bioinformatic analyses often places a challenge in configuration and maintenance of a workflow, and consequently limits the research reproducibility. Furthermore, scalability needs to be considered to handle the growing amount of data due to increase in sequence output and the scale of project. Here, we describe eDNAFlow, a fully automated workflow that employs a number of state-of-the-art applications to process eDNA data from raw sequences (single-end or paired-end) to generation of curated and noncurated zero-radius operational taxonomic units (ZOTUs) and their abundance tables. This pipeline is based on Nextflow and Singularity which enable a scalable, portable and reproducible workflow using software containers on a local computer, clouds and high-performance computing (HPC) clusters. Finally, we present an in-house Python script to assign taxonomy to ZOTUs based on user specified thresholds for assigning lowest common ancestor (LCA). We demonstrate the utility and efficiency of the pipeline using an example of a published coral diversity biomonitoring study. Our results were congruent with the aforementioned study. The scalability of the pipeline is also demonstrated through analysis of a large data set containing 154 samples. To our knowledge, this is the first automated bioinformatic pipeline for eDNA analysis using two powerful tools: Nextflow and Singularity. This pipeline addresses two major challenges in the analysis of eDNA data; scalability and reproducibility.

eDNA metabarcoding as a biomonitoring tool for marine protected areas

Monitoring of marine protected areas (MPAs) is critical for marine ecosystem management, yet current protocols rely on SCUBA-based visual surveys that are costly and time consuming, limiting their scope and effectiveness. Environmental DNA (eDNA) metabarcoding is a promising alternative for marine ecosystem monitoring, but more direct comparisons to visual surveys are needed to understand the strengths and limitations of each approach. This study compares fish communities inside and outside the Scorpion State Marine Reserve off Santa Cruz Island, CA using eDNA metabarcoding and underwater visual census surveys. Results from eDNA captured 76% (19/25) of fish species and 95% (19/20) of fish genera observed during pairwise underwater visual census. Species missed by eDNA were due to the inability of MiFish 12S barcodes to differentiate species of rockfishes (Sebastes, n = 4) or low site occupancy rates of crevice-dwelling Lythrypnus gobies. However, eDNA detected an additional 23 fish species not recorded in paired visual surveys, but previously reported from prior visual surveys, highlighting the sensitivity of eDNA. Significant variation in eDNA signatures by location (50 m) and site (~1000 m) demonstrates the sensitivity of eDNA to address key questions such as community composition inside and outside MPAs. Results demonstrate the utility of eDNA metabarcoding for monitoring marine ecosystems, providing an important complementary tool to visual methods.

Lake Sedimentary DNA Research on Past Terrestrial and Aquatic Biodiversity: Overview and Recommendations

The use of lake sedimentary DNA to track the long-term changes in both terrestrial and aquatic biota is a rapidly advancing field in paleoecological research. Although largely applied nowadays, knowledge gaps remain in this field and there is therefore still research to be conducted to ensure the reliability of the sedimentary DNA signal. Building on the most recent literature and seven original case studies, we synthesize the state-of-the-art analytical procedures for effective sampling, extraction, amplification, quantification and/or generation of DNA inventories from sedimentary ancient DNA (sedaDNA) via high-throughput sequencing technologies. We provide recommendations based on current knowledge and best practises.

At Palmyra Atoll, the fish-community environmental DNA signal changes across habitats but not with tides

At Palmyra Atoll, the environmental DNA (eDNA) signal on tidal sand flats was associated with fish biomass density and captured 98%-100% of the expected species diversity there. Although eDNA spilled over across habitats, species associated with reef habitat contributed more eDNA to reef sites than to sand-flat sites, and species associated with sand-flat habitat contributed more eDNA to sand-flat sites than to reef sites. Tides did not disrupt the sand-flat habitat signal. At least 25 samples give a coverage >97.5% at this diverse, tropical, marine system.

Environmental DNA metabarcoding reveals winners and losers of global change in coastal waters

Studies of the ecological effects of global change often focus on one or a few species at a time. Consequently, we know relatively little about the changes underway at real-world scales of biological communities, which typically have hundreds or thousands of interacting species. Here, we use COI mtDNA amplicons from monthly samples of environmental DNA to survey 221 planktonic taxa along a gradient of temperature, salinity, dissolved oxygen and carbonate chemistry in nearshore marine habitat. The result is a high-resolution picture of changes in ecological communities using a technique replicable across a wide variety of ecosystems. We estimate community-level differences associated with time, space and environmental variables, and use these results to forecast near-term community changes due to warming and ocean acidification. We find distinct communities in warmer and more acidified conditions, with overall reduced richness in diatom assemblages and increased richness in dinoflagellates. Individual taxa finding more suitable habitat in near-future waters are more taxonomically varied and include the ubiquitous coccolithophore Emiliania huxleyi and the harmful dinoflagellate Alexandrium sp. These results suggest foundational changes for nearshore food webs under near-future conditions.

PUMAA: A Platform for Accessible Microbiome Analysis in the Undergraduate Classroom

Improvements in high-throughput sequencing makes targeted amplicon analysis an ideal method for the study of human and environmental microbiomes by undergraduates. Multiple bioinformatics programs are available to process and interpret raw microbial diversity datasets, and the choice of programs to use in curricula is largely determined by student learning goals. Many of the most commonly used microbiome bioinformatics platforms offer end-to-end data processing and data analysis using a command line interface (CLI), but the downside for novice microbiome researchers is the steep learning curve often required. Alternatively, some sequencing providers include processing of raw data and taxonomy assignments as part of their pipelines. This, when coupled with available web-based or graphical user interface (GUI) analysis and visualization tools, eliminates the need for students or instructors to have extensive CLI experience. However, lack of universal data formats can make integration of these tools challenging. For example, tools for upstream and downstream analyses frequently use multiple different data formats which then require writing custom scripts or hours of manual work to make the files compatible. Here, we describe a microbial ecology bioinformatics curriculum that focuses on data analysis, visualization, and statistical reasoning by taking advantage of existing web-based and GUI tools. We created the Program for Unifying Microbiome Analysis Applications (PUMAA), which solves the problem of inconsistent files by formatting the output files from several raw data processing programs to seamlessly transition to a suite of GUI programs for analysis and visualization of microbiome taxonomic and inferred functional profiles. Additionally, we created a series of tutorials to accompany each of the microbiome analysis curricular modules. From pre- and post-course surveys, students in this curriculum self-reported conceptual and confidence gains in bioinformatics and data analysis skills. Students also demonstrated gains in biologically relevant statistical reasoning based on rubric-guided evaluations of open-ended survey questions and the Statistical Reasoning in Biology Concept Inventory. The PUMAA program and associated analysis tutorials enable students and researchers with no computational experience to effectively analyze real microbiome datasets to investigate real-world research questions.

MARES, a replicable pipeline and curated reference database for marine eukaryote metabarcoding

The use of DNA metabarcoding to characterise the biodiversity of environmental and community samples has exploded in recent years. However, taxonomic inferences from these studies are contingent on the quality and completeness of the sequence reference database used to characterise sample species-composition. In response, studies often develop custom reference databases to improve species assignment. The disadvantage of this approach is that it limits the potential for database re-use, and the transferability of inferences across studies. Here, we present the MARine Eukaryote Species (MARES) reference database for use in marine metabarcoding studies, created using a transparent and reproducible pipeline. MARES includes all COI sequences available in GenBank and BOLD for marine taxa, unified into a single taxonomy. Our pipeline facilitates the curation of sequences, synonymization of taxonomic identifiers used by different repositories, and formatting these data for use in taxonomic assignment tools. Overall, MARES provides a benchmark COI reference database for marine eukaryotes, and a standardised pipeline for (re)producing reference databases enabling integration and fair comparison of marine DNA metabarcoding results.

Measurement(s)	DNA
Technology Type(s)	bioinformatics analysis
Factor Type(s)	DNA sequence
Sample Characteristic - Organism	Eukaryota
Sample Characteristic - Environment	marine environment

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12324122

Experimental DNA Demethylation Associates with Changes in Growth and Gene Expression of Oak Tree Seedlings

Epigenetic modifications such as DNA methylation, where methyl groups are added to cytosine base pairs, have the potential to impact phenotypic variation and gene expression, and could influence plant response to changing environments. One way to test this impact is through the application of chemical demethylation agents, such as 5-Azacytidine, which inhibit DNA methylation and lead to a partial reduction in DNA methylation across the genome. In this study, we treated 5-month-old seedlings of the tree, Quercus lobata, with foliar application of 5-Azacytidine to test whether a reduction in genome-wide methylation would cause differential gene expression and change phenotypic development. First, we demonstrate that demethylation treatment led to 3–6% absolute reductions and 6.7–43.2% relative reductions in genome-wide methylation across CG, CHG, and CHH sequence contexts, with CHH showing the strongest relative reduction. Seedlings treated with 5-Azacytidine showed a substantial reduction in new growth, which was less than half that of control seedlings. We tested whether this result could be due to impact of the treatment on the soil microbiome and found minimal differences in the soil microbiome between two groups, although with limited sample size. We found no significant differences in leaf fluctuating asymmetry (i.e., deviations from bilateral symmetry), which has been found in other studies. Nonetheless, treated seedlings showed differential expression of a total of 23 genes. Overall, this study provides initial evidence that DNA methylation is involved in gene expression and phenotypic variation in seedlings and suggests that removal of DNA methylation affects plant development.