16S ribosomal RNA currently represents the most important target of study in bacterial ecology. Its use for the description of bacterial diversity is, however, limited by the presence of variable copy numbers in bacterial genomes and sequence variation within closely related taxa or within a genome. Here we use the information from sequenced bacterial genomes to explore the variability of 16S rRNA sequences and copy numbers at various taxonomic levels and apply it to estimate bacterial genome and DNA abundances. In total, 7,081 16S rRNA sequences were in silico extracted from 1,690 available bacterial genomes (1-15 per genome). While there are several phyla containing low 16S rRNA copy numbers, in certain taxa, e.g., the Firmicutes and Gammaproteobacteria, the variation is large. Genome sizes are more conserved at all tested taxonomic levels than 16S rRNA copy numbers. Only a minority of bacterial genomes harbors identical 16S rRNA gene copies, and sequence diversity increases with increasing copy numbers. While certain taxa harbor dissimilar 16S rRNA genes, others contain sequences common to multiple species. Sequence identity clusters (often termed operational taxonomic units) thus provide an imperfect representation of bacterial taxa of a certain phylogenetic rank. We have demonstrated that the information on 16S rRNA copy numbers and genome sizes of genome-sequenced bacteria may be used as an estimate for the closest related taxon in an environmental dataset to calculate alternative estimates of the relative abundance of individual bacterial taxa in environmental samples. Using an example from forest soil, this procedure would increase the abundance estimates of Acidobacteria and decrease these of Firmicutes. Using the currently available information, alternative estimates of bacterial community composition may be obtained in this way if the variation of 16S rRNA copy numbers among bacteria is considered.

Laboratory of Environmental Microbiology Institute of Microbiology of the Academy of Sciences of the Czech Republic Praha Czech Republic

See more in PubMed

Hugenholtz P, Goebel BM, Pace NR (1998) Impact of culture-independent studies on the emerging phylogenetic view of bacterial diversity. J Bacteriol 180: 4765–4774. PubMed PMC

Daubin V, Moran NA, Ochman H (2003) Phylogenetics and the cohesion of bacterial genomes. Science 301: 829–832. PubMed

Head IM, Saunders JR, Pickup RW (1998) Microbial evolution, diversity, and ecology: A decade of ribosomal RNA analysis of uncultivated microorganisms. Microb Ecol 35: 1–21. PubMed

Klappenbach JA, Saxman PR, Cole JR, Schmidt TM (2001) rrndb: the Ribosomal RNA Operon Copy Number Database. Nuc Acid Res 29: 181–184. PubMed PMC

Acinas SG, Marcelino LA, Klepac-Ceraj V, Polz MF (2004) Divergence and redundancy of 16S rRNA sequences in genomes with multiple rrn operons. J Bacteriol 186: 2629–2635. PubMed PMC

Eichorst SA, Breznak JA, Schmidt TM (2007) Isolation and Characterization of Soil Bacteria That Define Terriglobus gen. nov., in the Phylum Acidobacteria. Appl Environ Microbiol 73: 2708–2717. PubMed PMC

Klappenbach JA, Dunbar JM, Schmidt TM (2000) rRNA Operon Copy Number Reflects Ecological Strategies of Bacteria. Appl Environ Microbiol 66: 1328–1333. PubMed PMC

Hashimoto JQ, Stevenson BS, Schmidt TM (2003) Rates and consequences of recombination between rRNA operons. J Bacteriol 185: 966–972. PubMed PMC

Wang Y, Zwang ZS, Ramanan N (1997) The actinomycete Thermobispora bispora contains two distinct types of transcriptionally active 16S rRNA genes. J Bacteriol 179: 3270–3276. PubMed PMC

Yap WH, Zhang ZS, Wang Y (1999) Distinct types of rRNA operons exist in the genome of the actinomycete Thermomonospora chromogena and evidence for horizontal transfer of an entire rRNA operon. J Bacteriol 181: 5201–5209. PubMed PMC

Pei AY, Oberdorf WE, Nossa CW, Agarwal A, Chokshi P, et al. (2010) Diversity of 16S rRNA Genes within Individual Prokaryotic Genomes. Appl Environ Microbiol 76: 3886–3897. PubMed PMC

Farrelly V, Rainey FA, Stackebrandt E (1995) Effect of Genome Size and rrn Gene Copy Number on PCR Amplification of 16S Ribosomal RNA Genes from a Mixture of Bacterial Species. Appl Environ Microbiol 61: 2798–2801. PubMed PMC

Amend AS, Seifert KA, Bruns TD (2010) Quantifying microbial communities with 454 pyrosequencing: does read abundance count? Mol Ecol 19: 5555–5565. PubMed

Baldrian P, Větrovský T, Cajthaml T, Dobiášová P, Petránková M, et al. (2013) Estimation of fungal biomass in forest litter and soil. Fungal Ecol in press doi:10.1016/j.funeco.2012.10.002. DOI

Baldrian P, Kolarik M, Stursova M, Kopecky J, Valaskova V, et al. (2012) Active and total microbial communities in forest soil are largely different and highly stratified during decomposition. ISME J 6: 248–258. PubMed PMC

Dowd S, Callaway T, Wolcott R, Sun Y, McKeehan T, et al. (2008) Evaluation of the bacterial diversity in the feces of cattle using 16S rDNA bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP). BMC Microbiol 8: 125. PubMed PMC

Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol 28: 2731–2739. PubMed PMC

Li W, Godzik A (2006) Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22: 1658–1659. PubMed

Reeder J, Knight R (2010) Rapidly denoising pyrosequencing amplicon reads by exploiting rank-abundance distributions. Nat Meth 7: 668–669. PubMed PMC

Edgar RC (2010) Search and clustering orders of magnitude faster than BLAST. Bioinformatics btq461. PubMed

Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, et al. (1997) Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nuc Acid Res 25: 3389–3402. PubMed PMC

Candela M, Vitali B, Matteuzzi D, Brigidi P (2004) Evaluation of the rrn operon copy number in Bifidobacterium using real-time PCR. Lett Appl Microbiol 38: 229–232. PubMed

Lee CM, Sieo CC, Abdullah N, Ho YW (2008) Estimation of 16S rRNA gene copy number in several probiotic Lactobacillus strains isolated from the gastrointestinal tract of chicken. FEMS Microbiol Lett 287: 136–141. PubMed PMC

Fox GE, Wisotzkey JD, Jurtshuk P (1992) How close is close: 16S rRNA sequence identity may not be sufficient to guarantee species identity. Int J Systematic Bacteriol 42: 166–170. PubMed

Stackebrandt E, Goebel BM (1994) A Place for DNA-DNA Reassociation and 16S Ribosomal RNA Sequence Analysis in the Present Species Definition in Bacteriology. Int J Syst Bacteriol 44: 846–849.

Lauber CL, Hamady M, Knight R, Fierer N (2009) Pyrosequencing-Based Assessment of Soil pH as a Predictor of Soil Bacterial Community Structure at the Continental Scale. Appl Environ Microbiol 75: 5111–5120. PubMed PMC

Sukenik A, Kaplan-Levy RN, Welch JM, Post AF (2012) Massive multiplication of genome and ribosomes in dormant cells (akinetes) of Aphanizomenon ovalisporum (Cyanobacteria). ISME J 6: 670–679. PubMed PMC

Ohtani N, Tomita M, Itaya M (2010) An Extreme Thermophile, Thermus thermophilus, Is a Polyploid Bacterium. J Bacteriol 192: 5499–5505. PubMed PMC

Bresler V, Fishelson L (2003) Polyploidy and polyteny in the gigantic eubacterium Epulopiscium fishelsoni. Marine Biol 143: 17–21.

Dahllof I, Baillie H, Kjelleberg S (2000) rpoB-based microbial community analysis avoids limitations inherent in 16S rRNA gene intraspecies heterogeneity. Appl Environ Microbiol 66: 3376–3380. PubMed PMC

Santos SR, Ochman H (2004) Identification and phylogenetic sorting of bacterial lineages with universally conserved genes and proteins. Environ Microbiol 6: 754–759. PubMed

Konstantinidis KT, Ramette A, Tiedje JM (2006) Toward a more robust assessment of intraspecies diversity, using fewer genetic markers. Appl Environ Microbiol 72: 7286–7293. PubMed PMC

von Mering C, Hugenholtz P, Raes J, Tringe SG, Doerks T, et al. (2007) Quantitative phylogenetic assessment of microbial communities in diverse environments. Science 315: 1126–1130. PubMed

Schellenberg J, Links MG, Hill JE, Hemmingsen SM, Peters GA, et al... (2011) Pyrosequencing of Chaperonin-60 (cpn60) Amplicons as a Means of Determining Microbial Community Composition. In: Kwon YM, Ricke SC, editors. High-Throughput Next Generation Sequencing: Methods and Application. Totowa: Humana Press Inc. pp. 143–158. PubMed

Case RJ, Boucher Y, Dahllof I, Holmstrom C, Doolittle WF, et al. (2007) Use of 16S rRNA and rpoB genes as molecular markers for microbial ecology studies. Appl Environ Microbiol 73: 278–288. PubMed PMC

Roux S, Enault F, Bronner G, Debroas D (2011) Comparison of 16S rRNA and protein-coding genes as molecular markers for assessing microbial diversity (Bacteria and Archaea) in ecosystems. FEMS Microbiol Ecol 78: 617–628. PubMed

Vos M, Quince C, Pijl AS, de Hollander M, Kowalchuk GA (2012) A Comparison of rpoB and 16S rRNA as Markers in Pyrosequencing Studies of Bacterial Diversity. PLOS One 7: e30600. PubMed PMC

Host traits rather than migration and molting strategies explain feather bacterial load in Palearctic passerines

Impact of Wood Age on Termite Microbial Assemblages

The functional evolution of termite gut microbiota

Comparison of Actinobacteria communities from human-impacted and pristine karst caves

Determination of Butyrate Synthesis Capacity in Gut Microbiota: Quantification of but Gene Abundance by qPCR in Fecal Samples

Ecological Divergence Within the Enterobacterial Genus Sodalis: From Insect Symbionts to Inhabitants of Decomposing Deadwood

The Oral Microbiome in Periodontal Health

R/G Value-A Numeric Index of Individual Periodontal Health and Oral Microbiome Dynamics

Explorative Meta-Analysis of 417 Extant Archaeal Genomes to Predict Their Contribution to the Total Microbiome Functionality

16S rRNA Gene Copy Number Normalization Does Not Provide More Reliable Conclusions in Metataxonomic Surveys

Termites Are Associated with External Species-Specific Bacterial Communities

Characteristics of microbial community of soil subjected to industrial production of antibiotics

Phenotypic and Genomic Analysis of Clostridium beijerinckii NRRL B-598 Mutants With Increased Butanol Tolerance

Fragment of the aspartyl-tRNA synthetase applicable as a shared classification and phylogenetic marker in particular representatives of the order Lactobacillales

The concept of operational taxonomic units revisited: genomes of bacteria that are regarded as closely related are often highly dissimilar

Suitability and setup of next-generation sequencing-based method for taxonomic characterization of aquatic microbial biofilm

The threonine-tRNA ligase gene region is applicable in classification, typing, and phylogenetic analysis of bifidobacteria

Seasonal and Sexual Differences in the Microbiota of the Hoopoe Uropygial Secretion

Gene encoding the CTP synthetase as an appropriate molecular tool for identification and phylogenetic study of the family Bifidobacteriaceae

Mucin Cross-Feeding of Infant Bifidobacteria and Eubacterium hallii