Comprehensive Assessment of 16S rRNA Gene Amplicon Sequencing for Microbiome Profiling across Multiple Habitats

The 16S rRNA gene works as a rapid and effective marker for the identification of microorganisms in complex communities; hence, a huge number of microbiomes have been surveyed by 16S amplicon-based sequencing. The resolution of the 16S rRNA gene is always considered only at the genus level; however, it has not been verified on a wide range of microbes yet. To fully explore the ability and potential of the 16S rRNA gene in microbial profiling, here, we propose Qscore, a comprehensive method to evaluate the performance of amplicons by integrating the amplification rate, multitier taxonomic annotation, sequence type, and length. Our in silico assessment by a "global view" of 35,889 microbe species across multiple reference databases summarizes the optimal sequencing strategy for 16S short reads. On the other hand, since microbes are unevenly distributed according to their habitats, we also provide the recommended configuration for 16 typical ecosystems based on the Qscores of 157,390 microbiomes in the Microbiome Search Engine (MSE). Detailed data simulation further proves that the 16S amplicons produced with Qscore-suggested parameters exhibit high precision in microbiome profiling, which is close to that of shotgun metagenomes under CAMI metrics. Therefore, by reconsidering the precision of 16S-based microbiome profiling, our work not only enables the high-quality reusability of massive sequence legacy that has already been produced but is also significant for guiding microbiome studies in the future. We have implemented the Qscore as an online service at to parse the recommended sequencing strategy for specific habitats or expected microbial structures.