Shotgun metagenomes and multiple primer pair-barcode combinations of amplicons reveal biases in metabarcoding analyses of fungi

Rapid development of high-throughput (HTS) molecular identification methods has revolutionized our knowledge about taxonomic diversity and ecology of fungi. However, PCR-based methods exhibit multiple technical shortcomings that may bias our understanding of the fungal kingdom. This study was initiated to quantify potential biases in fungal community ecology by comparing the relative performance of amplicon-free shotgun metagenomics and amplicons of nine primer pairs over seven nuclear ribosomal DNA (rDNA) regions often used in metabarcoding analyses. The internal transcribed spacer (ITS) bar codes ITS1 and ITS2 provided greater taxonomic and functional resolution and richness of operational taxonomic units (OTUs) at the 97% similarity threshold compared to barcodes located within the ribosomal small subunit (SSU) and large subunit (LSU) genes. All barcode-primer pair combinations provided consistent results in ranking taxonomic richness and recovering the importance of floristic variables in driving fungal community composition in soils of Papua New Guinea. The choice of forward primer explained up to 2.0% of the variation in OTU-level analysis of the ITS1 and ITS2 barcode data sets. Across the whole data set, barcode-primer pair combination explained 37.6–38.1% of the variation, which surpassed any environmental signal. Overall, the metagenomics data set recovered a similar taxonomic overview, but resulted in much lower fungal rDNA sequencing depth, inability to infer OTUs, and high

[1]  J. Koricheva,et al.  Tree diversity and species identity effects on soil fungi, protists and animals are context dependent , 2015, The ISME Journal.

[2]  J. Bengtsson-Palme,et al.  ITS1: a DNA barcode better than ITS2 in eukaryotes? , 2015, Molecular ecology resources.

[3]  R. Henrik Nilsson,et al.  A Comprehensive, Automatically Updated Fungal ITS Sequence Dataset for Reference-Based Chimera Control in Environmental Sequencing Efforts , 2015, Microbes and environments.

[4]  T. Kull,et al.  Temporal patterns of orchid mycorrhizal fungi in meadows and forests as revealed by 454 pyrosequencing. , 2015, The New phytologist.

[5]  N. Ishii,et al.  PCR primers for assessing community structure of aquatic fungi including Chytridiomycota and Cryptomycota , 2015 .

[6]  A. Jumpponen,et al.  Scraping the bottom of the barrel: are rare high throughput sequences artifacts? , 2015 .

[7]  S. Tringe,et al.  High-Throughput Metagenomic Technologies for Complex Microbial Community Analysis: Open and Closed Formats , 2015, mBio.

[8]  L. Tedersoo,et al.  Global biogeography of Alnus-associated Frankia actinobacteria. , 2014, The New phytologist.

[9]  R. Henrik Nilsson,et al.  Global diversity and geography of soil fungi , 2014, Science.

[10]  Falk Hildebrand,et al.  Correction: LotuS: an efficient and user-friendly OTU processing pipeline , 2014, Microbiome.

[11]  Francisco M. Cornejo-Castillo,et al.  Metagenomic 16S rDNA Illumina tags are a powerful alternative to amplicon sequencing to explore diversity and structure of microbial communities. , 2014, Environmental microbiology.

[12]  H. Jacquemyn,et al.  Impact of primer choice on characterization of orchid mycorrhizal communities using 454 pyrosequencing , 2014, Molecular ecology resources.

[13]  Bart Lievens,et al.  Comparison and Validation of Some ITS Primer Pairs Useful for Fungal Metabarcoding Studies , 2014, PloS one.

[14]  Jack W. McFarland,et al.  A first comprehensive census of fungi in soil reveals both hyperdiversity and fine-scale niche partitioning , 2014 .

[15]  M. Hejda,et al.  Comparison of commonly used primer sets for evaluating arbuscular mycorrhizal fungal communities: Is there a universal solution? , 2014 .

[16]  L. Tedersoo,et al.  Lineages of ectomycorrhizal fungi revisited: Foraging strategies and novel lineages revealed by sequences from belowground , 2013 .

[17]  Michael Weiss,et al.  Towards a unified paradigm for sequence‐based identification of fungi , 2013, Molecular ecology.

[18]  R. Henrik Nilsson,et al.  Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data , 2013 .

[19]  T. Mattes,et al.  Phylogenetic detection of novel Cryptomycota in an Iowa (United States) aquifer and from previously collected marine and freshwater targeted high-throughput sequencing sets. , 2013, Environmental microbiology.

[20]  Qing Yang,et al.  Ultra-deep sequencing enables high-fidelity recovery of biodiversity for bulk arthropod samples without PCR amplification , 2013, GigaScience.

[21]  P. Kirk,et al.  ITS1 versus ITS2 as DNA metabarcodes for fungi , 2013, Molecular ecology resources.

[22]  Pelin Yilmaz,et al.  The SILVA ribosomal RNA gene database project: improved data processing and web-based tools , 2012, Nucleic Acids Res..

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

[24]  Hélène Touzet,et al.  SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data , 2012, Bioinform..

[25]  Zhengwei Zhu,et al.  CD-HIT: accelerated for clustering the next-generation sequencing data , 2012, Bioinform..

[26]  C. Murat,et al.  Fungal Diversity Is Not Determined by Mineral and Chemical Differences in Serpentine Substrates , 2012, PloS one.

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

[28]  M. Hartmann,et al.  Megraft: a software package to graft ribosomal small subunit (16S/18S) fragments onto full-length sequences for accurate species richness and sequencing depth analysis in pyrosequencing-length metagenomes and similar environmental datasets. , 2012, Research in microbiology.

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

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

[31]  Anthony Levasseur,et al.  Multiple markers pyrosequencing reveals highly diverse and host-specific fungal communities on the mangrove trees Avicennia marina and Rhizophora stylosa. , 2012, FEMS microbiology ecology.

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

[33]  M. Blackwell.  The fungi: 1, 2, 3 ... 5.1 million species? , 2011, American journal of botany.

[34]  T. Fennell,et al.  Analyzing and minimizing PCR amplification bias in Illumina sequencing libraries , 2011, Genome Biology.

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

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

[37]  Wolfgang Maier,et al.  Current state and perspectives of fungal DNA barcoding and rapid identification procedures , 2010, Applied Microbiology and Biotechnology.

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

[39]  C. Schadt,et al.  Soil Microbial Community Responses to Multiple Experimental Climate Change Drivers , 2009, Applied and Environmental Microbiology.

[40]  Sergei L. Kosakovsky Pond,et al.  Windshield splatter analysis with the Galaxy metagenomic pipeline. , 2009, Genome research.

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

[42]  Erik Kristiansson,et al.  The ITS region as a target for characterization of fungal communities using emerging sequencing technologies. , 2009, FEMS microbiology letters.

[43]  J. Cavender-Bares,et al.  The merging of community ecology and phylogenetic biology. , 2009, Ecology letters.

[44]  L. Tedersoo,et al.  Strong host preference of ectomycorrhizal fungi in a Tasmanian wet sclerophyll forest as revealed by DNA barcoding and taxon-specific primers. , 2008, The New phytologist.

[45]  M. McCormick,et al.  Internal transcribed spacer primers and sequences for improved characterization of basidiomycetous orchid mycorrhizas. , 2008, The New phytologist.

[46]  J. Tiedje,et al.  Naïve Bayesian Classifier for Rapid Assignment of rRNA Sequences into the New Bacterial Taxonomy , 2007, Applied and Environmental Microbiology.

[47]  Paul T Rygiewicz,et al.  Fungal-specific PCR primers developed for analysis of the ITS region of environmental DNA extracts , 2005, BMC Microbiology.

[48]  Christopher W. Schadt,et al.  Seasonal Dynamics of Previously Unknown Fungal Lineages in Tundra Soils , 2003, Science.

[49]  Takahiro Kanagawa,et al.  Bias and artifacts in multitemplate polymerase chain reactions (PCR). , 2003, Journal of bioscience and bioengineering.

[50]  Pedro Jordano,et al.  GEOGRAPHIC PATTERNS IN PLANT–POLLINATOR MUTUALISTIC NETWORKS , 2002 .

[51]  Piero Carninci,et al.  Balanced-size and long-size cloning of full-length, cap-trapped cDNAs into vectors of the novel lambda-FLC family allows enhanced gene discovery rate and functional analysis. , 2001, Genomics.

[52]  D. Bhattacharya,et al.  Widespread occurrence of spliceosomal introns in the rDNA genes of ascomycetes. , 2000, Molecular biology and evolution.

[53]  R. J. Hartin,et al.  PCR Primers That Amplify Fungal rRNA Genes from Environmental Samples , 2000, Applied and Environmental Microbiology.

[54]  J. Palmer,et al.  Investigating Deep Phylogenetic Relationships among Cyanobacteria and Plastids by Small Subunit rRNA Sequence Analysis 1 , 1999, The Journal of eukaryotic microbiology.

[55]  S. Sanche,et al.  Rapid Identification of Fungi by Using the ITS2 Genetic Region and an Automated Fluorescent Capillary Electrophoresis System , 1999, Journal of Clinical Microbiology.

[56]  R. Vilgalys,et al.  Phylogenetic relationships among coprinoid taxa and allies based on data from restriction site mapping of nuclear rDNA , 1994 .

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

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