Over 2.5 million COI sequences in GenBank and growing

Open Data In Neurophysiology: Advancements, Solutions & Challenges

Across the life sciences, an ongoing effort over the last 50 years has made data and methods more reproducible and transparent. This openness has led to transformative insights and vastly accelerated scientific progress1,2. For example, structural biology3 and genomics4,5 have undertaken systematic collection and publication of protein sequences and structures over the past half-century, and these data have led to scientific breakthroughs that were unthinkable when data collection first began (e.g.6). We believe that neuroscience is poised to follow the same path, and that principles of open data and open science will transform our understanding of the nervous system in ways that are impossible to predict at the moment. To this end, new social structures along with active and open scientific communities are essential7 to facilitate and expand the still limited adoption of open science practices in our field8. Unified by shared values of openness, we set out to organize a symposium for Open Data in Neuroscience (ODIN) to strengthen our community and facilitate transformative neuroscience research at large. In this report, we share what we learned during this first ODIN event. We also lay out plans for how to grow this movement, document emerging conversations, and propose a path toward a better and more transparent science of tomorrow.

The development of multiplex PCR assays for the rapid identification of multiple Saccostrea species, and their practical applications in restoration and aquaculture

BACKGROUND

The ecology and biology of oysters (Ostreidae) across the tropics is poorly understood. Morphological plasticity and shared characteristics among oysters have resulted in the misidentification of species, creating challenges for understanding basic species-specific biological information that is required for restoration and aquaculture. Genetic barcoding has proven essential for accurate species identification and understanding species geographic ranges. To reduce the costs of molecular species identification we developed multiplex assays using the cytochrome c oxidase subunit I (COI or cox1) barcoding gene for the rapid identification of five species of oysters within the genus Saccostrea that are commonly found in Queensland, Australia: Saccostrea glomerata, Saccostrea lineage B, Saccostrea lineage F, Saccostrea lineage G, and Saccostrea spathulata (lineage J).

RESULTS

Multiplex assays were successful in species-specific amplification of targeted species. The practical application of these primers was tested on wild spat collected from a pilot restoration project in Moreton Bay, Queensland, with identified species (S. glomerata, lineage B and lineage G) validated by Sanger sequencing. DNA sampling by extraction of oyster pallial fluid was also tested on adult oysters collected from the Noosa estuary in Queensland to assess whether oysters were able to be identified non-destructively. DNA concentrations as low as 1 ng/ μL still amplified in most cases, allowing for identification, and mortality at 6 weeks post pallial fluid collection was low (3 out of 104 sampled oysters).

CONCLUSION

These multiplex assays will be essential tools for species identification in future studies, and we successfully demonstrate their practical application in both restoration and aquaculture contexts in Queensland. The multiplex assays developed in this study outline easily replicable methods for the development of additional species-specific primer sets for the rapid identification of other species of Saccostrea found across the Indo-Pacific, which will be instrumental in unravelling the taxonomic ambiguities within this genus in tropical regions.

Data science and its future in large neuroscience collaborations

The rise of large scientific collaborations in neuroscience requires systematic, scalable, and reliable data management. How this is best done in practice remains an open question. To address this, we conducted a data science survey among currently active U19 grants, funded through the NIH's BRAIN Initiative. The survey was answered by both data science liaisons and Principal Investigators, speaking for ~500 researchers across 21 nation-wide collaborations. We describe the tools, technologies, and methods currently in use, and identify several shortcomings of current data science practice. Building on this survey, we develop plans and propose policies to improve data collection, use, publication, reuse and training in the neuroscience community.

Metabarcoding of zooplankton communities of Dianchi Lake based on the mitochondrial cytochrome oxidase subunit 1 gene

Freshwater lakes as an essential component of the ecosystem, provide ecological resources in addition to economic source for humans. Under recent climate change scenario, preserving the biodiversity of freshwater ecosystems is crucial. This study aimed to characterize the diversity of zooplankton communities in Dianchi Lake, located in Kunming Municipality, Yunnan Province, China, using Illumina high-throughput sequencing of the cytochrome oxidase subunit 1 (COI) gene marker. A total of 18 water samples were collected including 16 from the outer sea area of Dianchi Lake: 4 from the east (E1-4), 4 from the west (W1-4), 4 from the south (S1-4), and 4 from the north (N1-4), and: 2 from the Caohai area (C1-2) as research sites. All environmental parameters including pH, ammonium (NH4+), total nitrogen (TN), total phosphorus (TP), chlorophyll a content (CHLA) were found to be insignificant (p > 0.05), except for chemical oxygen demand (COD) and transparency (T), which were found to be significant (p < 0.05). Alpha diversity indices including ACE, Chao1, Shannon, and Simpson showed non-significant differences (p > 0.05), indicating no variation in the richness of zooplankton communities at different locations of Dianchi Lake. However, principal coordinate analysis (PCoA) showed that most of the samples from East, West, and South groups were close to each other, showing more similarities among them, while Caohai and North group samples were distant from each other, showing more differences with other groups. Rotifera, Arthropoda, and Chordata were the top three phyla, while Keratella, Macrothrix, and Brachionus were the dominant genera. Mantel test analysis showed that COD and transparency were important environmental factors that shaped the Rotifera community structure of Dianchi Lake. In conclusion, this study provides insights on conserving the diversity of zooplankton communities in Dianchi Lake, especially by controlling COD and maintaining water transparency, in order to preserve its ecological resources and economic significance.

Towards a Canary Islands barcode database for soil biodiversity: revealing cryptic and unrecorded mite species diversity within insular soils

Soil arthropod diversity contributes to a high proportion of the total biodiversity on Earth. However, most soil arthropods are still undescribed, hindering our understanding of soil functioning and global biodiversity estimations. Inventorying soil arthropods using conventional taxonomical approaches is particularly difficult and costly due to the great species richness, abundance and local-scale heterogeneity of mesofauna communities and the poor taxonomic background knowledge of most lineages. To alleviate this situation, we have designed and implemented a molecular barcoding framework adapted to soil fauna. This pipeline includes different steps, starting with a morphology-based selection of specimens which are imaged. Then, DNA is extracted non-destructively. Both images and voucher specimens are used to assign a taxonomic identification, based on morphology that is further checked for consistency with molecular information. Using this procedure, we studied 239 specimens of mites from the Canary Islands including representatives of Mesostigmata, Sarcoptiformes and Trombidiformes, of which we recovered barcode sequences for 168 specimens that were morphologically identified to 49 species, with nine specimens that could only be identified at the genus or family levels. Multiple species delimitation analyses were run to compare molecular delimitations with morphological identifications, including ASAP, mlPTP, BINs and 3% and 8% genetic distance thresholds. Additionally, a species-level search was carried out at the Biodiversity Databank of the Canary Islands (BIOTA) to evaluate the number of species in our dataset that were not previously recorded in the archipelago. In parallel, a sequence-level search of our sequences was performed against BOLD Systems. Our results reveal that multiple morphologically identified species correspond to different molecular lineages, which points to significant levels of unknown cryptic diversity within the archipelago. In addition, we evidenced that multiple species in our dataset constituted new records for the Canary Islands fauna and that the information for these lineages within online genetic repositories is very incomplete. Our study represents the first systematic effort to catalogue the soil arthropod mesofauna of the Canary Islands and establishes the basis for the Canary Islands Soil Biodiversity barcode database. This resource will constitute a step forward in the knowledge of these arthropods in a region of special interest.

Creating Virtuous Cycles for DNA Barcoding: A Case Study in Science Innovation, Entrepreneurship, and Diplomacy

This chapter on the history of the DNA barcoding enterprise attempts to set the stage for the more scholarly contributions in this volume by addressing the following questions. How did the DNA barcoding enterprise begin? What were its goals, how did it develop, and to what degree are its goals being realized? We have taken a keen interest in the barcoding movement and its relationship to taxonomy, collections, and biodiversity informatics more broadly considered. This chapter integrates our two different perspectives on barcoding. DES was the Executive Secretary of the Consortium for the Barcode of Life from 2004 to 2017, with the mission to support the success of DNA barcoding without being directly involved in generating barcode data. RDMP viewed barcoding as an important entry into the landscape of biodiversity data, with many potential linkages to other components of that landscape. We also saw it as a critical step toward the era of international genomic research that was sure to follow. Like the Mercury Program that paved the way for lunar landings by the Apollo Program, we saw DNA barcoding as the proving grounds for the interdisciplinary and international cooperation that would be needed for success of whole-genome research.

taxalogue: a toolkit to create comprehensive CO1 reference databases

Background

Taxonomic identification through DNA barcodes gained considerable traction through the invention of next-generation sequencing and DNA metabarcoding. Metabarcoding allows for the simultaneous identification of thousands of organisms from bulk samples with high taxonomic resolution. However, reliable identifications can only be achieved with comprehensive and curated reference databases. Therefore, custom reference databases are often created to meet the needs of specific research questions. Due to taxonomic inconsistencies, formatting issues, and technical difficulties, building a custom reference database requires tremendous effort. Here, we present taxalogue, an easy-to-use software for creating comprehensive and customized reference databases that provide clean and taxonomically harmonized records. In combination with extensive geographical filtering options, taxalogue opens up new possibilities for generating and testing evolutionary hypotheses.

Methods

taxalogue collects DNA sequences from several online sources and combines them into a reference database. Taxonomic incongruencies between the different data sources can be harmonized according to available taxonomies. Dereplication and various filtering options are available regarding sequence quality or metadata information. taxalogue is implemented in the open-source Ruby programming language, and the source code is available at https://github.com/nwnoll/taxalogue. We benchmark four reference databases by sequence identity against eight queries from different localities and trapping devices. Subsamples from each reference database were used to compare how well another one is covered.

Results

taxalogue produces reference databases with the best coverage at high identities for most tested queries, enabling more accurate, reliable predictions with higher certainty than the other benchmarked reference databases. Additionally, the performance of taxalogue is more consistent while providing good coverage for a variety of habitats, regions, and sampling methods. taxalogue simplifies the creation of reference databases and makes the process reproducible and transparent. Multiple available output formats for commonly used downstream applications facilitate the easy adoption of taxalogue in many different software pipelines. The resulting reference databases improve the taxonomic classification accuracy through high coverage of the query sequences at high identities.

Taxonomy and Biological Control: New Challenges in an Old Relationship

Biological control and taxonomy are continuously developing fields with remarkable impacts on society. At least 80 years of literature have documented this relationship, which remains essentially the same in its mutualistic nature, as well as in its major challenges. From the perspective of Brazilian taxonomists, we discuss the impacts of important scientific and social developments that directly affect research in these areas, posing new challenges for this lasting relationship. The increasing restrictions and concerns regarding the international transit of organisms require improvements in research related to risk assessment for exotic biological control agents and also stimulate prospecting within the native biota. In our view, this is a positive situation that can foster a closer relationship between taxonomists and applied entomologists, as well as local surveys and taxonomic studies that are necessary before new programs and agents can be implemented. We discuss the essential role of molecular biology in this context, as an iconic example of the synergy between applied sciences and natural history. As our society comes to need safer and more sustainable solutions for food security and the biodiversity crisis, scientific progress will build upon this integration, where biological control and taxonomy play an essential role.

Diet of a threatened rattlesnake (eastern massasauga) revealed by DNA metabarcoding

Characterizing the diet of imperiled species using minimally invasive methods is crucial to understanding their ecology and conservation requirements. Here, we apply a DNA metabarcoding approach to study the diet of the eastern massasauga rattlesnake (Sistrurus catenatus), a Federally Threatened snake found throughout the Great Lakes region. Eighty-three fecal samples collected across 10 different massasauga populations located in Michigan, USA, were sequenced, with 70 samples containing prey DNA. We used universal metazoan primers and developed a host-specific oligonucleotide blocker to characterize their diet. We identified at least 12 different prey species, with eastern massasaugas exhibiting opportunistic feeding and a strong preference towards small mammals. Meadow voles (Microtus pennsylvanicus) were the most common prey item (70% of diet) followed by the northern short-tailed shrew (Blarina brevicauda) and masked shrew (Sorex cinereus; 15.7% of diet each), along with occasional bird and snake prey. Adult individuals exhibited a more generalized diet, consuming a larger number of prey taxa on average. Younger snakes consumed a smaller variety of prey items and tended to consume smaller-sized mammals such as masked shrews (Sorex cinereus) and northern short-tailed shrews (Blarina brevicauda). We conclude that small mammals are a crucial part of eastern massasauga rattlesnake diet and recommend this be taken into consideration when conservation strategies are developed. The methods developed in this study can be applied to other reptile species, providing an accurate, minimally invasive, and thorough diet assessment for at-risk reptile species.

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).

Maximizing the reliability and the number of species assignments in metabarcoding studies using a curated regional library and a public repository

Biodiversity assessments relying on DNA have increased rapidly over the last decade. However, the reliability of taxonomic assignments in metabarcoding studies is variable and affected by the reference databases and the assignment methods used. Species level assignments are usually considered as reliable using regional libraries but unreliable using public repositories. In this study, we aimed to test this assumption for metazoan species detected in the Gulf of St. Lawrence in the Northwest Atlantic. We first created a regional library (GSL-rl) by data mining COI barcode sequences from BOLD, and included a reliability ranking system for species assignments. We then estimated 1) the accuracy and precision of the public repository NCBI-nt for species assignments using sequences from the regional library and 2) compared the detection and reliability of species assignments of a metabarcoding dataset using either NCBI-nt or the regional library and popular assignment methods. With NCBI-nt and sequences from the regional library, the BLAST-LCA (least common ancestor) method was the most precise method for species assignments, but the accuracy was higher with the BLAST-TopHit method (&gt;80% over all taxa, between 70% and 90% amongst taxonomic groups). With the metabarcoding dataset, the reliability of species assignments was greater using GSL-rl compared to NCBI-nt. However, we also observed that the total number of reliable species assignments could be maximized using both GSL-rl and NCBI-nt with different optimized assignment methods. The use of a two-step approach for species assignments, i.e., using a regional library and a public repository, could improve the reliability and the number of detected species in metabarcoding studies.

Demographic changes and life-history strategies predict the genetic diversity in crabs

Uncovering what predicts genetic diversity (GD) within species can help us access the status of populations and their evolutionary potential. Traits related to effective population size show a proportional association to GD, but evidence supports life-history strategies and habitat as the drivers of GD variation. Instead of investigating highly divergent taxa, focusing on one group could help to elucidate the factors influencing the GD. Additionally, most empirical data is based on vertebrate taxa; therefore, we might be missing novel patterns of GD found in neglected invertebrate groups. Here, we investigated the predictors of the GD in crabs (Brachyura) by compiling the most comprehensive cytochrome c oxidase subunit I (COI) available. Eight predictor variables were analysed across 150 species (16 992 sequences) using linear models (multiple linear regression) and comparative methods (PGLS). Our results indicate that population size fluctuation represents the most critical trait predicting GD, with species that have undergone bottlenecks followed by population expansion showing lower GD. Egg size, pelagic larval duration and habitat might play a role probably because of their association with how species respond to disturbances. Ultimately, K-strategists that have undergone bottlenecks are the species showing lower GD. Some variables do not show an association with GD as expected, most likely due to the taxon-specific role of some predictors, which should be considered in further investigations and generalizations. This work highlights the complexity underlying the predictors of GD and adds results from a marine invertebrate group to the current understanding of this topic.

Phylogenetic analyses reveal bat communities in Northwestern Mexico harbor a high diversity of novel cryptic ectoparasite species

Parasites are integral parts of ecosystem function and important drivers of evolutionary processes. Characterizing ectoparasite diversity is fundamental to studies of host-parasite interactions, evolution, and conservation, and also for understanding emerging disease threats for some vector borne pathogens. With more than 1400 species, bats represent the second most speciose mammalian clade, but their ectoparasite fauna are poorly known for most species. We sequenced mitochondrial Cytochrome Oxidase C subunit I and nuclear 18S ribosomal gene fragments, and used Bayesian phylogenetic analyses to characterize ectoparasite taxon identity and diversity for 17 species of parasitized bats sampled along the Baja California peninsula and in Northwestern Mexico. The sequence data revealed multiple novel lineages of bat bugs (Cimicidae), flies (Nycteribiidae and Streblidae), and ticks (Argasidae). Within families, the new linages showed more than 10% sequence divergence, which is consistent with separation at least at the species level. Both families of bat flies showed host specificity, particularly on Myotis species. We also identified new records for the Baja peninsula of one tick (Carios kelleyi), and of five Streblid bat fly species. One Nycteribiid bat fly haplotype from Pallid bat (Antrozous pallidus) hosts was found throughout the peninsula, suggesting potential long distance co-dispersal with hosts. Different bat bug and tick communities were found in the north and south of the peninsula. This study is the first systematic survey of bat ectoparasites in the Baja California peninsula, revealing novel lineages that are highly genetically differentiated from other parts of North America. For some ectoparasite species, haplotype distributions may reflect patterns of bat migration. This work is a first step in characterizing ectoparasite diversity over the Baja California peninsula, and understanding how ecological and evolutionary interactions shape bat ectoparasite communities among host species in different parts of their ranges.

The first DNA barcode library of Chironomidae from the Tibetan Plateau with an evaluation of the status of the public databases

The main aim of this study was to curate a COI barcode library of Chironomidae from the Tibetan Plateau (TP) as an essential supplement to the public database. Another aim is to evaluate the current status of the public database of Chironomidae in aspects of taxonomic coverage, geographic representation, barcode quality, and efficiency for molecular identification, the Tibetan Plateau, China. In this study, 512 individuals of Chironomidae from the TP were identified based on morphological taxonomy and barcode analysis. The metadata of public records of Chironomidae were downloaded from the BOLD, and the quality of the public barcodes was ranked using the BAGS program. The reliability of the public library for molecular identification was evaluated with the newly curated library using the BLAST method. The newly curated library comprised 159 barcode species of 54 genera, of which 58.4% of species were likely new to science. There were great gaps in the taxonomic coverage and geographic representation in the public database, and only 29.18% of barcodes were identified at the species level. The quality of the public database was of concern, with only 20% of species being determined as concordant between BINs and morphological species. The accuracy of molecular identification using the public database was poor, and about 50% of matched barcodes could be correctly identified at the species level at the identity threshold of 97%. Based on these data, some recommendations are included here for improving barcoding studies on Chironomidae. The species richness of Chironomidae from the TP is much higher than ever recorded. Barcodes from more taxonomic groups and geographic regions are urgently needed to fill the great gap in the current public database of Chironomidae. Users should take caution when public databases are adopted as reference libraries for the taxonomic assignment.

Spatial assortment of soil organisms supports the size-plasticity hypothesis

The size-plasticity hypothesis posits that larger size organisms are less plastic in their metabolic rates and, therefore, are more strongly environmental-filtered than smaller organisms. Many studies have supported this hypothesis by evaluating the relative roles of environmental filtration and dispersal for different taxonomic groups of soil organisms. Most observations are made at large spatial scales, which are assumed to have a wide array of varying habitats. However, since urbanization causes habitat fragmentation at smaller regional scales, testing the size-plasticity hypothesis at this scale would help better understand the spatial assortment of urban soil organisms which, in turn, would help to develop improved management and conservation strategies for urban soil health. Here, we used DNA metabarcoding on five groups of soil biota (bacteria, fungi, protists, nematodes, and invertebrates) to assess the relative importance of dispersal and environmental filters to examine the size-plasticity hypothesis at this spatial scale in an urban environment. We observed strong distance-decay of community similarities associated with higher levels of stochastic changes in bacteria, nematode, and protist communities but not fungal or invertebrate communities. Bacterial communities occupied the widest niche followed by protists and nematodes, potentially because of their higher dispersal abilities compared to the larger soil organisms. Null deviation of communities varied with taxonomic groups where bacteria and nematodes were mainly driven by homogenizing dispersal, protists and fungi by drift, and soil invertebrates by environmental selection. We further identified a small percentage of locally-adapted taxa (2.1%) that could be focal taxa for conservation and restoration efforts by, for example, restoring their habitats and enhancing their regional connectivity. These results support the size-plasticity hypothesis at the relatively unexplored regional scale in an urbanization context, and provide new information for improving urban soil health and sustainable city models.

Dietary preferences and feeding strategies of Colombian highland woolly monkeys

Primates are very selective in the foods they include in their diets with foraging strategies that respond to spatial and temporal changes in resource availability, distribution and quality. Colombian woolly monkeys (Lagothrix lagotricha lugens), one of the largest primate species in the Americas, feed mainly on fruits, but they also eat a high percentage of arthropods. This differs from closely related Atelid species that supplement their diet with leaves. In an 11 month study, we investigated the foraging strategies of this endemic monkey and assessed how resource availability affects dietary selection. Using behavioural, phenological, arthropod sampling and metabarcoding methods, we recorded respectively foraging time, forest productivity, arthropod availability in the forest and arthropod consumption. Scat samples and capturing canopy substrates (i.e. moss, bromeliads, aerial insects) were used for assigning arthropod taxonomy. The most important resource in the diet was fruits (54%), followed by arthropods (28%). Resource availability predicted feeding time for arthropods but not for fruits. Further, there was a positive relationship between feeding time on fruits and arthropods, suggesting that eating both resources during the same periods might work as an optimal strategy to maximize nutrient intake. Woolly monkeys preferred and avoided some fruit and arthropod items available in their home range, choosing a wide variety of arthropods. Geometrid moths (Lepidoptera) were the most important and consistent insects eaten over time. We found no differences in the type of arthropods adults and juveniles ate, but adults invested more time foraging for this resource, especially in moss. Although woolly monkeys are generalist foragers, they do not select their food items randomly or opportunistically.

Is species richness mediated by functional and genetic divergence? A global analysis in birds

Unravelling why species richness shows such dramatic spatial variation is an ongoing challenge. Common to many theories is that increasing species richness (e.g. with latitude) requires a compensatory trade-off on an axis of species' ecology. Spatial variation in species richness may also affect genetic diversity if large numbers of coexisting, related species result in smaller population sizes.Here, we test whether increasing species richness results in differential occupation of morphospace by the constituent species, or decreases species' genetic diversity. We test for two potential mechanisms of morphological accommodation: denser packing in ecomorphological space, and expansion of the space. We then test whether species differ in their nucleotide diversity depending on allopatry or sympatry with relatives, indicative of potential genetic consequences of coexistence that would reduce genetic diversity in sympatry. We ask these questions in a spatially explicit framework, using a global database of avian functional trait measurements in combination with >120,000 sequences downloaded from GenBank.We find that higher species richness within families is not systematically correlated with either packing in morphological space or overdispersion but, at the Class level, we find a general positive relationship between packing and species richness, but that points sampled in the tropics have comparatively greater packing than temperate ones relative to their species richness. We find limited evidence that geographical co-occurrence with closely related species or tropical distributions decreases nucleotide diversity of nuclear genes; however, this requires further analysis.Our results suggest that avian families can accumulate species regionally with minimal tradeoffs or cost, implying that external biotic factors do not limit species richness. Read the free Plain Language Summary for this article on the Journal blog.

Exploring Large-Scale Patterns of Genetic Variation in the COI Gene among Insecta: Implications for DNA Barcoding and Threshold-Based Species Delimitation Studies

The genetic variation in the COI gene has had a great effect on the final results of species delimitation studies. However, little research has comprehensively investigated the genetic divergence in COI among Insecta. The fast-growing COI data in BOLD provide an opportunity for the comprehensive appraisal of the genetic variation in COI among Insecta. We calculated the K2P distance of 64,414 insect species downloaded from BOLD. The match ratios of the clustering analysis, based on different thresholds, were also compared among 4288 genera (35,068 species). The results indicate that approximately one-quarter of the species of Insecta showed high intraspecific genetic variation (>3%), and a conservative estimate of this proportion ranges from 12.05% to 22.58%. The application of empirical thresholds (e.g., 2% and 3%) in the clustering analysis may result in the overestimation of the species diversity. If the minimum interspecific genetic distance of the congeneric species is greater than or equal to 2%, it is possible to avoid overestimating the species diversity on the basis of the empirical thresholds. In comparison to the fixed thresholds, the “threshOpt” and “localMinima” algorithms are recommended for the provision of a reference threshold for threshold-based species delimitation studies.

Water eDNA metabarcoding is effective in detecting non-native species in marinas, but detection errors still hinder its use for passive monitoring

Marinas are high-priority targets for marine non-indigenous species (NIS), where they compose a large portion of the biofouling communities. The practicality of water samples collection makes environmental DNA (eDNA) metabarcoding an interesting tool for routine NIS surveys. Here the effectiveness of water-eDNA-metabarcoding to identify biofouling NIS, in 10 marinas from western France, was examined. Morphological identification of specimens collected in quadrats brought out 18 sessile benthic NIS beneath floating pontoons. Water-eDNA-metabarcoding detected two thirds of them, failing to detect important NIS. However, sampling and bioinformatics filtering steps can be optimized to identify more species. In addition, this method allowed the detection of additional NIS from neighboring micro-habitats. Caution should, however, be taken when reporting putative novel NIS, because of errors in species assignment. This work highlights that water-eDNA-metabarcoding is effective for active (targeted) NIS surveys and could be significantly improved for its further use in marine NIS passive surveys.

An evaluation of errors in the mitochondrial COI sequences of Hydrachnidia (Acari, Parasitengona) in public databases

Public molecular databases are fundamental tools for modern taxonomic studies whose usefulness rely on the soundness of the data within them. Here, we study potential errors that can arise along the data pipeline from sampling, specimen identification and molecular processing (digestion, amplification and sequencing) to the submission of sequences to these databases by using the DNA sequences of Hydrachnidia (Acari, Parasitengona) as a case study. Our results indicate that molecular information is available for only about 3% of the Hydrachnidia species known to date; yet, within this small percentage, errors are present in almost 5% of the species analyzed (0.5% of the sequences and almost 11% of the genera). This study underscores the scarcity of genetic data available for Hydrachnidia, but also that the proportion of errors in DNA sequences is relatively small. Even so, it highlights the danger associated with using DNA sequences from public databases, particularly for species identification, and reinforces the need for greater quality control measures and/or protocols to avoid an intensification of errors in the (post) genomics era. Finally, our study emphasizes that potential errors may also reveal cryptic diversity within a species.

Revision and annotation of DNA barcode records for marine invertebrates: report of the 8th iBOL conference hackathon

The accuracy of specimen identification through DNA barcoding and metabarcoding relies on reference libraries containing records with reliable taxonomy and sequence quality. The considerable growth in barcode data requires stringent data curation, especially in taxonomically difficult groups such as marine invertebrates. A major effort in curating marine barcode data in the Barcode of Life Data Systems (BOLD) was undertaken during the 8th International Barcode of Life Conference (Trondheim, Norway, 2019). Major taxonomic groups (crustaceans, echinoderms, molluscs, and polychaetes) were reviewed to identify those which had disagreement between Linnaean names and Barcode Index Numbers (BINs). The records with disagreement were annotated with four tags: a) MIS-ID (misidentified, mislabeled, or contaminated records), b) AMBIG (ambiguous records unresolved with the existing data), c) COMPLEX (species names occurring in multiple BINs), and d) SHARE (barcodes shared between species). A total of 83,712 specimen records corresponding to 7,576 species were reviewed and 39% of the species were tagged (7% MIS-ID, 17% AMBIG, 14% COMPLEX, and 1% SHARE). High percentages (&gt;50%) of AMBIG tags were recorded in gastropods, whereas COMPLEX tags dominated in crustaceans and polychaetes. The high proportion of tagged species reflects either flaws in the barcoding workflow (e.g., misidentification, cross-contamination) or taxonomic difficulties (e.g., synonyms, undescribed species). Although data curation is essential for barcode applications, such manual attempts to examine large datasets are unsustainable and automated solutions are extremely desirable.

Apparent density, trypanosome infection rates and host preference of tsetse flies in the sleeping sickness endemic focus of northwestern Uganda

Background

African trypanosomiasis, caused by protozoa of the genus Trypanosoma and transmitted by the tsetse fly, is a serious parasitic disease of humans and animals. Reliable data on the vector distribution, feeding preference and the trypanosome species they carry is pertinent to planning sustainable control strategies.

Methodology

We deployed 109 biconical traps in 10 villages in two districts of northwestern Uganda to obtain information on the apparent density, trypanosome infection status and blood meal sources of tsetse flies. A subset (272) of the collected samples was analyzed for detection of trypanosomes species and sub-species using a nested PCR protocol based on primers amplifying the Internal Transcribed Spacer (ITS) region of ribosomal DNA. 34 blood-engorged adult tsetse midguts were analyzed for blood meal sources by sequencing of the mitochondrial cytochrome c oxidase 1 (COI) and cytochrome b (cytb) genes.

Results

We captured a total of 622 Glossina fuscipes fuscipes tsetse flies (269 males and 353 females) in the two districts with apparent density (AD) ranging from 0.6 to 3.7 flies/trap/day (FTD). 10.7% (29/272) of the flies were infected with one or more trypanosome species. Infection rate was not significantly associated with district of origin (Generalized linear model (GLM), χ² = 0.018, P = 0.895, df = 1, n = 272) and sex of the fly (χ2 = 1.723, P = 0.189, df = 1, n = 272). However, trypanosome infection was highly significantly associated with the fly’s age based on wing fray category (χ² = 22.374, P < 0.001, df = 1, n = 272), being higher among the very old than the young tsetse. Nested PCR revealed several species of trypanosomes: T. vivax (6.62%), T. congolense (2.57%), T. brucei and T. simiae each at 0.73%. Blood meal analyses revealed five principal vertebrate hosts, namely, cattle (Bos taurus), humans (Homo sapiens), Nile monitor lizard (Varanus niloticus), African mud turtle (Pelusios chapini) and the African Savanna elephant (Loxodonta africana).

Conclusion

We found an infection rate of 10.8% in the tsetse sampled, with all infections attributed to trypanosome species that are causative agents for AAT. However, more verification of this finding using large-scale passive and active screening of human and tsetse samples should be done. Cattle and humans appear to be the most important tsetse hosts in the region and should be considered in the design of control interventions.

Profiling of RNA Viruses in Biting Midges (Ceratopogonidae) and Related Diptera from Kenya Using Metagenomics and Metabarcoding Analysis

Vector-borne diseases (VBDs) cause enormous health burden worldwide, as they account for more than 17% of all infectious diseases and over 700,000 deaths each year. A significant number of these VBDs are caused by RNA virus pathogens. Here, we used metagenomics and metabarcoding analysis to characterize RNA viruses and their insect hosts among biting midges from Kenya. We identified a total of 15 phylogenetically distinct insect-specific viruses. These viruses fall into six families, with one virus falling in the recently proposed negevirus taxon. The six virus families include Partitiviridae, Iflaviridae, Tombusviridae, Solemoviridae, Totiviridae, and Chuviridae. In addition, we identified many insect species that were possibly associated with the identified viruses. Ceratopogonidae was the most common family of midges identified. Others included Chironomidae and Cecidomyiidae. Our findings reveal a diverse RNA virome among Kenyan midges that includes previously unknown viruses. Further, metabarcoding analysis based on COI (cytochrome c oxidase subunit 1 mitochondrial gene) barcodes reveal a diverse array of midge species among the insects used in the study. Successful application of metagenomics and metabarcoding methods to characterize RNA viruses and their insect hosts in this study highlights a possible simultaneous application of these two methods as cost-effective approaches to virus surveillance and host characterization. IMPORTANCE The majority of the viruses that currently cause diseases in humans and animals are RNA viruses, and more specifically arthropod-transmitted viruses. They cause diseases such as dengue, West Nile infection, bluetongue disease, Schmallenberg disease, and yellow fever, among others. Several sequencing investigations have shown us that a diverse array of RNA viruses among insect vectors remain unknown. Some of these could be ancient lineages that could aid in comprehensive studies on RNA virus evolution. Such studies may provide us with insights into the evolution of the currently pathogenic viruses. Here, we applied metagenomics to field-collected midges and we managed to characterize several RNA viruses, where we recovered complete and nearly complete genomes of these viruses. We also characterized the insect host species that are associated with these viruses. These results add to the currently known diversity of RNA viruses among biting midges as well as their associated insect hosts.

Excluding spatial sampling bias does not eliminate oversplitting in DNA-based species delimitation analyses

DNA barcoding and DNA-based species delimitation are major tools in DNA taxonomy. Sampling has been a central debate in this context, because the geographical composition of samples affects the accuracy and performance of DNA barcoding. Performance of complex DNA-based species delimitation is to be tested under simpler conditions in absence of geographic sampling bias. Here, we present an empirical dataset sampled from a single locality in a Southeast-Asian biodiversity hotspot (Laos: Phou Pan mountain). We investigate the performance of various species delimitation approaches on a megadiverse assemblage of herbivorous chafer beetles (Coleoptera: Scarabaeidae) to infer whether species delimitation suffers in the same way from exaggerate infraspecific variation despite the lack of geographic genetic variation that led to inconsistencies between entities from DNA-based and morphology-based species inference in previous studies. For this purpose, a 658 bp fragment of the mitochondrial cytochrome c oxidase subunit 1 (cox1) was analyzed for a total of 186 individuals of 56 morphospecies. Tree-based and distance-based species delimitation methods were used. All approaches showed a rather limited match ratio (max. 77%) with morphospecies. Poisson tree process (PTP) and statistical parsimony network analysis (TCS) prevailingly over-splitted morphospecies, while 3% clustering and Automatic Barcode Gap Discovery (ABGD) also lumped several species into one entity. ABGD revealed the highest congruence between molecular operational taxonomic units (MOTUs) and morphospecies. Disagreements between morphospecies and MOTUs have to be explained by historically acquired geographic genetic differentiation, incomplete lineage sorting, and hybridization. The study once again highlights how important morphology still is in order to correctly interpret the results of molecular species delimitation.

Just keep it simple? Benchmarking the accuracy of taxonomy assignment software in metabarcoding studies

How do you put a name on an unknown piece of DNA? From microbes to mammals, high-throughput metabarcoding studies provide a more objective view of natural communities, overcoming many of the inherent limitations of traditional field surveys and microscopy-based observations (Deiner et al., 2017). Taxonomy assignment is one of the most critical aspects of any metabarcoding study, yet this important bioinformatics task is routinely overlooked. Biodiversity surveys and conservation efforts often depend on formal species inventories: the presence (or absence) of species, and the number of individuals reported across space and time. However, computational workflows applied in eukaryotic metabarcoding studies were originally developed for use with bacterial/archaeal data sets, where microbial researchers rely on one conserved locus (nuclear 16S rRNA) and have access to vast databases with good coverage across most prokaryotic lineages - a situation not mirrored in most multicellular taxa. In this issue of Molecular Ecology Resources, Hleap et al. (2021) carry out an extensive benchmarking exercise focused on taxonomy assignment strategies for eukaryotic metabarcoding studies utilizing the mitochondrial Cytochrome C oxidase I marker gene (COI). They assess the performance and accuracy of software tools representing diverse methodological approaches: from "simple" strategies based on sequence similarity and composition, to model-based phylogenetic and probabilistic classification tools. Contrary to popular assumptions, less complex approaches (BLAST and the QIIME2 feature classifier) consistently outperformed more sophisticated mathematical algorithms and were highly accurate for assigning taxonomy at higher levels (e.g. family). Lower-level assignments at the genus and species level still pose significant challenge for most existing algorithms, and sparse eukaryotic reference databases further limit software performance. This study illuminates current best practices for metabarcoding taxonomy assignments, and underscores the need for community-driven efforts to expand taxonomic and geographic representation in reference DNA barcode databases.

The Future of DNA Barcoding: Reflections from Early Career Researchers

Over the last two decades, the use of DNA barcodes has transformed our ability to identify and assess life on our planet. Both strengths and weaknesses of the method have been exemplified through thousands of peer-reviewed scientific articles. Given the novel sequencing approaches, currently capable of generating millions of reads at low cost, we reflect on the questions: What will the future bring for DNA barcoding? Will identification of species using short, standardized fragments of DNA stand the test of time? We present reflected opinions of early career biodiversity researchers in the form of a SWOT analysis and discuss answers to these questions.

Biases in bulk: DNA metabarcoding of marine communities and the methodology involved

With the growing anthropogenic pressure on marine ecosystems, the need for efficient monitoring of biodiversity grows stronger. DNA metabarcoding of bulk samples is increasingly being implemented in ecosystem assessments and is more cost-efficient and less time-consuming than monitoring based on morphology. However, before raw sequences are obtained from bulk samples, a profound number of methodological choices must be made. Here, we critically review the recent methods used for metabarcoding of marine bulk samples (including benthic, plankton and diet samples) and indicate how potential biases can be introduced throughout sampling, preprocessing, DNA extraction, marker and primer selection, PCR amplification and sequencing. From a total of 64 studies evaluated, our recommendations for best practices include to (a) consider DESS as a fixative instead of ethanol, (b) use the DNeasy PowerSoil kit for any samples containing traces of sediment, (c) not limit the marker selection to COI only, but preferably include multiple markers for higher taxonomic resolution, (d) avoid touchdown PCR profiles, (e) use a fixed annealing temperature for each primer pair when comparing across studies or institutes, (f) use a minimum of three PCR replicates, and (g) include both negative and positive controls. Although the implementation of DNA metabarcoding still faces several technical complexities, we foresee wide-ranging advances in the near future, including improved bioinformatics for taxonomic assignment, sequencing of longer fragments and the use of whole-genome information. Despite the bulk of biases involved in metabarcoding of bulk samples, if appropriate controls are included along the data generation process, it is clear that DNA metabarcoding provides a valuable tool in ecosystem assessments.

Biodiversity at the global scale: the synthesis continues

Traditionally, the generation and use of biodiversity data and their associated specimen objects have been primarily the purview of individuals and small research groups. While deposition of data and specimens in herbaria and other repositories has long been the norm, throughout most of their history, these resources have been accessible only to a small community of specialists. Through recent concerted efforts, primarily at the level of national and international governmental agencies over the last two decades, the pace of biodiversity data accumulation has accelerated, and a wider array of biodiversity scientists has gained access to this massive accumulation of resources, applying them to an ever-widening compass of research pursuits. We review how these new resources and increasing access to them are affecting the landscape of biodiversity research in plants today, focusing on new applications across evolution, ecology, and other fields that have been enabled specifically by the availability of these data and the global scope that was previously beyond the reach of individual investigators. We give an overview of recent advances organized along three lines: broad-scale analyses of distributional data and spatial information, phylogenetic research circumscribing large clades with comprehensive taxon sampling, and data sets derived from improved accessibility of biodiversity literature. We also review synergies between large data resources and more traditional data collection paradigms, describe shortfalls and how to overcome them, and reflect on the future of plant biodiversity analyses in light of increasing linkages between data types and scientists in our field.

Profile hidden Markov model sequence analysis can help remove putative pseudogenes from DNA barcoding and metabarcoding datasets

BACKGROUND

Pseudogenes are non-functional copies of protein coding genes that typically follow a different molecular evolutionary path as compared to functional genes. The inclusion of pseudogene sequences in DNA barcoding and metabarcoding analysis can lead to misleading results. None of the most widely used bioinformatic pipelines used to process marker gene (metabarcode) high throughput sequencing data specifically accounts for the presence of pseudogenes in protein-coding marker genes. The purpose of this study is to develop a method to screen for nuclear mitochondrial DNA segments (nuMTs) in large COI datasets. We do this by: (1) describing gene and nuMT characteristics from an artificial COI barcode dataset, (2) show the impact of two different pseudogene removal methods on perturbed community datasets with simulated nuMTs, and (3) incorporate a pseudogene filtering step in a bioinformatic pipeline that can be used to process Illumina paired-end COI metabarcode sequences. Open reading frame length and sequence bit scores from hidden Markov model (HMM) profile analysis were used to detect pseudogenes.

RESULTS

Our simulations showed that it was more difficult to identify nuMTs from shorter amplicon sequences such as those typically used in metabarcoding compared with full length DNA barcodes that are used in the construction of barcode libraries. It was also more difficult to identify nuMTs in datasets where there is a high percentage of nuMTs. Existing bioinformatic pipelines used to process metabarcode sequences already remove some nuMTs, especially in the rare sequence removal step, but the addition of a pseudogene filtering step can remove up to 5% of sequences even when other filtering steps are in place.

CONCLUSIONS

Open reading frame length filtering alone or combined with hidden Markov model profile analysis can be used to effectively screen out apparent pseudogenes from large datasets. There is more to learn from COI nuMTs such as their frequency in DNA barcoding and metabarcoding studies, their taxonomic distribution, and evolution. Thus, we encourage the submission of verified COI nuMTs to public databases to facilitate future studies.

High-throughput sequencing on preservative ethanol is effective at jointly examining infraspecific and taxonomic diversity, although bioinformatics pipelines do not perform equally

High-throughput sequencing of amplicons (HTSA) has been proposed as an effective approach to evaluate taxonomic and genetic diversity at the same time. However, there are still uncertainties as to how the results produced by different bioinformatics treatments impact the conclusions drawn on biodiversity and population genetics indices.We evaluated the ability of six bioinformatics pipelines to recover taxonomic and genetic diversity from HTSA data obtained from controlled assemblages. To that end, 20 assemblages were produced using 354 colonies of Botrylloides spp., sampled in the wild in ten marinas around Brittany (France). We used DNA extracted from preservative ethanol (ebDNA) after various time of storage (3, 6, and 12 months), and from a bulk of preserved specimens (bulkDNA). DNA was amplified with primers designed for targeting this ascidian genus. Results obtained from HTSA data were compared with Sanger sequencing on individual zooids (i.e., individual barcoding).Species identification and relative abundance determined with HTSA data from either ebDNA or bulkDNA were similar to those obtained with traditional individual barcoding. However, after 12 months of storage, the correlation between HTSA and individual-based data was lower than after shorter durations. The six bioinformatics pipelines were able to depict accurately the genetic diversity using standard population genetics indices (HS and FST), despite producing false positives and missing rare haplotypes. However, they did not perform equally and dada2 was the only pipeline able to retrieve all expected haplotypes.This study showed that ebDNA is a nondestructive alternative for both species identification and haplotype recovery, providing storage does not last more than 6 months before DNA extraction. Choosing the bioinformatics pipeline is a matter of compromise, aiming to retrieve all true haplotypes while avoiding false positives. We here recommend to process HTSA data using dada2, including a chimera-removal step. Even if the possibility to use multiplexed primer sets deserves further investigation to expand the taxonomic coverage in future similar studies, we showed that primers targeting a particular genus allowed to reliably analyze this genus within a complex community.

DNA barcode analysis of the endangered green turtle (Chelonia mydas) in Mexico1

Technological and analytical advances to study evolutionary biology, ecology, and conservation of green turtles (Chelonia mydas) are realized through molecular approaches including DNA barcoding. We characterized the usefulness of COI DNA barcodes in green turtles in Mexico to better understand genetic divergence and other genetic parameters of this species. We analyzed 63 sequences, including 25 from green turtle field specimens collected from the Gulf of Mexico and from the Mexican Pacific and 38 already present in the Barcode of Life Data Systems (BOLD). A total of 13 haplotypes were identified with four novel haplotypes from the Pacific Ocean and three novel haplotypes from the Atlantic Ocean. Intraspecific distance values among COI gene sequences by two different models were 0.01, demonstrating that there is not a subdivision for green turtle species. Otherwise, the interspecific distance interval ranged from 0.07 to 0.13, supporting a clear subdivision among all sea turtle species. Haplotype and total nucleotide diversity values of the COI gene reflect a medium genetic diversity average. Green turtles of the Mexican Pacific showed common haplotypes to some Australian and Chinese turtles, but different from the haplotypes of the Mexican Atlantic. COI analysis revealed new haplotypes and confirmed that DNA barcodes were useful for evaluation of the population diversity of green turtles in Mexico.

Molecular and morphological revision of small Myotinae from the Himalayas shed new light on the poorly known genus Submyotodon (Chiroptera: Vespertilionidae)

The systematics status of the constituent species of the M. mystacinus morphogroup in the Himalayan region has long been marred by uncertainty. Lack of integrative studies combining morphological and genetic data from specimens recently collected in this region has hampered our understanding of cryptic variations in this complex taxonomic group. To address this issue, new material from the Himalayan region of India and Nepal was obtained and vouchered specimens in the holdings of various museums were also re-examined. As comparative material, a large series of relevant specimens from South and Southeast Asia were also included in this revision. Using a combination of multivariate analysis of craniodental characters and molecular reconstructions, we critically evaluated the systematic position of the small Myotinae in the Himalayas. We establish that M. nipalensis forms a very distinct lineage (which also includes the recently described M. annatessae) and refute previous taxonomic suggestions that it is related to M. davidii. Our study also conclusively proved the common occurrence of the poorly known genus Submyotodon in the Himalayan region (Afghanistan, Pakistan, India, Nepal and China) and evidenced species-level divergences within that genus. Submyotodon species share nyctalodont or semi-nyctalodont lower molar configuration with few other small and unrelated Myotinae from Asia suggesting that these unusual dental characters are homoplasious in this subfamily. We also noticed a very confused taxonomic situation associated with many DNA sequences of Asian Myotis deposited in public repositories and call for possibilities of better data curation.

Molecular identification of Brachygenys and Haemulon species (Perciformes: Haemulidae) from the Brazilian coast

Abstract The fishes of the Haemulidae family are currently allocated to 19 genera with a worldwide distribution in the tropical and subtropical waters of the world’s oceans. Brachygenys and Haemulon are important genera of reef fish in Brazil, as they occur in large shoals, which are both ecologically and commercially valuable. This study identified the Brazilian species of the genera Brachygenys and Haemulon using DNA barcodes. While we found only a single lineage in Brachygenys chrysargyrea, Haemulon melanurum, H. parra, and H. squamipinna, more than one molecular operational taxonomic unit (MOTU) was identified in H. atlanticus, H. aurolineatum, and H. plumieri, indicating the possible existence of discrete populations or cryptic species.

Efficient COI barcoding using high throughput single-end 400 bp sequencing

BACKGROUND

Over the last decade, the rapid development of high-throughput sequencing platforms has accelerated species description and assisted morphological classification through DNA barcoding. However, the current high-throughput DNA barcoding methods cannot obtain full-length barcode sequences due to read length limitations (e.g. a maximum read length of 300 bp for the Illumina's MiSeq system), or are hindered by a relatively high cost or low sequencing output (e.g. a maximum number of eight million reads per cell for the PacBio's SEQUEL II system).

RESULTS

Pooled cytochrome c oxidase subunit I (COI) barcodes from individual specimens were sequenced on the MGISEQ-2000 platform using the single-end 400 bp (SE400) module. We present a bioinformatic pipeline, HIFI-SE, that takes reads generated from the 5' and 3' ends of the COI barcode region and assembles them into full-length barcodes. HIFI-SE is written in Python and includes four function modules of filter, assign, assembly and taxonomy. We applied the HIFI-SE to a set of 845 samples (30 marine invertebrates, 815 insects) and delivered a total of 747 fully assembled COI barcodes as well as 70 Wolbachia and fungi symbionts. Compared to their corresponding Sanger sequences (72 sequences available), nearly all samples (71/72) were correctly and accurately assembled, including 46 samples that had a similarity score of 100% and 25 of ca. 99%.

CONCLUSIONS

The HIFI-SE pipeline represents an efficient way to produce standard full-length barcodes, while the reasonable cost and high sensitivity of our method can contribute considerably more DNA barcodes under the same budget. Our method thereby advances DNA-based species identification from diverse ecosystems and increases the number of relevant applications.

Metabarcoding From Microbes to Mammals: Comprehensive Bioassessment on a Global Scale

Global biodiversity loss is unprecedented, and threats to existing biodiversity are growing. Given pervasive global change, a major challenge facing resource managers is a lack of scalable tools to rapidly and consistently measure Earth's biodiversity. Environmental genomic tools provide some hope in the face of this crisis, and DNA metabarcoding, in particular, is a powerful approach for biodiversity assessment at large spatial scales. However, metabarcoding studies are variable in their taxonomic, temporal, or spatial scope, investigating individual species, specific taxonomic groups, or targeted communities at local or regional scales. With the advent of modern, ultra-high throughput sequencing platforms, conducting deep sequencing metabarcoding surveys with multiple DNA markers will enhance the breadth of biodiversity coverage, enabling comprehensive, rapid bioassessment of all the organisms in a sample. Here, we report on a systematic literature review of 1,563 articles published about DNA metabarcoding and summarize how this approach is rapidly revolutionizing global bioassessment efforts. Specifically, we quantify the stakeholders using DNA metabarcoding, the dominant applications of this technology, and the taxonomic groups assessed in these studies. We show that while DNA metabarcoding has reached global coverage, few studies deliver on its promise of near-comprehensive biodiversity assessment. We then outline how DNA metabarcoding can help us move toward real-time, global bioassessment, illustrating how different stakeholders could benefit from DNA metabarcoding. Next, we address barriers to widespread adoption of DNA metabarcoding, highlighting the need for standardized sampling protocols, experts and computational resources to handle the deluge of genomic data, and standardized, open-source bioinformatic pipelines. Finally, we explore how technological and scientific advances will realize the promise of total biodiversity assessment in a sample—from microbes to mammals—and unlock the rich information genomics exposes, opening new possibilities for merging whole-system DNA metabarcoding with (1) abundance and biomass quantification, (2) advanced modeling, such as species occupancy models, to improve species detection, (3) population genetics, (4) phylogenetics, and (5) food web and functional gene analysis. While many challenges need to be addressed to facilitate widespread adoption of environmental genomic approaches, concurrent scientific and technological advances will usher in methods to supplement existing bioassessment tools reliant on morphological and abiotic data. This expanded toolbox will help ensure that the best tool is used for the job and enable exciting integrative techniques that capitalize on multiple tools. Collectively, these new approaches will aid in addressing the global biodiversity crisis we now face.

Harnessing the power of eDNA metabarcoding for the detection of deep-sea fishes

The deep ocean is the largest biome on Earth and faces increasing anthropogenic pressures from climate change and commercial fisheries. Our ability to sustainably manage this expansive habitat is impeded by our poor understanding of its inhabitants and by the difficulties in surveying and monitoring these areas. Environmental DNA (eDNA) metabarcoding has great potential to improve our understanding of this region and to facilitate monitoring across a broad range of taxa. Here, we evaluate two eDNA sampling protocols and seven primer sets for elucidating fish diversity from deep sea water samples. We found that deep sea water samples (> 1400 m depth) had significantly lower DNA concentrations than surface or mid-depth samples necessitating a refined protocol with a larger sampling volume. We recovered significantly more DNA in large volume water samples (1.5 L) filtered at sea compared to small volume samples (250 mL) held for lab filtration. Furthermore, the number of unique sequences (exact sequence variants; ESVs) recovered per sample was higher in large volume samples. Since the number of ESVs recovered from large volume samples was less variable and consistently high, we recommend the larger volumes when sampling water from the deep ocean. We also identified three primer sets which detected the most fish taxa but recommend using multiple markers due the variability in detection probabilities and taxonomic resolution among fishes for each primer set. Overall, fish diversity results obtained from metabarcoding were comparable to conventional survey methods. While eDNA sampling and processing need be optimized for this unique environment, the results of this study demonstrate that eDNA metabarcoding can facilitate biodiversity surveys in the deep ocean, require less dedicated survey effort per unit identification, and are capable of simultaneously providing valuable information on other taxonomic groups.

BAGS: An automated Barcode, Audit & Grade System for DNA barcode reference libraries

Biodiversity studies greatly benefit from molecular tools, such as DNA metabarcoding, which provides an effective identification tool in biomonitoring and conservation programmes. The accuracy of species-level assignment, and consequent taxonomic coverage, relies on comprehensive DNA barcode reference libraries. The role of these libraries is to support species identification, but accidental errors in the generation of the barcodes may compromise their accuracy. Here, we present an R-based application, Barcode, Audit & Grade System (BAGS) (https://github.com/tadeu95/BAGS), that performs automated auditing and annotation of cytochrome c oxidase subunit I (COI) sequences libraries, for a given taxonomic group of animals, available in the Barcode of Life Data System (BOLD). This is followed by implementing a qualitative ranking system that assigns one of five grades (A to E) to each species in the reference library, according to the attributes of the data and congruency of species names with sequences clustered in barcode index numbers (BINs). Our goal is to allow researchers to obtain the most useful and reliable data, highlighting and segregating records according to their congruency. Different tests were performed to perceive its usefulness and limitations. BAGS fulfils a significant gap in the current landscape of DNA barcoding research tools by quickly screening reference libraries to gauge the congruence status of data and facilitate the triage of ambiguous data for posterior review. Thereby, BAGS has the potential to become a valuable addition in forthcoming DNA metabarcoding studies, in the long term contributing to globally improve the quality and reliability of the public reference libraries.

Drivers of tropical soil invertebrate community composition and richness across tropical secondary forests using DNA metasystematics

Tropical forests are fundamental ecosystems, essential for providing terrestrial primary productivity, global nutrient cycling, and biodiversity. Despite their importance, tropical forests are currently threatened by deforestation and associated activities. Moreover, tropical regions are now mostly represented by secondary forest regrowth, with half of the remaining tropical forests as secondary forest. Soil invertebrates are an important component to the functioning and biodiversity of these soil ecosystems. However, it remains unclear how these past land-use activities and subsequent secondary forest developments have altered the soil invertebrate communities and any potential ecological consequences associated with this. DNA metabarcoding offers an effective approach to rapidly monitor soil invertebrate communities under different land-use practices and within secondary forests. In this study, we used DNA metabarcoding to detect community-based patterns of soil invertebrate composition across a primary forest, a 23-year-old secondary forest, and a 33-year-old secondary forest and the associated soil environmental drivers of the soil invertebrate community structure in the Maquenque National Wildlife Refuge of Costa Rica (MNWR). We also used a species contribution analysis (SIMPER) to determine which soil invertebrate groups may be an indication of these soils reaching a pre-disturbed state such as a primary forest. We found that the soil invertebrate community composition at class, order, family, and ESV level were mostly significantly different across that habitats. We also found that the primary forest had a greater richness of soil invertebrates compared to the 23-year-old and 33-year-old secondary forest. Moreover, a redundancy analysis indicated that soil moisture influenced soil invertebrate community structure and explained up to 22% of the total variation observed in the community composition across the habitats; whereas soil invertebrate richness was structured by soil microbial biomass carbon (C) (Cmic) and explained up to 52% of the invertebrate richness across the primary and secondary forests. Lastly, the SIMPER analysis revealed that Naididae, Entomobryidae, and Elateridae could be important indicators of soil and forest recuperation in the MNWR. This study adds to the increasing evidence that soil invertebrates are intimately linked with the soil microbial biomass carbon (Cmic) and that even after 33 years of natural regrowth of a forest, these land use activities can still have persisting effects on the overall composition and richness of the soil invertebrate communities.

Widespread Torix Rickettsia in New Zealand amphipods and the use of blocking primers to rescue host COI sequences

Endosymbionts and intracellular parasites are common in arthropod hosts. As a consequence, (co)amplification of untargeted bacterial sequences has been occasionally reported as a common problem in DNA barcoding. While identifying amphipod species with universal COI primers, we unexpectedly detected rickettsial endosymbionts belonging to the Torix group. To map the distribution and diversity of Rickettsia species among amphipod hosts, we conducted a nationwide molecular screening of seven families of New Zealand freshwater amphipods. In addition to uncovering a diversity of Torix Rickettsia species across multiple amphipod populations from three different families, our research indicates that: (1) detecting Torix Rickettsia with universal primers is not uncommon, (2) obtaining 'Rickettsia COI sequences' from many host individuals is highly likely when a population is infected, and (3) obtaining 'host COI' may not be possible with a conventional PCR if an individual is infected. Because Rickettsia COI is highly conserved across diverse host taxa, we were able to design blocking primers that can be used in a wide range of host species infected with Torix Rickettsia. We propose the use of blocking primers to circumvent problems caused by unwanted amplification of Rickettsia and to obtain targeted host COI sequences for DNA barcoding, population genetics, and phylogeographic studies.

The mitochondrial genomes of 11 aquatic macroinvertebrate species from Cyprus

Aquatic macroinvertebrates are often identified, based on morphology, but molecular approaches like DNA barcoding, metabarcoding and metagenomics are increasingly used for species identification. These approaches require the availability of DNA references deposited in public databases. Here we report the mitochondrial genomes of 11 aquatic macroinvertebrate species from Cyprus, a European Union island country in the Mediterranean. Only three species could be provisionally assigned to a binomial species name, highlighting the current lack of molecular references for aquatic macroinvertebrates from Cyprus.

Graphical Abstract

A total crapshoot? Evaluating bioinformatic decisions in animal diet metabarcoding analyses

Metabarcoding studies provide a powerful approach to estimate the diversity and abundance of organisms in mixed communities in nature. While strategies exist for optimizing sample and sequence library preparation, best practices for bioinformatic processing of amplicon sequence data are lacking in animal diet studies. Here we evaluate how decisions made in core bioinformatic processes, including sequence filtering, database design, and classification, can influence animal metabarcoding results. We show that denoising methods have lower error rates compared to traditional clustering methods, although these differences are largely mitigated by removing low-abundance sequence variants. We also found that available reference datasets from GenBank and BOLD for the animal marker gene cytochrome oxidase I (COI) can be complementary, and we discuss methods to improve existing databases to include versioned releases. Taxonomic classification methods can dramatically affect results. For example, the commonly used Barcode of Life Database (BOLD) Classification API assigned fewer names to samples from order through species levels using both a mock community and bat guano samples compared to all other classifiers (vsearch-SINTAX and q2-feature-classifier's BLAST + LCA, VSEARCH + LCA, and Naive Bayes classifiers). The lack of consensus on bioinformatics best practices limits comparisons among studies and may introduce biases. Our work suggests that biological mock communities offer a useful standard to evaluate the myriad computational decisions impacting animal metabarcoding accuracy. Further, these comparisons highlight the need for continual evaluations as new tools are adopted to ensure that the inferences drawn reflect meaningful biology instead of digital artifacts.

Italian odonates in the Pandora's box: A comprehensive DNA barcoding inventory shows taxonomic warnings at the Holarctic scale

The Odonata are considered among the most endangered freshwater faunal taxa. Their DNA-based monitoring relies on validated reference data sets that are often lacking or do not cover important biogeographical centres of diversification. This study presents the results of a DNA barcoding campaign on Odonata, based on the standard 658-bp 5' end region of the mitochondrial COI gene, involving the collection of 812 specimens (409 of which barcoded) from peninsular Italy and its main islands (328 localities), belonging to all the 88 species (31 Zygoptera and 57 Anisoptera) known from the country. Additional BOLD and GenBank data from Holarctic samples expanded the data set to 1,294 DNA barcodes. A multi-approach species delimitation analysis involving two distance (OT and ABGD) and four tree-based (PTP, MPTP, GMYC and bGMYC) methods was used to explore these data. Of the 88 investigated morphospecies, 75 (85%) unequivocally corresponded to distinct molecular operational units, whereas the remaining ones were classified as 'warnings' (i.e. showing a mismatch between morphospecies assignment and DNA-based species delimitation). These results are in contrast with other DNA barcoding studies on Odonata showing up to 95% of identification success. The species causing warnings were grouped into three categories depending on if they showed low, high or mixed genetic divergence patterns. The analysis of haplotype networks revealed unexpected intraspecific complexity at the Italian, Palearctic and Holarctic scale, possibly indicating the occurrence of cryptic species. Overall, this study provides new insights into the taxonomy of odonates and a valuable basis for future DNA and eDNA-based monitoring studies.

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

Identifying Leafhopper Targets for Controlling Aster Yellows in Carrots and Celery

Aster yellows phytoplasma (Candidatus Phytoplasma asteris) is a multi-host plant pathogen and is transmitted by at least 24 leafhopper species. Pathogen management is complex and requires a thorough understanding of vector dynamics. In the American Midwest, aster yellows is of great concern for vegetable farmers who focus on controlling one vector, Macrosteles quadrilineatus-the aster leafhopper. However, vegetable-associated leafhopper communities can be diverse. To investigate whether additional species are important aster yellows vectors, we surveyed leafhopper communities at commercial celery and carrot farms in Michigan from 2018 to 2019 and conducted real-time PCR to determine infection status. Leafhoppers were collected within crop fields and field edges and identified with DNA barcoding. Overall, we collected 5049 leafhoppers, with the most abundant species being M. quadrilineatus (57%) and Empoasca fabae-the potato leafhopper (23%). Our results revealed the most abundant aster yellows vector in Michigan in both crops is M. quadrilineatus, but we also found that E. fabae may be a potential vector for this pathogen. While several taxa reside in and near these crops, we did not find strong evidence that they contribute to phytoplasma infection. These findings indicate that M. quadrilineatus should be the primary target for controlling this pathogen.

coil: an R package for cytochrome c oxidase I (COI) DNA barcode data cleaning, translation, and error evaluation

Biological conclusions based on DNA barcoding and metabarcoding analyses can be strongly influenced by the methods utilized for data generation and curation, leading to varying levels of success in the separation of biological variation from experimental error. The 5′ region of cytochrome c oxidase subunit I (COI-5P) is the most common barcode gene for animals, with conserved structure and function that allows for biologically informed error identification. Here, we present coil ( https://CRAN.R-project.org/package=coil ), an R package for the pre-processing and frameshift error assessment of COI-5P animal barcode and metabarcode sequence data. The package contains functions for placement of barcodes into a common reading frame, accurate translation of sequences to amino acids, and highlighting insertion and deletion errors. The analysis of 10 000 barcode sequences of varying quality demonstrated how coil can place barcode sequences in reading frame and distinguish sequences containing indel errors from error-free sequences with greater than 97.5% accuracy. Package limitations were tested through the analysis of COI-5P sequences from the plant and fungal kingdoms as well as the analysis of potential contaminants: nuclear mitochondrial pseudogenes and Wolbachia COI-5P sequences. Results demonstrated that coil is a strong technical error identification method but is not reliable for detecting all biological contaminants.

DNA metabarcoding reveals metacommunity dynamics in a threatened boreal wetland wilderness

The complexity and natural variability of ecosystems present a challenge for reliable detection of change due to anthropogenic influences. This issue is exacerbated by necessary trade-offs that reduce the quality and resolution of survey data for assessments at large scales. The Peace-Athabasca Delta (PAD) is a large inland wetland complex in northern Alberta, Canada. Despite its geographic isolation, the PAD is threatened by encroachment of oil sands mining in the Athabasca watershed and hydroelectric dams in the Peace watershed. Methods capable of reliably detecting changes in ecosystem health are needed to evaluate and manage risks. Between 2011 and 2016, aquatic macroinvertebrates were sampled across a gradient of wetland flood frequency, applying both microscope-based morphological identification and DNA metabarcoding. By using multispecies occupancy models, we demonstrate that DNA metabarcoding detected a much broader range of taxa and more taxa per sample compared to traditional morphological identification and was essential to identifying significant responses to flood and thermal regimes. We show that family-level occupancy masks high variation among genera and quantify the bias of barcoding primers on the probability of detection in a natural community. Interestingly, patterns of community assembly were nearly random, suggesting a strong role of stochasticity in the dynamics of the metacommunity. This variability seriously compromises effective monitoring at local scales but also reflects resilience to hydrological and thermal variability. Nevertheless, simulations showed the greater efficiency of metabarcoding, particularly at a finer taxonomic resolution, provided the statistical power needed to detect change at the landscape scale.

From metabarcoding to metaphylogeography: separating the wheat from the chaff

Metabarcoding is by now a well-established method for biodiversity assessment in terrestrial, freshwater, and marine environments. Metabarcoding data sets are usually used for α- and β-diversity estimates, that is, interspecies (or inter-MOTU [molecular operational taxonomic unit]) patterns. However, the use of hypervariable metabarcoding markers may provide an enormous amount of intraspecies (intra-MOTU) information-mostly untapped so far. The use of cytochrome oxidase (COI) amplicons is gaining momentum in metabarcoding studies targeting eukaryote richness. COI has been for a long time the marker of choice in population genetics and phylogeographic studies. Therefore, COI metabarcoding data sets may be used to study intraspecies patterns and phylogeographic features for hundreds of species simultaneously, opening a new field that we suggest to name metaphylogeography. The main challenge for the implementation of this approach is the separation of erroneous sequences from true intra-MOTU variation. Here, we develop a cleaning protocol based on changes in entropy of the different codon positions of the COI sequence, together with co-occurrence patterns of sequences. Using a data set of community DNA from several benthic littoral communities in the Mediterranean and Atlantic seas, we first tested by simulation on a subset of sequences a two-step cleaning approach consisting of a denoising step followed by a minimal abundance filtering. The procedure was then applied to the whole data set. We obtained a total of 563 MOTUs that were usable for phylogeographic inference. We used semiquantitative rank data instead of read abundances to perform AMOVAs and haplotype networks. Genetic variability was mainly concentrated within samples, but with an important between seas component as well. There were intergroup differences in the amount of variability between and within communities in each sea. For two species, the results could be compared with traditional Sanger sequence data available for the same zones, giving similar patterns. Our study shows that metabarcoding data can be used to infer intra- and interpopulation genetic variability of many species at a time, providing a new method with great potential for basic biogeography, connectivity and dispersal studies, and for the more applied fields of conservation genetics, invasion genetics, and design of protected areas.