High-resolution ISR amplicon sequencing reveals personalized oral microbiome

Chiranjit Mukherjee, Clifford J. Beall, Ann L. Griffen, Eugene J. Leys

Sequencing of the 16S rRNA gene has been the standard for studying the composition of microbial communities. While it allows identification of bacteria at the level of species, this method does not usually provide sufficient information to resolve communities at the sub-species level. Species-level resolution is not adequate for studies of transmission or stability or for exploring subspecies variation in disease association. Strain level analysis using whole metagenome shotgun sequencing has significant limitations that can make it unsuitable for large-scale studies. Achieving sufficient depth of sequencing can be cost-prohibitive, and even with adequate coverage, deconvoluting complex communities such as the oral microbiota is computationally very challenging. Thus, there is a need for high-resolution, yet cost-effective, high-throughput methods for characterizing microbial communities. Significant improvement in resolution for amplicon-based bacterial community analysis was achieved by combining amplicon sequencing of a high-diversity marker gene, the ribosomal 16-23S intergenic spacer region (ISR), with a probabilistic error modeling based denoising algorithm, DADA2. The resolving power of this new approach was compared to that of both standard and high-resolution 16S-based approaches using a set of longitudinal subgingival plaque samples. The ISR strategy resulted in a 5.2-fold increase in community resolution compared to reference-based 16S rRNA gene analysis and showed 100% accuracy in predicting the correct source of a clinical sample. Individuals' microbial communities were highly personalized, and although they exhibited some drift in membership and levels over time, that difference was always smaller than the differences between any two subjects, even after 1 year. The construction of an ISR database from publicly available genomic sequences allowed us to explore genomic variation within species, resulting in the identification of multiple variants of the ISR for most species. The ISR approach resulted in significantly improved resolution of communities and revealed a highly personalized human oral microbiota that was stable over 1 year. Multiple ISR types were observed for all species examined, demonstrating a high level of subspecies variation in the oral microbiota. The approach is high-throughput, high-resolution yet cost-effective, allowing subspecies-level community fingerprinting at a cost comparable to that of 16S rRNA gene amplicon sequencing. It will be useful for a range of applications that require high-resolution identification of organisms, including microbial tracking, community fingerprinting, and potentially for identification of virulence-associated strains.