Mapping the evolutionary path towards multi-drug resistance in the pandemic Escherichia coli ST131 lineage

Abstract Escherichia coli sequence type (ST) 131 is the most widely studied genetic lineage of the species to date, originally identified in the early 2000s as an increasingly common cause of human urinary tract and bloodstream infections worldwide. This lineage is subdivided into four extant main subclades A, B, C1 and C2 that exhibit distinct features in terms of invasiveness, antibiotic resistance and transmissibility. However, the evolutionary pathway from the generally susceptible ST131-B to the drug-resistant ST131-C clades remains poorly mapped. To fill this knowledge gap, we analysed in detail human clinical isolates obtained in Vietnam, designated as belonging to the generally neglected minor clade ST131-B0. We sequenced them using both short- and long-read technology, and combined these data with a recently published high-resolution genomic collection to provide further insight into the evolutionary process and its timeline. Extensive genomic analyses established ST131-B0 as an intermediary progenitor in the evolutionary path leading from ST131-B to the ST131-C clades, associated with multiple type I pili switches, as well as the loss and gain of specific chromosomal genes representing diverse core functions such as metabolism, transcription, DNA binding and type II toxin-antitoxin systems. Furthermore, all Vietnamese isolates of ST131-B0 unprecedentedly harboured bla CTX-M genes encoding extended-spectrum β-lactamases, a trait dominant in ST131-C clades and not previously observed in ST131-B0. Our study supports the hypothesis that the ST131-C clades have gradually evolved from ST131-B by reducing the host range with better adaptation to colonising humans under selective conditions.