Dataset Information

High-resolution microbial community reconstruction by integrating short reads from multiple 16S rRNA regions.

ABSTRACT: The emergence of massively parallel sequencing technology has revolutionized microbial profiling, allowing the unprecedented comparison of microbial diversity across time and space in a wide range of host-associated and environmental ecosystems. Although the high-throughput nature of such methods enables the detection of low-frequency bacteria, these advances come at the cost of sequencing read length, limiting the phylogenetic resolution possible by current methods. Here, we present a generic approach for integrating short reads from large genomic regions, thus enabling phylogenetic resolution far exceeding current methods. The approach is based on a mapping to a statistical model that is later solved as a constrained optimization problem. We demonstrate the utility of this method by analyzing human saliva and Drosophila samples, using Illumina single-end sequencing of a 750 bp amplicon of the 16S rRNA gene. Phylogenetic resolution is significantly extended while reducing the number of falsely detected bacteria, as compared with standard single-region Roche 454 Pyrosequencing. Our approach can be seamlessly applied to simultaneous sequencing of multiple genes providing a higher resolution view of the composition and activity of complex microbial communities.

SUBMITTER: Amir A

PROVIDER: S-EPMC3905898 | biostudies-literature | 2013 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

High-resolution microbial community reconstruction by integrating short reads from multiple 16S rRNA regions.

Amir Amnon A Zeisel Amit A Zuk Or O Elgart Michael M Stern Shay S Shamir Ohad O Turnbaugh Peter J PJ Soen Yoav Y Shental Noam N

Nucleic acids research 20131107 22

The emergence of massively parallel sequencing technology has revolutionized microbial profiling, allowing the unprecedented comparison of microbial diversity across time and space in a wide range of host-associated and environmental ecosystems. Although the high-throughput nature of such methods enables the detection of low-frequency bacteria, these advances come at the cost of sequencing read length, limiting the phylogenetic resolution possible by current methods. Here, we present a generic a ...[more]

PMID: 24214960

Similar Datasets

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

High-resolution microbial community reconstruction by integrating short reads from multiple 16S rRNA regions.

Publications

High-resolution microbial community reconstruction by integrating short reads from multiple 16S rRNA regions.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure