Dataset Information

Advanced computational algorithms for microbial community analysis using massive 16S rRNA sequence data.

ABSTRACT: With the aid of next-generation sequencing technology, researchers can now obtain millions of microbial signature sequences for diverse applications ranging from human epidemiological studies to global ocean surveys. The development of advanced computational strategies to maximally extract pertinent information from massive nucleotide data has become a major focus of the bioinformatics community. Here, we describe a novel analytical strategy including discriminant and topology analyses that enables researchers to deeply investigate the hidden world of microbial communities, far beyond basic microbial diversity estimation. We demonstrate the utility of our approach through a computational study performed on a previously published massive human gut 16S rRNA data set. The application of discriminant and topology analyses enabled us to derive quantitative disease-associated microbial signatures and describe microbial community structure in far more detail than previously achievable. Our approach provides rigorous statistical tools for sequence-based studies aimed at elucidating associations between known or unknown organisms and a variety of physiological or environmental conditions.

SUBMITTER: Sun Y

PROVIDER: S-EPMC3001099 | biostudies-literature | 2010 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Advanced computational algorithms for microbial community analysis using massive 16S rRNA sequence data.

Sun Yijun Y Cai Yunpeng Y Mai Volker V Farmerie William W Yu Fahong F Li Jian J Goodison Steve S

Nucleic acids research 20101006 22

With the aid of next-generation sequencing technology, researchers can now obtain millions of microbial signature sequences for diverse applications ranging from human epidemiological studies to global ocean surveys. The development of advanced computational strategies to maximally extract pertinent information from massive nucleotide data has become a major focus of the bioinformatics community. Here, we describe a novel analytical strategy including discriminant and topology analyses that enab ...[more]

PMID: 20929878

Similar Datasets

Project description:BackgroundSequencing of the PCR-amplified 16S rRNA gene has become a common approach to microbial community investigations in the fields of human health and environmental sciences. This approach, however, is difficult when the amount of DNA is too low to be amplified by standard PCR. Nested PCR can be employed as it can amplify samples with DNA concentration several-fold lower than standard PCR. However, potential biases with nested PCRs that could affect measurement of community structure have received little attention.ResultsIn this study, we used 17 DNAs extracted from vaginal swabs and 12 DNAs extracted from stool samples to study the influence of nested PCR amplification of the 16S rRNA gene on the estimation of microbial community structure using Illumina MiSeq sequencing. Nested and standard PCR methods were compared on alpha- and beta-diversity metrics and relative abundances of bacterial genera. The effects of number of cycles in the first round of PCR (10 vs. 20) and microbial diversity (relatively low in vagina vs. high in stool) were also investigated. Vaginal swab samples showed no significant difference in alpha diversity or community structure between nested PCR and standard PCR (one round of 40 cycles). Stool samples showed significant differences in alpha diversity (except Shannon's index) and relative abundance of 13 genera between nested PCR with 20 cycles in the first round and standard PCR (P<0.01), but not between nested PCR with 10 cycles in the first round and standard PCR. Operational taxonomic units (OTUs) that had low relative abundance (sum of relative abundance <0.167) accounted for most of the distortion (>27% of total OTUs in stool).ConclusionsNested PCR introduced bias in estimated diversity and community structure. The bias was more significant for communities with relatively higher diversity and when more cycles were applied in the first round of PCR. We conclude that nested PCR could be used when standard PCR does not work. However, rare taxa detected by nested PCR should be validated by other technologies.

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.

Project description:BackgroundFor decades, 16S ribosomal RNA sequencing has been the primary means for identifying the bacterial species present in a sample with unknown composition. One of the most widely used tools for this purpose today is the QIIME (Quantitative Insights Into Microbial Ecology) package. Recent results have shown that the newest release, QIIME 2, has higher accuracy than QIIME, MAPseq, and mothur when classifying bacterial genera from simulated human gut, ocean, and soil metagenomes, although QIIME 2 also proved to be the most computationally expensive. Kraken, first released in 2014, has been shown to provide exceptionally fast and accurate classification for shotgun metagenomics sequencing projects. Bracken, released in 2016, then provided users with the ability to accurately estimate species or genus relative abundances using Kraken classification results. Kraken 2, which matches the accuracy and speed of Kraken 1, now supports 16S rRNA databases, allowing for direct comparisons to QIIME and similar systems.MethodsFor a comprehensive assessment of each tool, we compare the computational resources and speed of QIIME 2's q2-feature-classifier, Kraken 2, and Bracken in generating the three main 16S rRNA databases: Greengenes, SILVA, and RDP. For an evaluation of accuracy, we evaluated each tool using the same simulated 16S rRNA reads from human gut, ocean, and soil metagenomes that were previously used to compare QIIME, MAPseq, mothur, and QIIME 2. We evaluated accuracy based on the accuracy of the final genera read counts assigned by each tool. Finally, as Kraken 2 is the only tool providing per-read taxonomic assignments, we evaluate the sensitivity and precision of Kraken 2's per-read classifications.ResultsFor both the Greengenes and SILVA database, Kraken 2 and Bracken are up to 100 times faster at database generation. For classification, using the same data as previous studies, Kraken 2 and Bracken are up to 300 times faster, use 100x less RAM, and generate results that more accurate at 16S rRNA profiling than QIIME 2's q2-feature-classifier.ConclusionKraken 2 and Bracken provide a very fast, efficient, and accurate solution for 16S rRNA metataxonomic data analysis. Video Abstract.

Dataset Information

Advanced computational algorithms for microbial community analysis using massive 16S rRNA sequence data.

Publications

Advanced computational algorithms for microbial community analysis using massive 16S rRNA sequence data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets