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 Chia; Marina R Walther-Antonio

doi:10.1128/mSphere.00560-17

Systematic Bias Introduced by Genomic DNA Template Dilution in 16S rRNA Gene-Targeted Microbiota Profiling in Human Stool Homogenates

mSphere. 2018 Mar 14;3(2):e00560-17. doi: 10.1128/mSphere.00560-17. eCollection 2018 Mar-Apr.

Authors

Francesco Multinu¹, Sean C Harrington^#², Jun Chen^#^{3

2}, Patricio R Jeraldo², Stephen Johnson³, Nicholas Chia^{2

4

5}, Marina R Walther-Antonio^{1

2

5}

Affiliations

¹ Department of Obstetrics and Gynecology, Mayo Clinic, Rochester, Minnesota, USA.
² Center for Individualized Medicine, Mayo Clinic, Rochester, Minnesota, USA.
³ Division of Biomedical Statistics and Informatics, Mayo Clinic, Rochester, Minnesota, USA.
⁴ Department of Physiology and Biomedical Engineering, Mayo Clinic, Rochester, Minnesota, USA.
⁵ Division of Surgery Research, Mayo Clinic, Rochester, Minnesota, USA.

^# Contributed equally.

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. IMPORTANCE The genomic DNA input for stool samples utilized for microbiome composition has not been determined. In this study, we determined the reliable threshold level under which conclusions drawn from the data may be compromised. We also determined the type of microbial bias introduced by less-than-ideal genomic input.

Keywords: 16S rRNA; gDNA; microbiome.

Abstract

Grants and funding