Illumina metabarcoding of a soil fungal community

Abstract Next generation metabarcoding is becoming an indispensable tool in fungal community ecology. Here we tested Illumina metabarcoding, a method that generates shorter reads but achieves deeper sequencing than 454 metabarcoding approaches. We found that paired-end Illumina MiSeq data cover the full ITS1 in many fungal lineages and are suitable for environmental fungal community assessment. There was substantial read loss during data cleanup (78.6%), which, however, did not impede the analyses, because of the large number of initial sequences (over 4Mio). We observed a high stochasticity in individual PCR reactions. Comparing three repeated sets of PCRs products showed that 58.5% of the total fungal operational taxonomic units (OTUs) found were not recovered by any single set of PCR reactions. Similarly, comparing three annealing temperatures showed that 63.6% of all fungal OTUs were not recovered using any single annealing temperature. These findings suggest that sampling of soil fungal communities is more exhaustive, if we combine repeated PCR products, and PCR products generated at various annealing temperatures. To analyze the above issues we sampled 16 soil cores along a 270 cm transect in a meadow. In total we recovered 3320 fungal OTUs (based on a 95% similarity threshold). Distance decay analysis indicated that community similarity decreased slightly with geographical distance.

[1]  Thomas J. White,et al.  PCR protocols: a guide to methods and applications. , 1990 .

[2]  G. Kovács,et al.  Variation in the nrDNA ITS sequences of some powdery mildew species: do routine molecular identification procedures hide valuable information? , 2011, European Journal of Plant Pathology.

[3]  J. Larsen,et al.  Molecular characterization of root-associated fungal communities in relation to health status of Pisum sativum using barcoded pyrosequencing , 2012, Plant and Soil.

[4]  P. White,et al.  The distance decay of similarity in biogeography and ecology , 1999 .

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

[6]  William A. Walters,et al.  Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms , 2012, The ISME Journal.

[7]  P. Meinicke,et al.  Fungal soil communities in a young transgenic poplar plantation form a rich reservoir for fungal root communities , 2012, Ecology and evolution.

[8]  Jean M. Macklaim,et al.  Microbiome Profiling by Illumina Sequencing of Combinatorial Sequence-Tagged PCR Products , 2010, PLoS ONE.

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

[10]  F. Martin,et al.  454 Pyrosequencing analyses of forest soils reveal an unexpectedly high fungal diversity. , 2009, The New phytologist.

[11]  S. Sørensen,et al.  454‐sequencing reveals stochastic local reassembly and high disturbance tolerance within arbuscular mycorrhizal fungal communities , 2012 .

[12]  A. Jumpponen,et al.  Massively parallel 454‐sequencing of fungal communities in Quercus spp. ectomycorrhizas indicates seasonal dynamics in urban and rural sites , 2010, Molecular ecology.

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

[14]  D. Lindner,et al.  Intragenomic variation in the ITS rDNA region obscures phylogenetic relationships and inflates estimates of operational taxonomic units in genus Laetiporus , 2011, Mycologia.

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

[16]  H. Seliger,et al.  PCR protocols — A guide to methods and applications , 1990 .

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

[18]  Rita Sipos,et al.  Effect of primer mismatch, annealing temperature and PCR cycle number on 16S rRNA gene-targetting bacterial community analysis. , 2007, FEMS microbiology ecology.

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

[20]  N. Kyrpides,et al.  Direct Comparisons of Illumina vs. Roche 454 Sequencing Technologies on the Same Microbial Community DNA Sample , 2012, PloS one.

[21]  Marti J. Anderson,et al.  Much ado about nothings: using zero similarity points in distance-decay curves. , 2011, Ecology.

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

[23]  Daniel G. Brown,et al.  PANDAseq: paired-end assembler for illumina sequences , 2012, BMC Bioinformatics.

[24]  M. Desprez-Loustau,et al.  Spatial variability of phyllosphere fungal assemblages: genetic distance predominates over geographic distance in a European beech stand (Fagus sylvatica) , 2012 .

[25]  C. Quince,et al.  Sample richness and genetic diversity as drivers of chimera formation in nSSU metagenetic analyses , 2012, Nucleic acids research.

[26]  M. Zobel,et al.  Large-scale parallel 454 sequencing reveals host ecological group specificity of arbuscular mycorrhizal fungi in a boreonemoral forest. , 2009, The New phytologist.

[27]  Peter J. Schaap,et al.  Molecular characterization of the , 1997 .

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

[29]  C. Dytham,et al.  Distinct seasonal assemblages of arbuscular mycorrhizal fungi revealed by massively parallel pyrosequencing. , 2011, The New phytologist.

[30]  L. Tedersoo,et al.  Fine scale distribution of ectomycorrhizal fungi and roots across substrate layers including coarse woody debris in a mixed forest. , 2003, The New phytologist.

[31]  L. H. Liow,et al.  Ecology and evolution. , 2014, Ecology and evolution.

[32]  S. Schuster,et al.  Integrative analysis of environmental sequences using MEGAN4. , 2011, Genome research.

[33]  Vanja Klepac-Ceraj,et al.  PCR-Induced Sequence Artifacts and Bias: Insights from Comparison of Two 16S rRNA Clone Libraries Constructed from the Same Sample , 2005, Applied and Environmental Microbiology.

[34]  R. H. Nilsson,et al.  Soil fungal community structure along a soil health gradient in pea fields examined using deep amplicon sequencing , 2012 .

[35]  John L. Spouge,et al.  Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi , 2012, Proceedings of the National Academy of Sciences.

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

[37]  K. I. Ugland,et al.  Changes in the root‐associated fungal communities along a primary succession gradient analysed by 454 pyrosequencing , 2012, Molecular ecology.

[38]  M. Weiß,et al.  Intragenomic variation of fungal ribosomal genes is higher than previously thought. , 2008, Molecular biology and evolution.

[39]  M. Fukui,et al.  Optimization of Annealing Temperature To Reduce Bias Caused by a Primer Mismatch in Multitemplate PCR , 2001, Applied and Environmental Microbiology.

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