Fine-Scale Structure Analysis Shows Epidemic Patterns of Clonal Complex 95, a Cosmopolitan Escherichia coli Lineage Responsible for Extraintestinal Infection

One day in Denmark: whole-genome sequence-based analysis of Escherichia coli isolates from clinical settings

BACKGROUND

WGS can potentially be routinely used in clinical microbiology settings, especially with the increase in sequencing accuracy and decrease in cost. Escherichia coli is the most common bacterial species analysed in those settings, thus fast and accurate diagnostics can lead to reductions in morbidity, mortality and healthcare costs.

OBJECTIVES

To evaluate WGS for diagnostics and surveillance in a collection of clinical E. coli; to examine the pool of antimicrobial resistance (AMR) determinants circulating in Denmark and the most frequent STs; and to evaluate core-genome MLST (cgMLST) and SNP-based clustering approaches for detecting genetically related isolates.

METHODS

We analysed the genomes of 699 E. coli isolates collected throughout all Danish Clinical Microbiology Laboratories. We used rMLST and KmerFinder for species identification, ResFinder for prediction of AMR, and PlasmidFinder for plasmid identification. We used Center for Genomic Epidemiology MLST, cgMLSTFinder and CSI Phylogeny to perform typing and clustering analysis.

RESULTS

Genetic AMR determinants were detected in 56.2% of isolates. We identified 182 MLSTs, most frequently ST-69, ST-73, ST-95 and ST-131. Using a maximum 15-allele difference as the threshold for genetic relatedness, we identified 23 clusters. SNP-based phylogenetic analysis within clusters revealed from 0 to 13 SNPs, except two cases with 111 and 461 SNPs.

CONCLUSIONS

WGS data are useful to characterize clinical E. coli isolates, including predicting AMR profiles and subtyping in concordance with surveillance data. We have shown that it is possible to adequately cluster isolates through a cgMLST approach, but it remains necessary to define proper interpretative criteria.

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.

HlyF, an underestimated virulence factor of uropathogenic Escherichia coli

OBJECTIVES

Urinary tract infections (UTIs) are primarily caused by uropathogenic Escherichia coli (UPEC). This study aims to elucidate the role of the virulence factor HlyF in the epidemiology and pathophysiology of UTIs and investigate the dissemination of plasmids carrying the hlyF gene.

METHODS

An epidemiological analysis was conducted on a representative collection of 225 UPEC strains isolated from community-acquired infections. Selected hlyF+ strains were fully sequenced using a combination of Illumina and Nanopore technologies. To investigate the impact of HlyF, a murine model of UTI was utilized to compare clinical signs, bacterial loads in the bladder, kidney, and spleen, onset of bacteraemia, and inflammation through cytokine quantification among wild-type hlyF+ strains, isogenic mutants, and complemented mutants.

RESULTS

Our findings demonstrate that 20% of UPEC encode the HlyF protein. These hlyF+ UPEC strains exhibited enhanced virulence, frequently leading to pyelonephritis accompanied by bloodstream infections. Unlike typical UPEC strains, hlyF+ UPEC strains demonstrate a broader phylogroup distribution and possess a unique array of virulence factors and antimicrobial resistance genes, primarily carried by ColV-like plasmids. In the murine UTI model, expression of HlyF was linked to the UPECs' capacity to induce urosepsis and elicit an exacerbated inflammatory response, setting them apart from typical UPEC strains.

DISCUSSION

Overall, our results strongly support the notion that HlyF serves as a significant virulence factor for UPECs, and the dissemination of ColV-like plasmids encoding HlyF warrants further investigation.

Genomics Analysis to Identify Multiple Genetic Determinants That Drive the Global Transmission of the Pandemic ST95 Lineage of Extraintestinal Pathogenic Escherichia coli (ExPEC)

Extraintestinal pathogenic Escherichia coli (ExPEC) is a pathogen that causes host extraintestinal diseases. The ST95 E. coli lineage is one of the dominant ExPEC lineages in humans and poultry. In this study, we took advantage of extensive E. coli genomes available through public open-access databases to construct a detailed understanding of the phylogeny and evolution of ST95. We used a high variability of accessory genomes to highlight the diversity and dynamic traits of ST95. Isolates from diverse hosts and geographic sources were randomly located on the phylogenetic tree, which suggested that there is no host specificity for ST95. The time-scaled phylogeny showed that ST95 is an ancient and long-lasting lineage. The virulence genes, resistance genes, and pathogenicity islands (PAIs) were characterized in ST95 pan-genomes to provide novel insights into the pathogenicity and multidrug resistance (MDR) genotypes. We found that a pool of large plasmids drives virulence and MDR. Based on the unique genes in the ST95 pan-genome, we designed a novel multiplex PCR reaction to rapidly detect ST95. Overall, our study addressed a gap in the current understanding of ST95 ExPEC genomes, with significant implications for recognizing the success and spread of ST95.

EnteroBase: hierarchical clustering of 100 000s of bacterial genomes into species/subspecies and populations

The definition of bacterial species is traditionally a taxonomic issue while bacterial populations are identified by population genetics. These assignments are species specific, and depend on the practitioner. Legacy multilocus sequence typing is commonly used to identify sequence types (STs) and clusters (ST Complexes). However, these approaches are not adequate for the millions of genomic sequences from bacterial pathogens that have been generated since 2012. EnteroBase (http://enterobase.warwick.ac.uk) automatically clusters core genome MLST allelic profiles into hierarchical clusters (HierCC) after assembling annotated draft genomes from short-read sequences. HierCC clusters span core sequence diversity from the species level down to individual transmission chains. Here we evaluate HierCC's ability to correctly assign 100 000s of genomes to the species/subspecies and population levels for Salmonella, Escherichia, Clostridoides, Yersinia, Vibrio and Streptococcus. HierCC assignments were more consistent with maximum-likelihood super-trees of core SNPs or presence/absence of accessory genes than classical taxonomic assignments or 95% ANI. However, neither HierCC nor ANI were uniformly consistent with classical taxonomy of Streptococcus. HierCC was also consistent with legacy eBGs/ST Complexes in Salmonella or Escherichia and with O serogroups in Salmonella. Thus, EnteroBase HierCC supports the automated identification of and assignment to species/subspecies and populations for multiple genera. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Molecular and metabolic characteristics of wastewater associated Escherichia coli strains

We previously characterized the genetic diversity of Escherichia coli strains isolated from septic tanks in the Canberra region, Australia. In this study, we used repetitive element palindromic (REP) PCR fingerprinting to identify dominant REP-types belonging to phylogroups A and B1 strains across septic tanks. Subsequently, 76 E. coli strains were selected for whole-genome sequencing and phenotype microarrays. Comparative genome analysis was performed to compare septic tank E. coli genomes with a collection of 433 E. coli isolates from different hosts and freshwater. Clonal complexes (CCs) 10 (n = 15) and 399 (n = 10) along with sequence type (ST) 401 (n = 9) were the common lineages in septic tanks. CC10 strains have been detected from animal hosts and freshwater, whereas CC399 and ST401 strains appeared to be associated with septic tanks as they were uncommon in isolates from other sources. Comparative genome analysis revealed that CC399 and ST401 were genetically distinct from other isolates and carried an abundance of niche-specific traits involved in environmental adaptation. These strains also showed distinct metabolic characteristics, such as the ability to utilize pectin, which may provide a fitness advantage under nutrient-limited conditions. The results of this study characterized the adaptive mechanisms allowing E. coli to persist in wastewater.

Characterization of Escherichia coli harboring colibactin genes (clb) isolated from beef production and processing systems

Certain strains of Escherichia coli possess and express the toxin colibactin (Clb) which induces host mutations identical to the signature mutations of colorectal cancer (CRC) that lead to tumorigenic lesions. Since cattle are a known reservoir of several Enterobacteriaceae including E. coli, this study screened for clb amongst E. coli isolated from colons of cattle-at-harvest (entering beef processing facility; n = 1430), across the beef processing continuum (feedlot to finished subprimal beef; n = 232), and in ground beef (n = 1074). Results demonstrated that clb⁺E. coli were present in cattle and beef. Prevalence of clb⁺E. coli from colonic contents of cattle and ground beef was 18.3% and 5.5%, respectively. clb⁺E. coli were found susceptible to commonly used meat processing interventions. Whole genome sequencing of 54 bovine and beef clb⁺ isolates showed clb occurred in diverse genetic backgrounds, most frequently in phylogroup B1 (70.4%), MLST 1079 (42.6%), and serogroup O49 (40.7%).

Prevalence, diversity and genetic structure of Escherichia coli isolates from septic tanks

The present study investigated the diversity and genetic structure of Escherichia coli isolates from 100 septic tanks in the Canberra region, Australia. The physicochemical characteristics of the septic tanks were determined to examine the extent to which environmental factors might influence E. coli prevalence, diversity and population structure. The results of this study indicated that the temperature of the septic tank could explain some of the variation observed in the number of E. coli isolates recovered per septic tank, whereas pH was an important driver of E. coli diversity. Conductivity, pH and household size had a significant impact on E. coli population structure, and household size significantly affected the probability of detecting human-associated E. coli lineages [sequence types (STs) 69, 73, 95 and 131] in septic tanks. Phylogroup A and B1 strains were not randomly distributed among septic tanks, and the strong negative association between them may indicate intraspecific competition. The findings of this study suggest that the combination of environmental factors and intraspecific interactions may influence the distribution and genetic structure of E. coli in the environment.

F Plasmid Lineages in Escherichia coli ST95: Implications for Host Range, Antibiotic Resistance, and Zoonoses

Escherichia coli sequence type 95 (ST95) is an extraintestinal pathogenic E. coli (ExPEC) renowned for its ability to cause significant morbidity and mortality in humans and poultry. A core genome analysis of 668 ST95 isolates generated 10 clades (A to J), 5 of which are reported here for the first time. F plasmid replicon sequence typing showed that almost a third (178/668 [27%]) of the collection carry pUTI89 (F29:B10) and were restricted to clade A and a sublineage of clade B. In contrast, almost half (328/668 [49%]) of the collection across multiple clades harbor ColV plasmids (multiple F types). Strikingly, ST95 lineages with pUTI89 were almost exclusively from humans, while ColV⁺ ST95 lineages were sourced from poultry and humans. Clade I was notable because it comprises temporally and geographically matched ColV⁺ isolates sourced from human and retail poultry meat, suggesting interspecies transmission via food. Clade F contained ST95 isolates of bovine origin, none of which carried ColV or pUTI89 plasmids. Remarkably, an analysis of a cohort of 34,176 E. coli isolates comprising 2,570 sequence types mirrored what was observed in ST95: (i) pUTI89 was overwhelmingly linked to E. coli sourced from humans but almost entirely absent from 13,027 E. coli isolates recovered from poultry, pigs, and cattle, and (ii) E. coli isolates harboring ColV plasmids were from multiple sources, including humans, poultry, and swine. Overall, our data suggest that F plasmids influence E. coli host range, clade structure, and zoonotic potential in ST95 and ExPEC more broadly.

IMPORTANCEE. coli ST95 is one of five dominant ExPEC lineages globally and noted for causing urinary tract and bloodstream infections and neonatal meningitis in humans and colibacillosis in poultry. Using high-resolution phylogenomics, we show that F replicon sequence type is linked to ST95 clade structure and zoonotic potential. Specifically, human centric ST95 clades overwhelmingly harbor F29:B10 (pUTI89) plasmids, while clades carrying both human- and poultry-sourced isolates are typically ColV⁺ with multiple replicon types. Importantly, several clades identified clonal ColV⁺ ST95 isolates from human and poultry sources, but clade I, which housed temporally and spatially matched isolates, provided the most robust evidence. Notably, patterns of association of F replicon types with E. coli host were mirrored within a diverse collection of 34,176 E. coli genomes. Our studies indicate that the role of food animals as a source of human ExPEC disease is complex and warrants further investigation.

Genomic and Phenotypic Analysis of Heat and Sanitizer Resistance in Escherichia coli from Beef in Relation to the Locus of Heat Resistance

The locus of heat resistance (LHR) can confer heat resistance to Escherichia coli to various extents. This study investigated the phylogenetic relationships and the genomic and phenotypic characteristics of E. coli with or without LHR recovered from beef by direct plating or from enrichment broth at 42°C. LHR-positive E. coli isolates (n = 24) were subjected to whole-genome sequencing by short and long reads. LHR-negative isolates (n = 18) from equivalent sources as LHR-positive isolates were short-read sequenced. All isolates were assessed for decimal reduction time at 60°C (D_60°C) and susceptibility to the sanitizers E-SAN and Perox-E. Selected isolates were evaluated for growth at 42°C. The LHR-positive and -negative isolates were well separated on the core genome tree, with 22/24 positive isolates clustering into three clades. Isolates within clade 1 and 2, despite their different D_60°C values, were clonal, as determined by subtyping (multilocus sequence typing [MLST], core genome MLST, and serotyping). Isolates within each clade are of one serotype. The LHR-negative isolates were genetically diverse. The LHR-positive isolates had a larger (P < 0.001) median genome size by 0.3 Mbp (5.0 versus 4.7 Mbp) and overrepresentation of genes related to plasmid maintenance, stress response, and cryptic prophages but underrepresentation of genes involved in epithelial attachment and virulence. All LHR-positive isolates harbored a chromosomal copy of LHR, and all clade 2 isolates had an additional partial copy of LHR on conjugative plasmids. The growth rates at 42°C were 0.71 ± 0.02 and 0.65 ± 0.02 log(OD) h^-1 for LHR-positive and -negative isolates, respectively. No meaningful difference in sanitizer susceptibility was noted between LHR-positive and -negative isolates. IMPORTANCE Resistant bacteria are serious food safety and public health concerns. Heat resistance conferred by the LHR varies largely among different strains of E. coli. The findings in this study show that genomic background and composition of LHR, in addition to the presence of LHR, play an important role in the degree of heat resistance in E. coli and that strains with certain genetic backgrounds are more likely to acquire and maintain the LHR. Also, caution should be exercised when recovering E. coli at elevated temperatures, as the presence of LHR may confer growth advantages to some strains. Interestingly, the LHR-harboring strains seem to have evolved further from their primary animal host to adapt to their secondary habitat, as reflected by fewer genes involved in virulence and epithelial attachment. The phylogenetic relationships among the isolates point toward multiple mechanisms for acquisition of LHR by E. coli, likely prior to its being deposited on meat.

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

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

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

Background

Escherichia coli is the leading cause of bloodstream infections, associated with a significant mortality. Recent genomic analyses revealed that few clonal lineages are involved in bloodstream infections and captured the emergence of some of them. However, data on within sequence type (ST) population genetic structure evolution are rare.

Methods

We compared whole genome sequences of 912 E. coli isolates responsible for bloodstream infections from two multicenter clinical trials that were conducted in the Paris area, France, 12 years apart, in teaching hospitals belonging to the same institution (“Assistance Publique-Hôpitaux de Paris”). We analyzed the strains at different levels of granularity, i.e., the phylogroup, the ST complex (STc), and the within STc clone taking into consideration the evolutionary history, the resistance, and virulence gene content as well as the antigenic diversity of the strains.

Results

We found a mix of stability and changes overtime, depending on the level of comparison. Overall, we observed an increase in antibiotic resistance associated to a restricted number of genetic determinants and in strain plasmidic content, whereas phylogroup distribution and virulence gene content remained constant. Focusing on STcs highlighted the pauci-clonality of the populations, with only 11 STcs responsible for more than 73% of the cases, dominated by five STcs (STc73, STc131, STc95, STc69, STc10). However, some STcs underwent dramatic variations, such as the global pandemic STc131, which replaced the previously predominant STc95. Moreover, within STc131, 95 and 69 genomic diversity analysis revealed a highly dynamic pattern, with reshuffling of the population linked to clonal replacement sometimes coupled with independent acquisitions of virulence factors such as the pap gene cluster bearing a papGII allele located on various pathogenicity islands. Additionally, STc10 exhibited huge antigenic diversity evidenced by numerous O:H serotype/fimH allele combinations, whichever the year of isolation.

Conclusions

Altogether, these data suggest that the bloodstream niche is occupied by a wide but specific phylogenetic diversity and that highly specialized extra-intestinal clones undergo frequent turnover at the within ST level. Additional worldwide epidemiological studies overtime are needed in different geographical and ecological contexts to assess how generalizable these data are.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00892-0.

Insights into the acquisition of the pks island and production of colibactin in the Escherichia coli population

The pks island codes for the enzymes necessary for synthesis of the genotoxin colibactin, which contributes to the virulence of Escherichia coli strains and is suspected of promoting colorectal cancer. From a collection of 785 human and bovine E. coli isolates, we identified 109 strains carrying a highly conserved pks island, mostly from phylogroup B2, but also from phylogroups A, B1 and D. Different scenarios of pks acquisition were deduced from whole genome sequence and phylogenetic analysis. In the main scenario, pks was introduced and stabilized into certain sequence types (STs) of the B2 phylogroup, such as ST73 and ST95, at the asnW tRNA locus located in the vicinity of the yersiniabactin-encoding High Pathogenicity Island (HPI). In a few B2 strains, pks inserted at the asnU or asnV tRNA loci close to the HPI and occasionally was located next to the remnant of an integrative and conjugative element. In a last scenario specific to B1/A strains, pks was acquired, independently of the HPI, at a non-tRNA locus. All the pks-positive strains except 18 produced colibactin. Sixteen strains contained mutations in clbB or clbD, or a fusion of clbJ and clbK and were no longer genotoxic but most of them still produced low amounts of potentially active metabolites associated with the pks island. One strain was fully metabolically inactive without pks alteration, but colibactin production was restored by overexpressing the ClbR regulator. In conclusion, the pks island is not restricted to human pathogenic B2 strains and is more widely distributed in the E. coli population, while preserving its functionality.

Evolutionary Dynamics Based on Comparative Genomics of Pathogenic Escherichia coli Lineages Harboring Polyketide Synthase (pks) Island

The genotoxin colibactin is a secondary metabolite produced by the polyketide synthase (pks) island harbored by extraintestinal pathogenic E. coli (ExPEC) and other members of the Enterobacteriaceae that has been increasingly reported to have critical implications in human health. The present study entails a high-throughput whole-genome comparison and phylogenetic analysis of such pathogenic E. coli isolates to gain insights into the patterns of distribution, horizontal transmission, and evolution of the island. For the current study, 23 pks-positive ExPEC genomes were newly sequenced, and their virulome and resistome profiles indicated a preponderance of virulence encoding genes and a reduced number of genes for antimicrobial resistance. In addition, 4,090 E. coli genomes from the public domain were also analyzed for large-scale screening for pks-positive genomes, out of which a total of 530 pks-positive genomes were studied to understand the subtype-based distribution pattern(s). The pks island showed a significant association with the B2 phylogroup (82.2%) and a high prevalence in sequence type 73 (ST73; n = 179) and ST95 (n = 110) and the O6:H1 (n = 110) serotype. Maximum-likelihood (ML) phylogeny of the core genome and intergenic regions (IGRs) of the ST95 model data set, which was selected because it had both pks-positive and pks-negative genomes, displayed clustering in relation to their carriage of the pks island. Prevalence patterns of genes encoding RM systems in the pks-positive and pks-negative genomes were also analyzed to determine their potential role in pks island acquisition and the maintenance capability of the genomes. Further, the maximum-likelihood phylogeny based on the core genome and pks island sequences from 247 genomes with an intact pks island demonstrated horizontal gene transfer of the island across sequence types and serotypes, with few exceptions. This study vitally contributes to understanding of the lineages and subtypes that have a higher propensity to harbor the pks island-encoded genotoxin with possible clinical implications.

Horizontally acquired papGII-containing pathogenicity islands underlie the emergence of invasive uropathogenic Escherichia coli lineages

Escherichia coli is the leading cause of urinary tract infection, one of the most common bacterial infections in humans. Despite this, a genomic perspective is lacking regarding the phylogenetic distribution of isolates associated with different clinical syndromes. Here, we present a large-scale phylogenomic analysis of a spatiotemporally and clinically diverse set of 907 E. coli isolates, including 722 uropathogenic E. coli (UPEC) isolates. A genome-wide association approach identifies the (P-fimbriae-encoding) papGII locus as the key feature distinguishing invasive UPEC, defined as isolates associated with severe UTI, i.e., kidney infection (pyelonephritis) or urinary-source bacteremia, from non-invasive UPEC, defined as isolates associated with asymptomatic bacteriuria or bladder infection (cystitis). Within the E. coli population, distinct invasive UPEC lineages emerged through repeated horizontal acquisition of diverse papGII-containing pathogenicity islands. Our findings elucidate the molecular determinants of severe UTI and have implications for the early detection of this pathogen.

The population genetics of pathogenic Escherichia coli

Escherichia coli is a commensal of the vertebrate gut that is increasingly involved in various intestinal and extra-intestinal infections as an opportunistic pathogen. Numerous pathotypes that represent groups of strains with specific pathogenic characteristics have been described based on heterogeneous and complex criteria. The democratization of whole-genome sequencing has led to an accumulation of genomic data that render possible a population phylogenomic approach to the emergence of virulence. Few lineages are responsible for the pathologies compared with the diversity of commensal strains. These lineages emerged multiple times during E. coli evolution, mainly by acquiring virulence genes located on mobile elements, but in a specific chromosomal phylogenetic background. This repeated emergence of stable and cosmopolitan lineages argues for an optimization of strain fitness through epistatic interactions between the virulence determinants and the remaining genome.

Phylogenetic background and habitat drive the genetic diversification of Escherichia coli

Escherichia coli is mostly a commensal of birds and mammals, including humans, where it can act as an opportunistic pathogen. It is also found in water and sediments. We investigated the phylogeny, genetic diversification, and habitat-association of 1,294 isolates representative of the phylogenetic diversity of more than 5,000 isolates from the Australian continent. Since many previous studies focused on clinical isolates, we investigated mostly other isolates originating from humans, poultry, wild animals and water. These strains represent the species genetic diversity and reveal widespread associations between phylogroups and isolation sources. The analysis of strains from the same sequence types revealed very rapid change of gene repertoires in the very early stages of divergence, driven by the acquisition of many different types of mobile genetic elements. These elements also lead to rapid variations in genome size, even if few of their genes rise to high frequency in the species. Variations in genome size are associated with phylogroup and isolation sources, but the latter determine the number of MGEs, a marker of recent transfer, suggesting that gene flow reinforces the association of certain genetic backgrounds with specific habitats. After a while, the divergence of gene repertoires becomes linear with phylogenetic distance, presumably reflecting the continuous turnover of mobile element and the occasional acquisition of adaptive genes. Surprisingly, the phylogroups with smallest genomes have the highest rates of gene repertoire diversification and fewer but more diverse mobile genetic elements. This suggests that smaller genomes are associated with higher, not lower, turnover of genetic information. Many of these genomes are from freshwater isolates and have peculiar traits, including a specific capsule, suggesting adaptation to this environment. Altogether, these data contribute to explain why epidemiological clones tend to emerge from specific phylogenetic groups in the presence of pervasive horizontal gene transfer across the species.

Previous large scale studies on the evolution of E. coli focused on clinical isolates emphasizing virulence and antibiotic resistance in medically important lineages. Yet, most E. coli strains are either human commensals or not associated with humans at all. Here, we analyzed a large collection of non-clinical isolates of the species to assess the mechanisms of gene repertoire diversification in the light of isolation sources and phylogeny. We show that gene repertoires evolve so rapidly by the high turnover of mobile genetic elements that epidemiologically indistinguishable strains can be phenotypically extremely heterogeneous, illustrating the velocity of bacterial adaptation and the importance of accounting for the information on the whole genome at the epidemiological scale. Phylogeny and habitat shape the genetic diversification of E. coli to similar extents. Surprisingly, freshwater strains seem specifically adapted to this environment, breaking the paradigm that E. coli environmental isolates are systematically fecal contaminations. As a consequence, the evolution of this species is also shaped by environmental habitats, and it may diversify by acquiring genes and mobile elements from environmental bacteria (and not just from gut bacteria). This may facilitate the acquisition of virulence factors and antibiotic resistance in the strains that become pathogenic.

Study of a classification algorithm for AIEC identification in geographically distinct E. coli strains

Adherent-invasive Escherichia coli (AIEC) have been extensively implicated in Crohn’s disease pathogenesis. Currently, AIEC is identified phenotypically, since no molecular marker specific for AIEC exists. An algorithm based on single nucleotide polymorphisms was previously presented as a potential molecular tool to classify AIEC/non-AIEC, with 84% accuracy on a collection of 50 strains isolated in Girona (Spain). Herein, our aim was to determine the accuracy of the tool using AIEC/non-AIEC isolates from different geographical origins and extraintestinal pathogenic E. coli (ExPEC) strains. The accuracy of the tool was significantly reduced (61%) when external AIEC/non-AIEC strains from France, Chile, Mallorca (Spain) and Australia (82 AIEC, 57 non-AIEC and 45 ExPEC strains in total) were included. However, the inclusion of only the ExPEC strains showed that the tool was fairly accurate at differentiating these two close pathotypes (84.6% sensitivity; 79% accuracy). Moreover, the accuracy was still high (81%) for those AIEC/non-AIEC strains isolated from Girona and Mallorca (N = 63); two collections obtained from independent studies but geographically close. Our findings indicate that the presented tool is not universal since it would be only applicable for strains from similar geographic origin and demonstrates the need to include strains from different origins to validate such tools.

Clonal Structure, Virulence Factor-encoding Genes and Antibiotic Resistance of Escherichia coli, Causing Urinary Tract Infections and Other Extraintestinal Infections in Humans in Spain and France during 2016

Escherichia coli is the main pathogen responsible for extraintestinal infections. A total of 196 clinical E. coli consecutively isolated during 2016 in Spain (100 from Lucus Augusti hospital in Lugo) and France (96 from Beaujon hospital in Clichy) were characterized. Phylogroups, clonotypes, sequence types (STs), O:H serotypes, virulence factor (VF)-encoding genes and antibiotic resistance were determined. Approximately 10% of the infections were caused by ST131 isolates in both hospitals and approximately 60% of these infections were caused by isolates belonging to only 10 STs (ST10, ST12, ST58, ST69, ST73, ST88, ST95, ST127, ST131, ST141). ST88 isolates were frequent, especially in Spain, while ST141 isolates significantly predominated in France. The 23 ST131 isolates displayed four clonotypes: CH40-30, CH40-41, CH40-22 and CH40-298. Only 13 (6.6%) isolates were carriers of extended-spectrum beta-lactamase (ESBL) enzymes. However, 37.2% of the isolates were multidrug-resistant (MDR). Approximately 40% of the MDR isolates belonged to only four of the dominant clones (B2-CH40-30-ST131, B2-CH40-41-ST131, C-CH4-39-ST88 and D-CH35-27-ST69). Among the remaining MDR isolates, two isolates belonged to B2-CH14-64-ST1193, i.e., the new global emergent MDR clone. Moreover, a hybrid extraintestinal pathogenic E.coli (ExPEC)/enteroaggregative isolate belonging to the A-CH11-54-ST10 clone was identified.

The EnteroBase user's guide, with case studies on Salmonella transmissions, Yersinia pestis phylogeny, and Escherichia core genomic diversity

EnteroBase is an integrated software environment that supports the identification of global population structures within several bacterial genera that include pathogens. Here, we provide an overview of how EnteroBase works, what it can do, and its future prospects. EnteroBase has currently assembled more than 300,000 genomes from Illumina short reads from Salmonella, Escherichia, Yersinia, Clostridioides, Helicobacter, Vibrio, and Moraxella and genotyped those assemblies by core genome multilocus sequence typing (cgMLST). Hierarchical clustering of cgMLST sequence types allows mapping a new bacterial strain to predefined population structures at multiple levels of resolution within a few hours after uploading its short reads. Case Study 1 illustrates this process for local transmissions of Salmonella enterica serovar Agama between neighboring social groups of badgers and humans. EnteroBase also supports single nucleotide polymorphism (SNP) calls from both genomic assemblies and after extraction from metagenomic sequences, as illustrated by Case Study 2 which summarizes the microevolution of Yersinia pestis over the last 5000 years of pandemic plague. EnteroBase can also provide a global overview of the genomic diversity within an entire genus, as illustrated by Case Study 3, which presents a novel, global overview of the population structure of all of the species, subspecies, and clades within Escherichia.

Diverse Fluoroquinolone Resistance Plasmids From Retail Meat E. coli in the United States

Fluoroquinolones are used to treat serious bacterial infections, including those caused by Escherichia coli and Salmonella enterica. The emergence of plasmid-mediated quinolone resistance (PMQR) represent a new challenge to the successful treatment of Gram-negative infections. As part of a long-term strategy to generate a reference database of closed plasmids from antimicrobial resistant foodborne bacteria, we performed long-read sequencing of 11 E. coli isolates from retail meats that were non-susceptible to ciprofloxacin. Each of the isolates had PMQR genes, including qnrA1, qnrS1, and qnrB19. The four qnrB19 genes were carried on two distinct ColE-type plasmids among isolates from pork chop and ground turkey and were identical to plasmids previously identified in Salmonella. Seven other plasmids differed from any other sequences in GenBank and comprised IncF and IncR plasmids that ranged in size from 48 to 180 kb. These plasmids also contained different combinations of resistance genes, including those conferring resistance to beta-lactams, macrolides, sulfonamides, tetracycline, and heavy metals. Although relatively few isolates have PMQR genes, the identification of diverse plasmids in multiple retail meat sources suggests the potential for further spread of fluoroquinolone resistance, including through co-selection. These results highlight the value of long-read sequencing in characterizing antimicrobial resistance genes of public health concern.

Genotypic analysis of uropathogenic Escherichia coli to understand factors that impact the prevalence of β-lactam-resistant urinary tract infections in a community

OBJECTIVES

To determine the contribution of specific uropathogenic Escherichia coli (UPEC) lineages, drug resistance genes, and plasmid incompatibility/replicon (Inc) groups to the prevalence of β-lactam-resistant urinary tract infections (UTIs) in a university community.

METHODS

Urine samples were consecutively collected and cultured over a 2-year period from patients presenting to a university health centre with symptoms of UTI. Isolated UPEC were subtyped by multilocus sequence typing and fimH typing, and tested by PCR and sequencing for β-lactamase genes and plasmid Inc groups.

RESULTS

Among 273 UPEC isolates, 85 (31%) were ampicillin-resistant (AMP-R) and 188 (69%) were susceptible to all β-lactam drugs (AMP-S). Six lineages accounted for two-thirds of the isolates: ST95 (21%), ST69 (11%), ST420 (11%), ST73 (10%), ST127 (8%), and ST404 (3%). ST69 and ST404 were associated with AMP-R (P=0.003, P=0.0005), while ST420 and ST127 were associated with AMP-S (P<0.0001, P=0.027). ST95 contained four fimH types; the ST95/f-6 sublineage was more frequently identified among the AMP-R population (P = 0.009), while the ST95/f-47 sublineage was more frequently identified among the AMP-S population (P=0.007). The most common β-lactamase gene was bla_TEM, which was identified in 81 (95%) AMP-R isolates. IncFIB, IncFIA, and IncB/O type plasmids were the most commonly identified types, and were associated with β-lactam resistance (P<0.001 for all).

CONCLUSIONS

These observations indicate that the prevalence of β-lactam-resistant UTIs in this community was largely determined by a limited set of circulating UPEC STs and sublineages, carrying TEM β-lactamase genes that were likely encoded on one of three Inc type plasmids.

Rapid and Simple Universal Escherichia coli Genotyping Method Based on Multiple-Locus Variable-Number Tandem-Repeat Analysis Using Single-Tube Multiplex PCR and Standard Gel Electrophoresis

We developed a multiplex PCR method based on multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) that was designed for the rapid typing of Escherichia coli and Shigella isolates. The method amplifies seven VNTRs and does not require a sequencing capillary or fluorescent dyes. The amplification products are simply loaded on a standard agarose gel for electrophoresis, and the banding patterns are analyzed visually. We evaluated the method on 220 strains belonging to different collections: the E. coli reference (ECOR) collection (n = 72), O1:K1 isolates causing neonatal meningitis (n = 38), extended-spectrum beta-lactamase-producing fecal isolates belonging to the worldwide sequence type 131 (ST131) clone (n = 38), Shiga toxin-producing E. coli (STEC) isolates of serogroups O157:H7 (n = 21) and O26 (n = 16, 8 of which belonged to an outbreak), 27 Shigella isolates (22 Shigella sonnei isolates, including 5 epidemic strains), and 8 reference strains. The performances were compared to those of multilocus sequence typing (MLST), the DiversiLab automated repetitive element palindromic PCR (REP-PCR), pulsed-field gel electrophoresis (PFGE), and whole-genome sequencing (WGS). We found 66 different profiles among the isolates in the ECOR collection. Among the clonal group O1:K1 isolates, 14 different profiles were identified. For the 37 STEC isolates, we found 23 profiles, with 1 corresponding to the 8 epidemic strains. We found 19 profiles among the 27 Shigella isolates, with 1 corresponding to the epidemic strain. The method was able to recognize strains of the ST131 clone and to distinguish the O16 and O25b serogroups and identified 15 different MLVA types among them. This method allows the simple, fast, and inexpensive typing of E. coli/Shigella isolates that can be carried out in any laboratory equipped for molecular biology and has a discriminatory power superior to that of MLST and DiversiLab REP-PCR but slightly lower than that of PFGE.IMPORTANCE Fast typing methods that can easily and accurately distinguish clonal groups and unrelated isolates are of particular interest for microbiologists confronted with outbreaks or performing epidemiological studies. Highly discriminatory universal methods, like PFGE, optical mapping, or WGS, are expensive and/or time-consuming. MLST is useful for phylogeny but is less discriminatory and requires sequencing facilities. PCR methods, which are fast and easy to perform, also have drawbacks. Random PCRs and REP-PCR are universal but lack reproducibility. Other PCR methods may lack the discriminatory power to differentiate isolates during outbreaks. MLVA combines the advantages of PCR methods with a high discriminatory power but in its standard form requires sequencing capillary electrophoresis. The method that we have developed combines the advantages of standard PCR (simple, fast, and inexpensive) with the high discriminatory power of MLVA and permits the typing of all E. coli isolates (either intestinal or extraintestinal pathogenic isolates as well as commensal isolates).

Genetic structure, antimicrobial resistance and frequency of human associated Escherichia coli sequence types among faecal isolates from healthy dogs and cats living in Canberra, Australia

Extraintestinal pathogenic Escherichia coli (ExPEC) cause clinical infections in humans. Understanding the evolution and dissemination of ExPEC strains via potential reservoirs is important due to associated morbidity, health care costs and mortality. To further understanding this survey has examined isolates recovered from the faeces of 221 healthy dogs and 427 healthy cats. The distribution of phylogroups varied with host species, and depended on whether the animal was living in a shelter or a home. The human associated STs 69, 73, 95, 131 and 127 were prevalent, with 30.5% of cat isolates and 10.3% of dog isolates representing these ExPEC sequence types. Resistance to the antibiotics ampicillin and tetracycline was common, but resistance to other antimicrobials was negligible.

Genome sequencing of strains of the most prevalent clonal group of O1:K1:H7 Escherichia coli that causes neonatal meningitis in France

Background

To describe the temporal dynamics, molecular characterization, clinical and ex vivo virulence of emerging O1:K1 neonatal meningitis Escherichia coli (NMEC) strains of Sequence Type complex (STc) 95 in France.

The national reference center collected NMEC strains and performed whole genome sequencing (WGS) of O1:K1 STc95 NMEC strains for phylogenetic and virulence genes content analysis. Data on the clinical and biological features of patients were also collected. Ex vivo virulence was assessed using the Dictyostelium discoideum amoeba model.

Results

Among 250 NMEC strains collected between 1998 and 2015, 38 belonged to O1:K1 STc95. This clonal complex was the most frequently collected after 2004, representing up to 25% of NMEC strains in France. Phylogenetic analysis demonstrated that most (74%) belonged to a cluster designated D-1, characterized by the adhesin FimH30. There is no clinical data to suggest that this cluster is more pathogenic than its counterparts, although it is highly predominant and harbors a large repertoire of extraintestinal virulence factors, including a pS88-like plasmid. Ex vivo virulence model showed that this cluster was generally less virulent than STc95 reference strains of O45_S88:H7 and O18:H7 serotypes. However, the model showed differences between several subclones, although they harbor the same known virulence determinants.

Conclusions

The emerging clonal group O1:K1 STc95 of NMEC strains is mainly composed of a cluster with many virulence factors but of only moderate virulence. Whether its emergence is due to its ability to colonize the gut thanks to FimH30 or pS88-like plasmid remains to be determined.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1376-4) contains supplementary material, which is available to authorized users.

Discovery of mcr-1-Mediated Colistin Resistance in a Highly Virulent Escherichia coli Lineage

Resistance to last-line polymyxins mediated by the plasmid-borne mobile colistin resistance gene (mcr-1) represents a new threat to global human health. Here we present the complete genome sequence of an mcr-1-positive multidrug-resistant Escherichia coli strain (MS8345). We show that MS8345 belongs to serotype O2:K1:H4, has a large 241,164-bp IncHI2 plasmid that carries 15 other antibiotic resistance genes (including the extended-spectrum β-lactamase bla _CTX-M-1) and 3 putative multidrug efflux systems, and contains 14 chromosomally encoded antibiotic resistance genes. MS8345 also carries a large ColV-like virulence plasmid that has been associated with E. coli bacteremia. Whole-genome phylogeny revealed that MS8345 clusters within a discrete clade in the sequence type 95 (ST95) lineage, and MS8345 is very closely related to the highly virulent O45:K1:H4 clone associated with neonatal meningitis. Overall, the acquisition of a plasmid carrying resistance to colistin and multiple other antibiotics in this virulent E. coli lineage is concerning and might herald an era where the empirical treatment of ST95 infections becomes increasingly more difficult.IMPORTANCE Escherichia coli ST95 is a globally disseminated clone frequently associated with bloodstream infections and neonatal meningitis. However, the ST95 lineage is defined by low levels of drug resistance amongst clinical isolates, which normally provides for uncomplicated treatment options. Here, we provide the first detailed genomic analysis of an E. coli ST95 isolate that has both high virulence potential and resistance to multiple antibiotics. Using the genome, we predicted its virulence and antibiotic resistance mechanisms, which include resistance to last-line antibiotics mediated by the plasmid-borne mcr-1 gene. Finding an ST95 isolate resistant to nearly all antibiotics that also has a high virulence potential is of major clinical importance and underscores the need to monitor new and emerging trends in antibiotic resistance development in this important global lineage.

Molecular Epidemiology of Extraintestinal Pathogenic Escherichia coli

Extraintestinal pathogenic Escherichia coli (ExPEC) are important pathogens in humans and certain animals. Molecular epidemiological analyses of ExPEC are based on structured observations of E. coli strains as they occur in the wild. By assessing real-world phenomena as they occur in authentic contexts and hosts, they provide an important complement to experimental assessment. Fundamental to the success of molecular epidemiological studies are the careful selection of subjects and the use of appropriate typing methods and statistical analysis. To date, molecular epidemiological studies have yielded numerous important insights into putative virulence factors, host-pathogen relationships, phylogenetic background, reservoirs, antimicrobial-resistant strains, clinical diagnostics, and transmission pathways of ExPEC, and have delineated areas in which further study is needed. The rapid pace of discovery of new putative virulence factors and the increasing awareness of the importance of virulence factor regulation, expression, and molecular variation should stimulate many future molecular epidemiological investigations. The growing sophistication and availability of molecular typing methodologies, and of the new computational and statistical approaches that are being developed to address the huge amounts of data that whole genome sequencing generates, provide improved tools for such studies and allow new questions to be addressed.

Clonal Diversity, Virulence Potential and Antimicrobial Resistance of Escherichia coli Causing Community Acquired Urinary Tract Infection in Switzerland

Objectives: The aim of this study was to assess the clonal structure, virulence potential and antibiotic susceptibility of uropathogenic Escherichia coli (UPEC) isolates causing community acquired urinary tract infection (CAUTI) in unselected primary care patients in Switzerland. Methods: We performed multilocus sequence typing, virulence factor determination, and phenotypic and genotypic antimicrobial resistance testing on 44 non-duplicate UPEC isolates. Results: Twenty-seven different sequence types (STs) were identified. Major UPEC clones were represented by 19 (43.2%) of the isolates, including E. coli ST131, ST69 (both 13.6%), ST73 (6.8%), ST10 (4.5%), ST127, ST140, (both 2.3%). Five (11.4%) isolates belonged to ST141. Aggregate virulence factor (VF) scores were highest among isolates belonging to ST127 and ST141. Overall, 50% of the isolates were susceptible to all 12 antimicrobials tested, and all isolates remained susceptible to fosfomycin and nitrofurantoin. Resistance to sulfamethoxazole and ciprofloxacin were found in 31.8, and 15.9% of the isolates, respectively. Plasmid-mediated resistance genes were detected in ST69 and ST131 and included aac(6')-Ib-cr (2.3% of all isolates) blaCTX-M-14 and blaCTX-M-15 (9%), and mph(A) (13.6%). None of the isolates tested positive for mcr-1 or mcr-2. Conclusions: Our results show that CAUTI in Switzerland is caused by a wide variety of UPEC STs for which fosfomycin remains a good treatment option. We suggest that ST141 is an emerging clone associated with UTI in the community, and warrants closer attention. Moreover, the high rate of E. coli harboring mph(A) from patients without a history of antimicrobial therapy or hospitalization indicates that UPEC is an important reservoir for mph(A).

Diversity of the auxotrophic requirements in natural isolates of Escherichia coli

Isolates of Escherichia coli, except Shigella, are generally prototrophic; they do not require any growth factors to grow in mineral medium. However, a nicotinic acid requirement is common among B2 phylogroup STc95 O18 E. coli clone strains. Nicotinic acid is a precursor of nicotinamide adenine dinucleotide (NAD), an essential molecule that plays central role in cellular metabolism. The defect in NAD synthesis of these strains is due to alterations in de novo biosynthesis pathway nadB gene. Here, by studying growth on minimal medium with glycolytic (glucose) or gluconeogenic (pyruvate or succinate) substrates as the carbon supply in a large panel of E. coli natural isolates representative of the species diversity, we identify an absolute nicotinic acid requirement in non-STc95 strains due in one case to a nadA inactivation. The growth on glucose medium of some extraintestinal pathogenic E. coli strains belonging to various non-O18 B2 phylogroup STc95 clones is restored either by aspartate or nicotinate, demonstrating that the nicotinic acid requirement can also be due to an intracellular aspartate depletion. The auxotrophic requirements depend on the carbon source available in the environment. Moreover, some strains prototrophic in glucose medium become auxotrophic in succinate medium, and conversely, some strainsauxotrophic in glucose medium become prototrophic in succinate medium. Finally, a partial depletion of intracellular aspartate can be observed in some prototrophic strains belonging to various phylogroups. The observed more or less significant depletion according to isolates may be due to differences in tricarboxylic acid cycle enzyme activities. These metabolic defects could be involved in the adaptation of E. coli to its various niches.