A novel barcoded nanopore sequencing workflow of high-quality, full-length bacterial 16S amplicons for taxonomic annotation of bacterial isolates and complex microbial communities

Due to recent improvements, Nanopore sequencing has become a promising method for experiments relying on amplicon sequencing. We describe a flexible workflow to generate and annotate high-quality, full-length 16S rDNA amplicons. We evaluated it for two applications, namely, (i) identification of bacterial isolates and (ii) species-level profiling of microbial communities. We assessed the identification of single bacterial isolates by sequencing, using a set of barcoded full-length 16S rRNA gene primer pairs (pair A), on 47 isolates encompassing multiple genera and compared those results with matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS)-based identification. Species-level community profiling was tested with two sets of barcoded full-length 16S primer pairs (A and B) and compared to the results obtained with shotgun Illumina sequencing using 27 stool samples. We developed a Nextflow pipeline to retain high-quality reads and taxonomically annotate them. We found high agreement between our workflow and MALDI-TOF data for isolate identification (positive predictive value = 0.90, Cramer's V = 0.857, and Theil's U = 0.316). For species-level community profiling, we found strong correlations (r(s) > 0.6) of alpha diversity indices between the two primer sets and Illumina sequencing. At the community level, we found significant but small differences when comparing sequencing techniques. Finally, we found a moderate to strong correlation when comparing the relative abundances of individual species (average r(s) = 0.6 and 0.533 for primers A and B). Despite identified shortcomings, the proposed workflow enabled accurate identification of single bacterial isolates and prominent features in microbial communities, making it a worthwhile alternative to MALDI-TOF MS and Illumina sequencing.