A comparison of approaches to scaffolding multiple regions along the 16S rRNA gene for improved resolution

Abstract

MotivationFull length, high resolution 16s rRNA marker gene sequencing has been challenging historically. Short amplicons provide high accuracy reads with widely available equipment, at the cost of taxonomic resolution. One recent proposal has been to reconstruct multiple amplicons along the full-length marker gene, however no barcode-free computationally tractable approach for this is available. To address this gap, we present Sidle (SMURF Implementation Done to acceLerate Efficiency), an implementation of the Short MUltiple Reads Framework algorithm with a novel tree building approach to reconstruct rRNA genes from individually amplified regions. ResultsUsing simulated and real data, we compared Sidle to two other approaches of leveraging multiple gene region data. We found that Sidle had the least bias in non-phylogenetic alpha diversity, feature-based measures of beta diversity, and the reconstruction of individual clades. With a curated database, Sidle also provided the most precise species-level resolution. Availability and ImplementationSidle is available under a BSD 3 license from https://github.com/jwdebelius/q2-sidle