Systematic longitudinal survey of invasive Escherichia coli in England demonstrates a stable population structure only transiently disturbed by the emergence of ST131

Genomic epidemiology of third-generation cephalosporin-resistant Escherichia coli from companion animals and human infections in Europe

In high-income countries, dogs and cats are often considered members of the family. Because of this proximity, it has been suggested that pets and humans might exchange bacterial species from their gut microbiota, with multidrug resistant bacteria being of particular concern. The aim of this study was to compare the genomes of third-generation cephalosporin-resistant (3GC-R) Escherichia coli responsible for human and pet infections in Europe. Whole-genome sequencing data from 3GC-R E. coli isolated from clinical samples of humans, dogs and cats, and published in eight European studies were re-analyzed using bioinformatics tools. The acquired genes responsible for 3GC-R were identified. The sequence type (ST) of all genomes were assessed by multilocus sequence typing. Alpha and beta diversities were measured within and between the two populations. We included genomes of 1327 3GC-R E. coli isolated from humans and animals with 109 (8.2 %) being responsible for infections in dogs and cat, and 1218 (91.8 %) responsible for human infections. Alpha diversity analysis suggested greater diversity within ST and 3GC-R genes in the animal population. Beta diversity analysis by principal coordinate analysis separated animal and human strains. ST131 was more abundant in human strains (43.4 %) than in animal strains (14.7 %) (p < 0.001). Six STs, including ST372, were identified almost exclusively in 3GC-R E. coli from animal origin. The bla CTX-M-15 gene was more frequent in humans (49.24 %) than in companion animals (17.9 %) (p < 0.001). The resistance genes bla CMY-2 (30.8 %) and bla CTX-M-1 (15.4 %) were more frequent in E. coli isolated from pets (p < 0.001). We found that populations of 3GC-R E. coli responsible for human and pet infections in Europe do not overlap. Although it cannot rule out occasional transmission of bacteria between pets and humans within a household, it suggests that dogs and cats are not a major source of human infection with this antibiotic-resistant pathogen.

Geographical variation in the incidence of colorectal cancer and urinary tract cancer is associated with population exposure to colibactin-producing Escherichia coli

Biomedical research has implicated the bacterial metabolite colibactin as a causal risk factor for several cancer types, in particular, colorectal cancer. Colibactin has been known to drive tumorigenesis by inducing double-strand breaks in the DNA of epithelial cells exposed to colibactin-producing bacteria. Some phylogroup B2 Escherichia coli secrete colibactin during interbacterial warfare, concomitantly exposing the host to an increasing risk of DNA damage. This Personal View reviews the current knowledge about the cancer-colibactin interface and summarises metagenomics-based and population-genomics-based surveys to show that the prevalence of dominant colibactin-producing lineages of E coli varies considerably across geographical regions. The prevalence is further strongly associated with the age-standardised incidences of colorectal cancer, bladder cancer, and prostate cancer, suggesting that the degree of colibactin exposure in a population might contribute to the geographical variation of these cancers. Our observations provide a strong impetus for further research and the development of novel interventions to reduce the risks for colibactin-related cancers.

A comparative study of the molecular characteristics of human uropathogenic Escherichia coli collected from two hospitals in Norway and France in 2019

BACKGROUND

Escherichia coli is the primary aetiological agent of urinary tract infections (UTIs), a frequent cause of antibiotic prescription. Over the years, E. coli has become increasingly resistant towards several antimicrobial drugs used to treat UTIs.

OBJECTIVES

To compare the population structures and antimicrobial resistance profiles by molecular characterization of longitudinally collected E. coli from urine samples at two university hospitals in Norway and France.

METHOD

Each month in 2019, the first 10 E. coli identified in urine cultures of hospitalized and primary healthcare (PHC) patients were collected at both locations. WGS was performed to detect acquired antimicrobial resistance genes (ARGs) and plasmids and to determine phylogenetic relationships between isolates.

RESULTS

A total of 478 isolates were included, revealing a high diversity of STs among the isolates. However, the same 10 STs (ST73, ST131, ST69, ST95, ST12, ST141, ST127, ST10, ST58 and ST404) were dominant in both countries, constituting 55.9% (n = 267) of all isolates. Seven of the 10 prevalent STs belonged to phylogroup B2. The number of ARGs varied from 0 (n = 229) to 17 (n = 2). The mean ARG was higher in the French (2.9) than in the Norwegian cohort (2.1, P < 0.001), as well as in the hospitalized (3.0) compared with the PHC cohort (2.0, P = 0.007).

CONCLUSIONS

Despite the differences in patterns of ARGs between the Norwegian and French cohorts, the E. coli population structure was surprisingly highly conserved, suggesting that the distribution of the most common lineages largely depends on factors other than antibiotic use and resistance.

Genetic Evolution of Antibiotic Resistance and Virulence Genes in Escherichia coli Isolates from a Chinese Hospital over a 12-Year Period

Escherichia coli is a significant pathogen capable of inducing a variety of infections in both human and animal hosts. Despite its clinical significance, there is a lack of longitudinal research aimed at elucidating the genomic attributes that facilitate antimicrobial resistance and clonal dissemination in this bacterium. In this study, we investigated the genetic dynamics of antibiotic resistance and virulence factors within a collection of 137 E. coli isolates retrieved from a Chinese hospital over a 12-year period. Notably, a substantial increase in resistance to various antibiotics, including broad-spectrum β-lactams, aminoglycosides, and quinolones, was observed. Additionally, our study revealed the acquisition of diverse antibiotic resistance and virulence genes across different sequence types (STs). Among the STs, ST131 emerged as the most prevalent, exhibiting a high level of multidrug resistance. In contrast, ST73 and ST12 demonstrated a higher prevalence of virulence genes, suggestive of a potential trade-off between antibiotic resistance and virulence. What's more, we identified significant intra-clonal diversification and convergence of antibiotic resistance and virulence traits within the dominant ST131 group. These findings underscore the importance of longitudinal studies in understanding the evolution of bacterial pathogens and the necessity for ongoing research to inform public health strategies.

Plasmid-driven strategies for clone success in Escherichia coli

Escherichia coli is the most widely studied microbe in history, but the population structure and evolutionary trends of its extrachromosomal elements known as plasmids remain poorly delineated. Here we used long-read technology to high-resolution sequence the entire plasmidome and the corresponding host chromosomes from an unbiased longitudinal survey covering two decades and over 2000 E. coli isolates. We find that some plasmids have persisted in lineages even for centuries, demonstrating strong plasmid-lineage associations. Our analysis provides a detailed map of recent vertical and horizontal evolutionary events involving plasmids with key antibiotic resistance, competition and virulence determinants. We present genomic evidence of both chromosomal and plasmid-driven success strategies adopted by distant lineages by independently inheriting the same genomic elements. Further, we use in vitro experiments to verify the importance of key bacteriocin-producing plasmids for clone success. Our study has general implications for understanding plasmid biology and bacterial evolutionary strategies.

Bloodstream infections: mechanisms of pathogenesis and opportunities for intervention

Bloodstream infections (BSIs) are common in hospitals, often life-threatening and increasing in prevalence. Microorganisms in the blood are usually rapidly cleared by the immune system and filtering organs but, in some cases, they can cause an acute infection and trigger sepsis, a systemic response to infection that leads to circulatory collapse, multiorgan dysfunction and death. Most BSIs are caused by bacteria, although fungi also contribute to a substantial portion of cases. Escherichia coli, Staphylococcus aureus, coagulase-negative Staphylococcus, Klebsiella pneumoniae and Candida albicans are leading causes of BSIs, although their prevalence depends on patient demographics and geographical region. Each species is equipped with unique factors that aid in the colonization of initial sites and dissemination and survival in the blood, and these factors represent potential opportunities for interventions. As many pathogens become increasingly resistant to antimicrobials, new approaches to diagnose and treat BSIs at all stages of infection are urgently needed. In this Review, we explore the prevalence of major BSI pathogens, prominent mechanisms of BSI pathogenesis, opportunities for prevention and diagnosis, and treatment options.

Major F plasmid clusters are linked with ColV and pUTI89-like marker genes in bloodstream isolates of Escherichia coli

BACKGROUND

F plasmids are abundant in E. coli, carrying a variety of genetic cargo involved in fitness, pathogenicity, and antimicrobial resistance. ColV and pUTI89-like plasmids have drawn attention for their potential roles in various forms of extra-intestinal pathogenicity. However, the rates of their carriage and the overall diversity of F plasmids in E. coli bloodstream infections (BSI E. coli) remain unknown.

METHODS

We performed a t-SNE-based cluster analysis of predicted F plasmids from a collection of 4711 BSI E. coli draft genomes to describe their diversity and abundance. We also screened them for markers of ColV and pUTI89-like plasmids, F plasmid replicon sequence types (RST) and E. coli sequence types (ST) to understand how genetic features were related to plasmid clusters.

RESULTS

Predicted F plasmids in BSI E. coli draft genomes were embedded within five major clusters based on a model of complete F plasmid sequences. Nearly half of the clustered sequences belonged to two major clusters, which were associated with ColV and pUTI89-like marker genes, respectively. Genomes from the ColV cluster featured F2:A-:B1 and F24:A-B1 RSTs in association with ST95, ST58 and ST88, whilst the pUTI89-like cluster was mostly F29:A-:B10 linked to ST73, ST69, ST95 and ST131. Plasmids associated with different lineages of ST131 formed additional major clusters, whilst F51:A-:B10 plasmids in ST73 were also common.

CONCLUSIONS

ColV and pUTI89-like plasmid markers are predominant in BSI E. coli that carry F plasmids. These markers are associated with distinct clusters of plasmids across diverse sequence types of E. coli. We hypothesise that their abundance in BSI E. coli is partially driven by carriage of backbone genes previously shown to contribute to virulence in models of bloodstream infection. Their carriage by pandemic E. coli STs suggests clonal expansion also plays a role in their success in BSI. Ecological pathways via which these plasmids evolve, and spread are likely to be distinct as other studies show ColV is strongly associated with poultry and food animal production, whereas pUTI89-like plasmids appear to be mostly human-restricted.

Quantifying the Evolutionary Dynamics of Structure and Content in Closely Related E. coli Genomes

Bacterial genomes primarily diversify via gain, loss, and rearrangement of genetic material in their flexible accessory genome. Yet the dynamics of accessory genome evolution are very poorly understood, in contrast to the core genome where diversification is readily described by mutations and homologous recombination. Here, we tackle this problem for the case of very closely related genomes. We comprehensively describe genome evolution within n=222 genomes of Escherichia coli ST131, which likely shared a common ancestor around 100 years ago. After removing putative recombinant diversity, the total length of the phylogeny is 6,000 core genome substitutions. Within this diversity, we find 22 modifications to core genome synteny and estimate around 2,000 structural changes within the accessory genome, i.e. one structural change for every three core genome substitutions. Sixty-three percent of loci with structural diversity could be resolved into individual gain and loss events with 10-fold more gains than losses, demonstrating a dominance of gains due to insertion sequences and prophage integration. Our results suggest the majority of synteny changes and insertions in our dataset are likely deleterious and only persist for a short time before being removed by purifying selection.

Genomic and growth fitness study of extended-spectrum β-lactamase-producing Escherichia coli from bloodstream infections after introduction of a national 4C antimicrobial stewardship policy in Scotland

BACKGROUND

Extended-spectrum β-lactamase-producing Escherichia coli remains a major cause of hospital-acquired bloodstream infections in countries with high antimicrobial stewardship compliance.

METHODS

Isolates from bloodstream infections that occurred between 2010 and 2020 in a tertiary-level hospital in North-East Scotland soon after introduction of the '4C' antimicrobial stewardship policy were analysed for phylogenetic structure, antimicrobial resistance, plasmid, and virulence gene carriage. Growth fitness was measured in kinetic assays. Non-metric-multidimensional-scaling was used to evaluate clonal relationships, antimicrobial resistance, and virulence profiles in early and later years after the 4C policy introduction. Clonal and fitness trends over the study period were determined by generalised additive modelling. The relationship between clonal type, antimicrobial resistance, and fitness was evaluated by linear regression.

RESULTS

Three hierarchical phylogenetic clusters were identified, with the most dominant cluster (O25:H4/fimH30) including all, and nearly exclusively, Clade C ST131 isolates as well as minor non-ST131 sequence types. The prevalence of ST131 was largely stable over the study period. Resistance to aminoglycosides and aztreonam in ST131 was lower (P = 0.019 and P = 0.004, respectively) during later years (2016-2020) by 28% on average compared to early years soon after 4C policy implementation (2010-2014). Carriage of virulence factors involved in bacterial adaptation was higher in ST131 compared to non-ST131 but mostly stable in early vs. later years. Growth fitness of ST131 was lower than non-ST131 and declined steadily in later years (P < 0.0001).

CONCLUSIONS

Despite stable virulence factor carriage, population structure, and resistance to cephalosporins, we show increased susceptibility of ST131 to aminoglycosides and aztreonam and concurrent fitness decline years after the introduction of the 4C policy.

Comparative genomic analysis of ESBL-selected and non-selected Escherichia coli in Australian wastewater: Elucidating differences in diversity, antimicrobial resistance, and virulence profiles

Extended-spectrum β-lactamases (ESBLs)-producing E. coli have been proposed as an indicator bacterium for antimicrobial resistance (AMR) surveillance within a OneHealth framework. However, it is important to understand the effects and potential biases ESBL-selection has on E. coli populations. Utilising whole genome sequencing, this study compared 80 ESBL-selected E. coli isolates with 201 non-selected isolates from Australian wastewater. The findings revealed significant variations between these cohorts in genetic diversity, AMR profiles, and carriage of virulence-associated genes (VAGs), plasmids, and the transmissible Locus of Stress Tolerance (tLST), a genomic island that imparts resistance to extreme heat and chlorination. The study highlights the predominance of certain sequence types (STs), particularly ST131 (75 % clade A), in ESBL-selected isolates (40 % vs 2 %) and overall the ESBL-selected isolates were largely multidrug-resistant (MDR), predominantly carrying genes for resistance to aminoglycosides, extended-spectrum β-lactams, fluoroquinolone, macrolides, sulphonamides/trimethoprim, and tetracyclines. The ESBLs identified were almost exclusively blaCTX-M genes, most commonly blaCTX-M-15 > blaCTX-M-27 > blaCTX-M-14. These were predominately carried on IncF plasmids or chromosomally (always ISEcp1 associated), in equal numbers. In contrast, 80 % of non-selected isolates carried no acquired ARGs, and none carried blaCTX-M genes. In both cohorts, extraintestinal pathogenic E. coli (ExPEC) was the dominate pathotype (35 % total) with few (4 % total) intestinal pathogenic E. coli pathotypes identified (aEPEC > ETEC > EAEC). Nevertheless, some clinically important genes were only identified in the non-selected group, namely tigecycline-resistance gene tet(X4) and AmpC ESBL blaCMY-2. Additionally, the presence of tLST, associated with higher metal resistance gene carriage (Ag, As, Cu, Hg, Ni), in a substantial portion of non-selected isolates (20 % vs 0 %), underscores environmental pressures shaping bacterial populations in wastewater ecosystems. These insights are important for developing comprehensive, less biased genomic surveillance strategies to understand and manage public health threats posed by pathogenic E. coli and AMR.

Exploring uncatalogued genetic variation in antimicrobial resistance gene families in Escherichia coli: an observational analysis

BACKGROUND

Antimicrobial resistance (AMR) in Escherichia coli is a global problem associated with substantial morbidity and mortality. AMR-associated genes are typically annotated based on similarity to variants in a curated reference database, with the implicit assumption that uncatalogued genetic variation within these is phenotypically unimportant. In this study, we evaluated the performance of the AMRFinder tool and, subsequently, the potential for discovering new AMR-associated gene families and characterising variation within existing ones to improve genotype-to-susceptibility phenotype predictions in E coli.

METHODS

In this cross-sectional study of international genome sequence data, we assembled a global dataset of 9001 E coli sequences from five publicly available data collections predominantly deriving from human bloodstream infections from: Norway, Oxfordshire (UK), Thailand, the UK, and Sweden. 8555 of these sequences had linked antibiotic susceptibility data. Raw reads were assembled using Shovill and AMR genes (relevant to amoxicillin-clavulanic acid, ampicillin, ceftriaxone, ciprofloxacin, gentamicin, piperacillin-tazobactam, and trimethoprim) extracted using the National Center for Biotechnology Information AMRFinder tool (using both default and strict [100%] coverage and identity filters). We assessed the predictive value of the presence of these genes for predicting resistance or susceptibility against US Food and Drug Administration thresholds for major and very major errors. Mash was used to calculate the similarity between extracted genes using Jaccard distances. We empirically reclustered extracted gene sequences into AMR-associated gene families (≥70% match) and antibiotic-resistance genes (ARGs; 100% match) and categorised these according to their frequency in the dataset. Accumulation curves were simulated and correlations between gene frequency in the Oxfordshire and other datasets calculated using the Spearman coefficient. Firth regression was used to model the association between the presence of blaTEM-1 variants and amoxicillin-clavulanic acid or piperacillin-tazobactam resistance, adjusted for the presence of other relevant ARGs.

FINDINGS

The performance of the AMRFinder database for genotype-to-phenotype predictions using strict 100% identity and coverage thresholds did not meet US Food and Drug Administration thresholds for any of the seven antibiotics evaluated. Relaxing filters to default settings improved sensitivity with a specificity cost. For all antibiotics, most explainable resistance was associated with the presence of a small number of genes. There was a proportion of resistance that could not be explained by known ARGs; this ranged from 75·1% for amoxicillin-clavulanic acid to 3·4% for ciprofloxacin. Only 18 199 (51·5%) of the 35 343 ARGs detected had a 100% identity and coverage match in the AMRFinder database. After empirically reclassifying genes at 100% nucleotide sequence identity, we identified 1042 unique ARGs, of which 126 (12·1%) were present ten times or more, 313 (30·0%) were present between two and nine times, and 603 (57·9%) were present only once. Simulated accumulation curves revealed that discovery of new (100% match) ARGs present more than once in the dataset plateaued relatively quickly, whereas new singleton ARGs were discovered even after many thousands of isolates had been included. We identified a strong correlation (Spearman coefficient 0·76 [95% CI 0·73-0·80], p<0·0001) between the number of times an ARG was observed in Oxfordshire and the number of times it was seen internationally, with ARGs that were observed six times in Oxfordshire always being found elsewhere. Finally, using the example of blaTEM-1, we showed that uncatalogued variation, including synonymous variation, is associated with potentially important phenotypic differences; for example, two common, uncatalogued blaTEM-1 alleles with only synonymous mutations compared with the known reference were associated with reduced resistance to amoxicillin-clavulanic acid (adjusted odds ratio 0·58 [95% CI 0·35-0·95], p=0·031) and piperacillin-tazobactam (0·50 [95% CI 0·29-0·82], p=0·005).

INTERPRETATION

We highlight substantial uncatalogued genetic variation with respect to known ARGs, although a relatively small proportion of these alleles are repeatedly observed in a large international dataset suggesting strong selection pressures. The current approach of using fuzzy matching for ARG detection, ignoring the unknown effects of uncatalogued variation, is unlikely to be acceptable for future clinical deployment. The association of synonymous mutations with potentially important phenotypic differences suggests that relying solely on amino acid-based gene detection to predict resistance is unlikely to be sufficient. Finally, the inability to explain all resistance using existing knowledge highlights the importance of new target gene discovery.

FUNDING

National Institute for Health and Care Research, Wellcome, and UK Medical Research Council.

The complete genome sequence of five pre-2013 Escherichia coli sequence type (ST)1193 strains reveals insights into an emerging pathogen

Fluoroquinolone-resistant Escherichia coli sequence type (ST)1193 is a profound, emerging lineage associated with systemic, urinary tract and neonatal infections. Humans, companion animals and the environment are reservoirs for ST1193, which has been disseminated globally. Following its detection in 2007, ST1193 has been identified repeatedly amongst fluoroquinolone-resistant clones in Australia. However, despite the growing importance of ST1193, only three complete genomes are published in the literature, none of which are from Australia. Here we expand on the available ST1193 resources with the complete genomes of five ST1193 strains sequenced using Oxford Nanopore Technologies and Illumina. Using in silico genotyping, we found that all strains were multi-drug resistant, including resistances to fluoroquinolones and cephalosporins. In vitro antibiotic susceptibility testing mostly correlated with individual genotypes. The exception was MS8320, which had additional in vitro resistance to piperacillin/tazobactam, ampicillin/sulbactam, cefazolin and doripenem (carbapenem). Further investigation identified seven additional copies of an IS26 transposable unit carrying a bla TEM-1B beta-lactamase gene, suggesting this tandem amplification is associated with extended resistance phenotypes. Uropathogenicity factors, including three separate siderophore-encoding loci, were conserved in chromosomal and plasmid regions. Using all complete genomes, we further elucidated the recombination events surrounding the previously described K5/K1 capsular locus switch. Phenotypic confirmation of differing capsules in Australian ST1193 strains, coupled with genetic analysis revealing insertions downstream of the capsular locus, underscored the genetic distinctions between K5 and K1 capsule encoding strains. This study provides five new reference ST1193 genomes from Australia. These include the earliest complete K5-capsule ST1193 genomes on record (collected 2007), alongside our reference genome (MS10858), a clinical isolate obtained early during the ST1193 expansion and representative of the predominant K1-associated clade. These findings lay the foundations for further genomic and molecular analyses that may help understand the underlying reasons for the rapid global expansion of ST1193.

PanKA: Leveraging population pangenome to predict antibiotic resistance

Machine learning has the potential to be a powerful tool in the fight against antimicrobial resistance (AMR), a critical global health issue. Machine learning can identify resistance mechanisms from DNA sequence data without prior knowledge. The first step in building a machine learning model is a feature extraction from sequencing data. Traditional methods like single nucleotide polymorphism (SNP) calling and k-mer counting yield numerous, often redundant features, complicating prediction and analysis. In this paper, we propose PanKA, a method using the pangenome to extract a concise set of relevant features for predicting AMR. PanKA not only enables fast model training and prediction but also improves accuracy. Applied to the Escherichia coli and Klebsiella pneumoniae bacterial species, our model is more accurate than conventional and state-of-the-art methods in predicting AMR.

A microbiological and genomic perspective of globally collected Escherichia coli from adults hospitalized with invasive E. coli disease

OBJECTIVES

Escherichia coli can cause infections in the urinary tract and in normally sterile body sites leading to invasive E. coli disease (IED), including bacteraemia and sepsis, with older populations at increased risk. We aimed to estimate the theoretical coverage rate by the ExPEC4V and 9V vaccine candidates. In addition, we aimed at better understanding the diversity of E. coli isolates, including their genetic and phenotypic antimicrobial resistance (AMR), sequence types (STs), O-serotypes and the bacterial population structure.

METHODS

Blood and urine culture E. coli isolates (n = 304) were collected from hospitalized patients ≥60 years (n = 238) with IED during a multicentric, observational study across three continents. All isolates were tested for antimicrobial susceptibility, O-serotyped, whole-genome sequenced and bioinformatically analysed.

RESULTS

A large diversity of STs and of O-serotypes were identified across all centres, with O25b-ST131, O6-ST73 and O1-ST95 being the most prevalent types. A total of 45.4% and 64.7% of all isolates were found to have an O-serotype covered by the ExPEC4V and ExPEC9V vaccine candidates, respectively. The overall frequency of MDR was 37.4% and ST131 was predominant among MDR isolates. Low in-patient genetic variability was observed in cases where multiple isolates were collected from the same patient.

CONCLUSIONS

Our results highlight the predominance of MDR O25b-ST131 E. coli isolates across diverse geographic areas. These findings provide further baseline data on the theoretical coverage of novel vaccines targeting E. coli associated with IED in older adults and their associated AMR levels.

CELEBRIMBOR: core and accessory genes from metagenomes

MOTIVATION

Metagenome-Assembled Genomes (MAGs) or Single-cell Amplified Genomes (SAGs) are often incomplete, with sequences missing due to errors in assembly or low coverage. This presents a particular challenge for the identification of true gene frequencies within a microbial population, as core genes missing in only a few assemblies will be mischaracterized by current pangenome approaches.

RESULTS

Here, we present CELEBRIMBOR, a Snakemake pangenome analysis pipeline which uses a measure of genome completeness to automatically adjust the frequency threshold at which core genes are identified, enabling accurate core gene identification in MAGs and SAGs.

AVAILABILITY AND IMPLEMENTATION

CELEBRIMBOR is published under open source Apache 2.0 licence at https://github.com/bacpop/CELEBRIMBOR and is available as a Docker container from this repository. Supplementary material is available in the online version of the article.

Evolution of extended-spectrum β-lactamase-producing ST131 Escherichia coli at a single hospital over 15 years

Escherichia coli multi-locus sequence type ST131 is a globally distributed pandemic lineage that causes multidrug-resistant extra-intestinal infections. ST131 E. coli frequently produce extended-spectrum β-lactamases (ESBLs), which confer resistance to many β-lactam antibiotics and make infections difficult to treat. We sequenced the genomes of 154 ESBL-producing E. coli clinical isolates belonging to the ST131 lineage from patients at the University of Pittsburgh Medical Center (UPMC) between 2004 and 2018. Isolates belonged to the well described ST131 clades A (8%), B (3%), and C (89%). Time-dated phylogenetic analysis estimated that the most recent common ancestor (MRCA) for all clade C isolates emerged around 1989, consistent with previous studies. We identified multiple genes potentially under selection in clade C, including the cell wall assembly gene ftsI, the LPS biosynthesis gene arnC, and the yersiniabactin uptake receptor fyuA. Diverse ESBL-encoding genes belonging to the blaCTX-M, blaSHV, and blaTEM families were identified; these genes were found at varying numbers of loci and in variable numbers of copies across isolates. Analysis of ESBL flanking regions revealed diverse mobile elements that varied by ESBL type. Overall, our findings show that ST131 subclade C dominated among patients and uncover possible signals of ongoing adaptation within this ST131 lineage.

Population structure and antibiotic resistance of swine extraintestinal pathogenic Escherichia coli from China

Extraintestinal Pathogenic Escherichia coli (ExPEC) pose a significant threat to human and animal health. However, the diversity and antibiotic resistance of animal ExPEC, and their connection to human infections, remain largely unexplored. The study performs large-scale genome sequencing and antibiotic resistance testing of 499 swine-derived ExPEC isolates from China. Results show swine ExPEC are phylogenetically diverse, with over 80% belonging to phylogroups B1 and A. Importantly, 15 swine ExPEC isolates exhibit genetic relatedness to human-origin E. coli strains. Additionally, 49 strains harbor toxins typical of enteric E. coli pathotypes, implying hybrid pathotypes. Notably, 97% of the total strains are multidrug resistant, including resistance to critical human drugs like third- and fourth-generation cephalosporins. Correspondingly, genomic analysis unveils prevalent antibiotic resistance genes (ARGs), often associated with co-transfer mechanisms. Furthermore, analysis of 20 complete genomes illuminates the transmission pathways of ARGs within swine ExPEC and to human pathogens. For example, the transmission of plasmids co-harboring fosA3, blaCTX-M-14, and mcr-1 genes between swine ExPEC and human-origin Salmonella enterica is observed. These findings underscore the importance of monitoring and controlling ExPEC infections in animals, as they can serve as a reservoir of ARGs with the potential to affect human health or even be the origin of pathogens infecting humans.

Deep sequencing of Escherichia coli exposes colonisation diversity and impact of antibiotics in Punjab, Pakistan

Multi-drug resistant (MDR) E. coli constitute a major public health burden globally, reaching the highest prevalence in the global south yet frequently flowing with travellers to other regions. However, our comprehension of the entire genetic diversity of E. coli colonising local populations remains limited. We quantified this diversity, its associated antimicrobial resistance (AMR), and assessed the impact of antibiotic use by recruiting 494 outpatients and 423 community dwellers in the Punjab province, Pakistan. Rectal swab and stool samples were cultured on CLED agar and DNA extracted from plate sweeps was sequenced en masse to capture both the genetic and AMR diversity of E. coli. We assembled 5,247 E. coli genomes from 1,411 samples, displaying marked genetic diversity in gut colonisation. Compared with high income countries, the Punjabi population generally showed a markedly different distribution of genetic lineages and AMR determinants, while use of antibiotics elevated the prevalence of well-known globally circulating MDR clinical strains. These findings implicate that longitudinal multi-regional genomics-based surveillance of both colonisation and infections is a prerequisite for developing mechanistic understanding of the interplay between ecology and evolution in the maintenance and dissemination of (MDR) E. coli.

Detecting co-selection through excess linkage disequilibrium in bacterial genomes

Population genomics has revolutionized our ability to study bacterial evolution by enabling data-driven discovery of the genetic architecture of trait variation. Genome-wide association studies (GWAS) have more recently become accompanied by genome-wide epistasis and co-selection (GWES) analysis, which offers a phenotype-free approach to generating hypotheses about selective processes that simultaneously impact multiple loci across the genome. However, existing GWES methods only consider associations between distant pairs of loci within the genome due to the strong impact of linkage-disequilibrium (LD) over short distances. Based on the general functional organisation of genomes it is nevertheless expected that majority of co-selection and epistasis will act within relatively short genomic proximity, on co-variation occurring within genes and their promoter regions, and within operons. Here, we introduce LDWeaver, which enables an exhaustive GWES across both short- and long-range LD, to disentangle likely neutral co-variation from selection. We demonstrate the ability of LDWeaver to efficiently generate hypotheses about co-selection using large genomic surveys of multiple major human bacterial pathogen species and validate several findings using functional annotation and phenotypic measurements. Our approach will facilitate the study of bacterial evolution in the light of rapidly expanding population genomic data.

High-risk Escherichia coli clones that cause neonatal meningitis and association with recrudescent infection

Neonatal meningitis is a devastating disease associated with high mortality and neurological sequelae. Escherichia coli is the second most common cause of neonatal meningitis in full-term infants (herein NMEC) and the most common cause of meningitis in preterm neonates. Here, we investigated the genomic relatedness of a collection of 58 NMEC isolates spanning 1974-2020 and isolated from seven different geographic regions. We show NMEC are comprised of diverse sequence types (STs), with ST95 (34.5%) and ST1193 (15.5%) the most common. No single virulence gene profile was conserved in all isolates; however, genes encoding fimbrial adhesins, iron acquisition systems, the K1 capsule, and O antigen types O18, O75, and O2 were most prevalent. Antibiotic resistance genes occurred infrequently in our collection. We also monitored the infection dynamics in three patients that suffered recrudescent invasive infection caused by the original infecting isolate despite appropriate antibiotic treatment based on antibiogram profile and resistance genotype. These patients exhibited severe gut dysbiosis. In one patient, the causative NMEC isolate was also detected in the fecal flora at the time of the second infection episode and after treatment. Thus, although antibiotics are the standard of care for NMEC treatment, our data suggest that failure to eliminate the causative NMEC that resides intestinally can lead to the existence of a refractory reservoir that may seed recrudescent infection.

Whole-genome characterisation of Escherichia coli isolates from patients with bacteraemia presenting with sepsis or septic shock in Spain: a multicentre cross-sectional study

BACKGROUND

Escherichia coli is the most frequent cause of bloodstream infections (BSIs). About one-third of patients with BSIs due to E coli develop sepsis or shock. The objective of this study is to characterise the microbiological features of E coli blood isolates causing sepsis or septic shock to provide exploratory information for future diagnostic, preventive, or therapeutic interventions.

METHODS

E coli blood isolates from a multicentre cross-sectional study of patients older than 14 years presenting with sepsis or septic shock (according to the Third International Consensus Definitions for Sepsis and Septic Shock criteria) from hospitals in Spain between Oct 4, 2016, and Oct 15, 2017, were studied by whole-genome sequencing. Phylogroups, sequence types (STs), serotype, FimH types, antimicrobial resistance (AMR) genes, pathogenicity islands, and virulence factors were identified. Susceptibility testing was performed by broth microdilution. The main outcome of this study was the characterisation of the E coli blood isolates in terms of population structure by phylogroups, groups (group 1: phylogroups B2, F, and G; group 2: A, B1, and C; group 3: D), and STs and distribution by geographical location and bloodstream infection source. Other outcomes were virulence score and prevalence of virulence-associated genes, pathogenicity islands, AMR, and AMR-associated genes. Frequencies were compared using χ² or Fisher's exact tests, and continuous variables using the Mann-Whitney test, with Bonferroni correction for multiple comparisons.

FINDINGS

We analysed 224 isolates: 140 isolates (63%) were included in phylogenetic group 1, 52 (23%) in group 2, and 32 (14%) in group 3. 85 STs were identified, with four comprising 44% (n=98) of the isolates: ST131 (38 [17%]), ST73 (25 [11%]), ST69 (23 [10%]), and ST95 (12 [5%]). No significant differences in phylogroup or ST distribution were found according to geographical areas or source of bloodstream infection, except for ST95, which was more frequent in urinary tract infections than in other sources (11 [9%] of 116 vs 1 [1%] of 108, p=0·0045). Median virulence score was higher in group 1 (median 25·0 [IQR 20·5-29·0) than in group 2 (median 14·5 [9·0-20·0]; p<0·0001) and group 3 (median 21 [16·5-23·0]; p<0·0001); prevalence of several pathogenicity islands was higher in group 1. No significant differences were found between phylogenetic groups in proportions of resistance to antibiotics. ST73 had higher median virulence score (32 [IQR 29-35]) than the other predominant clones (median range 21-28). Some virulence genes and pathogenicity islands were significantly associated with each ST. ST131 isolates had higher prevalence of AMR and a higher proportion of AMR genes, notably blaCTX-M-15 and blaOXA-1.

INTERPRETATION

In this exploratory study, the population structure of E coli causing sepsis or shock was similar to previous studies that included all bacteraemic isolates. Virulence genes, pathogenicity islands, and AMR genes were not randomly distributed among phylogroups or STs. These results provide a comprehensive characterisation of invasive E coli isolates causing severe response syndrome. Future studies are required to determine the contribution of these microbiological factors to severe clinical presentation and worse outcomes in patients with E coli bloodstream infection.

FUNDING

Instituto de Salud Carlos III.

Anthropogenic contamination sources drive differences in antimicrobial-resistant Escherichia coli in three urban lakes

Antimicrobial resistance (AMR) is an ever-present threat to the treatment of infectious diseases. However, the potential relevance of this phenomenon in environmental reservoirs still raises many questions. Detection of antimicrobial-resistant bacteria in the environment is a critical aspect for understanding the prevalence of resistance outside of clinical settings, as detection in the environment indicates that resistance is likely already widespread. We isolated antimicrobial-resistant Escherichia coli from three urban waterbodies over a 15-month time series, determined their antimicrobial susceptibilities, investigated their population structure, and identified genetic determinants of resistance. We found that E. coli populations at each site were composed of different dominant phylotypes and showed distinct patterns of antimicrobial and multidrug resistance, despite close geographic proximity. Many strains that were genome-sequenced belonged to sequence types of international concern, particularly the ST131 clonal complex. We found widespread resistance to clinically important antimicrobials such as amoxicillin, cefotaxime, and ciprofloxacin, but found that all strains were susceptible to amikacin and the last-line antimicrobials meropenem and fosfomycin. Resistance was most often due to acquirable antimicrobial resistance genes, while chromosomal mutations in gyrA, parC, and parE conferred resistance to quinolones. Whole-genome analysis of a subset of strains further revealed the diversity of the population of E. coli present, with a wide array of AMR and virulence genes identified, many of which were present on the chromosome, including blaCTX-M. Finally, we determined that environmental persistence, transmission between sites, most likely mediated by wild birds, and transfer of mobile genetic elements likely contributed significantly to the patterns observed.IMPORTANCEA One Health perspective is crucial to understand the extent of antimicrobial resistance (AMR) globally, and investigation of AMR in the environment has been increasing in recent years. However, most studies have focused on waterways that are directly polluted by sewage, industrial manufacturing, or agricultural activities. Therefore, there remains a lack of knowledge about more natural, less overtly impacted environments. Through phenotypic and genotypic investigation of AMR in Escherichia coli, this study adds to our understanding of the extent and patterns of resistance in these types of environments, including over a time series, and showed that complex biotic and abiotic factors contribute to the patterns observed. Our study further emphasizes the importance of incorporating the surveillance of microbes in freshwater environments in order to better comprehend potential risks for both human and animal health and how the environment may serve as a sentinel for potential future clinical infections.

Generalizability of machine learning in predicting antimicrobial resistance in E. coli: a multi-country case study in Africa

BACKGROUND

Antimicrobial resistance (AMR) remains a significant global health threat particularly impacting low- and middle-income countries (LMICs). These regions often grapple with limited healthcare resources and access to advanced diagnostic tools. Consequently, there is a pressing need for innovative approaches that can enhance AMR surveillance and management. Machine learning (ML) though underutilized in these settings, presents a promising avenue. This study leverages ML models trained on whole-genome sequencing data from England, where such data is more readily available, to predict AMR in E. coli, targeting key antibiotics such as ciprofloxacin, ampicillin, and cefotaxime. A crucial part of our work involved the validation of these models using an independent dataset from Africa, specifically from Uganda, Nigeria, and Tanzania, to ascertain their applicability and effectiveness in LMICs.

RESULTS

Model performance varied across antibiotics. The Support Vector Machine excelled in predicting ciprofloxacin resistance (87% accuracy, F1 Score: 0.57), Light Gradient Boosting Machine for cefotaxime (92% accuracy, F1 Score: 0.42), and Gradient Boosting for ampicillin (58% accuracy, F1 Score: 0.66). In validation with data from Africa, Logistic Regression showed high accuracy for ampicillin (94%, F1 Score: 0.97), while Random Forest and Light Gradient Boosting Machine were effective for ciprofloxacin (50% accuracy, F1 Score: 0.56) and cefotaxime (45% accuracy, F1 Score:0.54), respectively. Key mutations associated with AMR were identified for these antibiotics.

CONCLUSION

As the threat of AMR continues to rise, the successful application of these models, particularly on genomic datasets from LMICs, signals a promising avenue for improving AMR prediction to support large AMR surveillance programs. This work thus not only expands our current understanding of the genetic underpinnings of AMR but also provides a robust methodological framework that can guide future research and applications in the fight against AMR.

Detection and characterization of potentially hybrid enteroaggregative Escherichia coli (EAEC) strains isolated from urinary tract infection

Uropathogenic Escherichia coli (UPEC) have the potential to receive the virulence markers of intestinal pathotypes and transform into various important hybrid pathotypes. This study aimed to investigate the frequency and characteristics of hybrid enteroaggregative E. coli (EAEC)/UPEC strains. Out of 202 UPEC strains, nine (4.5%) were detected as hybrid EAEC/UPEC. These strains carried one to four iron uptake systems. Among nine investigated pathogenicity islands (PAIs), PAI IV536, PAI II536, and PAI ICFT073 were found in 9 (100%), 3 (33.3%), and 1 (11.1%) strains, respectively. The chuA and sitA genes were detected in 5 (55.5%) and 3 (33.3%) hybrid strains, respectively. Six hybrid strains were found to be typical extraintestinal pathogenic E. coli (ExPEC) according to their virulence traits. Most of the hybrid strains belonged to the phylogenetic group E (6/9). Among the hybrid strains, seven (7/9) were able to form biofilm and adhere to cells; however, only two strains penetrated into the HeLa cells. Our findings reveal some of the virulence characteristics of hybrid strains that lead to fitness and infection in the urinary tract. These strains, with virulence factors of intestinal and non-intestinal pathotypes, may become emerging pathogens in clinical settings; therefore, further studies are needed to reveal their pathogenicity mechanisms and so that preventive measures can be taken.

Emergence and dissemination of epidemic-causing OXA-244 carbapenemase-producing Escherichia coli ST38 through hospital sewage in Norway, 2020-2022

BACKGROUND

Population-based sewage surveillance has emerged as a promising approach for studying the prevalence of antibiotic resistance in pathogens.

AIM

To determine the temporal prevalence of cefotaxime-resistant Escherichia coli in sewage from five sewage treatment plants located in Bergen city, to determine whether ESBL- and carbapenemase-producing E. coli are consistently disseminated in the receiving environment through sewage.

METHOD

A total of 569 cefotaxime-resistant E. coli were isolated over a period of 19 months (August 2020 to February 2022) using ECC CHROMagar™ plates from 82 samples, antibiotic sensitivity profiles were determined, using Sensititre™ plates. The draft genome sequences were determined, using Illumina MiSeq-based sequencing. Complete genome sequences were determined, using Oxford Nanopore-based sequencing.

FINDINGS

All 569 strains obtained from influent (N=461) and effluent (N=108) were multi-drug resistant. Most of the sequenced strains (52 of 61) carried blaCTX-M-15 (38.5%) and blaCTX-M-27 (34.6%). The most prevalent sequence types (STs) for ESBL-carrying strains were ST131 (32.8%) and ST38 (21.3%). All CTX-M-27-carrying ST131 strains belonged to clade A or C1, while CTX-M-15-harbouring strains were present in all the clades. Five OXA-244-producing ST38 strains, genetically similar to epidemic-causing strains from Western Norway, France and the Netherlands, were isolated only from raw and treated sewage of the treatment plant receiving hospital sewage.

CONCLUSION

This is the first study showing persistent dissemination of OXA-244-producing ST38 clones through sewage in Norway, demonstrating that hospital sewage is the likely source of OXA-244-producing ST38 clones reaching the receiving environment.

The plasmidome associated with Gram-negative bloodstream infections: A large-scale observational study using complete plasmid assemblies

Plasmids carry genes conferring antimicrobial resistance and other clinically important traits, and contribute to the rapid dissemination of such genes. Previous studies using complete plasmid assemblies, which are essential for reliable inference, have been small and/or limited to plasmids carrying antimicrobial resistance genes (ARGs). In this study, we sequenced 1,880 complete plasmids from 738 isolates from bloodstream infections in Oxfordshire, UK. The bacteria had been originally isolated in 2009 (194 isolates) and 2018 (368 isolates), plus a stratified selection from intervening years (176 isolates). We demonstrate that plasmids are largely, but not entirely, constrained to a single host species, although there is substantial overlap between species of plasmid gene-repertoire. Most ARGs are carried by a relatively small number of plasmid groups with biological features that are predictable. Plasmids carrying ARGs (including those encoding carbapenemases) share a putative 'backbone' of core genes with those carrying no such genes. These findings suggest that future surveillance should, in addition to tracking plasmids currently associated with clinically important genes, focus on identifying and monitoring the dissemination of high-risk plasmid groups with the potential to rapidly acquire and disseminate these genes.

A convergent evolutionary pathway attenuating cellulose production drives enhanced virulence of some bacteria

Bacteria adapt to selective pressure in their immediate environment in multiple ways. One mechanism involves the acquisition of independent mutations that disable or modify a key pathway, providing a signature of adaptation via convergent evolution. Extra-intestinal pathogenic Escherichia coli (ExPEC) belonging to sequence type 95 (ST95) represent a global clone frequently associated with severe human infections including acute pyelonephritis, sepsis, and neonatal meningitis. Here, we analysed a publicly available dataset of 613 ST95 genomes and identified a series of loss-of-function mutations that disrupt cellulose production or its modification in 55.3% of strains. We show the inability to produce cellulose significantly enhances ST95 invasive infection in a rat model of neonatal meningitis, leading to the disruption of intestinal barrier integrity in newborn pups and enhanced dissemination to the liver, spleen and brain. Consistent with these observations, disruption of cellulose production in ST95 augmented innate immune signalling and tissue neutrophil infiltration in a mouse model of urinary tract infection. Mutations that disrupt cellulose production were also identified in other virulent ExPEC STs, Shigella and Salmonella, suggesting a correlative association with many Enterobacteriaceae that cause severe human infection. Together, our findings provide an explanation for the emergence of hypervirulent Enterobacteriaceae clones.

Genomic epidemiology reveals geographical clustering of multidrug-resistant Escherichia coli ST131 associated with bacteraemia in Wales

Antibiotic resistance is a significant global public health concern. Uropathogenic Escherichia coli sequence type (ST)131, a widely prevalent multidrug-resistant clone, is frequently associated with bacteraemia. This study investigates third-generation cephalosporin resistance in bloodstream infections caused by E. coli ST131. From 2013-2014 blood culture surveillance in Wales, 142 E. coli ST131 genomes were studied alongside global data. All three major ST131 clades were represented across Wales, with clade C/H30 predominant (n = 102/142, 71.8%). Consistent with global findings, Welsh strains of clade C/H30 contain β-lactamase genes from the blaCTX-M-1 group (n = 65/102, 63.7%), which confer resistance to third-generation cephalosporins. Most Welsh clade C/H30 genomes belonged to sub-clade C2/H30Rx (58.3%). A Wales-specific sub-lineage, named GB-WLS.C2, diverged around 1996-2000. An introduction to North Wales around 2002 led to a localised cluster by 2009, depicting limited genomic diversity within North Wales. This investigation emphasises the value of genomic epidemiology, allowing the detection of genetically similar strains in local areas, enabling targeted and timely public health interventions.

Modulation of multidrug-resistant clone success in Escherichia coli populations: a longitudinal, multi-country, genomic and antibiotic usage cohort study

BACKGROUND

The effect of antibiotic usage on the success of multidrug-resistant (MDR) clones in a population remains unclear. With this genomics-based molecular epidemiology study, we aimed to investigate the contribution of antibiotic use to Escherichia coli clone success, relative to intra-strain competition for colonisation and infection.

METHODS

We sequenced all the available E coli bloodstream infection isolates provided by the British Society for Antimicrobial Chemotherapy (BSAC) from 2012 to 2017 (n=718) and combined these with published data from the UK (2001-11; n=1090) and Norway (2002-17; n=3254). Defined daily dose (DDD) data from the European Centre for Disease Prevention and Control (retrieved on Sept 21, 2021) for major antibiotic classes (β-lactam, tetracycline, macrolide, sulfonamide, quinolone, and non-penicillin β-lactam) were used together with sequence typing, resistance profiling, regression analysis, and non-neutral Wright-Fisher simulation-based modelling to enable systematic comparison of resistance levels, clone success, and antibiotic usage between the UK and Norway.

FINDINGS

Sequence type (ST)73, ST131, ST95, and ST69 accounted for 892 (49·3%) of 1808 isolates in the BSAC collection. In the UK, the proportion of ST69 increased between 2001-10 and 2011-17 (p=0·0004), whereas the proportions of ST73 and ST95 did not vary between periods. ST131 expanded quickly after its emergence in 2003 and its prevalence remained consistent throughout the study period (apart from a brief decrease in 2009-10). The extended-spectrum β-lactamase (ESBL)-carrying, globally disseminated MDR clone ST131-C2 showed overall greater success in the UK (154 [56·8%] of 271 isolates in 2003-17) compared with Norway (51 [18·3%] of 278 isolates in 2002-17; p<0·0001). DDD data indicated higher total use of antimicrobials in the UK, driven mainly by the class of non-penicillin β-lactams, which were used between 2·7-times and 5·1-times more in the UK per annum (ratio mean 3·7 [SD 0·8]). This difference was associated with the higher success of the MDR clone ST131-C2 (pseudo-R2 69·1%). A non-neutral Wright-Fisher model replicated the observed expansion of non-MDR and MDR sequence types under higher DDD regimes.

INTERPRETATION

Our study indicates that resistance profiles of contemporaneously successful clones can vary substantially, warranting caution in the interpretation of correlations between aggregate measures of resistance and antibiotic usage. Our study further suggests that in countries with low-to-moderate use of antibiotics, such as the UK and Norway, the extent of non-penicillin β-lactam use modulates rather than determines the success of widely disseminated MDR ESBL-carrying E coli clones. Detailed understanding of underlying causal drivers of success is important for improved control of resistant pathogens.

FUNDING

Trond Mohn Foundation, Marie Skłodowska-Curie Actions, European Research Council, Royal Society, and Wellcome Trust.

Genetic features of antimicrobial drug-susceptible extraintestinal pathogenic Escherichia coli pandemic sequence type 95

IMPORTANCE

Despite the increasing prevalence of antibiotic-resistant Escherichia coli strains that cause urinary tract and bloodstream infections, a major pandemic lineage of extraintestinal pathogenic E. coli (ExPEC) ST95 has a comparatively low frequency of drug resistance. We compared the genomes of 1,749 ST95 isolates to identify genetic features that may explain why most strains of ST95 resist becoming drug-resistant. Identification of such genomic features could contribute to the development of novel strategies to prevent the spread of antibiotic-resistant genes and devise new measures to control antibiotic-resistant infections.

Evolution of extended-spectrum β-lactamase-producing ST131 Escherichia coli at a single hospital over 15 years

Escherichia coli belonging to sequence type ST131 constitute a globally distributed pandemic lineage that causes multidrug-resistant extra-intestinal infections. ST131 E. coli frequently produce extended-spectrum β-lactamases (ESBLs), which confer resistance to many β-lactam antibiotics and make infections difficult to treat. We sequenced the genomes of 154 ESBL-producing E. coli clinical isolates belonging to the ST131 lineage from patients at the University of Pittsburgh Medical Center (UPMC) between 2004 and 2018. Isolates belonged to the well described ST131 clades A (8%), B (3%), C1 (33%), and C2 (54%). An additional four isolates belonged to another distinct subclade within clade C and encoded genomic characteristics that have not been previously described. Time-dated phylogenetic analysis estimated that the most recent common ancestor (MRCA) for all clade C isolates from UPMC emerged around 1989, consistent with previous studies. We identified multiple genes potentially under selection in clade C, including the cell wall assembly gene ftsI, the LPS biosynthesis gene arnC, and the yersiniabactin uptake receptor fyuA. Diverse ESBL genes belonging to the blaCTX-M, blaSHV, and blaTEM families were identified; these genes were found at varying numbers of loci and in variable numbers of copies across isolates. Analysis of ESBL flanking regions revealed diverse mobile elements that varied by ESBL type. Overall, our findings show that ST131 subclades C1 and C2 dominated and were stably maintained among patients in the same hospital and uncover possible signals of ongoing adaptation within the clade C ST131 lineage.

ProQ binding to small RNA RyfA promotes virulence and biofilm formation in avian pathogenic Escherichia coli

Avian pathogenic Escherichia coli (APEC) is a notable subpathotype of the nonhuman extraintestinal pathogenic E. coli (ExPEC). Recognized as an extraintestinal foodborne pathogen, the zoonotic potential of APEC/ExPEC allows for cross-host transmission via APEC-contaminated poultry meat and eggs. ProQ, an RNA binding protein, is evolutionarily conserved in E. coli. However, its regulatory roles in the biofilm formation and virulence of APEC/ExPEC have not been explored. In this study, proQ deletion in the APEC strain FY26 significantly compromised its biofilm-forming ability. Furthermore, animal tests and cellular infection experiments showed that ProQ depletion significantly attenuated APEC virulence, thereby diminishing its capacity for bloodstream infection and effective adherence to and persistence within host cells. Transcriptome analysis revealed a decrease in the transcription level of the small RNA (sRNA) RyfA in the mutant FY26ΔproQ, suggesting a direct interaction between the sRNA RyfA and ProQ. This interaction might indicate that sRNA RyfA is a novel ProQ-associated sRNA. Moreover, the direct binding of ProQ to the sRNA RyfA was crucial for APEC biofilm formation, pathogenicity, adhesion, and intracellular survival. In conclusion, our findings provide detailed insight into the interaction between ProQ and sRNA RyfA and deepen our understanding of the regulatory elements that dictate APEC virulence and biofilm development. Such insights are instrumental in developing strategies to counteract APEC colonization within hosts and impede APEC biofilm establishment on food surfaces.

Dissemination of ESBL-producing E. coli ST131 through wastewater and environmental water in Switzerland

The E. coli lineage ST131 is a major cause of multidrug-resistant urinary tract and bloodstream infections worldwide. Recently emerged ST131 sublineages spread globally within few years, but their dissemination routes are incompletely understood. In this study, we investigate the potential role of wastewater and surface water in the spread of extended-spectrum β-lactamase (ESBL)-producing ST131. Streams, lakes, and two wastewater treatment plants (WWTPs) in the canton of Zug, Switzerland, were consecutively sampled over 1.5 years. ST131 was detected in 38% of the samples taken downstream (1-5 km) of WWTP discharge sites, but usually absent in water bodies distant from urban areas or WWTP discharge. Specific strains were repeatedly isolated (≤5 pairwise cgSNP distance) from wastewater or river sites downstream of effluent discharge, indicating their repeated entry or persistence in WWTPs in large concentrations. Genetic characterization of the ESBL-producing water isolates revealed a predominance of clades A and C1 and an emerging ciprofloxacin-resistant sublineage with mutations in quinolone resistance determining regions (QRDR) within clade A. Multiple isolates belonged to internationally circulating sublineages, including C1-M27 and papGII + sublineages with chromosomally encoded ESBLs. This study demonstrates that the clinically relevant E. coli lineage ST131 pollutes river ecosystems, representing a significant challenge to public health and to technologies to minimize their entry into the water environment.

Parallel loss of type VI secretion systems in two multi-drug-resistant Escherichia coli lineages

The repeated emergence of multi-drug-resistant (MDR) Escherichia coli clones is a threat to public health globally. In recent work, drug-resistant E. coli were shown to be capable of displacing commensal E. coli in the human gut. Given the rapid colonization observed in travel studies, it is possible that the presence of a type VI secretion system (T6SS) may be responsible for the rapid competitive advantage of drug-resistant E. coli clones. We employed large-scale genomic approaches to investigate this hypothesis. First, we searched for T6SS genes across a curated dataset of over 20 000 genomes representing the full phylogenetic diversity of E. coli. This revealed large, non-phylogenetic variation in the presence of T6SS genes. No association was found between T6SS gene carriage and MDR lineages. However, multiple clades containing MDR clones have lost essential structural T6SS genes. We characterized the T6SS loci of ST410 and ST131 and identified specific recombination and insertion events responsible for the parallel loss of essential T6SS genes in two MDR clones.

Occurrence of Major Human Extraintestinal Pathogenic Escherichia coli Sequence Types Among Diarrheic Pet Animals: A Potential Public Health Threat

Background: Extraintestinal pathogenic Escherichia coli (ExPEC) has become a mounting public health concern. The present study was conducted to address the role of diarrheic pet animals as potential reservoirs for major human ExPEC sequence types (STs). Materials and Methods: Rectal swabs were collected from 145 diarrheic pet animals (75 dogs and 70 cats). Samples were processed for isolation and identification of E. coli by culture methods. Afterward, ExPEC isolates were identified on a molecular basis through detection of ExPEC phylogroups (B2 and D) coupled with carriage of two or more of the virulence genes associated with ExPEC (papAH, papC, sfa/focDE, afa/draBC, iutA, and kpsMT II). ExPEC STs 131, 73, 69, and 95 were identified among ExPEC isolates by quadruplex PCR and tested for their antimicrobial susceptibility. Eventually, two isolates underwent gene sequencing for the phylogenetic analysis. Results: Of 145 pet animals, 16 (11%) E. coli strains were identified as ExPEC, in which 15 (10.3%) isolates belonged to phylogroup B2 and 1 (0.69%) strain belonged to phylogroup D. The major human ExPEC STs were detected in 13 (9%) isolates, whereas the prevalence rates were 5.3% and 12.9% for dogs and cats, respectively. The isolation rates of ExPEC STs were 4.8%, 2.8%, 0.69%, and 0.69% for ST73, ST131, ST95, and ST69, respectively. Regarding the prevalence of virulence genes among ExPEC STs, the most prevalent ones were papC and sfa/focDE (92.3%), followed by papAH (76.9%), iutA (53.8%), afa/draBC (30.8%), and kpsMT II (30.8%). Moreover, 38.5% of the obtained human ExPEC STs were multidrug resistant. The phylogenetic analysis of two ExPEC ST73 gene sequences showed high genetic relatedness to those isolated from humans in different countries. Conclusions: The fecal carriage of major human ExPEC STs among diarrheic dogs and cats poses a potential zoonotic hazard with serious public health implications.

Phylogenomic analysis of a global collection of Escherichia coli ST38: evidence of interspecies and environmental transmission?

IMPORTANCE

Extraintestinal pathogenic Escherichia coli (ExPEC) sequence type (ST) 38 is one of the top 10 human pandemic lineages. Although a major cause of urinary tract and blood stream infections, ST38 has been poorly characterized from a global phylogenomic perspective. A comprehensive genome-scale analysis of 925 ST38 isolate genomes identified two broad ancestral clades and linkage of discrete ST38 clusters with specific bla CTX-M variants. In addition, the clades and clusters carry important virulence genes, with diverse but poorly characterized plasmids. Numerous putative interhost and environment transmission events were identified here by the presence of ST38 clones (defined as isolates with ≤35 SNPs) within humans, companion animals, food sources, urban birds, wildlife, and the environment. A small cluster of international ST38 clones from diverse sources, likely representing progenitors of a hospital outbreak that occurred in Brisbane, Australia, in 2017, was also identified. Our study emphasizes the importance of characterizing isolate genomes derived from nonhuman sources and geographical locations, without any selection bias.

Identification and Evaluation of Pathogenic Genes (traT, hly, aer, pap, and fimH) and Antibiotic Resistance Genes (blaTEM, blaSHV, and blaCTX) in Escherichia coli in Patients Referred to Gonabad Hospitals, Iran

Background

Urinary tract infection (UTI) is one of the common bacterial infections. Escherichia coli is the most common cause of UTI. In this research, the prevalence of several virulence factors and beta-lactam resistance genes was investigated.

Methods

One hundred E. coli isolates were collected from patients' specimens with UTI referred to Allame-Bohlol Gonabadi hospital. Polymerase chain reaction (PCR) was performed to identify five pathogenic genes (fimH, aer, pap, hly, traT) and three antibiotic resistance genes (blaTEM, blaCTX, blaSHV).

Results

The frequencies of blaSHV, blaTEM and blaCTX beta-lactamase genes among extended-spectrum-beta-lactamases (ESBLs) positive isolates were 11.1%, 48.1%, and 93.3%, respectively. A significant number of isolates were resistant to the most commonly used antibiotics.

Conclusion

Pathogenic genes may also increase the severity, progression, and expansion of urinary tract infections. Therefore, identifying these genes as critical controllers of illness can use for better manage the treatment.

The clinical, genomic, and microbiological profile of invasive multi-drug resistant Escherichia coli in a major teaching hospital in the United Kingdom

Escherichia coli is a ubiquitous component of the human gut microbiome, but is also a common pathogen, causing around 40, 000 bloodstream infections (BSI) in the United Kingdom (UK) annually. The number of E. coli BSI has increased over the last decade in the UK, and emerging antimicrobial resistance (AMR) profiles threaten treatment options. Here, we combined clinical, epidemiological, and whole genome sequencing data with high content imaging to characterise over 300 E. coli isolates associated with BSI in a large teaching hospital in the East of England. Overall, only a limited number of sequence types (ST) were responsible for the majority of organisms causing invasive disease. The most abundant (20 % of all isolates) was ST131, of which around 90 % comprised the pandemic O25b:H4 group. ST131-O25b:H4 isolates were frequently multi-drug resistant (MDR), with a high prevalence of extended spectrum β-lactamases (ESBL) and fluoroquinolone resistance. There was no association between AMR phenotypes and the source of E. coli bacteraemia or whether the infection was healthcare-associated. Several clusters of ST131 were genetically similar, potentially suggesting a shared transmission network. However, there was no clear epidemiological associations between these cases, and they included organisms from both healthcare-associated and non-healthcare-associated origins. The majority of ST131 isolates exhibited strong binding with an anti-O25b antibody, raising the possibility of developing rapid diagnostics targeting this pathogen. In summary, our data suggest that a restricted set of MDR E. coli populations can be maintained and spread across both community and healthcare settings in this location, contributing disproportionately to invasive disease and AMR.

Accurate and fast graph-based pangenome annotation and clustering with ggCaller

Bacterial genomes differ in both gene content and sequence mutations, which underlie extensive phenotypic diversity, including variation in susceptibility to antimicrobials or vaccine-induced immunity. To identify and quantify important variants, all genes within a population must be predicted, functionally annotated, and clustered, representing the "pangenome." Despite the volume of genome data available, gene prediction and annotation are currently conducted in isolation on individual genomes, which is computationally inefficient and frequently inconsistent across genomes. Here, we introduce the open-source software graph-gene-caller (ggCaller). ggCaller combines gene prediction, functional annotation, and clustering into a single workflow using population-wide de Bruijn graphs, removing redundancy in gene annotation and resulting in more accurate gene predictions and orthologue clustering. We applied ggCaller to simulated and real-world bacterial data sets containing hundreds or thousands of genomes, comparing it to current state-of-the-art tools. ggCaller has considerable speed-ups with equivalent or greater accuracy, particularly with data sets containing complex sources of error, such as assembly contamination or fragmentation. ggCaller is also an important extension to bacterial genome-wide association studies, enabling querying of annotated graphs for functional analyses. We highlight this application by functionally annotating DNA sequences with significant associations to tetracycline and macrolide resistance in Streptococcus pneumoniae, identifying key resistance determinants that were missed when using only a single reference genome. ggCaller is a novel bacterial genome analysis tool with applications in bacterial evolution and epidemiology.

Distribution of virulence-associated genes, antibiotic resistance and phylogenetic groups in Escherichia coli isolated from domestic and racing pigeons

Despite a lot of information about virulence and resistance of Escherichia coli (E. coli) in poultry, very limited data are currently available on its occurrence in pigeon isolates, although this poses a threat to human and animal health. Therefore, this study was conducted to explore the phylogenetic classification, antibiotic sensitivity, and virulence factors in E. coli isolated from cloacal swabs of domestic pigeons bred for meat (n = 47) and racing pigeons (n = 44). The most frequent phylogroup in racing pigeons was E (36, 82.00%), unlike domestic pigeons (B2- 19, 40.00%). The most abundant iron uptake system in both groups of bird was feoB (racing = 40, 90.90%; domestic = 44, 93.61%). The presence of ibeA (52, 57.10%) and kpsMTII (46, 50.50%) genes was detected in more than half of all strains belonging exclusively to phylogroups B2, D, E, F, clade I. Antibiotic resistance was higher in racing pigeons. All racing pigeon isolates were resistant to tetracycline and trimethoprim + sulphonamide. Resistance to ciprofloxacin was determined in three isolates (6.38%) of domestic and 33 isolates (75%) of racing pigeons. Aminoglycosides and β-lactamases resistance were also recorded. One of the important detected phenotypic mechanisms of resistance occurring in isolates from racing pigeons was AGL AAC(6´)I. Our study confirms that healthy pigeons are a reservoir of antibiotic-resistant E. coli containing an arsenal of virulence factors, thus capable of potentially causing infection. Pigeons with the option to fly to multiple places can transfer virulent and resistant bacteria. Direct contact with pigeons and their faeces and the contamination of water and food pose a threat of infection to humans and other animal species.

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.

coli ESBL Study Group, Sundsfjord A, Gravningen K, Samuelsen Ø. Community carriage of ESBL-producing Escherichia coli and Klebsiella pneumoniae: a cross-sectional study of risk factors and comparative genomics of carriage and clinical isolates

The global prevalence of infections caused by extended-spectrum β-lactamase-producing Enterobacterales (ESBL-E) is increasing, and for Escherichia coli, observations indicate that this is partly driven by community-onset cases. The ESBL-E population structure in the community is scarcely described, and data on risk factors for carriage are conflicting. Here, we report the prevalence and population structure of fecal ESBL-producing E. coli and Klebsiella pneumoniae (ESBL-Ec/Kp) in a general adult population, examine risk factors, and compare carriage isolates with contemporary clinical isolates. Fecal samples obtained from 4,999 participants (54% women) ≥40 years in the seventh survey of the population-based Tromsø Study, Norway (2015, 2016), were screened for ESBL-Ec/Kp. In addition, we included 118 ESBL-Ec clinical isolates from the Norwegian surveillance program in 2014. All isolates were whole-genome sequenced. Risk factors associated with carriage were analyzed using multivariable logistic regression. ESBL-Ec gastrointestinal carriage prevalence was 3.3% [95% confidence interval (CI) 2.8%-3.9%, no sex difference] and 0.08% (0.02%-0.20%) for ESBL-Kp. For ESBL-Ec, travel to Asia was the only independent risk factor (adjusted odds ratio 3.46, 95% CI 2.18-5.49). E. coli ST131 was most prevalent in both collections. However, the ST131 proportion was significantly lower in carriage (24%) versus clinical isolates (58%, P < 0.001). Carriage isolates were genetically more diverse with a higher proportion of phylogroup A (26%) than clinical isolates (5%, P < 0.001), indicating that ESBL gene acquisition occurs in a variety of E. coli lineages colonizing the gut. STs commonly related to extraintestinal infections were more frequent in clinical isolates also carrying a higher prevalence of antimicrobial resistance, which could indicate clone-associated pathogenicity.IMPORTANCEESBL-Ec and ESBL-Kp are major pathogens in the global burden of antimicrobial resistance. However, there is a gap in knowledge concerning the bacterial population structure of human ESBL-Ec/Kp carriage isolates in the community. We have examined ESBL-Ec/Kp isolates from a population-based study and compared these to contemporary clinical isolates. The large genetic diversity of carriage isolates indicates frequent ESBL gene acquisition, while those causing invasive infections are more clone dependent and associated with a higher prevalence of antibiotic resistance. The knowledge of factors associated with ESBL carriage helps to identify patients at risk to combat the spread of resistant bacteria within the healthcare system. Particularly, previous travel to Asia stands out as a major risk factor for carriage and should be considered in selecting empirical antibiotic treatment in critically ill patients.

ESBL-Producing and Non-ESBL-Producing Escherichia coli Isolates from Urinary Tract Differ in Clonal Distribution, Virulence Gene Content and Phylogenetic Group

Purpose

The objectives of this study are to determine the differences in clonality, virulence gene (VG) content and phylogenetic group between non extended-spectrum beta-lactamase-producing E. coli (non-ESBL-EC) and ESBL-EC isolates from urine.

Patients and Methods

This study characterized a total of 100 clinical E. coli isolates consecutively obtained from the inpatients hospitalized in The First Affiliated Hospital of Ningbo University in China by polymerase-chain reaction (PCR).

Results

Phylogenetic group B2 was found to be the most prevalent in both ESBL-EC and non-ESBL-EC group. Among 100 clinical isolates, the count of acquired virulence genes in group B2 was found to be significantly higher than that in group A, B1, and D (p <0.001). Additionally, the presence of content within virulence genes (the total number of virulence genes detected per isolate) in B2 of non-ESBL-EC and ESBL-EC showed a significant difference (p<0.001). ST131 was detected exclusively in ESBL-EC, while ST95 and ST73 were the main sequence types in non-ESBL-EC.

Conclusion

Our study demonstrated the different distribution of MLST, phylogenetic group in ESBL-EC and non-ESBL-EC group. The inverse association between beta-lactamase resistance and VG content performed in this study should get a lot more attention. At the same time, we should also be wary of the appearance of non-ESBL-EC isolates of group B2 harboring more virulence genes which will lead to high pathogenicity.

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.

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.

Genomic analysis of extended-spectrum beta-lactamase (ESBL) producing Escherichia coli colonising adults in Blantyre, Malawi reveals previously undescribed diversity

Escherichia coli is one of the most prevalent Gram-negative species associated with drug resistant infections. Strains that produce extended-spectrum beta-lactamases (ESBLs) or carbapenemases are both particularly problematic and disproportionately impact resource limited healthcare settings where last-line antimicrobials may not be available. A large number of E. coli genomes are now available and have allowed insights into pathogenesis and epidemiology of ESBL E. coli but genomes from sub-Saharan Africa (sSA) are significantly underrepresented. To reduce this gap, we investigated ESBL-producing E. coli colonising adults in Blantyre, Malawi to assess bacterial diversity and AMR determinants and to place these isolates in the context of the wider population structure. We performed short-read whole-genome sequencing of 473 colonising ESBL E. coli isolated from human stool and contextualised the genomes with a previously curated multi-country collection of 10 146 E. coli genomes and sequence type (ST)-specific collections for our three most commonly identified STs. These were the globally successful ST131, ST410 and ST167, and the dominant ESBL genes were bla CTX-M, mirroring global trends. However, 37 % of Malawian isolates did not cluster with any isolates in the curated multicountry collection and phylogenies were consistent with locally spreading monophyletic clades, including within the globally distributed, carbapenemase-associated B4/H24RxC ST410 lineage. A single ST2083 isolate in this collection harboured a carbapenemase gene. Long read sequencing demonstrated the presence of a globally distributed ST410-associated carbapenemase carrying plasmid in this isolate, which was absent from the ST410 strains in our collection. We conclude there is a risk that carbapenem resistance in E. coli could proliferate rapidly in Malawi under increasing selection pressure, and that both ongoing antimicrobial stewardship and genomic surveillance are critical as local carbapenem use increases.

The evolution of colistin resistance increases bacterial resistance to host antimicrobial peptides and virulence

Antimicrobial peptides (AMPs) offer a promising solution to the antibiotic resistance crisis. However, an unresolved serious concern is that the evolution of resistance to therapeutic AMPs may generate cross-resistance to host AMPs, compromising a cornerstone of the innate immune response. We systematically tested this hypothesis using globally disseminated mobile colistin resistance (MCR) that has been selected by the use of colistin in agriculture and medicine. Here, we show that MCR provides a selective advantage to Escherichia coli in the presence of key AMPs from humans and agricultural animals by increasing AMP resistance. Moreover, MCR promotes bacterial growth in human serum and increases virulence in a Galleria mellonella infection model. Our study shows how the anthropogenic use of AMPs can drive the accidental evolution of resistance to the innate immune system of humans and animals. These findings have major implications for the design and use of therapeutic AMPs and suggest that MCR may be difficult to eradicate, even if colistin use is withdrawn.

A Cross-Sectional Study of Companion Animal-Derived Multidrug-Resistant Escherichia coli in Hangzhou, China

Antimicrobial resistance poses a challenge to global public health, and companion animals could serve as the reservoir for antimicrobial-resistant bacteria. However, the prevalence of antimicrobial-resistant bacteria, especially multidrug-resistant (MDR) bacteria, and the associated risk factors from companion animals are partially understood. Here, we aim to investigate the prevalence of MDR Escherichia coli, as an indicator bacterium, in pet cats and dogs in Hangzhou, China, and evaluate the factors affecting the prevalence of MDR E. coli. The proportion of pets carrying MDR E. coli was 35.77% (49/137), i.e., 40.96% (34/83) for dogs and 27.28% (15/54) for cats. Isolates resistant to trimethoprim-sulfamethoxazole (49.40% and 44.44%), amoxicillin-clavulanic acid (42.17% and 38.89%), and nalidixic acid (40.96% and 35.19%) were the most prevalent in dogs and cats. Interestingly, comparable prevalence of MDR E. coli was observed in pet dogs and cats regardless of the health condition and the history of antibiotic use. Genetic diversity analysis indicates a total of 86 sequencing types (23 clonal complexes), with ST12 being the most dominant. Further genomic investigation of a carbapenem-resistant E. coli ST410 isolate reveals abundant antimicrobial-resistance genes and a plasmid-borne carbapenemase gene (NDM-5) flanked by insertion sequences of IS91 and IS31, suggesting the plasmid and insertion sequences may be involved in carbapenem-resistance dissemination. These data show that companion animal-derived MDR bacteria could threaten public health, and further regulation and supervision of antimicrobial use in pet clinics should be established in China. IMPORTANCE MDR Escherichia coli are considered a global threat because of the decreasing options for antimicrobial therapy. Companion animals could be a reservoir of MDR E. coli, and the numbers of pets and households owning pets in China are booming. However, the prevalence and risk factors of MDR E. coli carriage in Chinese pets were rarely studied. Here, we investigated the prevalence of MDR E. coli in pets in Hangzhou, one of the leading cities with the most established pet market in China, and explored the factors that affected the prevalence. Our findings showed high prevalences of MDR E. coli in pet dogs and cats regardless of their health condition and the history of antibiotic use, suggesting their potential role of public health risk. A call-to-action for improved regulation of antimicrobial use in companion animal is needed in China.

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.