Comparing accuracy and biases of DNA metabarcoding, hybridization capture, and metagenomic sequencing for quantifying herbivore diets
Data files
Jul 10, 2026 version files 68.16 KB
-
MER_code.R
11.21 KB
-
plant_data.csv
2.04 KB
-
README.md
2.50 KB
-
WGS_HYB_MET_RRA_RGC.csv
52.41 KB
Abstract
DNA metabarcoding using relative read abundance (RRA) is commonly applied to estimate herbivore diet composition, yet its quantitative accuracy remains uncertain. We assessed taxonomic resolution and quantitative performance of RRA from DNA metabarcoding compared to metagenomic sequencing and hybridization capture, using deer scats from feeding trials and recreated diet samples using plant tissues. All methods recovered plant composition in recreated diets (R² = 0.59-0.82), indicating accurate scaling with biomass in the absence of digestion, with only minor bias from amplicon length in DNA metabarcoding. In contrast, RRA from scat samples performed poorly (R² < 0.01) across all methods largely due to differential plant digestibility. Correcting for digestibility, measured with acid detergent lignin and acid-insoluble ash, was strongly supported in mixed-effects models and improved prediction of dietary composition, although species-level variation remained. For metagenomic sequencing and hybridization capture, we also evaluated Relative Genome Coverage (RGC), a novel relative abundance metric quantifying the proportion of each plant’s chloroplast genome covered by mapped reads, normalized for genome length. RGC further improved correlations in recreated diets (R² = 0.82-0.84) and, with hybridization capture, largely overcame digestibility-related biases in scat samples (R² = 0.57) without correction. When such corrections are infeasible, hybridization capture with uncorrected RGC may achieve higher quantitative accuracy in scat samples. Our results provide practical guidance for improving molecular herbivore diet analysis and highlight the importance of accounting for digestion-related biases.
Dataset DOI: 10.5061/dryad.np5hqc06j
Description of the data and file structure
Data and R code used in "Comparing accuracy and biases of DNA metabarcoding, hybridization capture, and metagenomic sequencing for quantifying herbivore diets" Eriksson et al., (2025) Molecular Ecology Resources
Files and variables
File: MER_code.R
Description: 'MER_code.R' contains R code for all models. The two CSV files included in the deposit are the inputs for all analyses.
File: plant_data.csv
Description: Plant covariates
Variables
- Scientific: scientific name of each plant species and type of commercial pellet
- Common: common name of each plant species, commercial pellet types referred to as 'Pellet'
- Plant_part: part of the plant that was fed to the deer (e.g., leaves, fruit)
- Amplicon_GC: GC content (%) of trnL amplicon
- Amplicon_length: length of trnL amplicon (bp)
- Copies: chloroplast copy number
- 95%CI: 95% confidence interval of copy number
- corr_NDF: neutral detergent fiber (%)
- corr_ADL: acid detergent lignin (%)
- corr_AIA: acid-insoluble ash (%)
- corr_DMD: dry matter digestibility (%)
- Chloroplast_length: length of chloroplast genome (bp)
- Chloroplast_GC: GC content (%) of chloroplast genome
File: WGS_HYB_MET_RRA_RGC.csv
Description: Relative abundance data
Variables
- sample_id: Sample identification number
- total_read_count: Total number of sequence reads per sample
- count: Per base coverage per species ('NA' for DNA metabarcoding which did not use RGC)
- read_count: Number of sequence read per species
- scientific name: Scientific name of each plant species and pellet
- common: Common name of each plant species and pellet type
- diet: Feeding trial diet number
- sample_type: Deer scat or recreated diet sample
- propalldays: Proportion consumed by the deer or proportion used in recreated tissue mixes
- read_rra: Relative read abundance (RRA)
- norm_cov_rra: Relative genome coverage (RGC) ('NA' for DNA metabarcoding)
- method: DNA metabarcoding, Hybridization Capture, Metagenomic sequencing
- analysis_method: Method used for taxonomic assignment
- genus: Genus level identification
Code/software
All statistical analyses were conducted in R version 4.4.1. Code and required packages are included within the submission.
