Summary DNA metabarcoding holds great promise for the assessment of macroinvertebrates in stream ecosystems. However, few large‐scale studies have compared the performance of DNA metabarcoding with that of routine morphological identification. We performed metabarcoding using four primer sets on macroinvertebrate samples from 18 stream sites across Finland. The samples were collected in 2013 and identified based on morphology as part of a Finnish stream monitoring program. Specimens were morphologically classified, following standardised protocols, to the lowest taxonomic level for which identification was feasible in the routine national monitoring. DNA metabarcoding identified more than twice the number of taxa than the morphology‐based protocol, and also yielded a higher taxonomic resolution. For each sample, we detected more taxa by metabarcoding than by the morphological method, and all four primer sets exhibited comparably good performance. Sequence read abundance and the number of specimens per taxon (a proxy for biomass) were significantly correlated in each sample, although the adjusted R2 values were low. With a few exceptions, the ecological status assessment metrics calculated from morphological and DNA metabarcoding datasets were similar. Given the recent reduction in sequencing costs, metabarcoding is currently approximately as expensive as morphology‐based identification. Using samples obtained in the field, we demonstrated that DNA metabarcoding can achieve comparable assessment results to current protocols relying on morphological identification. Thus, metabarcoding represents a feasible and reliable method to identify macroinvertebrates in stream bioassessment, and offers powerful advantage over morphological identification in providing identification for taxonomic groups that are unfeasible to identify in routine protocols. To unlock the full potential of DNA metabarcoding for ecosystem assessment, however, it will be necessary to address key problems with current laboratory protocols and reference databases.

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...

Species diversity of metazoan bulk samples can be rapidly assessed using cytochrome c oxidase I (COI) metabarcoding. However, in some applications often only degraded DNA is available, e.g. from poorly conserved museum specimens, environmental DNA (eDNA) filtered from water or gut content analyses. Here universal primer sets targeting only a short COI fragment are advantageous, as they often can still amplify short DNA fragments. Using PrimerMiner, we optimised two universal primer sets targeting freshwater macroinvertebrates based on NCBI and BOLD reference sequences. The fwh1 and fwh2 primer sets targeting a 178 and 205 bp region were tested in vitro by sequencing previously used freshwater macroinvertebrate mock communities as well as three monitoring samples from Romanian streams of unknown composition. They were further evaluated in silico for their suitability to amplify other insect groups. The fwh1 primer sets showed the most consistent amplification in silico and in vitro, detecting 92% of the taxa present in the mock communities, and allowing clear differentiation between the three macroinvertebrate communities from the Romanian streams. In silico analysis indicates that the short primers are likely to perform well even for non-freshwater insects. Comparing the performance of the new fwh1 primer sets to a highly degenerate primer set targeting a longer fragment (BF2+BR2) revealed that detection efficiency is slightly lower for the new primer set. Nevertheless, the shorter new primer pairs might be useful for studies that have to rely on degraded or poorly conserved DNA and thus be of importance for biomonitoring, conservation biological or molecular ecological studies. Furthermore, our study highlights the need for in silico evaluation of primer sets in order to detect design errors in primers (fwhR2) and find optimal universal primer sets for the target taxa of interest.

15Background. DNA metabarcoding is used to generate species composition data for entire communities. However, 16

Species diversity of metazoan bulk samples can be rapidly assessed using cytochrome c oxidase I (COI) metabarcoding. However, in some applications often only degraded DNA is available, e.g. from poorly conserved museum specimens, environmental DNA (eDNA) filtered from water or gut content analyses. Here universal primer sets targeting only a short COI fragment are advantageous, as they often can still amplify short DNA fragments. Using communities from the Romanian streams. In silico analysis indicates that the short primers are likely to perform well even for non-freshwater insects. Comparing the performance of the new fwh1 primer sets to a highly degenerate primer set targeting a longer fragment (BF2+BR2) revealed that detection efficiency is slightly lower for the new primer set.Nevertheless, the shorter new primer pairs might be useful for studies that have to rely on degraded or poorly conserved DNA and thus be of importance for biomonitoring, conservation biological or molecular ecological studies. Furthermore, our study highlights the need for in silico evaluation of primer sets in order to detect design errors in primers (fwhR2) and find optimal universal primer sets for the target taxa of interest.