Dataset Information

The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses.

ABSTRACT: 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.

SUBMITTER: Vetrovsky T

PROVIDER: S-EPMC3583900 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:A curated murine oral microbiome database to be used as a reference for mouse-based studies has been constructed using a combination of bacterial culture, 16S rRNA gene amplicon, and whole-genome sequencing. The database comprises a collection of nearly full-length 16S rRNA gene sequences from cultured isolates and draft genomes from representative taxa collected from a range of sources, including specific-pathogen-free laboratory mice, wild Mus musculus domesticus mice, and formerly wild wood mouse Apodemus sylvaticus At present, it comprises 103 mouse oral taxa (MOT) spanning four phyla-Firmicutes, Proteobacteria, Actinobacteria, and Bacteroidetes-including 12 novel undescribed species-level taxa. The key observations from this study are (i) the low diversity and predominantly culturable nature of the laboratory mouse oral microbiome and (ii) the identification of three major murine-specific oral bacterial lineages, namely, Streptococcus danieliae (MOT10), Lactobacillus murinus (MOT93), and Gemella species 2 (MOT43), which is one of the novel, still-unnamed taxa. Of these, S. danieliae is of particular interest, since it is a major component of the oral microbiome from all strains of healthy and periodontally diseased laboratory mice, as well as being present in wild mice. It is expected that this well-characterized database should be a useful resource for in vitro experimentation and mouse model studies in the field of oral microbiology.IMPORTANCE Mouse model studies are frequently used in oral microbiome research, particularly to investigate diseases such as periodontitis and caries, as well as other related systemic diseases. We have reported here the details of the development of a curated reference database to characterize the oral microbial community in laboratory and some wild mice. The genomic information and findings reported here can help improve the outcomes and accuracy of host-microbe experimental studies that use murine models to understand health and disease. Work is also under way to make the reference data sets publicly available on a web server to enable easy access and downloading for researchers across the world.

Project description:BackgroundAccurate identification of bacterial communities is crucial for research applications, diagnostics, and clinical interventions. Although 16S ribosomal RNA (rRNA) gene sequencing is a widely employed technique for bacterial taxonomic classification, it often results in misclassified or unclassified bacterial taxa. This study sought to refine the full-length 16S rRNA gene sequencing protocol using the MinION sequencer, focusing on the V1-V9 regions. Our methodological enquiry examined several factors, including the number of PCR amplification cycles, choice of primers and Taq polymerase, and specific sequence databases and workflows employed. We used a microbial standard comprising eight bacterial strains (five gram-positive and three gram-negative) in known proportions as a validation control.ResultsBased on the MinION protocol, we employed the microbial standard as the DNA template for the 16S rRNA gene amplicon sequencing procedure. Our analysis showed that an elevated number of PCR amplification cycles introduced PCR bias, and the selection of Taq polymerase and primer sets significantly affected the subsequent analysis. Bacterial identification at genus level demonstrated Pearson correlation coefficients ranging from 0.73 to 0.79 when assessed using BugSeq, Kraken-Silva and EPI2ME-16S workflows. Notably, the EPI2ME-16S workflow exhibited the highest Pearson correlation with the microbial standard, minimised misclassification, and increased alignment accuracy. At the species taxonomic level, the BugSeq workflow was superior, with a Pearson correlation coefficient of 0.92.ConclusionsThese findings emphasise the importance of careful selection of PCR settings and a well-structured analytical framework for 16S rRNA full-length gene sequencing. The results showed a robust correlation between the predicted and observed bacterial abundances at both the genus and species taxonomic levels, making these findings applicable across diverse research contexts and with clinical utility for reliable pathogen identification.

Project description:Two consortia (Consortium A and Consortium B) that can use 1,4-dioxane (a groundwater contaminant of emerging concern) as the sole carbon source were enriched from Rice University (Houston, TX, USA) campus soil. Phylogenetic analysis by 16S rRNA sequencing revealed the dominant genus in both of the consortia is Mycobacterium (56% in Consortium A and 49% in Consortium B). The predominance of Mycobacterium spp, in these consortia support the notion that this is an important and commonly encountered genus of dioxane degraders. Among other genera present that make at least 2% of these consortia, only Afipia encompasses a strain (i.e., Afipia sp. D1) that was reported to degrade dioxane as sole carbon and energy source. A nested PCR analysis using two degenerate primers to target the hydroxylase alpha subunit of groups 3 to 6 SDIMOs was performed to gain insights into which enzymes were responsible for dioxane degradation by these consortia. The purified products obtained from the second PCR run were sequenced and compared to genes databases (NCBI) encompassing all of the currently reported SDIMOs. The dominant SDIMO genes in Consortium A corresponded to a group-6 putative propane monooxygenase-like SDIMO (98.8%); while in Consortium B, SDIMO genes from both groups 5 (47.3%) and 6 (51.9%) were observed. In both consortia, the relative abundance of thmA/dxmA gene was negligible (0.03%), which is consistent with the negative amplification of these genes as verified in qPCR. Overall, the high relative abundance of group-6 putative propane monooxygenases in our two consortia suggests the novel finding that group 6-SDIMOs could also play an important role in dioxane degradation. This underscores the need for further research on genes and enzymes involved in dioxane biodegradation to develop novel biomarkers that can be useful for forensic analysis and performance assessment of bioremediation and natural attenuation at dioxane-impacted sites. DNA was extracted from bacteria biomass harvested in exponential growth phase, when half or more of the added dioxane (100 mg/L) was consumed. Total DNA extractions were performed using the UltraClean® Microbial DNA Isolation Kit (MO BIO, Carlsbad, CA, USA) according to the manufacturer’s protocol. The V4 region of the 16S rRNA gene was amplified by PCR using the forward 515F and reverse 806R primers. Sequencing was performed at MR DNA (www.mrdnalab.com, Shallowater, TX, USA) by Illumina MiSeq paired-end sequencing (approximately 2×300 bp as the read length). Sequence data were processed using MR DNA analysis pipeline. Operational taxonomic units (OTUs) were defined by clustering at 3% divergence (97% similarity). Final OTUs were taxonomically classified using BLASTn against the RDPII (http://rdp.cme.msu.edu) and NCBI (www.ncbi.nlm.nih.gov) databases.Previously designed degenerate primers NVC57, NVC58, NVC65 and NVC66 to target conserved regions in the soluble di-iron monooxygenases (SDIMO) alpha subunit gene (groups 3 to 6) were used to examine the presence and diversity of SDIMO genes in these two consortia. A nested PCR strategy was used to increase the PCR product yield. In the first run, the PCR mixture contained 1 µL of NVC65 and NVC58 primer mixture (10 µM), 20 ng of the extracted genomic DNA, 12.5 µL of KAPA HiFi HotStart ReadyMix (2X) (KAPA Biosystems, Wilmington, MA, USA), and nuclease-free water to yield a total volume of 25 µL. PCR was performed in a Bio-Rad Thermal Cycler (Bio-Rad, Hercules, CA, USA) with the following temperature profile: initial denaturation (94°C, 5 min), then 29 amplification cycles (94°C for 30 s, 55°C for 30 s, 72°C for 1 min per kb) and a final extension (72°C for 5 min). The length of the PCR products in the first run was checked by 1% agarose gel and DNA bands of the correct size (1100 bp) were excised and purified. 20 ng of the purified PCR product was used as the DNA template in the second run, with the second set of primers (NVC57 and NVC66). The purified product (420 bp) from the second PCR was sent to MR DNA (www.mrdnalab.com, Shallowater, TX, USA) for Illumina MiSeq paired-end sequencing (approximately 2×300 bp as the read length). Sequence data were processed using MR DNA analysis pipeline. Operational taxonomic units (OTUs) were defined by clustering at 3% divergence (97% similarity). A database including all of the currently reported SDIMO genes on NCBI was created and used to taxonomically classify the final OTUs.

Dataset Information

The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets