Easy phylotyping of Escherichia coli via the EzClermont web app and command-line tool

Multi-omics strategy reveals potential role of antimicrobial resistance and virulence factor genes responsible for Simmental diarrheic calves caused by Escherichia coli

Escherichia coli (E. coli) is reported to be an important pathogen associated with calf diarrhea. Antibiotic resistance genes (ARGs) and virulence factor genes (VFGs) pose a considerable threat to both animal and human health. However, little is known about the characterization of ARGs and VFGs presented in the gut microbiota of diarrheic calves caused by E. coli. In this study, we used multi-omics strategy to analyze the ARG and VFG profiles of Simmental calves with diarrhea caused by E. coli K99. We found that gut bacterial composition and their microbiome metabolic functions varied greatly in diarrheic calves compared to healthy calves. In total, 175 ARGs were identified, and diarrheal calves showed a significantly higher diversity and abundance of ARGs than healthy calves. Simmental calves with diarrhea showed higher association of VFGs with pili function, curli assembly, and ferrienterobactin transport of E. coli. Co-occurrence patterns based on Pearson correlation analysis revealed that E. coli had a highly significant (P < 0.0001) correlation coefficient (>0.8) with 16 ARGs and 7 VFGs. Metabolomics analysis showed that differentially expressed metabolites in Simmental calves with diarrhea displayed a high correlation with the aforementioned ARGs and VFGs. Phylotype analysis of E. coli genomes showed that the predominant phylogroup B1 in diarrheic Simmental calves was associated with 10 ARGs and 3 VFGs. These findings provide an overview of the diversity and abundance of the gut microbiota in diarrheic calves caused by E. coli and pave the way for further studies on the mechanisms of antibiotic resistance and virulence in the calves affected with diarrhea.IMPORTANCESimmental is a well-recognized beef cattle breed worldwide. They also suffer significant economic losses due to diarrhea. In this study, fecal metagenomic analysis was applied to characterize the antibiotic resistance gene (ARG) and virulence factor gene (VFG) profiles of diarrheic Simmental calves. We identified key ARGs and VFGs correlated with Escherichia coli isolated from Simmental calves. Additionally, metabolomics analysis showed that differentially expressed metabolites in Simmental calves with diarrhea displayed a high correlation with the aforementioned ARGs and VFGs. Our findings provide an insight into the diversity and abundance of the gut microbiota in diarrheic calves caused by Escherichia coli and pave the way for further studies on the mechanisms of antibiotic resistance and virulence in the diarrheal calves from cattle hosts.

Genomic surveillance of antimicrobial-resistant Escherichia coli in fecal sludge and sewage in Uganda

The global increase of antimicrobial resistance (AMR) is a major public health concern. An effective AMR surveillance tool is needed to track the emergence and spread of AMR. Wastewater surveillance has been proposed as a resource-efficient tool for monitoring AMR carriage in the community. Here, we performed genomic surveillance of antimicrobial-resistant Escherichia coli obtained from fecal sludge and sewage in Uganda to gain insights into E. coli epidemiology and AMR burden in the underlying population. Selective media containing different antibiotic combinations (cefotaxime, ciprofloxacin, cefotaxime + ciprofloxacin + gentamicin) were used to obtain antimicrobial-resistant E. coli from fecal sludge and sewage. Short-read sequencing was performed for the obtained isolates, and a subset of isolates (selected from predominant sequence types (STs)) was also subjected to long-read sequencing. Genomic analysis of the obtained E. coli isolates (n = 181) revealed the prevalence of clonal complex 10, including ST167 (n = 43), ST10 (n = 28), ST1284 (n = 17), and ST617 (n = 4), in both fecal sludge and sewage, irrespective of antibiotics used for selection. We also detected global high-risk clones ST1193 (n = 10) and ST131 (n = 2 clade A, n = 3 subclade C1-M27, and n = 1 subclade C2). Diverse AMR determinants, including extended-spectrum β-lactamase genes (mostly blaCTX-M-15) and mutations in gyrA and parC, were identified. Analysis of the completed genomes revealed that diverse IncF plasmids and chromosomal integration were the major contributors to the spread of AMR genes in the predominant STs. This study showed that a combination of sewage surveillance (or fecal sludge surveillance) and whole-genome sequencing can be a powerful tool for monitoring AMR carriage in the underlying population.

Comprehensive probiogenomics analysis of the commensal Escherichia coli CEC15 as a potential probiotic strain

BACKGROUND

Probiotics have gained attention for their potential maintaining gut and immune homeostasis. They have been found to confer protection against pathogen colonization, possess immunomodulatory effects, enhance gut barrier functionality, and mitigate inflammation. However, a thorough understanding of the unique mechanisms of effects triggered by individual strains is necessary to optimize their therapeutic efficacy. Probiogenomics, involving high-throughput techniques, can help identify uncharacterized strains and aid in the rational selection of new probiotics. This study evaluates the potential of the Escherichia coli CEC15 strain as a probiotic through in silico, in vitro, and in vivo analyses, comparing it to the well-known probiotic reference E. coli Nissle 1917. Genomic analysis was conducted to identify traits with potential beneficial activity and to assess the safety of each strain (genomic islands, bacteriocin production, antibiotic resistance, production of proteins involved in host homeostasis, and proteins with adhesive properties). In vitro studies assessed survival in gastrointestinal simulated conditions and adhesion to cultured human intestinal cells. Safety was evaluated in BALB/c mice, monitoring the impact of E. coli consumption on clinical signs, intestinal architecture, intestinal permeability, and fecal microbiota. Additionally, the protective effects of both strains were assessed in a murine model of 5-FU-induced mucositis.

RESULTS

CEC15 mitigates inflammation, reinforces intestinal barrier, and modulates intestinal microbiota. In silico analysis revealed fewer pathogenicity-related traits in CEC15, when compared to Nissle 1917, with fewer toxin-associated genes and no gene suggesting the production of colibactin (a genotoxic agent). Most predicted antibiotic-resistance genes were neither associated with actual resistance, nor with transposable elements. The genome of CEC15 strain encodes proteins related to stress tolerance and to adhesion, in line with its better survival during digestion and higher adhesion to intestinal cells, when compared to Nissle 1917. Moreover, CEC15 exhibited beneficial effects on mice and their intestinal microbiota, both in healthy animals and against 5FU-induced intestinal mucositis.

CONCLUSIONS

These findings suggest that the CEC15 strain holds promise as a probiotic, as it could modulate the intestinal microbiota, providing immunomodulatory and anti-inflammatory effects, and reinforcing the intestinal barrier. These findings may have implications for the treatment of gastrointestinal disorders, particularly some forms of diarrhea.

Identification of plasmids in avian-associated Escherichia coli using nanopore and illumina sequencing

BACKGROUND

Avian pathogenic Escherichia coli (APEC) are the causative agents of colibacillosis in chickens, a disease which has significant economic impact on the poultry industry. Large plasmids detected in APEC are known to contribute to strain diversity for pathogenicity and antimicrobial resistance, but there could be other plasmids that are missed in standard analysis. In this study, we determined the impact of sequencing and assembly factors for the detection of plasmids in an E. coli whole genome sequencing project.

RESULTS

Hybrid assembly (Illumina and Nanopore) combined with plasmid DNA extractions allowed for detection of the greatest number of plasmids in E. coli, as detected by MOB-suite software. In total, 79 plasmids were identified in 19 E. coli isolates. Hybrid assemblies were robust and consistent in quality regardless of sequencing kit used or if long reads were filtered or not. In contrast, long read only assemblies were more variable and influenced by sequencing and assembly parameters. Plasmid DNA extractions allowed for the detection of physically smaller plasmids, but when averaged over 19 isolates did not significantly change the overall number of plasmids detected.

CONCLUSIONS

Hybrid assembly can be reliably used to detect plasmids in E. coli, especially if researchers are focused on large plasmids containing antimicrobial resistance genes and virulence factors. If the goal is comprehensive detection of all plasmids, particularly if smaller sized vectors are desired for biotechnology applications, the addition of plasmid DNA extractions to hybrid assemblies is prudent. Long read sequencing is sufficient to detect many plasmids in E. coli, however, it is more prone to errors when expanded to analyze a large number of isolates.

Genomic analysis of phylogroup D Escherichia coli strains using novel de-novo reference-based guided assembly

Escherichia coli are highly diverse bacteria with different pathogenic types, serotypes and phylogenetic types/phylotypes. In recent years, infections with E. coli have increased worldwide and so has the emergence of antibiotic resistant strains. In the present study we have assembled, annotated and analysed genome sequences of three strains of the phylogroup D of E. coli. These strains were isolated from the river Yamuna, a prominent anthropogenic urban river of northern India. These strains showed varied antibiotic susceptibilities, one was susceptible to all the antibiotics tested except ampicillin while of the other two, one was multi-β-lactam resistant and the other was multi-drug resistant (resistant to multiple β-lactams, fluoroquinolones and kanamycin). The short-sequence reads were assembled into contigs using the de-novo approach and further, scaffolding of contigs was performed by using the best reference genome for a particular isolate which resulted in a significant increase in the N50 value of each assembly. The bioinformatics assembly approach used in this study could be easily applied to study other bacterial genomes.

Commensal Escherichia coli Strains of Bovine Origin Competitively Mitigated Escherichia coli O157:H7 in a Gnotobiotic Murine Intestinal Colonization Model with or without Physiological Stress

Cattle are a primary reservoir of enterohemorrhagic Escherichia coli (EHEC) O157:H7. Currently, there are no effective methods of eliminating this important zoonotic pathogen from cattle, and colonization resistance in relation to EHEC O157:H7 in cattle is poorly understood. We developed a gnotobiotic EHEC O157:H7 murine model to examine aspects of the cattle pathogen-microbiota interaction, and to investigate competitive suppression of EHEC O157:H7 by 18 phylogenetically distinct commensal E. coli strains of bovine origin. As stress has been suggested to influence enteric colonization by EHEC O157:H7 in cattle, corticosterone administration (±) to incite a physiological stress response was included as an experimental variable. Colonization of the intestinal tract (IT) of mice by the bovine EHEC O157:H7 strain, FRIK-2001, mimicked characteristics of bovine IT colonization. In this regard, FRIK-2001 successfully colonized the IT and temporally incited minimal impacts on the host relative to other EHEC O157:H7 strains, including on the renal metabolome. The presence of the commensal E. coli strains decreased EHEC O157:H7 densities in the cecum, proximal colon, and distal colon. Moreover, histopathologic changes and inflammation markers were reduced in the distal colon of mice inoculated with commensal E. coli strains (both propagated separately and communally). Although stress induction affected the behavior of mice, it did not influence EHEC O157:H7 densities or disease. These findings support the use of a gnotobiotic murine model of enteric bovine EHEC O157:H7 colonization to better understand pathogen-host-microbiota interactions toward the development of effective on-farm mitigations for EHEC O157:H7 in cattle, including the identification of bacteria capable of competitively colonizing the IT.

Genomic analysis of multidrug-resistant Escherichia coli strains carrying the mcr-1 gene recovered from pigs in Lima-Peru

Antibiotic resistance is a current problem that significantly impacts overall health. The dissemination of antibiotic resistance genes (ARGs) to urban areas primarily occurs through ARG-carrying bacteria present in the gut microbiota of animals raised in intensive farming settings, such as pig production. Hence, this study aimed to isolate and analyzed 87 Escherichia coli strains from pig fecal samples obtained from intensive farms in Lima Department. The isolates were subjected to Kirby-Bauer-Disk Diffusion Test and PCR for mcr-1 gene identification. Disk-diffusion assay revealed a high level of resistance among these isolates to oxytetracycline, ampicillin, cephalothin, chloramphenicol, ciprofloxacin, and doxycycline. PCR analysis identified the mcr-1 gene in 8% (7/87) E. coli isolates. Further, whole genome sequencing was conducted on 17 isolates, including multidrug resistance (MDR) E. coli and/or mcr-1 gene carriers. This analysis unveiled a diverse array of ARGs. Alongside the mcr-1 gene, the bla_CTX-M55 gene was particularly noteworthy as it confers resistance to third generation cephalosporins, including ceftriaxone. MDR E. coli genomes exhibited other ARGs encoding resistance to fosfomycin (fosA3), quinolones (qnrB19, qnrS1, qnrE1), tetracyclines (tetA, tetB, tetD, tetM), sulfonamides (sul1, sul2, sul3), amphenicols (cmlA1, floR), lincosamides (inuE), as well as various aminoglycoside resistance genes. Additionally, Multi Locus Sequence Typing (MLST) revealed a high diversity of E. coli strains, including ST10, a pandemic clone. This information provides evidence of the dissemination of highly significant ARGs in public health. Therefore, it is imperative to implement measures aimed at mitigating and preventing the transmission of MDR bacteria carrying ARGs to urban environments.

Dominance of Escherichia coli sequence types ST73, ST95, ST127 and ST131 in Australian urine isolates: a genomic analysis of antimicrobial resistance and virulence linked to F plasmids

Extraintestinal pathogenic Escherichia coli (ExPEC) are the most frequent cause of urinary tract infections (UTIs) globally. Most studies of clinical E. coli isolates are selected based on their antimicrobial resistance (AMR) phenotypes; however, this selection bias may not provide an accurate portrayal of which sequence types (STs) cause the most disease. Here, whole genome sequencing (WGS) was performed on 320 E. coli isolates from urine samples sourced from a regional hospital in Australia in 2006. Most isolates (91%) were sourced from patients with UTIs and were not selected based on any AMR phenotypes. No significant differences were observed in AMR and virulence genes profiles across age sex, and uro-clinical syndromes. While 88 STs were identified, ST73, ST95, ST127 and ST131 dominated. F virulence plasmids carrying senB-cjrABC (126/231; 55%) virulence genes were a feature of this collection. These senB-cjrABC+ plasmids were split into two categories: pUTI89-like (F29:A-:B10 and/or >95 % identity to pUTI89) (n=73) and non-pUTI89-like (n=53). Compared to all other plasmid replicons, isolates with pUTI89-like plasmids carried fewer antibiotic resistance genes (ARGs), whilst isolates with senB-cjrABC+/non-pUTI89 plasmids had a significantly higher load of ARGs and class 1 integrons. F plasmids were not detected in 89 genomes, predominantly ST73. Our phylogenomic analyses identified closely related isolates from the same patient associated with different pathologies and evidence of strain-sharing events involving isolates sourced from companion and wild animals.

Comparative Genomic Analysis of ST131 Subclade C2 of ESBL-Producing E. coli Isolates from Patients with Recurrent and Sporadic Urinary Tract Infections

The global emergence of extended-spectrum beta-lactamase-producing Escherichia coli (ESBL-E. coli), mainly causing urinary tract infections (UTI), is a major threat to human health. ESBL-E. coli sequence type (ST) 131 is the dominating clone worldwide, especially its subclade C2. Patients developing recurrent UTI (RUTI) due to ST131 subclade C2 appear to have an increased risk of recurrent infections. We have thus compared the whole genome of ST131 subclade C2 isolates from 14 patients with RUTI to those from 14 patients with sporadic UTI (SUTI). We aimed to elucidate if isolates causing RUTI can be associated with specific genomic features. Paired isolates from patients with RUTI were identical, presenting 2-18 single nucleotide polymorphism (SNP) differences for all six patients investigated. Comparative genomic analyses, including virulence factors, antibiotic resistance, pangenome and SNP analyses did not find any pattern associated with isolates causing RUTI. Despite extensive whole genome analyses, an increased risk of recurrences seen in patients with UTI due to ST131 subclade C2 isolates could not be explained by bacterial genetic differences in the two groups of isolates. Hence, additional factors that could aid in identifying bacterial properties contributing to the increased risk of RUTI due to ESBL-E. coli ST131 subclade C2 remains to be explored.

Whole Genome Sequencing of Avian Pathogenic Escherichia coli Causing Bacterial Chondronecrosis and Osteomyelitis in Australian Poultry

Bacterial chondronecrosis with osteomyelitis (BCO) impacts animal welfare and productivity in the poultry industry worldwide, yet it has an understudied pathogenesis. While Avian Pathogenic Escherichia coli (APEC) are known to be one of the main causes, there is a lack of whole genome sequence data, with only a few BCO-associated APEC (APEC_BCO) genomes available in public databases. In this study, we conducted an analysis of 205 APEC_BCO genome sequences to generate new baseline phylogenomic knowledge regarding the diversity of E. coli sequence types and the presence of virulence associated genes (VAGs). Our findings revealed the following: (i) APEC_BCO are phylogenetically and genotypically similar to APEC that cause colibacillosis (APEC_colibac), with globally disseminated APEC sequence types ST117, ST57, ST69, and ST95 being predominate; (ii) APEC_BCO are frequent carriers of ColV-like plasmids that carry a similar set of VAGs as those found in APEC_colibac. Additionally, we performed genomic comparisons, including a genome-wide association study, with a complementary collection of geotemporally-matched genomes of APEC from multiple cases of colibacillosis (APEC_colibac). Our genome-wide association study found no evidence of novel virulence loci unique to APEC_BCO. Overall, our data indicate that APEC_BCO and APEC_colibac are not distinct subpopulations of APEC. Our publication of these genomes substantially increases the available collection of APEC_BCO genomes and provides insights for the management and treatment strategies of lameness in poultry.

Functional Genomics of Gastrointestinal Escherichia coli Isolated from Patients with Cancer and Diarrhea

We describe the epidemiology and clinical characteristics of 29 patients with cancer and diarrhea in whom Enteroaggregative Escherichia coli (EAEC) was initially identified by GI BioFire panel multiplex. E. coli strains were successfully isolated from fecal cultures in 14 of 29 patients. Six of the 14 strains were identified as EAEC and 8 belonged to other diverse E. coli groups of unknown pathogenesis. We investigated these strains by their adherence to human intestinal organoids, cytotoxic responses, antibiotic resistance profile, full sequencing of their genomes, and annotation of their functional virulome. Interestingly, we discovered novel and enhanced adherence and aggregative patterns for several diarrheagenic pathotypes that were not previously seen when co-cultured with immortalized cell lines. EAEC isolates displayed exceptional adherence and aggregation to human colonoids compared not only to diverse GI E. coli , but also compared to prototype strains of other diarrheagenic E. coli . Some of the diverse E. coli strains that could not be classified as a conventional pathotype also showed an enhanced aggregative and cytotoxic response. Notably, we found a high carriage rate of antibiotic resistance genes in both EAEC strains and diverse GI E. coli isolates and observed a positive correlation between adherence to colonoids and the number of metal acquisition genes carried in both EAEC and the diverse E. coli strains. This work indicates that E. coli from cancer patients constitute strains of remarkable pathotypic and genomic divergence, including strains of unknown disease etiology with unique virulomes. Future studies will allow for the opportunity to re-define E. coli pathotypes with greater diagnostic accuracy and into more clinically relevant groupings.

Escherichia cryptic clade I is an emerging source of human intestinal pathogens

BACKGROUND

Within the genus Escherichia, several monophyletic clades other than the traditionally defined species have been identified. Of these, cryptic clade I (C-I) appears to represent a subspecies of E. coli, but due to the difficulty in distinguishing it from E. coli sensu stricto, the population structure and virulence potential of C-I are unclear.

RESULTS

We defined a set of true C-I strains (n = 465), including a Shiga toxin 2a (Stx2a)-producing isolate from a patient with bloody diarrhoea identified by the retrospective analyses using a C-I-specific detection system. Through genomic analysis of 804 isolates from the cryptic clades, including these C-I strains, we revealed their global population structures and the marked accumulation of virulence genes and antimicrobial resistance genes in C-I. In particular, half of the C-I strains contained hallmark virulence genes of Stx-producing E. coli (STEC) and/or enterotoxigenic E. coli (ETEC). We also found the host-specific distributions of virulence genes, which suggests bovines as the potential source of human infections caused by STEC- and STEC/ETEC hybrid-type C-I strains, as is known in STEC.

CONCLUSIONS

Our findings demonstrate the emergence of human intestinal pathogens in C-I lineage. To better understand the features of C-I strains and their infections, extensive surveillance and larger population studies of C-I strains are needed. The C-I-specific detection system developed in this study will be a powerful tool for screening and identifying C-I strains.

Survey of clinical and commensal Escherichia coli from commercial broilers and turkeys, with emphasis on high-risk clones using APECTyper

Molecular characterization of avian pathogenic Escherichia coli (APEC) is challenging due to the complex nature of its associated disease, colibacillosis, in poultry. Numerous efforts have been made toward defining APEC, and it is becoming clear that certain clonal backgrounds are predictive of an avian E. coli isolate's virulence potential. Thus, APEC can be further differentiated as high-risk APEC based upon their clonal background's virulence potential. However, less clear is the degree of overlap between clinical isolates of differing bird type, and between clinical and gastrointestinal isolates. This study aimed to determine genomic similarities and differences between such populations, comparing commercial broiler vs. turkey isolates, and clinical vs. gastrointestinal isolates. Differences were observed in Clermont phylogenetic groups between isolate populations, with B2 as the dominant group in turkey clinical isolates and G as the dominant group in broiler clinical isolates. Nearly all clinical isolates were classified as APEC using a traditional gene-based typing scheme, whereas 53.4% and 44.1% of broiler and turkey gastrointestinal isolates were classified as APEC, respectively. High-risk APEC were identified among 31.0% and 46.9% of broiler and turkey clinical isolates, compared with 5.7% and 2.9% of broiler and turkey gastrointestinal isolates. As found in previous studies, no specific known virulence or fitness gene sets were identified which universally differentiate between clinical and gastrointestinal isolates. This study further demonstrates the utility of a hybrid APEC typing approach, considering both plasmid content and clonal background, for the identification of dominant and highly virulent APEC clones in poultry production.

Strain-level bacterial typing directly from patient samples using optical DNA mapping

BACKGROUND

Identification of pathogens is crucial to efficiently treat and prevent bacterial infections. However, existing diagnostic techniques are slow or have a too low resolution for well-informed clinical decisions.

METHODS

In this study, we have developed an optical DNA mapping-based method for strain-level bacterial typing and simultaneous plasmid characterisation. For the typing, different taxonomical resolutions were examined and cultivated pure Escherichia coli and Klebsiella pneumoniae samples were used for parameter optimization. Finally, the method was applied to mixed bacterial samples and uncultured urine samples from patients with urinary tract infections.

RESULTS

We demonstrate that optical DNA mapping of single DNA molecules can identify Escherichia coli and Klebsiella pneumoniae at the strain level directly from patient samples. At a taxonomic resolution corresponding to E. coli sequence type 131 and K. pneumoniae clonal complex 258 forming distinct groups, the average true positive prediction rates are 94% and 89%, respectively. The single-molecule aspect of the method enables us to identify multiple E. coli strains in polymicrobial samples. Furthermore, by targeting plasmid-borne antibiotic resistance genes with Cas9 restriction, we simultaneously identify the strain or subtype and characterize the corresponding plasmids.

CONCLUSION

The optical DNA mapping method is accurate and directly applicable to polymicrobial and clinical samples without cultivation. Hence, it has the potential to rapidly provide comprehensive diagnostics information, thereby optimizing early antibiotic treatment and opening up for future precision medicine management.

The enterohemorrhagic Escherichia coli insertion sequence-excision enhancer protein is a DNA polymerase with microhomology-mediated end-joining activity

Bacterial genomes contain an abundance of transposable insertion sequence (IS) elements that are essential for genome evolution and fitness. Among them, IS629 is present in most strains of enterohemorrhagic Escherichia coli O157 and accounts for many polymorphisms associated with gene inactivation and/or genomic deletions. The excision of IS629 from the genome is promoted by IS-excision enhancer (IEE) protein. Despite IEE has been identified in the most pathogenic serotypes of E. coli, its biochemical features that could explain its role in IS excision are not yet understood. We show that IEE is present in >30% of all available E. coli genome assemblies, and is highly conserved and very abundant within enterohemorrhagic, enteropathogenic and enterotoxigenic genomes. In vitro analysis of the recombinant protein from E. coli O157:H7 revealed the presence of a Mn2+-dependent error-prone DNA polymerase activity in its N-terminal archaeo-eukaryotic primase (AEP) domain able to promote dislocations of the primer and template strands. Importantly, IEE could efficiently perform in vitro an end-joining reaction of 3'-single-strand DNA overhangs with ≥4 bp of homology requiring both the N-terminal AEP and C-terminal helicase domains. The proposed role for IEE in the novel IS excision mechanism is discussed.

Escherichia coli O80 in Healthy Cattle: Absence of Shigatoxigenic and Enteropathogenic E. coli O80:H2 and (Phylo) Genomics of Non-Clonal Complex 165 E. coli O80

The origin of human and calf infections by Shigatoxigenic (STEC) and enteropathogenic (EPEC) Escherichia coli O80:H2 is still unknown. The aim of this study was to identify E. coli O80 in healthy cattle with an emphasis on melibiose non-fermenting E. coli O80:H2. Faecal materials collected from 149 bulls at 1 slaughterhouse and 194 cows on 9 farms were tested with O80 antigen-encoding gene PCR after overnight growth in enrichment broths. The 53 O80 PCR-positive broths were streaked on different (semi-)selective agar plates. Five E. coli colonies from 3 bulls and 11 from 2 cows tested positive with the O80 PCR, but no melibiose non-fermenting E. coli was isolated. However, these 16 E. coli O80 were negative with PCR targeting the fliC_H2, eae, stx1, stx2 and hlyF genes and were identified by WGS to serotypes and sequence types O80:H6/ST8619 and O80:H45/ST4175. They were phylogenetically related to E. coli O80:H6 and O80:H45 isolated from different animal species in different countries, respectively, but neither to STEC and EPEC O80:H2/ST301, nor to other serotypes of the clonal complex 165. As a conclusion, healthy adult cattle were not identified as a source of contamination of humans and calves by STEC or EPEC O80:H2.

A one-year genomic investigation of Escherichia coli epidemiology and nosocomial spread at a large US healthcare network

BACKGROUND

Extra-intestinal pathogenic Escherichia coli (ExPEC) are a leading cause of bloodstream and urinary tract infections worldwide. Over the last two decades, increased rates of antibiotic resistance in E. coli have been reported, further complicating treatment. Worryingly, specific lineages expressing extended-spectrum β-lactamases (ESBLs) and fluoroquinolone resistance have proliferated and are now considered a serious threat. Obtaining contemporary information on the epidemiology and prevalence of these circulating lineages is critical for containing their spread globally and within the clinic.

METHODS

Whole-genome sequencing (WGS), phylogenetic analysis, and antibiotic susceptibility testing were performed for a complete set of 2075 E. coli clinical isolates collected from 1776 patients at a large tertiary healthcare network in the USA between October 2019 and September 2020.

RESULTS

The isolates represented two main phylogenetic groups, B2 and D, with six lineages accounting for 53% of strains: ST-69, ST-73, ST-95, ST-131, ST-127, and ST-1193. Twenty-seven percent of the primary isolates were multidrug resistant (MDR) and 5% carried an ESBL gene. Importantly, 74% of the ESBL-E.coli were co-resistant to fluoroquinolones and mostly belonged to pandemic ST-131 and emerging ST-1193. SNP-based detection of possible outbreaks identified 95 potential transmission clusters totaling 258 isolates (12% of the whole population) from ≥ 2 patients. While the proportion of MDR isolates was enriched in the set of putative transmission isolates compared to sporadic infections (35 vs 27%, p = 0.007), a large fraction (61%) of the predicted outbreaks (including the largest cluster grouping isolates from 12 patients) were caused by the transmission of non-MDR clones.

CONCLUSION

By coupling in-depth genomic characterization with a complete sampling of clinical isolates for a full year, this study provides a rare and contemporary survey on the epidemiology and spread of E. coli in a large US healthcare network. While surveillance and infection control efforts often focus on ESBL and MDR lineages, our findings reveal that non-MDR isolates represent a large burden of infections, including those of predicted nosocomial origins. This increased awareness is key for implementing effective WGS-based surveillance as a routine technology for infection control.

Occurrence and Genomic Characterization of mcr-1-Harboring Escherichia coli Isolates from Chicken and Pig Farms in Lima, Peru

Resistance to colistin generated by the mcr-1 gene in Enterobacteriaceae is of great concern due to its efficient worldwide spread. Despite the fact that the Lima region has a third of the Peruvian population and more than half of the national pig and poultry production, there are no reports of the occurrence of the mcr-1 gene in Escherichia coli isolated from livestock. In the present work, we studied the occurrence of E. coli carrying the mcr-1 gene in chicken and pig farms in Lima between 2019 and 2020 and described the genomic context of the mcr-1 gene. We collected fecal samples from 15 farms in 4 provinces of Lima including the capital Lima Metropolitana and recovered 341 E. coli isolates. We found that 21.3% (42/197) and 12.5% (18/144) of the chicken and pig strains were mcr-1-positive by PCR, respectively. The whole genome sequencing of 14 mcr-1-positive isolates revealed diverse sequence types (e.g., ST48 and ST602) and the presence of other 38 genes that confer resistance to 10 different classes of antibiotics, including beta-lactamase bla_CTX-M-55. The mcr-1 gene was located on diverse plasmids belonging to the IncI2 and IncHI1A:IncHI1B replicon types. A comparative analysis of the plasmids showed that they contained the mcr-1 gene within varied structures (mikB-mcr1-pap2, ISApl1-mcr1-pap2, and Tn6330). To the best of our knowledge, this is the first attempt to study the prevalence of the mcr-1 gene in livestock in Peru, revealing its high occurrence in pig and chicken farms. The genetic diversity of mcr-1-positive strains suggests a complex local epidemiology calling for a coordinated surveillance under the One-Health approach that includes animals, retail meat, farmers, hospitals and the environment to effectively detect and limit the spread of colistin-resistant bacteria.

Virulence Potential of Escherichia coli Isolates from a Beef Farm in Puerto Rico

Sixteen Escherichia coli isolates were obtained from fecal matter from a beef farm in Puerto Rico. Isolates were whole-genome sequenced for in silico characterization, including pathotype characterization, virulence, and plasmid identification. The results of the draft genomes identified potential pathogenic E. coli strains from beef cattle in Puerto Rico.

Escherichia coli in different animal feces: phylotypes and virulence genes

In this study, it was aimed to determine the phylogroups of Escherichia coli isolates from horse, cat, dog, sheep, cattle, and chicken feces samples and to investigate some important virulence genes of the isolates. For this purpose, a total of 600 feces samples, 100 from each animal species, were used as material. For the isolation of E.coli, feces samples were directly inoculated on MacConkey agar. The identification of the isolates was performed via phenotypic tests and species-specific multiplex Polymerase Chain Reaction (mPCR) method. PCR methods were used to phylotype E.coli isolates and to investigate virulence genes (bfpA, eaeA, LT, ST, Stx1, and Stx2). Of the total 600 E.coli isolates recovered in this study, 120 (20%), 269 (44.8%), 58 (9.7%), 19 (3.2%), 35 (5.8%), 56 (9.3%), 31 (5.2%), and 12 (2%) were identified as phylogroup A, B1, B2, C, D, E, F, and Escherichia clade I, respectively. While the virulence gene was detected in 149 (24.8%) E.coli isolates, no virulence gene was detected in 451 (75.2%) isolates. According to the analysis results, the most determined virulence gene was Stx1, while the least determined virulence gene was LT. In conclusion, in this study, when both the animal species and the number of E.coli isolates examined are considered, the data obtained are of great importance in epidemiological terms. However, the detection of virulence genes in 13.5% among phylogroup A, B1, and C isolates with commensal characteristics suggest that these isolates may show pathogenic characteristics with the virulence genes they contain.

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'.

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

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

A comparative genomics approach for identifying genetic factors in Escherichia coli isolates associated with bovine diseases

AIMS

Escherichia coli are ubiquitously present bacterial pathogens that cause septicaemia, diarrhoea and other clinical illness in farm animals. Many pathogen factors can be associated with disease conditions. Currently, studies inferring E. coli genetic factors associated with infection in bovines are limited. Hence, the present study envisaged to determine the pathogen genetic factors associated with bovine disease conditions.

METHOD AND RESULTS

The comparative genomic analysis involved genome sequence data of 135 diseased and 145 healthy bovine origin E. coli strains. Phylogroups A and C, as well as pathotypes ExPEC and EPEC, were found to have a strong connection with bovine disease strains. STEC strains, including EHEC, seem to play a less important role in bovine disease. Sequence types (STs) predominant among strains from diarrhoeal origin were ST 301 (CC 165) and ST 342. Correlation of core genome phylogeny with accessory gene based clustering, phylogroups and pathotypes indicated lineage specific virulence factors mostly associated with disease conditions.

CONCLUSIONS

Comparative genomic analysis was applied to infer genetic factors significant in bovine disease origin E. coli strains. Isolates from bovine disease origin were enriched for the phylogroups A and C, and for the pathotypes ExPEC and EPEC. However, there was minimal evidence of STEC involvement. The study also indicated predominant genetic lineages and virulence genes (pap, sfa and afa) associated with disease origin strains. SIGNIFICANCE AND IMPACT OF STUDY;: The study revealed significant pathotypes, phylgroups, serotypes and sequence types associated with bovine disease conditions. These identified genetic factors can be applied for disease diagnosis, implementing vaccine and therapeutic measures. In addition, E. coli isolates from the bovine species revealed a complex pattern of disease epidemiology.

Characterization of Adherent-Invasive Escherichia coli (AIEC) Outer Membrane Proteins Provides Potential Molecular Markers to Screen Putative AIEC Strains

Adherent-invasive E. coli (AIEC) is a pathotype associated with the etiopathogenesis of Crohn's disease (CD), albeit with an as-yet unclear role. The main pathogenic mechanisms described for AIEC are adherence to epithelial cells, invasion of epithelial cells, and survival and replication within macrophages. A few virulence factors have been described as participating directly in these phenotypes, most of which have been evaluated only in AIEC reference strains. To date, no molecular markers have been identified that can differentiate AIEC from other E. coli pathotypes, so these strains are currently identified based on the phenotypic characterization of their pathogenic mechanisms. The identification of putative AIEC molecular markers could be beneficial not only from the diagnostic point of view but could also help in better understanding the determinants of AIEC pathogenicity. The objective of this study was to identify molecular markers that contribute to the screening of AIEC strains. For this, we characterized outer membrane protein (OMP) profiles in a group of AIEC strains and compared them with the commensal E. coli HS strain. Notably, we found a set of OMPs that were present in the AIEC strains but absent in the HS strain. Moreover, we developed a PCR assay and performed phylogenomic analyses to determine the frequency and distribution of the genes coding for these OMPs in a larger collection of AIEC and other E. coli strains. As result, it was found that three genes (chuA, eefC, and fitA) are widely distributed and significantly correlated with AIEC strains, whereas they are infrequent in commensal and diarrheagenic E. coli strains (DEC). Additional studies are needed to validate these markers in diverse strain collections from different geographical regions, as well as investigate their possible role in AIEC pathogenicity.

Corrigendum: Extended-Spectrum Beta-Lactamase Producing-Escherichia coli Isolated From Irrigation Waters and Produce in Ecuador

[This corrects the article DOI: 10.3389/fmicb.2021.709418.].

High Prevalence of bla CTXM-1/IncI1-Iγ/ST3 Plasmids in Extended-Spectrum β-Lactamase-Producing Escherichia coli Isolates Collected From Domestic Animals in Guadeloupe (French West Indies)

Extended-spectrum β-lactamase-producing Enterobacteriaceae (ESBL-E) have been classified in the group of resistant bacteria of highest priority. We determined the prevalence of ESBL-E collected in feces from household and shelter pets in Guadeloupe (French West Indies). A single rectal swab was taken from 125 dogs and 60 cats between June and September 2019. The prevalence of fecal carriage of ESBL-E was 7.6% (14/185, 95% CI: 4.2-12.4), within the range observed worldwide. The only risk factor associated with a higher prevalence of ESBL-E rectal carriage was a stay in a shelter, suggesting that refuges could be hotspots for their acquisition. All but one (Klebsiella pneumoniae from a cat) were Escherichia coli. We noted the presence of a bla _CTX-M-1/IncI1-Iγ/sequence type (ST3) plasmid in 11 ESBL-producing E. coli isolates belonging to ST328 (n = 6), ST155 (n = 4) and ST953 (n = 1). A bla _CTX-M-15 gene was identified in the three remaining ESBL-E isolates. The bla _CTX-M-1 and most of the antimicrobial resistance genes were present in a well-conserved large conjugative IncI1-Iγ/ST3 plasmid characterized by two accessory regions containing antibiotic resistance genes. The plasmid has been detected worldwide in E. coli isolates from humans and several animal species, such as food-producing animals, wild birds and pets, and from the environment. This study shows the potential role of pets as a reservoir of antimicrobial-resistant bacteria or genes for humans and underlines the importance of basic hygiene measures by owners of companion animals.

Escherichia coli and Pseudomonas aeruginosa Isolated From Urine of Healthy Bovine Have Potential as Emerging Human and Bovine Pathogens

The study of livestock microbiota has immediate benefits for animal health as well as mitigating food contamination and emerging pathogens. While prior research has indicated the gastrointestinal tract of cattle as the source for many zoonoses, including Shiga-toxin producing Escherichia coli and antibiotic resistant bacteria, the bovine urinary tract microbiota has yet to be thoroughly investigated. Here, we describe 5 E. coli and 4 Pseudomonas aeruginosa strains isolated from urine of dairy Gyr cattle. While both species are typically associated with urinary tract infections and mastitis, all of the animals sampled were healthy. The bovine urinary strains were compared to E. coli and P. aeruginosa isolates from other bovine samples as well as human urinary samples. While the bovine urinary E. coli isolates had genomic similarity to isolates from the gastrointestinal tract of cattle and other agricultural animals, the bovine urinary P. aeruginosa strains were most similar to human isolates suggesting niche adaptation rather than host adaptation. Examination of prophages harbored by these bovine isolates revealed similarity with prophages within distantly related E. coli and P. aeruginosa isolates from the human urinary tract. This suggests that related urinary phages may persist and/or be shared between mammals. Future studies of the bovine urinary microbiota are needed to ascertain if E. coli and P. aeruginosa are resident members of this niche and/or possible sources for emerging pathogens in humans.

Genetic and Ecological Diversity of Escherichia coli and Cryptic Escherichia Clades in Subtropical Aquatic Environments

Escherichia coli not only inhabit the large intestines of human and warm-blooded animals but could also persist in the external environment. However, current knowledge was largely based on host-associated strains. Moreover, cryptic Escherichia clades that were often misidentified as E. coli by conventional diagnostic methods were discovered. Failure to distinguish them from E. coli sensu stricto could lead to inaccurate conclusions about the population genetics of E. coli. Based on seven housekeeping genes, we determine the genetic and ecological diversity of E. coli and cryptic clades as they occupy aquatic habitats with different characteristics and human impact levels in subtropical Hong Kong. Contrary to previous reports, clade II was the most abundant cryptic lineage co-isolated with E. coli, being especially abundant in relatively pristine subtropical aquatic environments. The phylogenetically distinct cryptic clades and E. coli showed limited recombination and significant genetic divergence. Analyses indicated that these clade II strains were ecologically differentiated from typical E. coli; some may even represent novel environmental Escherichia clades that were closely related to the original clade II strains of fecal origins. E. coli of diverse origins exhibited clonality amidst divergent genotypes STs, echoing other studies in that recombination in housekeeping genes was insufficient to disrupt phylogenetic signals of the largely clonal E. coli. Notably, environmental E. coli were less diverse than fecal isolates despite contributing many new alleles and STs. Finally, we demonstrated that human activities influenced the distribution of E. coli and clade II in a small aquatic continuum. Moving from relatively pristine sites toward areas with higher human disturbance, the abundance of clade II isolates and new E. coli genotypes reduces, while E. coli bearing class I integrons and belonging to CCs of public health concern accumulates. Altogether, this work revealed the new genetic diversity of E. coli and cryptic clades embedded in selected subtropical aquatic habitats, especially relatively pristine sites, which will aid a more thorough understanding of the extent of their genetic and functional variations in relation to diverse habitats with varied conditions.

Persistence and Dissemination Capacities of a BlaNDM-5-Harboring IncX-3 Plasmid in Escherichia coli Isolated from an Urban River in Montpellier, France

To investigate the capacities of persistence and dissemination of bla_NDM-5 within Escherichia coli and in aquatic environment, we characterized E. coli (sequence type 636) strains B26 and B28 isolated one month apart from the same urban river in Montpellier, France. The two isolates carried a pTsB26 plasmid, which sized 45,495 Kb, harbored bla_NDM-5 gene and belonged to IncX-3 incompatibility group. pTsB26 was conjugative in vitro at high frequency, it was highly stable after 400 generations and it exerted no fitness cost on its host. bla_NDM-5harboring plasmids are widely dispersed in E. coli all around the world, with no lineage specialization. The genomic comparison between B26 and B28 stated that the two isolates probably originated from the same clone, suggesting the persistence of pTsB26 in an E. coli host in aquatic environment.

Specificities and Commonalities of Carbapenemase-Producing Escherichia coli Isolated in France from 2012 to 2015

Carbapenemase-producing Escherichia coli (CP-Ec) represents a major public health threat with a risk of dissemination in the community as has occurred for lineages producing extended-spectrum β-lactamases. To characterize the extent of CP-Ec spread in France, isolates from screening and infection samples received at the French National Reference Center (F-NRC) laboratory for carbapenemase-producing Enterobacterales were investigated. A total of 691 CP-Ec isolates collected between 2012 and 2015 and 22 isolates collected before 2012 were fully sequenced. Analysis of their genome sequences revealed some disseminating multidrug-resistant (MDR) lineages frequently acquiring diverse carbapenemase genes mainly belonging to clonal complex 23 (CC23) (sequence type 410 [ST410]) and CC10 (ST10 and ST167) and sporadic isolates, including rare ST131 isolates (n = 17). However, the most represented sequence type (ST) was ST38 (n = 92) with four disseminated lineages carrying bla_OXA-48-like genes inserted in the chromosome. Globally, the most frequent carbapenemase gene (n = 457) was bla_OXA-48. It was also less frequently associated with MDR isolates being the only resistance gene in 119 isolates. Thus, outside the ST38 clades, its acquisition was frequently sporadic with no sign of dissemination, reflecting the circulation of the IncL plasmid pOXA-48 in France and its high frequency of conjugation. In contrast, bla_OXA-181 and bla_NDM genes were often associated with the evolution of MDR E. coli lineages characterized by mutations in ftsI and ompC.

IMPORTANCE Carbapenemase-producing Escherichia coli (CP-Ec) might be difficult to detect, as MICs can be very low. However, their absolute number and their proportion among carbapenem-resistant Enterobacterales have been increasing, as reported by WHO and national surveillance programs. This suggests a still largely uncharacterized community spread of these isolates. Here, we have characterized the diversity and evolution of CP-Ec isolated in France before 2016. We show that carbapenemase genes are associated with a wide variety of E. coli genomic backgrounds and a small number of dominant phylogenetic lineages. In a significant proportion of CP-Ec, the most frequent carbapenemase gene bla_OXA-48, was detected in isolates lacking any other resistance gene, reflecting the dissemination of pOXA-48 plasmids, likely in the absence of any antibiotic pressure. In contrast, carbapenemase gene transfer may also occur in multidrug-resistant E. coli, ultimately giving rise to at-risk lineages encoding carbapenemases with a high potential of dissemination.

metaSNV v2: detection of SNVs and subspecies in prokaryotic metagenomes

SUMMARY

Taxonomic analysis of microbial communities is well supported at the level of species and strains. However, species can contain significant phenotypic diversity and strains are rarely widely shared across global populations. Stratifying the diversity between species and strains can identify 'subspecies', which are a useful intermediary. High-throughput identification and profiling of subspecies is not yet supported in the microbiome field. Here, we use an operational definition of subspecies based on single nucleotide variant (SNV) patterns within species to identify and profile subspecies in metagenomes, along with their distinctive SNVs and genes. We incorporate this method into metaSNV v2, which extends existing SNV-calling software to support further SNV interpretation for population genetics. These new features support microbiome analyses to link SNV profiles with host phenotype or environment and niche-specificity. We demonstrate subspecies identification in marine and fecal metagenomes. In the latter, we analyze 70 species in 7524 adult and infant subjects, supporting a common subspecies population structure in the human gut microbiome and illustrating some limits in subspecies calling.

AVAILABILITY AND IMPLEMENTATION

Source code, documentation, tutorials and test data are available at https://github.com/metasnv-tool/metaSNV and https://metasnv.embl.de.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Extended-Spectrum Beta-Lactamase Producing-Escherichia coli Isolated From Irrigation Waters and Produce in Ecuador

In cities across the globe, the majority of wastewater – that includes drug resistant and pathogenic bacteria among other contaminants – is released into streams untreated. This water is often subsequently used for irrigation of pastures and produce. This use of wastewater-contaminated streams allows antibiotic-resistant bacteria to potentially cycle back to humans through agricultural products. In this study, we investigated the prevalence of extended-spectrum β-lactamase (ESBL)-producing Escherichia coli isolated from produce and irrigation water across 17 provinces of Ecuador. A total of 117 vegetable samples, 119 fruit samples, and 38 irrigation water samples were analyzed. Results showed that 11% of the samples were positive for E. coli including 11 irrigation water samples (29%), and samples of 13 vegetables (11%), and 11 fruits (9%). Among the 165 E. coli isolates cultured, 96 (58%) had the ESBL phenotype, and 58% of ESBL producing E. coli came from irrigation water samples, 11% from vegetables, and 30% from fruits. The bla_CTX–M–55, bla_{CTX–M 65}, and bla_{CTX–M 15} genes were the most frequently found gene associated with the ESBL phenotype and coincided with the bla_CTX–M alleles associated with human infections in Ecuador. Three isolates had the mcr-1 gene which is responsible for colistin resistance. This report provides evidence of the potential role of irrigation water in the growing antimicrobial resistance crisis in Ecuador.

Emergence and dissemination of antimicrobial resistance in Escherichia coli causing bloodstream infections in Norway in 2002-17: a nationwide, longitudinal, microbial population genomic study

BACKGROUND

The clonal diversity underpinning trends in multidrug resistant Escherichia coli causing bloodstream infections remains uncertain. We aimed to determine the contribution of individual clones to resistance over time, using large-scale genomics-based molecular epidemiology.

METHODS

This was a longitudinal, E coli population, genomic, cohort study that sampled isolates from 22 512 E coli bloodstream infections included in the Norwegian surveillance programme on resistant microbes (NORM) from 2002 to 2017. 15 of 22 laboratories were able to share their isolates, and the first 22·5% of isolates from each year were requested. We used whole genome sequencing to infer the population structure (PopPUNK), and we investigated the clade composition of the dominant multidrug resistant clonal complex (CC)131 using genetic markers previously reported for sequence type (ST)131, effective population size (BEAST), and presence of determinants of antimicrobial resistance (ARIBA, PointFinder, and ResFinder databases) over time. We compared these features between the 2002-10 and 2011-17 time periods. We also compared our results with those of a longitudinal study from the UK done between 2001 and 2011.

FINDINGS

Of the 3500 isolates requested from the participating laboratories, 3397 (97·1%) were received, of which 3254 (95·8%) were successfully sequenced and included in the analysis. A significant increase in the number of multidrug resistant CC131 isolates from 71 (5·6%) of 1277 in 2002-10 to 207 (10·5%) of 1977 in 2011-17 (p<0·0001), was the largest clonal expansion. CC131 was the most common clone in extended-spectrum β-lactamase (ESBL)-positive isolates (75 [58·6%] of 128) and fluoroquinolone non-susceptible isolates (148 [39·2%] of 378). Within CC131, clade A increased in prevalence from 2002, whereas the global multidrug resistant clade C2 was not observed until 2007. Multiple de-novo acquisitions of both blaCTX-M ESBL-encoding genes in clades A and C1 and gain of phenotypic fluoroquinolone non-susceptibility across the clade A phylogeny were observed. We estimated that exponential increases in the effective population sizes of clades A, C1, and C2 occurred in the mid-2000s, and in clade B a decade earlier. The rate of increase in the estimated effective population size of clade A (Ne=3147) was nearly ten-times that of C2 (Ne=345), with clade A over-represented in Norwegian CC131 isolates (75 [27·0%] of 278) compared with the UK study (8 [5·4%] of 147 isolates).

INTERPRETATION

The early and sustained establishment of predominantly antimicrobial susceptible CC131 clade A isolates, relative to multidrug resistant clade C2 isolates, suggests that resistance is not necessary for clonal success. However, even in the low antibiotic use setting of Norway, resistance to important antimicrobial classes has rapidly been selected for in CC131 clade A isolates. This study shows the importance of genomic surveillance in uncovering the complex ecology underlying multidrug resistance dissemination and competition, which have implications for the design of strategies and interventions to control the spread of high-risk multidrug resistant clones.

FUNDING

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

Predominance of CTX-M-15-producing ST131 strains among ESBL-producing Escherichia coli isolated from asylum seekers in the Netherlands

Objectives

Numerous studies show increased prevalence of MDR bacteria amongst asylum seekers, but data on the molecular profiles of such strains are limited. We aimed to evaluate the molecular profiles of ESBL-producing Escherichia coli (ESBL-E. coli) strains isolated from asylum seekers and investigate their phylogenetic relatedness.

Methods

WGS data of ESBL-E. coli isolates from asylum seekers, retrieved from 1 January to 31 December 2016, were analysed to assess MLST STs, fim types, phylogroups and resistance genes. Fifty-two ESBL-E. coli isolates from the Dutch–German border region were used for genome comparison purposes as a control group.

Results

Among 112 ESBL-E. coli isolates from asylum seekers, originating mostly from Syria (n = 40) and Iraq (n = 15), the majority belonged to ST131 (21.4%) and ST10 (17.0%). The predominant gene for β-lactam resistance was bla_CTX-M-15 (67.9%), followed by the often co-detected bla_TEM-1B (39.3%). No mcr or carbapenemase genes were detected. The majority of the strains belonged to phylogroups B2 (38.4%) and A (32.1%), carrying fimH27 (25%) and fimH30 (19.6%). A core genome MLST minimum spanning tree did not reveal clusters containing strains from the asylum seekers and the control group. Five clusters were formed within the asylum seeker group, by strains isolated from people originating from different countries.

Conclusions

The most frequently isolated clones in this study were isolated on a regular basis within the Dutch population before the increase in the asylum seeker population. No mcr- or carbapenemase-producing clones were detected among the asylum seeker population. Minor clustering was observed amongst the asylum seeker strains.

Diversity and Adaptations of Escherichia coli Strains: Exploring the Intestinal Community in Crohn's Disease Patients and Healthy Individuals

Crohn's disease (CD) is characterized by a chronic, progressive inflammation across the gastrointestinal tract with a series of exacerbations and remissions. A significant factor in the CD pathogenesis is an imbalance in gut microbiota composition, particularly the prevalence of Escherichia coli. In the present study, the genomes of sixty-three E. coli strains from the gut of patients with CD and healthy subjects were sequenced. In addition, eighteen E. coli-like metagenome-assembled genomes (MAGs) were reconstructed from the shotgun-metagenome sequencing data of fecal samples. The comparative analysis revealed the similarity of E. coli genomes regardless of the origin of the strain. The strains exhibited similar genetic patterns of virulence, antibiotic resistance, and bacteriocin-producing systems. The study showed antagonistic activity of E. coli strains and the metabolic features needed for their successful competition in the human gut environment. These observations suggest complex bacterial interactions within the gut which may affect the host and cause intestinal damage.

Comparative Genomics of Emerging Lineages and Mobile Resistomes of Contemporary Broiler Strains of Salmonella Infantis and E. coli

Introduction

Commensal and pathogenic strains of multidrug-resistant (MDR) Escherichia coli and non-typhoid strains of Salmonella represent a growing foodborne threat from foods of poultry origin. MDR strains of Salmonella Infantis and E. coli are frequently isolated from broiler chicks and the simultaneous presence of these two enteric bacterial species would potentially allow the exchange of mobile resistance determinants.

Objectives

In order to understand possible genomic relations and to obtain a first insight into the potential interplay of resistance genes between enteric bacteria, we compared genomic diversity and mobile resistomes of S. Infantis and E. coli from broiler sources.

Results

The core genome MLST analysis of 56 S. Infantis and 90 E. coli contemporary strains revealed a high genomic heterogeneity of broiler E. coli. It also allowed the first insight into the genomic diversity of the MDR clone B2 of S. Infantis, which is endemic in Hungary. We also identified new MDR lineages for S. Infantis (ST7081 and ST7082) and for E. coli (ST8702 and ST10088). Comparative analysis of antibiotic resistance genes and plasmid types revealed a relatively narrow interface between the mobile resistomes of E. coli and S. Infantis. The mobile resistance genes tet(A), aadA1, and sul1 were identified at an overall high prevalence in both species. This gene association is characteristic to the plasmid pSI54/04 of the epidemic clone B2 of S. Infantis. Simultaneous presence of these genes and of IncI plasmids of the same subtype in cohabitant caecal strains of E. coli and S. Infantis suggests an important role of these plasmid families in a possible interplay of resistance genes between S. Infantis and E. coli in broilers.

Conclusion

This is the first comparative genomic analysis of contemporary broiler strains of S. Infantis and E. coli. The diversity of mobile resistomes suggests that commensal E. coli could be potential reservoirs of resistance for S. Infantis, but so far only a few plasmid types and mobile resistance genes could be considered as potentially exchangeable between these two species. Among these, IncI1 plasmids could make the greatest contribution to the microevolution and genetic interaction between E. coli and S. Infantis.

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

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

A comprehensive and high-quality collection of Escherichia coli genomes and their genes

Escherichia coli is a highly diverse organism that includes a range of commensal and pathogenic variants found across a range of niches and worldwide. In addition to causing severe intestinal and extraintestinal disease, E. coli is considered a priority pathogen due to high levels of observed drug resistance. The diversity in the E. coli population is driven by high genome plasticity and a very large gene pool. All these have made E. coli one of the most well-studied organisms, as well as a commonly used laboratory strain. Today, there are thousands of sequenced E. coli genomes stored in public databases. While data is widely available, accessing the information in order to perform analyses can still be a challenge. Collecting relevant available data requires accessing different sources, where data may be stored in a range of formats, and often requires further manipulation and processing to apply various analyses and extract useful information. In this study, we collated and intensely curated a collection of over 10 000 E. coli and Shigella genomes to provide a single, uniform, high-quality dataset. Shigella were included as they are considered specialized pathovars of E. coli. We provide these data in a number of easily accessible formats that can be used as the foundation for future studies addressing the biological differences between E. coli lineages and the distribution and flow of genes in the E. coli population at a high resolution. The analysis we present emphasizes our lack of understanding of the true diversity of the E. coli species, and the biased nature of our current understanding of the genetic diversity of such a key pathogen.

Genomic Survey of E. coli From the Bladders of Women With and Without Lower Urinary Tract Symptoms

Urinary tract infections (UTIs) are one of the most common human bacterial infections. While UTIs are commonly associated with colonization by Escherichia coli, members of this species also have been found within the bladder of individuals with no lower urinary tract symptoms (no LUTS), also known as asymptomatic bacteriuria. Prior studies have found that both uropathogenic E. coli (UPEC) strains and E. coli isolates that are not associated with UTIs encode for virulence factors. Thus, the reason(s) why E. coli sometimes causes UTI-like symptoms remain(s) elusive. In this study, the genomes of 66 E. coli isolates from adult female bladders were sequenced. These isolates were collected from four cohorts, including women: (1) without lower urinary tract symptoms, (2) overactive bladder symptoms, (3) urgency urinary incontinence, and (4) a clinical diagnosis of UTI. Comparative genomic analyses were conducted, including core and accessory genome analyses, virulence and motility gene analyses, and antibiotic resistance prediction and testing. We found that the genomic content of these 66 E. coli isolates does not correspond with the participant's symptom status. We thus looked beyond the E. coli genomes to the composition of the entire urobiome and found that the presence of E. coli alone was not sufficient to distinguish between the urobiomes of individuals with UTI and those with no LUTS. Because E. coli presence, abundance, and genomic content appear to be weak predictors of UTI status, we hypothesize that UTI symptoms associated with detection of E. coli are more likely the result of urobiome composition.