Establishing arthropod community composition using metabarcoding: Surprising inconsistencies between soil samples and preservative ethanol and homogenate from Malaise trap catches

Using DNA barcoding to identify high-priority taxa (Hymenoptera: Ichneumonidae) from Great Smoky Mountains National Park

The All Taxa Biodiversity Inventory (ATBI) in Great Smoky Mountains National Park (GSMNP) seeks to document every species of living thing in the park. The ATBI is decades in progress, yet some taxa remain virtually untouched by taxonomists. Such "high priority" taxa include the hyper-diverse parasitoid wasp family Ichneumonidae. Despite the positive and multifaceted effects ichneumonids have on their environment, only a small percentage of those collected in the park have been identified as species, mostly to their complex morphology and overwhelming diversity. Recently, DNA barcoding has transformed biodiversity inventories, streamlining the process to be more rapid and efficient. To test the effectiveness of barcoding 20 + year-old specimens of Ichneumonidae and catalog new records for GSMNP, COI was amplified from 95 ichneumonid morphospecies collected from Andrew's Bald, NC. Species identifications were confirmed morphologically. Eighty-one ichneumonids generated sequence data, representing 16 subfamilies and 44 genera. The subfamily Oxytorinae is newly recorded from GSMNP, along with 10 newly recorded genera and 23 newly recorded species across Ichneumonidae. These results contribute significantly to the ATBI by adding new park records for a high-priority taxon and demonstrate the effectiveness of applying DNA barcoding to samples in long-term storage or those lacking immediate taxonomic expertise.

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.

Effects of livestock on arthropod biodiversity in Iberian holm oak savannas revealed by metabarcoding

Increasing food production while avoiding negative impacts on biodiversity constitutes one of the main challenges of our time. Traditional silvopastoral systems like Iberian oak savannas ("dehesas") set an example, where free-range livestock has been reared for centuries while preserving a high natural value. Nevertheless, factors decreasing productivity need to be addressed, one being acorn losses provoked by pest insects. An increased and focalized grazing by livestock on infested acorns would kill the larvae inside and decrease pest numbers, but increased livestock densities could have undesired side effects on ground arthropod communities as a whole. We designed an experimental setup including areas under trees with livestock exclosures of different ages (short-term: 1-year exclusion, long-term: 10-year exclusion), along with controls (continuous grazing), using DNA metabarcoding (mitochondrial markers COI and 16S) to rapidly assess arthropod communities' composition. Livestock removal quickly increased grass cover and arthropod taxonomic richness and diversity, which was already higher in short-term (1-year exclosures) than beneath the canopies of control trees. Interestingly, arthropod diversity was not highest at long-term exclosures (≥10 years), although their community composition was the most distinct. Also, regardless of treatment, we found that functional diversity strongly correlated with the vegetation structure, being higher at trees beneath which there was higher grass cover and taller herbs. Overall, the taxonomic diversity peak at short term exclosures would support the intermediate disturbance hypothesis, which relates it with the higher microhabitat heterogeneity at moderately disturbed areas. Thus, we propose a rotatory livestock management in dehesas: plots with increased grazing should co-exist with temporal short-term exclosures. Ideally, a few long-term excluded areas should be also kept for the singularity of their arthropod communities. This strategy would make possible the combination of biological pest control and arthropod conservation in Iberian dehesas.

Towards holistic insect monitoring: species discovery, description, identification and traits for all insects

Holistic insect monitoring needs scalable techniques to overcome taxon biases, determine species abundances, and gather functional traits for all species. This requires that we address taxonomic impediments and the paucity of data on abundance, biomass and functional traits. We here outline how these data deficiencies could be addressed at scale. The workflow starts with large-scale barcoding (megabarcoding) of all specimens from mass samples obtained at biomonitoring sites. The barcodes are then used to group the specimens into molecular operational taxonomic units that are subsequently tested/validated as species with a second data source (e.g. morphology). New species are described using barcodes, images and short diagnoses, and abundance data are collected for both new and described species. The specimen images used for species discovery then become the raw material for training artificial intelligence identification algorithms and collecting trait data such as body size, biomass and feeding modes. Additional trait data can be obtained from vouchers by using genomic tools developed by molecular ecologists. Applying this pipeline to a few samples per site will lead to greatly improved insect monitoring regardless of whether the species composition of a sample is determined with images, metabarcoding or megabarcoding. This article is part of the theme issue 'Towards a toolkit for global insect biodiversity monitoring'.

Three steps towards comparability and standardization among molecular methods for characterizing insect communities

Molecular methods are currently some of the best-suited technologies for implementation in insect monitoring. However, the field is developing rapidly and lacks agreement on methodology or community standards. To apply DNA-based methods in large-scale monitoring, and to gain insight across commensurate data, we need easy-to-implement standards that improve data comparability. Here, we provide three recommendations for how to improve and harmonize efforts in biodiversity assessment and monitoring via metabarcoding: (i) we should adopt the use of synthetic spike-ins, which will act as positive controls and internal standards; (ii) we should consider using several markers through a multiplex polymerase chain reaction (PCR) approach; and (iii) we should commit to the publication and transparency of all protocol-associated metadata in a standardized fashion. For (i), we provide a ready-to-use recipe for synthetic cytochrome c oxidase spike-ins, which enable between-sample comparisons. For (ii), we propose two gene regions for the implementation of multiplex PCR approaches, thereby achieving a more comprehensive community description. For (iii), we offer guidelines for transparent and unified reporting of field, wet-laboratory and dry-laboratory procedures, as a key to making comparisons between studies. Together, we feel that these three advances will result in joint quality and calibration standards rather than the current laboratory-specific proof of concepts. This article is part of the theme issue 'Towards a toolkit for global insect biodiversity monitoring'.

From DNA to diagnostics: A case study using macroinvertebrate metabarcoding to assess the effectiveness of restoration measures in a Dutch stream

Stream ecosystems are under pressure due to multiple stressors. Restoration measures can halt further degradation and improve their ecological status. However, assessment of the effectiveness of the implemented measures is often insufficient because of logistic and financial constraints. DNA-metabarcoding has been proposed to scale up sample processing, although its application as a diagnostic tool has received less attention. The aim of our study was to evaluate if DNA-metabarcoding of stream macroinvertebrates can be used to compute a stressor-specific index to assess the effectiveness of a stream restoration project. For this purpose, we sampled the upstream, restored, and downstream section of a recently restored lowland stream in the Netherlands. At each site, we applied three different methods of macroinvertebrate identification: morphological identification of bulk samples (morphology), DNA-metabarcoding of the same bulk samples (DNA) and metabarcoding of eDNA extracted from the water (eDNA). First, we compared the community composition identified by each method. The communities identified by morphology and DNA were highly similar, whereas the communities generated by the eDNA differed. Second, we analysed whether the identification methods could be used to assess the effectiveness of the restoration project, focussing on a stressor-specific index for flow as the restoration measures aimed at improving flow conditions. Both the morphology and bulk DNA samples indicated improved flow conditions in the restored section of the stream (i.e., less stress from the reduction or absence of flow than in the unrestored sections). Contrary, the eDNA-water samples did not differentiate the amount of stress throughout the catchment, although applying recent developments in eDNA sampling could lead to more robust results. In conclusion, this study forms proof of concept that DNA from bulk samples can be utilized to assess the effectiveness of restoration measures, showing the added value of this approach for water managers.

Novel foods, old issues: Metabarcoding revealed mislabeling in insect-based products sold by e-commerce on the EU market

Insects intended for human consumption are considered Novel Foods according to EU legislation. marketed in form of powders, bars, snacks are increasingly available on the EU market, especially on e-commerce. The commercial form and the way of distribution make IBPs particularly prone to mislabeling. Literature concerning the mislabeling occurrence in IBPs is extremely scarce. In this study, 46 processed IBPs were collected on nine EU e-commerce platforms (e-CO) to be authenticated by metabarcoding. A 200 bp region from 16S rRNA gene was used as molecular target. Sequencing data were processed using DADA2 R package, and sequences were taxonomically assigned through BLAST analysis against GenBank. Procedural blanks and positive controls were included in the analysis, and threshold values were established to filter the final data. The mislabeling rate (i. e. the mismatch between the species declared on the IBP label and the species identified by metabarcoding) was calculated. Overall, a high mislabeling rate (33.3 %) was observed, although this percentage is influenced by the e-CO platform and the insect species, with A. domesticus particularly involved. The use of species not listed in authorized Novel Food (e. g. Gryllus locorojo), and/or the partial replacement of high value species with lower value species was highlighted for the first time in processed IBPs. The presence of insect pests was also detected. Metabarcoding was confirmed as an effective tool for IBPs authentication. Also, outcomes from this study can provide useful data on the main issues involving the EU IBPs' market, that can represent an incentive to reinforce both official controls and FBO's self-controls on these poorly investigated products.

Diversity of Parasitoid Wasps and Comparison of Sampling Strategies in Rice Fields Using Metabarcoding

A comprehensive and precise evaluation of Arthropoda diversity in agricultural landscapes can enhance biological pest control strategies. We used Malaise traps and sweep nets to collect insects from three double-cropping paddy fields. DNA was extracted from the ethanol preservative of the Malaise traps and from tissue samples of selected parasitoid wasps. This was followed by amplification using DNA barcoding primers to prepare high-throughput sequencing libraries. We annotated a total of 4956 operational taxonomic units (OTUs), encompassing 174 genera and 32 families of parasitoid wasps. The ethanol filter method efficiently captured a wide range of information. However, the method has low resolution and may result in a reduced estimate of species abundance. Additional insect species were also identified in the parasitoid samples. This suggests that high throughput sequencing from adult parasitoid wasps can also detect host species, enabling a better understanding of host species and providing insights into food webs.

Inferring the ecological and evolutionary determinants of community genetic diversity

Understanding the relative contributions of ecological and evolutionary processes to the structuring of ecological communities is needed to improve our ability to predict how communities may respond to future changes in an increasingly human-modified world. Metabarcoding methods make it possible to gather population genetic data for all species within a community, unlocking a new axis of data to potentially unveil the origins and maintenance of biodiversity at local scales. Here, we present a new eco-evolutionary simulation model for investigating community assembly dynamics using metabarcoding data. The model makes joint predictions of species abundance, genetic variation, trait distributions and phylogenetic relationships under a wide range of parameter settings (e.g. high speciation/low dispersal or vice versa) and across a range of community states, from pristine and unmodified to heavily disturbed. We first demonstrate that parameters governing metacommunity and local community processes leave detectable signatures in simulated biodiversity data axes. Next, using a simulation-based machine learning approach we show that neutral and non-neutral models are distinguishable and that reasonable estimates of several model parameters within the local community can be obtained using only community-scale genetic data, while phylogenetic information is required to estimate those describing metacommunity dynamics. Finally, we apply the model to soil microarthropod metabarcoding data from the Troodos mountains of Cyprus, where we find that communities in widespread forest habitats are structured by neutral processes, while high-elevation and isolated habitats act as an abiotic filter generating non-neutral community structure. We implement our model within the ibiogen R package, a package dedicated to the investigation of island, and more generally community-scale, biodiversity using community-scale genetic data.

Persisting roadblocks in arthropod monitoring using non-destructive metabarcoding from collection media of passive traps

Background

Broad-scale monitoring of arthropods is often carried out with passive traps (e.g., Malaise traps) that can collect thousands of specimens per sample. The identification of individual specimens requires time and taxonomic expertise, limiting the geographical and temporal scale of research and monitoring studies. DNA metabarcoding of bulk-sample homogenates has been found to be faster, efficient and reliable, but the destruction of samples prevents a posteriori validation of species occurrences and relative abundances. Non-destructive metabarcoding of DNA extracted from collection medium has been applied in a limited number of studies, but further tests of efficiency are required with different trap types and collection media to assess the consistency of the method.

Methods

We quantified the detection rate of arthropod species when applying non-destructive DNA metabarcoding with a short (127-bp) fragment of mitochondrial COI on two combinations of passive traps and collection media: (1) water with monopropylene glycol (H2O-MPG) used in window-flight traps (WFT, 53 in total); (2) ethanol with monopropylene glycol (EtOH-MPG) used in Malaise traps (MT, 27 in total). We then compared our results with those obtained for the same samples using morphological identification (for WFTs) or destructive metabarcoding of bulk homogenate (for MTs). This comparison was applied as part of a larger study of arthropod species richness in silver fir (Abies alba Mill., 1759) stands across a range of climate-induced tree dieback levels and forest management strategies.

Results

Of the 53 H2O-MPG samples from WFTs, 16 produced no metabarcoding results, while the remaining 37 samples yielded 77 arthropod MOTUs in total, of which none matched any of the 343 beetle species morphologically identified from the same traps. Metabarcoding of 26 EtOH-MPG samples from MTs detected more arthropod MOTUs (233) than destructive metabarcoding of homogenate (146 MOTUs, 8 orders), of which 71 were shared MOTUs, though MOTU richness per trap was similar between treatments. While we acknowledge the failure of metabarcoding from WFT-derived collection medium (H2O-MPG), the treatment of EtOH-based Malaise trapping medium remains promising. We conclude however that DNA metabarcoding from collection medium still requires further methodological developments and cannot replace homogenate metabarcoding as an approach for arthropod monitoring. It can be used nonetheless as a complementary treatment when enhancing the detection of soft-bodied arthropods like spiders and Diptera.

Decontamination protocols affect the internal microbiota of ticks

Studies on the microbiota of ticks have promoted hypotheses about the combined effects of the bacterial community, its functional contributions to the tick's physiology or probable competition effects with some tick-borne pathogens. However, knowledge on the origin of the microbiota of newly hatched larvae is missing. This study aimed to elucidate the source(s) of the microbiota in unfed tick larvae, addressing the composition of the "core microbiota" and the best ways to decontaminate eggs for microbiota studies. We applied laboratory degree bleach washes and/or ultraviolet light treatments on engorged Rhipicephalus australis females and/or their eggs. No significant effects of these treatments on the reproductive parameters of females and the hatching rates of eggs were observed. However, the different treatments did show striking effects on the composition of the microbiota. The results indicated that bleach washes disrupted the internal tick microbiota in females, implying that bleach may have entered the tick and subsequently affected the microbiota. Furthermore, the analyses of results demonstrated that the ovary is a main source of tick microbiota, while the contribution of Gené's organ (a part of the female reproductive system that secretes a protective wax coat onto tick eggs) or the male's spermatophore requires further investigation. Further studies are needed to identify best practice protocols for the decontamination of ticks for microbiota studies.

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.

Future of DNA-based insect monitoring

Insects are crucial for ecosystem health but climate change and pesticide use are driving massive insect decline. To mitigate this loss, we need new and effective monitoring techniques. Over the past decade there has been a shift to DNA-based techniques. We describe key emerging techniques for sample collection. We suggest that the selection of tools should be broadened, and that DNA-based insect monitoring data need to be integrated more rapidly into policymaking. We argue that there are four key areas for advancement, including the generation of more complete DNA barcode databases to interpret molecular data, standardisation of molecular methods, scaling up of monitoring efforts, and integrating molecular tools with other technologies that allow continuous, passive monitoring based on images and/or laser imaging, detection, and ranging (LIDAR).

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.

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.

Emerging technologies revolutionise insect ecology and monitoring

Insects are the most diverse group of animals on Earth, but their small size and high diversity have always made them challenging to study. Recent technological advances have the potential to revolutionise insect ecology and monitoring. We describe the state of the art of four technologies (computer vision, acoustic monitoring, radar, and molecular methods), and assess their advantages, current limitations, and future potential. We discuss how these technologies can adhere to modern standards of data curation and transparency, their implications for citizen science, and their potential for integration among different monitoring programmes and technologies. We argue that they provide unprecedented possibilities for insect ecology and monitoring, but it will be important to foster international standards via collaboration.

Comparison of traditional and DNA metabarcoding samples for monitoring tropical soil arthropods (Formicidae, Collembola and Isoptera)

The soil fauna of the tropics remains one of the least known components of the biosphere. Long-term monitoring of this fauna is hampered by the lack of taxonomic expertise and funding. These obstacles may potentially be lifted with DNA metabarcoding. To validate this approach, we studied the ants, springtails and termites of 100 paired soil samples from Barro Colorado Island, Panama. The fauna was extracted with Berlese-Tullgren funnels and then either sorted with traditional taxonomy and known, individual DNA barcodes ("traditional samples") or processed with metabarcoding ("metabarcoding samples"). We detected 49 ant, 37 springtail and 34 termite species with 3.46 million reads of the COI gene, at a mean sequence length of 233 bp. Traditional identification yielded 80, 111 and 15 species of ants, springtails and termites, respectively; 98%, 37% and 100% of these species had a Barcode Index Number (BIN) allowing for direct comparison with metabarcoding. Ants were best surveyed through traditional methods, termites were better detected by metabarcoding, and springtails were equally well detected by both techniques. Species richness was underestimated, and faunal composition was different in metabarcoding samples, mostly because 37% of ant species were not detected. The prevalence of species in metabarcoding samples increased with their abundance in traditional samples, and seasonal shifts in species prevalence and faunal composition were similar between traditional and metabarcoding samples. Probable false positive and negative species records were reasonably low (13-18% of common species). We conclude that metabarcoding of samples extracted with Berlese-Tullgren funnels appear suitable for the long-term monitoring of termites and springtails in tropical rainforests. For ants, metabarcoding schemes should be complemented by additional samples of alates from Malaise or light traps.

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.

DNA Metabarcoding Enables High-Throughput Detection of Spotted Wing Drosophila (Drosophila suzukii) Within Unsorted Trap Catches

The spotted wing drosophila (Drosophila suzukii, Matsumara) is a rapidly spreading global pest of soft and stone fruit production. Due to the similarity of many of its life stages to other cosmopolitan drosophilids, surveillance for this pest is currently bottlenecked by the laborious sorting and morphological identification of large mixed trap catches. DNA metabarcoding presents an alternative high-throughput sequencing (HTS) approach for multi-species identification, which may lend itself ideally to rapid and scalable diagnostics of D. suzukii within unsorted trap samples. In this study, we compared the qualitative (identification accuracy) and quantitative (bias toward each species) performance of four metabarcoding primer pairs on D. suzukii and its close relatives. We then determined the sensitivity of a non-destructive metabarcoding assay (i.e., which retains intact specimens) by spiking whole specimens of target species into mock communities of increasing specimen number, as well as 29 field-sampled communities from a cherry and a stone fruit orchard. Metabarcoding successfully detected D. suzukii and its close relatives Drosophila subpulchrella and Drosophila biarmipes in the spiked communities with an accuracy of 96, 100, and 100% respectively, and identified a further 57 non-target arthropods collected as bycatch by D. suzukii surveillance methods in a field scenario. While the non-destructive DNA extraction retained intact voucher specimens, dropouts of single species and entire technical replicates suggests that these protocols behave more similarly to environmental DNA than homogenized tissue metabarcoding and may require increased technical replication to reliably detect low-abundance taxa. Adoption of high-throughput metabarcoding assays for screening bulk trap samples could enable a substantial increase in the geographic scale and intensity of D. suzukii surveillance, and thus likelihood of detecting a new introduction. Trap designs and surveillance protocols will, however, need to be optimized to adequately preserve specimen DNA for molecular identification.

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.

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.

Coming of age for COI metabarcoding of whole organism community DNA: Towards bioinformatic harmonisation

Metabarcoding of DNA extracted from community samples of whole organisms (whole organism community DNA, wocDNA) is increasingly being applied to terrestrial, marine and freshwater metazoan communities to provide rapid, accurate and high resolution data for novel molecular ecology research. The growth of this field has been accompanied by considerable development that builds on microbial metabarcoding methods to develop appropriate and efficient sampling and laboratory protocols for whole organism metazoan communities. However, considerably less attention has focused on ensuring bioinformatic methods are adapted and applied comprehensively in wocDNA metabarcoding. In this study we examined over 600 papers and identified 111 studies that performed COI metabarcoding of wocDNA. We then systematically reviewed the bioinformatic methods employed by these papers to identify the state‐of‐the‐art. Our results show that the increasing use of wocDNA COI metabarcoding for metazoan diversity is characterised by a clear absence of bioinformatic harmonisation, and the temporal trends show little change in this situation. The reviewed literature showed (i) high heterogeneity across pipelines, tasks and tools used, (ii) limited or no adaptation of bioinformatic procedures to the nature of the COI fragment, and (iii) a worrying underreporting of tasks, software and parameters. Based upon these findings we propose a set of recommendations that we think the metabarcoding community should consider to ensure that bioinformatic methods are appropriate, comprehensive and comparable. We believe that adhering to these recommendations will improve the long‐term integrative potential of wocDNA COI metabarcoding for biodiversity science.

Wet grinding of invertebrate bulk samples – a scalable and cost-efficient protocol for metabarcoding and metagenomics

Most metabarcoding protocols for invertebrate bulk samples start with sample homogenisation, followed by DNA extraction, amplification of a specific marker region, and sequencing. Many of the above-mentioned laboratory steps have been verified thoroughly and best practice strategies exist, yet, no clear recommendation for the basis of almost all metabarcoding studies exists: the homogenisation of samples itself. Two different categories of devices are typically used for homogenisation: bead mills or blenders. Both have upsides and downsides. Bead mills rely on single-use plastics and therefore produce a lot of waste and are expensive. In addition to that, processing times can go up to 30 minutes making them unsuitable for large-scale studies. Blenders can handle larger sample volumes in a shorter time, and be cleaned – yet suffer from an increased risk of cross-contamination. We aimed to develop a fast, robust, cheap, and reliable sample homogenisation protocol that overcomes limitations of both approaches, i.e. does not produce difficult to discard waste and avoid single-use plastics while reducing overall costs. We tested the performance of the new protocol using six size-sorted Malaise trap samples and six unsorted stream macroinvertebrate kick-net samples. We used 14 replicates per sample and included many negative controls at different steps of the protocol to quantify the impacts of i) insufficient homogenisation and ii) cross-contamination. Our results show that 3-min homogenisation is sufficient to recover about 80% of OTUs per sample in each replicate and that a non-hazardous DIY cleaning solution provides an effective and efficient way of cleaning. The improvements of the protocol in terms of speed, ease of handling, an overall reduction of costs as well as the documented reliability and robustness make it an important candidate for sample homogenisation after sampling in particular for large-scale and regulatory metabarcoding but also metagenomics biodiversity assessments and monitoring.

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.

All together now: Limitations and recommendations for the simultaneous analysis of all eukaryotic soil sequences

The soil environment contains a large, but historically underexplored, reservoir of biodiversity. Sequencing prokaryotic marker genes has become commonplace for the discovery and characterization of soil bacteria and archaea. Increasingly, this approach is also applied to eukaryotic marker genes to characterize the diversity and distribution of soil eukaryotes. However, understanding the properties and limitations of eukaryotic marker sequences is essential for correctly analysing, interpreting, and synthesizing the resulting data. Here, we illustrate several biases from sequencing data that affect measurements of biodiversity that arise from variation in morphology, taxonomy and phylogeny between organisms, as well as from sampling designs. We recommend analytical approaches to overcome these limitations, and outline how the benchmarking and standardization of sequencing protocols may improve the comparability of the data.

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.

Meta-analysis shows that environmental DNA outperforms traditional surveys, but warrants better reporting standards

Decades of environmental DNA (eDNA) method application, spanning a wide variety of taxa and habitats, has advanced our understanding of eDNA and underlined its value as a tool for conservation practitioners. The general consensus is that eDNA methods are more accurate and cost-effective than traditional survey methods. However, they are formally approved for just a few species globally (e.g., Bighead Carp, Silver Carp, Great Crested Newt). We conducted a meta-analysis of studies that directly compare eDNA with traditional surveys to evaluate the assertion that eDNA methods are consistently "better."Environmental DNA publications for multiple species or single macro-organism detection were identified using the Web of Science, by searching "eDNA" and "environmental DNA" across papers published between 1970 and 2020. The methods used, focal taxa, habitats surveyed, and quantitative and categorical results were collated and analyzed to determine whether and under what circumstances eDNA outperforms traditional surveys.Results show that eDNA methods are cheaper, more sensitive, and detect more species than traditional methods. This is, however, taxa-dependent, with amphibians having the highest potential for detection by eDNA survey. Perhaps most strikingly, of the 535 papers reviewed just 49 quantified the probability of detection for both eDNA and traditional survey methods and studies were three times more likely to give qualitative statements of performance. Synthesis and applications: The results of this meta-analysis demonstrate that where there is a direct comparison, eDNA surveys of macro-organisms are more accurate and efficient than traditional surveys. This conclusion, however, is based on just a fraction of available eDNA papers as most do not offer this granularity. We recommend that conclusions are substantiated with comparable and quantitative data. Where a direct comparison has not been made, we caution against the use of qualitative statements about relative performance. This consistency and rigor will simplify how the eDNA research community tracks methods-based advances and will also provide greater clarity for conservation practitioners. To this end suggest reporting standards for eDNA studies.

Connecting high-throughput biodiversity inventories: Opportunities for a site-based genomic framework for global integration and synthesis

High-throughput sequencing (HTS) is increasingly being used for the characterization and monitoring of biodiversity. If applied in a structured way, across broad geographical scales, it offers the potential for a much deeper understanding of global biodiversity through the integration of massive quantities of molecular inventory data generated independently at local, regional and global scales. The universality, reliability and efficiency of HTS data can potentially facilitate the seamless linking of data among species assemblages from different sites, at different hierarchical levels of diversity, for any taxonomic group and regardless of prior taxonomic knowledge. However, collective international efforts are required to optimally exploit the potential of site-based HTS data for global integration and synthesis, efforts that at present are limited to the microbial domain. To contribute to the development of an analogous strategy for the nonmicrobial terrestrial domain, an international symposium entitled "Next Generation Biodiversity Monitoring" was held in November 2019 in Nicosia (Cyprus). The symposium brought together evolutionary geneticists, ecologists and biodiversity scientists involved in diverse regional and global initiatives using HTS as a core tool for biodiversity assessment. In this review, we summarize the consensus that emerged from the 3-day symposium. We converged on the opinion that an effective terrestrial Genomic Observatories network for global biodiversity integration and synthesis should be spatially led and strategically united under the umbrella of the metabarcoding approach. Subsequently, we outline an HTS-based strategy to collectively build an integrative framework for site-based biodiversity data generation.

Standards and Best Practices for Monitoring and Benchmarking Insects

Benchmark studies of insect populations are increasingly relevant and needed amid accelerating concern about insect trends in the Anthropocene. The growing recognition that insect populations may be in decline has given rise to a renewed call for insect population monitoring by scientists, and a desire from the broader public to participate in insect surveys. However, due to the immense diversity of insects and a vast assortment of data collection methods, there is a general lack of standardization in insect monitoring methods, such that a sudden and unplanned expansion of data collection may fail to meet its ecological potential or conservation needs without a coordinated focus on standards and best practices. To begin to address this problem, we provide simple guidelines for maximizing return on proven inventory methods that will provide insect benchmarking data suitable for a variety of ecological responses, including occurrence and distribution, phenology, abundance and biomass, and diversity and species composition. To track these responses, we present seven primary insect sampling methods—malaise trapping, light trapping, pan trapping, pitfall trappings, beating sheets, acoustic monitoring, and active visual surveys—and recommend standards while highlighting examples of model programs. For each method, we discuss key topics such as recommended spatial and temporal scales of sampling, important metadata to track, and degree of replication needed to produce rigorous estimates of ecological responses. We additionally suggest protocols for scalable insect monitoring, from backyards to national parks. Overall, we aim to compile a resource that can be used by diverse individuals and organizations seeking to initiate or improve insect monitoring programs in this era of rapid change.

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.

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.

In silico and empirical evaluation of twelve metabarcoding primer sets for insectivorous diet analyses

During the most recent decade, environmental DNA metabarcoding approaches have been both developed and improved to minimize the biological and technical biases in these protocols. However, challenges remain, notably those relating to primer design. In the current study, we comprehensively assessed the performance of ten COI and two 16S primer pairs for eDNA metabarcoding, including novel and previously published primers. We used a combined approach of in silico, in vivo-mock community (33 arthropod taxa from 16 orders), and guano-based analyses to identify primer sets that would maximize arthropod detection and taxonomic identification, successfully identify the predator (bat) species, and minimize the time and financial costs of the experiment. We focused on two insectivorous bat species that live together in mixed colonies: the greater horseshoe bat (Rhinolophus ferrumequinum) and Geoffroy's bat (Myotis emarginatus). We found that primer degeneracy is the main factor that influences arthropod detection in silico and mock community analyses, while amplicon length is critical for the detection of arthropods from degraded DNA samples. Our guano-based results highlight the importance of detecting and identifying both predator and prey, as guano samples can be contaminated by other insectivorous species. Moreover, we demonstrate that amplifying bat DNA does not reduce the primers' capacity to detect arthropods. We therefore recommend the simultaneous identification of predator and prey. Finally, our results suggest that up to one-third of prey occurrences may be unreliable and are probably not of primary interest in diet studies, which may decrease the relevance of combining several primer sets instead of using a single efficient one. In conclusion, this study provides a pragmatic framework for eDNA primer selection with respect to scientific and methodological constraints.

Comparing diversity levels in environmental samples: DNA sequence capture and metabarcoding approaches using 18S and COI genes

Environmental DNA studies targeting multiple taxa using metabarcoding provide remarkable insights into levels of species diversity in any habitat. The main drawbacks are the presence of primer bias and difficulty in identifying rare species. We tested a DNA sequence-capture method in parallel with the metabarcoding approach to reveal possible advantages of one method over the other. Both approaches were performed using the same eDNA samples and the same 18S and COI regions, followed by high throughput sequencing. Metabarcoded eDNA libraries were PCR amplified with one primer pair from 18S and COI genes. DNA sequence-capture libraries were enriched with 3,639 baits targeting the same gene regions. We tested amplicon sequence variants (ASVs) and operational taxonomic units (OTUs) in silico approaches for both markers and methods, using for this purpose the metabarcoding data set. ASVs methods uncovered more species for the COI gene, whereas the opposite occurred for the 18S gene, suggesting that clustering reads into OTUs could bias diversity richness especially using 18S with relaxed thresholds. Additionally, metabarcoding and DNA sequence-capture recovered 80%-90% of the control sample species. DNA sequence-capture was 8x more expensive, nonetheless it identified 1.5x more species for COI and 13x more genera for 18S than metabarcoding. Both approaches offer reliable results, sharing ca. 40% species and 72% families and retrieve more taxa when nuclear and mitochondrial markers are combined. eDNA metabarcoding is quite well established and low-cost, whereas DNA-sequence capture for biodiversity assessment is still in its infancy, is more time-consuming but provides more taxonomic assignments.

Completing Linnaeus's inventory of the Swedish insect fauna: Only 5,000 species left?

Despite more than 250 years of taxonomic research, we still have only a vague idea about the true size and composition of the faunas and floras of the planet. Many biodiversity inventories provide limited insight because they focus on a small taxonomic subsample or a tiny geographic area. Here, we report on the size and composition of the Swedish insect fauna, thought to represent roughly half of the diversity of multicellular life in one of the largest European countries. Our results are based on more than a decade of data from the Swedish Taxonomy Initiative and its massive inventory of the country's insect fauna, the Swedish Malaise Trap Project The fauna is considered one of the best known in the world, but the initiative has nevertheless revealed a surprising amount of hidden diversity: more than 3,000 new species (301 new to science) have been documented so far. Here, we use three independent methods to analyze the true size and composition of the fauna at the family or subfamily level: (1) assessments by experts who have been working on the most poorly known groups in the fauna; (2) estimates based on the proportion of new species discovered in the Malaise trap inventory; and (3) extrapolations based on species abundance and incidence data from the inventory. For the last method, we develop a new estimator, the combined non-parametric estimator, which we show is less sensitive to poor coverage of the species pool than other popular estimators. The three methods converge on similar estimates of the size and composition of the fauna, suggesting that it comprises around 33,000 species. Of those, 8,600 (26%) were unknown at the start of the inventory and 5,000 (15%) still await discovery. We analyze the taxonomic and ecological composition of the estimated fauna, and show that most of the new species belong to Hymenoptera and Diptera groups that are decomposers or parasitoids. Thus, current knowledge of the Swedish insect fauna is strongly biased taxonomically and ecologically, and we show that similar but even stronger biases have distorted our understanding of the fauna in the past. We analyze latitudinal gradients in the size and composition of known European insect faunas and show that several of the patterns contradict the Swedish data, presumably due to similar knowledge biases. Addressing these biases is critical in understanding insect biomes and the ecosystem services they provide. Our results emphasize the need to broaden the taxonomic scope of current insect monitoring efforts, a task that is all the more urgent as recent studies indicate a possible worldwide decline in insect faunas.

High-throughput sequencing for community analysis: the promise of DNA barcoding to uncover diversity, relatedness, abundances and interactions in spider communities

Large-scale studies on community ecology are highly desirable but often difficult to accomplish due to the considerable investment of time, labor and, money required to characterize richness, abundance, relatedness, and interactions. Nonetheless, such large-scale perspectives are necessary for understanding the composition, dynamics, and resilience of biological communities. Small invertebrates play a central role in ecosystems, occupying critical positions in the food web and performing a broad variety of ecological functions. However, it has been particularly difficult to adequately characterize communities of these animals because of their exceptionally high diversity and abundance. Spiders in particular fulfill key roles as both predator and prey in terrestrial food webs and are hence an important focus of ecological studies. In recent years, large-scale community analyses have benefitted tremendously from advances in DNA barcoding technology. High-throughput sequencing (HTS), particularly DNA metabarcoding, enables community-wide analyses of diversity and interactions at unprecedented scales and at a fraction of the cost that was previously possible. Here, we review the current state of the application of these technologies to the analysis of spider communities. We discuss amplicon-based DNA barcoding and metabarcoding for the analysis of community diversity and molecular gut content analysis for assessing predator-prey relationships. We also highlight applications of the third generation sequencing technology for long read and portable DNA barcoding. We then address the development of theoretical frameworks for community-level studies, and finally highlight critical gaps and future directions for DNA analysis of spider communities.

The Swedish Malaise Trap Project: A 15 Year Retrospective on a Countrywide Insect Inventory

The Swedish Malaise Trap Project (SMTP) is one of the most ambitious insect inventories ever attempted. The project was designed to target poorly known insect groups across a diverse range of habitats in Sweden. The field campaign involved the deployment of 73 Malaise traps at 55 localities across the country for three years (2003-2006). Over the past 15 years, the collected material has been hand sorted by trained technicians into over 300 taxonomic fractions suitable for expert attention. The resulting collection is a tremendous asset for entomologists around the world, especially as we now face a desperate need for baseline data to evaluate phenomena like insect decline and climate change. Here, we describe the history, organisation, methodology and logistics of the SMTP, focusing on the rationale for the decisions taken and the lessons learned along the way. The SMTP represents one of the early instances of community science applied to large-scale inventory work, with a heavy reliance on volunteers in both the field and the laboratory. We give estimates of both staff effort and volunteer effort involved. The project has been funded by the Swedish Taxonomy Initiative; in total, the inventory has cost less than 30 million SEK (approximately 3.1 million USD). Based on a subset of the samples, we characterise the size and taxonomic composition of the SMTP material. Several different extrapolation methods suggest that the material comprises around 20 million specimens in total. The material is dominated by Diptera (75% of the specimens) and Hymenoptera (15% of specimens). Amongst the Diptera, the dominant groups are Chironomidae (37% of specimens), Sciaridae (15%), Phoridae (13%), Cecidomyiidae (9.5%) and Mycetophilidae (9.4%). Within Hymenoptera, the major groups are Ichneumonidae (44% of specimens), Diaprioidea (19%), Braconidae (9.6%), Platygastroidea (8.5%) and Chalcidoidea (7.9%). The taxonomic composition varies with latitude and season. Several Diptera and Hymenoptera groups are more common in non-summer samples (collected from September to April) and in the North, while others show the opposite pattern. About 1% of the total material has been processed and identified by experts so far. This material represents over 4,000 species. One third of these had not been recorded from Sweden before and almost 700 of them are new to science. These results reveal the large amounts of taxonomic work still needed on Palaearctic insect faunas. Based on the SMTP experiences, we discuss aspects of planning and conducting future large-scale insect inventory projects using mainly traditional approaches in relation to more recent approaches that rely on molecular techniques.

Establishing arthropod community composition using metabarcoding: Surprising inconsistencies between soil samples and preservative ethanol and homogenate from Malaise trap catches

DNA metabarcoding allows the analysis of insect communities faster and more efficiently than ever before. However, metabarcoding can be conducted through several approaches, and the consistency of results across methods has rarely been studied. We compare the results obtained by DNA metabarcoding of the same communities using two different markers – COI and 16S – and three different sampling methods: (a) homogenized Malaise trap samples (homogenate), (b) preservative ethanol from the same samples, and (c) soil samples. Our results indicate that COI and 16S offer partly complementary information on Malaise trap samples, with each marker detecting a significant number of species not detected by the other. Different sampling methods offer highly divergent estimates of community composition. The community recovered from preservative ethanol of Malaise trap samples is significantly different from that recovered from homogenate. Small and weakly sclerotized insects tend to be overrepresented in ethanol while strong and large taxa are overrepresented in homogenate. For soil samples, highly degenerate COI primers pick up large amounts of nontarget DNA and only 16S provides adequate analyses of insect diversity. However, even with 16S, very little overlap in molecular operational taxonomic unit (MOTU) content was found between the trap and the soil samples. Our results demonstrate that none of the tested sampling approaches is satisfactory on its own. For instance, DNA extraction from preservative ethanol is not a valid replacement for destructive bulk extraction but a complement. In future metabarcoding studies, both should ideally be used together to achieve comprehensive representation of the target community.