StrainR2 accurately deconvolutes strain-level abundances in synthetic microbial communities

Abstract

BackgroundSynthetic microbial communities offer an opportunity to conduct reductionist research in tractable model systems. However, deriving abundances of highly related strains within these communities is currently unreliable. 16S rRNA gene sequencing does not resolve abundance at the strain level, standard methods for analysis of shotgun metagenomic sequencing do not account for ambiguous mapping between closely related strains, and other methods such as quantitative PCR (qPCR) scale poorly and are resource prohibitive for complex communities. We present StrainR2, which utilizes shotgun metagenomic sequencing paired with a k-mer-based normalization strategy to provide high accuracy strain-level abundances for all members of a synthetic community, provided their genomes. ResultsBoth in silico, and using sequencing data derived from gnotobiotic mice colonized with a synthetic fecal microbiota, StrainR2 resolves strain abundances with greater accuracy than other tools utilizing shotgun metagenomic sequencing reads and can resolve complex mixtures of highly related strains. Through experimental validation and benchmarking, we demonstrate that StrainR2s accuracy is comparable to that of qPCR on a subset of strains resolved using absolute quantification. Further, it is capable of scaling to communities of hundreds of strains and efficiently utilizes memory being capable of running both on personal computers and high-performance computing nodes. ConclusionsUsing shotgun metagenomic sequencing reads is a viable method for determining accurate strain-level abundances in synthetic communities using StrainR2.