- Research article
- Open Access
ITS as an environmental DNA barcode for fungi: an in silico approach reveals potential PCR biases
© Bellemain et al; licensee BioMed Central Ltd. 2010
Received: 9 April 2010
Accepted: 9 July 2010
Published: 9 July 2010
During the last 15 years the internal transcribed spacer (ITS) of nuclear DNA has been used as a target for analyzing fungal diversity in environmental samples, and has recently been selected as the standard marker for fungal DNA barcoding. In this study we explored the potential amplification biases that various commonly utilized ITS primers might introduce during amplification of different parts of the ITS region in samples containing mixed templates ('environmental barcoding'). We performed in silico PCR analyses with commonly used primer combinations using various ITS datasets obtained from public databases as templates.
Some of the ITS primers, such as ITS1-F, were hampered with a high proportion of mismatches relative to the target sequences, and most of them appeared to introduce taxonomic biases during PCR. Some primers, e.g. ITS1-F, ITS1 and ITS5, were biased towards amplification of basidiomycetes, whereas others, e.g. ITS2, ITS3 and ITS4, were biased towards ascomycetes. The assumed basidiomycete-specific primer ITS4-B only amplified a minor proportion of basidiomycete ITS sequences, even under relaxed PCR conditions. Due to systematic length differences in the ITS2 region as well as the entire ITS, we found that ascomycetes will more easily amplify than basidiomycetes using these regions as targets. This bias can be avoided by using primers amplifying ITS1 only, but this would imply preferential amplification of 'non-dikarya' fungi.
We conclude that ITS primers have to be selected carefully, especially when used for high-throughput sequencing of environmental samples. We suggest that different primer combinations or different parts of the ITS region should be analyzed in parallel, or that alternative ITS primers should be searched for.
Molecular identification through DNA barcoding of fungi has, during the last 15-20 years, become an integrated and essential part of fungal ecology research and has provided new insights into the diversity and ecology of many different groups of fungi (reviewed by [1–4]). Molecular identification has made it possible to study the ecology of fungi in their dominant but inconspicuous mycelial stage and not only by means of fruiting bodies. Interest in sequenced-based analysis of environmental samples ('environmental barcoding') has increased in the past decade as it allows to study abundance and species richness of fungi at a high rate and more reliably than conventional biotic surveys (e.g. [5–10]). The internal transcribed spacer (ITS) of nuclear DNA (nrDNA) is the preferred DNA barcoding marker both for the identification of single taxa and mixed environmental templates ('environmental DNA barcoding'). It has recently been proposed as the official primary barcoding marker for fungi (Deliberation of 37 mycologists from 12 countries at the Smithsonian's Conservation and Research Centre, Front Royal, Virginia, May 2007). More than 100 000 fungal ITS sequences generated by conventional Sanger sequencing are deposited in the International Nucleotide Sequence Databases and/or other databases , providing a large reference material for identification of fungal taxa. However, these data are to some extent hampered by misidentifications or technical errors such as mixing of DNA templates or sequencing errors . Furthermore, a large amount of partial ITS sequences generated by next-generation sequencing has recently been deposited in public sequence databases.
The aim of this study was to analyse the biases commonly used ITS primers might introduce during PCR amplification. First, we addressed to what degree the various primers mismatch with the target sequence and whether the mismatches are more widespread in some taxonomic groups. Second, we considered the length variation in the amplified products, in relation to taxonomic group, to assess amplification biases during real (in vitro) PCR amplification, as shorter DNA fragments are preferentially amplified from environmental samples containing DNA from a mixture of different species . Finally, we analyzed to what degree the various primers co-amplify plants, which often co-occur in environmental samples. For these purposes we performed in silico PCR using various primer combinations on target sequences retrieved from EMBL databases as well as subset databases using the bioinformatic tool EcoPCR . In order to better simulate real PCR conditions, we allowed a maximum of 0 to 3 mismatches except for the 2 last bases of each primer and we assessed the melting temperature (Tm) for each primer in relation to primer mismatches.
Compilation of datasets
The EcoPCR package contains a set of bioinformatics tools developed at the Laboratoire d'Ecologie Alpine, Grenoble, France (, freely available at http://www.grenoble.prabi.fr/trac/ecoPCR). The package is composed of four pieces of software, namely 'ecoPCRFormat', 'ecoFind', 'ecoPCR' and 'ecoGrep'. Briefly, EcoPCR is based on the pattern matching algorithm agrep  and selects sequences from a database that match (exhibit similarity to) two PCR primers. The user can specify (1) which database the given primers should be tested against, and (2) the primer sequences. Different options allow specification of the minimum and maximum amplification length, the maximum count of mismatched positions between each primer and the target sequence (excluding the two bases on the 3'end of each primer), and restriction of the search to given taxonomic groups. The ecoPCR output contains, for each target sequence, amplification length, melting temperature (Tm), taxonomic information as well as the number of mismatched positions for each strand.
First, we retrieved from EMBL sequences from fungi in the following categories: 'standard', 'Genome sequence scan', 'High Throughput Genome sequencing', 'Whole Genome Sequence' from ftp://ftp.ebi.ac.uk/pub/databases/embl/release/ (release embl_102, January 2010) to create our initial database. It corresponds to 1,212,954 sequences including approximately 79,500 ITS sequences (estimated from EMBL SRS website requesting for fungi sequences annotated with 'ITS' or 'Internal Transcribed Spacer'). These ITS entries refer to more than 10,800 taxa. This database hereafter referred to as the "fungi database" was compiled using EcoPCRFormat.
To assess the specificity of the primers to fungi, we used the plant database from EMBL (release embl_102, January 2010 from ftp://ftp.ebi.ac.uk/pub/databases/embl/release/) to run amplifications using the same primers as for fungi. This database, hereafter referred to as the "plant database", contained 1,253,565 sequences, including approximately 65,000 ITS sequences (estimated from EMBL SRS website requesting for viridiplantae sequences annotated with 'ITS' or 'Internal Transcribed Spacer'). These ITS entries refer to more than 6,100 taxa. This database was also compiled using EcoPCRFormat.
As there are relatively few sequences submitted to public databases covering the entire ITS region as well as the commonly used universal primer sites in the flanking SSU and LSU regions, we created three subset datasets covering either ITS1, ITS2 or the entire ITS region. From the initial fungi database, we compiled three subset databases (hereafter referred to as subset 1, 2, and 3) by in silico amplification (see below) of target sequences using the following primer pairs: NS7-ITS2 (dataset 1, focused on ITS1 region), ITS5-ITS4 (dataset 2, including both ITS1 and ITS2 regions) and ITS3-LR3 (dataset 3, focused on ITS2 region). To simulate relatively stringent PCR conditions, a single mismatch between each primer and the template was allowed except in the 2 bases of the 3' primer end. These three subsets were then compiled using EcoPCRFormat and included 1291, 5924 and 2459 partial nrDNA sequences, respectively.
In silico amplification and primer specificity to fungi
Using EcoPCR, we ran in silico amplifications from both the fungi and the plant databases using various commonly used primer combinations, to assess the number of amplifications and the specificity of the primers to fungi. For each amplification, we allowed from 0 to 3 mismatches between each primer and the template (excluding mismatches in the 2 bases of the 3' primer end) in order to simulate different stringency conditions of PCRs. Secondly, from the three subsets, we amplified sequences using different internal primer combinations in order to evaluate the various primers (Figure 1). From dataset 1 we used the primer combinations ITS1-F-ITS2, ITS5-ITS2 and ITS1-ITS2. From dataset 2 we used the combinations ITS1-ITS4 (amplifying both ITS1 and ITS2 introns), ITS3-ITS4 and ITS5-ITS2. From dataset 3 we used the combinations ITS3-ITS4 and ITS3-ITS4B. During these virtual PCRs we also allowed from 0 to 3 mismatches between each primer and the template, except in the 2 bases of the 3' primer end.
Assessing the degree of primer mismatches and Tm
For all in silico internal amplifications from each subset, we assessed the proportion of sequences retrieved when allowing for 0 to 3 mismatches between each primer and the template. For the amplifications from each subset, we used an external primer (one of the primers used to create the subset) and an internal primer. Therefore, for each analysis, we assessed the proportion of sequences including mismatches for the internal primer only. The primer pair ITS5-ITS2 was evaluated both for subset 1 and subset 2, with the focus on ITS5 for subset 1 and on ITS2 for subset 2 (as those primers correspond to internal primers within their respective subsets). Similarly, the primer pair ITS3-ITS4 was evaluated both for subsets 2 and 3, with the focus on ITS3 in subset 2 and ITS4 in subset 3. The primer ITS1 was evaluated both for subset 1 (with the combination ITS1-ITS2) and for subset 2 (with the combination ITS1-ITS4) as ITS2 and ITS4 were used as external primers in subsets 1 and 2, respectively.
To assess whether certain taxonomic groups were more prone to mismatches, we assessed the proportion of sequences including one mismatch for each of the three taxonomic groups 'ascomycetes', 'basidiomycetes' and 'non-dikarya' (the latter is a highly polyphyletic group including e.g. Blastocladiomycota, Chytridiomycota, Glomeromycota and Zygomycota ). We also assessed the Tm for each primer based on the analyses from internal amplifications, allowing a single mismatch. The Tm is defined as the temperature at which half of the DNA strands are in the double-helical state and half are in the "random-coil" states. The strength of hybridization between the primers and the template affects Tm. It is therefore informative to assess how Tm decreases as the number of mismatches increases, i.e. with less stringent PCR conditions. Tm was calculated in ecoPCR based on a thermodynamic nearest neighbor model . Exact computation was performed following .
Assessing bias in amplification length relative to taxonomic group
To further assess the taxonomic bias introduced by the use of the different primer pairs, we separated the amplified sequences from selected analyses into the groups 'ascomycetes', 'basidomycetes' and 'non-dikarya' based on their taxonomic identification number, using the ecoGrep tool. These selected analyses were (1) the three subsets, and (2) all internal amplifications within each subset with one mismatch allowed. The amplification length was reported for each analysis.
Relative amplification of different primer combinations from the fungi and plant databases
Number of plant and fungi ITS sequences amplified in silico from EMBL fungal and plant databases, using the various primer combinations and allowing none to three mismatches.
Fungal ITS sequences
Plant ITS sequences
Number of mismatches *
Primer mismatches in sequence subsets
Allowing one mismatch increased the proportion of amplified sequences from 36% to 91.6% for the commonly used primer ITS1-F, implying that more than half of the amplified sequences included one mismatch. ITS5 amplified the highest proportion of the sequences when allowing for a single mismatch (97.5%), and less than 10% of the sequences in each taxonomic group included one mismatch. The primer ITS1, on the other hand, only amplified 56.8% and 65.9% of the sequences from subsets one and two, respectively, when allowing no mismatches. Allowing three mismatches, ITS1 was still only able to amplify 92% of the sequences in subsets one and two. Allowing no mismatches, the complementary primers ITS2 and ITS3 amplified 79.4% and 77.3% of all sequences respectively, in subset 2. Allowing one mismatch, these numbers increased to 87.5 and 90%, respectively. Primer ITS4 amplified 74.9% of all sequences in subset 3 and this proportion only increased to 93.7% when allowing three mismatches. The assumed basidiomycete-specific primer ITS4-B amplified only 5.6% of the sequences in subset 3 under strict conditions (corresponding to 46% of the basidiomycetes sequences, see below) and up to 14.9% allowing 3 mismatches. However, about half of the sequences included a mismatch when a single mismatch was allowed.
Taxonomic bias for different primers
Percentage of sequences amplified in silico, allowing one mismatch, from ascomycetes, basidiomycetes and 'non-Dikarya' with different primer combinations and using the three sequence subsets 1-3 (see Material and Methods) as templates.
Melting temperature (Tm) of each primer according to the number of mismatches allowed between the primer and the target sequence.
Number of mismatches allowed
Taxonomic bias relative to length of the amplified region
Number (percentages) of sequences and amplified in each of the most common Basidiomycete groups, from the original subset3 and from the amplification of ITS3-ITS4-B from subset3, allowing no or 3 mismatches.
Other categories *
Our in silico analyses further indicate that most of the primers will introduce a taxonomic bias due to higher levels of mismatches in certain taxonomic groups. When allowing one mismatch (corresponding to rather stringent PCR conditions) we found that the primer pairs ITS1-F, ITS1 and ITS5 preferentially amplified basidiomycetes whereas the primer pairs ITS2, ITS3 and ITS4 preferentially amplified ascomycetes. This type of bias must also be considered before selecting primer pairs for a given study. Also in molecular surveys of protistan and prokaryotic diversity, it has been documented that different 16S primers target different parts of the diversity [32–34].
In addition, our results clearly demonstrate that basidiomycetes, on average, have significantly longer amplicon sequences than ascomycetes both for the whole ITS region, and the ITS2 region. This fact probably also introduces taxonomic bias during PCR amplification of environmental samples, since shorter fragments are more readily amplified compared to longer ones. In several studies, it has been demonstrated that a greater proportion of the diversity can be detected with short target sequences compared to longer ones [35, 36]. Hence, using the ITS2 region or the whole ITS region, a higher number of the ascomycetes will probably be targeted compared to basidiomycetes. This bias could be avoided by using primers amplifying ITS1 only, but this would imply a preferential amplification of the 'non-dikarya' fungi.
The in silico method used here allowed for the assessment of different parameters for commonly used ITS primers, including the length amplicons generated, taxonomic biases, and the consequences of primer mismatches. The results provide novel insights into the relative performance of commonly used ITS primer pairs. Our analyses suggest that studies using these ITS primers to retrieve the entire fungal diversity from environmental samples including mixed templates should use lower annealing temperatures than the recommended Tm to allow for primer mismatches. A high Tm has been used in most studies, which likely biases the inferred taxonomic composition and diversity. However, one has to find a balance between allowing some mismatches and avoiding non-specific binding in other genomic regions, which can also be a problem.
Considering the different types of biases (specificity to fungi; mismatches; length; taxonomy), we suggest that different primer combinations targeting different parts of the ITS region should be analyzed in parallel. When dealing with single culture isolates compared to environmental samples, the choice of a primer pair to amplify ITS is less problematic because there is no 'competition' between DNA fragments of different taxonomic groups/lengths, and the DNA quality is generally higher.
This study also illustrates potential benefits of using a bioinformatics approach before selecting primer pairs for a given study. We nevertheless emphasize that an in silico analysis does not necessarily reflect the performance of the primers in vitro, since there are many other PCR parameters such as ITS copy number, amplification program, and salt and primer concentration in the PCR mix that cannot easily be simulated. This study should therefore be followed up by in vitro PCR analyses of the fungal ITS primers where biases are measured based on sequence output, although it will be a huge task to control and check for all types of biases that might be involved. We are currently performing further bioinformatics analyses using the tool 'ecoPrimer' (http://www.grenoble.prabi.fr/trac/ecoPrimers; Riaz et al. unpublished) to identify the most appropriate barcoding primers within the ITS region and other regions, with the intent of determining whether new ITS primers, such as those recently published by Martin and Rygiewicz , should replace the currently used ones.
Eva Bellemain was funded by the Natural History Museum, University of Oslo and this work has been initiated as part of the BarFrost project (Barcoding of permafrost samples). We are thankful to four anonymous reviewers for constructive comments and to Marie Davey for helping to improve the style of written English.
- Anderson I, Cairney J: Diversity and ecology of soil fungal communities: increased understanding through the application of molecular techniques. Environmental Microbiology. 2004, 6 (8): 769-779. 10.1111/j.1462-2920.2004.00675.x.View ArticlePubMedGoogle Scholar
- Chase M, Fay M: Barcoding of plants and fungi. Science. 2009, 325: 682-683. 10.1126/science.1176906.View ArticlePubMedGoogle Scholar
- Horton T, Bruns T: The molecular revolution in ectomycorrhizal ecology: Peeking into the black box. Molecular Ecology. 2001, 10: 1855-1871. 10.1046/j.0962-1083.2001.01333.x.View ArticlePubMedGoogle Scholar
- Seiffert K: Progress toward DNA barcoding of fungi. Molecular Ecology Resources. 2009, 9 (Suppl 1): 83-89. 10.1111/j.1755-0998.2009.02635.x.View ArticleGoogle Scholar
- Freeman K, Martin A, Karki D, Lynch R, Mitter M, Meyer A, Longcore J, Simmons D, Schmidt S: Evidence that chytrids dominate fungal communities in high-elevation soils. Proceeding of the National Academy of Sciences USA. 2009, 106 (43): 18315-18320. 10.1073/pnas.0907303106.View ArticleGoogle Scholar
- Frohlich-Nowoisky J, Pickergill D, Despres V, Poschl U: High diversity of fungi in air particulate matter. Proceeding of the National Academy of Sciences USA. 2009, 106: 12814-12819. 10.1073/pnas.0811003106.View ArticleGoogle Scholar
- Lindahl B, Ihrmark K, Boberg J, Trumbore S, Högberg P, Stenlid J, Finlay R: Spatial separation of litter decomposition and mycorrhizal nitrogen uptake in a boreal forest. New Phytologist. 2007, 173: 611-620. 10.1111/j.1469-8137.2006.01936.x.View ArticlePubMedGoogle Scholar
- O'Brien H, Parrent J, KJackson J, Moncalvo J, Vilgalys R: Fungal community analysis by large-scale sequencing of environmental samples. Applied and Environmental Microbiology. 2005, 71 (9): 5544-5550. 10.1128/AEM.71.9.5544-5550.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Pickles B, Genney D, Potts J, Lennon J, Andersonand I, Alexander I: Spatial and temporal ecology of Scots pine ectomycorrhizas. New Phytologist. 2010, 10.1111/j.1469-8137.2010.03204.xGoogle Scholar
- Zinger L, Coissac E, Choler P, Geremia R: Assessment of microbial communities by graph partitioning in a study of soil fungi in two alpine meadows. Applied and Environmental Microbiology. 2009, 75: 5863-5870. 10.1128/AEM.00748-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Nilsson R, Ryberg M, Abarenkov K, Sjökvist E, Kristiansson E: The ITS region as a target for characterization of fungal communities using emerging sequencing technologies. FEMS Microbiology Letters. 2009, 296: 97-101. 10.1111/j.1574-6968.2009.01618.x.View ArticlePubMedGoogle Scholar
- Nilsson R, Ryberg M, Kristiansson E, Abarenkov K, Larsson K, Koljalg U: Taxonomic reliability of DNA sequences in public sequence databases: a fungal perspective. PLoS One. 2006, 1 (1): e59-10.1371/journal.pone.0000059.PubMed CentralView ArticlePubMedGoogle Scholar
- Vilgalys R, Gonzalez D: Organisation of ribosomal DNA in the basidiomycete Thanatephorus praticola. Current Genetics. 1990, 18: 277-280. 10.1007/BF00318394.View ArticlePubMedGoogle Scholar
- Buée M, Reich M, Murat C, Morin E, Nilsson R, Uroz S, Martin F: 454 Pyrosequencing analyses of forest soils reveal an unexpectedly high fungal diversity. New Phytologist. 2009, 2: 449-456. 10.1111/j.1469-8137.2009.03003.x.View ArticleGoogle Scholar
- Ghannoum M, Jurevic R, Mukherjee P, Cui F, Sikaroodi M, Naqvi A, Gillevet P: Characterization of the Oral Fungal Microbiome (Mycobiome) in Healthy Individuals. PLoS Pathogens. 2010, 6 (1): e1000713-10.1371/journal.ppat.1000713.PubMed CentralView ArticlePubMedGoogle Scholar
- Jumpponen A, Jones K: Massively parallel 454-sequencing of Quercus macrocarpa phyllosphere fungal communities indicates reduced richness and diversity in urban environments. New Phytologist. 2009, 184: 438-448. 10.1111/j.1469-8137.2009.02990.x.View ArticlePubMedGoogle Scholar
- Jumpponen A, Jones K, Mattox J, Yeage C: Massively parallel 454-sequencing of Quercus spp. ectomycorrhizosphere indicates differences in fungal community composition richness, and diversity among urban and rural environments. Molecular Ecology. 2010,Google Scholar
- Gardes M, Bruns T: ITS primers with enhanced specificity for basidiomycetes - application to the identification of mycorrhizae and rusts. Molecular Ecology. 1993, 2 (2): 113-118. 10.1111/j.1365-294X.1993.tb00005.x.View ArticlePubMedGoogle Scholar
- White T, Bruns T, Lee S, Taylor J: Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. PCR-protocols a guide to methods and applications. Edited by: Innis MA, Gelfand DH, Sninski JJ, White TJ. 1990, San Diego: Academic press, 315-322.Google Scholar
- Martin K, Rygiewicz P: Fugal-specific primers developed for analysis of the ITS region of environmental DNA extracts. BMC Microbiology. 2005, 5: 28-10.1186/1471-2180-5-28.PubMed CentralView ArticlePubMedGoogle Scholar
- Kirk PM, Cannon PF, David JC, Stalpers J: Ainsworth and Bisby's Dictionary of the Fungi. 2001, Wallingford UK: CAB International, 9Google Scholar
- Deagle B, Eveson J, Jarman S: Quantification of damage in DNA recovered from highly degraded samples - a case study on DNA in faeces. Frontiers in Zoology. 2006, 3: 11-10.1186/1742-9994-3-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Ficetola GF, Coissac E, Zundel S, Riaz T, Shehzad W, Bessière J, Taberlet P, Pompanon F: An In silico approach for the evaluation of DNA barcodes. BMC Genomics.Google Scholar
- Wu S, Mamber U: Agrep- a fast approximate pattern matching tool. Proceedings of the Winter 1992 USENIX Conference San Francisco USA. Berkeley. 1992, 153-162.Google Scholar
- James T: Reconstructing the early evolution of Fungi using a six-gene phylogeny. Nature. 2006, 443: 818-822. 10.1038/nature05110.View ArticlePubMedGoogle Scholar
- SantaLucia JJ, Hicks D: The thermodynamics of DNA structural motifs. Annual Review of Biophysics and Biomolecular Structure. 2004, 33: 415-440. 10.1146/annurev.biophys.32.110601.141800.View ArticlePubMedGoogle Scholar
- Duitama J, Kumar D, Hemphill E, Khan M, Mandoiu I, Nelson C: Primerhunter: a primer design tool for pcr-based virus subtype identification. Nucleic Acids research. 2009, 37 (8): 2483-2492. 10.1093/nar/gkp073.PubMed CentralView ArticlePubMedGoogle Scholar
- Peay K, Kennedy P, Davies S, Tan S, Bruns T: Potential link between plant and fungal distributions in a dipterocarp rainforest: community and phylogenetic structure of tropical ectomycorrhizal fungi across a plant and soil ecotone. New Phytologist. 2010, 185: 529-542. 10.1111/j.1469-8137.2009.03075.x.View ArticlePubMedGoogle Scholar
- Harris D: Can you bank on GenBank?. Trends in Ecology and Evolution. 2003, 18 (7): 317-319. 10.1016/S0169-5347(03)00150-2.View ArticleGoogle Scholar
- Landeweert R, Leeflang P, Kuyper T, Hoffland E, Rosling A, Wernars K, Smit E: Molecular identification of ectomycorrhizal mycelium in soil horizons. Applied and Environmental Microbiology. 2003, 69 (1): 10.1128/AEM.69.1.327-333.2003. DOI: 10.1128/AEM.1169.1121.1327-1333.2003Google Scholar
- Robinson C, Szaro T, Izzo A, Anderson I, Parkin P, Bruns T: Spatial distribution of fungal communities in a coastal graasland soil. Soil Biology and Biochemistry. 2009, 41: 414-416. 10.1016/j.soilbio.2008.10.021.View ArticleGoogle Scholar
- Hong S, Bunge J, Leslin C, S J, Epstein S: Polymerase chain reaction primers miss half of rRNA microbial diversity. The ISME shopping. 2009, 3: 1365-1373. 10.1038/ismej.2009.89.View ArticleGoogle Scholar
- Jeon S, Bunge J, Leslin C, Stoeck T, Hong S, Epstein S: Environmental rRNA inventories miss over half of protistan diversity. BMC Microbiology. 2008, 8: 222-10.1186/1471-2180-8-222.PubMed CentralView ArticlePubMedGoogle Scholar
- Sipos R, Szekely A, Palatinszky M, Revesz M, K M, Nikolausz M: Effect of primer mismatch annealing temperature and PCR cycle number on 16S rRNA gene -targetting bacterial community analysis. FEMS Microbiology Ecology. 2007, 60: 341-350. 10.1111/j.1574-6941.2007.00283.x.View ArticlePubMedGoogle Scholar
- Engelbrektson A, Kunin V, Wrighton K, Zvenigorodsky N, Chen F, Ochman H, Hugenholtz P: Experimental factors affecting PCR-based estimates of microbial species richness and evenness. The International Society for Microbial Ecology Journal. 2010, doi:10.1038/ismej.2009.153Google Scholar
- Huber J, Morrison H, SM H, Neal P, Sogin M, Welch D: Effect of PCR amplicon size on assessments of clone library microbial diversity and community structure. Environmental Microbiology. 2009, 11 (5): 1292-1302. 10.1111/j.1462-2920.2008.01857.x.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.