Ion Torrent PGM as Tool for Fungal Community Analysis: A Case Study of Endophytes in Eucalyptus grandis Reveals High Taxonomic Diversity

The Kingdom Fungi adds substantially to the diversity of life, but due to their cryptic morphology and lifestyle, tremendous diversity, paucity of formally described specimens, and the difficulty in isolating environmental strains into culture, fungal communities are difficult to characterize. This is especially true for endophytic communities of fungi living in healthy plant tissue. The developments in next generation sequencing technologies are, however, starting to reveal the true extent of fungal diversity. One of the promising new technologies, namely semiconductor sequencing, has thus far not been used in fungal diversity assessments. In this study we sequenced the internal transcribed spacer 1 (ITS1) nuclear encoded ribosomal RNA of the endophytic community of the economically important tree, Eucalyptus grandis, from South Africa using the Ion Torrent Personal Genome Machine (PGM). We determined the impact of various analysis parameters on the interpretation of the results, namely different sequence quality parameter settings, different sequence similarity cutoffs for clustering and filtering of databases for removal of sequences with incomplete taxonomy. Sequence similarity cutoff values only had a marginal effect on the identified family numbers, whereas different sequence quality filters had a large effect (89 vs. 48 families between least and most stringent filters). Database filtering had a small, but statistically significant, effect on the assignment of sequences to reference sequences. The community was dominated by Ascomycota, and particularly by families in the Dothidiomycetes that harbor well-known plant pathogens. The study demonstrates that semiconductor sequencing is an ideal strategy for environmental sequencing of fungal communities. It also highlights some potential pitfalls in subsequent data analyses when using a technology with relatively short read lengths.

[1]  C. Mitchell,et al.  Release of invasive plants from fungal and viral pathogens , 2003, Nature.

[2]  K. O’Donnell,et al.  Two divergent intragenomic rDNA ITS2 types within a monophyletic lineage of the fungus Fusarium are nonorthologous. , 1997, Molecular phylogenetics and evolution.

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

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

[5]  Michael J. Wingfield,et al.  Botryosphaeria dothidea endophytic in Eucalyptus grandis and Eucalyptus nitens in South Africa , 1996 .

[6]  A. V. Kachalkin,et al.  Endophytic yeast fungi in plant storage tissues , 2010, Biology Bulletin.

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

[8]  Jens Stoye,et al.  Bacterial Community Shift in Treated Periodontitis Patients Revealed by Ion Torrent 16S rRNA Gene Amplicon Sequencing , 2012, PloS one.

[9]  R. Henrik Nilsson,et al.  Progress in molecular and morphological taxon discovery in Fungi and options for formal classification of environmental sequences , 2011 .

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

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

[12]  Jaysheel D. Bhavsar,et al.  Metagenomics: Read Length Matters , 2008, Applied and Environmental Microbiology.

[13]  John C. Wooley,et al.  Metagenomics: Facts and Artifacts, and Computational Challenges , 2010, Journal of Computer Science and Technology.

[14]  C. Rosa,et al.  Biodiversity and ecophysiology of yeasts , 2006 .

[15]  John C. Wooley,et al.  A Primer on Metagenomics , 2010, PLoS Comput. Biol..

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

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

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

[19]  K. Jones,et al.  Seasonally dynamic fungal communities in the Quercus macrocarpa phyllosphere differ between urban and nonurban environments. , 2010, The New phytologist.

[20]  M. Wingfield,et al.  Co-occurring species of Teratosphaeria on Eucalyptus , 2009, Persoonia.

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

[22]  Andy F. S. Taylor,et al.  The UNITE database for molecular identification of fungi--recent updates and future perspectives. , 2010, The New phytologist.

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

[24]  Marcus J. Claesson,et al.  Comparison of two next-generation sequencing technologies for resolving highly complex microbiota composition using tandem variable 16S rRNA gene regions , 2010, Nucleic acids research.

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

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

[27]  Anthony O'Donnell,et al.  Microbial 16S rRNA Ion Tag and community metagenome sequencing using the Ion Torrent (PGM) Platform. , 2012, Journal of microbiological methods.

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

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

[30]  K. Katoh,et al.  MAFFT version 5: improvement in accuracy of multiple sequence alignment , 2005, Nucleic acids research.

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

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

[33]  O. Petrini,et al.  A comparative study of fungal endophytes in leaves, xylem and bark of Eucalyptus in Australia and England , 2011 .

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

[35]  Tor Carlsen,et al.  Employing 454 amplicon pyrosequencing to reveal intragenomic divergence in the internal transcribed spacer rDNA region in fungi , 2013, Ecology and evolution.

[36]  Robert Samson,et al.  Indoor fungal composition is geographically patterned and more diverse in temperate zones than in the tropics , 2010, Proceedings of the National Academy of Sciences.

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

[38]  G. B. Golding,et al.  Are similarity- or phylogeny-based methods more appropriate for classifying internal transcribed spacer (ITS) metagenomic amplicons? , 2011, The New phytologist.

[39]  T. Burgess,et al.  A diverse assemblage of Botryosphaeriaceae infect Eucalyptus in native and non-native environments , 2009 .

[40]  Bernard P. Puc,et al.  An integrated semiconductor device enabling non-optical genome sequencing , 2011, Nature.

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

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

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

[44]  Pedro W. Crous,et al.  Phylogenetic reassessment of Mycosphaerella spp. and their anamorphs occurring on Eucalyptus. II. , 2006, Studies in mycology.

[45]  H. Swerdlow,et al.  A tale of three next generation sequencing platforms: comparison of Ion Torrent, Pacific Biosciences and Illumina MiSeq sequencers , 2012, BMC Genomics.

[46]  L. Øvreås,et al.  Microbial diversity and function in soil: from genes to ecosystems. , 2002, Current opinion in microbiology.

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

[48]  T. Glenn.  Field guide to next‐generation DNA sequencers , 2011, Molecular ecology resources.

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

[50]  D. Hawksworth.  The magnitude of fungal diversity: the 1.5 million species estimate revisited * * Paper presented at , 2001 .

[51]  T. Burgess,et al.  Mycosphaerella and Teratosphaeria diseases of Eucalyptus; easily confused and with serious consequences , 2011, Fungal Diversity.