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.