Strain Identification and Quantitative Analysis in Microbial Communities

Microbiology has long studied the ways in which subtle genetic differences between closely related microbial strains can have profound impacts on their phenotypes and those of their surrounding environments and communities. Despite the growth in high-throughput microbial community profiling, however, such strain-level differences remain challenging to detect. Once detected, few quantitative approaches have been well-validated for associating strain variants from microbial communities with phenotypes of interest, such as medication usage, treatment efficacy, host environment, or health. First, the term “strain” itself is not used consistently when defining a highly-resolved taxonomic or genomic unit from within a microbial community. Second, computational methods for identifying such strains directly from shotgun metagenomics are difficult, with several possible reference- and assembly-based approaches available, each with different sensitivity/specificity tradeoffs. Finally, statistical challenges exist in using any of the resulting strain profiles for downstream analyses, which can include strain tracking, phylogenetic analysis, or genetic association studies. We provide an in depth discussion of recently available computational tools that can be applied for this task, as well as statistical models and gaps in performing and interpreting any of these three main types of studies using strain-resolved shotgun metagenomic profiling of microbial communities.