Estimating intraspecific genetic diversity from community DNA metabarcoding data

Abstract

BackgroundDNA metabarcoding is used to generate species composition data for entire communities. However, sequencing errors in high-throughput sequencing instruments are fairly common, usually requiring reads to be clustered into operational taxonomic units (OTUs), losing information on intraspecific diversity in the process. While Cytochrome c oxidase subunit I (COI) haplotype information is limited in resolving intraspecific diversity it is nevertheless often useful e.g. in a phylogeographic context, helping to formulate hypotheses on taxon distribution and dispersal.MethodsThis study combines sequence denoising strategies, normally applied in microbial research, with additional abundance-based filtering to extract haplotype information from freshwater macroinvertebrate metabarcoding datasets. This novel approach was added to the R package “JAMP” and can be applied to COI amplicon datasets. We tested our haplotyping method by sequencing (i) a single-species mock community composed of 31 individuals with 15 different haplotypes spanning three orders of magnitude in biomass and (ii) 18 monitoring samples each amplified with four different primer sets and two PCR replicates.ResultsWe detected all 15 haplotypes of the single specimens in the mock community with relaxed filtering and denoising settings. However, up to 480 additional unexpected haplotypes remained in both replicates. Rigorous filtering removes most unexpected haplotypes, but also can discard expected haplotypes mainly from the small specimens. In the monitoring samples, the different primer sets detected 177–200 OTUs, each containing an average of 2.40–3.30 haplotypes per OTU. The derived intraspecific diversity data showed population structures that were consistent between replicates and similar between primer pairs but resolution depended on the primer length. A closer look at abundant taxa in the dataset revealed various population genetic patterns, e.g. the stonefly Taeniopteryx nebulosa and the caddisfly Hydropsyche pellucidula showed a distinct north–south cline with respect to haplotype distribution, while the beetle Oulimnius tuberculatus and the isopod Asellus aquaticus displayed no clear population pattern but differed in genetic diversity.DiscussionWe developed a strategy to infer intraspecific genetic diversity from bulk invertebrate metabarcoding data. It needs to be stressed that at this point this metabarcoding-informed haplotyping is not capable of capturing the full diversity present in such samples, due to variation in specimen size, primer bias and loss of sequence variants with low abundance. Nevertheless, for a high number of species intraspecific diversity was recovered, identifying potentially isolated populations and taxa for further more detailed phylogeographic investigation. While we are currently lacking large-scale metabarcoding datasets to fully take advantage of our new approach, metabarcoding-informed haplotyping holds great promise for biomonitoring efforts that not only seek information about species diversity but al...

From metabarcoding to metaphylogeography: separating the wheat from the chaff

“…Both the simulation approach and the analysis of the real data set pointed to an α‐value of 5 in the denoising step, which was also the optimal value selected in Elbrecht et al. (). Whether this value can be taken as a general rule of thumb or not will require analyses of more data sets.…”

“…In this study, we suggest an operational approach based on the stabilization of the entropy ratio to guide the cleaning procedures. Both the simulation approach and the analysis of the real data set pointed to an avalue of 5 in the denoising step, which was also the optimal value selected in Elbrecht et al (2018a). Whether this value can be taken as a general rule of thumb or not will require analyses of more data sets.…”

“…Finally, even if sequences retained were mostly correct, they can still include a number of nontarget variants due to heteroplasmy or numts (Elbrecht et al 2018a). However, numts tend to accumulate mutations resulting in stop codons (Song et al 2008).…”

Invasion genetics from eDNA and thousands of larvae: A targeted metabarcoding assay that distinguishes species and population variation of zebra and quagga mussels

“…Although those studies detected dreissenids (De Ventura et al, 2017;Egan et al, 2015;Gingera et al, 2017;Ram et al, 2011), their methodology was unable to differentiate between zebra and quagga mussel species using a single assay and could not address population-level variation (e.g., haplotypes). Targeted metabarcoding assays with high-throughput sequencing (HTS) and bioinformatic processing can be designed to simultaneously resolve entire communities to species level, as well as distinguish population genetic variation (Elbrecht, Vamos, Steinke, & Leese, 2018;Parsons, Everett, Dahlheim, & Park, 2018;Sigsgaard et al, 2016;Stepien, Snyder, & Elz, 2019).…”

Cobble community DNA as a tool to monitor patterns of biodiversity within kelp forest ecosystems

“…For the first time in history, abalone fishing has been closed state‐wide since 2017, and monitoring for recovery through recruitment will be critical. Further, utilizing metabarcoding‐based haplotypes to investigate population genetic patterns allows identifying the presence and distribution of genetically distinguishable groups (Elbrecht, Vamos, Steinke, & Leese, ).…”