A simple method that enhances minority species detection in the microbiota: 16S metagenome-DRIP (Deeper Resolution using an Inhibitory Primer)

Abstract

Aim: 16S rRNA gene-based microbiota analyses (16S metagenomes) using next-generation sequencing (NGS) technologies are widely used to examine the microbial community composition in environmental samples. However, the sequencing capacity of NGS is sometimes insufficient to cover the whole microbial community, especially when analyzing soil and fecal microbiotas. This limitation may have hampered the detection of minority species that potentially affect microbiota formation and structure.

Methods: We developed a simple method, termed 16S metagenome-DRIP (Deeper Resolution using an Inhibitory Primer), that not only enhances minority species detection but also increases the accuracy of their abundance estimation. The method relies on the inhibition of normal amplicon formation of the 16S rRNA gene of a target major (abundant) species during the first PCR step. The addition of a biotinylated primer that is complementary to the variable sequence of the V3-V4 region of the target species inhibits a normal amplification process to form an aberrant short amplicon. The fragment is then captured by streptavidin beads for removal from the reaction mixture, and the resulting mixture is utilized for the second PCR with barcode-tag primers. Thus, this method only requires two additional experimental procedures to the conventional 16S metagenome analysis. A proof-of-concept experiment was first conducted using a mock sample consisting of the genomes of 14 bacterial species. Then, the method was applied to infant fecal samples using a Bifidobacterium-specific inhibitory primer (n = 11).

Results: As a result, the reads assigned to the family Bifidobacteriaceae decreased on average from 16,657 to 1718 per sample without affecting the total read counts (36,073 and 34,778 per sample for the conventional and DRIP methods, respectively). Furthermore, the minority species detection rate increased with neither affecting Bray-Curtis dissimilarity calculated by omitting the target Bifidobacterium species (median: 0.049) nor changing the relative abundances of the non-target species. While 115 amplicon sequence variants (ASVs) were unique to the conventional method, 208 ASVs were uniquely detected for the DRIP method. Moreover, the abundance estimation for minority species became more accurate, as revealed thorough comparison with the results of quantitative PCR analysis.

Conclusion: The 16S metagenome-DRIP method serves as a useful technique to grasp a deeper and more accurate microbiota composition when combined with conventional 16S metagenome analysis methods.

Keywords

16S metagenome-DRIP / α-Type DNA polymerase / biotinylated inhibitory primer / minority species / streptavidin-beads

References

[1] Paine RT.A note on trophic complexity and community stability.The American Naturalist1969;103:91-3

[2] Power ME,Estes JA.Challenges in the quest for keystones.BioScience1996;46:609-20

[3] Ze X,Louis P.Ruminococcus bromii is a keystone species for the degradation of resistant starch in the human colon.ISME J2012;6:1535-43 PMCID:PMC3400402

[4] Centanni M,Butts CA.Bifidobacterium pseudolongum in the ceca of rats fed hi-maize starch has characteristics of a keystone species in bifidobacterial blooms.Appl Environ Microbiol2018;84:e00547-18 PMCID:PMC6052269

[5] Yamada C,Sakanaka M.Molecular insight into evolution of symbiosis between breast-fed infants and a member of the human gut microbiome bifidobacterium longum.Cell Chem Biol2017;24:515-524.e5

[6] Raynaud X.Spatial ecology of bacteria at the microscale in soil.PLoS One2014;9:e87217 PMCID:PMC3905020

[7] Sender R,Milo R.Are we really vastly outnumbered?.Cell2016;164:337-40

[8] Odamaki T,Sugahara H.Age-related changes in gut microbiota composition from newborn to centenarian: a cross-sectional study.BMC Microbiol2016;16:90 PMCID:PMC4879732

[9] Gotoh A,Sakanaka M.Sharing of human milk oligosaccharides degradants within bifidobacterial communities in faecal cultures supplemented with Bifidobacterium bifidum.Sci Rep2018;8:13958 PMCID:PMC6143587

[10] Untergasser A,Koressaar T.Primer3 - new capabilities and interfaces.Nucleic Acids Res2012;40:e115 PMCID:PMC3424584

[11] Ye J,Zaretskaya I,Rozen S.Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction.BMC Bioinformatics2012;13:134 PMCID:PMC3412702

[12] Langmead B.Fast gapped-read alignment with Bowtie 2.Nat Methods2012;9:357-9 PMCID:PMC3322381

[13] Callahan BJ,Rosen MJ,Johnson AJ.DADA2: high-resolution sample inference from illumina amplicon data.Nat Methods2016;13:581-3 PMCID:PMC4927377

[14] Clifford RJ,Prestwood J.Detection of bacterial 16S rRNA and identification of four clinically important bacteria by real-time PCR.PLoS One2012;7:e48558 PMCID:PMC3490953

[15] R Core Team. R: a language and environment for statistical computing. Vienna, Austria; 2021. Available from: https://www.r-project.org/ [Last accessed on 25 May 2022]

[16] Harrell Jr FE. Hmisc: Harrell Miscellaneous; 2021. Available from: https://cran.r-project.org/package=Hmisc [Last accessed on 25 May 2022]

[17] Oksanen J,Friendly M. vegan: Community Ecology Package; 2020. Available from: https://cran.r-project.org/package=vegan [Last accessed on 25 May 2022]

[18] Nishijima S,Oshima K.The gut microbiome of healthy Japanese and its microbial and functional uniqueness.DNA Res2016;23:125-33 PMCID:PMC4833420

[19] Sakanaka M,Gotoh A.Evolutionary adaptation in fucosyllactose uptake systems supports bifidobacteria-infant symbiosis.Sci Adv2019;5:eaaw7696 PMCID:PMC6713505

[20] Gonzalez JM,Belda-Ferre P.Amplification by PCR artificially reduces the proportion of the rare biosphere in microbial communities.PLoS One2012;7:e29973 PMCID:PMC3256211