Genomovar-level resolution reveals rapid pathotype switching and genomovar-specific disease potential in diarrheagenic Escherichia coli populations in northern Ecuador
TL;DR
These findings demonstrate that, although there seems to be within-lineage (phylogroup) conservation of virulence, pathotype identity is a labile state defined by horizontally acquired virulence genes at the genomovar level, and that the genomovar framework provides a biologically meaningful unit for linking intra-species diversity to pathogenic potential and outbreaks.
Abstract
Diarrheagenic Escherichia coli (DEC) pathotypes are commonly defined by molecular detection of discrete virulence genes, yet how quickly these diagnostic genes emerge and move among co-circulating lineages remain unclear. Here, we classified 248 whole-genome-sequenced E. coli isolates from the EcoZUR case-control study in northern Ecuador into intra-species genomovar units using the recently described 99.5% ANI threshold. This framework exposed cryptic population structure, revealing that single sequence types, representing identical multilocus sequence types (MLST), can harbor multiple distinct genomovars. Within individual genomovars, we observed a few cases of different pathotypes among isolates showing ∼99.7% ANI (and many such cases between genomovars). Coupled with synteny and phylogeny analyses that revealed pervasive incongruences between pathotype-diagnostic virulence genes and the core genome, these findings suggest recent horizontal gene transfer as the primary driver of pathotype evolution. Virulence gene profiling further revealed that accessory virulence repertoires are hierarchically structured by phylogroup across pathotypes, with genomovars assigned to phylogroups B2 and D exhibiting more conserved virulence architectures than those in phylogroup A and B1. Among DAEC isolates specifically, the B2- and D-associated genomovars showed elevated diarrhea-association rates relative to their phylogroup A counterparts. Rare virulence genes, including Type VI secretion systems, further distinguished diarrhea-associated from asymptomatic genomovars. These findings demonstrate that, although there seems to be within-lineage (phylogroup) conservation of virulence, pathotype identity is a labile state defined by horizontally acquired virulence genes at the genomovar level, and that the genomovar framework provides a biologically meaningful unit for linking intra-species diversity to pathogenic potential and outbreaks. Importance Efforts to diagnose diarrheal Escherichia coli infections depend on our ability to reliably identify which strain is dangerous, a task that for decades has rested on sorting strains into pathotypes defined by a few virulence genes. Whether these labels mark stable lineages or fleeting states is hard to judge with traditional typing methods such as Sequence Types (STs), which group together isolates with identical sequences in a handful of housekeeping loci. Using a recently defined genome-wide threshold (99.5% ANI), we resolved isolates into fine-scale genomovars and found that individual STs often conceal multiple distinct genomovars, some carrying conflicting virulence repertoires. At this resolution, the loci that define pathotypes are gained and lost far faster than the core genome diverges, and a genomovar’s genomic background shapes its association with disease. Genomovars therefore complement MLST with the resolution needed to interpret genomic surveillance data and to build robust diagnostic and public-health frameworks.