Systematic bias introduced by genomic DNA template dilution in 16S rRNA gene-targeted microbiota profiling in human stool homogenates

Francesco Multinu, Sean C. Harrington, Jun Chen, Patricio R. Jeraldo, Stephen Johnson, Nicholas D Chia, Marina Walther-Antonio

Research output: Contribution to journalArticle

6 Citations (Scopus)

Abstract

Variability in representation of microbial communities can be caused by differences in microbial composition or artifacts introduced at sample collection or processing. Alterations in community representation introduced by variations in starting DNA concentrations have not been systematically investigated in stool samples. The goal of this study was to evaluate the effect of the genomic DNA (gDNA) concentration in the resulting 16S rRNA gene library composition and compare its effect to other sample processing variables in homogenized human fecal material. Compared to a gDNA input of 1 ng/μl, inputs of ≤1.6 × 10-3 ng/μl resulted in a marked decrease in the concentration of the 16S rRNA gene amplicon (P < 0.001). Low gDNA concentrations (≤1.6 × 10-3 ng/μl) were also associated with a decrease (P < 0.001) in the number of operational taxonomic units and significant divergence in β-diversity profiles (unweighted UniFrac distance, P < 0.001), as characterized by an overestimation of Proteobacteria and underestimation of Firmicutes. Even a gDNA concentration of 4 × 10-2 ng/μl showed a significant impact on the β-diversity profile (unweighted UniFrac distance, P < 0.03). Overall, the gDNA concentration explained 22.4% to 38.1% of the microbiota variation based on various β-diversity measures (P < 0.001). By comparison, the DNA extraction methods and PCR volumes tested did not significantly affect the microbial composition profile, and the PCR cycling method explained less than 3.7% of the microbiota variation (weighted UniFrac distance, P < 0.03). The 16S rRNA gene yield and the microbial community representation of human homogenized stool samples are significantly altered by gDNA template concentrations of ≤1.6 × 10-3 ng/μl. In addition, data from studies with a gDNA input of ≤4 × 10-2 ng/μl should be interpreted with caution.

Original languageEnglish (US)
Article numbere00560-17
JournalmSphere
Volume3
Issue number2
DOIs
StatePublished - Mar 1 2018

Fingerprint

Microbiota
rRNA Genes
DNA
Proteobacteria
Polymerase Chain Reaction
Gene Library
Artifacts

Keywords

  • 16S rRNA
  • GDNA
  • Microbiome

ASJC Scopus subject areas

  • Microbiology
  • Molecular Biology

Cite this

Systematic bias introduced by genomic DNA template dilution in 16S rRNA gene-targeted microbiota profiling in human stool homogenates. / Multinu, Francesco; Harrington, Sean C.; Chen, Jun; Jeraldo, Patricio R.; Johnson, Stephen; Chia, Nicholas D; Walther-Antonio, Marina.

In: mSphere, Vol. 3, No. 2, e00560-17, 01.03.2018.

Research output: Contribution to journalArticle

@article{316c509cd43244d2afdedb491dbe2546,
title = "Systematic bias introduced by genomic DNA template dilution in 16S rRNA gene-targeted microbiota profiling in human stool homogenates",
abstract = "Variability in representation of microbial communities can be caused by differences in microbial composition or artifacts introduced at sample collection or processing. Alterations in community representation introduced by variations in starting DNA concentrations have not been systematically investigated in stool samples. The goal of this study was to evaluate the effect of the genomic DNA (gDNA) concentration in the resulting 16S rRNA gene library composition and compare its effect to other sample processing variables in homogenized human fecal material. Compared to a gDNA input of 1 ng/μl, inputs of ≤1.6 × 10-3 ng/μl resulted in a marked decrease in the concentration of the 16S rRNA gene amplicon (P < 0.001). Low gDNA concentrations (≤1.6 × 10-3 ng/μl) were also associated with a decrease (P < 0.001) in the number of operational taxonomic units and significant divergence in β-diversity profiles (unweighted UniFrac distance, P < 0.001), as characterized by an overestimation of Proteobacteria and underestimation of Firmicutes. Even a gDNA concentration of 4 × 10-2 ng/μl showed a significant impact on the β-diversity profile (unweighted UniFrac distance, P < 0.03). Overall, the gDNA concentration explained 22.4{\%} to 38.1{\%} of the microbiota variation based on various β-diversity measures (P < 0.001). By comparison, the DNA extraction methods and PCR volumes tested did not significantly affect the microbial composition profile, and the PCR cycling method explained less than 3.7{\%} of the microbiota variation (weighted UniFrac distance, P < 0.03). The 16S rRNA gene yield and the microbial community representation of human homogenized stool samples are significantly altered by gDNA template concentrations of ≤1.6 × 10-3 ng/μl. In addition, data from studies with a gDNA input of ≤4 × 10-2 ng/μl should be interpreted with caution.",
keywords = "16S rRNA, GDNA, Microbiome",
author = "Francesco Multinu and Harrington, {Sean C.} and Jun Chen and Jeraldo, {Patricio R.} and Stephen Johnson and Chia, {Nicholas D} and Marina Walther-Antonio",
year = "2018",
month = "3",
day = "1",
doi = "10.1128/mSphere.00560-17",
language = "English (US)",
volume = "3",
journal = "mSphere",
issn = "2379-5042",
publisher = "American Society for Microbiology",
number = "2",

}

TY - JOUR

T1 - Systematic bias introduced by genomic DNA template dilution in 16S rRNA gene-targeted microbiota profiling in human stool homogenates

AU - Multinu, Francesco

AU - Harrington, Sean C.

AU - Chen, Jun

AU - Jeraldo, Patricio R.

AU - Johnson, Stephen

AU - Chia, Nicholas D

AU - Walther-Antonio, Marina

PY - 2018/3/1

Y1 - 2018/3/1

N2 - Variability in representation of microbial communities can be caused by differences in microbial composition or artifacts introduced at sample collection or processing. Alterations in community representation introduced by variations in starting DNA concentrations have not been systematically investigated in stool samples. The goal of this study was to evaluate the effect of the genomic DNA (gDNA) concentration in the resulting 16S rRNA gene library composition and compare its effect to other sample processing variables in homogenized human fecal material. Compared to a gDNA input of 1 ng/μl, inputs of ≤1.6 × 10-3 ng/μl resulted in a marked decrease in the concentration of the 16S rRNA gene amplicon (P < 0.001). Low gDNA concentrations (≤1.6 × 10-3 ng/μl) were also associated with a decrease (P < 0.001) in the number of operational taxonomic units and significant divergence in β-diversity profiles (unweighted UniFrac distance, P < 0.001), as characterized by an overestimation of Proteobacteria and underestimation of Firmicutes. Even a gDNA concentration of 4 × 10-2 ng/μl showed a significant impact on the β-diversity profile (unweighted UniFrac distance, P < 0.03). Overall, the gDNA concentration explained 22.4% to 38.1% of the microbiota variation based on various β-diversity measures (P < 0.001). By comparison, the DNA extraction methods and PCR volumes tested did not significantly affect the microbial composition profile, and the PCR cycling method explained less than 3.7% of the microbiota variation (weighted UniFrac distance, P < 0.03). The 16S rRNA gene yield and the microbial community representation of human homogenized stool samples are significantly altered by gDNA template concentrations of ≤1.6 × 10-3 ng/μl. In addition, data from studies with a gDNA input of ≤4 × 10-2 ng/μl should be interpreted with caution.

AB - Variability in representation of microbial communities can be caused by differences in microbial composition or artifacts introduced at sample collection or processing. Alterations in community representation introduced by variations in starting DNA concentrations have not been systematically investigated in stool samples. The goal of this study was to evaluate the effect of the genomic DNA (gDNA) concentration in the resulting 16S rRNA gene library composition and compare its effect to other sample processing variables in homogenized human fecal material. Compared to a gDNA input of 1 ng/μl, inputs of ≤1.6 × 10-3 ng/μl resulted in a marked decrease in the concentration of the 16S rRNA gene amplicon (P < 0.001). Low gDNA concentrations (≤1.6 × 10-3 ng/μl) were also associated with a decrease (P < 0.001) in the number of operational taxonomic units and significant divergence in β-diversity profiles (unweighted UniFrac distance, P < 0.001), as characterized by an overestimation of Proteobacteria and underestimation of Firmicutes. Even a gDNA concentration of 4 × 10-2 ng/μl showed a significant impact on the β-diversity profile (unweighted UniFrac distance, P < 0.03). Overall, the gDNA concentration explained 22.4% to 38.1% of the microbiota variation based on various β-diversity measures (P < 0.001). By comparison, the DNA extraction methods and PCR volumes tested did not significantly affect the microbial composition profile, and the PCR cycling method explained less than 3.7% of the microbiota variation (weighted UniFrac distance, P < 0.03). The 16S rRNA gene yield and the microbial community representation of human homogenized stool samples are significantly altered by gDNA template concentrations of ≤1.6 × 10-3 ng/μl. In addition, data from studies with a gDNA input of ≤4 × 10-2 ng/μl should be interpreted with caution.

KW - 16S rRNA

KW - GDNA

KW - Microbiome

UR - http://www.scopus.com/inward/record.url?scp=85045219775&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85045219775&partnerID=8YFLogxK

U2 - 10.1128/mSphere.00560-17

DO - 10.1128/mSphere.00560-17

M3 - Article

C2 - 29564398

AN - SCOPUS:85045219775

VL - 3

JO - mSphere

JF - mSphere

SN - 2379-5042

IS - 2

M1 - e00560-17

ER -