Comparison and Validation of Some ITS Primer Pairs Useful for Fungal Metabarcoding Studies
Identifieur interne : 000228 ( Pmc/Corpus ); précédent : 000227; suivant : 000229Comparison and Validation of Some ITS Primer Pairs Useful for Fungal Metabarcoding Studies
Auteurs : Michiel Op De Beeck ; Bart Lievens ; Pieter Busschaert ; Stéphan Declerck ; Jaco Vangronsveld ; Jan V. ColpaertSource :
- PLoS ONE [ 1932-6203 ] ; 2014.
Abstract
Current metabarcoding studies aiming to characterize microbial communities generally rely on the amplification and sequencing of relatively short DNA regions. For fungi, the internal transcribed spacer (ITS) region in the ribosomal RNA (rRNA) operon has been accepted as the formal fungal barcode. Despite an increasing number of fungal metabarcoding studies, the amplification efficiency of primers is generally not tested prior to their application in metabarcoding studies. Some of the challenges that metabarcoding primers should overcome efficiently are the amplification of target DNA strands in samples rich in non-target DNA and environmental pollutants, such as humic acids, that may have been co-extracted with DNA. In the current study, three selected primer pairs were tested for their suitability as fungal metabarcoding primers. The selected primer pairs include two primer pairs that have been frequently used in fungal metabarcoding studies (ITS1F/ITS2 and ITS3/ITS4) and a primer pair (ITS86F/ITS4) that has been shown to efficiently amplify the ITS2 region of a broad range of fungal taxa in environmental soil samples. The selected primer pairs were evaluated in a 454 amplicon pyrosequencing experiment, real-time PCR (qPCR) experiments and
Url:
DOI: 10.1371/journal.pone.0097629
PubMed: 24933453
PubMed Central: 4059633
Links to Exploration step
PMC:4059633Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Comparison and Validation of Some ITS Primer Pairs Useful for Fungal Metabarcoding Studies</title><author><name sortKey="Op De Beeck, Michiel" sort="Op De Beeck, Michiel" uniqKey="Op De Beeck M" first="Michiel" last="Op De Beeck">Michiel Op De Beeck</name><affiliation><nlm:aff id="aff1"><addr-line>Centre for Environmental Sciences, Hasselt University, Hasselt, Limburg, Belgium</addr-line></nlm:aff></affiliation></author><author><name sortKey="Lievens, Bart" sort="Lievens, Bart" uniqKey="Lievens B" first="Bart" last="Lievens">Bart Lievens</name><affiliation><nlm:aff id="aff2"><addr-line>Department of Microbial and Molecular Systems (M2S), Catholic University of Leuven, Antwerp, Belgium</addr-line></nlm:aff></affiliation></author><author><name sortKey="Busschaert, Pieter" sort="Busschaert, Pieter" uniqKey="Busschaert P" first="Pieter" last="Busschaert">Pieter Busschaert</name><affiliation><nlm:aff id="aff2"><addr-line>Department of Microbial and Molecular Systems (M2S), Catholic University of Leuven, Antwerp, Belgium</addr-line></nlm:aff></affiliation></author><author><name sortKey="Declerck, Stephan" sort="Declerck, Stephan" uniqKey="Declerck S" first="Stéphan" last="Declerck">Stéphan Declerck</name><affiliation><nlm:aff id="aff3"><addr-line>Earth & Life Institute, Catholic University of Louvain, Louvain-la-Neuve, Walloon Brabant, Belgium</addr-line></nlm:aff></affiliation></author><author><name sortKey="Vangronsveld, Jaco" sort="Vangronsveld, Jaco" uniqKey="Vangronsveld J" first="Jaco" last="Vangronsveld">Jaco Vangronsveld</name><affiliation><nlm:aff id="aff1"><addr-line>Centre for Environmental Sciences, Hasselt University, Hasselt, Limburg, Belgium</addr-line></nlm:aff></affiliation></author><author><name sortKey="Colpaert, Jan V" sort="Colpaert, Jan V" uniqKey="Colpaert J" first="Jan V." last="Colpaert">Jan V. Colpaert</name><affiliation><nlm:aff id="aff1"><addr-line>Centre for Environmental Sciences, Hasselt University, Hasselt, Limburg, Belgium</addr-line></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">24933453</idno><idno type="pmc">4059633</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4059633</idno><idno type="RBID">PMC:4059633</idno><idno type="doi">10.1371/journal.pone.0097629</idno><date when="2014">2014</date><idno type="wicri:Area/Pmc/Corpus">000228</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000228</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Comparison and Validation of Some ITS Primer Pairs Useful for Fungal Metabarcoding Studies</title><author><name sortKey="Op De Beeck, Michiel" sort="Op De Beeck, Michiel" uniqKey="Op De Beeck M" first="Michiel" last="Op De Beeck">Michiel Op De Beeck</name><affiliation><nlm:aff id="aff1"><addr-line>Centre for Environmental Sciences, Hasselt University, Hasselt, Limburg, Belgium</addr-line></nlm:aff></affiliation></author><author><name sortKey="Lievens, Bart" sort="Lievens, Bart" uniqKey="Lievens B" first="Bart" last="Lievens">Bart Lievens</name><affiliation><nlm:aff id="aff2"><addr-line>Department of Microbial and Molecular Systems (M2S), Catholic University of Leuven, Antwerp, Belgium</addr-line></nlm:aff></affiliation></author><author><name sortKey="Busschaert, Pieter" sort="Busschaert, Pieter" uniqKey="Busschaert P" first="Pieter" last="Busschaert">Pieter Busschaert</name><affiliation><nlm:aff id="aff2"><addr-line>Department of Microbial and Molecular Systems (M2S), Catholic University of Leuven, Antwerp, Belgium</addr-line></nlm:aff></affiliation></author><author><name sortKey="Declerck, Stephan" sort="Declerck, Stephan" uniqKey="Declerck S" first="Stéphan" last="Declerck">Stéphan Declerck</name><affiliation><nlm:aff id="aff3"><addr-line>Earth & Life Institute, Catholic University of Louvain, Louvain-la-Neuve, Walloon Brabant, Belgium</addr-line></nlm:aff></affiliation></author><author><name sortKey="Vangronsveld, Jaco" sort="Vangronsveld, Jaco" uniqKey="Vangronsveld J" first="Jaco" last="Vangronsveld">Jaco Vangronsveld</name><affiliation><nlm:aff id="aff1"><addr-line>Centre for Environmental Sciences, Hasselt University, Hasselt, Limburg, Belgium</addr-line></nlm:aff></affiliation></author><author><name sortKey="Colpaert, Jan V" sort="Colpaert, Jan V" uniqKey="Colpaert J" first="Jan V." last="Colpaert">Jan V. Colpaert</name><affiliation><nlm:aff id="aff1"><addr-line>Centre for Environmental Sciences, Hasselt University, Hasselt, Limburg, Belgium</addr-line></nlm:aff></affiliation></author></analytic><series><title level="j">PLoS ONE</title><idno type="eISSN">1932-6203</idno><imprint><date when="2014">2014</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>Current metabarcoding studies aiming to characterize microbial communities generally rely on the amplification and sequencing of relatively short DNA regions. For fungi, the internal transcribed spacer (ITS) region in the ribosomal RNA (rRNA) operon has been accepted as the formal fungal barcode. Despite an increasing number of fungal metabarcoding studies, the amplification efficiency of primers is generally not tested prior to their application in metabarcoding studies. Some of the challenges that metabarcoding primers should overcome efficiently are the amplification of target DNA strands in samples rich in non-target DNA and environmental pollutants, such as humic acids, that may have been co-extracted with DNA. In the current study, three selected primer pairs were tested for their suitability as fungal metabarcoding primers. The selected primer pairs include two primer pairs that have been frequently used in fungal metabarcoding studies (ITS1F/ITS2 and ITS3/ITS4) and a primer pair (ITS86F/ITS4) that has been shown to efficiently amplify the ITS2 region of a broad range of fungal taxa in environmental soil samples. The selected primer pairs were evaluated in a 454 amplicon pyrosequencing experiment, real-time PCR (qPCR) experiments and <italic>in silico</italic> analyses. Results indicate that experimental evaluation of primers provides valuable information that could aid in the selection of suitable primers for fungal metabarcoding studies. Furthermore, we show that the ITS86F/ITS4 primer pair outperforms other primer pairs tested in terms of <italic>in silico</italic> primer efficiency, PCR efficiency, coverage, number of reads and number of species-level operational taxonomic units (OTUs) obtained. These traits push the ITS86F/ITS4 primer pair forward as highly suitable for studying fungal diversity and community structures using DNA metabarcoding.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Stahl, Da" uniqKey="Stahl D">DA Stahl</name></author><author><name sortKey="Lane, Dj" uniqKey="Lane D">DJ Lane</name></author><author><name sortKey="Olsen, Gj" uniqKey="Olsen G">GJ Olsen</name></author><author><name sortKey="Pace, Nr" uniqKey="Pace N">NR Pace</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hugenholtz, P" uniqKey="Hugenholtz P">P Hugenholtz</name></author><author><name sortKey="Pace, Nr" uniqKey="Pace N">NR Pace</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Margulies, M" uniqKey="Margulies M">M Margulies</name></author><author><name sortKey="Egholm, M" uniqKey="Egholm M">M Egholm</name></author><author><name sortKey="Altman, We" uniqKey="Altman W">WE Altman</name></author><author><name sortKey="Attyia, S" uniqKey="Attyia S">S Attyia</name></author><author><name sortKey="Bader, Js" uniqKey="Bader J">JS Bader</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sogin, Ml" uniqKey="Sogin M">ML Sogin</name></author><author><name sortKey="Morrison, Hg" uniqKey="Morrison H">HG Morrison</name></author><author><name sortKey="Huber, Ja" uniqKey="Huber J">JA Huber</name></author><author><name sortKey="Welch, Dw" uniqKey="Welch D">DW Welch</name></author><author><name sortKey="Huse, Sm" uniqKey="Huse S">SM Huse</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Buee, M" uniqKey="Buee M">M Buée</name></author><author><name sortKey="Reich, M" uniqKey="Reich M">M Reich</name></author><author><name sortKey="Murat, C" uniqKey="Murat C">C Murat</name></author><author><name sortKey="Morin, E" uniqKey="Morin E">E Morin</name></author><author><name sortKey="Nilsson, Rh" uniqKey="Nilsson R">RH Nilsson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Taberlet, P" uniqKey="Taberlet P">P Taberlet</name></author><author><name sortKey="Coissac, E" uniqKey="Coissac E">E Coissac</name></author><author><name sortKey="Pompanon, F" uniqKey="Pompanon F">F Pompanon</name></author><author><name sortKey="Brochmann, C" uniqKey="Brochmann C">C Brochmann</name></author><author><name sortKey="Willerslev, E" uniqKey="Willerslev E">E Willerslev</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hebert, Pdn" uniqKey="Hebert P">PDN Hebert</name></author><author><name sortKey="Cywinska, A" uniqKey="Cywinska A">A Cywinska</name></author><author><name sortKey="Ball, Sl" uniqKey="Ball S">SL Ball</name></author><author><name sortKey="Dewaard, Jr" uniqKey="Dewaard J">JR deWaard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Juste, A" uniqKey="Juste A">A Justé</name></author><author><name sortKey="Thomma, Bphj" uniqKey="Thomma B">BPHJ Thomma</name></author><author><name sortKey="Lievens, B" uniqKey="Lievens B">B Lievens</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Seifert, Ka" uniqKey="Seifert K">KA Seifert</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Begerow, D" uniqKey="Begerow D">D Begerow</name></author><author><name sortKey="Nilsson, H" uniqKey="Nilsson H">H Nilsson</name></author><author><name sortKey="Unterseher, M" uniqKey="Unterseher M">M Unterseher</name></author><author><name sortKey="Maier, W" uniqKey="Maier W">W Maier</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schoch, Cl" uniqKey="Schoch C">CL Schoch</name></author><author><name sortKey="Seifert, Ka" uniqKey="Seifert K">KA Seifert</name></author><author><name sortKey="Huhndorf, S" uniqKey="Huhndorf S">S Huhndorf</name></author><author><name sortKey="Robert, V" uniqKey="Robert V">V Robert</name></author><author><name sortKey="Spouge, Jl" uniqKey="Spouge J">JL Spouge</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Porter, Tm" uniqKey="Porter T">TM Porter</name></author><author><name sortKey="Golding, Gb" uniqKey="Golding G">GB Golding</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Haas, Bj" uniqKey="Haas B">BJ Haas</name></author><author><name sortKey="Gevers, D" uniqKey="Gevers D">D Gevers</name></author><author><name sortKey="Earl, Am" uniqKey="Earl A">AM Earl</name></author><author><name sortKey="Feldgarden, M" uniqKey="Feldgarden M">M Feldgarden</name></author><author><name sortKey="Ward, Dv" uniqKey="Ward D">DV Ward</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kosch, Ta" uniqKey="Kosch T">TA Kosch</name></author><author><name sortKey="Summers, K" uniqKey="Summers K">K Summers</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Polz, Mf" uniqKey="Polz M">MF Polz</name></author><author><name sortKey="Cavanaugh, Cm" uniqKey="Cavanaugh C">CM Cavanaugh</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jumpponen, A" uniqKey="Jumpponen A">A Jumpponen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Engelbrektson, A" uniqKey="Engelbrektson A">A Engelbrektson</name></author><author><name sortKey="Kunin, V" uniqKey="Kunin V">V Kunin</name></author><author><name sortKey="Wrighton, Kc" uniqKey="Wrighton K">KC Wrighton</name></author><author><name sortKey="Zvenigorodsky, N" uniqKey="Zvenigorodsky N">N Zvenigorodsky</name></author><author><name sortKey="Chen, F" uniqKey="Chen F">F Chen</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Gardes, M" uniqKey="Gardes M">M Gardes</name></author><author><name sortKey="Bruns, Td" uniqKey="Bruns T">TD Bruns</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jumpponen, A" uniqKey="Jumpponen A">A Jumpponen</name></author><author><name sortKey="Jones, Kl" uniqKey="Jones K">KL Jones</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Amend, As" uniqKey="Amend A">AS Amend</name></author><author><name sortKey="Seifert, Ka" uniqKey="Seifert K">KA Seifert</name></author><author><name sortKey="Bruns, Td" uniqKey="Bruns T">TD Bruns</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ghannoum, Ma" uniqKey="Ghannoum M">MA Ghannoum</name></author><author><name sortKey="Jurevic, Rj" uniqKey="Jurevic R">RJ Jurevic</name></author><author><name sortKey="Mukherjee, Pk" uniqKey="Mukherjee P">PK Mukherjee</name></author><author><name sortKey="Cui, F" uniqKey="Cui F">F Cui</name></author><author><name sortKey="Sikaroodi, M" uniqKey="Sikaroodi M">M Sikaroodi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jumpponen, A" uniqKey="Jumpponen A">A Jumpponen</name></author><author><name sortKey="Jones, Kl" uniqKey="Jones K">KL Jones</name></author><author><name sortKey="Mattox, D" uniqKey="Mattox D">D Mattox</name></author><author><name sortKey="Yaege, C" uniqKey="Yaege C">C Yaege</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tedersoo, L" uniqKey="Tedersoo L">L Tedersoo</name></author><author><name sortKey="Nilsson, Rh" uniqKey="Nilsson R">RH Nilsson</name></author><author><name sortKey="Abarenkov, K" uniqKey="Abarenkov K">K Abarenkov</name></author><author><name sortKey="Jairus, T" uniqKey="Jairus T">T Jairus</name></author><author><name sortKey="Sadam, A" uniqKey="Sadam A">A Sadam</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hibbett, Ds" uniqKey="Hibbett D">DS Hibbett</name></author><author><name sortKey="Ohman, A" uniqKey="Ohman A">A Ohman</name></author><author><name sortKey="Glotzer, D" uniqKey="Glotzer D">D Glotzer</name></author><author><name sortKey="Nuhn, M" uniqKey="Nuhn M">M Nuhn</name></author><author><name sortKey="Kirk, P" uniqKey="Kirk P">P Kirk</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Vancov, T" uniqKey="Vancov T">T Vancov</name></author><author><name sortKey="Keen, B" uniqKey="Keen B">B Keen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Turenne, Cy" uniqKey="Turenne C">CY Turenne</name></author><author><name sortKey="Sanche, Se" uniqKey="Sanche S">SE Sanche</name></author><author><name sortKey="Hoban, Dj" uniqKey="Hoban D">DJ Hoban</name></author><author><name sortKey="Karlowsky, Ja" uniqKey="Karlowsky J">JA Karlowsky</name></author><author><name sortKey="Kabani, Am" uniqKey="Kabani A">AM Kabani</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Carlsen, T" uniqKey="Carlsen T">T Carlsen</name></author><author><name sortKey="Aas, Ab" uniqKey="Aas A">AB Aas</name></author><author><name sortKey="Lindner, D" uniqKey="Lindner D">D Lindner</name></author><author><name sortKey="Vr Lstad, T" uniqKey="Vr Lstad T">T Vrålstad</name></author><author><name sortKey="Schumacher, T" uniqKey="Schumacher T">T Schumacher</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schloss, Pd" uniqKey="Schloss P">PD Schloss</name></author><author><name sortKey="Westcott, Sl" uniqKey="Westcott S">SL Westcott</name></author><author><name sortKey="Ryabin, T" uniqKey="Ryabin T">T Ryabin</name></author><author><name sortKey="Hall, Jr" uniqKey="Hall J">JR Hall</name></author><author><name sortKey="Hartmann, M" uniqKey="Hartmann M">M Hartmann</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Blaalid, R" uniqKey="Blaalid R">R Blaalid</name></author><author><name sortKey="Kumar, S" uniqKey="Kumar S">S Kumar</name></author><author><name sortKey="Nilsson, Rh" uniqKey="Nilsson R">RH Nilsson</name></author><author><name sortKey="Abarenkov, K" uniqKey="Abarenkov K">K Abarenkov</name></author><author><name sortKey="Kirk, Pm" uniqKey="Kirk P">PM Kirk</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Abarenkov, K" uniqKey="Abarenkov K">K Abarenkov</name></author><author><name sortKey="Tedersoo, L" uniqKey="Tedersoo L">L Tedersoo</name></author><author><name sortKey="Nilsson, Rh" uniqKey="Nilsson R">RH Nilsson</name></author><author><name sortKey="Vellak, K" uniqKey="Vellak K">K Vellak</name></author><author><name sortKey="Saar, I" uniqKey="Saar I">I Saar</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="K Ljalg, U" uniqKey="K Ljalg U">U Kõljalg</name></author><author><name sortKey="Larsson, Kh" uniqKey="Larsson K">KH Larsson</name></author><author><name sortKey="Abarenkov, K" uniqKey="Abarenkov K">K Abarenkov</name></author><author><name sortKey="Nilsson, Rh" uniqKey="Nilsson R">RH Nilsson</name></author><author><name sortKey="Alexander, Ij" uniqKey="Alexander I">IJ Alexander</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nakamura, Y" uniqKey="Nakamura Y">Y Nakamura</name></author><author><name sortKey="Cochrane, G" uniqKey="Cochrane G">G Cochrane</name></author><author><name sortKey="Karsch Mizrachi, I" uniqKey="Karsch Mizrachi I">I Karsch-Mizrachi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Walters, Wa" uniqKey="Walters W">WA Walters</name></author><author><name sortKey="Caporaso, Jg" uniqKey="Caporaso J">JG Caporaso</name></author><author><name sortKey="Lauber, Cl" uniqKey="Lauber C">CL Lauber</name></author><author><name sortKey="Berg Lyons, D" uniqKey="Berg Lyons D">D Berg-Lyons</name></author><author><name sortKey="Fierer, N" uniqKey="Fierer N">N Fierer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ihrmark, K" uniqKey="Ihrmark K">K Ihrmark</name></author><author><name sortKey="Bodeker, Itm" uniqKey="Bodeker I">ITM Bödeker</name></author><author><name sortKey="Cruz Martinez, K" uniqKey="Cruz Martinez K">K Cruz-Martinez</name></author><author><name sortKey="Friberg, H" uniqKey="Friberg H">H Friberg</name></author><author><name sortKey="Kubartova, A" uniqKey="Kubartova A">A Kubartova</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Toju, H" uniqKey="Toju H">H Toju</name></author><author><name sortKey="Tanabe, As" uniqKey="Tanabe A">AS Tanabe</name></author><author><name sortKey="Yamamoto, S" uniqKey="Yamamoto S">S Yamamoto</name></author><author><name sortKey="Sato, H" uniqKey="Sato H">H Sato</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Martin, Kj" uniqKey="Martin K">KJ Martin</name></author><author><name sortKey="Rygiewicz, Pt" uniqKey="Rygiewicz P">PT Rygiewicz</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Egger, Kn" uniqKey="Egger K">KN Egger</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Bazzicalupo, Al" uniqKey="Bazzicalupo A">AL Bazzicalupo</name></author><author><name sortKey="Balint, M" uniqKey="Balint M">M Bálint</name></author><author><name sortKey="Schmitt, I" uniqKey="Schmitt I">I Schmitt</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lefever, S" uniqKey="Lefever S">S Lefever</name></author><author><name sortKey="Pattyn, F" uniqKey="Pattyn F">F Pattyn</name></author><author><name sortKey="Hellemans, J" uniqKey="Hellemans J">J Hellemans</name></author><author><name sortKey="Vandesompele, J" uniqKey="Vandesompele J">J Vandesompele</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bellemain, E" uniqKey="Bellemain E">E Bellemain</name></author><author><name sortKey="Carlsen, T" uniqKey="Carlsen T">T Carlsen</name></author><author><name sortKey="Brochmann, C" uniqKey="Brochmann C">C Brochmann</name></author><author><name sortKey="Coissac, E" uniqKey="Coissac E">E Coissac</name></author><author><name sortKey="Taberlet, P" uniqKey="Taberlet P">P Taberlet</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nilsson, Rh" uniqKey="Nilsson R">RH Nilsson</name></author><author><name sortKey="Ryberg, M" uniqKey="Ryberg M">M Ryberg</name></author><author><name sortKey="Abarenkov, K" uniqKey="Abarenkov K">K Abarenkov</name></author><author><name sortKey="Sjokvist, E" uniqKey="Sjokvist E">E Sjökvist</name></author><author><name sortKey="Kristiansson, E" uniqKey="Kristiansson E">E Kristiansson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Monard, C" uniqKey="Monard C">C Monard</name></author><author><name sortKey="Gantner, S" uniqKey="Gantner S">S Gantner</name></author><author><name sortKey="Stenlid, J" uniqKey="Stenlid J">J Stenlid</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Feinstein, Lm" uniqKey="Feinstein L">LM Feinstein</name></author><author><name sortKey="Sul, Wj" uniqKey="Sul W">WJ Sul</name></author><author><name sortKey="Blackwood, Cb" uniqKey="Blackwood C">CB Blackwood</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lindahl, Bd" uniqKey="Lindahl B">BD Lindahl</name></author><author><name sortKey="Nilsson, Rh" uniqKey="Nilsson R">RH Nilsson</name></author><author><name sortKey="Tedersoo, L" uniqKey="Tedersoo L">L Tedersoo</name></author><author><name sortKey="Abarenkov, K" uniqKey="Abarenkov K">K Abarenkov</name></author><author><name sortKey="Carlsen, T" uniqKey="Carlsen T">T Carlsen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Davey, Ml" uniqKey="Davey M">ML Davey</name></author><author><name sortKey="Currah, Rs" uniqKey="Currah R">RS Currah</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Munzenberger, B" uniqKey="Munzenberger B">B Münzenberger</name></author><author><name sortKey="Schneider, B" uniqKey="Schneider B">B Schneider</name></author><author><name sortKey="Nilsson, Rh" uniqKey="Nilsson R">RH Nilsson</name></author><author><name sortKey="Bubner, B" uniqKey="Bubner B">B Bubner</name></author><author><name sortKey="Larsson, Kh" uniqKey="Larsson K">KH Larsson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ashkannejhad, S" uniqKey="Ashkannejhad S">S Ashkannejhad</name></author><author><name sortKey="Horton, Tr" uniqKey="Horton T">TR Horton</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS One</journal-id><journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">plosone</journal-id><journal-title-group><journal-title>PLoS ONE</journal-title></journal-title-group><issn pub-type="epub">1932-6203</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">24933453</article-id><article-id pub-id-type="pmc">4059633</article-id><article-id pub-id-type="publisher-id">PONE-D-14-00713</article-id><article-id pub-id-type="doi">10.1371/journal.pone.0097629</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biology and Life Sciences</subject><subj-group><subject>Ecology</subject><subj-group><subject>Plant Ecology</subject><subj-group><subject>Plant-Environment Interactions</subject></subj-group></subj-group><subj-group><subject>Biodiversity</subject><subject>Community Ecology</subject><subject>Microbial Ecology</subject></subj-group></subj-group><subj-group><subject>Microbiology</subject><subj-group><subject>Plant Microbiology</subject></subj-group></subj-group><subj-group><subject>Mycology</subject></subj-group><subj-group><subject>Organisms</subject><subj-group><subject>Fungi</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Ecology and Environmental Sciences</subject><subj-group><subject>Soil Science</subject><subj-group><subject>Soil Ecology</subject></subj-group></subj-group></subj-group></article-categories><title-group><article-title>Comparison and Validation of Some ITS Primer Pairs Useful for Fungal Metabarcoding Studies</article-title><alt-title alt-title-type="running-head">Primers Used in Fungal Metabarcoding Studies</alt-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Op De Beeck</surname><given-names>Michiel</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Lievens</surname><given-names>Bart</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Busschaert</surname><given-names>Pieter</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author"><name><surname>Declerck</surname><given-names>Stéphan</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author"><name><surname>Vangronsveld</surname><given-names>Jaco</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author"><name><surname>Colpaert</surname><given-names>Jan V.</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Centre for Environmental Sciences, Hasselt University, Hasselt, Limburg, Belgium</addr-line></aff><aff id="aff2"><label>2</label><addr-line>Department of Microbial and Molecular Systems (M2S), Catholic University of Leuven, Antwerp, Belgium</addr-line></aff><aff id="aff3"><label>3</label><addr-line>Earth & Life Institute, Catholic University of Louvain, Louvain-la-Neuve, Walloon Brabant, Belgium</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>Neilan</surname><given-names>Brett</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1"><addr-line>University of New South Wales, Australia</addr-line></aff><author-notes><corresp id="cor1">* E-mail: <email>jan.colpaert@uhasselt.be</email></corresp><fn fn-type="conflict"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="con"><p>Conceived and designed the experiments: MODB JV JVC. Performed the experiments: MODB. Analyzed the data: MODB PB. Contributed reagents/materials/analysis tools: PB BL SD. Wrote the paper: MODB.</p></fn></author-notes><pub-date pub-type="collection"><year>2014</year></pub-date><pub-date pub-type="epub"><day>16</day><month>6</month><year>2014</year></pub-date><volume>9</volume><issue>6</issue><elocation-id>e97629</elocation-id><history><date date-type="received"><day>17</day><month>1</month><year>2014</year></date><date date-type="accepted"><day>21</day><month>4</month><year>2014</year></date></history><permissions><copyright-year>2014</copyright-year><copyright-holder>Op De Beeck et al</copyright-holder><license><license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p></license></permissions><abstract><p>Current metabarcoding studies aiming to characterize microbial communities generally rely on the amplification and sequencing of relatively short DNA regions. For fungi, the internal transcribed spacer (ITS) region in the ribosomal RNA (rRNA) operon has been accepted as the formal fungal barcode. Despite an increasing number of fungal metabarcoding studies, the amplification efficiency of primers is generally not tested prior to their application in metabarcoding studies. Some of the challenges that metabarcoding primers should overcome efficiently are the amplification of target DNA strands in samples rich in non-target DNA and environmental pollutants, such as humic acids, that may have been co-extracted with DNA. In the current study, three selected primer pairs were tested for their suitability as fungal metabarcoding primers. The selected primer pairs include two primer pairs that have been frequently used in fungal metabarcoding studies (ITS1F/ITS2 and ITS3/ITS4) and a primer pair (ITS86F/ITS4) that has been shown to efficiently amplify the ITS2 region of a broad range of fungal taxa in environmental soil samples. The selected primer pairs were evaluated in a 454 amplicon pyrosequencing experiment, real-time PCR (qPCR) experiments and <italic>in silico</italic> analyses. Results indicate that experimental evaluation of primers provides valuable information that could aid in the selection of suitable primers for fungal metabarcoding studies. Furthermore, we show that the ITS86F/ITS4 primer pair outperforms other primer pairs tested in terms of <italic>in silico</italic> primer efficiency, PCR efficiency, coverage, number of reads and number of species-level operational taxonomic units (OTUs) obtained. These traits push the ITS86F/ITS4 primer pair forward as highly suitable for studying fungal diversity and community structures using DNA metabarcoding.</p></abstract><funding-group><funding-statement>The study was made possible by funding of the Research Foundation Flanders (FWO). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement></funding-group><counts><page-count count="11"></page-count></counts></article-meta></front><body><sec id="s1"><title>Introduction</title><p>Until the late 1980s, microbial ecologists and taxonomists relied on culturing and morphological and physiological characteristics to describe microbial communities and members thereof. In the last two decades, DNA sequencing has revolutionized the way microbial communities are being characterized <xref rid="pone.0097629-Stahl1" ref-type="bibr">[1]</xref>, <xref rid="pone.0097629-Hugenholtz1" ref-type="bibr">[2]</xref>. In addition, since the introduction of pyrosequencing by Margulies et al. <xref rid="pone.0097629-Margulies1" ref-type="bibr">[3]</xref>, characterization of microbial communities has undergone a second revolution as this technology (used by e.g. Sogin et al. <xref rid="pone.0097629-Sogin1" ref-type="bibr">[4]</xref> and Buée et al. <xref rid="pone.0097629-Bue1" ref-type="bibr">[5]</xref>) enables detailed microbial community characterization at greater sequencing depth than was deemed possible via cloning and Sanger sequencing. A number of next-generation sequencing technologies now enable researchers to identify a large number of organisms from environmental samples using relatively short DNA sequences. This molecular identification method has been termed metabarcoding <xref rid="pone.0097629-Taberlet1" ref-type="bibr">[6]</xref>. Nevertheless, whatever sequencing technology is used, DNA metabarcoding generally depends on the amplification of barcode regions using taxon-specific primers <xref rid="pone.0097629-Hebert1" ref-type="bibr">[7]</xref>. Such primers need to be universal enough to cover a large group of taxa (e.g. the fungal kingdom), but at the same time have to result in amplicons that are variable enough to efficiently distinguish between closely related species or to identify operational taxonomic units (OTUs) <xref rid="pone.0097629-Hebert1" ref-type="bibr">[7]</xref>, <xref rid="pone.0097629-Just1" ref-type="bibr">[8]</xref>. For fungi and oomycetes, the internal transcribed spacer region (ITS; spanning the ITS1, 5.8S and ITS2 regions) in the ribosomal RNA (rRNA) operon has been recognized as the formal DNA barcoding region <xref rid="pone.0097629-Seifert1" ref-type="bibr">[9]</xref>–<xref rid="pone.0097629-Schoch1" ref-type="bibr">[11]</xref>.</p><p>The full ITS region in fungi has an average length of 500 and 600 base pairs (bp) for ascomycetes and basidiomycetes, respectively, and an average length of 600 bp across all fungal lineages <xref rid="pone.0097629-Porter1" ref-type="bibr">[12]</xref>. As current 454 amplicon pyrosequencing (using Roche's Genome Sequencer FLX (GS-FLX) instrument and Titanium chemistry) generates read lengths averaging 450 bp, it is impossible to span the entire ITS region in a single run. Even with recent advances in sequencing technologies that enable sequencing across the entire ITS region, it will probably remain desirable for fungal metabarcoding studies to exclude the 5.8S region of the rRNA operon. The inclusion of conserved regions in DNA sequences is known to increase the risk of chimera formation during PCR <xref rid="pone.0097629-Haas1" ref-type="bibr">[13]</xref>. Therefore, generally, either the ITS1 or the ITS2 region is used in ecological studies aiming at the characterization of fungal communities.</p><p>Primers that will be used in metabarcoding studies should be able to efficiently amplify their target DNA regions in the presence of high concentrations of non-target DNA and contaminants, such as humic acids, that may have been co-extracted with DNA <xref rid="pone.0097629-Kosch1" ref-type="bibr">[14]</xref>. Therefore, <italic>in silico</italic> testing of primers is expected to result in an incomplete picture of how primers will behave during amplification of DNA extracted from environmental samples. Comparing the amplification efficiency and robustness of primers used in metabarcoding studies is important because differences in primer efficiency may result in strong biases in favour of more easily amplifiable sequences during PCR reactions, potentially influencing our view on fungal communities <xref rid="pone.0097629-Polz1" ref-type="bibr">[15]</xref>–<xref rid="pone.0097629-Engelbrektson1" ref-type="bibr">[17]</xref>. Moreover, a primer set that covers a large proportion of the species that compose a community of interest and that produces a reliable outcome is desired as ecological metabarcoding studies typically rely on a single primer pair to map microbial diversity.</p><p>The most commonly used primers in fungal ecology for sequence-based fungal identification at the species level were published by White et al. <xref rid="pone.0097629-White1" ref-type="bibr">[18]</xref>: ITS1, ITS2, ITS3 and ITS4, and by Gardes and Bruns <xref rid="pone.0097629-Gardes1" ref-type="bibr">[19]</xref>: ITS1F and ITS4B. Whereas the primers developed by White et al. <xref rid="pone.0097629-White1" ref-type="bibr">[18]</xref> had a broad spectrum, ITS1F and ITS4B were developed to be specific for fungi and basidiomycetes respectively <xref rid="pone.0097629-Gardes1" ref-type="bibr">[19]</xref>. ITS1F is most frequently combined with ITS2 to amplify the ITS1 region of the fungal rRNA operon and ITS3 is usually combined with ITS4 to amplify the ITS2 region. These primer pairs have been used in many branches of mycological research in the past twenty years and are popular tools in recent fungal community studies as well <xref rid="pone.0097629-Bue1" ref-type="bibr">[5]</xref>, <xref rid="pone.0097629-Jumpponen2" ref-type="bibr">[20]</xref>–<xref rid="pone.0097629-Tedersoo1" ref-type="bibr">[24]</xref> (also reviewed in Hibbett et al. <xref rid="pone.0097629-Hibbett1" ref-type="bibr">[25]</xref>).</p><p>The aim of the current study was to evaluate the amplification efficiency of these established primer pairs and to compare them to a selected primer pair (ITS86F/ITS4) that has been shown to specifically and efficiently amplify ITS sequences from a broad range of fungal taxa in human blood samples as well as in environmental soil samples <xref rid="pone.0097629-Vancov1" ref-type="bibr">[26]</xref>.</p></sec><sec sec-type="materials|methods" id="s2"><title>Materials and methods</title><sec id="s2a"><title>Study site and soil sampling</title><p>A pioneer pine forest on a stabilised sand dune in the northern part of Limburg, Belgium (Hechtel-Eksel: 51°7′33″N, 5°22′22″E) was selected to obtain samples for this study. The study site is not freely accessible. To gain access to this study site, please contact the responsible authorities (<xref ref-type="supplementary-material" rid="pone.0097629.s003">Table S1</xref>). The soil in this study site is a dry sandy soil without a litter layer, poor in organic matter and slightly acidic. The average organic carbon content for this site is 0.7% and the average pH is 4.7. The pioneer vegetation at the study site is dominated by young Scots Pine trees (<italic>Pinus sylvestris</italic> L.), mosses and lichens, with only few grasses and heather shrubs (<italic>Calluna vulgaris</italic> (L.) Hull). Tree ages at the time of sampling ranged from one to five years. The region has an average annual rainfall of 800 mm per square meter and the average annual temperature is 10°C (Royal Meteorological Institute, Ukkel, Belgium).</p><p>Soil samples for fungal community characterization were collected in November 2009. Samples were collected at a depth of 0 to 20 cm using a soil corer with a diameter of 1 cm. Four replicate soil samples were collected within a distance of ten centimetres from each other for seven sampling locations. Each sampling location was chosen close to a three to five year old pine tree randomly selected in the field. Selected pine trees were at least 20 m apart from each other. The 4 replicate soil samples were pooled for each sampling location, resulting in a total of seven pooled samples. Samples were sealed in plastic bags and tightly closed to prevent desiccation during transportation. Upon arrival in the lab, soil samples were sieved using a 2 mm sieve to homogenize the sample and remove roots, large pieces of organic matter and stones. Samples were subsequently stored at −80°C until DNA was extracted. No protected species were sampled during the study.</p></sec><sec id="s2b"><title>DNA extraction, PCR amplification and pyrosequencing</title><p>Approximately 250 mg of soil was used for each DNA extraction. DNA was extracted in quadruplicate from each pooled sample using the UltraClean Soil DNA Isolation Kit according to the manufacturer's protocol (MoBio, Carlsbad, CA, USA). This resulted in four replicates for each of seven pooled soil samples. Subsequently, amplicon libraries were created using barcode-tagged primers for the primer pairs ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4 (<xref ref-type="table" rid="pone-0097629-t001">Table 1</xref>). Both forward and reverse primers were synthesized with a tail containing the Roche 454 pyrosequencing adaptors and a sample-specific 10 bp barcode (multiplex identifiers: MIDs) <xref rid="pone.0097629-Carlsen1" ref-type="bibr">[28]</xref> enabling sorting out the obtained sequences after sequencing (Roche Applied Science, Mannheim, Germany). Fusion primers were designed according to the scheme provided in <xref ref-type="supplementary-material" rid="pone.0097629.s004">Table S2</xref>.</p><table-wrap id="pone-0097629-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0097629.t001</object-id><label>Table 1</label><caption><title>Primers used in the current study.</title></caption><alternatives><graphic id="pone-0097629-t001-1" xlink:href="pone.0097629.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Primer name</td><td align="left" rowspan="1" colspan="1">Primer sequence (5′-3′)</td><td align="left" rowspan="1" colspan="1">rRNA operon binding site</td><td align="left" rowspan="1" colspan="1">Reference</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">ITS1F (F)</td><td align="left" rowspan="1" colspan="1"><named-content content-type="gene">CTTGGTCATTTAGAGGAAGTAA</named-content></td><td align="left" rowspan="1" colspan="1">Small subunit</td><td align="left" rowspan="1" colspan="1"><xref rid="pone.0097629-Gardes1" ref-type="bibr">[19]</xref></td></tr><tr><td align="left" rowspan="1" colspan="1">ITS2 (R)</td><td align="left" rowspan="1" colspan="1"><named-content content-type="gene">GCTGCGTTCTTCATCGATGC</named-content></td><td align="left" rowspan="1" colspan="1">5.8S</td><td align="left" rowspan="1" colspan="1"><xref rid="pone.0097629-White1" ref-type="bibr">[18]</xref></td></tr><tr><td align="left" rowspan="1" colspan="1">ITS3 (F)</td><td align="left" rowspan="1" colspan="1"><named-content content-type="gene">GCATCGATGAAGAACGCAGC</named-content></td><td align="left" rowspan="1" colspan="1">5.8S</td><td align="left" rowspan="1" colspan="1"><xref rid="pone.0097629-White1" ref-type="bibr">[18]</xref></td></tr><tr><td align="left" rowspan="1" colspan="1">ITS4 (R)</td><td align="left" rowspan="1" colspan="1"><named-content content-type="gene">TCCTCCGCTTATTGATATGC</named-content></td><td align="left" rowspan="1" colspan="1">Large subunit</td><td align="left" rowspan="1" colspan="1"><xref rid="pone.0097629-White1" ref-type="bibr">[18]</xref></td></tr><tr><td align="left" rowspan="1" colspan="1">ITS86F (F)</td><td align="left" rowspan="1" colspan="1"><named-content content-type="gene">GTGAATCATCGAATCTTTGAA</named-content></td><td align="left" rowspan="1" colspan="1">5.8S</td><td align="left" rowspan="1" colspan="1"><xref rid="pone.0097629-Turenne1" ref-type="bibr">[27]</xref></td></tr><tr><td align="left" rowspan="1" colspan="1">ITS86R (R)</td><td align="left" rowspan="1" colspan="1"><named-content content-type="gene">TTCAAAGATTCGATGATTCA<underline>G</underline></named-content></td><td align="left" rowspan="1" colspan="1">5.8S</td><td align="left" rowspan="1" colspan="1"><xref rid="pone.0097629-Vancov1" ref-type="bibr">[26]</xref></td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt101"><label></label><p>Primers are indicated as forward (F) or reverse (R). ITS86R contains a wrong base at the 3′ end. The <underline>G</underline> should be replaced by a C (see <xref ref-type="sec" rid="s4">Discussion</xref>).</p></fn></table-wrap-foot></table-wrap><p>DNA samples were amplified using a Techne TC-5000 thermocycler (Bibby Scientific Limited, Staffordshire, UK) under the following conditions: initial denaturation at 95°C for 2 minutes, followed by 40 cycles of denaturation at 95°C for 30 seconds, annealing at 55°C for 30 seconds and extension at 72°C during 1 minute; a final extension phase was performed at 72°C during 10 minutes. Reactions were carried out in 25 µl reaction volumes using the FastStart High Fidelity PCR System (Roche Applied Science, Mannheim, Germany). Each reaction contained 2.75 µl FastStart 10× reaction buffer, 1.8 mM MgCl, 0.2 mM dNTP mix, 0.4 µM of each primer, 1.25 U FastStart HiFi polymerase and 5 ng template DNA (as measured by a Nanodrop spectrophotometer).</p><p>Amplified DNA was cleared from PCR primers and primer dimers using the Agencourt AMpure XP System according to the manufacturer's protocol (Beckman Coulter, Brea, CA, USA). Finally, purified dsDNA was quantified with the Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA) and a Fluostar Omega plate reader (BMG Labtech, Ortenberg, Germany) and subsequently pooled in equimolar concentrations. The resulting amplicon pool, containing all 84 samples, was sequenced on one fourth of a Pico Titer Plate on a Roche Genome Sequencer FLX System using Titanium chemistry (Roche Applied Science, Mannheim, Germany) according to the manufacturer's instructions.</p></sec><sec id="s2c"><title>Bioinformatics processing</title><p>The standard flowgram format (SFF) file that resulted from the interpreted flowgrams was deposited in the NCBI Sequence Read Archive under accession number SRP026207 (SRA, <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/Traces/sra">http://www.ncbi.nlm.nih.gov/Traces/sra</ext-link>). From the original SFF file, three separate quality and fasta files were created with a custom Python script according to the three primer pairs used (<xref ref-type="supplementary-material" rid="pone.0097629.s004">Table S2</xref>). Further analyses were carried out in Mothur 1.31.2 on the individual fastq and fasta files <xref rid="pone.0097629-Schloss1" ref-type="bibr">[29]</xref>. Quality trimming in Mothur removed reads shorter than 200 bases, reads longer than 600 bases, reads with homopolymers longer than 8 bases and reads containing ambiguous bases. Reads were trimmed when the average Phred quality score dropped below 35 over a window of 50 bases. Next, sequences were compared to each other and duplicate sequences were replaced by a single sequence, while archiving the abundance data of the unique sequences. Subsequently, unique reads were checked for chimeric sequences using the Uchime tool in Mothur followed by their removal from the datasets. Unique reads were aligned with the pairwise alignment tool in Mothur. Finally, species-level OTUs were defined based on a 97% sequence similarity level, which is within the range of intraspecific ITS sequence similarity <xref rid="pone.0097629-Blaalid1" ref-type="bibr">[30]</xref>. In order to further remove potential sequencing errors from the analysis, global singletons (i.e. OTUs represented by only a single sequence over an entire dataset) were removed <xref rid="pone.0097629-Tedersoo1" ref-type="bibr">[24]</xref>.</p><p>Because the primer pairs resulted in different amounts of reads per sample, the number of reads per sample were rarefied to 200 reads per sample. Samples for which less than 200 reads were obtained were removed from the dataset. For ITS1F/ITS2 14 of 28 samples were removed. For ITS3/ITS4 4 samples were removed and for ITS86F/ITS4 no samples were removed. Inter-sample rarefaction curves were constructed based on 10,000 iterations. Subsequently, intra-sample diversity, richness and Good's coverage estimates were calculated in Mothur 1.31.2 based on 10,000 iterations. BLAST searches for a representative sequence of each OTU (as determined by Mothur) were conducted using the PlutoF v2.0 massBLASTer online tool <xref rid="pone.0097629-Abarenkov1" ref-type="bibr">[31]</xref>. Reads were blasted against the UNITE <xref rid="pone.0097629-Kljalg1" ref-type="bibr">[32]</xref> and INSD <xref rid="pone.0097629-Nakamura1" ref-type="bibr">[33]</xref> databases. Resulting HTML files were combined with the abundance data obtained in Mothur using a custom Python script. This script also acquired the names of species or genera that resemble Latin binomials with the highest BLAST score, avoiding unidentified OTUs in the databases to be seen as best BLAST hits. Unidentified OTUs were indicated as “not applicable (NA)”.</p></sec><sec id="s2d"><title>Quantitative real-time PCR</title><p>To evaluate the performance of the primer pairs amplifying target DNA from a heterogeneous pool of DNA in environmental samples, all primer pairs were tested in a qPCR set-up. A 2-fold dilution series (1∶1 to 1∶64) was made from twelve DNA samples (ranging from 5 ng µl<sup>−1</sup> to 78 pg µl<sup>−1</sup>, including one no-template control (NTC) for each sample). Amplification was performed in optical 96-well plates using a 7500 Fast Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) and SYBR Green chemistry. PCR conditions were as follows: initial denaturation at 95°C for two minutes, followed by 40 cycles of 95°C (30 s), 55°C (30 s) and 72°C (60 s) and a final extension phase at 72°C for 10 minutes followed by the generation of a dissociation curve to verify amplification specificity. These qPCR conditions were chosen to mimic the PCR conditions used during the PCR step prior to emPCR and amplicon pyrosequencing. Reactions contained 2.5 µL template DNA, 5 µL 2× Fast SYBR Green Master Mix (Applied Biosystems, Foster City, CA, USA), 0.3 µl forward and reverse primers (3.3 µM each) and 1.9 µL nuclease-free H<sub>2</sub>O in a total volume of 10 µL. PCR efficiencies (E) were calculated as E = (10<sup>−1/slope</sup>−1)×100.</p><p>To assess a potential PCR-bias at the phylum level, DNA was extracted from 15 pure cultures provided by the Mycothèque de l'Université Catholique de Louvain (BCCM/MUCL) including 5 basidiomycetes (<italic>Lentinula edodes</italic> (MUCL 44827), <italic>Agrocybe praecox</italic> (MUCL 46727), <italic>Coniophora marmorata</italic> (MUCL 39471), <italic>Suillus luteus</italic> (UH-Slu-LM8-n1) and <italic>Antrodia vaillantii</italic> (MUCL 54533)), 5 ascomycetes (<italic>Cladosporium cladosporioides</italic> (MUCL 53652), <italic>Cryptosporiopsis radicicola</italic> (MUCL 53485), <italic>Monilinia laxa</italic> (MUCL 30841), <italic>Arthroderma otae</italic> (MUCL 39756) and <italic>Galactomyces geotrichum</italic> (MUCL 52377)), 2 glomeromycetes (<italic>Rhizophagus clareus</italic> (MUCL 46238) and <italic>Rhizophagus</italic> sp. (MUCL 41833)) and 3 zygomycetes (<italic>Mortierella verticillata</italic> (MUCL 9658), <italic>Absidia corymbifera</italic> (MUCL 38907) and <italic>Mucor hiemalis</italic> (MUCL 15439), also see <xref ref-type="supplementary-material" rid="pone.0097629.s006">Table S4</xref>). DNA was extracted from cultures using the DNeasy Plant Mini Kit according to the manufacturer's instructions (Qiagen, Venlo, Netherlands). DNA concentrations extracted from pure cultures used for qPCR ranged from 5 ng µl<sup>−1</sup> to 20 ng µl<sup>−1</sup>. PCR bias at the phylum level was tested according to the qPCR protocol described above.</p></sec><sec id="s2e"><title><italic>In silico</italic> evaluation of primer pairs</title><p>To evaluate the primer-to-target mismatches <italic>in silico</italic>, primers were tested with PrimerProspector 1.0.1 <xref rid="pone.0097629-Walters1" ref-type="bibr">[34]</xref> against sequences downloaded from NCBI. Three sets of sequences were downloaded from NCBI containing only full-length fungal 5.8S, 18S and 28S sequences. Duplicate sequences were removed using Mothur 1.31.2. ITS1F was tested against 3,748 18S rDNA sequences. ITS2, ITS3 and ITS86F were tested against 4,421 5.8S rDNA sequences. ITS4 was tested against 4,270 28S rDNA sequences. For comparison, also all primers investigated by Ihrmark et al. <xref rid="pone.0097629-Ihrmark1" ref-type="bibr">[35]</xref> and Toju et al. <xref rid="pone.0097629-Toju1" ref-type="bibr">[36]</xref> were tested <xref rid="pone.0097629-White1" ref-type="bibr">[18]</xref>–<xref rid="pone.0097629-Gardes1" ref-type="bibr">[19]</xref>, <xref rid="pone.0097629-Turenne1" ref-type="bibr">[27]</xref>, <xref rid="pone.0097629-Ihrmark1" ref-type="bibr">[35]</xref>–<xref rid="pone.0097629-Egger1" ref-type="bibr">[38]</xref>. All tests were performed as described by Walters et al. <xref rid="pone.0097629-Walters1" ref-type="bibr">[34]</xref> using standard settings. Primer scores were calculated based on the following formula: weighted score = non-3′ mismatches ×0.40+3′ mismatches ×1.00+non-3′ gaps ×1.00+3′ gaps ×3.00. An additional penalty score of 3.00 was assigned if the final 3′ base of a primer had a mismatch with its target sequence <xref rid="pone.0097629-Walters1" ref-type="bibr">[34]</xref>.</p></sec><sec id="s2f"><title>Statistical analysis</title><p>Statistical analyses were conducted in R 2.13.0 (The R Foundation for Statistical Computing, Vienna, Austria). Normal distributions of the residuals of models were checked with the Shapiro-Wilk test, while homoscedasticity of variances was analysed using either Bartlett's or the Fligner-Killeen test. Depending on the distribution of the estimated parameters, either ANOVA or the Kruskal-Wallis Rank Sum Test was used to check for significant differences in variances of parameters. Two-by-two comparisons were conducted using either Tukey's Honest Significant Differences tests or Pairwise Wilcoxon Rank Sum Tests. Poisson corrections were implemented for abundance data. Distributions of ratios were compared with Pearson's Chi-squared tests. Non-metric multi-dimensional scaling (NMDS) was performed using the Vegan 2.0 - 8 package in R.</p></sec></sec><sec id="s3"><title>Results</title><sec id="s3a"><title>Parametrical analysis of 454 amplicon pyrosequencing data</title><p>For the three tested primer pairs, GS-FLX sequencing of the amplicon libraries generated a total of 151,650 reads. For a read to be successfully assigned to a sample, we required that both the forward and the reverse MIDs and primers were identified in a read with no more than one erroneous base in the MIDs and no more than two erroneous bases in the primer sequences. Based on the primer and MID sequences, 65,133 reads were assigned to their respective sample and 86,517 reads remained unassigned. The average length of reads assigned to either ITS1F/ITS2, ITS3/ITS4 or ITS86F/ITS4 prior to quality checking and trimming was 314, 331 and 369 bp respectively (excluding primers). The average read length of the unassigned reads was 116 bp (including primers, data not shown).</p><p>Rarefaction curves were constructed showing the rarefied number of OTUs defined at a 97% sequence similarity threshold relative to the number of samples (<xref ref-type="fig" rid="pone-0097629-g001">Fig. 1</xref>). These results indicate that, on average, a higher OTU richness and a better coverage of the fungal community can be expected for the ITS86F/ITS4 and ITS3/ITS4 primer pairs. The lowest OTU richness and coverage was predicted for the ITS1F/ITS2 primer pair. As most rarefaction curves tended towards saturation, the sequencing depth was assumed to be sufficient to retrieve the most abundant fungal OTUs in analysed soil samples that are detectable by the respective primers and 454 amplicon pyrosequencing.</p><fig id="pone-0097629-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0097629.g001</object-id><label>Figure 1</label><caption><title>Rarefaction curves for each of the three primer pairs used in this study: ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4.</title><p>In these graphs, the number of samples is plotted against the rarefied number of operational taxonomic units (OTUs) that were created based on a 97% sequence similarity cut-off value.</p></caption><graphic xlink:href="pone.0097629.g001"></graphic></fig><p>To compare primer pair performance in the 454 amplicon pyrosequencing experiment, averages of the number of reads were calculated across replicates (four replicates per sample) and samples (seven samples) for each primer pair. The average number of reads per sample obtained by ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4 after quality trimming differed significantly (p<0.01) and primer pairs yielded on average (± standard error) 356 (±26), 523 (±43) and 797 (±34) high quality reads per sample, respectively (<xref ref-type="fig" rid="pone-0097629-g002">Fig. 2A</xref>). The average number of OTUs found for each primer pair at a 97% sequence similarity threshold (observed OTU richness) also differed significantly (p<0.01). The highest OTU richness was observed for ITS86F/ITS4 with an average of 62 OTUs per sample (min = 42; max = 106). ITS1F/ITS2 yielded on average 32 OTUs per sample (min = 15; max = 60), whereas ITS3/ITS4 resulted in an average of 50 OTUs per sample (min = 27; max = 76) (<xref ref-type="fig" rid="pone-0097629-g002">Fig. 2B</xref>). Diversity was estimated with the inverse Simpson index. The inverse Simpson index differed significantly between ITS86F/ITS4 and ITS1F/ITS2, whereas with ITS1F/ITS2 a lower diversity was found than with ITS86F/ITS4 (p = 0.04). However, no significant differences were found between ITS3/ITS4 and ITS1F/ITS2 or between ITS3/ITS4 and ITS86F/ITS4 (p = 0.31 and p = 0.53, respectively) (<xref ref-type="fig" rid="pone-0097629-g002">Fig. 2C</xref>). The average Good's coverage per sample obtained for ITS1F/ITS2 was 96.8% (min = 93.8%, max = 98.9%), whereas the average Good's coverage obtained for ITS3/ITS4 and ITS86F/ITS4 was 96.5% (min = 93.2%, max = 99.0%) and 97.5% (min = 95.3%, max = 99.6%) respectively (<xref ref-type="fig" rid="pone-0097629-g002">Fig. 2D</xref>). Significant differences in Good's coverage were found between ITS3/ITS4 and ITS86F/ITS4 (p<0.01). However, no significant differences were found between ITS1F/ITS2 and ITS3/ITS4 (p = 0.81) or between ITS1F/ITS2 and ITS86F/ITS4 (p = 0.31).</p><fig id="pone-0097629-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0097629.g002</object-id><label>Figure 2</label><caption><title>Parametrical comparison between the three primer pairs used in this study (ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4).</title><p>A. Average number of sequences obtained after quality trimming. B. Average number of operational taxonomic units (OTUs), based on a 97% sequence similarity cut-off value. C. Average inverse Simpson index. D. Average Good's coverage. Averages were calculated across replicates (four) and samples (seven) for each primer pair. Differences at the 95% significance level are indicated with an asterisk “*”.</p></caption><graphic xlink:href="pone.0097629.g002"></graphic></fig></sec><sec id="s3b"><title>Community similarity compared between primer pairs</title><p>To compare the fungal community characterized with ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4 at the species and phylum level, a representative sequence of each OTU (as selected by Mothur) was blasted against the UNITE and INSD databases using the massBLASTer tool in PlutoF v2.0 <xref rid="pone.0097629-Abarenkov1" ref-type="bibr">[31]</xref>. Relative frequency distributions of the obtained species-level OTUs and phyla were analysed with chi-squared tests for the different primer pairs, based on the average abundances across replicates (four) and samples (seven). Representative reads of OTUs that could not be coupled to an accession of either the UNITE or INSD databases were considered as unidentified OTUs (indicated as not applicable “NA” in <xref ref-type="supplementary-material" rid="pone.0097629.s005">Table S3</xref>). A total of 51 unidentified OTUs were found of which 50 were found with ITS86F/ITS4 and 1 with ITS3/ITS4. BLAST scores and corresponding E-values for all OTUs can be found in <xref ref-type="supplementary-material" rid="pone.0097629.s005">Table S3</xref>. At the species level, differences were observed between the fungal communities identified by the three primer pairs studied (p<0.01). To give an idea of the fungal communities identified by each primer pair, pie charts displaying the top ten most abundant OTUs were constructed covering 68%, 62% and 64% of all sequences obtained with ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4, respectively (<xref ref-type="fig" rid="pone-0097629-g003">Fig. 3</xref>). Using the ITS1F/ITS2 primer pair (targeting the ITS1 region) a total of 183 OTUs across all samples were observed, with the most abundant OTUs corresponding to <italic>Sistotrema</italic> sp. Fr. (27%), <italic>Rhizopogon luteolus</italic> Fr. (9%), <italic>Wilcoxina mikolae</italic> (Chin S. Yang & H.E. Wilcox) Chin S. Yang & Korf (8%), <italic>Cladophialophora minutissima</italic> M.L. Davey & Currah (7%), and <italic>Capronia</italic> sp. Sacc. (5%) (<xref ref-type="fig" rid="pone-0097629-g003">Fig. 3A</xref>). The primer pairs ITS3/ITS4 and ITS86F/ITS4 (targeting the ITS2 region) identified 333 and 346 OTUs across all samples, respectively. In line with ITS1F/ITS2, the fungal communities identified with ITS3/ITS4 and ITS86F/ITS4 were also dominated by <italic>Sistotrema</italic> sp. (21-19%), but the subdominant OTUs were not exactly the same (<xref ref-type="fig" rid="pone-0097629-g003">Fig. 3B,C</xref>). Interesting to note is that the soil samples were dominated by ectomycorrhizal and ericoid mycorrhizal fungi and mycobionts from lichens. Based on field observations, we assumed that the fungal community in the pioneer forest that was sampled in this study would be relatively species poor compared to old forest soils <xref rid="pone.0097629-Bue1" ref-type="bibr">[5]</xref> and that biotrophic fungi would dominate over saprotrophic ones. These assumptions were confirmed by all three primer pairs (<xref ref-type="fig" rid="pone-0097629-g003">Fig. 3</xref>). At the phylum level, differences in community composition were found between all primer pairs tested (p<0.01 for all comparisons) (<xref ref-type="fig" rid="pone-0097629-g004">Fig. 4</xref>). Nevertheless, the majority of OTUs identified by all tested primer pairs belonged to the phyla Ascomycota (56% to 71%), followed by Basidiomycota (14% to 17%). A minority of OTUs identified, belonged to the Zygomycota (3% to 4%), Chytridiomycota (3% to 4%) and Glomeromycota (0% to 3%) (<xref ref-type="fig" rid="pone-0097629-g004">Fig. 4</xref>).</p><fig id="pone-0097629-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0097629.g003</object-id><label>Figure 3</label><caption><title>Relative abundance for the top ten most abundant species-level operational taxonomic units (OTUs), based on a 97% sequence similarity cut-off value, obtained for each of the three primer pairs studied (ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4).</title><p>Reads that did not result in a BLAST hit against the UNITE or INSD databases were indicated as “not applicable (NA)”. Ecological functions of OTUs are indicated between brackets behind the OTU identities (ECM: ectomycorrhizal, ERM: ericoid mycorrhizal, SAP: saprotrophic, LICH: lichenized, END: endophytic). OTUs not belonging to the top ten most abundant OTUs were pooled in the category “Remaining taxa”. OTUs that appear exclusively in a single chart are indicated in grayscale. OTUs that can be found in multiple pie charts are indicated in colour. OTU abundance scores were averaged across replicates (four) and samples (seven). A. ITS1F/ITS2. B. ITS3/ITS4. C. ITS86F/ITS4.</p></caption><graphic xlink:href="pone.0097629.g003"></graphic></fig><fig id="pone-0097629-g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0097629.g004</object-id><label>Figure 4</label><caption><title>Relative number of OTUs belonging to different fungal phyla. OTUs that could not be assigned to a phylum were grouped together under “not applicable (NA)”.</title><p>Averages were calculated across replicates (four) and samples (seven). A. ITS1F/ITS2. B. ITS3/ITS4. C. ITS86F/ITS4.</p></caption><graphic xlink:href="pone.0097629.g004"></graphic></fig></sec><sec id="s3c"><title>Repeatability of metabarcoding results</title><p>The repeatability of the molecular identification of fungal OTUs from environmental samples was compared between the three tested primer pairs to assess their experimental robustness. Replicates of samples were compared for each primer pair using NMDS with Bray-Curtis dissimilarities. In this analysis, samples with a similar OTU-composition will have smaller Bray-Curtis distances than samples with more dissimilar OTU compositions. In general, for all three primer pairs, replicates from the same sample grouped closely together (especially for ITS3/ITS4) (<xref ref-type="supplementary-material" rid="pone.0097629.s001">Fig. S1</xref>). Hence, the results of molecular identification of fungal OTUs are fairly consistent between replicated samples using the current experimental set-up. In order to test the possibility that some OTUs are missed in metabarcoding analyses based on the amplification and sequencing of target DNA from a single DNA extraction, results from the four replicated DNA extractions of the same sample were compared (<xref ref-type="supplementary-material" rid="pone.0097629.s002">Fig. S2</xref>). This assessment was performed for the four most abundant OTUs, representing <italic>Sistotrema</italic> sp., <italic>Rhizopogon luteolus</italic>, <italic>Cladophialophora minutissima</italic> and <italic>Wilcoxina mikolae</italic>. From <xref ref-type="supplementary-material" rid="pone.0097629.s002">Fig. S2</xref>, it is clear that in some replicated extractions of the same sample abundant OTUs can be missed. These results indicate that PCR amplification and sequencing can best be performed on multiple DNA extractions from the same environmental sample that are pooled prior to PCR in order to obtain an accurate picture of a fungal community.</p></sec><sec id="s3d"><title>Efficiency of primer pairs studied</title><p>To test the amplification efficiency of the three primer pairs in a heterogeneous pool of DNA (environmental sample) a qPCR experiment was conducted. More specifically, a 2-fold dilution series, ranging from 1∶1 to 1∶64 dilutions of twelve randomly selected DNA samples were amplified with randomly selected ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4 primers with MIDs and 454 adaptors attached. For ITS1F/ITS2, exponential amplification was obtained between 24 and 32 PCR cycles for ten out of twelve samples (data not shown). For two samples no exponential amplification phase was obtained within 40 cycles with this primer pair. ITS3/ITS4 showed exponential amplification after 22 to 36 cycles for all twelve samples, whereas ITS86F/ITS4 already showed an exponential amplification phase after 20 to 31 cycles for all samples (data not shown). Average PCR efficiencies (± standard error) were calculated to be 76% (±4%) for ITS3/ITS4, 82% (±5%) for ITS1F/ITS2 and 97% (±6%) for ITS86F/ITS4 (<xref ref-type="table" rid="pone-0097629-t002">Table 2</xref>).</p><table-wrap id="pone-0097629-t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0097629.t002</object-id><label>Table 2</label><caption><title>Average PCR amplification efficiencies obtained for twelve environmental DNA samples using quantitative real-time PCR.</title></caption><alternatives><graphic id="pone-0097629-t002-2" xlink:href="pone.0097629.t002"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Primer pair</td><td align="left" rowspan="1" colspan="1">ITS1F_ITS2</td><td align="left" rowspan="1" colspan="1">ITS3_ITS4</td><td align="left" rowspan="1" colspan="1">ITS86F_ITS4</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Average (%)</td><td align="left" rowspan="1" colspan="1">82</td><td align="left" rowspan="1" colspan="1">76</td><td align="left" rowspan="1" colspan="1">97</td></tr><tr><td align="left" rowspan="1" colspan="1">Standard error (%)</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">6</td></tr><tr><td align="left" rowspan="1" colspan="1">Minimum (%)</td><td align="left" rowspan="1" colspan="1">64</td><td align="left" rowspan="1" colspan="1">67</td><td align="left" rowspan="1" colspan="1">78</td></tr><tr><td align="left" rowspan="1" colspan="1">Maximum (%)</td><td align="left" rowspan="1" colspan="1">97</td><td align="left" rowspan="1" colspan="1">103</td><td align="left" rowspan="1" colspan="1">120</td></tr></tbody></table></alternatives></table-wrap></sec><sec id="s3e"><title>Phylum-level PCR bias</title><p>qPCR amplification efficiency did not significantly differ between primer pairs tested (ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4), nor between phyla (Ascomycota, Basidiomycota, Glomeromycota and Zygomycota) (<xref ref-type="fig" rid="pone-0097629-g005">Fig. 5</xref>). Two-way ANOVA resulted in p = 0.14 for phylum and p = 0.59 for primer pair. Primer pair - rDNA target combinations with poor PrimerProspector scores tended to have slightly lower PCR efficiencies, but these differences were not significant. Strains of species used for this experiment and PCR efficiencies can be found in <xref ref-type="supplementary-material" rid="pone.0097629.s006">Table S4</xref> and <xref ref-type="supplementary-material" rid="pone.0097629.s007">Table S5</xref>, respectively.</p><fig id="pone-0097629-g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0097629.g005</object-id><label>Figure 5</label><caption><title>Phylum-level PCR bias assessed using qPCR.</title><p>Average PCR efficiencies were calculated for each phylum using 5 basidiomycetes, 5 ascomycetes, 2 glomeromycetes and 3 zygomycetes. Error bars represent standard errors. No significant differences between primer pairs and phyla were found at the 95% significance level.</p></caption><graphic xlink:href="pone.0097629.g005"></graphic></fig></sec><sec id="s3f"><title><italic>In silico</italic> evaluation of primers</title><p>In a final analysis, the primer-to-target mismatches of the three primer pairs used in this study were evaluated with PrimerProspector <xref rid="pone.0097629-Walters1" ref-type="bibr">[34]</xref>. PrimerProspector was used to calculate a score for each primer based on mismatches between primers and target DNA sequences. The closer the score of a primer is to 0, the fewer mismatches were detected between primers and target sequences. The average scores (± standard error) for primers used in our study were: ITS1F = 4.55 (±0.05), ITS2 = 0.70 (±0.03), ITS3 = 0.58 (±0.03), ITS4 = 3.96 (±0.04) and ITS86F = 0.52 (±0.02) (<xref ref-type="table" rid="pone-0097629-t003">Table 3</xref>). Moreover, it was found that 44% of the tested sequences had a mismatch with the last base at the 3′ end of primer ITS1F. This particular mismatch between the last base at the 3′ end of a primer sequence and a target sequence occurred with only 9%, 4%, 16% and 3% of the tested sequences for ITS2, ITS3, ITS4 and ITS86F respectively (<xref ref-type="table" rid="pone-0097629-t003">Table 3</xref>). For comparison, also the primers suggested by Ihrmark et al. <xref rid="pone.0097629-Ihrmark1" ref-type="bibr">[35]</xref> and Toju et al. <xref rid="pone.0097629-Toju1" ref-type="bibr">[36]</xref> were tested with PrimerProspector. Also in this analysis, ITS86F was found to have the best primer score of all tested primers (<xref ref-type="supplementary-material" rid="pone.0097629.s008">Table S6</xref>).</p><table-wrap id="pone-0097629-t003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0097629.t003</object-id><label>Table 3</label><caption><title><xref ref-type="sec" rid="s3">Results</xref> of <italic>in silico</italic> testing of primers using PrimerProspector 1.0.1 <xref rid="pone.0097629-Walters1" ref-type="bibr">[34]</xref>.</title></caption><alternatives><graphic id="pone-0097629-t003-3" xlink:href="pone.0097629.t003"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col><col align="center" span="1"></col></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Primer</td><td align="left" rowspan="1" colspan="1">Number of sequences tested</td><td align="left" rowspan="1" colspan="1">3′ end base mismatch (%)</td><td align="left" rowspan="1" colspan="1">Average score ± SE</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">ITS1F (F)</td><td align="left" rowspan="1" colspan="1">3748</td><td align="left" rowspan="1" colspan="1">44%</td><td align="left" rowspan="1" colspan="1">4.6±0.05</td></tr><tr><td align="left" rowspan="1" colspan="1">ITS2 (R)</td><td align="left" rowspan="1" colspan="1">4421</td><td align="left" rowspan="1" colspan="1">9%</td><td align="left" rowspan="1" colspan="1">0.7±0.03</td></tr><tr><td align="left" rowspan="1" colspan="1">ITS3 (F)</td><td align="left" rowspan="1" colspan="1">4421</td><td align="left" rowspan="1" colspan="1">4%</td><td align="left" rowspan="1" colspan="1">0.6±0.03</td></tr><tr><td align="left" rowspan="1" colspan="1">ITS4 (R)</td><td align="left" rowspan="1" colspan="1">4270</td><td align="left" rowspan="1" colspan="1">16%</td><td align="left" rowspan="1" colspan="1">4.0±0.04</td></tr><tr><td align="left" rowspan="1" colspan="1">ITS86F (F)</td><td align="left" rowspan="1" colspan="1">4421</td><td align="left" rowspan="1" colspan="1">3%</td><td align="left" rowspan="1" colspan="1">0.0±0.00</td></tr></tbody></table></alternatives><table-wrap-foot><fn id="nt102"><label></label><p>Primers are indicated as forward (F) or reverse (R). Average PrimerProspector scores are shown ± standard errors (SE).</p></fn></table-wrap-foot></table-wrap></sec></sec><sec id="s4"><title>Discussion</title><p>Amplification and sequencing of short, standard DNA regions (metabarcoding) is becoming an increasingly popular tool for the characterization of fungal communities. Nevertheless, in most fungal metabarcoding studies, primers are generally used without being tested for their efficiency to amplify heterogeneous DNA pools, which may affect our view on studied fungal communities. Whereas the most commonly used primers in fungal metabarcoding studies were designed in the 90 s for species identification of a limited number of focal species, environmental metabarcoding studies generally aim to characterize diverse communities in environmental samples. Hence, primers used for fungal metabarcoding should be able to amplify a broad range of target DNA sequences in a sample that is also rich in non-target DNA and that may contain environmental contaminants <xref rid="pone.0097629-Kennedy1" ref-type="bibr">[39]</xref>. Even though recent efforts have resulted in new primers that could amplify a large proportion of target fungal DNA sequences <xref rid="pone.0097629-Ihrmark1" ref-type="bibr">[35]</xref>, <xref rid="pone.0097629-Toju1" ref-type="bibr">[36]</xref>, an experimental evaluation of PCR efficiency and primer performance should be performed on real environmental samples.</p><p>Initially, also ITS1F/ITS86R was included in our study design, but this primer pair was discarded from the study as no amplification was obtained in exploratory PCR and gel-electrophoresis tests. A plausible explanation for this failure can be found in the fact that the reverse primer (ITS86R) used and reported by Turenne et al. <xref rid="pone.0097629-Turenne1" ref-type="bibr">[27]</xref> and Vancov and Keen <xref rid="pone.0097629-Vancov1" ref-type="bibr">[26]</xref> contains an incorrect base at the 3′ end of the primer sequence. In order to be the perfect reverse complement of ITS86F, the sequence of ITS86R should be <named-content content-type="gene">5′-TTCAAAGATTCGATGATTCA<underline>C</underline>-3′</named-content>, and not <named-content content-type="gene">5′-TTCAAAGATTCGATGATTCA<underline>G</underline>-3′</named-content> as reported. GS-FLX sequencing of the amplicon pool resulted in 151,650 raw reads prior to quality trimming. Of these reads, 65,133 were assigned to their respective sample and 86,517 reads remained unassigned. The unassigned reads were investigated manually revealing that the majority were primer sequences probably resulting from primer dimers in our sequenced amplicon pool. Most likely, these primer dimers were not sufficiently removed during post-PCR clean-up steps.</p><p>Rarefaction curves were constructed for each primer pair (<xref ref-type="fig" rid="pone-0097629-g001">Fig. 1</xref>). These rarefaction curves indicate that the highest rarefied OTU richness and best coverage of the fungal community can be expected for the ITS86F/ITS4 and ITS3/ITS4 primer pairs. The average observed number of reads and the average observed number of OTUs (derived from these reads at a 97% sequence identity cut-off) indeed were highest for the ITS86F/ITS4 primer pair (797 reads and 62 OTUs on average per sample) and the ITS3/ITS4 primer pair (523 reads, 50 OTUs) and were much lower for the ITS1F/ITS2 primer pair (356 reads and 32 OTUs) (<xref ref-type="fig" rid="pone-0097629-g002">Fig. 2</xref>). The average observed diversity per sample, estimated by the inverse Simpson index, did not differ between ITS3/ITS4 and ITS86F/ITS4, but was significantly lower for ITS1F/ITS2 (<xref ref-type="fig" rid="pone-0097629-g002">Fig. 2</xref>). Overall, the low number of OTUs per sample found in the current study, are in sharp contrast with the more than 1000 OTUs per gram of forest soil found by Buée et al. <xref rid="pone.0097629-Bue1" ref-type="bibr">[5]</xref> based on amplification with the ITS1F/ITS2 primer pair. This difference in richness may be explained by the fact that pioneer forests probably contain relatively fewer fungal species compared to old forest soils <xref rid="pone.0097629-Bue1" ref-type="bibr">[5]</xref>. Additionally, overestimation or underestimation of species richness can also originate from data handling and analysis <xref rid="pone.0097629-Bazzicalupo1" ref-type="bibr">[40]</xref>. Based on the <italic>in silico</italic> performance and high Good's coverage calculated for ITS86F, it can be expected that the 62 OTUs found on average per sample by the ITS86F/ITS4 primer pair is close to the real species richness for the pioneer ecosystem growing on stabilised sand dunes which was studied here. The 50 OTUs per sample found by ITS3/ITS4 and the 32 OTUs found by ITS1F/ITS2, are probably underestimations due to a more narrow primer spectrum and/or lower PCR efficiencies. The fact that a high Good's coverage was found for the ITS1F/ITS2 primer pair despite a low observed OTU richness indicates that this primer pair is unable to multiply the ITS1 region of a large number of fungi. This is also supported by the <italic>in silico</italic> analysis. In this analysis, ITS1F was shown to have the poorest primer score of 4.6 and its sequence was shown to have a mismatch at the final base at the 3′ end of the primer (having a detrimental effect on amplification efficiency <xref rid="pone.0097629-Lefever1" ref-type="bibr">[41]</xref>) with no less than 44% of the tested fungal sequences (<xref ref-type="table" rid="pone-0097629-t003">Table 3</xref>). The large number of mismatches between the ITS1F primer and its target sequences was previously also addressed by Bellemain et al. <xref rid="pone.0097629-Bellemain1" ref-type="bibr">[42]</xref> and Ihrmark et al. <xref rid="pone.0097629-Ihrmark1" ref-type="bibr">[35]</xref>. In comparison, the ITS4 primer was given a score of 4.0 and was found to have a primer-to-target mismatch at the 3′ end of the primer with only 16% of the tested sequences. For the ITS2, ITS3 and ITS86F primers a score of 0.7, 0.6 and 0.0 was obtained respectively (<xref ref-type="table" rid="pone-0097629-t003">Table 3</xref>). These primers were shown to have a mismatch at the 3′ end of the primer with only 9%, 4% and 4% of the tested sequences, respectively (<xref ref-type="table" rid="pone-0097629-t003">Table 3</xref>), illustrating their broad amplification potential. Furthermore, our <italic>in silico</italic> analyses indicated that the primers suggested by Ihrmark et al. <xref rid="pone.0097629-Ihrmark1" ref-type="bibr">[35]</xref> and Toju et al. <xref rid="pone.0097629-Toju1" ref-type="bibr">[36]</xref> had more mismatches to their respective target sequences than ITS86F.</p><p>To test how these parametrical differences would translate to amplification efficiency during PCR amplification preceding emulsion PCR (emPCR) and pyrosequencing, a first qPCR experiment was conducted. To this end, DNA was extracted from 12 soil samples and amplified with the same primer pairs used in the pyrosequencing experiment. The calculated PCR efficiencies were 82% for ITS1F/ITS2, 76% for ITS3/ITS4 and 97% for ITS86F/ITS4 (<xref ref-type="table" rid="pone-0097629-t002">Table 2</xref>). From these PCR efficiencies, it is clear that ITS86F/ITS4 amplified its target ITS regions with greater efficiency than the other two primer pairs. Contrary to our expectations from the <italic>in silico</italic> analysis, ITS3/ITS4 obtained a lower efficiency than the ITS1F/ITS2 primer pair. This could be explained by the fact that also other factors determine the amplification efficiency of PCR reactions beside binding and dissociation of primers to their target DNA sequences. Such factors include the temperature-dependent properties of target DNA sequences and primer sequences in the PCR mixture, the temperature-dependent behaviour of the used polymerase enzyme mixtures, the use of ROX as an endogenous reference dye, etc. <xref rid="pone.0097629-Kennedy1" ref-type="bibr">[39]</xref>. Alternatively, the range of target sequences that ITS1F and ITS2 may bind to during PCR amplification is smaller, but the sequences that do get bound by these primers are amplified efficiently.</p><p>To see whether differences in amplification efficiency between primer pairs would also be reflected in the identities of the OTUs identified in the 454 amplicon pyrosequencing experiment, a representative read for each OTU was blasted against the UNITE and INSD databases and the BLAST hits with the highest score and a species or genus name were used to reconstruct the fungal community for each primer pair (<xref ref-type="fig" rid="pone-0097629-g003">Fig. 3</xref>). According to all three primer pairs, the soil fungal community was dominated by an OTU corresponding to <italic>Sistotrema</italic> sp. Additionally, all primer sets produced a number of OTUs that were commonly identified by all primer pairs (<xref ref-type="fig" rid="pone-0097629-g003">Fig. 3</xref>). The community identified by the three tested primer pairs still differed significantly, however. These differences confirm the finding that targeting either the ITS1 or the ITS2 region may result in different pictures of the fungal communities at the OTU level, as was previously assessed by both <italic>in silico</italic><xref rid="pone.0097629-Nilsson1" ref-type="bibr">[43]</xref> and sequencing studies <xref rid="pone.0097629-Bazzicalupo1" ref-type="bibr">[40]</xref>, <xref rid="pone.0097629-Monard1" ref-type="bibr">[44]</xref>. In addition, it was found that primers targeting the same ITS region do not necessarily result in the same OTU composition (<xref ref-type="fig" rid="pone-0097629-g003">Fig. 3</xref>), highlighting the importance of primer choice in a given study. However, it needs to be noted that in comparative studies, it has been shown that an ecological signal can be much stronger than the differences in community composition originating from primer choice <xref rid="pone.0097629-Monard1" ref-type="bibr">[44]</xref>.</p><p>At the phylum level, significant differences between ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4 were found as well (<xref ref-type="fig" rid="pone-0097629-g004">Fig. 4</xref>). Although in varying proportions, all three primer pairs identified more OTUs belonging to ascomycetes (70%, 71% and 56% respectively) than basidiomycetes (17%, 14% and 15%), but also Chytridiomycota (3%, 4%, 4%), Glomeromycota (0%, 3% and 2%) and Zygomycota (3%, 3% and 4%) were detected (<xref ref-type="fig" rid="pone-0097629-g004">Fig. 4</xref>). This might suggest that more ascomycetes were present in the soil at the time of investigation. However, amplification of DNA from ascomycetes may be favoured relative to amplification of DNA from basidiomycetes as the ITS sequences for ascomycetes are generally shorter than basidiomycete ITS sequences (this is especially true for the ITS2 region <xref rid="pone.0097629-Porter1" ref-type="bibr">[12]</xref>) and amplification of shorter DNA fragments is favoured during PCR. Whereas in previous <italic>in silico</italic> analyses indeed a phylum-level bias was expected for some of the primers used <xref rid="pone.0097629-Bellemain1" ref-type="bibr">[42]</xref>, no such bias was found in the current study based on experimental data derived from qPCR of DNA extracted from 15 fungal species belonging to the major fungal phyla (<xref ref-type="fig" rid="pone-0097629-g005">Fig. 5</xref>).</p><p>Whatever the aim of a metabarcoding study, results obtained from metabarcoding need to be reliable. To assess the repeatability of the fungal metabarcoding experiment, we analysed four replicate DNA extractions of seven soil samples separately. The analysis of all replicates of samples revealed that replicated analysis of the same sample with a specific primer pair generally results in similar fungal community compositions (<xref ref-type="supplementary-material" rid="pone.0097629.s001">Fig. S1</xref>). This is especially true for the ITS3/ITS4 and ITS86F/ITS4 primer pairs as their replicated samples clustered nicely together. However, this is less true for the ITS1F/ITS2 primer pair, where replicates of samples tend to have greater projected Bray-Curtis distances (<xref ref-type="supplementary-material" rid="pone.0097629.s001">Fig. S1</xref>). Moreover, we have shown that it is possible to miss certain OTUs, even abundant ones, when one sequences amplicon pools that are constructed from a single DNA extraction (<xref ref-type="supplementary-material" rid="pone.0097629.s002">Fig. S2</xref>). It is therefore advisable to extract DNA from environmental samples in multiple replicates, pool the eluates and perform PCR and sequencing on the DNA from the mixed eluate. This observation is in line with other studies performed previously, demonstrating that at least three replicated extractions are required to obtain a DNA pool that is representative for the microbial community present in a given soil sample <xref rid="pone.0097629-Feinstein1" ref-type="bibr">[45]</xref>, <xref rid="pone.0097629-Lindahl1" ref-type="bibr">[46]</xref>.</p><p>Apart from the technical issues that were addressed in this study, our data also provided a glimpse at the fungal community present in the studied site. Based on field observations of above-ground basidiocarps, we assumed that pioneer pine forests in the Campine region in Belgium are dominated by biotrophic species (mostly lichens, ectomycorrhizal and ericoid mycorrhizal fungi) over saprotrophic species. All three primer pairs confirmed this assumption, but they found different fungal OTUs to be dominant. According to the results obtained with ITS1F/ITS2, the fungal community in the studied site was dominated by OTUs corresponding to <italic>Sistotrema</italic> sp. (27%), followed by <italic>Rhizopogon luteolus</italic> (9%), <italic>Wilcoxina mikolae</italic> (8%) and <italic>Cladophialophora minutissima</italic> (7%) (<xref ref-type="fig" rid="pone-0097629-g003">Fig. 3</xref>) <xref rid="pone.0097629-Davey1" ref-type="bibr">[47]</xref>. These OTUs were also found to be very important members of the studied community according to ITS3/ITS4 and ITS86F/ITS4 as they appeared in the top ten of the most abundant OTUs found by both primer pairs, although in varying proportions (<xref ref-type="fig" rid="pone-0097629-g003">Fig. 3</xref>). <italic>Sistotrema</italic> sp., likely an important member of our studied ecosystem, was recently shown to be polyphyletic, containing both ectomycorrhizal and saprotrophic taxa <xref rid="pone.0097629-Mnzenberger1" ref-type="bibr">[48]</xref>. The reads that were found in the current study correspond to <italic>Sistotrema</italic> strains that were sampled from ectomycorrhizal root tips of <italic>Pinus contorta</italic> Dougl. growing on coastal sand dunes <xref rid="pone.0097629-Ashkannejhad1" ref-type="bibr">[49]</xref>. This genus provides a fine example of the power of molecular tools, such as DNA metabarcoding, to draw attention to ecologically important, cryptic fungal species. Based on field observations alone (basidiocarps observations and root tip morphotying), we never expected this genus to be so abundant in this pioneer ecosystem.</p><sec id="s4a"><title>Concluding remarks</title><p>In many fungal metabarcoding studies universal primers from previous phylogenetic or ecological studies are used without first performing an evaluation of their spectrum and performance for high-throughput sequencing, potentially resulting in a biased description of fungal communities. Whereas <italic>in silico</italic> PCR analyses on sequences retrieved from sequence databases may suggest promising primers <xref rid="pone.0097629-Ihrmark1" ref-type="bibr">[35]</xref>, <xref rid="pone.0097629-Toju1" ref-type="bibr">[36]</xref>, we showed that an experimental set-up to evaluate their usefulness in practice provides complementary information on the actual performance of the primers for high-throughput sequencing of environmental samples. Indeed, here we demonstrated that the choice of primers has a significant impact on how fungal communities are translated into OTU communities and subsequent data analysis. As such, before setting up large-scale sequencing experiments, we recommend to first test a number of promising primer pairs, e.g. selected with <italic>in silico</italic> analyses, under real PCR conditions for a subset of the samples under investigation. In case an in-depth characterization of a fungal community is desired, the use of more than one primer pair is advisable. We also showed that quantitative real-time PCR, evaluating the efficiency of selected primer pairs, may help in selecting the most efficient primer pairs. After all, using primer pairs that are not very efficient in amplifying DNA from an environmental sample will undoubtedly result in a low number of reads, and consequently in biased community descriptions.</p><p>In this study, the primer pair ITS86F/ITS4, which amplifies the ITS2 region of the fungal rRNA operon, was shown to be the most suitable primer pair for the characterization of fungal communities with metabarcoding. This primer pair not only resulted in superior amplification efficiency leading to a significantly higher number of reads, but also yielded a high number of OTUs belonging to different phyla. In addition, this primer pair resulted in a robust amplification reaction for the broadest range of samples and across replicated extractions.</p></sec></sec><sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title><supplementary-material content-type="local-data" id="pone.0097629.s001"><label>Figure S1</label><caption><p><bold>NMDS comparing community dissimilarities (based on Bray-Curtis distances) between replicates of samples.</bold></p><p>(PDF)</p></caption><media xlink:href="pone.0097629.s001.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0097629.s002"><label>Figure S2</label><caption><p><bold>Bar charts displaying the relative number of reads identified by the three primer pairs studied (ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4) assigned to one of the four most abundant OTUs.</bold></p><p>(PDF)</p></caption><media xlink:href="pone.0097629.s002.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0097629.s003"><label>Table S1</label><caption><p><bold>Contact addresses for access to the study site.</bold></p><p>(PDF)</p></caption><media xlink:href="pone.0097629.s003.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0097629.s004"><label>Table S2</label><caption><p><bold>Primer design used in the current 454 amplicon pyrosequencing experiment.</bold></p><p>(PDF)</p></caption><media xlink:href="pone.0097629.s004.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0097629.s005"><label>Table S3</label><caption><p><bold>Complete list of OTUs identified in the 454 amplicon pyrosequencing experiment and corresponding BLAST scores and E-values.</bold></p><p>(PDF)</p></caption><media xlink:href="pone.0097629.s005.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0097629.s006"><label>Table S4</label><caption><p><bold>Species used to assess PCR bias of the studied primer pairs (ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4).</bold></p><p>(PDF)</p></caption><media xlink:href="pone.0097629.s006.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0097629.s007"><label>Table S5</label><caption><p><bold>PCR efficiencies for the amplification of the ITS region of the fungal rRNA operon with the studied primer pairs (ITS1F/ITS2, ITS3/ITS4 and ITS86F/ITS4).</bold></p><p>(PDF)</p></caption><media xlink:href="pone.0097629.s007.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pone.0097629.s008"><label>Table S6</label><caption><p><bold>PrimerProspector scores of primers used in the current study and investigated by Ihrmark et al. </bold><xref rid="pone.0097629-Ihrmark1" ref-type="bibr">[35]</xref><bold> and Toju et al. </bold><xref rid="pone.0097629-Toju1" ref-type="bibr">[36]</xref><bold>.</bold></p><p>(PDF)</p></caption><media xlink:href="pone.0097629.s008.pdf"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><ack><p>The authors would like to thank Remans T. for assisting in the qPCR experiments and Waud M. for performing the pyrosequencing run.</p></ack><ref-list><title>References</title><ref id="pone.0097629-Stahl1"><label>1</label><mixed-citation publication-type="journal"><name><surname>Stahl</surname><given-names>DA</given-names></name>, <name><surname>Lane</surname><given-names>DJ</given-names></name>, <name><surname>Olsen</surname><given-names>GJ</given-names></name>, <name><surname>Pace</surname><given-names>NR</given-names></name> (<year>1984</year>) <article-title>Analysis of hydrothermal vent-associated symbionts by ribosomal-RNA sequences</article-title>. <source>Science</source><volume>224</volume>: <fpage>409</fpage>–<lpage>411</lpage><pub-id pub-id-type="pmid">17741220</pub-id></mixed-citation></ref><ref id="pone.0097629-Hugenholtz1"><label>2</label><mixed-citation publication-type="journal"><name><surname>Hugenholtz</surname><given-names>P</given-names></name>, <name><surname>Pace</surname><given-names>NR</given-names></name> (<year>1996</year>) <article-title>Identifying microbial diversity in the natural environment: a molecular phylogenetic approach</article-title>. <source>Trends Biotechnol</source><volume>14</volume>: <fpage>190</fpage>–<lpage>197</lpage><pub-id pub-id-type="pmid">8663938</pub-id></mixed-citation></ref><ref id="pone.0097629-Margulies1"><label>3</label><mixed-citation publication-type="journal"><name><surname>Margulies</surname><given-names>M</given-names></name>, <name><surname>Egholm</surname><given-names>M</given-names></name>, <name><surname>Altman</surname><given-names>WE</given-names></name>, <name><surname>Attyia</surname><given-names>S</given-names></name>, <name><surname>Bader</surname><given-names>JS</given-names></name>, <etal>et al</etal> (<year>2005</year>) <article-title>Genome sequencing in microfabricated high-density picolitre reactors</article-title>. <source>Nature</source><volume>437</volume>: <fpage>376</fpage>–<lpage>380</lpage><pub-id pub-id-type="pmid">16056220</pub-id></mixed-citation></ref><ref id="pone.0097629-Sogin1"><label>4</label><mixed-citation publication-type="journal"><name><surname>Sogin</surname><given-names>ML</given-names></name>, <name><surname>Morrison</surname><given-names>HG</given-names></name>, <name><surname>Huber</surname><given-names>JA</given-names></name>, <name><surname>Welch</surname><given-names>DW</given-names></name>, <name><surname>Huse</surname><given-names>SM</given-names></name>, <etal>et al</etal> (<year>2006</year>) <article-title>Microbial diversity in the deep sea and the underexplored “rare biosphere”</article-title>. <source>Proc Natl Acad Sci USA</source><volume>103</volume>: <fpage>12115</fpage>–<lpage>12120</lpage><pub-id pub-id-type="pmid">16880384</pub-id></mixed-citation></ref><ref id="pone.0097629-Bue1"><label>5</label><mixed-citation publication-type="journal"><name><surname>Buée</surname><given-names>M</given-names></name>, <name><surname>Reich</surname><given-names>M</given-names></name>, <name><surname>Murat</surname><given-names>C</given-names></name>, <name><surname>Morin</surname><given-names>E</given-names></name>, <name><surname>Nilsson</surname><given-names>RH</given-names></name>, <etal>et al</etal> (<year>2009</year>) <article-title>454 Pyrosequencing analyses of forest soils reveal an unexpectedly high fungal diversity</article-title>. <source>New Phytol</source><volume>184</volume>: <fpage>449</fpage>–<lpage>456</lpage><pub-id pub-id-type="pmid">19703112</pub-id></mixed-citation></ref><ref id="pone.0097629-Taberlet1"><label>6</label><mixed-citation publication-type="journal"><name><surname>Taberlet</surname><given-names>P</given-names></name>, <name><surname>Coissac</surname><given-names>E</given-names></name>, <name><surname>Pompanon</surname><given-names>F</given-names></name>, <name><surname>Brochmann</surname><given-names>C</given-names></name>, <name><surname>Willerslev</surname><given-names>E</given-names></name> (<year>2012</year>) <article-title>Towards next-generation biodiversity assessment using DNA metabarcoding</article-title>. <source>Mol Ecol</source><volume>8</volume>: <fpage>2045</fpage>–<lpage>2050</lpage><pub-id pub-id-type="pmid">22486824</pub-id></mixed-citation></ref><ref id="pone.0097629-Hebert1"><label>7</label><mixed-citation publication-type="journal"><name><surname>Hebert</surname><given-names>PDN</given-names></name>, <name><surname>Cywinska</surname><given-names>A</given-names></name>, <name><surname>Ball</surname><given-names>SL</given-names></name>, <name><surname>deWaard</surname><given-names>JR</given-names></name> (<year>2003</year>) <article-title>Biological identifications through DNA barcodes</article-title>. <source>Proc R Soc B - Biol Sci</source><volume>270</volume>: <fpage>313</fpage>–<lpage>321</lpage></mixed-citation></ref><ref id="pone.0097629-Just1"><label>8</label><mixed-citation publication-type="journal"><name><surname>Justé</surname><given-names>A</given-names></name>, <name><surname>Thomma</surname><given-names>BPHJ</given-names></name>, <name><surname>Lievens</surname><given-names>B</given-names></name> (<year>2008</year>) <article-title>Recent advances in molecular techniques to study microbial communities in food-associated matrices and processes</article-title>. <source>Food Microbiol</source><volume>25</volume>: <fpage>745</fpage>–<lpage>761</lpage><pub-id pub-id-type="pmid">18620966</pub-id></mixed-citation></ref><ref id="pone.0097629-Seifert1"><label>9</label><mixed-citation publication-type="journal"><name><surname>Seifert</surname><given-names>KA</given-names></name> (<year>2009</year>) <article-title>Progress towards DNA barcoding of fungi</article-title>. <source>Mol Ecol Resour</source><volume>9</volume>: <fpage>83</fpage>–<lpage>89</lpage><pub-id pub-id-type="pmid">21564968</pub-id></mixed-citation></ref><ref id="pone.0097629-Begerow1"><label>10</label><mixed-citation publication-type="journal"><name><surname>Begerow</surname><given-names>D</given-names></name>, <name><surname>Nilsson</surname><given-names>H</given-names></name>, <name><surname>Unterseher</surname><given-names>M</given-names></name>, <name><surname>Maier</surname><given-names>W</given-names></name> (<year>2010</year>) <article-title>Current state and perspectives of fungal DNA barcoding and rapid identification procedures</article-title>. <source>Appl Microbiol Biotechnol</source><volume>87</volume>: <fpage>99</fpage>–<lpage>108</lpage><pub-id pub-id-type="pmid">20405123</pub-id></mixed-citation></ref><ref id="pone.0097629-Schoch1"><label>11</label><mixed-citation publication-type="journal"><name><surname>Schoch</surname><given-names>CL</given-names></name>, <name><surname>Seifert</surname><given-names>KA</given-names></name>, <name><surname>Huhndorf</surname><given-names>S</given-names></name>, <name><surname>Robert</surname><given-names>V</given-names></name>, <name><surname>Spouge</surname><given-names>JL</given-names></name>, <etal>et al</etal> (<year>2012</year>) <article-title>Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi</article-title>. <source>Proc Natl Acad Sci USA</source><volume>109</volume>: <fpage>6241</fpage>–<lpage>6246</lpage><pub-id pub-id-type="pmid">22454494</pub-id></mixed-citation></ref><ref id="pone.0097629-Porter1"><label>12</label><mixed-citation publication-type="journal"><name><surname>Porter</surname><given-names>TM</given-names></name>, <name><surname>Golding</surname><given-names>GB</given-names></name> (<year>2011</year>) <article-title>Are similarity- or phylogeny-based methods more appropriate for classifying internal transcribed spacer (ITS) metagenomic amplicons?</article-title><source>New Phytol</source><volume>192</volume>: <fpage>775</fpage>–<lpage>782</lpage><pub-id pub-id-type="pmid">21806618</pub-id></mixed-citation></ref><ref id="pone.0097629-Haas1"><label>13</label><mixed-citation publication-type="journal"><name><surname>Haas</surname><given-names>BJ</given-names></name>, <name><surname>Gevers</surname><given-names>D</given-names></name>, <name><surname>Earl</surname><given-names>AM</given-names></name>, <name><surname>Feldgarden</surname><given-names>M</given-names></name>, <name><surname>Ward</surname><given-names>DV</given-names></name>, <etal>et al</etal> (<year>2011</year>) <article-title>Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons</article-title>. <source>Genome Res</source><volume>21</volume>: <fpage>494</fpage>–<lpage>504</lpage><pub-id pub-id-type="pmid">21212162</pub-id></mixed-citation></ref><ref id="pone.0097629-Kosch1"><label>14</label><mixed-citation publication-type="journal"><name><surname>Kosch</surname><given-names>TA</given-names></name>, <name><surname>Summers</surname><given-names>K</given-names></name> (<year>2013</year>) <article-title>Techniques for minimizing the effects of PCR inhibitors in the chytridiomycosis assay</article-title>. <source>Mol Ecol Resour</source><volume>13</volume>: <fpage>230</fpage>–<lpage>236</lpage><pub-id pub-id-type="pmid">23241137</pub-id></mixed-citation></ref><ref id="pone.0097629-Polz1"><label>15</label><mixed-citation publication-type="journal"><name><surname>Polz</surname><given-names>MF</given-names></name>, <name><surname>Cavanaugh</surname><given-names>CM</given-names></name> (<year>1998</year>) <article-title>Bias in template-to-product ratios in multitemplate PCR</article-title>. <source>Appl Environ Microbiol</source><volume>64</volume>: <fpage>3724</fpage>–<lpage>3730</lpage><pub-id pub-id-type="pmid">9758791</pub-id></mixed-citation></ref><ref id="pone.0097629-Jumpponen1"><label>16</label><mixed-citation publication-type="journal"><name><surname>Jumpponen</surname><given-names>A</given-names></name> (<year>2007</year>) <article-title>Soil fungal communities underneath willow canopies on a primary successional glacier forefront: rDNA sequence results can be affected by primer selection and chimeric data</article-title>. <source>Microb Ecol</source><volume>53</volume>: <fpage>233</fpage>–<lpage>246</lpage><pub-id pub-id-type="pmid">17106807</pub-id></mixed-citation></ref><ref id="pone.0097629-Engelbrektson1"><label>17</label><mixed-citation publication-type="journal"><name><surname>Engelbrektson</surname><given-names>A</given-names></name>, <name><surname>Kunin</surname><given-names>V</given-names></name>, <name><surname>Wrighton</surname><given-names>KC</given-names></name>, <name><surname>Zvenigorodsky</surname><given-names>N</given-names></name>, <name><surname>Chen</surname><given-names>F</given-names></name>, <etal>et al</etal> (<year>2010</year>) <article-title>Experimental factors affecting PCR-based estimates of microbial species richness and evenness</article-title>. <source>ISME J</source><volume>4</volume>: <fpage>642</fpage>–<lpage>647</lpage><pub-id pub-id-type="pmid">20090784</pub-id></mixed-citation></ref><ref id="pone.0097629-White1"><label>18</label><mixed-citation publication-type="book">White TJ, Bruns TD, Lee SB, Taylor JW (1990) Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ, editors. PCR protocols: a guide to methods and applications. United States: Academic Press. pp. 315–322.</mixed-citation></ref><ref id="pone.0097629-Gardes1"><label>19</label><mixed-citation publication-type="journal"><name><surname>Gardes</surname><given-names>M</given-names></name>, <name><surname>Bruns</surname><given-names>TD</given-names></name> (<year>1993</year>) <article-title>ITS primers with enhanced specificity for basidiomycetes – application to the identification of mycorrhizae and rusts</article-title>. <source>Mol Ecol</source><volume>2</volume>: <fpage>113</fpage>–<lpage>118</lpage><pub-id pub-id-type="pmid">8180733</pub-id></mixed-citation></ref><ref id="pone.0097629-Jumpponen2"><label>20</label><mixed-citation publication-type="journal"><name><surname>Jumpponen</surname><given-names>A</given-names></name>, <name><surname>Jones</surname><given-names>KL</given-names></name> (<year>2009</year>) <article-title>Massively parallel 454 sequencing indicates hyperdiverse fungal communities in temperate <italic>Quercus macrocarpa</italic> phyllosphere</article-title>. <source>New Phytol</source><volume>184</volume>: <fpage>438</fpage>–<lpage>448</lpage><pub-id pub-id-type="pmid">19674337</pub-id></mixed-citation></ref><ref id="pone.0097629-Amend1"><label>21</label><mixed-citation publication-type="journal"><name><surname>Amend</surname><given-names>AS</given-names></name>, <name><surname>Seifert</surname><given-names>KA</given-names></name>, <name><surname>Bruns</surname><given-names>TD</given-names></name> (<year>2010</year>) <article-title>Quantifying microbial communities with 454 pyrosequencing: does read abundance count?</article-title><source>Mol Ecol</source><volume>19</volume>: <fpage>5555</fpage>–<lpage>5565</lpage><pub-id pub-id-type="pmid">21050295</pub-id></mixed-citation></ref><ref id="pone.0097629-Ghannoum1"><label>22</label><mixed-citation publication-type="journal"><name><surname>Ghannoum</surname><given-names>MA</given-names></name>, <name><surname>Jurevic</surname><given-names>RJ</given-names></name>, <name><surname>Mukherjee</surname><given-names>PK</given-names></name>, <name><surname>Cui</surname><given-names>F</given-names></name>, <name><surname>Sikaroodi</surname><given-names>M</given-names></name>, <etal>et al</etal> (<year>2010</year>) <article-title>Characterization of the oral fungal microbiome (mycobiome) in healthy individuals</article-title>. <source>PLoS Pathog</source><volume>6</volume>: <fpage>e1000713</fpage><pub-id pub-id-type="pmid">20072605</pub-id></mixed-citation></ref><ref id="pone.0097629-Jumpponen3"><label>23</label><mixed-citation publication-type="journal"><name><surname>Jumpponen</surname><given-names>A</given-names></name>, <name><surname>Jones</surname><given-names>KL</given-names></name>, <name><surname>Mattox</surname><given-names>D</given-names></name>, <name><surname>Yaege</surname><given-names>C</given-names></name> (<year>2010</year>) <article-title>Massively parallel 454-sequencing of fungal communities in <italic>Quercus</italic> spp. ectomycorrhizas indicates seasonal dynamics in urban and rural sites</article-title>. <source>Mol Ecol</source><volume>19</volume>: <fpage>41</fpage>–<lpage>53</lpage><pub-id pub-id-type="pmid">20331769</pub-id></mixed-citation></ref><ref id="pone.0097629-Tedersoo1"><label>24</label><mixed-citation publication-type="journal"><name><surname>Tedersoo</surname><given-names>L</given-names></name>, <name><surname>Nilsson</surname><given-names>RH</given-names></name>, <name><surname>Abarenkov</surname><given-names>K</given-names></name>, <name><surname>Jairus</surname><given-names>T</given-names></name>, <name><surname>Sadam</surname><given-names>A</given-names></name>, <etal>et al</etal> (<year>2010</year>) <article-title>454 Pyrosequencing and Sanger sequencing of tropical mycorrhizal fungi provide similar results but reveal substantial methodological biases</article-title>. <source>New Phytol</source><volume>188</volume>: <fpage>291</fpage>–<lpage>301</lpage><pub-id pub-id-type="pmid">20636324</pub-id></mixed-citation></ref><ref id="pone.0097629-Hibbett1"><label>25</label><mixed-citation publication-type="journal"><name><surname>Hibbett</surname><given-names>DS</given-names></name>, <name><surname>Ohman</surname><given-names>A</given-names></name>, <name><surname>Glotzer</surname><given-names>D</given-names></name>, <name><surname>Nuhn</surname><given-names>M</given-names></name>, <name><surname>Kirk</surname><given-names>P</given-names></name>, <etal>et al</etal> (<year>2011</year>) <article-title>Progress in molecular and morphological taxon discovery in Fungi and options for formal classification of environmental sequences</article-title>. <source>Fungal Biol Rev</source><volume>25</volume>: <fpage>38</fpage>–<lpage>47</lpage></mixed-citation></ref><ref id="pone.0097629-Vancov1"><label>26</label><mixed-citation publication-type="journal"><name><surname>Vancov</surname><given-names>T</given-names></name>, <name><surname>Keen</surname><given-names>B</given-names></name> (<year>2009</year>) <article-title>Amplification of soil fungal community DNA using the ITS86F and ITS4 primers</article-title>. <source>FEMS Microbiol Lett</source><volume>296</volume>: <fpage>91</fpage>–<lpage>96</lpage><pub-id pub-id-type="pmid">19459948</pub-id></mixed-citation></ref><ref id="pone.0097629-Turenne1"><label>27</label><mixed-citation publication-type="journal"><name><surname>Turenne</surname><given-names>CY</given-names></name>, <name><surname>Sanche</surname><given-names>SE</given-names></name>, <name><surname>Hoban</surname><given-names>DJ</given-names></name>, <name><surname>Karlowsky</surname><given-names>JA</given-names></name>, <name><surname>Kabani</surname><given-names>AM</given-names></name> (<year>1999</year>) <article-title>Rapid identification of fungi by using the ITS2 genetic region and an automated fluorescent capillary electrophoresis system</article-title>. <source>J Clin Microbiol</source><volume>37</volume>: <fpage>1846</fpage>–<lpage>1851</lpage><pub-id pub-id-type="pmid">10325335</pub-id></mixed-citation></ref><ref id="pone.0097629-Carlsen1"><label>28</label><mixed-citation publication-type="journal"><name><surname>Carlsen</surname><given-names>T</given-names></name>, <name><surname>Aas</surname><given-names>AB</given-names></name>, <name><surname>Lindner</surname><given-names>D</given-names></name>, <name><surname>Vrålstad</surname><given-names>T</given-names></name>, <name><surname>Schumacher</surname><given-names>T</given-names></name>, <etal>et al</etal> (<year>2012</year>) <article-title>Don't make a mista(g)ke: is tag switching an overlooked source of error in amplicon pyrosequencing studies?</article-title><source>Fungal Ecol</source><volume>5</volume>: <fpage>747</fpage>–<lpage>749</lpage></mixed-citation></ref><ref id="pone.0097629-Schloss1"><label>29</label><mixed-citation publication-type="journal"><name><surname>Schloss</surname><given-names>PD</given-names></name>, <name><surname>Westcott</surname><given-names>SL</given-names></name>, <name><surname>Ryabin</surname><given-names>T</given-names></name>, <name><surname>Hall</surname><given-names>JR</given-names></name>, <name><surname>Hartmann</surname><given-names>M</given-names></name>, <etal>et al</etal> (<year>2009</year>) <article-title>Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities</article-title>. <source>Appl Environ Microbiol</source><volume>75</volume>: <fpage>7537</fpage>–<lpage>7541</lpage><pub-id pub-id-type="pmid">19801464</pub-id></mixed-citation></ref><ref id="pone.0097629-Blaalid1"><label>30</label><mixed-citation publication-type="journal"><name><surname>Blaalid</surname><given-names>R</given-names></name>, <name><surname>Kumar</surname><given-names>S</given-names></name>, <name><surname>Nilsson</surname><given-names>RH</given-names></name>, <name><surname>Abarenkov</surname><given-names>K</given-names></name>, <name><surname>Kirk</surname><given-names>PM</given-names></name>, <etal>et al</etal> (<year>2013</year>) <article-title>ITS1 versus ITS2 as DNA metabarcodes for fungi</article-title>. <source>Mol Ecol Resour</source><volume>13</volume>: <fpage>218</fpage>–<lpage>224</lpage><pub-id pub-id-type="pmid">23350562</pub-id></mixed-citation></ref><ref id="pone.0097629-Abarenkov1"><label>31</label><mixed-citation publication-type="journal"><name><surname>Abarenkov</surname><given-names>K</given-names></name>, <name><surname>Tedersoo</surname><given-names>L</given-names></name>, <name><surname>Nilsson</surname><given-names>RH</given-names></name>, <name><surname>Vellak</surname><given-names>K</given-names></name>, <name><surname>Saar</surname><given-names>I</given-names></name>, <etal>et al</etal> (<year>2010</year>) <article-title>PlutoF - a web based workbench for ecological and taxonomic research, with an online implementation for fungal ITS sequences</article-title>. <source>Evol Bioinform</source><volume>6</volume>: <fpage>189</fpage>–<lpage>196</lpage></mixed-citation></ref><ref id="pone.0097629-Kljalg1"><label>32</label><mixed-citation publication-type="journal"><name><surname>Kõljalg</surname><given-names>U</given-names></name>, <name><surname>Larsson</surname><given-names>KH</given-names></name>, <name><surname>Abarenkov</surname><given-names>K</given-names></name>, <name><surname>Nilsson</surname><given-names>RH</given-names></name>, <name><surname>Alexander</surname><given-names>IJ</given-names></name>, <etal>et al</etal> (<year>2005</year>) <article-title>UNITE: a database providing web-based methods for the molecular identification of ectomycorrhizal fungi</article-title>. <source>New Phytol</source><volume>166</volume>: <fpage>1063</fpage>–<lpage>1068</lpage><pub-id pub-id-type="pmid">15869663</pub-id></mixed-citation></ref><ref id="pone.0097629-Nakamura1"><label>33</label><mixed-citation publication-type="journal"><name><surname>Nakamura</surname><given-names>Y</given-names></name>, <name><surname>Cochrane</surname><given-names>G</given-names></name>, <name><surname>Karsch-Mizrachi</surname><given-names>I</given-names></name> (<year>2013</year>) <article-title>The international nucleotide sequence database collaboration</article-title>. <source>Nucl Acids Res</source><volume>41</volume>: <fpage>D21</fpage>–<lpage>D24</lpage><pub-id pub-id-type="pmid">23180798</pub-id></mixed-citation></ref><ref id="pone.0097629-Walters1"><label>34</label><mixed-citation publication-type="journal"><name><surname>Walters</surname><given-names>WA</given-names></name>, <name><surname>Caporaso</surname><given-names>JG</given-names></name>, <name><surname>Lauber</surname><given-names>CL</given-names></name>, <name><surname>Berg-Lyons</surname><given-names>D</given-names></name>, <name><surname>Fierer</surname><given-names>N</given-names></name>, <etal>et al</etal> (<year>2011</year>) <article-title>PrimerProspector: de novo design and taxonomic analysis of PCR primers</article-title>. <source>Bioinformatics</source><volume>27</volume>: <fpage>1159</fpage>–<lpage>1161</lpage><pub-id pub-id-type="pmid">21349862</pub-id></mixed-citation></ref><ref id="pone.0097629-Ihrmark1"><label>35</label><mixed-citation publication-type="journal"><name><surname>Ihrmark</surname><given-names>K</given-names></name>, <name><surname>Bödeker</surname><given-names>ITM</given-names></name>, <name><surname>Cruz-Martinez</surname><given-names>K</given-names></name>, <name><surname>Friberg</surname><given-names>H</given-names></name>, <name><surname>Kubartova</surname><given-names>A</given-names></name>, <etal>et al</etal> (<year>2012</year>) <article-title>New primers to amplify the fungal ITS2 region - evaluation by 454-sequencing of artificial and natural communities</article-title>. <source>FEMS Microbiol Ecol</source><volume>82</volume>: <fpage>666</fpage>–<lpage>677</lpage><pub-id pub-id-type="pmid">22738186</pub-id></mixed-citation></ref><ref id="pone.0097629-Toju1"><label>36</label><mixed-citation publication-type="journal"><name><surname>Toju</surname><given-names>H</given-names></name>, <name><surname>Tanabe</surname><given-names>AS</given-names></name>, <name><surname>Yamamoto</surname><given-names>S</given-names></name>, <name><surname>Sato</surname><given-names>H</given-names></name> (<year>2012</year>) <article-title>High-coverage ITS primers for the DNA-based identification of ascomycetes and basidiomycetes in environmental samples</article-title>. <source>PLoS One</source><volume>7</volume>: <fpage>e40863</fpage><pub-id pub-id-type="pmid">22808280</pub-id></mixed-citation></ref><ref id="pone.0097629-Martin1"><label>37</label><mixed-citation publication-type="journal"><name><surname>Martin</surname><given-names>KJ</given-names></name>, <name><surname>Rygiewicz</surname><given-names>PT</given-names></name> (<year>2005</year>) <article-title>Fungal-specific PCR primers developed for analysis of the ITS region of environmental DNA extracts</article-title>. <source>BMC Microbiol</source><volume>5</volume>: <fpage>28</fpage><pub-id pub-id-type="pmid">15904497</pub-id></mixed-citation></ref><ref id="pone.0097629-Egger1"><label>38</label><mixed-citation publication-type="journal"><name><surname>Egger</surname><given-names>KN</given-names></name> (<year>1995</year>) <article-title>Molecular analysis of ectomycorrhizal fungal communities</article-title>. <source>Canadian J Bot</source><volume>73</volume>: <fpage>1415</fpage>–<lpage>1415</lpage></mixed-citation></ref><ref id="pone.0097629-Kennedy1"><label>39</label><mixed-citation publication-type="book">Kennedy S, Oswald N (2011) PCR troubleshooting and optimization: the essential guide. United Kingdom: Caister Academic Press. 235 p.</mixed-citation></ref><ref id="pone.0097629-Bazzicalupo1"><label>40</label><mixed-citation publication-type="journal"><name><surname>Bazzicalupo</surname><given-names>AL</given-names></name>, <name><surname>Bálint</surname><given-names>M</given-names></name>, <name><surname>Schmitt</surname><given-names>I</given-names></name> (<year>2013</year>) <article-title>Comparison of ITS1 and ITS2 rDNA in 454 sequencing of hyperdiverse fungal communities</article-title>. <source>Fungal Ecol</source><volume>6</volume>: <fpage>102</fpage>–<lpage>109</lpage></mixed-citation></ref><ref id="pone.0097629-Lefever1"><label>41</label><mixed-citation publication-type="journal"><name><surname>Lefever</surname><given-names>S</given-names></name>, <name><surname>Pattyn</surname><given-names>F</given-names></name>, <name><surname>Hellemans</surname><given-names>J</given-names></name>, <name><surname>Vandesompele</surname><given-names>J</given-names></name> (<year>2013</year>) <article-title>Single-nucleotide polymorphisms and other mismatches reduce performance of quantitative PCR assays</article-title>. <source>Clin Chem</source><volume>59</volume>: <fpage>1470</fpage>–<lpage>1480</lpage><pub-id pub-id-type="pmid">24014836</pub-id></mixed-citation></ref><ref id="pone.0097629-Bellemain1"><label>42</label><mixed-citation publication-type="journal"><name><surname>Bellemain</surname><given-names>E</given-names></name>, <name><surname>Carlsen</surname><given-names>T</given-names></name>, <name><surname>Brochmann</surname><given-names>C</given-names></name>, <name><surname>Coissac</surname><given-names>E</given-names></name>, <name><surname>Taberlet</surname><given-names>P</given-names></name>, <etal>et al</etal> (<year>2010</year>) <article-title>ITS as an environmental DNA barcode for fungi: an <italic>in silico</italic> approach reveals potential PCR biases</article-title>. <source>BMC Microbiol</source><volume>10</volume>: <fpage>1</fpage>–<lpage>9</lpage><pub-id pub-id-type="pmid">20051107</pub-id></mixed-citation></ref><ref id="pone.0097629-Nilsson1"><label>43</label><mixed-citation publication-type="journal"><name><surname>Nilsson</surname><given-names>RH</given-names></name>, <name><surname>Ryberg</surname><given-names>M</given-names></name>, <name><surname>Abarenkov</surname><given-names>K</given-names></name>, <name><surname>Sjökvist</surname><given-names>E</given-names></name>, <name><surname>Kristiansson</surname><given-names>E</given-names></name> (<year>2009</year>) <article-title>The ITS region as a target for characterization of fungal communities using emerging sequencing technologies</article-title>. <source>FEMS Microbiol Lett</source><volume>296</volume>: <fpage>97</fpage>–<lpage>101</lpage><pub-id pub-id-type="pmid">19459974</pub-id></mixed-citation></ref><ref id="pone.0097629-Monard1"><label>44</label><mixed-citation publication-type="journal"><name><surname>Monard</surname><given-names>C</given-names></name>, <name><surname>Gantner</surname><given-names>S</given-names></name>, <name><surname>Stenlid</surname><given-names>J</given-names></name> (<year>2013</year>) <article-title>Utilizing ITS1 and ITS2 to study environmental fungal diversity using pyrosequencing</article-title>. <source>FEMS Microbiol Ecol</source><volume>84</volume>: <fpage>165</fpage>–<lpage>175</lpage><pub-id pub-id-type="pmid">23176677</pub-id></mixed-citation></ref><ref id="pone.0097629-Feinstein1"><label>45</label><mixed-citation publication-type="journal"><name><surname>Feinstein</surname><given-names>LM</given-names></name>, <name><surname>Sul</surname><given-names>WJ</given-names></name>, <name><surname>Blackwood</surname><given-names>CB</given-names></name> (<year>2009</year>) <article-title>Assessment of bias associated with incomplete extraction of microbial DNA from soil</article-title>. <source>Appl Environ Microbiol</source><volume>75</volume>: <fpage>5428</fpage>–<lpage>5433</lpage><pub-id pub-id-type="pmid">19561189</pub-id></mixed-citation></ref><ref id="pone.0097629-Lindahl1"><label>46</label><mixed-citation publication-type="journal"><name><surname>Lindahl</surname><given-names>BD</given-names></name>, <name><surname>Nilsson</surname><given-names>RH</given-names></name>, <name><surname>Tedersoo</surname><given-names>L</given-names></name>, <name><surname>Abarenkov</surname><given-names>K</given-names></name>, <name><surname>Carlsen</surname><given-names>T</given-names></name>, <etal>et al</etal> (<year>2013</year>) <article-title>Fungal community analysis by high-throughput sequencing of amplified markers - a user's guide</article-title>. <source>New Phytol</source><volume>199</volume>: <fpage>288</fpage>–<lpage>299</lpage><pub-id pub-id-type="pmid">23534863</pub-id></mixed-citation></ref><ref id="pone.0097629-Davey1"><label>47</label><mixed-citation publication-type="journal"><name><surname>Davey</surname><given-names>ML</given-names></name>, <name><surname>Currah</surname><given-names>RS</given-names></name> (<year>2007</year>) <article-title>A new species of <italic>Cladophialophora</italic> (hyphomycetes) from boreal and montane bryophytes</article-title>. <source>Mycol Res</source><volume>111</volume>: <fpage>106</fpage>–<lpage>116</lpage><pub-id pub-id-type="pmid">17169546</pub-id></mixed-citation></ref><ref id="pone.0097629-Mnzenberger1"><label>48</label><mixed-citation publication-type="journal"><name><surname>Münzenberger</surname><given-names>B</given-names></name>, <name><surname>Schneider</surname><given-names>B</given-names></name>, <name><surname>Nilsson</surname><given-names>RH</given-names></name>, <name><surname>Bubner</surname><given-names>B</given-names></name>, <name><surname>Larsson</surname><given-names>KH</given-names></name>, <etal>et al</etal> (<year>2012</year>) <article-title>Morphology, anatomy, and molecular studies of the ectomycorrhiza formed axenically by the fungus <italic>Sistotrema</italic> sp. (Basidiomycota)</article-title>. <source>Mycol Prog</source><volume>11</volume>: <fpage>817</fpage>–<lpage>826</lpage></mixed-citation></ref><ref id="pone.0097629-Ashkannejhad1"><label>49</label><mixed-citation publication-type="journal"><name><surname>Ashkannejhad</surname><given-names>S</given-names></name>, <name><surname>Horton</surname><given-names>TR</given-names></name> (<year>2006</year>) <article-title>Ectomycorrhizal ecology under primary succession on coastal sand dunes: interactions involving <italic>Pinus contorta</italic>, suilloid fungi and deer</article-title>. <source>New Phytol</source><volume>169</volume>: <fpage>345</fpage>–<lpage>354</lpage><pub-id pub-id-type="pmid">16411937</pub-id></mixed-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Bois/explor/CheneBelgiqueV2/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000228 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000228 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Bois
|area= CheneBelgiqueV2
|flux= Pmc
|étape= Corpus
|type= RBID
|clé= PMC:4059633
|texte= Comparison and Validation of Some ITS Primer Pairs Useful for Fungal Metabarcoding Studies
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:24933453" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a CheneBelgiqueV2
| This area was generated with Dilib version V0.6.27. |  |