A Next-Generation Sequencing Approach to River Biomonitoring Using Benthic Diatoms

Abstract: Diatoms are main bioindicators used to assess the ecological quality of rivers, but their identification is difficult and time-consuming. Next Generation Sequencing (NGS) can be used to study communities of microorganisms, so we carried out a test of the reliability of 454 pyrosequencing for estimating diatom inventories in environmental samples. We used small subunit ribosomal deoxyribonucleic acid (SSU rDNA), ribulose-1, 5-bisphosphate carboxylase (rbcL), and cytochrome oxidase I (COI) markers and examined reference libraries to define thresholds between the intra- and interspecific and intra- and intergeneric genetic distances. Based on tests of 1 mock community, we used a threshold of 99% identity for SSU rDNA and rbcL sequences to study freshwater diatoms at the species level. We applied 454 pyrosequencing to 4 contrasting environmental samples (with one in duplicate), assigned taxon names to environmental sequences, and compared the qualitative and quantitative molecular inventories to those obtained by microscopy. Species richness detected by microscopy was always higher than that detected by pyrosequencing. Some morphologically detected taxa may have been persistent frustules from dead cells. Some taxa detected by molecular analysis were not detected by morphology and vice versa. The main source of divergence appears to be inadequate taxonomic coverage in DNA reference libraries. Only a small percentage of species (but almost all genera) in morphological inventories were included in DNA reference libraries. DNA reference libraries contained a smaller percentage of species from tropical (27.1–38.1%) than from temperate samples (53.7–77.8%). Agreement between morphological and molecular inventories was better for species with relative abundance >1% than for rare species. The rbcL marker appeared to provide more reproducible results (94.9% species similarity between the 2 duplicates) and was very useful for molecular identification, but procedural standardization is needed. The water-quality ranking assigned to a site via the Pollution Sensitivity diatom index was the same whether calculated with molecular or morphological data. Pyrosequencing is a promising approach for detecting all species, even rare ones, once reference libraries have been developed.

[1]  Joel Closier,et al.  DIRAC: a community grid solution , 2008 .

[2]  I. Kaczmarska,et al.  Barcoding of diatoms: nuclear encoded ITS revisited. , 2010, Protist.

[3]  O Hammer-Muntz,et al.  PAST: paleontological statistics software package for education and data analysis version 2.09 , 2001 .

[4]  S. Weisberg,et al.  Performance of Two Southern California Benthic Community Condition Indices Using Species Abundance and Presence-Only Data: Relevance to DNA Barcoding , 2012, PloS one.

[5]  J. Chun,et al.  Assessment of soil fungal communities using pyrosequencing , 2010, The Journal of Microbiology.

[6]  Martin F. Polz,et al.  Bias in Template-to-Product Ratios in Multitemplate PCR , 1998, Applied and Environmental Microbiology.

[7]  R. Giblin-Davis,et al.  Linking operational clustered taxonomic units (OCTUs) from parallel ultra sequencing (PUS) to nematode species , 2010 .

[8]  T. Ohama,et al.  Phylogenetic analysis of diatom coxI genes and implications of a fluctuating GC content on mitochondrial genetic code evolution , 2000, Current Genetics.

[9]  F. Rimet,et al.  Biomonitoring river diatoms: Implications of taxonomic resolution , 2012 .

[10]  I. Kaczmarska,et al.  Barcoding diatoms: Is there a good marker? , 2009, Molecular ecology resources.

[11]  Maria E. Asplund,et al.  Quantification of Diatom and Dinoflagellate Biomasses in Coastal Marine Seawater Samples by Real-Time PCR , 2008, Applied and Environmental Microbiology.

[12]  L. Medlin,et al.  Evolution of the diatoms (Bacillariophyta). IV. A reconstruction of their age from small subunit rRNA coding regions and the fossil record. , 1996, Molecular phylogenetics and evolution.

[13]  D. Mann,et al.  An assessment of potential diatom "barcode" genes (cox1, rbcL, 18S and ITS rDNA) and their effectiveness in determining relationships in Sellaphora (Bacillariophyta). , 2007, Protist.

[14]  Piet F. M. Verdonschot,et al.  Uncertainty in Diatom Assessment: Sampling, Identification and Counting Variation , 2006, Hydrobiologia.

[15]  J. Gong,et al.  Extremely high copy numbers and polymorphisms of the rDNA operon estimated from single cell analysis of oligotrich and peritrich ciliates. , 2013, Protist.

[16]  R. Knight,et al.  Fast UniFrac: Facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data , 2009, The ISME Journal.

[17]  F. Rimet,et al.  Next‐generation sequencing to inventory taxonomic diversity in eukaryotic communities: a test for freshwater diatoms , 2013, Molecular ecology resources.

[18]  M. Meybeck.  Global analysis of river systems: from Earth system controls to Anthropocene syndromes. , 2003, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[19]  F. Sylvestre,et al.  Enregistrement de la variabilité hydroclimatique au cours des 6500 dernières années sur le littoral de Cabo Frio (Rio de Janeiro, Brésil) , 2006 .

[20]  R. A. Zampella,et al.  Distribution of Diatoms in Relation to Land Use and pH in Blackwater Coastal Plain Streams , 2007, Environmental management.

[21]  E. Antoine,et al.  ISSR as new markers for genetic characterization and evaluation of relationships among phytoplankton , 2004, Journal of Applied Phycology.

[22]  Marina G. Potapova,et al.  Diatoms from the genus Achnanthidium in flowing waters of the Appalachian Mountains (North America): Ecology, distribution and taxonomic notes , 2007 .

[23]  Ø. Hammer,et al.  PAST: PALEONTOLOGICAL STATISTICAL SOFTWARE PACKAGE FOR EDUCATION AND DATA ANALYSIS , 2001 .

[24]  S. Sabater,et al.  Determination of the biological diatom index (IBD NF T 90–354): results of an intercomparison exercise , 2002, Journal of Applied Phycology.

[25]  W. Stam,et al.  Quantification of algal endosymbionts (Symbiodinium) in coral tissue using real‐time PCR , 2009, Molecular ecology resources.

[26]  M. Kelly,et al.  A comparison of national approaches to setting ecological status boundaries in phytobenthos assessment for the European Water Framework Directive: results of an intercalibration exercise , 2009, Hydrobiologia.

[27]  D. Mann,et al.  DNA barcoding for species identification and discovery in diatoms , 2010 .

[28]  Frédéric Rimet,et al.  Recent views on river pollution and diatoms , 2012, Hydrobiologia.

[29]  D. Mann,et al.  Barcoding diatoms: exploring alternatives to COI-5P. , 2011, Protist.

[30]  D. Vaulot,et al.  Mapping of picoeucaryotes in marine ecosystems with quantitative PCR of the 18S rRNA gene. , 2005, FEMS microbiology ecology.

[31]  P. Hamilton,et al.  MORPHOLOGICAL AND ECOLOGICAL VARIATION WITHIN THE ACHNANTHIDIUM MINUTISSIMUM (BACILLARIOPHYCEAE) SPECIES COMPLEX 1 , 2007 .

[32]  M. Nei,et al.  MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. , 2011, Molecular biology and evolution.

[33]  T. Stoeck,et al.  Depicting more accurate pictures of protistan community complexity using pyrosequencing of hypervariable SSU rRNA gene regions. , 2011, Environmental microbiology.

[34]  M. Trauth,et al.  Hidden diversity in diatoms of Kenyan Lake Naivasha: a genetic approach detects temporal variation , 2012, Molecular ecology.

[35]  F. Sylvestre,et al.  A high-resolution diatom reconstruction between 21,000 and 17,4000 14C yr BP from the southern Bolivian Altiplano (18–23°S) , 2002 .

[36]  M. Sogin,et al.  The characterization of enzymatically amplified eukaryotic 16S-like rRNA-coding regions. , 1988, Gene.

[37]  A. Chao.  Estimating the population size for capture-recapture data with unequal catchability. , 1987, Biometrics.

[38]  G. Olsen,et al.  The small-subunit ribosomal RNA gene sequences from the hypotrichous ciliates Oxytricha nova and Stylonychia pustulata. , 1985, Molecular biology and evolution.

[39]  L. Tedersoo,et al.  454 Pyrosequencing and Sanger sequencing of tropical mycorrhizal fungi provide similar results but reveal substantial methodological biases. , 2010, The New phytologist.

[40]  R. Andersen,et al.  A MOLECULAR PHYLOGENY OF THE HETEROKONT ALGAE BASED ON ANALYSES OF CHLOROPLAST‐ENCODED rbcL SEQUENCE DATA 1 , 1997 .

[41]  Andrew P. Martin.  Phylogenetic Approaches for Describing and Comparing the Diversity of Microbial Communities , 2002, Applied and Environmental Microbiology.

[42]  Regine Jahn,et al.  Barcoding diatoms: evaluation of the V4 subregion on the 18S rRNA gene, including new primers and protocols , 2011, Organisms Diversity & Evolution.

[43]  Mehrdad Hajibabaei,et al.  Assessing biodiversity of a freshwater benthic macroinvertebrate community through non-destructive environmental barcoding of DNA from preservative ethanol , 2012, BMC Ecology.

[44]  François Pompanon,et al.  Persistence of Environmental DNA in Freshwater Ecosystems , 2011, PloS one.

[45]  S. Carpenter,et al.  Global Consequences of Land Use , 2005, Science.

[46]  Christophe Klopp,et al.  Assessment of replicate bias in 454 pyrosequencing and a multi-purpose read-filtering tool , 2011, BMC Research Notes.

[47]  Christus,et al.  A General Method Applicable to the Search for Similarities in the Amino Acid Sequence of Two Proteins , 2022 .

[48]  P. Hebert,et al.  Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species , 2003, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[49]  D. Baird,et al.  Environmental Barcoding: A Next-Generation Sequencing Approach for Biomonitoring Applications Using River Benthos , 2011, PloS one.

[50]  C. Dormann,et al.  Species abundance distributions and richness estimations in fungal metagenomics – lessons learned from community ecology , 2011, Molecular ecology.

[51]  Pedro Belda-Ferre,et al.  Amplification by PCR Artificially Reduces the Proportion of the Rare Biosphere in Microbial Communities , 2012, PloS one.

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

[53]  P. Dillon,et al.  The effect of excluding diatom taxa and reducing taxonomic resolution on multivariate analyses and stream bioassessment , 2009 .

[54]  Mehrdad Hajibabaei,et al.  Biomonitoring 2.0: a new paradigm in ecosystem assessment made possible by next‐generation DNA sequencing , 2012, Molecular ecology.

[55]  B. W. Sweeney,et al.  Can DNA barcodes of stream macroinvertebrates improve descriptions of community structure and water quality? , 2011, Journal of the North American Benthological Society.

[56]  Garrick Staples,et al.  TORQUE resource manager , 2006, SC.