Dataset Information

Stool sampling and DNA isolation kits affect DNA quality and bacterial composition following 16S rRNA gene sequencing using MiSeq Illumina platform.

ABSTRACT: Many studies correlate changes in human gut microbiome with the onset of various diseases, mostly by 16S rRNA gene sequencing. Setting up the optimal sampling and DNA isolation procedures is crucial for robustness and reproducibility of the results. We performed a systematic comparison of several sampling and DNA isolation kits, quantified their effect on bacterial gDNA quality and the bacterial composition estimates at all taxonomic levels. Sixteen volunteers tested three sampling kits. All samples were consequently processed by two DNA isolation kits. We found that the choice of both stool sampling and DNA isolation kits have an effect on bacterial composition with respect to Gram-positivity, however the isolation kit had a stronger effect than the sampling kit. The proportion of bacteria affected by isolation and sampling kits was larger at higher taxa levels compared to lower taxa levels. The PowerLyzer PowerSoil DNA Isolation Kit outperformed the QIAamp DNA Stool Mini Kit mainly due to better lysis of Gram-positive bacteria while keeping the values of all the other assessed parameters within a reasonable range. The presented effects need to be taken into account when comparing results across multiple studies or computing ratios between Gram-positive and Gram-negative bacteria.

SUBMITTER: Videnska P

PROVIDER: S-EPMC6761292 | biostudies-literature | 2019 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Stool sampling and DNA isolation kits affect DNA quality and bacterial composition following 16S rRNA gene sequencing using MiSeq Illumina platform.

Videnska Petra P Smerkova Kristyna K Zwinsova Barbora B Popovici Vlad V Micenkova Lenka L Sedlar Karel K Budinska Eva E

Scientific reports 20190925 1

Many studies correlate changes in human gut microbiome with the onset of various diseases, mostly by 16S rRNA gene sequencing. Setting up the optimal sampling and DNA isolation procedures is crucial for robustness and reproducibility of the results. We performed a systematic comparison of several sampling and DNA isolation kits, quantified their effect on bacterial gDNA quality and the bacterial composition estimates at all taxonomic levels. Sixteen volunteers tested three sampling kits. All sam ...[more]

PMID: 31554833

Similar Datasets

Project description:In this meta-analysis, 17 rumen epithelial 16S rRNA gene Illumina MiSeq amplicon sequencing data sets were analyzed to identify a core rumen epithelial microbiota and core rumen epithelial OTUs shared between the different studies included. Sequences were quality-filtered and screened for chimeric sequences before performing closed-reference 97% OTU clustering, and de novo 97% OTU clustering. Closed-reference OTU clustering identified the core rumen epithelial OTUs, defined as any OTU present in ≥ 80% of the samples, while the de novo data was randomly subsampled to 10,000 reads per sample to generate phylum- and genus-level distributions and beta diversity metrics. 57 core rumen epithelial OTUs were identified including metabolically important taxa such as Ruminococcus, Butyrivibrio, and other Lachnospiraceae, as well as sulfate-reducing bacteria Desulfobulbus and Desulfovibrio. Two Betaproteobacteria OTUs (Neisseriaceae and Burkholderiaceae) were core rumen epithelial OTUs, in contrast to rumen content where previous literature indicates they are rarely found. Two core OTUs were identified as the methanogenic archaea Methanobrevibacter and Methanomethylophilaceae. These core OTUs are consistently present across the many variables between studies which include different host species, geographic region, diet, age, farm management practice, time of year, hypervariable region sequenced, and more. When considering only cattle samples, the number of core rumen epithelial OTUs expands to 147, highlighting the increased similarity within host species despite geographical location and other variables. De novo OTU clustering revealed highly similar rumen epithelial communities, predominated by Firmicutes, Bacteroidetes, and Proteobacteria at the phylum level which comprised 79.7% of subsampled sequences. The 15 most abundant genera represented an average of 54.5% of sequences in each individual study. These abundant taxa broadly overlap with the core rumen epithelial OTUs, with the exception of Prevotellaceae which were abundant, but not identified within the core OTUs. Our results describe the core and abundant bacteria found in the rumen epithelial environment and will serve as a basis to better understand the composition and function of rumen epithelial communities.

Project description:As new sequencing technologies become cheaper and older ones disappear, laboratories switch vendors and platforms. Validating the new setups is a crucial part of conducting rigorous scientific research. Here we report on the reliability and biases of performing bacterial 16S rRNA gene amplicon paired-end sequencing on the MiSeq Illumina platform. We designed a protocol using 50 barcode pairs to run samples in parallel and coded a pipeline to process the data. Sequencing the same sediment sample in 248 replicates as well as 70 samples from alkaline soda lakes, we evaluated the performance of the method with regards to estimates of alpha and beta diversity. Using different purification and DNA quantification procedures we always found up to 5-fold differences in the yield of sequences between individually barcodes samples. Using either a one-step or a two-step PCR preparation resulted in significantly different estimates in both alpha and beta diversity. Comparing with a previous method based on 454 pyrosequencing, we found that our Illumina protocol performed in a similar manner - with the exception for evenness estimates where correspondence between the methods was low. We further quantified the data loss at every processing step eventually accumulating to 50% of the raw reads. When evaluating different OTU clustering methods, we observed a stark contrast between the results of QIIME with default settings and the more recent UPARSE algorithm when it comes to the number of OTUs generated. Still, overall trends in alpha and beta diversity corresponded highly using both clustering methods. Our procedure performed well considering the precisions of alpha and beta diversity estimates, with insignificant effects of individual barcodes. Comparative analyses suggest that 454 and Illumina sequence data can be combined if the same PCR protocol and bioinformatic workflows are used for describing patterns in richness, beta-diversity and taxonomic composition.

Project description:Illumina's MiSeq has become the dominant platform for gene amplicon sequencing in microbial ecology studies; however, various technical concerns, such as reproducibility, still exist. To assess reproducibility, 16S rRNA gene amplicons from 18 soil samples of a reciprocal transplantation experiment were sequenced on an Illumina MiSeq. The V4 region of 16S rRNA gene from each sample was sequenced in triplicate with each replicate having a unique barcode. The average OTU overlap, without considering sequence abundance, at a rarefaction level of 10,323 sequences was 33.4±2.1% and 20.2±1.7% between two and among three technical replicates, respectively. When OTU sequence abundance was considered, the average sequence abundance weighted OTU overlap was 85.6±1.6% and 81.2±2.1% for two and three replicates, respectively. Removing singletons significantly increased the overlap for both (~1-3%, p<0.001). Increasing the sequencing depth to 160,000 reads by deep sequencing increased OTU overlap both when sequence abundance was considered (95%) and when not (44%). However, if singletons were not removed the overlap between two technical replicates (not considering sequence abundance) plateaus at 39% with 30,000 sequences. Diversity measures were not affected by the low overlap as α-diversities were similar among technical replicates while β-diversities (Bray-Curtis) were much smaller among technical replicates than among treatment replicates (e.g., 0.269 vs. 0.374). Higher diversity coverage, but lower OTU overlap, was observed when replicates were sequenced in separate runs. Detrended correspondence analysis indicated that while there was considerable variation among technical replicates, the reproducibility was sufficient for detecting treatment effects for the samples examined. These results suggest that although there is variation among technical replicates, amplicon sequencing on MiSeq is useful for analyzing microbial community structure if used appropriately and with caution. For example, including technical replicates, removing spurious sequences and unrepresentative OTUs, using a clustering method with a high stringency for OTU generation, estimating treatment effects at higher taxonomic levels, and adapting the unique molecular identifier (UMI) and other newly developed methods to lower PCR and sequencing error and to identify true low abundance rare species all can increase reproducibility.

Project description:BackgroundAdvances in sequencing technologies and bioinformatics have made the analysis of microbial communities almost routine. Nonetheless, the need remains to improve on the techniques used for gathering such data, including increasing throughput while lowering cost and benchmarking the techniques so that potential sources of bias can be better characterized.MethodsWe present a triple-index amplicon sequencing strategy to sequence large numbers of samples at significantly lower c ost and in a shorter timeframe compared to existing methods. The design employs a two-stage PCR protocol, incorpo rating three barcodes to each sample, with the possibility to add a fourth-index. It also includes heterogeneity spacers to overcome low complexity issues faced when sequencing amplicons on Illumina platforms.ResultsThe library preparation method was extensively benchmarked through analysis of a mock community in order to assess biases introduced by sample indexing, number of PCR cycles, and template concentration. We further evaluated the method through re-sequencing of a standardized environmental sample. Finally, we evaluated our protocol on a set of fecal samples from a small cohort of healthy adults, demonstrating good performance in a realistic experimental setting. Between-sample variation was mainly related to batch effects, such as DNA extraction, while sample indexing was also a significant source of bias. PCR cycle number strongly influenced chimera formation and affected relative abundance estimates of species with high GC content. Libraries were sequenced using the Illumina HiSeq and MiSeq platforms to demonstrate that this protocol is highly scalable to sequence thousands of samples at a very low cost.ConclusionsHere, we provide the most comprehensive study of performance and bias inherent to a 16S rRNA gene amplicon sequencing method to date. Triple-indexing greatly reduces the number of long custom DNA oligos required for library preparation, while the inclusion of variable length heterogeneity spacers minimizes the need for PhiX spike-in. This design results in a significant cost reduction of highly multiplexed amplicon sequencing. The biases we characterize highlight the need for highly standardized protocols. Reassuringly, we find that the biological signal is a far stronger structuring factor than the various sources of bias.

Dataset Information

Stool sampling and DNA isolation kits affect DNA quality and bacterial composition following 16S rRNA gene sequencing using MiSeq Illumina platform.

Publications

Stool sampling and DNA isolation kits affect DNA quality and bacterial composition following 16S rRNA gene sequencing using MiSeq Illumina platform.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets