Genomovar-level resolution reveals rapid pathotype switching and genomovar-specific disease potential in diarrheagenic Escherichia coli populations in northern Ecuador
Dorian J. FeisteliD, Kelsey J. JesseriD, Karen LevyiD, Gabriel A. TruebaiD, Konstantinos T. KonstantinidisiD
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.
02 /Specification & Classification
Click any tag to browse matching papers- Citations
- 0
- Organism
- Assay
- Research Focus
- Approach
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Functional Genomics of Gastrointestinal Escherichia coli Isolated from Patients with Cancer and Diarrhea
H. Carter, J. Clark, … P. C. Okhuysen·bioRxiv·May 31, 2023
- Strain intrinsic properties and environmental constraints together shape Escherichia coli dynamics and diversity over a twenty-year human gut time series
Bénédicte Condamine, Thibaut Morel‐Journel, … Olivier Clermont·The ISME JournalPeer-reviewed·Dec 10, 2024
- Enteroaggregative Escherichia clade I from Nigeria
Rotimi Ayodeji Dada, Olabisi Comfort Akinlabi, … Iruka N. Okeke·bioRxiv·Apr 22, 2026
- Genomic epidemiology of enteropathogenic Escherichia coli (EPEC) in southwestern Nigeria
Olabisi Comfort Akinlabi, Rotimi Ayodeji Dada, … Iruka N. Okeke·bioRxiv·Feb 14, 2025
- Captive green iguana is a reservoir of diarrheogenic Escherichia coli pathotypes
G. U. Bautista-Trujillo, F. A. Gutierrez-Miceli, … J. Gutierrez-Jimenez·bioRxiv·May 2, 2019
- The population genomics of increased virulence and antibiotic resistance in human commensal Escherichia coli over 30 years in France
Julie Marin, Olivier Clermont, … François Blanquart·Applied and Environmental MicrobiologyPeer-reviewed·Jul 18, 2022