A pipeline for metabarcoding and diet analysis from fecal samples developed for a small semi-aquatic mammal

Metabarcoding allows the genetic analysis of pooled samples of various sources. It is becoming popular in the study of animal diet, especially because it allows the analysis of the composition of feces without the need of handling animals. In this work, we studied the diet of the Pyrenean desman (Galemys pyrenaicus), a small semi-aquatic mammal endemic to the Iberian Peninsula and the Pyrenees, by sequencing COI minibarcodes from feces using next-generation sequencing techniques. For the identification of assembled sequences, we employed a tree-based identification method that used a reference tree of sequences of freshwater organisms. The comparison of freshly collected fecal samples and older samples showed that fresh samples produced significantly more sequencing reads. They also rendered more operational taxonomical units (OTUs), but not significantly. Our analyses of 41 samples identified 224 OTUs corresponding to species of the reference tree. Ephemeroptera, Diptera excl. Chironomidae, and Chironomidae were the most highly represented groups in terms of reads as well as samples. Other groups of freshwater organisms detected were Plecoptera, Trichoptera, Neuropteroida, Coleoptera, Crustacea, and Annelida. Our results are largely in line with previous morphological and genetic studies on the diet of the Pyrenean desman, but allowed the identification of a higher diversity of OTUs in each sample. Additionally, the bioinformatic pipeline we developed for deep sequencing of fecal samples will enable the quantitative analysis of the diet of this and other species, which can be highly useful to determine their ecological requirements.

[1]  P. Jaccard,et al.  Etude comparative de la distribution florale dans une portion des Alpes et des Jura , 1901 .

[2]  P. Jaccard.  THE DISTRIBUTION OF THE FLORA IN THE ALPINE ZONE.1 , 1912 .

[3]  C. E. SHANNON,et al.  A mathematical theory of communication , 1948, MOCO.

[4]  R. Vrijenhoek,et al.  DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. , 1994, Molecular marine biology and biotechnology.

[5]  E. Castién,et al.  Diet of Galemys Pyrenaicus (Geoffroy, 1811 ) in the North of the Iberian Peninsula , 1994 .

[6]  Thomas L. Madden,et al.  Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. , 1997, Nucleic acids research.

[7]  J A Eisen,et al.  Phylogenomics: improving functional predictions for uncharacterized genes by evolutionary analysis. , 1998, Genome research.

[8]  J. Hemingway,et al.  Can multiple‐copy sequences of prey DNA be detected amongst the gut contents of invertebrate predators? , 1999, Molecular ecology.

[9]  L. Koski,et al.  The Closest BLAST Hit Is Often Not the Nearest Neighbor , 2001, Journal of Molecular Evolution.

[10]  P. Boixader,et al.  Distribució del Galemys pyrenaicus (Insectivora, Talpidae) a Catalunya , 2001 .

[11]  K. Katoh,et al.  MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. , 2002, Nucleic acids research.

[12]  Jeremy R. deWaard,et al.  Biological identifications through DNA barcodes , 2003, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[13]  D. Janzen,et al.  Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[14]  J. Santamarina.  Feeding ecology of a vertebrate assemblage inhabiting a stream of NW Spain (Riobo; Ulla basin) , 1993, Hydrobiologia.

[15]  B. Deagle,et al.  Group‐specific polymerase chain reaction for DNA‐based analysis of species diversity and identity in dietary samples , 2004, Molecular ecology.

[16]  L. Waits,et al.  NONINVASIVE GENETIC SAMPLING TOOLS FOR WILDLIFE BIOLOGISTS: A REVIEW OF APPLICATIONS AND RECOMMENDATIONS FOR ACCURATE DATA COLLECTION , 2005 .

[17]  P. Hebert,et al.  DNA barcoding for effective biodiversity assessment of a hyperdiverse arthropod group: the ants of Madagascar , 2005, Philosophical Transactions of the Royal Society B: Biological Sciences.

[18]  P. Hebert,et al.  Biological identifications of mayflies (Ephemeroptera) using DNA barcodes , 2005, Journal of the North American Benthological Society.

[19]  D. Strayer.  Challenges for freshwater invertebrate conservation , 2006, Journal of the North American Benthological Society.

[20]  C. Vilà,et al.  Reliability of noninvasive genetic census of otters compared to field censuses , 2007, Conservation Genetics.

[21]  S. Johansen,et al.  Halibut mitochondrial genomes contain extensive heteroplasmic tandem repeat arrays involved in DNA recombination , 2008, BMC Genomics.

[22]  J. Landry,et al.  A universal DNA mini-barcode for biodiversity analysis , 2008, BMC Genomics.

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

[24]  P. Hebert,et al.  DNA barcoding: how it complements taxonomy, molecular phylogenetics and population genetics. , 2007, Trends in genetics : TIG.

[25]  Kazutaka Katoh,et al.  Improved accuracy of multiple ncRNA alignment by incorporating structural information into a MAFFT-based framework , 2008, BMC Bioinformatics.

[26]  S. Ratnasingham,et al.  BOLD : The Barcode of Life Data System (www.barcodinglife.org) , 2007 .

[27]  Alexander F. Auch,et al.  MEGAN analysis of metagenomic data. , 2007, Genome research.

[28]  Kathleen Marchal,et al.  Evaluation of time profile reconstruction from complex two-color microarray designs , 2008, BMC Bioinformatics.

[29]  Tu Bao Ho,et al.  Finding microRNA regulatory modules in human genome using rule induction , 2008, BMC Bioinformatics.

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

[31]  D. Strayer,et al.  Freshwater biodiversity conservation: recent progress and future challenges , 2010, Journal of the North American Benthological Society.

[32]  Robert C. Edgar,et al.  BIOINFORMATICS APPLICATIONS NOTE , 2001 .

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

[34]  Tao Jiang,et al.  SEED: efficient clustering of next-generation sequences , 2011, Bioinform..

[35]  C. Mora,et al.  How Many Species Are There on Earth and in the Ocean? , 2011, PLoS biology.

[36]  Denis Krompass,et al.  Performance, Accuracy, and Web Server for Evolutionary Placement of Short Sequence Reads under Maximum Likelihood , 2011, Systematic biology.

[37]  M. Némoz,et al.  A French Conservation Action Plan for the Pyrenean Desman Galemys pyrenaicus , 2011 .

[38]  J. Luque-Larena,et al.  New insights into social and space use behaviour of the endangered Pyrenean desman (Galemys pyrenaicus) , 2012, European Journal of Wildlife Research.

[39]  D. Baird,et al.  Building freshwater macroinvertebrate DNA-barcode libraries from reference collection material: formalin preservation vs specimen age , 2011, Journal of the North American Benthological Society.

[40]  P. Taberlet,et al.  Carnivore diet analysis based on next‐generation sequencing: application to the leopard cat (Prionailurus bengalensis) in Pakistan , 2012, Molecular ecology.

[41]  J. Castresana,et al.  Phylogeography and postglacial expansion of the endangered semi-aquatic mammal Galemys pyrenaicus , 2013, BMC Evolutionary Biology.

[42]  L. Orlando,et al.  Meta‐barcoding of ‘dirt’ DNA from soil reflects vertebrate biodiversity , 2012, Molecular ecology.

[43]  C. Wiuf,et al.  Monitoring endangered freshwater biodiversity using environmental DNA. , 2012, Molecular ecology.

[44]  W. E. Harris,et al.  An emergent science on the brink of irrelevance: a review of the past 8 years of DNA barcoding , 2012, Molecular ecology resources.

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

[46]  Frederick Albert Matsen IV,et al.  A Format for Phylogenetic Placements , 2012, PloS one.

[47]  N. Baeshen,et al.  Biological Identifications Through DNA Barcodes , 2012 .

[48]  Z. Nagy,et al.  Reliable DNA Barcoding Performance Proved for Species and Island Populations of Comoran Squamate Reptiles , 2013, PloS one.

[49]  Z. Nagy,et al.  DNA barcoding : a practical tool for fundamental and applied biodiversity research , 2013 .

[50]  Alexandros Stamatakis,et al.  RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies , 2014, Bioinform..

[51]  F. Altermatt,et al.  Transport Distance of Invertebrate Environmental DNA in a Natural River , 2014, PloS one.

[52]  J. Gómez‐Zurita,et al.  Automated DNA‐based plant identification for large‐scale biodiversity assessment , 2015, Molecular ecology resources.

[53]  S. Aulagnier,et al.  Integrating hydrological features and genetically validated occurrence data in occupancy modeling of an endemic and endangered semi-aquatic mammal species, Galemys pyrenaicus, in a Pyrenean catchment , 2015 .

[54]  M. Wilkinson,et al.  Quantitative evaluation of bias in PCR amplification and next-generation sequencing derived from metabarcoding samples , 2015, Analytical and Bioanalytical Chemistry.

[55]  D. P. Paula,et al.  Detection and decay rates of prey and prey symbionts in the gut of a predator through metagenomics , 2015, Molecular ecology resources.

[56]  S. Aulagnier,et al.  A new method to identify the endangered Pyrenean desman (Galemys pyrenaicus) and to study its diet, using next generation sequencing from faeces , 2015 .

[57]  J. Piñol,et al.  Universal and blocking primer mismatches limit the use of high‐throughput DNA sequencing for the quantitative metabarcoding of arthropods , 2015, Molecular ecology resources.

[58]  David L. Erickson,et al.  DNA metabarcoding illuminates dietary niche partitioning by African large herbivores , 2015, Proceedings of the National Academy of Sciences.

[59]  D. Hering,et al.  www.freshwaterecology.info – An online tool that unifies, standardises and codifies more than 20,000 European freshwater organisms and their ecological preferences ☆ , 2015 .

[60]  A. Hausmann,et al.  A comprehensive DNA barcode database for Central European beetles with a focus on Germany: adding more than 3500 identified species to BOLD , 2015, Molecular ecology resources.

[61]  A. Vogler,et al.  Comparing the effectiveness of metagenomics and metabarcoding for diet analysis of a leaf‐feeding monkey (Pygathrix nemaeus) , 2015, Molecular ecology resources.

[62]  Samuel Arvidsson,et al.  Species Identification in Malaise Trap Samples by DNA Barcoding Based on NGS Technologies and a Scoring Matrix , 2016, PloS one.

[63]  J. Lennon,et al.  Scaling laws predict global microbial diversity , 2016, Proceedings of the National Academy of Sciences.

[64]  B. Deagle,et al.  Quantitative DNA metabarcoding: improved estimates of species proportional biomass using correction factors derived from control material , 2016, Molecular ecology resources.

[65]  F. Glaw,et al.  DNA barcoding of crickets, katydids and grasshoppers (Orthoptera) from Central Europe with focus on Austria, Germany and Switzerland , 2017, Molecular ecology resources.

[66]  L. Buisson,et al.  Comparison of diet and prey selectivity of the Pyrenean desman and the Eurasian water shrew using next-generation sequencing methods , 2017 .

[67]  F. Leese,et al.  A DNA barcode library for Germany′s mayflies, stoneflies and caddisflies (Ephemeroptera, Plecoptera and Trichoptera) , 2017, Molecular ecology resources.

[68]  S. Aulagnier,et al.  Novel insights into the diet of the Pyrenean desman (Galemys Pyrenaicus) using next-generation sequencing molecular analyses , 2017, Journal of Mammalogy.