NGS metabarcoding proves successful for quantitative assessment of symbiont abundance: the case of feather mites on birds

Mosquito Identification From Bulk Samples Using DNA Metabarcoding: a Protocol to Support Mosquito-Borne Disease Surveillance in Canada

Approximately 80 species of mosquitoes (Diptera: Culicidae) have been documented in Canada. Exotic species such as Aedes albopictus (Skuse) (Diptera: Culicidae) are becoming established. Recently occurring endemic mosquito-borne diseases (MBD) in Canada including West-Nile virus (WNV) and Eastern Equine Encephalitis (EEE) are having significant public health impacts. Here we explore the use of DNA metabarcoding to identify mosquitoes from CDC light-trap collections from two locations in eastern Canada. Two primer pairs (BF2-BR2 and F230) were used to amplify regions of the cytochrome c oxidase subunit I (CO1) gene. High throughput sequencing was conducted using an Illumina MiSeq platform and GenBank-based species identification was applied using a QIIME 1.9 bioinformatics pipeline. From a site in southeastern Ontario, Canada, 26 CDC light trap collections of 72 to >300 individual mosquitoes were used to explore the capacity of DNA metabarcoding to identify and quantify captured mosquitoes. The DNA metabarcoding method identified 33 species overall while 24 species were identified by key. Using replicates from each trap, the dried biomass needed to identify the majority of species was determined to be 76 mg (equivalent to approximately 72 mosquitoes), and at least two replicates from the dried biomass would be needed to reliably detect the majority of species in collections of 144-215 mosquitoes and three replicates would be advised for collections with >215 mosquitoes. This study supports the use of DNA metabarcoding as a mosquito surveillance tool in Canada which can help identify the emergence of new mosquito-borne disease potential threats.

Culicidae-centric metabarcoding through targeted use of D2 ribosomal DNA primers

A practical limitation to many metabarcoding initiatives is that sampling methods tend to collect many non-target taxa, which become “amplicon noise” that can saturate Next Generation Sequencing results and lead to both financial and resource inefficiencies. An available molecular tool that can significantly decrease these non-target amplicons and decrease the need for pre-DNA-extraction sorting of bycatch is the design of PCR primers tailored to the taxa under investigation. We assessed whether the D2 extension segment of the 28S ribosomal operon can limit this shortcoming within the context of mosquito (Culicidae) monitoring. We designed PCR primers that are fully conserved across mosquitos and exclude from amplification most other taxa likely to be collected with current sampling apparatuses. We show that, given enough sequencing depth, D2 is an effective marker for the detection of mosquito sequences within mock genomic DNA pools. As few as 3,050 quality-filtered Illumina reads were able to recover all 17 species in a bulk pool containing as little as 0.2% of constituent DNA from single taxa. We also mixed these mosquito DNA pools with high concentrations of non-Culicidae bycatch DNA and show that the component mosquito species are generally still recoverable and faithful to their original relative frequencies. Finally, we show that there is little loss of fidelity in abundance parameters when pools from degraded DNA samples were sequenced using the D2 primers.

Understanding PCR Processes to Draw Meaningful Conclusions from Environmental DNA Studies

As environmental DNA (eDNA) studies have grown in popularity for use in ecological applications, it has become clear that their results differ in significant ways from those of traditional, non-PCR-based surveys. In general, eDNA studies that rely on amplicon sequencing may detect hundreds of species present in a sampled environment, but the resulting species composition can be idiosyncratic, reflecting species' true biomass abundances poorly or not at all. Here, we use a set of simulations to develop a mechanistic understanding of the processes leading to the kinds of results common in mixed-template PCR-based (metabarcoding) studies. In particular, we focus on the effects of PCR cycle number and primer amplification efficiency on the results of diversity metrics in sequencing studies. We then show that proportional indices of amplicon reads capture trends in taxon biomass with high accuracy, particularly where amplification efficiency is high (median correlation up to 0.97). Our results explain much of the observed behavior of PCR-based studies, and lead to recommendations for best practices in the field.

Unexpected bird-feather mite associations revealed by DNA metabarcoding uncovers a dynamic ecoevolutionary scenario

The high relevance of host-switching for the diversification of highly host-specific symbionts (i.e., those commonly inhabiting a single host species) demands a better understanding of host-switching dynamics at an ecological scale. Here, we used DNA metabarcoding to study feather mites on passerine birds in Spain, sequencing mtDNA (COI) for 25,540 individual mites (representing 64 species) from 1,130 birds (representing 71 species). Surprisingly, 1,228 (4.8%) mites from 84 (7.4%) birds were found on host species that were not the expected to be a host according to a recent bird-feather mite associations catalog. Unexpected associations were widespread across studied mite (40.6%) and bird (43.7%) species and showed smaller average infrapopulation sizes than typical associations. Unexpected mite species colonized hosts being distantly related to the set of their usual hosts, but with similar body size. The network of bird-mite associations was modular (i.e., some groups of bird and mite species tended to be more associated with each other than with the others), with 75.9% of the unexpected associations appearing within the module of the typical hosts of the mite species. Lastly, 68.4% of mite species found on unexpected hosts showed signatures of genetic differentiation, and we found evidence for reproduction or the potential for it in many of the unexpected associations. Results show host colonization as a common phenomenon even for these putatively highly host-specific symbionts. Thus, host-switching by feather mites, rather than a rare phenomenon, appears as a relatively frequent phenomenon shaped by ecological filters such as host morphology and is revealed as a fundamental component for a dynamic coevolutionary and codiversification scenario.

Biases in the metabarcoding of plant pathogens using rust fungi as a model system

Plant pathogens such as rust fungi (Pucciniales) are of global economic and ecological importance. This means there is a critical need to reliably and cost-effectively detect, identify, and monitor these fungi at large scales. We investigated and analyzed the causes of differences between next-generation sequencing (NGS) metabarcoding approaches and traditional DNA cloning in the detection and quantification of recognized species of rust fungi from environmental samples. We found significant differences between observed and expected numbers of shared rust fungal operational taxonomic units (OTUs) among different methods. However, there was no significant difference in relative abundance of OTUs that all methods were capable of detecting. Differences among the methods were mainly driven by the method's ability to detect specific OTUs, likely caused by mismatches with the NGS metabarcoding primers to some Puccinia species. Furthermore, detection ability did not seem to be influenced by differences in sequence lengths among methods, the most appropriate bioinformatic pipeline used for each method, or the ability to detect rare species. Our findings are important to future metabarcoding studies, because they highlight the main sources of difference among methods, and rule out several mechanisms that could drive these differences. Furthermore, strong congruity among three fundamentally different and independent methods demonstrates the promising potential of NGS metabarcoding for tracking important taxa such as rust fungi from within larger NGS metabarcoding communities. Our results support the use of NGS metabarcoding for the large-scale detection and quantification of rust fungi, but not for confirming the absence of species.

How quantitative is metabarcoding: A meta-analytical approach

Metabarcoding has been used in a range of ecological applications such as taxonomic assignment, dietary analysis and the analysis of environmental DNA. However, after a decade of use in these applications there is little consensus on the extent to which proportions of reads generated corresponds to the original proportions of species in a community. To quantify our current understanding, we conducted a structured review and meta‐analysis. The analysis suggests that a weak quantitative relationship may exist between the biomass and sequences produced (slope = 0.52 ± 0.34, p < 0.01), albeit with a large degree of uncertainty. None of the tested moderators, sequencing platform type, the number of species used in a trial or the source of DNA, were able to explain the variance. Our current understanding of the factors affecting the quantitative performance of metabarcoding is still limited: additional research is required before metabarcoding can be confidently utilized for quantitative applications. Until then, we advocate the inclusion of mock communities when metabarcoding as this facilitates direct assessment of the quantitative ability of any given study.

Enabling large-scale feather mite studies: an Illumina DNA metabarcoding pipeline

Feather mites are among the most common and diverse ectosymbionts of birds, yet basic questions such as the nature of their relationship remain largely unanswered. One reason for feather mites being understudied is that their morphological identification is often virtually impossible when using female or young individuals. Even for adult male specimens this task is tedious and requires advanced taxonomic expertise, thus hampering large-scale studies. In addition, molecular-based methods are challenging because the low DNA amounts usually obtained from these tiny mites do not reach the levels required for high-throughput sequencing. This work aims to overcome these issues by using a DNA metabarcoding approach to accurately identify and quantify the feather mite species present in a sample. DNA metabarcoding is a widely used molecular technique that takes advantage of high-throughput sequencing methodologies to assign the taxonomic identity to all the organisms present in a complex sample (i.e., a sample made up of multiple specimens that are hard or impossible to individualise). We present a high-throughput method for feather mite identification using a fragment of the COI gene as marker and Illumina Miseq technology. We tested this method by performing two experiments plus a field test over a total of 11,861 individual mites (5360 of which were also morphologically identified). In the first experiment, we tested the probability of detecting a single feather mite in a heterogeneous pool of non-conspecific individuals. In the second experiment, we made 2 × 2 combinations of species and studied the relationship between the proportion of individuals of a given species in a sample and the proportion of sequences retrieved to test whether DNA metabarcoding can reliably quantify the relative abundance of mites in a sample. Here we also tested the efficacy of degenerate primers (i.e., a mixture of similar primers that differ in one or several bases that are designed to increase the chance of annealing) and investigated the relationship between the number of mismatches and PCR success. Finally, we applied our DNA metabarcoding pipeline to a total of 6501 unidentified and unsorted feather mite individuals sampled from 380 European passerine birds belonging to 10 bird species (field test). Our results show that this proposed pipeline is suitable for correct identification and quantitative estimation of the relative abundance of feather mite species in complex samples, especially when dealing with a moderate number (> 30) of individuals per sample.

Environmental DNA detection of rare and invasive fish species in two Great Lakes tributaries

The extraction and characterization of DNA from aquatic environmental samples offers an alternative, noninvasive approach for the detection of rare species. Environmental DNA, coupled with PCR and next-generation sequencing ("metabarcoding"), has proven to be very sensitive for the detection of rare aquatic species. Our study used a custom-designed group-specific primer set and next-generation sequencing for the detection of three species at risk (Eastern Sand Darter, Ammocrypta pellucida; Northern Madtom, Noturus stigmosus; and Silver Shiner, Notropis photogenis), one invasive species (Round Goby, Neogobius melanostomus) and an additional 78 native species from two large Great Lakes tributary rivers in southern Ontario, Canada: the Grand River and the Sydenham River. Of 82 fish species detected in both rivers using capture-based and eDNA methods, our eDNA method detected 86.2% and 72.0% of the fish species in the Grand River and the Sydenham River, respectively, which included our four target species. Our analyses also identified significant positive and negative species co-occurrence patterns between our target species and other identified species. Our results demonstrate that eDNA metabarcoding that targets the fish community as well as individual species of interest provides a better understanding of factors affecting the target species spatial distribution in an ecosystem than possible with only target species data. Additionally, eDNA is easily implemented as an initial survey tool, or alongside capture-based methods, for improved mapping of species distribution patterns.

PCR cycles above routine numbers do not compromise high-throughput DNA barcoding results

High-throughput DNA barcoding has become essential in ecology and evolution, but some technical questions still remain. Increasing the number of PCR cycles above the routine 20-30 cycles is a common practice when working with old-type specimens, which provide little amounts of DNA, or when facing annealing issues with the primers. However, increasing the number of cycles can raise the number of artificial mutations due to polymerase errors. In this work, we sequenced 20 COI libraries in the Illumina MiSeq platform. Libraries were prepared with 40, 45, 50, 55, and 60 PCR cycles from four individuals belonging to four species of four genera of cephalopods. We found no relationship between the number of PCR cycles and the number of mutations despite using a nonproofreading polymerase. Moreover, even when using a high number of PCR cycles, the resulting number of mutations was low enough not to be an issue in the context of high-throughput DNA barcoding (but may still remain an issue in DNA metabarcoding due to chimera formation). We conclude that the common practice of increasing the number of PCR cycles should not negatively impact the outcome of a high-throughput DNA barcoding study in terms of the occurrence of point mutations.

Replacing Sanger with Next Generation Sequencing to improve coverage and quality of reference DNA barcodes for plants

We estimate the global BOLD Systems database holds core DNA barcodes (rbcL + matK) for about 15% of land plant species and that comprehensive species coverage is still many decades away. Interim performance of the resource is compromised by variable sequence overlap and modest information content within each barcode. Our model predicts that the proportion of species-unique barcodes reduces as the database grows and that 'false' species-unique barcodes remain >5% until the database is almost complete. We conclude the current rbcL + matK barcode is unfit for purpose. Genome skimming and supplementary barcodes could improve diagnostic power but would slow new barcode acquisition. We therefore present two novel Next Generation Sequencing protocols (with freeware) capable of accurate, massively parallel de novo assembly of high quality DNA barcodes of >1400 bp. We explore how these capabilities could enhance species diagnosis in the coming decades.

A quantitative protocol for DNA metabarcoding of springtails (Collembola)

We developed a novel protocol with superior quantitative analysis results for DNA metabarcoding of Collembola, a major soil microarthropod order. Degenerate PCR primers were designed for conserved regions in the mitochondrial cytochrome c oxidase subunit I (mtCOI) and 16S ribosomal RNA (mt16S) genes based on published collembolan mitogenomes. The best primer pair was selected based on its ability to amplify each gene, irrespective of the species. DNA was extracted from 10 natural communities sampled in a temperate forest (with typically 25-30 collembolan species per 10 soil samples) and 10 mock communities (with seven cultured collembolan species). The two gene regions were then amplified using the selected primers, ligated with adapters for 454 technology, and sequenced. Examination of the natural community samples showed that 32 and 36 operational taxonomic units (defined at a 90% sequence similarity threshold) were recovered from the mtCOI and mt16S data, respectively, which were comparable to the results of the microscopic identification of 25 morphospecies. Further, sequence abundances for each collembolan species from the mtCOI and mt16S data of the mock communities, after normalization by using a species as the internal control, showed good correlation with the number of individuals in the samples (R = 0.91-0.99), although relative species abundances within a mock community sample estimated from sequences were skewed from community composition in terms of the number of individuals or biomass of the species. Thus, this protocol enables the comparison of collembolan communities in a quantitative manner by metabarcoding.

Species mtDNA genetic diversity explained by infrapopulation size in a host-symbiont system

Understanding what shapes variation in genetic diversity among species remains a major challenge in evolutionary ecology, and it has been seldom studied in parasites and other host‐symbiont systems. Here, we studied mtDNA variation in a host‐symbiont non‐model system: 418 individual feather mites from 17 feather mite species living on 17 different passerine bird species. We explored how a surrogate of census size, the median infrapopulation size (i.e., the median number of individual parasites per infected host individual), explains mtDNA genetic diversity. Feather mite species genetic diversity was positively correlated with mean infrapopulation size, explaining 34% of the variation. As expected from the biology of feather mites, we found bottleneck signatures for most of the species studied but, in particular, three species presented extremely low mtDNA diversity values given their infrapopulation size. Their star‐like haplotype networks (in contrast with more reticulated networks for the other species) suggested that their low genetic diversity was the consequence of severe bottlenecks or selective sweeps. Our study shows for the first time that mtDNA diversity can be explained by infrapopulation sizes, and suggests that departures from this relationship could be informative of underlying ecological and evolutionary processes.