Disparities in second‐generation DNA metabarcoding results exposed with accessible and repeatable workflows

Different second‐generation sequencing technologies may have taxon‐specific biases when DNA metabarcoding prey in predator faeces. Our major objective was to examine differences in prey recovery from bat guano across two different sequencing workflows using the same faecal DNA extracts. We compared results between the Ion Torrent PGM and the Illumina MiSeq with similar library preparations and the same analysis pipeline. We focus on repeatability and provide an R Notebook in an effort towards transparency for future methodological improvements. Full documentation of each step enhances the accessibility of our analysis pipeline. We tagged DNA from insectivorous bat faecal samples, targeted the arthropod cytochrome c oxidase I minibarcode region and sequenced the product on both second‐generation sequencing platforms. We developed an analysis pipeline with a high operational taxonomic unit (OTU) clustering threshold (i.e., ≥98.5%) followed by copy number filtering to avoid merging rare but genetically similar prey into the same OTUs. With this workflow, we detected 297 unique prey taxa, of which 74% were identified at the species level. Of these, 104 (35%) prey OTUs were detected by both platforms, 176 (59%) OTUs were detected by the Illumina MiSeq system only, and 17 (6%) OTUs were detected using the Ion Torrent system only. Costs were similar between platforms but the Illumina MiSeq recovered six times more reads and four additional insect orders than did Ion Torrent. The considerations we outline are particularly important for long‐term ecological monitoring; a more standardized approach will facilitate comparisons between studies and allow faster recognition of changes within ecological communities.

[1]  B. D. Burks The Mayflies, or Ephemeroptera, of Illinois , 1975 .

[2]  T. Kunz,et al.  An evaluation of fecal analysis for determining food habits of insectivorous bats , 1983 .

[3]  R. Putman Facts from faeces , 1984 .

[4]  T. Kunz,et al.  Food habits analysis of insectivorous bats. , 1988 .

[5]  R. Wayne,et al.  Facts from feces revisited. , 1997, Trends in ecology & evolution.

[6]  G. Carvalho Advances in molecular ecology , 1998 .

[7]  L. Waits,et al.  Noninvasive genetic sampling: look before you leap. , 1999, Trends in ecology & evolution.

[8]  C. Boesch,et al.  Factors affecting the amount of genomic DNA extracted from ape faeces and the identification of an improved sample storage method , 2004, Molecular ecology.

[9]  James R. Knight,et al.  Genome sequencing in microfabricated high-density picolitre reactors , 2005, Nature.

[10]  S. Sheppard,et al.  Advances in molecular ecology: tracking trophic links through predator–prey food‐webs , 2005 .

[11]  TIM M. BLACKBURN,et al.  Reproducibility and Repeatability in Ecology , 2006 .

[12]  P. Hebert,et al.  bold: The Barcode of Life Data System (http://www.barcodinglife.org) , 2007, Molecular ecology notes.

[13]  R. Sikes,et al.  Guidelines of the American Society of Mammalogists for the Use of Wild Mammals in Research , 2007 .

[14]  J. Fullard,et al.  The effectiveness of katydid (Neoconocephalus ensiger) song cessation as antipredator defence against the gleaning bat Myotis septentrionalis , 2008, Behavioral Ecology and Sociobiology.

[15]  G. A. Feldhamer,et al.  Prey Consumed by Eight Species of Insectivorous Bats from Southern Illinois , 2009 .

[16]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[17]  王丽华,et al.  国际生命条形码计划—DNA Barcoding , 2009 .

[18]  P. Hebert,et al.  Species on the menu of a generalist predator, the eastern red bat (Lasiurus borealis): using a molecular approach to detect arthropod prey , 2009, Molecular ecology.

[19]  P. Taberlet,et al.  DNA barcoding for ecologists. , 2009, Trends in ecology & evolution.

[20]  R. Hensinger Peer Reviewers , 2010, Economic Botany.

[21]  D. Heath,et al.  Revolution in food web analysis and trophic ecology: diet analysis by DNA and stable isotope analysis , 2010, Molecular ecology.

[22]  William A. Walters,et al.  QIIME allows analysis of high-throughput community sequencing data , 2010, Nature Methods.

[23]  R. W. Carter,et al.  Conservation in the dark? The information used to support management decisions , 2010 .

[24]  Wes McKinney,et al.  Data Structures for Statistical Computing in Python , 2010, SciPy.

[25]  Gareth Jones,et al.  Taxon‐specific PCR for DNA barcoding arthropod prey in bat faeces , 2011, Molecular ecology resources.

[26]  John A Darling,et al.  From molecules to management: adopting DNA-based methods for monitoring biological invasions in aquatic environments. , 2011, Environmental research.

[27]  Kristine Bohmann,et al.  Molecular Diet Analysis of Two African Free-Tailed Bats (Molossidae) Using High Throughput Sequencing , 2011, PloS one.

[28]  Elizabeth L Clare,et al.  High-throughput sequencing offers insight into mechanisms of resource partitioning in cryptic bat species , 2011, Ecology and evolution.

[29]  P. Hebert,et al.  BOLD’s role in barcode data management and analysis: a response , 2011 .

[30]  P. Racey Ecological and Behavioral Methods for the Study of Bats , 2011 .

[31]  Matthew B. Jones,et al.  Challenges and Opportunities of Open Data in Ecology , 2011, Science.

[32]  D. Higgins,et al.  Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega , 2011, Molecular systems biology.

[33]  P. Taberlet,et al.  Who is eating what: diet assessment using next generation sequencing , 2012, Molecular ecology.

[34]  P. Taberlet,et al.  Towards next‐generation biodiversity assessment using DNA metabarcoding , 2012, Molecular ecology.

[35]  E. Chapman,et al.  Identification of prey of Myotis septentrionalis using DNA-based techniques , 2012 .

[36]  T. Dallman,et al.  Performance comparison of benchtop high-throughput sequencing platforms , 2012, Nature Biotechnology.

[37]  Mehrdad Hajibabaei,et al.  Next‐generation sequencing technologies for environmental DNA research , 2012, Molecular ecology.

[38]  R. Swihart,et al.  The Hardwood Ecosystem Experiment: a framework for studying responses to forest management , 2013 .

[39]  Kristine Bohmann,et al.  Second generation sequencing and morphological faecal analysis reveal unexpected foraging behaviour by Myotis nattereri (Chiroptera, Vespertilionidae) in winter , 2014, Frontiers in Zoology.

[40]  N. Wahlberg,et al.  Next Generation Sequencing of Fecal DNA Reveals the Dietary Diversity of the Widespread Insectivorous Predator Daubenton’s Bat (Myotis daubentonii) in Southwestern Finland , 2013, PloS one.

[41]  A. Kurta,et al.  Analysis of DNA from Feces to Identify Prey of Big Brown Bats (Eptesicus fuscus) Caught in Apple Orchards , 2013 .

[42]  Philip Hugenholtz,et al.  Shining a Light on Dark Sequencing: Characterising Errors in Ion Torrent PGM Data , 2013, PLoS Comput. Biol..

[43]  Rebecca A. Kalb,et al.  The hardwood ecosystem experiment: goals, design, and implementation , 2013 .

[44]  J. Kruschke Bayesian estimation supersedes the t test. , 2013, Journal of experimental psychology. General.

[45]  P. Taberlet,et al.  DNA metabarcoding multiplexing and validation of data accuracy for diet assessment: application to omnivorous diet , 2014, Molecular ecology resources.

[46]  Stephen J. Salipante,et al.  Performance Comparison of Illumina and Ion Torrent Next-Generation Sequencing Platforms for 16S rRNA-Based Bacterial Community Profiling , 2014, Applied and Environmental Microbiology.

[47]  James Haile,et al.  Who's for dinner? High‐throughput sequencing reveals bat dietary differentiation in a biodiversity hotspot where prey taxonomy is largely undescribed , 2014, Molecular ecology.

[48]  B. Siemers,et al.  An integrative approach to detect subtle trophic niche differentiation in the sympatric trawling bat species Myotis dasycneme and Myotis daubentonii , 2014, Molecular ecology.

[49]  A. Kurta,et al.  Species-level analysis of diets of two mormoopid bats from Puerto Rico , 2014 .

[50]  Elizabeth L Clare,et al.  An inordinate fondness for beetles? Variation in seasonal dietary preferences of night‐roosting big brown bats (Eptesicus fuscus) , 2014, Molecular ecology.

[51]  M. B. Fenton,et al.  Resource partitioning by insectivorous bats in Jamaica , 2014, Molecular ecology.

[52]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[53]  Frédéric Mahé,et al.  Swarm: robust and fast clustering method for amplicon-based studies , 2014, PeerJ.

[54]  L. Weyrich,et al.  Environmental metabarcodes for insects: in silico PCR reveals potential for taxonomic bias , 2014, Molecular ecology resources.

[55]  François Fabianek,et al.  The diet of Myotis lucifugus across Canada: assessing foraging quality and diet variability , 2014, Molecular ecology.

[56]  Elizabeth L Clare,et al.  Molecular detection of trophic interactions: emerging trends, distinct advantages, significant considerations and conservation applications , 2014, Evolutionary applications.

[57]  J. A. Encarnação,et al.  Importance of multi-dimensional analyses of resource partitioning in highly mobile species assemblages , 2015, Population Ecology.

[58]  P. Taberlet,et al.  Replication levels, false presences and the estimation of the presence/absence from eDNA metabarcoding data , 2015, Molecular ecology resources.

[59]  P. Lim,et al.  DNA metabarcoding of insects and allies: an evaluation of primers and pipelines , 2015, Bulletin of Entomological Research.

[60]  Kristine Bohmann,et al.  Tag jumps illuminated – reducing sequence‐to‐sample misidentifications in metabarcoding studies , 2015, Molecular ecology resources.

[61]  John G Kenny,et al.  A comprehensive benchmarking study of protocols and sequencing platforms for 16S rRNA community profiling , 2016, BMC Genomics.

[62]  Antton Alberdi,et al.  Trophic requirements beyond foraging habitats: The importance of prey source habitats in bat conservation , 2015 .

[63]  G. McCracken,et al.  Crop pests eaten by bats in organic pecan orchards , 2015 .

[64]  E. Clare,et al.  Dietary overlap and seasonality in three species of mormoopid bats from a tropical dry forest , 2015, Molecular ecology.

[65]  C. Quince,et al.  Insight into biases and sequencing errors for amplicon sequencing with the Illumina MiSeq platform , 2015, Nucleic acids research.

[66]  Stephan Lewandowsky,et al.  The Peer Reviewers' Openness Initiative: incentivizing open research practices through peer review , 2016, Royal Society Open Science.

[67]  Pierre Taberlet,et al.  How to limit false positives in environmental DNA and metabarcoding? , 2016, Molecular ecology resources.

[68]  Vanessa A. Mata,et al.  Female dietary bias towards large migratory moths in the European free-tailed bat (Tadarida teniotis) , 2016, Biology Letters.

[69]  Umer Zeeshan Ijaz,et al.  Illumina error profiles: resolving fine-scale variation in metagenomic sequencing data , 2016, BMC Bioinformatics.

[70]  Stephen R. Piccolo,et al.  Tools and techniques for computational reproducibility , 2016, GigaScience.

[71]  N. Wahlberg,et al.  What you need is what you eat? Prey selection by the bat Myotis daubentonii , 2016, Molecular ecology.

[72]  Jj Allaire,et al.  Dynamic Documents for R , 2016 .

[73]  Sujeevan Ratnasingham,et al.  Counting animal species with DNA barcodes: Canadian insects , 2016, Philosophical Transactions of the Royal Society B: Biological Sciences.

[74]  Frédéric J. J. Chain,et al.  The effects of parameter choice on defining molecular operational taxonomic units and resulting ecological analyses of metabarcoding data. , 2016, Genome.

[75]  Anna Krystalli,et al.  Striving for transparent and credible research: practical guidelines for behavioral ecologists , 2017, Behavioral ecology : official journal of the International Society for Behavioral Ecology.

[76]  G. Guillera‐Arroita,et al.  Dealing with false‐positive and false‐negative errors about species occurrence at multiple levels , 2017 .

[77]  L. Waits,et al.  Comparing morphological and molecular diet analyses and fecal DNA sampling protocols for a terrestrial carnivore , 2017 .

[78]  J. Boyles,et al.  Illuminating prey selection in an insectivorous bat community exposed to artificial light at night , 2018 .

[79]  Antton Alberdi,et al.  Scrutinizing key steps for reliable metabarcoding of environmental samples , 2018 .