Cited by 765Environmental DNA in Biodiversity StudiesInvertebrate Taxonomy and EcologyMicrobial Community Ecology and Physiology

The DOI matches its Crossref record; 59 of the 60 references it cites were found in a registry or an index; none of the 35 sources its claims cite is checked for support yet.

Can DNA-Based Ecosystem Assessments Quantify Species Abundance? Testing Primer Bias and Biomass—Sequence Relationships with an Innovative Metabarcoding Protocol

Open access

Abstract

Metabarcoding is an emerging genetic tool to rapidly assess biodiversity in ecosystems. It involves high-throughput sequencing of a standard gene from an environmental sample and comparison to a reference database. However, no consensus has emerged regarding laboratory pipelines to screen species diversity and infer species abundances from environmental samples. In particular, the effect of primer bias and the detection limit for specimens with a low biomass has not been systematically examined, when processing samples in bulk. We developed and tested a DNA metabarcoding protocol that utilises the standard cytochrome c oxidase subunit I (COI) barcoding fragment to detect freshwater macroinvertebrate taxa. DNA was extracted in bulk, amplified in a single PCR step, and purified, and the libraries were directly sequenced in two independent MiSeq runs (300-bp paired-end reads). Specifically, we assessed the influence of specimen biomass on sequence read abundance by sequencing 31 specimens of a stonefly species with known haplotypes spanning three orders of magnitude in biomass (experiment I). Then, we tested the recovery of 52 different freshwater invertebrate taxa of similar biomass using the same standard barcoding primers (experiment II). Each experiment was replicated ten times to maximise statistical power. The results of both experiments were consistent across replicates. We found a distinct positive correlation between species biomass and resulting numbers of MiSeq reads. Furthermore, we reliably recovered 83% of the 52 taxa used to test primer bias. However, sequence abundance varied by four orders of magnitudes between taxa despite the use of similar amounts of biomass. Our metabarcoding approach yielded reliable results for high-throughput assessments. However, the results indicated that primer efficiency is highly species-specific, which would prevent straightforward assessments of species abundance and biomass in a sample. Thus, PCR-based metabarcoding assessments of biodiversity should rely on presence-absence metrics.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

This study developed and tested a replicated COI metabarcoding protocol for freshwater macroinvertebrates, using bulk DNA extraction and Illumina MiSeq sequencing. It found a strong biomass–read relationship within a single species, but strong species-specific primer bias across taxa, concluding that PCR-based metabarcoding cannot reliably quantify species abundance or biomass and should be used for presence–absence assessments.

Key findings

  • In experiment I, all 31 *Dinocras cephalotes* specimens were recovered across ten replicates, including the smallest specimen at 0.023% (0.145 mg) of total biomass; specimen biomass and sequence abundance showed a highly significant positive linear correlation in all ten replicates (p<0.001, R² = 0.65).
  • In experiment II, 83% (43/52) of taxa were reliably recovered; 34 taxa were recovered in all ten replicates. Recovery by group was Ephemeroptera 7/8 (88%), Plecoptera 4/4 (100%), Trichoptera 13/15 (86%), Diptera 7/8 (88%), and Others 12/17 (71%).
  • Despite similar biomass, sequence abundance varied by four orders of magnitude among taxa, and there was no correlation between taxon biomass and number of recovered sequences. Primer efficiency was highly species-specific, preventing straightforward abundance or biomass inference.
  • LCO and HCO sequencing direction had small effects on final sequence abundance (2.7% and 4.4% differences in experiments I and II), though LCO was 34% more frequent in experiment I (p = 0.006) and no primer bias was detected in experiment II (p = 0.41).

Methods

  • Two controlled experiments with ten replicates each were sequenced on two independent MiSeq runs using 300-bp paired-end reads; experiment I used 31 *D. cephalotes* stonefly specimens spanning three orders of magnitude in biomass, and experiment II used 52 freshwater invertebrate taxa with similar tissue biomass.
  • DNA was extracted in bulk, amplified in a single PCR step with tagged fusion primers based on standard COI LCO1490/HCO2198, purified, and directly sequenced; experiment I used 1.4% and 20% tissue aliquots, while experiment II pooled roughly equal tissue amounts per taxon.
  • Bioinformatic analysis used demultiplexing, cutadapt primer removal, read concatenation, blastn against reference haplotypes for experiment I, and UPARSE OTU clustering at 97% identity with BOLD identification for experiment II; sequence counts were normalised to total sequence abundance and tissue weight.

Limitations

  • The paper states that PCR-based metabarcoding cannot accurately estimate species biomass or even abundance in diverse environmental samples because primer efficiency is highly species-specific; it recommends presence–absence metrics and notes that PCR-free approaches are needed for biomass estimates.
  • The paper states that only 83% of taxa were recovered with one universal COI primer pair, that barcode reference databases remain incomplete, that DNA-based assessments cannot assign biological traits, and that PCR replication effects were not evaluated due to a small number of replicates.
  • Evident from the text: the experiments are limited to freshwater macroinvertebrates from western Germany and a specific MiSeq/COI Folmer protocol; some morphotaxa were identified only to family or order and contained multiple MOTUs, so species-level amplification behaviour may be more variable than the morphotaxon-level summaries suggest.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1371/journal.pone.0130324

Cites

59 of 59 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Figures

Fig 1 Overview of the experimental setup of Experiment I and II. Two MiSeq runs were used to increase the reproducibility and reliability of our novel metabarcoding protocol. A Experiment I: Dinocras cephalotes specimens with different COI barcodes were used to determine the reproducibility of the protocol and the influence of biomass on sequence abundance. B Experiment II: Ten sets of 52 aquatic taxa were homogenised, and DNA was extracted and amplified to determine which taxa could be recovered with MiSeq using the "ready to load" primers developed in this study.
Fig 1 Overview of the experimental setup of Experiment I and II. Two MiSeq runs were used to increase the reproducibility and reliability of our novel metabarcoding protocol. A Experiment I: Dinocras cephalotes specimens with different COI barcodes were used to determine the reproducibility of the protocol and the influence of biomass on sequence abundance. B Experiment II: Ten sets of 52 aquatic taxa were homogenised, and DNA was extracted and amplified to determine which taxa could be recovered with MiSeq using the "ready to load" primers developed in this study.CC BY · PubMed Central
Fig 2 Results of Experiment I: Dependence of read abundance on specimen biomass. A Specimen weight (y-axis) is plotted against read abundance per specimen (x-axis) for all ten replicates. The linear regression (red line) was highly significant for all replicates with p<0.001. B Ratio of N5LCO/N5HCO sequence abundance for the mixed and four single DNA aliquots. The red line indicates the mean ratio.
Fig 2 Results of Experiment I: Dependence of read abundance on specimen biomass. A Specimen weight (y-axis) is plotted against read abundance per specimen (x-axis) for all ten replicates. The linear regression (red line) was highly significant for all replicates with p<0.001. B Ratio of N5LCO/N5HCO sequence abundance for the mixed and four single DNA aliquots. The red line indicates the mean ratio.CC BY · PubMed Central
Fig 3 Overview of taxa recovery in experiment II. Sequence abundances for the 52 morphologically identified taxa is shown in rows and the ten replicates used in the experiment in columns. Sequence abundance was normalised across the ten replicates and the amount of tissue used in each extraction. Sequence abundance of each specimens (morphotaxon) of the ten replicates is visualised by different shades of blue. If a field is i.e. half filled (50%) with the mid blue shade (= 1% of total sequences), the respective specimen represented 0.5% (50% of 1%) of the total sequences in that replicate. When no sequences or only a few sequences (below 0.003% of total abundance per replicate) were found for a specimen, it was scored as "No Hit," as indicated by an orange asterisk. On the right, K2P-corrected neighbour-joining (NJ) trees for each taxon, based on the most abundant sequence obtained for each specimen (calculated with MEGA6.06), are shown. MOTUs are defined by a 2% sequence difference based on the NJ tree.
Fig 3 Overview of taxa recovery in experiment II. Sequence abundances for the 52 morphologically identified taxa is shown in rows and the ten replicates used in the experiment in columns. Sequence abundance was normalised across the ten replicates and the amount of tissue used in each extraction. Sequence abundance of each specimens (morphotaxon) of the ten replicates is visualised by different shades of blue. If a field is i.e. half filled (50%) with the mid blue shade (= 1% of total sequences), the respective specimen represented 0.5% (50% of 1%) of the total sequences in that replicate. When no sequences or only a few sequences (below 0.003% of total abundance per replicate) were found for a specimen, it was scored as "No Hit," as indicated by an orange asterisk. On the right, K2P-corrected neighbour-joining (NJ) trees for each taxon, based on the most abundant sequence obtained for each specimen (calculated with MEGA6.06), are shown. MOTUs are defined by a 2% sequence difference based on the NJ tree.CC BY · PubMed Central

Verification

Claims the paper makes

44 Own result · 24 Own interpretation · 22 Not checked yet · 8 No basis given · 5 Own method

Verify all 35 references — free account

Show all 65 claims

  • Probably not a substantial claim. 23
  • It may state two separate claims as one. 23
38 of 38
  • The paper relies on this without support.

    Precise estimates of biomass with respect to specimen recovery and primer bias have not been performed.

    “Precise estimates of biomass with respect to specimen recovery and primer bias have not been performed.”
  • The paper relies on this without support.

    COI amplicons are typically similar in base composition.

    “COI amplicons are typically similar in base composition; therefore, three strategies were used to increase sequence diversity.”
  • The paper relies on this without support.

    Sanger sequencing requires that each specimen is processed individually in the laboratory, which is costly and extremely time-consuming for routine community assessments involving hundreds or thousands of specimens per sample.

    “However, Sanger sequencing requires that each specimen is processed individually in the laboratory, which is costly and extremely time-consuming for routine community assessments involving hundreds or thousands of specimens per sample.”
  • The paper relies on this without support.

    Benthic indicator organisms vary in biomass, and this variation depends on taxonomic group and life stage.

    “Benthic indicator organisms vary in biomass, and this variation depends on taxonomic group and life stage.”
  • The paper relies on this without support.

    Morphotaxa could often only be identified at family or order level, and several species per morphotaxon may occur in the sampled aquatic habitats even at the same habitat patch.

    “This was expected, as morphotaxa could often only be identified at family or order level, and several species per morphotaxon may occur in the sampled aquatic habitats even at the same habitat patch.”
  • The paper relies on this without support.

    Stochastic effects increase with reduced specimen biomass.

    “This result is expected because stochastic effects increase with reduced specimen biomass.”
  • The paper relies on this without support.

    Currently, the monitoring of freshwater ecosystems is based on abundance metrics, which cannot be generated using the metabarcoding solutions currently available.

    “Currently, the monitoring of freshwater ecosystems is based on abundance metrics, which cannot be generated using the metabarcoding solutions currently available.”
  • The paper relies on this without support.

    All assessment protocols rely on reference catalogues against which inventory data of a species from an ecosystem are compared.

    “All assessment protocols rely on reference catalogues against which inventory data of a species from an ecosystem are compared.”
  • The paper relies on this without support.

    Changes in species traits in a community (e.g. functional feeding groups, tolerance against pollution) are used as indicator values to evaluate the biological significance of inferred community changes.

    “In particular, changes in species traits in a community (e.g. functional feeding groups, tolerance against pollution) are used as indicator values to evaluate the biological significance of inferred community changes.”
  • The paper relies on this without support.

    DNA-based assessments cannot assign biological traits to species.

    “However, DNA-based assessments cannot assign biological traits to species.”
10 of 10

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

1 of 1

Citations

Cites

Show all 60 references

Cited by

Sorted by citations

Metadata

Authors
Elbrecht, V., & Leese, F.
Affiliations
Ruhr University Bochum (DE)
Corresponding author
Leese, Florian
Venue
PLOS ONE
Published
2015-07-08
Volume
10
Issue
7
Pages
e0130324
ISSN
1932-6203 (ISSN-L)
OpenAlex
W1551951435
Cited by
765
References
60
Citations by year
2026: 75 · 2025: 73 · 2024: 49 · 2023: 69 · 2022: 71 · 2021: 91 · 2020: 75 · 2019: 73 · 2018: 108 · 2017: 51 · 2016: 26 · 2015: 3
FWCI
56.82
Type
journal article
Language
en
Open access
Locations
  • PLoS ONE (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • Figshare (submittedVersion, cc-by)
Topics
  • Physical Sciences › Environmental Science › Ecology › Environmental DNA in Biodiversity Studies
  • Life Sciences › Agricultural and Biological Sciences › Ecology, Evolution, Behavior and Systematics › Invertebrate Taxonomy and Ecology
  • Physical Sciences › Environmental Science › Ecology › Microbial Community Ecology and Physiology
Keywords
DNA metabarcoding, primer bias, species abundance estimation, COI barcoding, freshwater macroinvertebrates, Plecoptera, MiSeq sequencing, biodiversity assessment
MeSH terms
Animals, Electron Transport Complex IV, Diptera, DNA, Fresh Water, Haplotypes, Species Specificity, Reproducibility of Results, Polymerase Chain Reaction, Ecosystem, DNA Primers, Biomass, Biodiversity, DNA Barcoding, Taxonomic, High-Throughput Nucleotide Sequencing, Odonata

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-10-06, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-06 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review