The effect of DNA degradation bias in passive sampling devices on metabarcoding studies of arthropod communities and their associated microbiota

Effectiveness of eDNA for monitoring riverine macroinvertebrates

Environmental DNA (eDNA) is a technique increasingly used for monitoring organisms in the natural environment including riverine macroinvertebrates. However, the effectiveness of eDNA for monitoring riverine macroinvertebrates compared with the more traditional method of sampling the organisms directly and identifying them via morphological analysis, has not been well established. Furthermore, the ability of the various gene markers and PCR primer sets to detect the full range of riverine invertebrate taxa has not been quantified. Here we conducted a meta-analysis of the available literature, to assess the effectiveness of eDNA sampling for detecting riverine macroinvertebrates compared with sampling for the organisms directly and applying morphological analysis. We found, on average, eDNA sampling, irrespective of the gene marker used, detected fewer riverine invertebrates than morphological sampling. The most effective PCR primer set for identifying taxa was mlCOIintF/jgHCO2198, (mlCOIintF- forward primer, jgHCO2198, - reverse primer). Regardless of the gene marker or primer sets used, however, many taxa were not detected by eDNA metabarcoding that were detected by sampling directly for these invertebrates, including over 100 members of Arthropoda. eDNA sampling failed to detect any species belonging to Nematoda, Platyhelminthes, Cnidaria or Nematomorpha and these markers applied for eDNA sampling in terrestrial systems also do not detect members of Nematoda. In addition to these issues, uncertainties relating to false positives from upstream DNA sources, the stability of DNA from different species, differences in the propensity for DNA release into the environment for different organisms, and lack of available sequence information for numerous taxa illustrates the use of eDNA is not yet applicable as a robust stand-alone method for the monitoring of riverine invertebrates. As a primary consideration, further methodological developments are needed to ensure eDNA captures some of the key freshwater taxa, notably taxa belonging to the phyla Arthropoda, Nematoda, Platyhelminthes, Cnidaria and Nematomorpha.

eDNA metabarcoding of avocado flowers: 'Hass' it got potential to survey arthropods in food production systems?

In the face of global biodiversity declines, surveys of beneficial and antagonistic arthropod diversity as well as the ecological services that they provide are increasingly important in both natural and agro-ecosystems. Conventional survey methods used to monitor these communities often require extensive taxonomic expertise and are time-intensive, potentially limiting their application in industries such as agriculture, where arthropods often play a critical role in productivity (e.g. pollinators, pests and predators). Environmental DNA (eDNA) metabarcoding of a novel substrate, crop flowers, may offer an accurate and high throughput alternative to aid in the detection of these managed and unmanaged taxa. Here, we compared the arthropod communities detected with eDNA metabarcoding of flowers, from an agricultural species (Persea americana-'Hass' avocado), with two conventional survey techniques: digital video recording (DVR) devices and pan traps. In total, 80 eDNA flower samples, 96 h of DVRs and 48 pan trap samples were collected. Across the three methods, 49 arthropod families were identified, of which 12 were unique to the eDNA dataset. Environmental DNA metabarcoding from flowers revealed potential arthropod pollinators, as well as plant pests and parasites. Alpha diversity levels did not differ across the three survey methods although taxonomic composition varied significantly, with only 12% of arthropod families found to be common across all three methods. eDNA metabarcoding of flowers has the potential to revolutionize the way arthropod communities are monitored in natural and agro-ecosystems, potentially detecting the response of pollinators and pests to climate change, diseases, habitat loss and other disturbances.

Haplotype-level metabarcoding of freshwater macroinvertebrate species: A prospective tool for population genetic analysis

Metabarcoding is a molecular-based tool capable of large quantity high-throughput species identification from bulk samples that is a faster and more cost-effective alternative to conventional DNA-sequencing approaches. Still, further exploration and assessment of the laboratory and bioinformatics strategies are required to unlock the potential of metabarcoding-based inference of haplotype information. In this study, we assessed the inference of freshwater macroinvertebrate haplotypes from metabarcoding data in a mock sample. We also examined the influence of DNA template concentration and PCR cycle on detecting true and spurious haplotypes. We tested this strategy on a mock sample containing twenty individuals from four species with known haplotypes based on the 658-bp Folmer region of the mitochondrial cytochrome c oxidase gene. We recovered fourteen zero-radius operational taxonomic units (zOTUs) of 421-bp length, with twelve zOTUs having a 100% match with the Sanger haplotype sequences. High-quality reads relatively increased with increasing PCR cycles, and the relative abundance of each zOTU was consistent for each cycle. This suggests that increasing the PCR cycles from 24 to 64 did not affect the relative abundance of each zOTU. As metabarcoding becomes more established and laboratory protocols and bioinformatic pipelines are continuously being developed, our study demonstrated the method's ability to infer intraspecific variability while highlighting the challenges that must be addressed before its eventual application for population genetic studies.

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.

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.

Accurate quantification of bacterial abundance in metagenomic DNAs accounting for variable DNA integrity levels

The quantification of the total microbial content in metagenomic samples is critical for investigating the interplay between the microbiome and its host, as well as for assessing the accuracy and precision of the relative microbial composition which can be strongly biased in low microbial biomass samples. In the present study, we demonstrate that digital droplet PCR (ddPCR) can provide accurate quantification of the total copy number of the 16S rRNA gene, the gene usually exploited for assessing total bacterial abundance in metagenomic DNA samples. Notably, using DNA templates with different integrity levels, as measured by the DNA integrity number (DIN), we demonstrated that 16S rRNA copy number quantification is strongly affected by DNA quality and determined a precise correlation between quantification underestimation and DNA degradation levels. Therefore, we propose an input DNA mass correction, according to the observed DIN value, which could prevent inaccurate quantification of 16S copy number in degraded metagenomic DNAs. Our results highlight that a preliminary evaluation of the metagenomic DNA integrity should be considered before performing metagenomic analyses of different samples, both for the assessment of the reliability of observed differential abundances in different conditions and to obtain significant functional insights.

Soil Sample Preservation Strategy Affects the Microbial Community Structure

Sample preservation is a critical procedure in any research that relies on molecular tools and is conducted in remote areas. Sample preservation options include low and room temperature storage, which require freezing equipment and specific buffering solutions, respectively. The aim of the present study was to investigate whether DNA/RNA Shield 1x from Zymo Research and DESS (Dimethyl sulfoxide, Ethylenediamine tetraacetic acid, Saturated Salt) solution performed similarly to snap freezing in liquid nitrogen. Soil samples were stored for 1 month in each of the buffers and without any solution at a range of temperatures: –20, +4, and +23°C. All treatments were compared to the “optimal treatment”, namely, snap freezing in liquid nitrogen. The quality and quantity of DNA were analyzed, and the microbial community structure was investigated in all samples. The results obtained indicated that the quantity and integrity of DNA was preserved well in all samples; however, the taxonomic distribution was skewed in samples stored without any solution at ambient temperatures, particularly when analyses were performed at lower taxonomic levels. Although both solutions performed equally well, sequencing output and OTU numbers in DESS-treated samples were closer to those snap frozen with liquid nitrogen. Furthermore, DNA/RNA Shield-stored samples performed better for the preservation of rare taxa.

Colonization of Beauveria bassiana 08F04 in root-zone soil and its biocontrol of cereal cyst nematode (Heterodera filipjevi)

Cereal cyst nematodes cause serious yield losses of wheat in Hunaghuai winter wheat growing region in China. Beauveria bassiana 08F04 isolated from the surface of cysts is a promising biological control agent for cereal cyst nematodes. As the colonization capacity is a crucial criteria to assess biocontrol effectiveness for a microbial agent candidate, we aimed to label B. bassiana 08F04 for efficient monitoring of colonization in the soil. The binary pCAM-gfp plasmid containing sgfp and hph was integrated into B. bassiana 08F04 using the Agrobacterium tumefaciens-mediated transformation. The transformation caused a significant change in mycelial and conidial yields, and in extracellular chitinase activity in some transformants. The cultural filtrates of some transformants also decreased acetylcholinesterase activity and the survival of Heterodera filipjevi second-stage juveniles relative to the wild-type strain. One transformant (G10) had a growth rate and biocontrol efficacy similar to the wild-type strain, so it was used for a pilot study of B. bassiana colonization conducted over 13 weeks. Real-time PCR results and CFU counts revealed that the population of G10 increased quickly over the first 3 weeks, then decreased slowly over the following 4 weeks before stabilizing. In addition, the application of wild-type B. bassiana 08F04 and transformant G10 significantly reduced the number of H. filipjevi females in roots by 64.4% and 60.2%, respectively. The results of this study have practical applications for ecological, biological and functional studies of B. bassiana 08F04 and for bionematicide registration.

Validation of COI metabarcoding primers for terrestrial arthropods

Metabarcoding can rapidly determine the species composition of bulk samples and thus aids biodiversity and ecosystem assessment. However, it is essential to use primer sets that minimize amplification bias among taxa to maximize species recovery. Despite this fact, the performance of primer sets employed for metabarcoding terrestrial arthropods has not been sufficiently evaluated. This study tests the performance of 36 primer sets on a mock community containing 374 insect species. Amplification success was assessed with gradient PCRs and the 21 most promising primer sets selected for metabarcoding. These 21 primer sets were also tested by metabarcoding a Malaise trap sample. We identified eight primer sets, mainly those including inosine and/or high degeneracy, that recovered more than 95% of the species in the mock community. Results from the Malaise trap sample were congruent with the mock community, but primer sets generating short amplicons produced potential false positives. Taxon recovery from both mock community and Malaise trap sample metabarcoding were used to select four primer sets for additional evaluation at different annealing temperatures (40-60 °C) using the mock community. The effect of temperature varied by primer pair but overall it only had a minor effect on taxon recovery. This study reveals the weak performance of some primer sets employed in past studies. It also demonstrates that certain primer sets can recover most taxa in a diverse species assemblage. Thus, based our experimental set up, there is no need to employ several primer sets targeting the same gene region. We identify several suitable primer sets for arthropod metabarcoding, and specifically recommend BF3 + BR2, as it is not affected by primer slippage and provides maximal taxonomic resolution. The fwhF2 + fwhR2n primer set amplifies a shorter fragment and is therefore ideal when targeting degraded DNA (e.g., from gut contents).

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.

Prospects and challenges of implementing DNA metabarcoding for high-throughput insect surveillance

Trap-based surveillance strategies are widely used for monitoring of invasive insect species, aiming to detect newly arrived exotic taxa as well as track the population levels of established or endemic pests. Where these surveillance traps have low specificity and capture non-target endemic species in excess of the target pests, the need for extensive specimen sorting and identification creates a major diagnostic bottleneck. While the recent development of standardized molecular diagnostics has partly alleviated this requirement, the single specimen per reaction nature of these methods does not readily scale to the sheer number of insects trapped in surveillance programmes. Consequently, target lists are often restricted to a few high-priority pests, allowing unanticipated species to avoid detection and potentially establish populations. DNA metabarcoding has recently emerged as a method for conducting simultaneous, multi-species identification of complex mixed communities and may lend itself ideally to rapid diagnostics of bulk insect trap samples. Moreover, the high-throughput nature of recent sequencing platforms could enable the multiplexing of hundreds of diverse trap samples on a single flow cell, thereby providing the means to dramatically scale up insect surveillance in terms of both the quantity of traps that can be processed concurrently and number of pest species that can be targeted. In this review of the metabarcoding literature, we explore how DNA metabarcoding could be tailored to the detection of invasive insects in a surveillance context and highlight the unique technical and regulatory challenges that must be considered when implementing high-throughput sequencing technologies into sensitive diagnostic applications.

Have the cake and eat it: Optimizing nondestructive DNA metabarcoding of macroinvertebrate samples for freshwater biomonitoring

DNA metabarcoding can contribute to improving cost-effectiveness and accuracy of biological assessments of aquatic ecosystems, but significant optimization and standardization efforts are still required to mainstream its application into biomonitoring programmes. In assessments based on freshwater macroinvertebrates, a key challenge is that DNA is often extracted from cleaned, sorted and homogenized bulk samples, which is time-consuming and may be incompatible with sample preservation requirements of regulatory agencies. Here, we optimize and evaluate metabarcoding procedures based on DNA recovered from 96% ethanol used to preserve field samples and thus including potential PCR inhibitors and nontarget organisms. We sampled macroinvertebrates at five sites and subsampled the preservative ethanol at 1 to 14 days thereafter. DNA was extracted using column-based enzymatic (TISSUE) or mechanic (SOIL) protocols, or with a new magnetic-based enzymatic protocol (BEAD), and a 313-bp COI fragment was amplified. Metabarcoding detected at least 200 macroinvertebrate taxa, including most taxa detected through morphology and for which there was a reference barcode. Better results were obtained with BEAD than SOIL or TISSUE, and with subsamples taken 7-14 than 1-7 days after sampling, in terms of DNA concentration and integrity, taxa diversity and matching between metabarcoding and morphology. Most variation in community composition was explained by differences among sites, with small but significant contributions of subsampling day and extraction method, and negligible contributions of extraction and PCR replication. Our methods enhance reliability of preservative ethanol as a potential source of DNA for macroinvertebrate metabarcoding, with a strong potential application in freshwater biomonitoring.

Nanopore sequencing of long ribosomal DNA amplicons enables portable and simple biodiversity assessments with high phylogenetic resolution across broad taxonomic scale

BACKGROUND

In light of the current biodiversity crisis, DNA barcoding is developing into an essential tool to quantify state shifts in global ecosystems. Current barcoding protocols often rely on short amplicon sequences, which yield accurate identification of biological entities in a community but provide limited phylogenetic resolution across broad taxonomic scales. However, the phylogenetic structure of communities is an essential component of biodiversity. Consequently, a barcoding approach is required that unites robust taxonomic assignment power and high phylogenetic utility. A possible solution is offered by sequencing long ribosomal DNA (rDNA) amplicons on the MinION platform (Oxford Nanopore Technologies).

FINDINGS

Using a dataset of various animal and plant species, with a focus on arthropods, we assemble a pipeline for long rDNA barcode analysis and introduce a new software (MiniBar) to demultiplex dual indexed Nanopore reads. We find excellent phylogenetic and taxonomic resolution offered by long rDNA sequences across broad taxonomic scales. We highlight the simplicity of our approach by field barcoding with a miniaturized, mobile laboratory in a remote rainforest. We also test the utility of long rDNA amplicons for analysis of community diversity through metabarcoding and find that they recover highly skewed diversity estimates.

CONCLUSIONS

Sequencing dual indexed, long rDNA amplicons on the MinION platform is a straightforward, cost-effective, portable, and universal approach for eukaryote DNA barcoding. Although bulk community analyses using long-amplicon approaches may introduce biases, the long rDNA amplicons approach signifies a powerful tool for enabling the accurate recovery of taxonomic and phylogenetic diversity across biological communities.

Toward accurate species-level metabarcoding of arthropod communities from the tropical forest canopy

Metabarcoding of arthropod communities can be used for assessing species diversity in tropical forests but the methodology requires validation for accurate and repeatable species occurrences in complex mixtures. This study investigates how the composition of ecological samples affects the accuracy of species recovery.Starting with field-collected bulk samples from the tropical canopy, the recovery of specimens was tested for subsets of different body sizes and major taxa, by assembling these subsets into increasingly complex composite pools. After metabarcoding, we track whether richness, diversity, and most importantly composition of any size class or taxonomic subset are affected by the presence of other subsets in the mixture.Operational taxonomic units (OTUs) greatly exceeded the number of morphospecies in most taxa, even under very stringent sequencing read filtering. There was no significant effect on the recovered OTU richness of small and medium-sized arthropods when metabarcoded alongside larger arthropods, despite substantial biomass differences in the mixture. The recovery of taxonomic subsets was not generally influenced by the presence of other taxa, although with some exceptions likely due to primer mismatches. Considerable compositional variation within size and taxon-based subcommunities was evident resulting in high beta-diversity among samples from within a single tree canopy, but this beta-diversity was not affected by experimental manipulation.We conclude that OTU recovery in complex arthropod communities, with sufficient sequencing depth and within reasonable size ranges, is not skewed by variable biomass of the constituent species. This could remove the need for time-intensive manual sorting prior to metabarcoding. However, there remains a chance of taxonomic bias, which may be primer-dependent. There will never be a panacea primer; instead, metabarcoding studies should carefully consider whether the aim is broadscale turnover, in which case these biases may not be important, or species lists, in which case separate PCRs and sequencing might be necessary. OTU number inflation remains an issue in metabarcoding and requires bioinformatic development, particularly in read filtering and OTU clustering, and/or greater use of species-identifying sequences generated outside of bulk sequencing.