Tissue storage and primer selection influence pyrosequencing‐based inferences of diversity and community composition of endolichenic and endophytic fungi

Next‐generation sequencing technologies have provided unprecedented insights into fungal diversity and ecology. However, intrinsic biases and insufficient quality control in next‐generation methods can lead to difficult‐to‐detect errors in estimating fungal community richness, distributions and composition. The aim of this study was to examine how tissue storage prior to DNA extraction, primer design and various quality‐control approaches commonly used in 454 amplicon pyrosequencing might influence ecological inferences in studies of endophytic and endolichenic fungi. We first contrast 454 data sets generated contemporaneously from subsets of the same plant and lichen tissues that were stored in CTAB buffer, dried in silica gel or freshly frozen prior to DNA extraction. We show that storage in silica gel markedly limits the recovery of sequence data and yields a small fraction of the diversity observed by the other two methods. Using lichen mycobiont sequences as internal positive controls, we next show that despite careful filtering of raw reads and utilization of current best‐practice OTU clustering methods, homopolymer errors in sequences representing rare taxa artificially increased estimates of richness c. 15‐fold in a model data set. Third, we show that inferences regarding endolichenic diversity can be improved using a novel primer that reduces amplification of the mycobiont. Together, our results provide a rationale for selecting tissue treatment regimes prior to DNA extraction, demonstrate the efficacy of reducing mycobiont amplification in studies of the fungal microbiomes of lichen thalli and highlight the difficulties in differentiating true information about fungal biodiversity from methodological artefacts.

[1]  C. Quince,et al.  Accurate determination of microbial diversity from 454 pyrosequencing data , 2009, Nature Methods.

[2]  P. Coley,et al.  Are tropical fungal endophytes hyperdiverse , 2000 .

[3]  Jolanta Miadlikowska,et al.  Host and geographic structure of endophytic and endolichenic fungi at a continental scale. , 2012, American journal of botany.

[4]  C. Kubicek,et al.  Identity, Diversity, and Molecular Phylogeny of the Endophytic Mycobiota in the Roots of Rare Wild Rice (Oryza granulate) from a Nature Reserve in Yunnan, China , 2009, Applied and Environmental Microbiology.

[5]  H. Friberg,et al.  New primers to amplify the fungal ITS2 region--evaluation by 454-sequencing of artificial and natural communities. , 2012, FEMS microbiology ecology.

[6]  F. Lutzoni,et al.  Phylogenetic relationships, host affinity, and geographic structure of boreal and arctic endophytes from three major plant lineages. , 2007, Molecular phylogenetics and evolution.

[7]  David L. Hawksworth,et al.  The fungal dimension of biodiversity: magnitude, significance, and conservation , 1991 .

[8]  M. Chase,et al.  Silica gel: An ideal material for field preservation of leaf samples for DNA studies , 1991 .

[9]  P. Keim,et al.  Specific detection of bacillus anthracis using a TaqMan mismatch amplification mutation assay. , 2005, BioTechniques.

[10]  M. Berbee,et al.  Culturing and direct DNA extraction find different fungi from the same ericoid mycorrhizal roots. , 2003, The New phytologist.

[11]  R. O’Hara,et al.  Host Genotype Shapes the Foliar Fungal Microbiome of Balsam Poplar (Populus balsamifera) , 2013, PloS one.

[12]  Ignazio Carbone,et al.  SNAP: workbench management tool for evolutionary population genetic analysis , 2005, Bioinform..

[13]  T. Bruns,et al.  Quantifying microbial communities with 454 pyrosequencing: does read abundance count? , 2010, Molecular ecology.

[14]  Scott T. Bates,et al.  A preliminary survey of lichen associated eukaryotes using pyrosequencing , 2011, The Lichenologist.

[15]  A. Arnold,et al.  Understanding the diversity of foliar endophytic fungi: progress, challenges, and frontiers , 2007 .

[16]  M. Hart,et al.  Differential effect of sample preservation methods on plant and arbuscular mycorrhizal fungal DNA. , 2010, Journal of microbiological methods.

[17]  A. Glenn,et al.  Exploring the evolutionary ecology of fungal endophytes in agricultural systems: using functional traits to reveal mechanisms in community processes , 2010, Evolutionary applications.

[18]  A. Arnold,et al.  Fungal endophytes: diversity and functional roles. , 2009, The New phytologist.

[19]  P. Green,et al.  Base-calling of automated sequencer traces using phred. I. Accuracy assessment. , 1998, Genome research.

[20]  Nathaniel L. Raizen Fungal endophyte diversity in foliage of native and cultivated Rhododendron species determined by culturing, ITS sequencing, and pyrosequencing , 2013 .

[21]  Rachel E Gallery,et al.  Diversity and evolutionary origins of fungi associated with seeds of a neotropical pioneer tree: a case study for analysing fungal environmental samples. , 2009, Mycological research.

[22]  T. Chujo,et al.  Fungal endophytes of grasses. , 2012, Current opinion in plant biology.

[23]  P. Coley,et al.  Culturing and direct PCR suggest prevalent host generalism among diverse fungal endophytes of tropical forest grasses , 2011, Mycologia.

[24]  K. R. Clarke,et al.  Non‐parametric multivariate analyses of changes in community structure , 1993 .

[25]  Susan M. Huse,et al.  Exploring Microbial Diversity and Taxonomy Using SSU rRNA Hypervariable Tag Sequencing , 2008, PLoS genetics.

[26]  K. Rodrigues The foliar fungal endophytes of the Amazonian palm Euterpe oleracea , 1994 .

[27]  G. May,et al.  Fungal-Fungal Associations Affect the Assembly of Endophyte Communities in Maize (Zea mays) , 2009, Microbial Ecology.

[28]  M. Sogin,et al.  Correction: Exploring Microbial Diversity and Taxonomy Using SSU rRNA Hypervariable Tag Sequencing , 2008, PLoS Genetics.

[29]  P. Fisher Survival and spread of the endophyte Stagonospora pteridiicola in Pteridium aquilinum, other ferns and some flowering plants. , 1996, The New phytologist.

[30]  A. Solow,et al.  Measuring biological diversity , 2006, Environmental and Ecological Statistics.

[31]  Douglas Ladd,et al.  Lichens of North America , 2002, Brittonia.

[32]  N. Pace,et al.  Differential amplification of rRNA genes by polymerase chain reaction , 1992, Applied and environmental microbiology.

[33]  Jolanta Miadlikowska,et al.  Community Analysis Reveals Close Affinities Between Endophytic and Endolichenic Fungi in Mosses and Lichens , 2010, Microbial Ecology.

[34]  Patrick D. Schloss,et al.  Reducing the Effects of PCR Amplification and Sequencing Artifacts on 16S rRNA-Based Studies , 2011, PloS one.

[35]  Kessy Abarenkov,et al.  Fungal community analysis by high-throughput sequencing of amplified markers – a user's guide , 2013, The New phytologist.

[36]  A. J. Shaw,et al.  Biogeographic and phylogenetic patterns in diversity of liverwort-associated endophytes. , 2008, American journal of botany.

[37]  Nicholas B. Simpson,et al.  Diverse Helotiales associated with the roots of three species of Arctic Ericaceae provide no evidence for host specificity. , 2011, The New phytologist.

[38]  Rytas Vilgalys,et al.  Fungal Community Analysis by Large-Scale Sequencing of Environmental Samples , 2005, Applied and Environmental Microbiology.

[39]  Tom O. Delmont,et al.  Accessing the Soil Metagenome for Studies of Microbial Diversity , 2010, Applied and Environmental Microbiology.

[40]  S. Jarman,et al.  Blocking primers to enhance PCR amplification of rare sequences in mixed samples – a case study on prey DNA in Antarctic krill stomachs , 2008, Frontiers in Zoology.

[41]  Satoshi Yamamoto,et al.  High-Coverage ITS Primers for the DNA-Based Identification of Ascomycetes and Basidiomycetes in Environmental Samples , 2012, PloS one.

[42]  P Green,et al.  Base-calling of automated sequencer traces using phred. II. Error probabilities. , 1998, Genome research.

[43]  K. Jones,et al.  Massively parallel 454 sequencing indicates hyperdiverse fungal communities in temperate Quercus macrocarpa phyllosphere. , 2009, The New phytologist.

[44]  James Long,et al.  TOPO TA is A-OK: a test of phylogenetic bias in fungal environmental clone library construction. , 2007, Environmental microbiology.

[45]  R. Fisher,et al.  The Relation Between the Number of Species and the Number of Individuals in a Random Sample of an Animal Population , 1943 .

[46]  Robert C. Edgar,et al.  MUSCLE: multiple sequence alignment with high accuracy and high throughput. , 2004, Nucleic acids research.

[47]  F. Lutzoni,et al.  Diversity and host range of foliar fungal endophytes: are tropical leaves biodiversity hotspots? , 2007, Ecology.

[48]  Rytas Vilgalys,et al.  Diversity and phylogenetic affinities of foliar fungal endophytes in loblolly pine inferred by culturing and environmental PCR , 2007, Mycologia.

[49]  T. Bruns,et al.  ITS primers with enhanced specificity for basidiomycetes ‐ application to the identification of mycorrhizae and rusts , 1993, Molecular ecology.

[50]  Jolanta Miadlikowska,et al.  A phylogenetic estimation of trophic transition networks for ascomycetous fungi: are lichens cradles of symbiotrophic fungal diversification? , 2009, Systematic biology.

[51]  Erik Kristiansson,et al.  An outlook on the fungal internal transcribed spacer sequences in GenBank and the introduction of a web-based tool for the exploration of fungal diversity. , 2009, The New phytologist.

[52]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[53]  Susan M. Huse,et al.  Accuracy and quality of massively parallel DNA pyrosequencing , 2007, Genome Biology.

[54]  Deepak Bhatnagar,et al.  Promising nucleic acid analogs and mimics: characteristic features and applications of PNA, LNA, and morpholino , 2006, Applied Microbiology and Biotechnology.

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

[56]  V. Kunin,et al.  Wrinkles in the rare biosphere: pyrosequencing errors can lead to artificial inflation of diversity estimates. , 2009, Environmental microbiology.

[57]  Guo,et al.  Endophytic fungi associated with lichens in Baihua mountain of Beijing, China , 2007 .

[58]  Z. Nagy A hands-on overview of tissue preservation methods for molecular genetic analyses , 2010, Organisms Diversity & Evolution.

[59]  R. Henrik Nilsson,et al.  Approaching the taxonomic affiliation of unidentified sequences in public databases – an example from the mycorrhizal fungi , 2005, BMC Bioinformatics.

[60]  D. Maddison,et al.  Mesquite: a modular system for evolutionary analysis. Version 2.6 , 2009 .

[61]  Susan M. Huse,et al.  Ironing out the wrinkles in the rare biosphere through improved OTU clustering , 2010, Environmental microbiology.

[62]  J. Tibbits,et al.  A rapid method for tissue collection and high-throughput isolation of genomic DNA from mature trees , 2006, Plant Molecular Biology Reporter.

[63]  Martin Hartmann,et al.  Introducing mothur: Open-Source, Platform-Independent, Community-Supported Software for Describing and Comparing Microbial Communities , 2009, Applied and Environmental Microbiology.

[64]  Ignazio Carbone,et al.  Mobyle SNAP Workbench: a web-based analysis portal for population genetics and evolutionary genomics , 2014, Bioinform..

[65]  K. Hyde,et al.  Detection and taxonomic placement of endophytic fungi within frond tissues of Livistona chinensis based on rDNA sequences. , 2001, Molecular phylogenetics and evolution.

[66]  William G. Mckendree,et al.  ESPRIT: estimating species richness using large collections of 16S rRNA pyrosequences , 2009, Nucleic acids research.

[67]  Rob Knight,et al.  UCHIME improves sensitivity and speed of chimera detection , 2011, Bioinform..

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

[69]  Susan M. Huse,et al.  Microbial diversity in the deep sea and the underexplored “rare biosphere” , 2006, Proceedings of the National Academy of Sciences.

[70]  R. Halvorsen,et al.  Amplicon‐pyrosequencing‐based detection of compositional shifts in bryophyte‐associated fungal communities along an elevation gradient , 2013, Molecular ecology.

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

[72]  Peter M. Vitousek,et al.  Fungal endophyte communities reflect environmental structuring across a Hawaiian landscape , 2012, Proceedings of the National Academy of Sciences.

[73]  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.

[74]  M. B. Couger,et al.  Phylogenetic diversity and community structure of anaerobic gut fungi (phylum Neocallimastigomycota) in ruminant and non-ruminant herbivores , 2010, The ISME Journal.

[75]  T. White Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics , 1990 .

[76]  R. Henrik Nilsson,et al.  An open source software package for automated extraction of ITS1 and ITS2 from fungal ITS sequences for use in high-throughput community assays and molecular ecology , 2010 .

[77]  D. J. Lodge,et al.  ENDOPHYTIC FUNGI OF MANILKARA BIDENTATA LEAVES IN PUERTO RICO , 1996 .

[78]  I. Schmitt,et al.  Comparison of ITS1 and ITS2 rDNA in 454 sequencing of hyperdiverse fungal communities , 2013 .

[79]  O. Petrini,et al.  Fungal endophytes from the leaves and twigs of Quercus ilex L. from England, Majorca and Switzerland. , 1994, The New phytologist.

[80]  S. Chao,et al.  High throughput tissue preparation for large‐scale genotyping experiments , 2008, Molecular ecology resources.

[81]  J. A. Johnson,et al.  Isolation of endophytic fungi from eastern larch (Larix laricina) leaves from New Brunswick, Canada , 1995 .

[82]  F. Graf,et al.  FUNGAL ENDOPHYTES OF DRYAS OCTOPETALA FROM A HIGH ARCTIC POLAR SEMIDESERT AND FROM THE SWISS ALPS , 1995 .

[83]  Pierre Taberlet,et al.  ITS as an environmental DNA barcode for fungi: an in silico approach reveals potential PCR biases , 2010, BMC Microbiology.

[84]  E. Herre,et al.  Ecological implications of anti-pathogen effects of tropical fungal endophytes and mycorrhizae. , 2007, Ecology.

[85]  R. Hrdličková,et al.  An improved method of DNA isolation from plants collected in the field and conserved in saturated NaCl/CTAB solution. , 2000 .