Sorting things out: Assessing effects of unequal specimen biomass on DNA metabarcoding

Environmental gradients and optimal fixation time revealed with DNA metabarcoding of benthic sample fixative

Assessments of biodiversity and ecosystem status can benefit from DNA metabarcoding as a means to streamline sample processing and specimen identification. Moreover, processing the fixation medium instead of the precious material introduces straightforward protocols that allow subsequent focus on certain organisms detected among the preserved specimens. In this study, we present a proof of concept via the analysis of freshwater invertebrate samples from the Tatra Mountain lakes (Slovakia). Besides highlighting a match between the lake-specific environmental conditions and the results of our fixative DNA metabarcoding, we observed an option to fine-tune the fixation time: to prefer two weeks over a day or a month. This effect emerged from the presence/absence of individual taxa rather than from coarse per-sample records of taxonomic richness, demonstrating that metabarcoding studies-and efforts to optimize their protocols-can use the robust metrics to explore even subtle trends. We also provide evidence that fixative DNA might better capture large freshwater species than terrestrial or meiofauna.

DNA and RNA Stability of Marine Microalgae in Cold-Stored Sediments and Its Implications in Metabarcoding Analyses

The ever-increasing applications of metabarcoding analyses for environmental samples demand a well-designed assessment of the stability of DNA and RNA contained in cells that are deposited or buried in marine sediments. We thus conducted a qPCR quantification of the DNA and RNA in the vegetative cells of three microalgae entrapped in facsimile marine sediments and found that >90% of DNA and up to 99% of RNA for all microalgal species were degraded within 60 days at 4 °C. A further examination of the potential interference of the relic DNA of the vegetative cells with resting cyst detection in sediments was performed via a metabarcoding analysis in artificial marine sediments spiked with the vegetative cells of two Kareniaceae dinoflagellates and the resting cysts of another three dinoflagellates. The results demonstrated a dramatic decrease in the relative abundances of the two Kareniaceae dinoflagellates in 120 days, while those of the three resting cysts increased dramatically. Together, our results suggest that a positive detection of microalgae via metabarcoding analysis in DNA or RNA extracted from marine sediments strongly indicates the presence of intact or viable cysts or spores due to the rapid decay of relic DNA/RNA. This study provides a solid basis for the data interpretation of metabarcoding surveys, particularly in resting cyst detection.

Black gold rush - Evaluating the efficiency of the Fractionator in separating Hymenoptera families in a meadow ecosystem over a two week period

In the face of insect decline, monitoring projects are launched widely to assess trends of insect populations. Collecting over long time periods results in large numbers of samples with thousands of individuals that are often just stored in freezers waiting to be further processed. As the time-consuming process of sorting and identifying specimens prevents taxonomists from working on mass samples, important information on species composition remains unknown and taxonomically neglected species remain undiagnosed. Size fractioning of bulk samples can improve sample handling and, thus, can help to overcome the taxonomic impediment. In this paper, we evaluate the efficiency of the fractionator in separating Hymenoptera families from a Malaise trap sample of a meadow ecosystem over a two week interval to make them available for further morphological identification. The fractionator system by Buffington and Gates (2008) was used to separate the sample in two size classes - a large (macro) and a small (micro) fraction - and Hymenoptera specimens were then counted and identified on family level. In total, 2,449 Hymenoptera specimens were found in the macro fraction and 3,016 in the micro fraction (5,465 specimens in total). For 24 out of 34 Hymenoptera families (71%), separation was significant. This study illustrates the efficiency of the fractionator and its potential to improve workflows dealing with specimen-rich Malaise trap samples.

THAPBI PICT-a fast, cautious, and accurate metabarcoding analysis pipeline

THAPBI PICT is an open source software pipeline for metabarcoding analysis of Illumina paired-end reads, including cases of multiplexing where more than one amplicon is amplified per DNA sample. Initially a Phytophthora ITS1 Classification Tool (PICT), we demonstrate using worked examples with our own and public data sets how, with appropriate primer settings and a custom database, it can be applied to other amplicons and organisms, and used for reanalysis of existing datasets. The core dataflow of the implementation is (i) data reduction to unique marker sequences, often called amplicon sequence variants (ASVs), (ii) dynamic thresholds for discarding low abundance sequences to remove noise and artifacts (rather than error correction by default), before (iii) classification using a curated reference database. The default classifier assigns a label to each query sequence based on a database match that is either perfect, or a single base pair edit away (substitution, deletion or insertion). Abundance thresholds for inclusion can be set by the user or automatically using per-batch negative or synthetic control samples. Output is designed for practical interpretation by non-specialists and includes a read report (ASVs with classification and counts per sample), sample report (samples with counts per species classification), and a topological graph of ASVs as nodes with short edit distances as edges. Source code available from https://github.com/peterjc/thapbi-pict/ with documentation including installation instructions.

Herbivorous dietary selection shown by hawfinch (Coccothraustes coccothraustes) within mixed woodland habitats

Knowledge of diet and dietary selectivity is vital, especially for the conservation of declining species. Accurately obtaining this information, however, is difficult, especially if the study species feeds on a wide range of food items within heterogeneous and inaccessible environments, such as the tree canopy. Hawfinches (Coccothraustes coccothraustes), like many woodland birds, are declining for reasons that are unclear. We investigated the possible role that dietary selection may have in these declines in the UK. Here, we used a combination of high-throughput sequencing of 261 hawfinch faecal samples assessed against tree occurrence data from quadrats sampled in three hawfinch population strongholds in the UK to test for evidence of selective foraging. This revealed that hawfinches show selective feeding and consume certain tree genera disproportionally to availability. Positive selection was shown for beech (Fagus), cherry (Prunus), hornbeam (Carpinus), maples (Acer) and oak (Quercus), while Hawfinch avoided ash (Fraxinus), birch (Betula), chestnut (Castanea), fir (Abies), hazel (Corylus), rowan (Sorbus) and lime (Tilia). This approach provided detailed information on hawfinch dietary choice and may be used to predict the effects of changing food resources on other declining passerines populations in the future.

Monitoring insect biodiversity and comparison of sampling strategies using metabarcoding: A case study in the Yanshan Mountains, China

Insects are the richest and most diverse group of animals and yet there remains a lack, not only of systematic research into their distribution across some key regions of the planet, but of standardized sampling strategies for their study. The Yanshan Mountains, being the boundary range between the Inner Mongolian Plateau and the North China Plain, present an indispensable piece of the insect biodiversity puzzle: both requiring systematic study and offering opportunities for the development of standardized methodologies. This is the first use of DNA metabarcoding to survey the insect biodiversity of the Yanshan Mountains. The study focuses on differences of community composition among samples collected via different methods and from different habitat types. In total, 74 bulk samples were collected from five habitat types (scrubland, woodland, wetland, farmland and grassland) using three collection methods (sweep netting, Malaise traps and light traps). After DNA extraction, PCR amplification, sequencing and diversity analysis were performed, a total of 7427 Operational Taxonomic Units (OTUs) at ≥97% sequence similarity level were delimited, of which 7083 OTUs were identified as belonging to Insecta. Orthoptera, Diptera, Coleoptera and Hemiptera were found to be the dominant orders according to community composition analysis. Nonmetric multidimensional scaling (NMDS) analysis based on Bray-Curtis distances revealed highly divergent estimates of insect community composition among samples differentiated by the collection method (R = .524802, p = .001), but nonsignificant difference among samples differentiated according to habitat (R = .051102, p = .078). The study therefore appears to indicate that the concurrent use of varied collection methods is essential to the accurate monitoring of insect biodiversity.

Metabarcoding reveals high diversity of benthic foraminifera linked to water masses circulation at coastal Svalbard

Arctic marine biodiversity is undergoing rapid changes due to global warming and modifications of oceanic water masses circulation. These changes have been demonstrated in the case of mega- and macrofauna, but much less is known about their impact on the biodiversity of smaller size organisms, such as foraminifera that represent a main component of meiofauna in the Arctic. Several studies analyzed the distribution and diversity of Arctic foraminifera. However, all these studies are based exclusively on the morphological identification of specimens sorted from sediment samples. Here, we present the first assessment of Arctic foraminifera diversity based on metabarcoding of sediment DNA samples collected in fjords and open sea areas in the Svalbard Archipelago. We obtained a total of 5,968,786 reads that represented 1384 amplicon sequence variants (ASVs). More than half of the ASVs (51.7%) could not be assigned to any group in the reference database suggesting a high genetic novelty of Svalbard foraminifera. The sieved and unsieved samples resolved comparable communities, sharing 1023 ASVs, comprising over 97% of reads. Our analyses show that the foraminiferal assemblage differs between the localities, with communities distinctly separated between fjord and open sea stations. Each locality was characterized by a specific assemblage, with only a small overlap in the case of open sea areas. Our study demonstrates a clear pattern of the influence of water masses on the structure of foraminiferal communities. The stations situated on the western coast of Svalbard that are strongly influenced by warm and salty Atlantic water (AW) are characterized by much higher diversity than stations in the northern and eastern part, where the impact of AW is less pronounced. This high diversity and specificity of Svalbard foraminifera associated with water mass distribution indicate that the foraminiferal metabarcoding data can be very useful for inferring present and past environmental conditions in the Arctic.

Designing a surveillance program for early detection of alien plants and insects in Norway

Naturalized species of alien plants and animals comprise < 3% of biodiversity recorded in Norway but have had major impacts on natural ecosystems through displacement of native species. Encroachment of alien species has been especially problematic for coastal sites close to transport facilities and urban areas with high density housing. The goal of our field project was to design and test a surveillance program for early detection of alien species of vascular plants and terrestrial insects at the first phase of establishment in natural areas. In our 3-year project (2018–2020), we sampled 60 study plots in three counties in the Oslofjord region of southern Norway. Study plots (6.25 ha) were selected by two criteria: manual selection based on expert opinion (27 plots) or by random selection based on weights from a hotspot analysis of occurrence of alien species (33 plots). Vascular plants were surveyed by two experienced botanists who found a total of 239 alien species of vascular plants in 95 rounds of surveys. Insects and other invertebrates were captured with a single Malaise trap per site, with 3–4 rounds of repeated sampling. We used DNA-metabarcoding to identify invertebrates based on DNA extractions from crushed insects or from the preservative media. Over 3500 invertebrate taxa were detected in 255 rounds of sampling. We recorded 20 alien species of known risk, and 115 species that were new to Norway, including several ‘doorknocker’ species identified by previous risk assessments. We modeled the probabilities of occupancy (ψ) and detection (p) with occupancy models with repeated visits by multiple observers (vascular plants) or multiple rounds of sampling (insects). The two probabilities covaried with risk category for alien organisms and both were low for species categorized as no known or low risk (range = 0.052–0.326) but were higher for species categorized as severe risk (range = 0.318–0.651). Selecting sites at random or manually did not improve the probability of finding novel alien species, but occupancy had a weak positive relationship with housing density for some categories of alien plants and insects. We used our empirical estimates to test alternative sampling designs that would minimize the combined variance of occupancy and detection (A-optimality criterion). Sampling designs with 8–10 visits per site were best for surveillance of new alien species if the probabilities of occupancy and detection were both low, and provided low conditional probabilities of site occupancy ($$\hat{\psi }_{condl}$$ ψ ^ condl ≤ 0.032) and a high probabilities of cumulative detection ($$\hat{p}*$$ p ^ ∗  ≥ 0.943). Our field results demonstrate that early detection is feasible as a key component of a national surveillance program based on early detection and rapid response.

Discovering marine biodiversity in the 21st century

We review the current knowledge of the biodiversity of the ocean as well as the levels of decline and threat for species and habitats. The lack of understanding of the distribution of life in the ocean is identified as a significant barrier to restoring its biodiversity and health. We explore why the science of taxonomy has failed to deliver knowledge of what species are present in the ocean, how they are distributed and how they are responding to global and regional to local anthropogenic pressures. This failure prevents nations from meeting their international commitments to conserve marine biodiversity with the results that investment in taxonomy has declined in many countries. We explore a range of new technologies and approaches for discovery of marine species and their detection and monitoring. These include: imaging methods, molecular approaches, active and passive acoustics, the use of interconnected databases and citizen science. Whilst no one method is suitable for discovering or detecting all groups of organisms many are complementary and have been combined to give a more complete picture of biodiversity in marine ecosystems. We conclude that integrated approaches represent the best way forwards for accelerating species discovery, description and biodiversity assessment. Examples of integrated taxonomic approaches are identified from terrestrial ecosystems. Such integrated taxonomic approaches require the adoption of cybertaxonomy approaches and will be boosted by new autonomous sampling platforms and development of machine-speed exchange of digital information between databases.

Fluctuating insect diversity, abundance and biomass across agricultural landscapes

Habitat destruction and deterioration of habitat quality caused a severe decline of biodiversity, such as insect diversity. In this study, we analyze insect diversity and biomass across agro-environments. We collected flying insects with 20 malaise traps across a landscape mosaic consisting of organic (eight traps) and conventional (four traps) farmland, as well as across agricultural land that has been recently converted from conventional to organic farming (eight traps). Sampling was conducted over 2 years, in 2019 and 2020, with in total 340 sampling events. We measured the dry weight of the captured organisms and identified species diversity by analyzing Operational Taxonomic Units (OTUs) and Barcode Index Numbers (BINs) via metabarcoding. The results obtained show temporal dynamics. The number of OTUs were always higher than the number of BINs. OTUs and BINs were moderately to highly correlated, while the number of OTUs and BINs were only moderately positively correlated with dry biomass. OTUs and BINs as well as biomass were highest in the recently transformed farmland if compared with pure organic and conventional farmland sites, which showed no significant differences in respect of insect diversity. OTU and BIN numbers but not the OTU/BIN ratio significantly decreased with increasing distance from the nearest forest fringe. The numbers of OTUs, BINs and the OTU/BIN proportion, as well as OTU and BIN/biomass proportions varied strongly over seasons, irrespective of agricultural practice. Based on our findings, we suggest to combine data on insect species richness and biomass measured over a period of time, to derive a largely complete and meaningful assessment of biodiversity for a specific region.

Toward global integration of biodiversity big data: a harmonized metabarcode data generation module for terrestrial arthropods

Metazoan metabarcoding is emerging as an essential strategy for inventorying biodiversity, with diverse projects currently generating massive quantities of community-level data. The potential for integrating across such data sets offers new opportunities to better understand biodiversity and how it might respond to global change. However, large-scale syntheses may be compromised if metabarcoding workflows differ from each other. There are ongoing efforts to improve standardization for the reporting of inventory data. However, harmonization at the stage of generating metabarcode data has yet to be addressed. A modular framework for harmonized data generation offers a pathway to navigate the complex structure of terrestrial metazoan biodiversity. Here, through our collective expertise as practitioners, method developers, and researchers leading metabarcoding initiatives to inventory terrestrial biodiversity, we seek to initiate a harmonized framework for metabarcode data generation, with a terrestrial arthropod module. We develop an initial set of submodules covering the 5 main steps of metabarcode data generation: (i) sample acquisition; (ii) sample processing; (iii) DNA extraction; (iv) polymerase chain reaction amplification, library preparation, and sequencing; and (v) DNA sequence and metadata deposition, providing a backbone for a terrestrial arthropod module. To achieve this, we (i) identified key points for harmonization, (ii) reviewed the current state of the art, and (iii) distilled existing knowledge within submodules, thus promoting best practice by providing guidelines and recommendations to reduce the universe of methodological options. We advocate the adoption and further development of the terrestrial arthropod module. We further encourage the development of modules for other biodiversity fractions as an essential step toward large-scale biodiversity synthesis through harmonization.

Evaluation of non-destructive DNA extraction protocols for insect metabarcoding: gentler and shorter is better

DNA metabarcoding can accelerate research on insect diversity, as it is cheap and fast compared to manual sorting and identification. Most metabarcoding protocols require homogenisation of the sample, preventing further work on the specimens. Mild digestion of the tissue by incubation in a lysis buffer has been proposed as an alternative, and, although some mild lysis protocols have already been presented, they have so far not been evaluated against each other. Here, we analyse how two mild lysis buffers (one more aggressive, one gentler in terms of tissue degradation), two different incubation times, and two DNA purification methods (a manual precipitation and an automated protocol) affect the accuracy of retrieving the true composition of mock communities using two mitochondrial markers (COI and 16S). We found that protocol-specific variation in concentration and purity of the DNA extracts produced had little effect on the recovery of species. However, the two lysis treatments differed in quantification of species abundances. Digestion in the gentler buffer and for a shorter time yielded better representation of original sample composition. Digestion in a more aggressive buffer or longer incubation time yielded lower alpha diversity values and increased differences between metabarcoding results and the true species-abundance distribution. We conclude that the details of non-destructive protocols can have a significant effect on metabarcoding performance. A short and mild lysis treatment appears the best choice for recovering the true composition of the sample. This not only improves accuracy, but also comes with a faster processing time than the other treatments.

Speeding up the detection of invasive bivalve species using environmental DNA: a Nanopore and Illumina sequencing comparison

Traditional detection of aquatic invasive species via morphological identification is often time-consuming and can require a high level of taxonomic expertise, leading to delayed mitigation responses. Environmental DNA (eDNA) detection approaches of multiple species using Illumina-based sequencing technology have been used to overcome these hindrances, but sample processing is often lengthy. More recently, portable nanopore sequencing technology has become available, which has the potential to make molecular detection of invasive species more widely accessible and substantially decrease sample turnaround times. However, nanopore-sequenced reads have a much higher error rate than those produced by Illumina platforms, which has so far hindered the adoption of this technology. We provide a detailed laboratory protocol and bioinformatic tools (msi package) to increase the reliability of nanopore sequencing to detect invasive species, and we test its application using invasive bivalves while comparing it with Illumina-based sequencing. We sampled water from sites with pre-existing bivalve occurrence and abundance data, and contrasting bivalve communities, in Italy and Portugal. Samples were extracted, amplified, and sequenced by the two platforms. The mean agreement between sequencing methods was 69% and the difference between methods was non-significant. The lack of detections of some species at some sites could be explained by their known low abundances. This is the first reported use of MinION to detect aquatic invasive species from eDNA samples.

Disentangling bias for non-destructive insect metabarcoding

A fast and reliable method for obtaining a species-level identification is a fundamental requirement for a wide range of activities, from plant protection and invasive species management to biodiversity assessments and ecological studies. For insects, novel molecular techniques such as DNA metabarcoding have emerged as a rapid alternative to traditional morphological identification, reducing the dependence on limited taxonomic experts. Until recently, molecular techniques have required a destructive DNA extraction, precluding the possibility of preserving voucher specimens for future studies, or species descriptions. Here we paired insect metabarcoding with two recent non-destructive DNA extraction protocols, to obtain a rapid and high-throughput taxonomic identification of diverse insect taxa while retaining a physical voucher specimen. The aim of this work was to explore how non-destructive extraction protocols impact the semi-quantitative nature of metabarcoding, which alongside species presence/absence also provides a quantitative, but biased, representation of their relative abundances. By using a series of mock communities representing each stage of a typical metabarcoding workflow we were able to determine how different morphological (i.e., insect biomass and exoskeleton hardness) and molecular traits (i.e., primer mismatch and amplicon GC%), interact with different protocol steps to introduce quantitative bias into non-destructive metabarcoding results. We discuss the relevance of taxonomic bias to metabarcoding identification of insects and potential approaches to account for it.

Aquatic Ecological Risk of Heavy-Metal Pollution Associated with Degraded Mining Landscapes of the Southern Africa River Basins: A Review

Africa accounts for nearly 30% of the discovered world’s mineral reserves, with half of the world’s platinum group metals deposits, 36% of gold, and 20% of cobalt being in Southern Africa (SA). The intensification of heavy-metal production in the SA region has exacerbated negative human and environmental health impacts. In recent years, mining waste generated from industrial and artisanal mining has significantly affected the ecological integrity of SA aquatic ecosystems due to the accelerated introduction and deposition of heavy metals. However, the extent to which heavy-metal pollution associated with mining has impacted the aquatic ecosystems has not been adequately documented, particularly during bioassessments. This review explores the current aquatic ecological impacts on the heavily mined river basins of SA. It also discusses the approaches to assessing the ecological risks, inherent challenges, and potential for developing an integrated ecological risk assessment protocol for aquatic systems in the region. Progress has been made in developing rapid bioassessment schemes (RBS) for SA aquatic ecosystems. Nevertheless, method integration, which also involves heavy-metal pollution monitoring and molecular technology, is necessary to overcome the current challenges of the standardisation of RBS protocols. Citizenry science will also encourage community and stakeholder involvement in sustainable environmental management in SA.

Unearthing soil arthropod diversity through DNA metabarcoding

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

Evaluating restoration trajectories using DNA metabarcoding of ground-dwelling and airborne invertebrates and associated plant communities

Invertebrates are important for restoration processes as they are key drivers of many landscape‐scale ecosystem functions; including pollination, nutrient cycling and soil formation. However, invertebrates are often overlooked in restoration monitoring because they are highly diverse, poorly described, and time‐consuming to survey, and require increasingly scarce taxonomic expertise to enable identification. DNA metabarcoding is a relatively new tool for rapid survey that is able to address some of these concerns, and provide information about the taxa with which invertebrates are interacting via food webs and habitat. Here, we evaluate how invertebrate communities may be used to determine ecosystem trajectories during restoration. We collected ground‐dwelling and airborne invertebrates across chronosequences of mine‐site restoration in three ecologically disparate locations in Western Australia and identified invertebrate and plant communities using DNA metabarcoding. Ground‐dwelling invertebrates showed the clearest restoration signals, with communities becoming more similar to reference communities over time. These patterns were weaker in airborne invertebrates, which have higher dispersal abilities and therefore less local fidelity to environmental conditions. Although we detected directional changes in community composition indicative of invertebrate recovery, patterns observed were inconsistent between study locations. The inclusion of plant assays allowed identification of plant species, as well as potential food sources and habitat. We demonstrate that DNA metabarcoding of invertebrate communities can be used to evaluate restoration trajectories. Testing and incorporating new monitoring techniques such as DNA metabarcoding is critical to improving restoration outcomes.

Measuring mass: variation among 3,161 species of Canadian Coleoptera and the prospects of a mass registry for all insects

Although biomass values are critical for diverse ecological and evolutionary analyses, they are unavailable for most insect species. Museum specimens have the potential to address this gap, but the variation introduced by sampling and preservation methods is uncertain. This study quantifies species-level variation in the body mass of Canadian Coleoptera based on the analysis of 3,744 specimens representing 3,161 Barcode Index Number (BIN) clusters. Employing the BIN system as a proxy for species allows the inclusion of groups where the taxonomic impediment prevents the assignment of specimens to a Linnaean species. By validating the reproducibility of measurements and evaluating the error introduced by operational complexities such as curatorial practice and the loss of body parts, this study demonstrates that museum specimens can speed the assembly of a mass registry. The results further indicate that congeneric species of Coleoptera generally have limited variation in mass, so a genus-level identification allows prediction of the body mass of species that have not been weighed or measured. Building on the present results, the construction of a mass registry for all insects is feasible.

Climate-induced forest dieback drives compositional changes in insect communities that are more pronounced for rare species

Species richness, abundance and biomass of insects have recently undergone marked declines in Europe. We metabarcoded 211 Malaise-trap samples to investigate whether drought-induced forest dieback and subsequent salvage logging had an impact on ca. 3000 species of flying insects in silver fir Pyrenean forests. While forest dieback had no measurable impact on species richness, there were significant changes in community composition that were consistent with those observed during natural forest succession. Importantly, most observed changes were driven by rare species. Variation was explained primarily by canopy openness at the local scale, and the tree-related microhabitat diversity and deadwood amount at landscape scales. The levels of salvage logging in our study did not explain compositional changes. We conclude that forest dieback drives changes in species assemblages that mimic natural forest succession, and markedly increases the risk of catastrophic loss of rare species through homogenization of environmental conditions.

Metabarcoding reveals massive species diversity of Diptera in a subtropical ecosystem

Diptera is often considered to be the richest insect group due to its great species diversity and broad ecological versatility. However, data on dipteran diversity from subtropical ecosystems have hitherto been scarce, due to the lack of studies conducted at an appropriate large scale. We investigated the diversity and composition of Diptera communities on Tianmu Mountain, Zhejiang, China, using DNA metabarcoding technology, and evaluated their dynamic responses to the effects of slope aspect, season, and altitudinal zone. A total of 5,092 operational taxonomic units (OTUs) were discovered and tentatively assigned to 72 dipteran families, including 2 family records new for China and 30 family records new for the locality. Cecidomyiidae, Sciaridae, and Phoridae were the predominant families, representing 53.6% of total OTUs, while 52 families include >95% unidentified and presumed undescribed species. We found that the community structure of Diptera was significantly affected by aspect, seasonality (month) and elevation, with richer diversity harbored in north-facing than south-facing slopes, and seasonality a more profound driver of community structure and diversity than elevation. Overall, massive species diversity of Diptera communities was discovered in this subtropical ecosystem of east China. The huge diversity of potentially undescribed species only revealed by metabarcoding now requires more detailed taxonomic study, as a step toward an evolutionary integration that accumulates information on species' geographic ranges, ecological traits, functional roles, and species interactions, and thus places the local communities in the context of the growing knowledge base of global biodiversity and its response to environmental change.

A comparison of genetically and morphometrically identified macroinvertebrate community index scores with implications for aquatic life use attainment

Measures of aquatic macroinvertebrate communities are often used to characterize water quality and indicate whether waterbodies are meeting management expectations. The accuracy of these measures depends on the skill and experience of the person identifying the macroinvertebrates, and obtaining these measures can be relatively expensive due to the time necessary for identification. Utilizing genetic identification of macroinvertebrate taxa has the potential to reduce the time of sample processing, identify a greater number of taxa, and increase the resolution of identification. We compared Colorado multi-metric index (MMI) scores from seven locations in the Big Thompson River, CO, based on genetic and morphometric identification and estimated the ability of MMI scores based on genetic identification to characterize aquatic life use attainment management thresholds. We found a significant linear relationship (p = 0.002, R² = 0.87) between MMI scores generated by genetic and morphological identification. MMIs support the following aquatic life use designations as defined by the Colorado Water Quality Control Commission: Impaired < 40, Attaining > 48, and Ambiguous 40-48. These values correspond to MMIs based on genetic identification as Impaired < 20, Attaining > 64, and Ambiguous = 21-63 based on the prediction interval of the regression equation. Our results suggest that using genetically identified macroinvertebrates to estimate MMI scores can provide some degree of certainty regarding aquatic life use designations, and while it may be inappropriate at the current time to entirely replace morphologically based biotic integrity measures with those based on molecular identification, there are opportunities in their use.

Semi-quantitative metabarcoding reveals how climate shapes arthropod community assembly along elevation gradients on Hawaii Island

Spatial variation in climatic conditions along elevation gradients provides an important backdrop by which communities assemble and diversify. Lowland habitats tend to be connected through time, whereas highlands can be continuously or periodically isolated, conditions that have been hypothesized to promote high levels of species endemism. This tendency is expected to be accentuated among taxa that show niche conservatism within a given climatic envelope. While species distribution modeling approaches have allowed extensive exploration of niche conservatism among target taxa, a broad understanding of the phenomenon requires sampling of entire communities. Species-rich groups such as arthropods are ideal case studies for understanding ecological and biodiversity dynamics along elevational gradients given their important functional role in many ecosystems, but community-level studies have been limited due to their tremendous diversity. Here, we develop a novel semi-quantitative metabarcoding approach that combines specimen counts and size-sorting to characterize arthropod community-level diversity patterns along elevational transects on two different volcanoes of the island of Hawai'i. We found that arthropod communities between the two transects became increasingly distinct compositionally at higher elevations. Resistance surface approaches suggest that climatic differences between sampling localities are an important driver in shaping beta-diversity patterns, though the relative importance of climate varies across taxonomic groups. Nevertheless, the climatic niche position of OTUs between transects was highly correlated, suggesting that climatic filters shape the colonization between adjacent volcanoes. Taken together, our results highlight climatic niche conservatism as an important factor shaping ecological assembly along elevational gradients and suggest topographic complexity as an important driver of diversification.

Assessing changes in stream macroinvertebrate communities across ecological gradients using morphological versus DNA metabarcoding approaches

Freshwater macroinvertebrates provide valuable indicators for biomonitoring ecosystem change in relation to natural and anthropogenic drivers. DNA metabarcoding is an efficient approach for estimating such indicators, but its results may differ from morphotaxonomic approaches traditionally used in biomonitoring. Here we test the hypothesis that despite differences in the number and identity of taxa recorded, both approaches may retrieve comparable patterns of community change, and detect similar ecological gradients influencing such changes. We compared results obtained with morphological identification at family level of macroinvertebrates collected at 80 streams under a Water Framework Directive biomonitoring program in Portugal, with results obtained with metabarcoding from the ethanol preserving the bulk samples, using either single (COI-M19BR2, 16S-Inse01, 18S-Euka02) or multiple markers. Metabarcoding recorded less families and different communities compared to morphotaxonomy, but community sensitivities to disturbance estimated with the IASPT index were more similar across approaches. Spatial variation in local community metrics and the factors influencing such variation were significantly correlated between morphotaxonomy and metabarcoding. After reducing random noise in the dissimilarity matrices, the spatial variation in community composition was also significantly correlated across methods. A dominant gradient of community change was consistently retrieved, and all methods identified a largely similar set of anthropogenic stressors strongly influencing such gradient. Overall, results confirm our initial hypothesis, suggesting that morphotaxonomy and metabarcoding can estimate consistent spatial patterns of community variation and their main drivers. These results are encouraging for macroinvertebrate biomonitoring using metabarcoding approaches, suggesting that they can be intercalibrated with morphotaxonomic approaches to recover equivalent spatial and temporal gradients of ecological change.

Pooling size sorted Malaise trap fractions to maximize taxon recovery with metabarcoding

Background

Small and rare specimens can remain undetected when metabarcoding is applied on bulk samples with a high specimen size heterogeneity. This is especially critical for Malaise trap samples, where most of the biodiversity is contributed by small taxa with low biomass. The separation of samples in different size fractions for downstream analysis is one possibility to increase detection of small and rare taxa. However, experiments systematically testing different size sorting approaches and subsequent proportional pooling of fractions are lacking, but would provide important information for the optimization of metabarcoding protocols. We set out to find a size sorting strategy for Malaise trap samples that maximizes taxonomic recovery but remains scalable and time efficient.

Methods

Three Malaise trap samples were sorted into four size classes using dry sieving. Each fraction was homogenized and lysed. The corresponding lysates were pooled to simulate unsorted samples. Pooling was additionally conducted in equal proportions and in four different proportions enriching the small size fraction of samples. DNA from the individual size classes as well as the pooled fractions was extracted and metabarcoded using the FwhF2 and Fol-degen-rev primer set. Additionally, alternative wet sieving strategies were explored.

Results

The small size fractions harboured the highest diversity and were best represented when pooling in favour of small specimens. Metabarcoding of unsorted samples decreases taxon recovery compared to size sorted samples. A size separation into only two fractions (below 4 mm and above) can double taxon recovery compared to not size sorting. However, increasing the sequencing depth 3- to 4-fold can also increase taxon recovery to levels comparable with size sorting, but remains biased towards biomass rich taxa in the sample.

Conclusion

We demonstrate that size fractionation of Malaise trap bulk samples can increase taxon recovery. While results show distinct patterns, the lack of statistical support due to the limited number of samples processed is a limitation. Due to increased speed and lower risk of cross-contamination as well as specimen damage we recommend wet sieving and proportional pooling of the lysates in favour of the small size fraction (80–90% volume). However, for large-scale projects with time constraints, increasing sequencing depth is an alternative solution.

Quantification of marine benthic communities with metabarcoding

DNA metabarcoding methods have been implemented in studies aimed at detecting and quantifying marine benthic biodiversity. In such surveys, universal barcodes are amplified and sequenced from environmental DNA. To quantify biodiversity with DNA metabarcoding, a relation between the number of DNA sequences of a species and its biomass and/or the abundance is required. However, this relationship is complicated by many factors, and it is often unknown. In this study, we validate estimates of biomass and abundance from molecular approaches with those from the traditional morphological approach. Abundance and biomass were quantified from 126 samples of benthic intertidal mudflat using traditional morphological approaches and compared with frequency of occurrence and relative read abundance estimates from a molecular approach. A relationship between biomass and relative read abundance was found for two widely dispersed annelid taxa (Pygospio and Scoloplos). None of the other taxons, however, showed such a relationship. We discuss how quantification of abundance and biomass using molecular approaches are hampered by the ecology of DNA i.e. all the processes that determine the amount of DNA in the environment, including the ecology of the benthic species as well as the compositional nature of sequencing data.

Relationship of insect biomass and richness with land use along a climate gradient

Recently reported insect declines have raised both political and social concern. Although the declines have been attributed to land use and climate change, supporting evidence suffers from low taxonomic resolution, short time series, a focus on local scales, and the collinearity of the identified drivers. In this study, we conducted a systematic assessment of insect populations in southern Germany, which showed that differences in insect biomass and richness are highly context dependent. We found the largest difference in biomass between semi-natural and urban environments (-42%), whereas differences in total richness (-29%) and the richness of threatened species (-56%) were largest from semi-natural to agricultural environments. These results point to urbanization and agriculture as major drivers of decline. We also found that richness and biomass increase monotonously with increasing temperature, independent of habitat. The contrasting patterns of insect biomass and richness question the use of these indicators as mutual surrogates. Our study provides support for the implementation of more comprehensive measures aimed at habitat restoration in order to halt insect declines.

Standardized high-throughput biomonitoring using DNA metabarcoding: Strategies for the adoption of automated liquid handlers

Reliable and comprehensive monitoring data are required to trace and counteract biodiversity loss. High-throughput metabarcoding using DNA extracted from community samples (bulk) or from water or sediment (environmental DNA) has revolutionized biomonitoring, given the capability to assess biodiversity across the tree of life rapidly with feasible effort and at a modest price. DNA metabarcoding can be upscaled to process hundreds of samples in parallel. However, while automated high-throughput analysis workflows are well-established in the medical sector, manual sample processing still predominates in biomonitoring laboratory workflows limiting the upscaling and standardization for routine monitoring applications. Here we present an automated, scalable, and reproducible metabarcoding workflow to extract DNA from bulk samples, perform PCR and library preparation on a liquid handler. Key features are the independent sample replication throughout the workflow and the use of many negative controls for quality assurance and quality control. We generated two datasets: i) a validation dataset consisting of 42 individual arthropod specimens of different species, and ii) a routine monitoring dataset consisting of 60 stream macroinvertebrate bulk samples. As a marker, we used the mitochondrial COI gene. Our results show that the developed single-deck workflow is free of laboratory-derived contamination and produces highly consistent results. Minor deviations between replicates are mostly due to stochastic differences for low abundant OTUs. Thus, we successfully demonstrated that robotic liquid handling can be used reliably from DNA extraction to final library preparation on a single deck, thereby substantially increasing throughput, reducing costs, and increasing data robustness for biodiversity assessments and monitoring.

Highly comparable metabarcoding results from MGI-Tech and Illumina sequencing platforms

With the developments in DNA nanoball sequencing technologies and the emergence of new platforms, there is an increasing interest in their performance in comparison with the widely used sequencing-by-synthesis methods. Here, we test the consistency of metabarcoding results from DNBSEQ-G400RS (DNA nanoball sequencing platform by MGI-Tech) and NovaSeq 6000 (sequencing-by-synthesis platform by Illumina) platforms using technical replicates of DNA libraries that consist of COI gene amplicons from 120 soil DNA samples. By subjecting raw sequencing data from both platforms to a uniform bioinformatics processing, we found that the proportion of high-quality reads passing through the filtering steps was similar in both datasets. Per-sample operational taxonomic unit (OTU) and amplicon sequence variant (ASV) richness patterns were highly correlated, but sequencing data from DNBSEQ-G400RS harbored a higher number of OTUs. This may be related to the lower dominance of most common OTUs in DNBSEQ data set (thus revealing higher richness by detecting rare taxa) and/or to a lower effective read quality leading to generation of spurious OTUs. However, there was no statistical difference in the ASV and post-clustered ASV richness between platforms, suggesting that additional denoising step in the ASV workflow had effectively removed the 'noisy' reads. Both OTU-based and ASV-based composition were strongly correlated between the sequencing platforms, with essentially interchangeable results. Therefore, we conclude that DNBSEQ-G400RS and NovaSeq 6000 are both equally efficient high-throughput sequencing platforms to be utilized in studies aiming to apply the metabarcoding approach, but the main benefit of the former is related to lower sequencing cost.

Dispersal limitations and long-term persistence drive differentiation from haplotypes to communities within a tropical sky-island: Evidence from community metabarcoding

Neutral theory proposes that dispersal stochasticity is one of the main drivers of local diversity. Haplotypes-level genetic variation can now be efficiently sampled from across whole communities, thus making it possible to test neutral predictions from the genetic to species-level diversity, and higher. However, empirical data is still limited, with the few studies to date coming from temperate latitudes. Here, we focus on a tropical mountain within the Transmexican Volcanic Belt to evaluate spatially fine-scale patterns of arthropod community assembly to understand the role of dispersal limitation and landscape features as drivers of diversity. We sampled whole-communities of arthropods for eight orders at a spatial scale ranging from 50 m to 19 km, using whole community metabarcoding. We explored multiple hierarchical levels, from individual haplotypes to lineages at 0.5, 1.5, 3, 5, and 7.5% similarity thresholds, to evaluate patterns of richness, turnover, and distance decay of similarity with isolation-by-distance and isolation-by-resistance (costs to dispersal given by landscape features) approaches. Our results showed that distance and altitude influence distance decay of similarity at all hierarchical levels. This holds for arthropod groups of contrasting dispersal abilities, but with different strength depending on the spatial scale. Our results support a model where local-scale differentiation mediated by dispersal constraints, combined with long-term persistence of lineages, is an important driver of diversity within tropical sky islands.

Effects of consumer surface sterilization on diet DNA metabarcoding data of terrestrial invertebrates in natural environments and feeding trials

DNA metabarcoding is an emerging tool used to quantify diet in environments and consumer groups where traditional approaches are unviable, including small-bodied invertebrate taxa. However, metabarcoding of small taxa often requires DNA extraction from full body parts (without dissection), and it is unclear whether surface contamination from body parts alters presumed diet presence or diversity.We examined four different measures of diet (presence, rarefied read abundance, richness, and species composition) for a terrestrial invertebrate consumer (the spider Heteropoda venatoria) both collected in its natural environment and fed an offered diet item in contained feeding trials using DNA metabarcoding of full body parts (opisthosomas). We compared diet from consumer individuals surface sterilized to remove contaminants in 10% commercial bleach solution followed by deionized water with a set of unsterilized individuals.We found that surface sterilization did not significantly alter any measure of diet for consumers in either a natural environment or feeding trials. The best-fitting model predicting diet detection in feeding trial consumers included surface sterilization, but this term was not statistically significant (β = -2.3, p-value = .07).Our results suggest that surface contamination does not seem to be a significant concern in this DNA diet metabarcoding study for consumers in either a natural terrestrial environment or feeding trials. As the field of diet DNA metabarcoding continues to progress into new environmental contexts with various molecular approaches, we suggest ongoing context-specific consideration of the possibility of surface contamination.

Issues of under-representation in quantitative DNA metabarcoding weaken the inference about diet of the tundra vole Microtus oeconomus

During the last decade, methods based on high-throughput sequencing such as DNA metabarcoding have opened up for a range of new questions in animal dietary studies. One of the major advantages of dietary metabarcoding resides in the potential to infer a quantitative relationship between sequence read proportions and biomass of ingested food. However, this relationship's robustness is highly dependent on the system under study, calling for case-specific assessments. Herbivorous small rodents often play important roles in the ecosystem, and the use of DNA metabarcoding for analyses of rodent diets is increasing. However, there has been no direct validation of the quantitative reliability of DNA metabarcoding for small rodents. Therefore, we used an experimental approach to assess the relationship between input plant biomass and sequence reads proportions from DNA metabarcoding in the tundra vole Microtus oeconomus. We found a weakly positive relationship between the number of high-throughput DNA sequences and the expected biomass proportions of food plants. The weak relationship was possibly caused by a systematic under-amplification of one of the three plant taxa fed. Generally, our results add to the growing evidence that case-specific validation studies are required to reliably make use of sequence read abundance as a proxy of relative food proportions in the diet.

The power of metabarcoding: Can we improve bioassessment and biodiversity surveys of stream macroinvertebrate communities?

Most stream bioassessment and biodiversity surveys are currently based on morphological identification of communities. However, DNA metabarcoding is emerging as a fast and cost-effective alternative for species identification. We compared both methods in a survey of benthic macroinvertebrate communities across 36 stream sites in northern Finland. We identified 291 taxa of which 62% were identified only by DNA metabarcoding. DNA metabarcoding produced extensive species level inventories for groups (Oligochaeta, Chironomidae, Simuliidae, Limoniidae and Limnephilidae), for which morphological identification was not feasible due to the high level of expertise needed. Metabarcoding also provided more insightful taxonomic information on the occurrence of three red-listed vulnerable or data deficient species, the discovery of two likely cryptic and potentially new species to Finland and species information of insect genera at an early larval stage that could not be separated morphologically. However, it systematically failed to reliably detect the occurrence of gastropods that were easily identified morphologically. The impact of mining on community structure could only be shown using DNA metabarcoding data which suggests that the finer taxonomic detail can improve detection of subtle impacts. Both methods generally exhibited similar strength of community-environment relationships, but DNA metabarcoding showed better performance with presence/absence data than with relative DNA sequence abundances. Our results suggest that DNA metabarcoding holds a promise for future anthropogenic impact assessments, although, in our case, the performance did not improve much from the morphological species identification. The key advantage of DNA metabarcoding lies in efficient biodiversity surveys, taxonomical studies and applications in conservation biology.

A bright idea-metabarcoding arthropods from light fixtures

Arthropod communities in buildings have not been extensively studied, although humans have always shared their homes with them. In this study we explored if arthropod DNA can be retrieved and metabarcoded from indoor environments through the collection of dead specimens in light fixtures to better understand what shapes arthropod diversity in our homes. Insects were collected from 45 light fixtures at the Centre for Biodiversity Genomics (CBG, Guelph, Canada), and by community scientists at 12 different residential homes in Southern Ontario. The CBG ground floor of the CBG showed the greatest arthropod diversity, especially in light fixtures that were continuously illuminated. The community scientist samples varied strongly by light fixture type, lightbulb used, time passed since lamp was last cleaned, and specimen size. In all cases, the majority of OTUs was not shared between samples even within the same building. This study demonstrates that light fixtures might be a useful resource to determine arthropod diversity in our homes, but individual samples are likely not representative of the full diversity.

Testing small-scale ecological gradients and intraspecific differentiation for hundreds of kelp forest species using haplotypes from metabarcoding

DNA metabarcoding has been increasingly used to detail distributions of hundreds of species. Most analyses focus on creating molecular operational taxonomic units (MOTUs) from complex mixtures of DNA sequences, but much less common is use of the sequence diversity within these MOTUs. Here we use the diversity of COI haplotypes within MOTUs from a California kelp forest to infer patterns of population abundance, dispersal and population history from 527 species of animals and algae from 106 samples of benthic habitats in Monterey Bay. Using haplotypes as a unit we show fine-grained differences of abundance across locations for 15 species, and marked aggregation from sample to sample for most of the common species of plants and animals. Previous analyses could not distinguish these patterns from artefacts of amplification or sequence bias. Our haplotype data also reveal strong population genetic differentiation over small spatial scales for 48 species of red algae, sponges and Bryozoa. Last, phylogenetic analysis of mismatch frequencies among haplotypes show a wide variety of demographic histories from recent expansions to long, stable population sizes. These analyses show that abundant, small-bodied marine species that are often overlooked in ecological surveys can have strikingly different patterns of ecological and genetic structure leading to population, ecological and perhaps adaptive differences between habitats. MOTU diversity data from the same sequencing efforts that generate species-level analyses can greatly increase the scope and value of metabarcoding studies.

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.

Unearthing the Potential of Soil eDNA Metabarcoding—Towards Best Practice Advice for Invertebrate Biodiversity Assessment

Metabarcoding has proven to be a powerful tool to assess ecological patterns and diversity from different habitats. Terrestrial invertebrate diversity is frequently based on bulk samples, which require comparatively high sampling effort. With environmental DNA (eDNA) metabarcoding, field sampling effort can be reduced while increasing the number of recovered organism groups. However, a proof of concept is missing for several invertebrate groups, hampering the development of best-practice advice for these groups. This study aims to provide recommendations on key aspects for the processing of soil samples, from sampling effort to choice of DNA extraction method and marker genes. This study uses eDNA metabarcoding as a tool for assessing invertebrate biodiversity in soil samples, specifically comparing two DNA extraction methods (with and without a lysis step) and two genes, 18S and COI markers. The results show that the choice of marker and DNA extraction method (including a lysis step) significantly affect species detection rates and concomitantly observed invertebrate community composition. Combining methods, by using larger amounts of starting material and including a lysis step resulted in an increase of invertebrate species numbers. Together, these methods improved the detection of species with known lower population densities and allowed the assessment of temporary mesofauna. Furthermore, the choice of marker significantly influenced the diversity levels found. The 18S marker allowed the detection of a higher number of annelid and nematode OTUs, while the COI marker was more suitable for detecting changes in arthropod community structure, especially at the species level. This study makes significant advances to the field of invertebrate biodiversity assessment, particularly using metabarcoding tools by addressing several methodological considerations that are key for accurate ecological appraisals.

Metabarcoding Malaise traps and soil eDNA reveals seasonal and local arthropod diversity shifts

Forest habitats host enormous diversity, but little is known about the seasonal turnover of arthropod species between the above- and below ground forest layers. In this study, we used metabarcoding approaches to uncover arthropod diversity in different forest types and seasons. Our study shows that metabarcoding soil eDNA and Malaise trap bulk samples can provide valuable insights into the phenology and life cycles of arthropods. We found major differences in arthropod species diversity between soil samples and Malaise traps, with only 11.8% species overlap. Higher diversity levels were found in Malaise traps in summer whereas soil samples showed a diversity peak in winter, highlighting the seasonal habitat preferences and life strategies of arthropods. We conclude that collecting time series of bulk arthropod samples and eDNA in the same locations provides a more complete picture of local arthropod diversity and turnover rates and may provide valuable information on climate induced phenological shifts for long-term monitoring.

Evaluating sediment and water sampling methods for the estimation of deep-sea biodiversity using environmental DNA

Despite representing one of the largest biomes on earth, biodiversity of the deep seafloor is still poorly known. Environmental DNA metabarcoding offers prospects for fast inventories and surveys, yet requires standardized sampling approaches and careful choice of environmental substrate. Here, we aimed to optimize the genetic assessment of prokaryote (16S), protistan (18S V4), and metazoan (18S V1-V2, COI) communities, by evaluating sampling strategies for sediment and aboveground water, deployed simultaneously at one deep-sea site. For sediment, while size-class sorting through sieving had no significant effect on total detected alpha diversity and resolved similar taxonomic compositions at the phylum level for all markers studied, it effectively increased the detection of meiofauna phyla. For water, large volumes obtained from an in situ pump (~ 6000 L) detected significantly more metazoan diversity than 7.5 L collected in sampling boxes. However, the pump being limited by larger mesh sizes (> 20 µm), only captured a fraction of microbial diversity, while sampling boxes allowed access to the pico- and nanoplankton. More importantly, communities characterized by aboveground water samples significantly differed from those characterized by sediment, whatever volume used, and both sample types only shared between 3 and 8% of molecular units. Together, these results underline that sediment sieving may be recommended when targeting metazoans, and aboveground water does not represent an alternative to sediment sampling for inventories of benthic diversity.

Quantitative assessment of fish larvae community composition in spawning areas using metabarcoding of bulk samples

Accurate assessment of larval community composition in spawning areas is essential for fisheries management and conservation but is often hampered by the cryptic nature of many larvae, which renders them difficult to identify morphologically. Metabarcoding is a rapid and cost-effective method to monitor early life stages for management and environmental impact assessment purposes but its quantitative capability is under discussion. We compared metabarcoding with traditional morphological identification to evaluate taxonomic precision and reliability of abundance estimates, using 332 fish larvae from multinet hauls (0-50 m depth) collected at 14 offshore sampling sites in the Irish and Celtic seas. To improve quantification accuracy (relative abundance estimates), the amount of tissue for each specimen was standardized and mitochondrial primers (12S gene) with conserved binding sites were used. Relative family abundance estimated from metabarcoding reads and morphological assessment were positively correlated, as well as taxon richness (R_S  = 0.81, P = 0.007) and diversity (R_S  = 0.90, P = 0.002). Spatial patterns of community composition did not differ significantly between metabarcoding and morphological assessments. Our results show that DNA metabarcoding of bulk tissue samples can be used to monitor changes in fish larvae abundance and community composition. This represents a feasible, efficient, and faster alternative to morphological methods that can be applied to terrestrial and aquatic habitats.

TaxonTableTools: A comprehensive, platform-independent graphical user interface software to explore and visualise DNA metabarcoding data

DNA metabarcoding is increasingly used as a tool to assess biodiversity in research and environmental management. Powerful analysis software exists to process raw data. However, the translation of sequence read data into biological information and downstream analyses may be difficult for end users with limited expertise in bioinformatics. Thus, there is a growing need for easy-to-use, graphical user interface (GUI) software to analyse and visualise DNA metabarcoding data. Here, we present TaxonTableTools (TTT), a new platform-independent GUI that aims to fill this gap by providing simple, reproducible analysis and visualisation workflows. At its base, TTT uses a "TaXon table", which is a data format that can be generated easily within TTT from two input files: a read table and a taxonomy table obtained using various published metabarcoding pipelines. TTT analysis and visualisation modules include Venn diagrams to compare taxon overlap among replicates, samples, or analysis methods. TTT analyses and visualises basic statistics, such as read proportion per taxon, as well as more sophisticated visualisations, such as interactive Krona charts for taxonomic data exploration. Various ecological analyses can be produced directly, including alpha or beta diversity estimates, and rarefaction analysis ordination plots. Metabarcoding data can be converted into formats required for traditional, taxonomy-based analyses performed by regulatory bioassessment programs. In addition, TTT is able to produce html-based interactive graphics that can be analysed in any web browser. The software comes with a manual and tutorial, is free and publicly available through GitHub (https://github.com/TillMacher/TaxonTableTools) or the Python package index (https://pypi.org/project/taxontabletools/).

Integration of DNA-Based Approaches in Aquatic Ecological Assessment Using Benthic Macroinvertebrates

Benthic macroinvertebrates are among the most used biological quality elements for assessing the condition of all types of aquatic ecosystems worldwide (i.e., fresh water, transitional, and marine). Current morphology-based assessments have several limitations that may be circumvented by using DNA-based approaches. Here, we present a comprehensive review of 90 publications on the use of DNA metabarcoding of benthic macroinvertebrates in aquatic ecosystems bioassessments. Metabarcoding of bulk macrozoobenthos has been preferentially used in fresh waters, whereas in marine waters, environmental DNA (eDNA) from sediment and bulk communities from deployed artificial structures has been favored. DNA extraction has been done predominantly through commercial kits, and cytochrome c oxidase subunit I (COI) has been, by far, the most used marker, occasionally combined with others, namely, the 18S rRNA gene. Current limitations include the lack of standardized protocols and broad-coverage primers, the incompleteness of reference libraries, and the inability to reliably extrapolate abundance data. In addition, morphology versus DNA benchmarking of ecological status and biotic indexes are required to allow general worldwide implementation and higher end-user confidence. The increased sensitivity, high throughput, and faster execution of DNA metabarcoding can provide much higher spatial and temporal data resolution on aquatic ecological status, thereby being more responsive to immediate management needs.

Adding DNA barcoding to stream monitoring protocols - What's the additional value and congruence between morphological and molecular identification approaches?

Although aquatic macroinvertebrates and freshwater fishes are important indicators for freshwater quality assessments, the morphological identification to species-level is often impossible and thus especially in many invertebrate taxa not mandatory during Water Framework Directive monitoring, a pragmatism that potentially leads to information loss. Here, we focus on the freshwater fauna of the River Sieg (Germany) to test congruence and additional value in taxa detection and taxonomic resolution of DNA barcoding vs. morphology-based identification in monitoring routines. Prior generated morphological identifications of juvenile fishes and aquatic macroinvertebrates were directly compared to species assignments using the identification engine of the Barcode of Life Data System. In 18% of the invertebrates morphology allowed only assignments to higher systematic entities, but DNA barcoding lead to species-level assignment. Dissimilarities between the two approaches occurred in 7% of the invertebrates and in 1% of the fishes. The 18 fish species were assigned to 20 molecular barcode index numbers, the 104 aquatic invertebrate taxa to 113 molecular entities. Although the cost-benefit analysis of both methods showed that DNA barcoding is still more expensive (5.30–8.60€ per sample) and time consuming (12.5h), the results emphasize the potential to increase taxonomic resolution and gain a more complete profile of biodiversity, especially in invertebrates. The provided reference DNA barcodes help building the foundation for metabarcoding approaches, which provide faster sample processing and more cost-efficient ecological status determination.

Assessing the impact of the threatened crucian carp (Carassius carassius) on pond invertebrate diversity: A comparison of conventional and molecular tools

Fishes stocked for recreation and angling can damage freshwater habitats and negatively impact biodiversity. The pond-associated crucian carp (Carassius carassius) is rare across Europe and is stocked for conservation management in England, but its impacts on pond biota are understudied. Freshwater invertebrates contribute substantially to aquatic biodiversity, encompassing many rare and endemic species, but their small size and high abundance complicate their assessment. Practitioners have employed sweep-netting and kick-sampling with microscopy (morphotaxonomy), but specimen size/quality and experience can bias identification. DNA and environmental DNA (eDNA) metabarcoding offer alternative means of invertebrate assessment. We compared invertebrate diversity in ponds (N = 18) with and without crucian carp using morphotaxonomic identification, DNA metabarcoding and eDNA metabarcoding. Five 2 L water samples and 3 min sweep-net samples were collected at each pond. Inventories produced by morphotaxonomic identification of netted samples, DNA metabarcoding of bulk tissue samples and eDNA metabarcoding of water samples were compared. Alpha diversity was greatest with DNA or eDNA metabarcoding, depending on whether standard or unbiased methods were considered. DNA metabarcoding reflected morphotaxonomic identification, whereas eDNA metabarcoding produced markedly different communities. These complementary tools should be combined for comprehensive invertebrate assessment. Crucian carp presence minimally reduced alpha diversity in ponds, but positively influenced beta diversity through taxon turnover (i.e., ponds with crucian carp contained different invertebrates to fishless ponds). Crucian carp presence contributes to landscape-scale invertebrate diversity, supporting continued conservation management in England. Our results show that molecular tools can enhance freshwater invertebrate assessment and facilitate development of more accurate and ecologically effective pond management strategies.

Modelling technical and biological biases in macroinvertebrate community assessment from bulk preservative using multiple metabarcoding markers

DNA metabarcoding from the ethanol used to store macroinvertebrate bulk samples is a convenient methodological option in molecular biodiversity assessment and biomonitoring of aquatic ecosystems, as it preserves specimens and reduces problems associated with sample sorting. However, this method may be affected by errors and biases, which need to be thoroughly quantified before it can be mainstreamed into biomonitoring programmes. Here, we used 80 unsorted macroinvertebrate samples collected in Portugal under a Water Framework Directive monitoring programme, to compare community diversity and taxonomic composition metrics estimated through morphotaxonomy versus metabarcoding from storage ethanol using three markers (COI‐M19BR2, 16S‐Inse01 and 18S‐Euka02) and a multimarker approach. A preliminary in silico analysis showed that the three markers were adequate for the target taxa, with detection failures related primarily to the lack of adequate barcodes in public databases. Metabarcoding of ethanol samples retrieved far less taxa per site (alpha diversity) than morphotaxonomy, albeit with smaller differences for COI‐M19BR2 and the multimarker approach, while estimates of taxa turnover (beta diversity) among sites were similar across methods. Using generalized linear mixed models, we found that after controlling for differences in read coverage across samples, the probability of detection of a taxon was positively related to its proportional abundance, and negatively so to the presence of heavily sclerotized exoskeleton (e.g., Coleoptera). Overall, using our experimental protocol with different template dilutions, the COI marker showed the best performance, but we recommend the use of a multimarker approach to detect a wider range of taxa in freshwater macroinvertebrate samples. Further methodological development and optimization efforts are needed to reduce biases associated with body armouring and rarity in some macroinvertebrate taxa.

Marine biomonitoring with eDNA: Can metabarcoding of water samples cut it as a tool for surveying benthic communities?

In the marine realm, biomonitoring using environmental DNA (eDNA) of benthic communities requires destructive direct sampling or the setting-up of settlement structures. Comparatively much less effort is required to sample the water column, which can be accessed remotely. In this study we assess the feasibility of obtaining information from the eukaryotic benthic communities by sampling the adjacent water layer. We studied two different rocky-substrate benthic communities with a technique based on quadrat sampling. We also took replicate water samples at four distances (0, 0.5, 1.5, and 20 m) from the benthic habitat. Using broad range primers to amplify a ca. 313 bp fragment of the cytochrome oxidase subunit I gene, we obtained a total of 3,543 molecular operational taxonomic units (MOTUs). The structure obtained in the two environments was markedly different, with Metazoa, Archaeplastida and Stramenopiles being the most diverse groups in benthic samples, and Hacrobia, Metazoa and Alveolata in the water. Only 265 MOTUs (7.5%) were shared between benthos and water samples and, of these, 180 (5.1%) were identified as benthic taxa that left their DNA in the water. Most of them were found immediately adjacent to the benthos, and their number decreased as we moved apart from the benthic habitat. It was concluded that water eDNA, even in the close vicinity of the benthos, was a poor proxy for the analysis of benthic structure, and that direct sampling methods are required for monitoring these complex communities via metabarcoding.

Nematodes as bioindicators of polluted sediments using metabarcoding and microscopic taxonomy

The use of bioindicator species is a widely applied approach to evaluate ecological conditions, and several indices have been designed for this purpose. To assess the impact of pollution, especially in sediments, a pollution-sensitive index based on nematodes, one of the most abundant and species-rich groups of metazoa, was developed. The NemaSPEAR[%] index in its original form relies on the morphological inspection of nematode species. The application of a morphologically based NemaSPEAR[%] at the genus-level was previously validated. The present study evaluated a NemaSPEAR[%] index based on metabarcoding of nematode communities and tested the potential of fragments from the 28S rDNA, 18S rDNA and cytochrome c oxidase subunit I (COI) genes. In general, molecular-based results tended to show a poorer condition than morphology-based results for the investigated sites. At the genus level, NemaSPEAR[%] values based on morphological data strongly correlated with those based on molecular data for both the 28S rDNA and the 18S rDNA gene fragments (R² = 0.86 and R² = 0.74, respectively). Within the dominant genera (>3%) identified by morphology, 68% were detected by at least one of the two ribosomal markers. At the species level, however, concordance was less pronounced, as there were several deviations of the molecular from the morphological data. These differences could mostly be attributed to shortcomings in the reference database used in the molecular-based assignments. Our pilot study shows that a molecularly based, genus-level NemaSPEAR[%] can be successfully applied to evaluate polluted sediment. Future studies need to validate this approach further, e.g. with bulk extractions of whole meiofaunal communities in order to circumvent time-consuming nematode isolation. Further database curation with abundant NemaSPEAR[%] species will also increase the applicability of this approach.