Phylogroup stability contrasts with high within sequence type complex dynamics of Escherichia coli bloodstream infection isolates over a 12-year period

Impressive pan-genomic diversity of E. coli from a wild animal community near urban development reflects human impacts

Human and domesticated animal waste infiltrates global freshwater, terrestrial, and marine environments, widely disseminating fecal microbes, antibiotics, and other chemical pollutants. Emerging evidence suggests that guts of wild animals are being invaded by our microbes, including Escherichia coli, which face anthropogenic selective pressures to gain antimicrobial resistance (AMR) and increase virulence. However, wild animal sources remain starkly under-represented among genomic sequence repositories. We sequenced whole genomes of 145 E. coli isolates from 55 wild and 13 domestic animal fecal samples, averaging 2 (ranging 1-7) isolates per sample, on a preserve imbedded in a human-dominated landscape in California Bay Area, USA, to assess AMR, virulence, and pan-genomic diversity. With single nucleotide polymorphism analyses we predict potential transmission routes. We illustrate the usefulness of E. coli to aid our understanding of and ability to surveil the emergence of zoonotic pathogens created by the mixing of human and wild bacteria in the environment.

Genomic insights into virulence, antimicrobial resistance, and adaptation acumen of Escherichia coli isolated from an urban environment

Populations of common commensal bacteria such as Escherichia coli undergo genetic changes by the acquisition of certain virulence and antimicrobial resistance (AMR) encoding genetic elements leading to the emergence of pathogenic strains capable of surviving in the previously uninhabited or protected niches. These bacteria are also reported to be prevalent in the environment where they survive by adopting various recombination strategies to counter microflora of the soil and water, under constant selection pressure(s). In this study, we performed molecular characterization, phenotypic AMR analysis, and whole genome sequencing (WGS) of E. coli (n = 37) isolated from soil and surface water representing the urban and peri-urban areas. The primary aim of this study was to understand the genetic architecture and pathogenic acumen exhibited by environmental E. coli. WGS-based analysis entailing resistome and virulome profiling indicated the presence of various virulence (adherence, iron uptake, and toxins) and AMR encoding genes, including blaNDM-5 in the environmental isolates. A majority of our isolates belonged to phylogroup B1 (73%). A few isolates in our collection were of sequence type(s) (ST) 58 and 224 that could have emerged recently as clonal lineages and might pose risk of infection/transmission. Mobile genetic elements (MGEs) such as plasmids (predominantly) of the IncF family, prophages, pipolins, and insertion elements such as IS1 and IS5 were also observed to exist, which may presumably aid in the propagation of genes encoding resistance against antimicrobial drugs. The observed high prevalence of MGEs associated with multidrug resistance in pathogenic E. coli isolates belonging to the phylogroup B1 underscores the need for extended surveillance to keep track of and prevent the transmission of the bacterium to certain vulnerable human and animal populations.

IMPORTANCE

Evolutionary patterns of E. coli bacteria convey that they evolve into highly pathogenic forms by acquiring fitness advantages, such as AMR, and various virulence factors through the horizontal gene transfer (HGT)-mediated acquisition of MGEs. However, limited research on the genetic profiles of environmental E. coli, particularly from India, hinders our understanding of their transition to pathogenic forms and impedes the adoption of a comprehensive approach to address the connection between environmentally dwelling E. coli populations and human and veterinary public health. This study focuses on high-resolution genomic analysis of the environmental E. coli isolates aiming to understand the genetic similarities and differences among isolates from different environmental niches and uncover the survival strategies employed by these bacteria to thrive in their surroundings. Our approach involved molecular characterization of environmental samples using PCR-based DNA fingerprinting and subsequent WGS analysis. This multidisciplinary approach is likely to provide valuable insights into the understanding of any potential spill-over to human and animal populations and locales. Investigating these environmental isolates has significant potential for developing epidemiological strategies against transmission and understanding niche-specific evolutionary patterns.

Temporal dynamics of genetically heterogeneous extended-spectrum cephalosporin-resistant Escherichia coli bloodstream infections

Extended-spectrum cephalosporin-resistant Escherichia coli (ESC-R-Ec) is an urgent public health threat with sequence type clonal complex 131 (STc131), phylogroup B2 strains being particularly concerning as the dominant cause of ESC-R-Ec infections. To address the paucity of recent ESC-R-Ec molecular epidemiology data in the United States, we used whole-genome sequencing (WGS) to fully characterize a large cohort of invasive ESC-R-Ec at a tertiary care cancer center in Houston, Texas, collected from 2016 to 2020. During the study time frame, there were 1,154 index E. coli bloodstream infections (BSIs) of which 389 (33.7%) were ESC-R-Ec. Using time series analyses, we identified a temporal dynamic of ESC-R-Ec distinct from ESC-susceptible E. coli (ESC-S-Ec), with cases peaking in the last 6 months of the calendar year. WGS of 297 ESC-R-Ec strains revealed that while STc131 strains accounted for ~45% of total BSIs, the proportion of STc131 strains remained stable across the study time frame with infection peaks driven by genetically heterogeneous ESC-R-Ec clonal complexes. bla CTX-M variants accounted for most β-lactamases conferring the ESC-R phenotype (89%; 220/248 index ESC-R-Ec), and amplification of bla CTX-M genes was widely detected in ESC-R-Ec strains, particularly in carbapenem non-susceptible, recurrent BSI strains. Bla CTX-M-55 was significantly enriched within phylogroup A strains, and we identified bla CTX-M-55 plasmid-to-chromosome transmission occurring across non-B2 strains. Our data provide important information regarding the current molecular epidemiology of invasive ESC-R-Ec infections at a large tertiary care cancer center and provide novel insights into the genetic basis of observed temporal variability for these clinically important pathogens. IMPORTANCE Given that E. coli is the leading cause of worldwide ESC-R Enterobacterales infections, we sought to assess the current molecular epidemiology of ESC-R-Ec using a WGS analysis of many BSIs over a 5-year period. We identified fluctuating temporal dynamics of ESC-R-Ec infections, which have also recently been identified in other geographical regions such as Israel. Our WGS data allowed us to visualize the stable nature of STc131 over the study period and demonstrate a limited but genetically diverse group of ESC-R-Ec clonal complexes are detected during infection peaks. Additionally, we provide a widespread assessment of β-lactamase gene copy number in ESC-R-Ec infections and delineate mechanisms by which such amplifications are achieved in a diverse array of ESC-R-Ec strains. These data suggest that serious ESC-R-Ec infections are driven by a diverse array of strains in our cohort and impacted by environmental factors suggesting that community-based monitoring could inform novel preventative measures.

The bacterial genetic determinants of Escherichia coli capacity to cause bloodstream infections in humans

Escherichia coli is both a highly prevalent commensal and a major opportunistic pathogen causing bloodstream infections (BSI). A systematic analysis characterizing the genomic determinants of extra-intestinal pathogenic vs. commensal isolates in human populations, which could inform mechanisms of pathogenesis, diagnostic, prevention and treatment is still lacking. We used a collection of 912 BSI and 370 commensal E. coli isolates collected in France over a 17-year period (2000-2017). We compared their pangenomes, genetic backgrounds (phylogroups, STs, O groups), presence of virulence-associated genes (VAGs) and antimicrobial resistance genes, finding significant differences in all comparisons between commensal and BSI isolates. A machine learning linear model trained on all the genetic variants derived from the pangenome and controlling for population structure reveals similar differences in VAGs, discovers new variants associated with pathogenicity (capacity to cause BSI), and accurately classifies BSI vs. commensal strains. Pathogenicity is a highly heritable trait, with up to 69% of the variance explained by bacterial genetic variants. Lastly, complementing our commensal collection with an older collection from 1980, we predict that pathogenicity continuously increased through 1980, 2000, to 2010. Together our findings imply that E. coli exhibit substantial genetic variation contributing to the transition between commensalism and pathogenicity and that this species evolved towards higher pathogenicity.

Epistatic interactions between the high pathogenicity island and other iron uptake systems shape Escherichia coli extra-intestinal virulence

The intrinsic virulence of extra-intestinal pathogenic Escherichia coli is associated with numerous chromosomal and/or plasmid-borne genes, encoding diverse functions such as adhesins, toxins, and iron capture systems. However, the respective contribution to virulence of those genes seems to depend on the genetic background and is poorly understood. Here, we analyze genomes of 232 strains of sequence type complex STc58 and show that virulence (quantified in a mouse model of sepsis) emerged in a sub-group of STc58 due to the presence of the siderophore-encoding high-pathogenicity island (HPI). When extending our genome-wide association study to 370 Escherichia strains, we show that full virulence is associated with the presence of the aer or sit operons, in addition to the HPI. The prevalence of these operons, their co-occurrence and their genomic location depend on strain phylogeny. Thus, selection of lineage-dependent specific associations of virulence-associated genes argues for strong epistatic interactions shaping the emergence of virulence in E. coli.

Evolutionary and functional history of the Escherichia coli K1 capsule

Escherichia coli is a leading cause of invasive bacterial infections in humans. Capsule polysaccharide has an important role in bacterial pathogenesis, and the K1 capsule has been firmly established as one of the most potent capsule types in E. coli through its association with severe infections. However, little is known about its distribution, evolution and functions across the E. coli phylogeny, which is fundamental to elucidating its role in the expansion of successful lineages. Using systematic surveys of invasive E. coli isolates, we show that the K1-cps locus is present in a quarter of bloodstream infection isolates and has emerged in at least four different extraintestinal pathogenic E. coli (ExPEC) phylogroups independently in the last 500 years. Phenotypic assessment demonstrates that K1 capsule synthesis enhances E. coli survival in human serum independent of genetic background, and that therapeutic targeting of the K1 capsule re-sensitizes E. coli from distinct genetic backgrounds to human serum. Our study highlights that assessing the evolutionary and functional properties of bacterial virulence factors at population levels is important to better monitor and predict the emergence of virulent clones, and to also inform therapies and preventive medicine to effectively control bacterial infections whilst significantly lowering antibiotic usage.

Bacterial DNA metabolism analysis by metagenomic next-generation sequencing (mNGS) after treatment of bloodstream infection

BACKGROUND

With the advent of metagenomic next-generation sequencing (mNGS), the time of DNA metabolism can be explored after bacteria be killed. In this study, we applied mNGS in investigation of the clearance profile of circulating bacteria DNA.

METHODS

All of the rabbits were injected with the inactivated Escherichia coli. Using mNGS, we analyzed serial samples of plasma collected from rabbits to detect clearance profile of circulating E. coli DNA.

RESULTS

In this study, we found that the of E. coli DNA could still be detected 6 h after injecting killed bacteria. The clearance half-lives associated with the 2 phases are 0.37 and 1.81 h. We also explored there is no correlation between the disease severity with the E. coli DNA reads in circulation.

CONCLUSIONS

After the bacteria were completely killed, their DNA could still be detected in the blood circulation. The metabolism of bacterial DNA in the circulation had two phases: fast and slow phases, and there were no correlations between the level of bacteria reads with the severity of patients' disease after the bacteria have been completely killed.

Effect of a karst system (France) on extended spectrum beta-lactamase (ESBL)-producing Escherichia coli

Karst aquifers are an important water resource worldwide particularly exposed to anthropogenic pollution, including antibiotic-resistance. The release of antibiotic-resistant bacterial pathogens in the environment is a major public health challenge worldwide. In this One Health study, we aimed to determine the effect of karst on antibiotic-resistant bacteria. For this purpose, we determined the concentrations of extended-spectrum β-lactamases-producing Escherichia coli (ESBL-Ec) for 92 weeks in a rural karst hydrosystem providing drinking water. ESBL-Ec isolates (n = 130) were sequenced by whole genome sequencing. We analysed the isolates at different levels of granularity, i.e., phylogroup, sequence type, presence of antibiotic-resistance genes, mutations conferring antibiotic-resistance, and virulence genes. The ESBL-Ec concentrations were spatially and temporally heterogeneous in the studied karst hydrosystem. ESBL-Ec isolates survived in the karst and their concentrations were mostly explained by the hydrodynamic of the hydrosystem. We demonstrate that the studied karst has no filtration effect on ESBL-Ec, either quantitatively (i.e., in the ESBL-Ec concentrations) or qualitatively (i.e., in the genetic characteristics of ESBL-Ec isolates).

Dynamic Periodicity of Extended-Spectrum Cephalosporin Resistant Bloodstream Escherichia coli Infections Driven by Non-CC131 Strains

Extended-spectrum cephalosporin resistant (ESC-R) Escherichia coli (ESC-R- Ec ) is an urgent public health threat with clonal complex (CC) 131, phylogroup B2 strains being particularly concerning as the dominant cause of ESC-R- Ec infections. To address the paucity of recent ESC-R- Ec molecular epidemiology data in the United States (US), we used whole genome sequencing (WGS) to fully characterize a large cohort of invasive ESC-R- Ec at a tertiary care cancer center in Houston, Texas collected from 2016-2020. During the study timeframe, there were 1154 index E. coli bloodstream infections (BSIs) of which 389 (33.7%) were ESC-R. Using time series analyses, we identified a temporal dynamic of ESC-R E. coli BSIs ( Ec -BSIs), distinct from ESC-susceptible Ec -BSIs, with cases peaking in the last 6 months of the calendar year. WGS of 297 ESC-R Ec -BSI strains revealed that while CC131 strains accounted for ∼45% of total infections, the proportion of CC131 strains remained stable across the time-period, and infection peaks were driven by genetically diverse, non-CC131 isolates. Bla _CTX-M variants accounted for most β-lactamases conferring the ESC-R phenotype (89%; 220/248 index ESC-R Ec -BSIs), and amplification of bla _CTX-M genes was widely detected in ESC-R Ec -BSI strains, particularly in carbapenem non-susceptible strains and in strains causing recurrent BSIs. Bla _{CTX- M-55} was significantly enriched within phylogroup A strains, and we identified bla _CTX-M-55 plasmid-to-chromosome transmission occurring across non-B2 strains. Our data provide important information regarding the current molecular epidemiology of invasive ESC-R E. coli and provide novel insights into the genetic basis of observed temporal variability for these clinically important pathogens.

IMPORTANCE

Given that E. coli is the leading cause of ESC-R Enterobacterales infection, we sought to assess the current molecular epidemiology of ESC-R- Ec using a WGS based analysis of many BSIs over a several year period. We identified a clear temporal dynamic of ESC-R- Ec infections, which has also recently been identified in other geographical regions such as Israel. Our WGS data allowed us to visualize the stable nature of CC131 over the study period and demonstrate that non-CC131 strains drove the infection peaks. Additionally, we provide the first widespread assessment of β-lactamase gene copy number in ESC-R- Ec infections and delineate mechanisms by which such amplifications are achieved in a diverse array of ESC-R- Ec strains. These data suggest that even for a tertiary care center, serious ESC-R- Ec infections are driven by a diverse array of strains and impacted by environmental factors suggesting that community-based monitoring could inform novel preventative measures.

Epidemiology, Clinical Features, and Antimicrobial Resistance of Invasive Escherichia Coli Disease in Patients Admitted in Tertiary Care Hospitals

Background

Invasive Escherichia coli disease (IED), including bloodstream infection, sepsis, and septic shock, can lead to high hospitalization and mortality rates. This multinational study describes the clinical profile of patients with IED in tertiary care hospitals.

Methods

We applied clinical criteria of systemic inflammatory response syndrome (SIRS), sepsis, or septic shock to patients hospitalized with culture-confirmed E coli from urine or a presumed sterile site. We assessed a proposed clinical case definition against physician diagnoses.

Results

Most patients with IED (N = 902) were adults aged ≥60 years (76.5%); 51.9%, 25.1%, and 23.0% of cases were community-acquired (CA), hospital-acquired (HA), and healthcare-associated (HCA), respectively. The urinary tract was the most common source of infection (52.3%). Systemic inflammatory response syndrome, sepsis, and septic shock were identified in 77.4%, 65.3%, and 14.1% of patients, respectively. Patients >60 years were more likely to exhibit organ dysfunction than those ≤60 years; this trend was not observed for SIRS. The case-fatality rate (CFR) was 20.0% (60-75 years, 21.5%; ≥75 years, 22.2%), with an increase across IED acquisition settings (HA, 28.3%; HCA, 21.7%; CA, 15.2%). Noticeably, 77.8% of patients initiated antibiotic use on the day of culture sample collection. A total of 65.6% and 40.8% of E coli isolates were resistant to ≥1 agent in ≥1 or ≥2 drug class(es). A 96.1% agreement was seen between the proposed clinical case definition and physician's diagnoses of IED.

Conclusions

This study contributes valuable, real-world data about IED severity. An accepted case definition could promote timely and accurate diagnosis of IED and inform the development of novel preventative strategies.

Genomic Characterization of ESBL/AmpC-Producing Escherichia coli in Stray Dogs Sheltered in Yangzhou, China

Purpose

Limited data are available on the prevalence and antimicrobial resistance of extended spectrum β-lactamase- (ESBL) and AmpC β-lactamase-producing Escherichia coli in stray dogs. We aimed to investigate the genomic characteristics of ESBL/AmpC-producing E. coli isolated from stray dogs sheltered in Yangzhou, China.

Methods

We collected 156 samples including 115 fecal swabs, 35 kennel floor swabs, two breeder hand and shoe sole swabs, and four feed samples. The isolates were tested for resistance by antimicrobial susceptibility testing and further analyzed for cefotaxime-resistant E. coli isolates by whole genome sequencing.

Results

We identified 80 cefotaxime-resistant E. coli isolates (51.3%), 59 isolates (73.8%) from feces and 21 (26.2%) from the environment. Whole-genome sequencing analysis showed that bla _CTX-M-15 (n=30) and bla _CTX-M-55 (n=29) were the most prevalent genotypes. Two isolates only carried the AmpC β-lactamase gene bla _CMY-2; one isolate had a combination of AmpC β-lactamase gene bla _DHA-1 and ESBL β-lactamase gene bla _CTX-M-14. Other important resistance genes such as bla _OXA-10, bla _TEM-1B, bla _TEM-135, bla _TEM-106, tet(A), qnrS1, qnrB4, and oqxAB were also detected. The serotype combination was highly abundant, with O10:H25 predominating (n=12). Most cefotaxime-resistant E. coli isolates belonged to phylogroup A (62.5%, n=50), followed by phylogroup B1 (26.3%, n=21). Thirty different sequence types (STs) and 27 distinct plasmid replicons were identified, among which ST2325 (n=12) and IncFII (n=38) was the most frequent ST and plasmid, respectively. ESBL/AmpC-producing isolates were divided into four major clades; clade IV was the primary lineage containing 37 isolates from feces and 13 from the environment. Three high-risk E. coli clone ST23 strains and one ST10 strain belonged to clades III and IV, respectively.

Conclusion

Our study provides a comprehensive overview of resistance profiles and genomic characteristics in ESBL/AmpC-producing E. coli and highlights the possible role of stray dogs as an antibiotic resistance gene reservoir.

An ST131 clade and a phylogroup A clade bearing an O101-like O-antigen cluster predominate among bloodstream Escherichia coli isolates from South-West Nigeria hospitals

Escherichia coli bloodstream infections are typically attributed to a limited number of lineages that carry virulence factors associated with invasiveness. In Nigeria, the identity of circulating clones is largely unknown and surveillance of their antimicrobial resistance has been limited. We verified and whole-genome sequenced 68 2016-2018 bloodstream E. coli isolates from three sentinel sites in South-Western Nigeria and susceptibility tested 67 of them. Resistance to antimicrobials commonly used in Nigeria was high, with 67 (100 %), 62 (92.5 %), 53 (79.1 %) and 37 (55.2 %) showing resistance to trimethoprim, ampicillin, ciprofloxacin and aminoglycosides, respectively. Thirty-five (51 %) isolates carried extended-spectrum β-lactamase genes and 32 (91 %) of these were multidrug resistant. All the isolates were susceptible to carbapenems and colistin. The strain set included globally disseminated high-risk clones from sequence type (ST)12 (2), ST131 (12) and ST648 (4). Twenty-three (33.8 %) of the isolates clustered within two clades. The first of these consisted of ST131 strains, comprising O16:H5 and O25:H4 sub-lineages. The second was an ST10-ST167 complex clade comprising strains carrying O-antigen and capsular genes of likely Klebsiella origin, identical to those of avian pathogenic E. coli Sanji, and serotyped in silico as O89, O101 or ONovel32, depending on the tool used. Four temporally associated ST90 strains from one sentinel were closely related enough to suggest that at least some of them represented a retrospectively detected outbreak cluster. Our data implicate a broad repertoire of E. coli isolates associated with bloodstream infections in South-West Nigeria. Continued genomic surveillance is valuable for tracking clones of importance and for outbreak identification.

The cnf1 gene is associated with an expanding Escherichia coli ST131 H30Rx/C2 subclade and confers a competitive advantage for gut colonization

Epidemiological projections point to acquisition of ever-expanding multidrug resistance (MDR) by Escherichia coli, a commensal of the digestive tract and a source of urinary tract pathogens. Bioinformatics analyses of a large collection of E. coli genomes from EnteroBase, enriched in clinical isolates of worldwide origins, suggest the Cytotoxic Necrotizing Factor 1 (CNF1)-toxin encoding gene, cnf1, is preferentially distributed in four common sequence types (ST) encompassing the pandemic E. coli MDR lineage ST131. This lineage is responsible for a majority of extraintestinal infections that escape first-line antibiotic treatment, with known enhanced capacities to colonize the gastrointestinal tract. Statistical projections based on this dataset point to a global expansion of cnf1-positive multidrug-resistant ST131 strains from subclade H30Rx/C2, accounting for a rising prevalence of cnf1-positive strains in ST131. Despite the absence of phylogeographical signals, cnf1-positive isolates segregated into clusters in the ST131-H30Rx/C2 phylogeny, sharing a similar profile of virulence factors and the same cnf1 allele. The suggested dominant expansion of cnf1-positive strains in ST131-H30Rx/C2 led us to uncover the competitive advantage conferred by cnf1 for gut colonization to the clinical strain EC131GY ST131-H30Rx/C2 versus cnf1-deleted isogenic strain. Complementation experiments showed that colon tissue invasion was compromised in the absence of deamidase activity on Rho GTPases by CNF1. Hence, gut colonization factor function of cnf1 was confirmed for another clinical strain ST131-H30Rx/C2. In addition, functional analysis of the cnf1-positive clinical strain EC131GY ST131-H30Rx/C2 and a cnf1-deleted isogenic strain showed no detectable impact of the CNF1 gene on bacterial fitness and inflammation during the acute phase of bladder monoinfection. Together these data argue for an absence of role of CNF1 in virulence during UTI, while enhancing gut colonization capacities of ST131-H30Rx/C2 and suggested expansion of cnf1-positive MDR isolates in subclade ST131-H30Rx/C2.

Characterizing Antimicrobial Resistance in Clinically Relevant Bacteria Isolated at the Human/Animal/Environment Interface Using Whole-Genome Sequencing in Austria

Antimicrobial resistance (AMR) is a public health issue attributed to the misuse of antibiotics in human and veterinary medicine. Since AMR surveillance requires a One Health approach, we sampled nine interconnected compartments at a hydrological open-air lab (HOAL) in Austria to obtain six bacterial species included in the WHO priority list of antibiotic-resistant bacteria (ARB). Whole genome sequencing-based typing included core genome multilocus sequence typing (cgMLST). Genetic and phenotypic characterization of AMR was performed for all isolates. Eighty-nine clinically-relevant bacteria were obtained from eight compartments including 49 E. coli, 27 E. faecalis, 7 K. pneumoniae and 6 E. faecium. Clusters of isolates from the same species obtained in different sample collection dates were detected. Of the isolates, 29.2% were resistant to at least one antimicrobial. E. coli and E. faecalis isolates from different compartments had acquired antimicrobial resistance genes (ARGs) associated with veterinary drugs such as aminoglycosides and tetracyclines, some of which were carried in conjugative and mobilizable plasmids. Three multidrug resistant (MDR) E. coli isolates were found in samples from field drainage and wastewater. Early detection of ARGs and ARB in natural and farm-related environments can identify hotspots of AMR and help prevent its emergence and dissemination along the food/feed chain.

The Population Genomics of Increased Virulence and Antibiotic Resistance in Human Commensal Escherichia coli over 30 Years in France

Escherichia coli is a commensal species of the lower intestine but is also a major pathogen causing intestinal and extraintestinal infections that is increasingly prevalent and resistant to antibiotics. Most studies on genomic evolution of E. coli used isolates from infections. Here, instead, we whole-genome sequenced a collection of 403 commensal E. coli isolates from fecal samples of healthy adult volunteers in France (1980 to 2010). These isolates were distributed mainly in phylogroups A and B2 (30% each) and belonged to 152 sequence types (STs), the five most frequent being ST10 (phylogroup A; 16.3%), ST73 and ST95 (phylogroup B2; 6.3 and 5.0%, respectively), ST69 (phylogroup D; 4.2%), and ST59 (phylogroup F; 3.9%), and 224 O:H serotypes. ST and serotype diversity increased over time. The O1, O2, O6, and O25 groups used in bioconjugate O-antigen vaccine against extraintestinal infections were found in 23% of the strains of our collection. The increase in frequency of virulence-associated genes and antibiotic resistance was driven by two evolutionary mechanisms. Evolution of virulence gene frequency was driven by both clonal expansion of STs with more virulence genes ("ST-driven") and increases in gene frequency within STs independent of changes in ST frequencies ("gene-driven"). In contrast, the evolution of resistance was dominated by increases in frequency within STs ("gene-driven"). This study provides a unique picture of the phylogenomic evolution of E. coli in its human commensal habitat over 30 years and will have implications for the development of preventive strategies. IMPORTANCE Escherichia coli is an opportunistic pathogen with the greatest burden of antibiotic resistance, one of the main causes of bacterial infections and an increasing concern in an aging population. Deciphering the evolutionary dynamics of virulence and antibiotic resistance in commensal E. coli is important to understand adaptation and anticipate future changes. The gut of vertebrates is the primary habitat of E. coli and probably where selection for virulence and resistance takes place. Unfortunately, most whole-genome-sequenced strains are isolated from pathogenic conditions. Here, we whole-genome sequenced 403 E. coli commensals isolated from healthy French subjects over a 30-year period. Virulence genes increased in frequency by both clonal expansion of clones carrying them and increases in frequency within clones, whereas resistance genes increased by within-clone increased frequency. Prospective studies of E. coli commensals should be performed worldwide to have a broader picture of evolution and adaptation of this species.

Convergence of virulence and antimicrobial resistance in increasingly prevalent Escherichia coli ST131 papGII+ sublineages

Escherichia coli lineage ST131 is an important cause of urinary tract and bloodstream infections worldwide and is highly resistant to antimicrobials. Specific ST131 lineages carrying invasiveness-associated papGII pathogenicity islands (PAIs) were previously described, but it is unknown how invasiveness relates to the acquisition of antimicrobial resistance (AMR). In this study, we analysed 1638 ST131 genomes and found that papGII+ isolates carry significantly more AMR genes than papGII-negative isolates, suggesting a convergence of virulence and AMR. The prevalence of papGII+ isolates among human clinical ST131 isolates increased dramatically since 2005, accounting for half of the recent E. coli bloodstream isolates. Emerging papGII+ lineages within clade C2 were characterized by a chromosomally integrated blaCTX-M-15 and the loss and replacement of F2:A1:B- plasmids. Convergence of virulence and AMR is worrying, and further dissemination of papGII+ ST131 lineages may lead to a rise in severe and difficult-to-treat extraintestinal infections.

Genomic analyses indicate that the ST131 lineage of E. coli, responsible for urinary and bloodstream infections globally, is evolving towards both increased virulence and increased resistance to antimicrobials.

Genomic and Temporal Trends in Canine ExPEC Reflect Those of Human ExPEC

Companion animals and humans are known to share extraintestinal pathogenic Escherichia coli (ExPEC), but the extent of E. coli sequence types (STs) that cause extraintestinal diseases in dogs is not well understood. Here, we generated whole-genome sequences of 377 ExPEC collected by the University of Melbourne Veterinary Hospital from dogs over an 11-year period from 2007 to 2017. Isolates were predominantly from urogenital tract infections (219, 58.1%), but isolates from gastrointestinal specimens (51, 13.5%), general infections (72, 19.1%), and soft tissue infections (34, 9%) were also represented. A diverse collection of 53 STs were identified, with 18 of these including at least five sequences. The five most prevalent STs were ST372 (69, 18.3%), ST73 (31, 8.2%), ST127 (22, 5.8%), ST80 (19, 5.0%), and ST58 (14, 3.7%). Apart from ST372, all of these are prominent human ExPEC STs. Other common ExPEC STs identified included ST12, ST131, ST95, ST141, ST963, ST1193, ST88, and ST38. Virulence gene profiles, antimicrobial resistance carriage, and trends in plasmid carriage for specific STs were generally reflective of those seen in humans. Many of the prominent STs were observed repetitively over an 11-year time span, indicating their persistence in the dogs in the community, which is most likely driven by household sharing of E. coli between humans and their pets. The case of ST372 as a dominant canine lineage observed sporadically in humans is flagged for further investigation. IMPORTANCE Pathogenic E. coli that causes extraintestinal infections (ExPEC) in humans and canines represents a significant burden in hospital and veterinary settings. Despite the obvious interrelationship between dogs and humans favoring both zoonotic and anthropozoonotic infections, whole-genome sequencing projects examining large numbers of canine-origin ExPEC are lacking. In support of anthropozoonosis, we found that most STs from canine infections are dominant human ExPEC STs (e.g., ST73, ST127, ST131) with similar genomic traits, such as plasmid carriage and virulence gene burden. In contrast, we identified ST372 as the dominant canine ST and a sporadic cause of infection in humans, supporting zoonotic transfer. Furthermore, we highlight that, as is the case in humans, STs in canine disease are consistent over time, implicating the gastrointestinal tract as the major community reservoir, which is likely augmented by exposure to human E. coli via shared diet and proximity.

Escherichia coli ST1193: Following in the Footsteps of E. coli ST131

Escherichia coli ST1193 is an emerging global multidrug (MDR) high-risk clone and an important cause of community-onset urinary and bloodstream infections. ST1193 is imitating E. coli ST131, the most successful MDR clone of all time. Both clones emerged in the early 1990s by acquiring quinolone resistance-determining region (QRDR) mutations, IncF plasmids, virulence factors, and type 1 pilus (fimH) recombination. They are the only MDR clones that are dominant among unselected E. coli populations. ST131 is the most frequent clone and ST1193 the second most frequent clone among fluoroquinolone/cephalosporin-resistant E. coli isolates. Both clones have played pivotal roles in the global spread of MDR E. coli. ST1193 originated from ST clonal complex 14 (STc14), is lactose nonfermenting, belongs to phylogenetic group B2, and contains the O type O75. Global ST1193 prevalence has been increasing since 2012, even replacing ST131 in certain regions. bla_CTX-M genes are rapidly expanding among ST1193 isolates, a scenario that occurred with ST131 during the 2000s. A validated PCR will enable global surveys to determine the extent of ST1193 among One Health E. coli isolates. The rapid emergence of ST1193 is concerning and is adding to the public health burden of MDR E. coli clones. Basic mechanistic, evolutionary, surveillance, and clinical studies are urgently required to investigate the success of ST1193. Such information will aid with management and prevention strategies. The medical community can ill afford to ignore the spread of another global successful MDR high-risk E. coli clone, especially one that is following in the footsteps of E. coli ST131.

Escherichia marmotae-a Human Pathogen Easily Misidentified as Escherichia coli

We hereby present the first descriptions of human-invasive infections caused by Escherichia marmotae, a recently described species that encompasses the former “Escherichia cryptic clade V.” We describe four cases, one acute sepsis of unknown origin, one postoperative sepsis after cholecystectomy, one spondylodiscitis, and one upper urinary tract infection. Cases were identified through unsystematic queries in a single clinical lab over 6 months. Through genome sequencing of the causative strains combined with available genomes from elsewhere, we demonstrate Es. marmotae to be a likely ubiquitous species containing genotypic virulence traits associated with Escherichia pathogenicity. The invasive isolates were scattered among isolates from a range of nonhuman sources in the phylogenetic analyses, thus indicating inherent virulence in multiple lineages. Pan genome analyses indicate that Es. marmotae has a large accessory genome and is likely to obtain ecologically advantageous traits, such as genes encoding antimicrobial resistance. Reliable identification might be possible by matrix-assisted laser desorption ionization–time of flight mass spectrometry (MALDI-TOF MS), but relevant spectra are missing in commercial databases. It can be identified through 16S rRNA gene sequencing. Escherichia marmotae could represent a relatively common human pathogen, and improved diagnostics will provide a better understanding of its clinical importance.

IMPORTANCEEscherichia coli is the most common pathogen found in blood cultures and urine and among the most important pathogenic species in the realm of human health. The notion that some of these isolates are not Es. coli but rather another species within the same genus may have implications for what Es. coli constitutes. We only recently have obtained methods to separate the two species, which means that possible differences in important clinical aspects, such as antimicrobial resistance rates, virulence, and phylogenetic structure, may exist. We believe that Es. marmotae as a common pathogen is new merely because we have not looked or bothered to distinguish between the thousands of invasive Escherichia passing through microbiological laboratories each day.

Genome wide association study of Escherichia coli bloodstream infection isolates identifies genetic determinants for the portal of entry but not fatal outcome

Escherichia coli is an important cause of bloodstream infections (BSI), which is of concern given its high mortality and increasing worldwide prevalence. Finding bacterial genetic variants that might contribute to patient death is of interest to better understand infection progression and implement diagnostic methods that specifically look for those factors. E. coli samples isolated from patients with BSI are an ideal dataset to systematically search for those variants, as long as the influence of host factors such as comorbidities are taken into account. Here we performed a genome-wide association study (GWAS) using data from 912 patients with E. coli BSI from hospitals in Paris, France. We looked for associations between bacterial genetic variants and three patient outcomes (death at 28 days, septic shock and admission to intensive care unit), as well as two portals of entry (urinary and digestive tract), using various clinical variables from each patient to account for host factors. We did not find any association between genetic variants and patient outcomes, potentially confirming the strong influence of host factors in influencing the course of BSI; we however found a strong association between the papGII operon and entrance of E. coli through the urinary tract, which demonstrates the power of bacterial GWAS when applied to actual clinical data. Despite the lack of associations between E. coli genetic variants and patient outcomes, we estimate that increasing the sample size by one order of magnitude could lead to the discovery of some putative causal variants. Given the wide adoption of bacterial genome sequencing of clinical isolates, such sample sizes may be soon available.

Reduced chlorhexidine susceptibility is associated with tetracycline resistance tet genes in clinical isolates of Escherichia coli

Chlorhexidine is a widely used antiseptic in hospital and community healthcare. Decreased susceptibility to this compound has been recently described in Klebsiella pneumoniae and Pseudomonas aeruginosa, together with cross-resistance to colistin. Surprisingly, few data are available for Escherichia coli, the main species responsible for community and healthcare-associated infections. In order to decipher chlorhexidine resistance mechanisms in E. coli, we studied both in vitro derived and clinical isolates through whole-genome sequence analysis. Comparison of strains grown in vitro under chlorhexidine pressure identified mutations in the gene mlaA coding for a phospholipid transport system. Phenotypic analyses of single-gene mutant from the Keio collection confirmed the role of this mutation in the decreased susceptibility to chlorhexidine. However, mutations in mlaA were not found in isolates from large clinical collections. In contrast, genome wide association studies (GWAS) showed that, in clinical strains, chlorhexidine reduced susceptibility was associated with the presence of tetA genes of class B coding for efflux pumps and located in a Tn10 transposon. Construction of recombinant strains in E. coli K-12 confirmed the role of tetA determinant in acquired resistance to both chlorhexidine and tetracycline. Our results reveal two different evolutionary paths leading to chlorhexidine decreased susceptibility: one restricted to in vitro evolution conditions and involving a retrograde phospholipid transport system; the other observed in clinical isolates associated with efflux pump TetA. None of these mechanisms provides cross-resistance to colistin. This work demonstrates the GWAS power to identify new resistance mechanisms in bacterial species.

A role for ColV plasmids in the evolution of pathogenic Escherichia coli ST58

Escherichia coli ST58 has recently emerged as a globally disseminated uropathogen that often progresses to sepsis. Unlike most pandemic extra-intestinal pathogenic E. coli (ExPEC), which belong to pathogenic phylogroup B2, ST58 belongs to the environmental/commensal phylogroup B1. Here, we present a pan-genomic analysis of a global collection of 752 ST58 isolates from diverse sources. We identify a large ST58 sub-lineage characterized by near ubiquitous carriage of ColV plasmids, which carry genes encoding virulence factors, and by a distinct accessory genome including genes typical of the Yersiniabactin High Pathogenicity Island. This sub-lineage includes three-quarters of all ExPEC sequences in our study and has a broad host range, although poultry and porcine sources predominate. By contrast, strains isolated from cattle often lack ColV plasmids. Our data indicate that ColV plasmid acquisition contributed to the divergence of the major ST58 sub-lineage, and different sub-lineages inhabit poultry, swine and cattle.

Recurring outbreaks by the same Escherichia coli ST10 clone in a broiler unit during 18 months

The current investigation aimed at characterizing the cause of multiple disease outbreaks in the same broiler production unit during a course of 18 months. The outbreaks had mortality rates of up to 22%. Escherichia coli was diagnosed as the responsible agent. Multiple-locus variable-number tandem-repeat analysis showed that all chicken isolates had identical band patterns. Core genome comparisons demonstrated that the 36 chicken isolates differed with maximum of nine nucleotides indicating that the same E. coli clone was responsible for all seven disease outbreaks despite adherence to the all-in-all production principle and rigorous cleaning and disinfection procedures.

Prevalence, Risk Factors, and Genetic Characterization of Extended-Spectrum Beta-Lactamase Escherichia coli Isolated From Healthy Pregnant Women in Madagascar

Antimicrobial resistance is a major public health concern worldwide affecting humans, animals and the environment. However, data is lacking especially in developing countries. Thus, the World Health Organization developed a One-Health surveillance project called Tricycle focusing on the prevalence of ESBL-producing Escherichia coli in humans, animals, and the environment. Here we present the first results of the human community component of Tricycle in Madagascar. From July 2018 to April 2019, rectal swabs from 492 pregnant women from Antananarivo, Mahajanga, Ambatondrazaka, and Toamasina were tested for ESBL-E. coli carriage. Demographic, sociological and environmental risk factors were investigated, and E. coli isolates were characterized (antibiotic susceptibility, resistance and virulence genes, plasmids, and genomic diversity). ESBL-E. coli prevalence carriage in pregnant women was 34% varying from 12% (Toamasina) to 65% (Ambatondrazaka). The main risk factor associated with ESBL-E. coli carriage was the rainy season (OR = 2.9, 95% CI 1.3-5.6, p = 0.009). Whole genome sequencing was performed on 168 isolates from 144 participants. bla CTX-M-15 was the most frequent ESBL gene (86%). One isolate was resistant to carbapenems and carried the bla NDM-5 gene. Most isolates belonged to commensalism associated phylogenetic groups A, B1, and C (90%) and marginally to extra-intestinal virulence associated phylogenetic groups B2, D and F (10%). Multi locus sequence typing showed 67 different sequence types gathered in 17 clonal complexes (STc), the most frequent being STc10/phylogroup A (35%), followed distantly by the emerging STc155/phylogroup B1 (7%), STc38/phylogroup D (4%) and STc131/phylogroup B2 (3%). While a wide diversity of clones has been observed, SNP analysis revealed several genetically close isolates (n = 34/168) which suggests human-to-human transmissions. IncY plasmids were found with an unusual prevalence (23%), all carrying a bla CTX-M-15. Most of them (85%) showed substantial homology (≥85%) suggesting a dissemination of IncY ESBL plasmids in Madagascar. This large-scale study reveals a high prevalence of ESBL-E. coli among pregnant women in four cities in Madagascar associated with warmth and rainfall. It shows the great diversity of E. coli disseminating throughout the country but also transmission of specific clones and spread of plasmids. This highlights the urgent need of public-health interventions to control antibiotic resistance in the country.

The E phylogroup of Escherichia coli is highly diverse and mimics the whole E. coli species population structure

To get a global picture of the population structure of the Escherichia coli phylogroup E, encompassing the O157:H7 EHEC lineage, we analysed the whole genome of 144 strains isolated from various continents, hosts and lifestyles and representative of the phylogroup diversity. The strains possess 4331 to 5440 genes with a core genome of 2771 genes and a pangenome of 33 722 genes. The distribution of these genes among the strains shows an asymmetric U-shaped distribution. E phylogenetic strains have the largest genomes of the species, partly explained by the presence of mobile genetic elements. Sixty-eight lineages were delineated, some of them exhibiting extra-intestinal virulence genes and being virulent in the mouse sepsis model. Except for the EHEC lineages and the reference EPEC, EIEC and ETEC strains, very few strains possess intestinal virulence genes. Most of the strains were devoid of acquired resistance genes, but eight strains possessed extended-spectrum beta-lactamase genes. Human strains belong to specific lineages, some of them being virulent and antibiotic-resistant [sequence type complexes (STcs) 350 and 2064]. The E phylogroup mimics all the features of the species as a whole, a phenomenon already observed at the STc level, arguing for a fractal population structure of E. coli.

Genome Informatics and Machine Learning-Based Identification of Antimicrobial Resistance-Encoding Features and Virulence Attributes in Escherichia coli Genomes Representing Globally Prevalent Lineages, Including High-Risk Clonal Complexes

Escherichia coli, a ubiquitous commensal/pathogenic member from the Enterobacteriaceae family, accounts for high infection burden, morbidity, and mortality throughout the world. With emerging multidrug resistance (MDR) on a massive scale, E. coli has been listed as one of the Global Antimicrobial Resistance and Use Surveillance System (GLASS) priority pathogens. Understanding the resistance mechanisms and underlying genomic features appears to be of utmost importance to tackle further spread of these multidrug-resistant superbugs. While a few of the globally prevalent sequence types (STs) of E. coli, such as ST131, ST69, ST405, and ST648, have been previously reported to be highly virulent and harboring MDR, there is no clarity if certain ST lineages have a greater propensity to acquire MDR. In this study, large-scale comparative genomics of a total of 5,653 E. coli genomes from 19 ST lineages revealed ST-wide prevalence patterns of genomic features, such as antimicrobial resistance (AMR)-encoding genes/mutations, virulence genes, integrons, and transposons. Interpretation of the importance of these features using a Random Forest Classifier trained with 11,988 genomic features from whole-genome sequence data identified ST-specific or phylogroup-specific signature proteins mostly belonging to different protein superfamilies, including the toxin-antitoxin systems. Our study provides a comprehensive understanding of a myriad of genomic features, ST-specific proteins, and resistance mechanisms entailing different lineages of E. coli at the level of genomes; this could be of significant downstream importance in understanding the mechanisms of AMR, in clinical discovery, in epidemiology, and in devising control strategies. IMPORTANCE With the leap in whole-genome data being generated, the application of relevant methods to mine biologically significant information from microbial genomes is of utmost importance to public health genomics. Machine-learning methods have been used not only to mine, curate, or classify the data but also to identify the relevant features that could be linked to a particular class/target. This is perhaps one of the pioneering studies that has attempted to classify a large repertoire of E. coli genome data sets (5,653 genomes) belonging to 19 different STs (including well-studied as well as understudied STs) using machine learning approaches. Important features identified by these approaches have revealed ST-specific signature proteins, which could be further studied to predict possible associations with the phenotypic profiles, thereby providing a better understanding of virulence and the resistance mechanisms among different clonal lineages of E. coli.