Management of DNA reference libraries for barcoding and metabarcoding studies with the R package refdb

DNA barcoding and metabarcoding are revolutionizing the study and survey of biodiversity. In order to assign taxonomic labels to the DNA sequence data retrieved, these methods are strongly dependent on comprehensive and accurate reference databases. Producing reliable databases linking biological sequences and taxonomic data can be—and often has been—done using mainstream tools such as spreadsheet software. However, spreadsheets quickly become insufficient when the amount of data increases to thousands of taxa and sequences to be matched, and validation operations become more complex and are error prone if done in a manual way. Thus, there is a clear need for providing scientists with user-friendly, reliable and powerful tools to manipulate and manage DNA reference databases in tractable, sound and efficient ways. Here, we introduce the R package refdb as an environment for semi-automatic and assisted construction of DNA reference libraries. The refdb package is a reference database manager offering a set of powerful functions to import, organize, clean, filter, audit and export the data. It is broadly applicable in metabarcoding data generally obtained in biodiversity and biomonitoring studies. We present the main features of the package and outline how refdb can speed up reference database generation, management and handling, and thus contribute to standardization and repeatability in barcoding and metabarcoding studies.     

crabs—A software program to generate curated reference databases for metabarcoding sequencing data

The measurement of biodiversity is an integral aspect of life science research. With the establishment of second- and third-generation sequencing technologies, an increasing amount of metabarcoding data is being generated as we seek to describe the extent and patterns of biodiversity in multiple contexts. The reliability and accuracy of taxonomically assigning metabarcoding sequencing data have been shown to be critically influenced by the quality and completeness of reference databases. Custom, curated, eukaryotic reference databases, however, are scarce, as are the software programs for generating them. Here, we present crabs (Creating Reference databases for Amplicon-Based Sequencing), a software package to create custom reference databases for metabarcoding studies. crabs includes tools to download sequences from multiple online repositories (i.e., NCBI, BOLD, EMBL, MitoFish), retrieve amplicon regions through in silico PCR analysis and pairwise global alignments, curate the database through multiple filtering parameters (e.g., dereplication, sequence length, sequence quality, unresolved taxonomy, inclusion/exclusion filter), export the reference database in multiple formats for immediate use in taxonomy assignment software, and investigate the reference database through implemented visualizations for diversity, primer efficiency, reference sequence length, database completeness and taxonomic resolution. crabs is a versatile tool for generating curated reference databases of user-specified genetic markers to aid taxonomy assignment from metabarcoding sequencing data. crabs can be installed via docker and is available for download as a conda package and via GitHub ( https://github.com/gjeunen/reference_database_creator ).     

Blind assessment of vertebrate taxonomic diversity across spatial scales by clustering environmental DNA metabarcoding sequences

Human activities impact all ecosystems on Earth, which urges scientists to better understand biodiversity changes across temporal and spatial scales. Environmental DNA (eDNA) metabarcoding is a promising non-invasive method to assess species composition in a wide range of ecosystems. Yet, this method requires the completeness of a reference database, i.e. a list of DNA sequences attached to each species of the regional pool, which is rarely met. As an alternative, molecular operational taxonomic units (MOTUs) can be extracted as clusters of sequences. However, the extent to which the diversity of MOTUs can predict the diversity of species across spatial scales is unknown. Here, we used 196 samples along the Rhone river (France) for which the reference database is complete to assess whether a blind eDNA approach can reliably predict the ground-truth number of species at different spatial scales. Using the 12S rDNA teleo primer, we curated and clustered 60 million sequences into MOTUs using a new assembled bioinformatic pipeline. We show that stringent quality filters were necessary to remove artefact noise, notably MOTUs present in a single PCR replicate, which represented 55% of MOTUs (103). Post-clustering cleaning also removed 19 additional erroneous MOTUs and only discarded one truly present species. We then show that the diversity of retained fish MOTUs accurately predicted the local (α, r = 0.98) and regional (γ) ground-truth species diversity (67 MOTUs versus 63 species), but also the species dissimilarity between samples (β-diversity, r = 0.98). This work paves the way towards extending the use of eDNA metabarcoding in community ecology and biogeography despite major gaps in genetic reference databases.     

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 .     

A DNA barcode library for 5,200 German flies and midges (Insecta: Diptera) and its implications for metabarcoding‐based biomonitoring

This study summarizes results of a DNA barcoding campaign on German Diptera, involving analysis of 45,040 specimens. The resultant DNA barcode library includes records for 2,453 named species comprising a total of 5,200 barcode index numbers (BINs), including 2,700 COI haplotype clusters without species‐level assignment, so called “dark taxa.” Overall, 88 out of 117 families (75%) recorded from Germany were covered, representing more than 50% of the 9,544 known species of German Diptera. Until now, most of these families, especially the most diverse, have been taxonomically inaccessible. By contrast, within a few years this study provided an intermediate taxonomic system for half of the German Dipteran fauna, which will provide a useful foundation for subsequent detailed, integrative taxonomic studies. Using DNA extracts derived from bulk collections made by Malaise traps, we further demonstrate that species delineation using BINs and operational taxonomic units (OTUs) constitutes an effective method for biodiversity studies using DNA metabarcoding. As the reference libraries continue to grow, and gaps in the species catalogue are filled, BIN lists assembled by metabarcoding will provide greater taxonomic resolution. The present study has three main goals: (a) to provide a DNA barcode library for 5,200 BINs of Diptera; (b) to demonstrate, based on the example of bulk extractions from a Malaise trap experiment, that DNA barcode clusters, labelled with globally unique identifiers (such as OTUs and/or BINs), provide a pragmatic, accurate solution to the “taxonomic impediment”; and (c) to demonstrate that interim names based on BINs and OTUs obtained through metabarcoding provide an effective method for studies on species‐rich groups that are usually neglected in biodiversity research projects because of their unresolved taxonomy.     

Taxonomy‐free molecular diatom index for high‐throughput eDNA biomonitoring

Current biodiversity assessment and biomonitoring are largely based on the morphological identification of selected bioindicator taxa. Recently, several attempts have been made to use eDNA metabarcoding as an alternative tool. However, until now, most applied metabarcoding studies have been based on the taxonomic assignment of sequences that provides reference to morphospecies ecology. Usually, only a small portion of metabarcoding data can be used due to a limited reference database and a lack of phylogenetic resolution. Here, we investigate the possibility to overcome these limitations using a taxonomy‐free approach that allows the computing of a molecular index directly from eDNA data without any reference to morphotaxonomy. As a case study, we use the benthic diatoms index, commonly used for monitoring the biological quality of rivers and streams. We analysed 87 epilithic samples from Swiss rivers, the ecological status of which was established based on the microscopic identification of diatom species. We compared the diatom index derived from eDNA data obtained with or without taxonomic assignment. Our taxonomy‐free approach yields promising results by providing a correct assessment for 77% of examined sites. The main advantage of this method is that almost 95% of OTUs could be used for index calculation, compared to 35% in the case of the taxonomic assignment approach. Its main limitations are under‐sampling and the need to calibrate the index based on the microscopic assessment of diatoms communities. However, once calibrated, the taxonomy‐free molecular index can be easily standardized and applied in routine biomonitoring, as a complementary tool allowing fast and cost‐effective assessment of the biological quality of watercourses.     

Species‐level biodiversity assessment using marine environmental DNA metabarcoding requires protocol optimization and standardization

DNA extraction from environmental samples (environmental DNA; eDNA) for metabarcoding‐based biodiversity studies is gaining popularity as a noninvasive, time‐efficient, and cost‐effective monitoring tool. The potential benefits are promising for marine conservation, as the marine biome is frequently under‐surveyed due to its inaccessibility and the consequent high costs involved. With increasing numbers of eDNA‐related publications have come a wide array of capture and extraction methods. Without visual species confirmation, inconsistent use of laboratory protocols hinders comparability between studies because the efficiency of target DNA isolation may vary. We determined an optimal protocol (capture and extraction) for marine eDNA research based on total DNA yield measurements by comparing commonly employed methods of seawater filtering and DNA isolation. We compared metabarcoding results of both targeted (small taxonomic group with species‐level assignment) and universal (broad taxonomic group with genus/family‐level assignment) approaches obtained from replicates treated with the optimal and a low‐performance capture and extraction protocol to determine the impact of protocol choice and DNA yield on biodiversity detection. Filtration through cellulose‐nitrate membranes and extraction with Qiagen's DNeasy Blood & Tissue Kit outperformed other combinations of capture and extraction methods, showing a ninefold improvement in DNA yield over the poorest performing methods. Use of optimized protocols resulted in a significant increase in OTU and species richness for targeted metabarcoding assays. However, changing protocols made little difference to the OTU and taxon richness obtained using universal metabarcoding assays. Our results demonstrate an increased risk of false‐negative species detection for targeted eDNA approaches when protocols with poor DNA isolation efficacy are employed. Appropriate optimization is therefore essential for eDNA monitoring to remain a powerful, efficient, and relatively cheap method for biodiversity assessments. For seawater, we advocate filtration through cellulose‐nitrate membranes and extraction with Qiagen's DNeasy Blood & Tissue Kit or phenol‐chloroform‐isoamyl for successful implementation of eDNA multi‐marker metabarcoding surveys.     

Towards next-generation biodiversity assessment using DNA metabarcoding

Virtually all empirical ecological studies require species identification during data collection. DNA metabarcoding refers to the automated identification of multiple species from a single bulk sample containing entire organisms or from a single environmental sample containing degraded DNA (soil, water, faeces, etc.). It can be implemented for both modern and ancient environmental samples. The availability of next-generation sequencing platforms and the ecologists' need for high-throughput taxon identification have facilitated the emergence of DNA metabarcoding. The potential power of DNA metabarcoding as it is implemented today is limited mainly by its dependency on PCR and by the considerable investment needed to build comprehensive taxonomic reference libraries. Further developments associated with the impressive progress in DNA sequencing will eliminate the currently required DNA amplification step, and comprehensive taxonomic reference libraries composed of whole organellar genomes and repetitive ribosomal nuclear DNA can be built based on the well-curated DNA extract collections maintained by standardized barcoding initiatives. The near-term future of DNA metabarcoding has an enormous potential to boost data acquisition in biodiversity research.     

Unlocking biodiversity and conservation studies in high‐diversity environments using environmental DNA (eDNA): A test with Guianese freshwater fishes

Determining the species compositions of local assemblages is a prerequisite to understanding how anthropogenic disturbances affect biodiversity. However, biodiversity measurements often remain incomplete due to the limited efficiency of sampling methods. This is particularly true in freshwater tropical environments that host rich fish assemblages, for which assessments are uncertain and often rely on destructive methods. Developing an efficient and nondestructive method to assess biodiversity in tropical freshwaters is highly important. In this study, we tested the efficiency of environmental DNA (eDNA) metabarcoding to assess the fish diversity of 39 Guianese sites. We compared the diversity and composition of assemblages obtained using traditional and metabarcoding methods. More than 7,000 individual fish belonging to 203 Guianese fish species were collected by traditional sampling methods, and ~17 million reads were produced by metabarcoding, among which ~8 million reads were assigned to 148 fish taxonomic units, including 132 fish species. The two methods detected a similar number of species at each site, but the species identities partially matched. The assemblage compositions from the different drainage basins were better discriminated using metabarcoding, revealing that while traditional methods provide a more complete but spatially limited inventory of fish assemblages, metabarcoding provides a more partial but spatially extensive inventory. eDNA metabarcoding can therefore be used for rapid and large‐scale biodiversity assessments, while at a local scale, the two approaches are complementary and enable an understanding of realistic fish biodiversity.     

Detection of rare and invasive freshwater fish species using eDNA pyrosequencing: Lake Iznik ichthyofauna revised

Assessment of fish biodiversity in freshwater environments is challenging, especially when rare species or species with low population densities exist. Environmental DNA is becoming a common tool in molecular ecology to detect target species found in the environment. Moreover, eDNA metabarcoding is now used to determine a complete list of target organisms without any prior knowledge on the species inhabiting the environment. This study is the first environmental DNA study designed to assess complete ichthyofauna of the largest lake in Marmara Region of Turkey. For this purpose, an eDNA metabarcoding approach enhanced with tagged primers according to sampling stations for a station specific species listing was used to revise the ichthyofauna of Lake Iznik. Results of pyrosequencing data indicate the presence of 23 species in the lake, five of which are reported for the first time. Short fragment of cytochrome b gene sequences amplified in this study were able to make identifications at species level and the eDNA metabarcoding approach was more cost effective and precise compared to conventional surveys. More molecular data from further studies will enhance the reference databases and eDNA metabarcoding could be used more efficiently as an important molecular tool in biodiversity assessment studies.     

DNA metabarcoding and the cytochrome c oxidase subunit I marker: not a perfect match

DNA metabarcoding enables efficient characterization of species composition in environmental DNA or bulk biodiversity samples, and this approach is making significant and unique contributions in the field of ecology. In metabarcoding of animals, the cytochrome c oxidase subunit I (COI) gene is frequently used as the marker of choice because no other genetic region can be found in taxonomically verified databases with sequences covering so many taxa. However, the accuracy of metabarcoding datasets is dependent on recovery of the targeted taxa using conserved amplification primers. We argue that COI does not contain suitably conserved regions for most amplicon-based metabarcoding applications. Marker selection deserves increased scrutiny and available marker choices should be broadened in order to maximize potential in this exciting field of research.     

DNA宏条形码技术在海洋浮游动物多样性和生态学研究中的应用

浮游动物是海洋生态系统的关键类群,其覆盖门类广泛,多样性高。传统形态鉴定技术需要检测人员具备专业的形态鉴定知识,且费时费力。宏条形码技术无需分离生物个体,而是提取拖网采集到的浮游动物混合样本的总DNA,或者水体中的环境DNA （eDNA）,依托高通量测序平台测序,能够实现对大规模样本快速、准确、经济的分析,在海洋浮游动物生态学研究中得到越来越广泛的应用。分析了DNA宏条形码技术常用的核糖体和线粒体分子标记,在浮游动物多样性和数量研究中的可靠性和不足,并给出在海洋浮游动物群落监测,食物关系分析及生物入侵早期预警等研究中的应用。未来,开发多基因片段组合条形码,发展完备的参考数据库及实现准确的量化研究是DNA宏条形码技术发展的重要方向。     

基于环境DNA-宏条形码技术的底栖动物监测及水质评价研究进展

底栖动物是淡水生态系统中物种多样性最高的类群,也是应用最广泛的水质监测指示生物之一。传统的底栖动物监测以形态学为基础,耗时费力,无法满足流域尺度大规模监测的需求。环境DNA-宏条形码技术是一种新兴的生物监测方法,其与传统方法相比优势在于采样方法简单、低成本、高灵敏度,不受生物样本和环境状况的影响,不依赖分类专家和鉴定资料,能够快速准确地对多个类群进行大规模、高通量的物种鉴定。然而,在实际应用中该方法的效果受诸多因素的影响,不同的方法、流程往往会产生差异较大的结果。鉴于此,着重分析总结了应用环境DNA-宏条形码技术监测底栖动物的关键影响因素,包括样品采集与处理流程、分子标记选择、引物设计、PCR偏好性、参考数据库的完整性及相应的优化。并基于此探讨了提高环境DNA-宏条形码技术在底栖动物监测效率和准确率的途径,以期为底栖动物环境DNA-宏条形码监测方案的制定提供可靠的参考。最后对该技术在底栖动物监测和水质评价中的最新发展方向进行了展望。     

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

The genomic revolution has fundamentally changed how we survey biodiversity on earth. High‐throughput sequencing (“HTS”) platforms now enable the rapid sequencing of DNA from diverse kinds of environmental samples (termed “environmental DNA” or “eDNA”). Coupling HTS with our ability to associate sequences from eDNA with a taxonomic name is called “eDNA metabarcoding” and offers a powerful molecular tool capable of noninvasively surveying species richness from many ecosystems. Here, we review the use of eDNA metabarcoding for surveying animal and plant richness, and the challenges in using eDNA approaches to estimate relative abundance. We highlight eDNA applications in freshwater, marine and terrestrial environments, and in this broad context, we distill what is known about the ability of different eDNA sample types to approximate richness in space and across time. We provide guiding questions for study design and discuss the eDNA metabarcoding workflow with a focus on primers and library preparation methods. We additionally discuss important criteria for consideration of bioinformatic filtering of data sets, with recommendations for increasing transparency. Finally, looking to the future, we discuss emerging applications of eDNA metabarcoding in ecology, conservation, invasion biology, biomonitoring, and how eDNA metabarcoding can empower citizen science and biodiversity education.     

Carrion fly‐derived DNA metabarcoding is an effective tool for mammal surveys: Evidence from a known tropical mammal community

Metabarcoding of vertebrate DNA derived from carrion flies has been proposed as a promising tool for biodiversity monitoring. To evaluate its efficacy, we conducted metabarcoding surveys of carrion flies on Barro Colorado Island (BCI), Panama, which has a well‐known mammal community, and compared our results against diurnal transect counts and camera trapping. We collected 1,084 flies in 29 sampling days, conducted metabarcoding with mammal‐specific (16S) and vertebrate‐specific (12S) primers, and sequenced amplicons on Illumina MiSeq. For taxonomic assignment, we compared blast with the new program protax , and we found that protax improved species identifications. We detected 20 mammal, four bird, and one lizard species from carrion fly metabarcoding, all but one of which are known from BCI. Fly metabarcoding detected more mammal species than concurrent transect counts (29 sampling days, 13 species) and concurrent camera trapping (84 sampling days, 17 species), and detected 67% of the number of mammal species documented by 8 years of transect counts and camera trapping combined, although fly metabarcoding missed several abundant species. This study demonstrates that carrion fly metabarcoding is a powerful tool for mammal biodiversity surveys and has the potential to detect a broader range of species than more commonly used methods.     

obitools: a unix‐inspired software package for DNA metabarcoding

DNA metabarcoding offers new perspectives in biodiversity research. This recently developed approach to ecosystem study relies heavily on the use of next‐generation sequencing (NGS) and thus calls upon the ability to deal with huge sequence data sets. The obitools package satisfies this requirement thanks to a set of programs specifically designed for analysing NGS data in a DNA metabarcoding context. Their capacity to filter and edit sequences while taking into account taxonomic annotation helps to set up tailor‐made analysis pipelines for a broad range of DNA metabarcoding applications, including biodiversity surveys or diet analyses. The obitools package is distributed as an open source software available on the following website: http://metabarcoding.org/obitools . A Galaxy wrapper is available on the GenOuest core facility toolshed: http://toolshed.genouest.org .     

eDNA metabarcoding as a new surveillance approach for coastal Arctic biodiversity

Because significant global changes are currently underway in the Arctic, creating a large‐scale standardized database for Arctic marine biodiversity is particularly pressing. This study evaluates the potential of aquatic environmental DNA (eDNA) metabarcoding to detect Arctic coastal biodiversity changes and characterizes the local spatio‐temporal distribution of eDNA in two locations. We extracted and amplified eDNA using two COI primer pairs from ~80 water samples that were collected across two Canadian Arctic ports, Churchill and Iqaluit, based on optimized sampling and preservation methods for remote regions surveys. Results demonstrate that aquatic eDNA surveys have the potential to document large‐scale Arctic biodiversity change by providing a rapid overview of coastal metazoan biodiversity, detecting nonindigenous species, and allowing sampling in both open water and under the ice cover by local northern‐based communities. We show that DNA sequences of ~50% of known Canadian Arctic species and potential invaders are currently present in public databases. A similar proportion of operational taxonomic units was identified at the species level with eDNA metabarcoding, for a total of 181 species identified at both sites. Despite the cold and well‐mixed coastal environment, species composition was vertically heterogeneous, in part due to river inflow in the estuarine ecosystem, and differed between the water column and tide pools. Thus, COI‐based eDNA metabarcoding may quickly improve large‐scale Arctic biomonitoring using eDNA, but we caution that aquatic eDNA sampling needs to be standardized over space and time to accurately evaluate community structure changes.     

DNA barcoding of economically important freshwater fish species from north‐central Nigeria uncovers cryptic diversity

This study examines the utility of morphology and DNA barcoding in species identification of freshwater fishes from north‐central Nigeria. We compared molecular data (mitochondrial cytochrome c oxidase subunit I (COI) sequences) of 136 de novo samples from 53 morphologically identified species alongside others in GenBank and BOLD databases. Using DNA sequence similarity‐based (≥97% cutoff) identification technique, 50 (94.30%) and 24 (45.30%) species were identified to species level using GenBank and BOLD databases, respectively. Furthermore, we identified cases of taxonomic problems in 26 (49.00%) morphologically identified species. There were also four (7.10%) cases of mismatch in DNA barcoding in which our query sequence in GenBank and BOLD showed a sequence match with different species names. Using DNA barcode reference data, we also identified four unknown fish samples collected from fishermen to species level. Our Neighbor‐joining (NJ) tree analysis recovers several intraspecific species clusters with strong bootstrap support (≥95%). Analysis uncovers two well‐supported lineages within Schilbe intermedius. The Bayesian phylogenetic analyses of Nigerian S. intermedius with others from GenBank recover four lineages. Evidence of genetic structuring is consistent with geographic regions of sub‐Saharan Africa. Thus, cryptic lineage diversity may illustrate species’ adaptive responses to local environmental conditions. Finally, our study underscores the importance of incorporating morphology and DNA barcoding in species identification. Although developing a complete DNA barcode reference library for Nigerian ichthyofauna will facilitate species identification and diversity studies, taxonomic revisions of DNA sequences submitted in databases alongside voucher specimens are necessary for a reliable taxonomic and diversity inventory.     

DNA from soil mirrors plant taxonomic and growth form diversity

Ecosystems across the globe are threatened by climate change and human activities. New rapid survey approaches for monitoring biodiversity would greatly advance assessment and understanding of these threats. Taking advantage of next-generation DNA sequencing, we tested an approach we call metabarcoding: high-throughput and simultaneous taxa identification based on a very short (usually <100 base pairs) but informative DNA fragment. Short DNA fragments allow the use of degraded DNA from environmental samples. All analyses included amplification using plant-specific versatile primers, sequencing and estimation of taxonomic diversity. We tested in three steps whether degraded DNA from dead material in soil has the potential of efficiently assessing biodiversity in different biomes. First, soil DNA from eight boreal plant communities located in two different vegetation types (meadow and heath) was amplified. Plant diversity detected from boreal soil was highly consistent with plant taxonomic and growth form diversity estimated from conventional above-ground surveys. Second, we assessed DNA persistence using samples from formerly cultivated soils in temperate environments. We found that the number of crop DNA sequences retrieved strongly varied with years since last cultivation, and crop sequences were absent from nearby, uncultivated plots. Third, we assessed the universal applicability of DNA metabarcoding using soil samples from tropical environments: a large proportion of species and families from the study site were efficiently recovered. The results open unprecedented opportunities for large-scale DNA-based biodiversity studies across a range of taxonomic groups using standardized metabarcoding approaches.