Species-level bacterial community profiling of the healthy sinonasal microbiome using Pacific Biosciences sequencing of full-length 16S rRNA genes

• Published in: bioRxiv
• Main Authors: Earl, Joshua P, Adappa, Nithin D, Krol, Jaroslaw, Bhat, Archana S, Balashov, Sergey, Ehrlich, Rachel L, Palmer, James N, Workman, Alan D, Blasetti, Mariel, Sen, Bhaswati, Hammond, Jocelyn, Cohen, Noam A, Ehrlich, Garth D, Mell, Joshua C
• Format: Paper
• Language: English
• Published: Cold Spring Harbor Cold Spring Harbor Laboratory Press 30.08.2018; Cold Spring Harbor Laboratory
• Edition: 1.3
• Subjects: Bacteria; Bioinformatics; Biopsy; Conserved sequence; DNA sequencing; Entropy; Genes; Microbiomes; Nose; Phylogeny; Pituitary; rRNA 16S; Sinus; Species; Taxonomy; 16S rRNA; circular consensus sequencing; database; sinonasal; long-read DNA sequencing; paranasal sinuses; microbiome
• ISSN: 2692-8205; 2692-8205
• DOI: 10.1101/338731

Background: Pan-bacterial 16S rRNA microbiome surveys performed with massively parallel DNA sequencing technologies have transformed community microbiological studies. Current 16S profiling methods, however, fail to provide sufficient taxonomic resolution and accuracy to adequately perform species-level associative studies for specific conditions. This is due to the amplification and sequencing of only short 16S rRNA gene regions, typically providing for only family- or genus-level taxonomy. Moreover, sequencing errors often inflate the number of taxa present. Pacific Biosciences' (PacBio's) long-read technology in particular suffers from high error rates per base. Herein we present a microbiome analysis pipeline that takes advantage of PacBio circular consensus sequencing (CCS) technology to sequence and error correct full-length bacterial 16S rRNA genes, which provides high-fidelity species-level microbiome data. Results: Analysis of a mock community with 20 bacterial species demonstrated 100% specificity and sensitivity. Examination of a 250-plus species mock community demonstrated correct species-level classification of >90% of taxa and relative abundances were accurately captured. The majority of the remaining taxa were demonstrated to be multiply, incorrectly, or incompletely classified. Using this methodology, we examined the microgeographic variation present among the microbiomes of six sinonasal sites, by both swab and biopsy, from the anterior nasal cavity to the sphenoid sinus from 12 subjects undergoing trans-sphenoidal hypophysectomy. We found greater variation among subjects than among sites within a subject, although significant within-individual differences were also observed. Propiniobacterium acnes (recently renamed Cutibacterium acnes [1]) was the predominant species throughout, but was found at distinct relative abundances by site. Conclusions: Our microbial composition analysis pipeline for single-molecule real-time 16S rRNA gene sequencing (MCSMRT, https://github.com/jpearl01/mcsmrt) overcomes deficits of standard marker gene based microbiome analyses by using CCS of entire 16S rRNA genes to provide increased taxonomic and phylogenetic resolution. Extensions of this approach to other marker genes could help refine taxonomic assignments of microbial species and improve reference databases, as well as strengthen the specificity of associations between microbial communities and dysbiotic states. Footnotes * Extensive revision to main text, and supplementary text, including: Amplicon Sequence Variant analysis via Minimum Entropy Decomposition Figures 1, 2, 4, and 8 extensively revised to include new analysis Figure 9 removed. Multiple additional supplementary figures detailing new analyses. Dr. Bhaswati Sen added to author list