The impact of primer design on amplicon-based metagenomic profiling accuracy: Detailed insights into bifidobacterial community structure

Abstract

Next Generation Sequencing (NGS) technologies have overcome the limitations of cultivation-dependent approaches and allowed detailed study of bacterial populations that inhabit the human body. The consortium of bacteria residing in the human intestinal tract, also known as the gut microbiota, impacts several physiological processes important for preservation of the health status of the host. The most widespread microbiota profiling method is based on amplification and sequencing of a variable portion of the 16S rRNA gene as a universal taxonomic marker among members of the Bacteria domain. Despite its popularity and obvious advantages, this 16S rRNA gene-based approach comes with some important limitations. In particular, the choice of the primer pair for amplification plays a major role in defining the accuracy of the reconstructed bacterial profiles. In the current study, we performed an in silico PCR using all currently described 16S rRNA gene-targeting primer pairs (PP) in order to assess their efficiency. Our results show that V3, V4, V5, and V6 were the optimal regions on which to design 16S rRNA metagenomic primers. In detail, PP39 (Probio_Uni/Probio_Rev), PP41 (341F/534R), and PP72 (970F/1050R) were the most suitable primer pairs with an amplification efficiency of >98.5%. Furthermore, the Bifidobacterium genus was examined as a test case for accurate evaluation of intra-genus performances at subspecies level. Intriguingly, the in silico analysis revealed that primer pair PP55 (527f/1406r) was unable to amplify the targeted region of any member of this bacterial genus, while several other primer pairs seem to rather inefficiently amplify the target region of the main bifidobacterial taxa. These results highlight that selection of a 16S rRNA gene-based PP should be done with utmost care in order to avoid biases in microbiota profiling results.