Single-gene long-read sequencing illuminates Escherichia coli strain dynamics in the human intestinal microbiome

Gut microbiome is of major interest due to its close relationship to health and disease. Bacteria usually vary in gene content, leading to functional variations within species, so resolution higher than species-level methods is needed for ecological and clinical relevance. We design a protocol to identify strains in selected species with high discrimination and in high numbers by amplicon sequencing of the flagellin gene. We apply the protocol to fecal samples from a human diet trial, targeting Escherichia coli. Across the 119 samples from 16 individuals, there are 1,532 amplicon sequence variants (ASVs), but only 32 ASVs are dominant in one or more fecal samples, despite frequent dominant strain turnover. Major strains in an intestine are found to be commonly accompanied by a large number of satellite cells, and many are identified as potential extraintestinal pathogens. The protocol could be used to track epidemics or investigate the intra- or inter-host diversity of pathogens.