Characterization of the cryptic Escherichia lineages: rapid identification and prevalence

Therapeutic Phage Candidates for Targeting Prevalent Sequence Types of Carbapenem-Resistant Escherichia coli

Carbapenem-resistant Escherichia coli (CREC) is a global threat to public health; therefore, alternative treatment options are urgently needed. Bacteriophages have emerged as promising candidates for combating CREC infections. This study aimed to investigate the genetic basis of phage sensitivity in CREC by evaluating carbapenem resistance among multidrug-resistant (MDR) E. coli isolated in Daegu, South Korea and analyzing their sequence types (STs) with phage susceptibility spectra. Among the 60 MDR E. coli isolates, 80.4% were identified as CREC, with 77.0% demonstrating resistance to imipenem and 66.6% to meropenem. Moreover, 70 lytic E. coli bacteriophages were isolated from hospital sewage water and evaluated against those 60 E. coli isolates. The phages exhibited lytic activity of 33%-60%, with average titers ranging from 5.6 × 1012 to 2.4 × 1013 PFU/mL (Plaque-Forming Unit). Furthermore, multilocus sequence typing (MLST) analysis of the bacterial isolates revealed 14 distinct STs, mostly belonging to ST131, ST410, and ST648. Notably, the phage susceptibility spectra of ST73, ST13003, ST648, ST2311, ST167, ST405, ST607, ST7962, and ST131 were significantly different. Thus, the isolated phages can effectively lyse CREC isolates, particularly those with clinically dominant STs. Conversely, ST410 exhibited a 14.2%-87.14% susceptibility spectrum, whereas ST1139, ST1487, ST10, and ST206 did not lyse, suggesting the presence of more resistant STs. Future studies are warranted to identify the reasons behind this resistance and address it. Ultimately, this study will aid in developing focused treatments to address these pressing global health issues.

Mobile genetic elements contributing to horizontal gene transfer of blaNDM among Escherichia coli in the community setting

OBJECTIVE

To investigate the distribution of carbapenem-resistant Enterobacterales (CRE) in the community and to describe the genomic characteristics.

METHODS

CRE screened from fecal samples in healthy people at the health examination center of a tertiary hospital in China underwent Whole genome sequencing (WGS) to analyze genotypic characteristics of CRE. The flanking DNA sequence of blaNDM-5 and mcr1.1 genes were analyzed by Gcluster software.

RESULTS

A total of 7187 fecal samples were screened, and CRE carriage was detected in 0.4 % of the sampled population. In total, 30 Escherichia coli, one Citrobacter freundii and one Klebsiella aerogene were screened. The 30 carbapenem-resistant Escherichia coli (CREC) isolates displayed slight resistance to amikacin (13.3 %) and aztreonam (20.0 %). All the CRE isolates contained blaNDM, and blaNDM-5 (84.4 %) was the most common one. B1 (n = 11) and A (n = 7) were predominant phylogroups. Furthermore, 34 distinct plasmid replicons, 67 different VFs, 22 distinct STs, 17 different FimH types, 26 O:H serotypes as well as 74 MGEs including 61 insertion sequences and 13 transposons were identified. The flanking DNA sequence analysis of blaNDM-5 and mcr1.1 genes indicates the key role of horizontal transfer of blaNDM-5 in the CRE development evidenced by diverse STs and phylogenetic tree.

CONCLUSION

E. coli was the most predominant CRE isolates in community setting, and blaNDM (blaNDM-5) was the main CHβL encoding genes. The high prevalence of ARGs was associated with high resistance to commonly used antimicrobials. Besides, the genetic diversity of these isolates suggested the key role of blaNDM horizontal transfer in the CRE development. Thus, active screening of blaNDM in communities is particularly important for the prevention and control of CRE.

Population structure and pathogen interaction of Escherichia coli in freshwater: Implications of land-use for water quality and public health in Aotearoa New Zealand

Freshwater samples (n = 199) were obtained from 41 sites with contrasting land-uses (avian, low impact, dairy, urban, sheep and beef, and mixed sheep, beef and dairy) and the E. coli phylotype of 3980 isolates (20 per water sample enrichment) was determined. Eight phylotypes were identified with B1 (48.04%), B2 (14.87%) and A (14.79%) the most abundant. Escherichia marmotae (n = 22), and Escherichia ruysiae (n = 1), were rare (0.68%) suggesting that these environmental strains are unlikely to confound water quality assessments. Phylotypes A and B1 were overrepresented in dairy and urban sites (p < 0.0001), whilst B2 were overrepresented in low impact sites (p < 0.0001). Pathogens ((Salmonella, Campylobacter, Cryptosporidium or Giardia) and the presence of diarrhoeagenic E. coli-associated genes (stx and eae) were detected in 89.9% (179/199) samples, including 80.5% (33/41) of samples with putative non-recent faecal inputs. Quantitative PCR to detect microbial source tracking targets from human, ruminant and avian contamination were concordant with land-use type and E. coli phylotype abundance. This study demonstrated that a potential recreational health risk remains where pathogens occurred in water samples with low E. coli concentration, potential non-recent faecal sources, low impact sites and where human, ruminant and avian faecal sources were absent.

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.

Development of a multiplex droplet digital PCR assay for simultaneous detection and quantification of Escherichia coli, E. marmotae, and E. ruysiae in water samples

Escherichia coli are widely used by water quality managers as Fecal Indicator Bacteria, but current quantification methods do not differentiate them from benign, environmental Escherichia species such as E. marmotae (formerly named cryptic clade V) or E. ruysiae (cryptic clades III and IV). Reliable and specific techniques for their identification are required to avoid confounding microbial water quality assessments. To address this, a multiplex droplet digital PCR (ddPCR) assay targeting lipB (E. coli and E. ruysiae) and bglC (E. marmotae) was designed. The ddPCR performance was assessed using in silico analysis; genomic DNA from 40 local, international, and reference strains of target and non-target coliforms; and spiked water samples in a range relevant to water quality managers (1 to 1000 cells/100 mL). Results were compared to an analogous quantitative PCR (qPCR) and the Colilert method. Both PCR assays showed excellent sensitivity with a limit of detection of 0.05 pg/μL and 0.005 pg/μl for ddPCR and qPCR respectively, and of quantification of 0.5 pg/μL of genomic DNA. The ddPCR allowed differentiation and quantification of three Escherichia species per run by amplitude multiplexing and showed a high concordance with concentrations measured by Colilert once proportional bias was accounted for. In silico specificity testing underlined the possibility to further detect and distinguish Escherichia cryptic clade VI. Finally, the applicability of the ddPCR was successfully tested on environmental water samples where E. marmotae and E. ruysiae potentially confound E. coli counts based on the Most Probable Number method, highlighting the utility of this novel ddPCR as an efficient and rapid discriminatory test to improve water quality assessments.

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

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

High-resolution genomic analysis to investigate the impact of the invasive brushtail possum (Trichosurus vulpecula) and other wildlife on microbial water quality assessments

Escherichia coli are routine indicators of fecal contamination in water quality assessments. Contrary to livestock and human activities, brushtail possums (Trichosurus vulpecula), common invasive marsupials in Aotearoa/New Zealand, have not been thoroughly studied as a source of fecal contamination in freshwater. To investigate their potential role, Escherichia spp. isolates (n = 420) were recovered from possum gut contents and feces and were compared to those from water, soil, sediment, and periphyton samples, and from birds and other introduced mammals collected within the Mākirikiri Reserve, Dannevirke. Isolates were characterized using E. coli-specific real-time PCR targeting the uidA gene, Sanger sequencing of a partial gnd PCR product to generate a gnd sequence type (gST), and for 101 isolates, whole genome sequencing. Escherichia populations from 106 animal and environmental sample enrichments were analyzed using gnd metabarcoding. The alpha diversity of Escherichia gSTs was significantly lower in possums and animals compared with aquatic environmental samples, and some gSTs were shared between sample types, e.g., gST535 (in 85% of samples) and gST258 (71%). Forty percent of isolates gnd-typed and 75% of reads obtained by metabarcoding had gSTs shared between possums, other animals, and the environment. Core-genome single nucleotide polymorphism (SNP) analysis showed limited variation between several animal and environmental isolates (<10 SNPs). Our data show at an unprecedented scale that Escherichia clones are shared between possums, other wildlife, water, and the wider environment. These findings support the potential role of possums as contributors to fecal contamination in Aotearoa/New Zealand freshwater. Our study deepens the current knowledge of Escherichia populations in under-sampled wildlife. It presents a successful application of high-resolution genomic methods for fecal source tracking, thereby broadening the analytical toolbox available to water quality managers. Phylogenetic analysis of isolates and profiling of Escherichia populations provided useful information on the source(s) of fecal contamination and suggest that comprehensive invasive species management strategies may assist in restoring not only ecosystem health but also water health where microbial water quality is compromised.

Occurrence and temporal distribution of extended-spectrum β-lactamase-producing Escherichia coli in clams from the Central Adriatic, Italy

The spread of extended-spectrum β-lactamase (ESBL)-producing Escherichia coli is a major public health issue. Bivalves are filter-feeder animals capable of bioaccumulating the microorganisms present in water. This physiological characteristic makes them both good indicators of environmental contamination and possible carriers of pathogenic bacteria, including those resistant to antimicrobials. The aim of this study was to investigate the occurrence of ESBL-producing E. coli in clams (n = 308) collected from harvesting areas of the Central Adriatic Sea between 2018 and 2019. ESBL- /class C β-lactamase (AmpC)- producing E. coli and Escherichia spp. were isolated by streaking over the surface of MacConkey agar plates supplemented with cefotaxime enriched broths of the initial shellfish suspension. E. coli and Escherichia spp. resistant to cefotaxime were screened for ESBL production by using the double disk synergy test. Susceptibility to different antimicrobials and confirmation of ESBL-production were determined by the minimum inhibitory concentration (MIC) test. Isolates were further characterized by whole genome sequencing (WGS) and bioinformatic analysis of genomes with different tools. Overall, ESBL-producing E. coli were isolated from 3% of the samples. Of 13 ESBL- and ESBL-/AmpC-producing Escherichia spp. (n = 11 E. coli, n = 1 E. marmotae, n = 1 E. ruysiae) isolates, 13 were resistant to ampicillin and cefotaxime, 9 to sulfamethoxazole, 6 to tetracycline and nalidixic acid, 4 to trimethoprim, and 3 to ceftazidime, cefoxitin, ciprofloxacin, and chloramphenicol. Moreover, the majority (8/11) of the ESBL-producing E. coli isolates were multidrug-resistant. WGS showed that the isolates predominantly carried the blaCTX-M-15 gene (3/11) and blaCTX-M-14 and blaCTX-M-1 (2/11 each). The AmpC β-lactamase CMY-2 was found in two isolates. Phylogroup A was the most prevalent (5/11), followed by phylogroups D (4/11), F (1/11), and B2 (1/11). Ten different sequence types (STs) were identified. Occurrence at sampling sites ranged between 0 and 27%. To identify associations between the occurrence of ESBL-producing E. coli and E. coli levels, samples were divided into two groups, with E. coli at >230 MPN/100 g and E. coli at ≤230 MPN/100 g. ESBL-producing E. coli isolates were significantly more commonly recovered in samples with higher E. coli levels (14%) than in those with lower levels of E. coli (2%). Moreover, the majority (3/4) of the potentially pathogenic strains were isolated in samples with higher E. coli levels. These findings provided evidence for the bacterial indicator of fecal contamination, E. coli, as an index organism for ESBL-producing E. coli isolates in bivalves.

Effect of Type of Mulch on Microbial Food Safety Risk on Cucumbers Irrigated with Contaminated Water

Mulches are used to block light and retain soil moisture which may affect the survival of bacterial pathogens on soil. This study examined the effectiveness of different types of mulches to minimize microbial risk from contaminated water used for irrigation of cucumbers. A production bed of 120 ft2 with 18 beds (30 ft long) covered with five different types of mulch (paper, paper with fertilizer incorporated (PF), maize-based mulch, biodegradable plastic covering, and conventional plastic) including three beds with no cover was planted with Dasher 2 Variety cucumber. Soil samples from each bed were collected for the first five weeks to examine natural E. coli and coliforms. Well water contaminated with or without nalidixic acid-resistant mutant of E. coli (8 Log CFU/mL) was used for drip irrigation for 7 days before harvesting. Prior to irrigation with contaminated water, naturally present E. coli and coliform in the soil samples with or without mulch were in the range of 3.45-3.78 Log CFU/g and 4.18-5.31 Log CFU/g, respectively. E. coli levels on cucumbers harvested from mulched beds and irrigated with contaminated irrigation water had significantly higher (P < 0.05) levels of E. coli as compared with samples from similar beds irrigated with noncontaminated water. However, Cucumber, harvested within each irrigation water quality were not significantly different regardless of the type of mulch with E. coli levels from 1.72 to 3.30 Log CFU/cm2 (contaminated water) and 0.28-1.86 Log CFU/cm2 (noncontaminated water). A significant die-off of inoculated E. coli was observed on cucumber within 3 days (>1.17 Log CFU/cm2) and >1.38 Log CFU/cm2 after 4 days. Beds with maize mulch were effective on minimizing E. coli contamination on cucumber from contaminated irrigation water.

Molecular characterization of highly prevalent Escherichia coli and Escherichia marmotae resistant to extended-spectrum cephalosporins in European starlings (Sturnus vulgaris) in Tunisia

European starlings are widespread migratory birds that have already been described as carrying bacteria resistant to extended-spectrum cephalosporins (ESC-R). These birds are well known in Tunisia because they spend the wintertime in this country and are hunted for human consumption. The goal of our study was to estimate the proportion of ESC-R in these birds and to characterize the collected isolates using whole-genome sequencing. Results showed that 21.5% (42/200) of the birds carried either an extended-spectrum beta-lactamase (ESBL) or an acquired AmpC gene. Diverse bla CTX-M genes were responsible for the ESBL phenotype, bla CTX-M-14 being the most prevalent, while only bla CMY-2 and one bla CMY-62 were found in AmpC-positive isolates. Likewise, different genetic determinants carried these resistance genes, including IncHI2, and IncF plasmids for bla CTX-M genes and IncI1 plasmids for bla CMY-2 genes. Three chromosomally encoded bla CTX-M-15 genes were also identified. Surprisingly, species identification revealed a large proportion (32.7%) of Escherichia marmotae isolates. This species is phenotypically indistinguishable from Escherichia coli and has obviously the same capacity to acquire ESC-R genes. Our data also strongly suggest that at least the IncHI2/pST3 plasmid can spread equally between E. coli and E. marmotae. Given the potential transmission routes between humans and animals, either by direct contact with dejections or through meat preparation, it is important to closely monitor antimicrobial resistance in European starlings in Tunisia and to set up further studies to identify the sources of contamination of these birds. IMPORTANCE The One Health concept highlighted knowledge gaps in the understanding of the transmission routes of resistant bacteria. A major interest was shown in wild migratory birds since they might spread resistant bacteria over long distances. Our study brings further evidence that wild birds, even though they are not directly submitted to antibiotic treatments, can be heavily contaminated by resistant bacteria. Our results identified numerous combinations of resistance genes, genetic supports, and bacterial clones that can spread vertically or horizontally and maintain a high level of resistance in the bird population. Some of these determinants are widespread in humans or animals (IncHI2/pST3 plasmids and pandemic clones), while some others are less frequent (atypical IncI1 plasmid and minor clones). Consequently, it is essential to be aware of the risks of transmission and to take all necessary measures to prevent the proportions of resistant isolates from increasing uncontrollably.

Escherichia ruysiae May Serve as a Reservoir of Antibiotic Resistance Genes across Multiple Settings and Regions

Gut colonization with multidrug-resistant Enterobacterales (MDR-Ent) has reached worrisome levels worldwide. In this context, Escherichia ruysiae is a recently described species mostly found in animals. However, its spread and impact on humans is poorly understood. A stool sample from a healthy individual living in India was screened for the presence of MDR-Ent using culture-based methods. Colonies were routinely identified using matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) and phenotypically characterized by broth microdilution. Illumina and Nanopore whole-genome sequencing (WGS) platforms were implemented to generate a complete assembly. E. ruysiae genomes deposited in international databases were used for a core genome phylogenetic analysis. An extended-spectrum β-lactamase (ESBL)-producing E. coli strain (S1-IND-07-A) was isolated from the stool. WGS confirmed that S1-IND-07-A was indeed E. ruysiae, belonged to sequence type 5792 (ST5792), core genome (cg) ST89059, serotype O13/O129-:H56-like, clade IV phylogroup, and possessed five virulence factors. A copy of blaCTX-M-15 and five other antimicrobial resistance genes (ARGs) were detected in a conjugative IncB/O/K/Z plasmid. A database search identified 70 further E. ruysiae strains from 16 countries (44, 15, and 11 strains isolated from animals, the environment, and humans, respectively). The core genome phylogeny revealed five major STs: ST6467, ST8084, ST2371, ST9287, and ST5792. Three out of the seventy strains possessed important ARGs: OTP1704 (blaCTX-M-14; ST6467), SN1013-18 (blaCTX-M-15; ST5792), and CE1758 (blaCMY-2; ST7531). These strains were of human, environmental, and wild animal origin, respectively. E. ruysiae may acquire clinically important ARGs and transmit them to other species. Due to its zoonotic potential, further efforts are needed to improve routine detection and surveillance across One Health settings. IMPORTANCE Escherichia ruysiae is a recently described species of the cryptic clades III and IV of the genus Escherichia and is commonly found in animals and the environment. This work highlights the zoonotic potential of E. ruysiae, as it has been shown to colonize the human intestinal tract. Importantly, E. ruysiae may be associated with conjugative plasmids carrying clinically relevant antibiotic resistance genes. Therefore, it is important to closely monitor this species. Overall, this study highlights the need for improved identification of Escherichia species and continued surveillance of zoonotic pathogens in One Health settings.

Under-Appreciated Phylogroup Diversity of Escherichia coli within and between Animals at the Urban-Wildland Interface

Wild animals have been implicated as reservoirs and even "melting pots" of pathogenic and antimicrobial-resistant bacteria of concern to human health. Though Escherichia coli is common among vertebrate guts and plays a role in the propagation of such genetic information, few studies have explored its diversity beyond humans nor the ecological factors that influence its diversity and distribution in wild animals. We characterized an average of 20 E. coli isolates per scat sample (n = 84) from a community of 14 wild and 3 domestic species. The phylogeny of E. coli comprises 8 phylogroups that are differentially associated with pathogenicity and antibiotic resistance, and we uncovered all of them in one small biological preserve surrounded by intense human activity. Challenging previous assumptions that a single isolate is representative of within-host phylogroup diversity, 57% of individual animals sampled carried multiple phylogroups simultaneously. Host species' phylogroup richness saturated at different levels across species and encapsulated vast within-sample and within-species variation, indicating that distribution patterns are influenced both by isolation source and laboratory sampling depth. Using ecological methods that ensure statistical relevance, we identify trends in phylogroup prevalence associated with host and environmental factors. The vast genetic diversity and broad distribution of E. coli in wildlife populations has implications for biodiversity conservation, agriculture, and public health, as well as for gauging unknown risks at the urban-wildland interface. We propose critical directions for future studies of the "wild side" of E. coli that will expand our understanding of its ecology and evolution beyond the human environment. IMPORTANCE To our knowledge, neither the phylogroup diversity of E. coli within individual wild animals nor that within an interacting multispecies community have previously been assessed. In doing so, we uncovered the globally known phylogroup diversity from an animal community on a preserve imbedded in a human-dominated landscape. We revealed that the phylogroup composition in domestic animals differed greatly from that in their wild counterparts, implying potential human impacts on the domestic animal gut. Significantly, many wild individuals hosted multiple phylogroups simultaneously, indicating the potential for strain-mixing and zoonotic spillback, especially as human encroachment into wildlands increases in the Anthropocene. We reason that due to extensive anthropogenic environmental contamination, wildlife is increasingly exposed to our waste, including E. coli and antibiotics. The gaps in the ecological and evolutionary understanding of E. coli thus necessitate a significant uptick in research to better understand human impacts on wildlife and the risk for zoonotic pathogen emergence.

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.

Genomic diversity of non-diarrheagenic fecal Escherichia coli from children in sub-Saharan Africa and south Asia and their relatedness to diarrheagenic E. coli

Escherichia coli is a frequent member of the healthy human gastrointestinal microbiota, as well as an important human pathogen. Previous studies have focused on the genomic diversity of the pathogenic E. coli and much remains unknown about the non-diarrheagenic E. coli residing in the human gut, particularly among young children in low and middle income countries. Also, gaining additional insight into non-diarrheagenic E. coli is important for understanding gut health as non-diarrheagenic E. coli can prevent infection by diarrheagenic bacteria. In this study we examine the genomic diversity of non-diarrheagenic fecal E. coli from male and female children with or without diarrhea from countries in sub-Saharan Africa and south Asia as part of the Global Enteric Multicenter Study (GEMS). We find that these E. coli exhibit considerable genetic diversity as they were identified in all E. coli phylogroups and an Escherichia cryptic clade. Although these fecal E. coli lack the characteristic virulence factors of diarrheagenic E. coli pathotypes, many exhibit remarkable genomic similarity to previously described diarrheagenic isolates with differences attributed to mobile elements. This raises an important question of whether these non-diarrheagenic fecal E. coli may have at one time possessed the mobile element-encoded virulence factors of diarrheagenic pathotypes or may have the potential to acquire these virulence factors.

Escherichia cryptic clade II through clade VIII: Rapid detection and prevalence in feces and surface water

Bacteria of the cryptic lineage of genus Escherichia, or Escherichia cryptic clades (cryptic clades), are phenotypically indistinguishable from Escherichia coli (E. coli) using standard biochemical tests. Except for clade I (C-I), cryptic clades were hypothetically believed to be environmental but not enteric. If so, they would hinder the interpretation of current E. coli-based water quality (fecal pollution) monitoring in the United States because environmental bacteria do not indicate the presence of harmful fecal material. This study was performed to develop a rapid method for the detection of cryptic clades and to investigate their potential impact on water quality monitoring. By whole-genome comparison, one gene, named ecc (Escherichiacryptic clades), was identified to be unique to C-II through C-VIII. An end-point polymerase chain reaction (PCR) method, eccPCR, was developed by targeting the ecc. The results of in-silico and wet tests demonstrated 100 % sensitivity and specificity of the eccPCR to detect C-II through C-VIII. Based on the EPA Method 1603, 519 presumptive E. coli isolates were obtained from the fecal samples of 13 different host species and 192 isolates from surface water samples taken at four locations in a watershed of mid-Missouri. As indicated by the eccPCR amplification, the overall prevalence of C-II through C-VIII in the presumptive E. coli isolates was estimated to be about 0.6 % in the fecal samples and about 1.6 % in the water samples. Therefore, the potential impact of cryptic clades on water quality monitoring may be limited if EPA Method 1603 is used. Furthermore, clades C-II through C-VIII in stream water samples were found repeatedly only at a single sampling site, but neither at the upstream sites nor five kilometers downstream of the site. The data do not support nor reject the environmental hypothesis about cryptic clades. Further study is needed to determine the implication of the observation.

Clonal Complexes 23, 10, 131 and 38 as Genetic Markers of the Environmental Spread of Extended-Spectrum β-Lactamase (ESBL)-Producing E. coli

In accordance with the global action plan on antimicrobial resistance adopted by the World Health Assembly in 2015, there is a need to develop surveillance programs for antimicrobial resistant bacteria. In this context, we have analyzed the clonal diversity of Extended-spectrum β-lactamase (ESBL)-producing Escherichia coli (E. coli) isolated from aquatic environments and human and food samples in Spain, with the aim of determining possible clonal complexes (CCs) that act as markers of the potential risk of transmission of these resistant bacteria. The phylogenetic groups, sequence types (STs) and CCs were determined by different Polymerase Chain Reaction (PCR) and Multilocus Sequence Typing (MLST) techniques. Phylogroup A was prevalent and was mainly present in food and water strains, while human strains were mostly associated with phylogroup B2. According to the observed prevalence in the different niches, CC23 and CC10 are proposed as markers of phylogroups A and C, related with the spread of blaCTX-M1 and blaCTX-M15 genes. Similarly, CC131 and CC38 could be associated to the dissemination of pathogenic strains (phylogroups B2 and D) carrying mainly blaCTX-M14 and blaCTX-M15 genes. Some strains isolated from wastewater treatment plants (WWTPs) showed identical profiles to those isolated from other environments, highlighting the importance that water acquires in the dissemination of bacterial resistance. In conclusion, the detection of these genetic markers in different environments could be considered as an alert in the spread of ESBL.

The microbial ecology of Escherichia coli in the vertebrate gut

Escherichia coli has a rich history as biology's 'rock star', driving advances across many fields. In the wild, E. coli resides innocuously in the gut of humans and animals but is also a versatile pathogen commonly associated with intestinal and extraintestinal infections and antimicrobial resistance-including large foodborne outbreaks such as the one that swept across Europe in 2011, killing 54 individuals and causing approximately 4000 infections and 900 cases of haemolytic uraemic syndrome. Given that most E. coli are harmless gut colonizers, an important ecological question plaguing microbiologists is what makes E. coli an occasionally devastating pathogen? To address this question requires an enhanced understanding of the ecology of the organism as a commensal. Here, we review how our knowledge of the ecology and within-host diversity of this organism in the vertebrate gut has progressed in the 137 years since E. coli was first described. We also review current approaches to the study of within-host bacterial diversity. In closing, we discuss some of the outstanding questions yet to be addressed and prospects for future research.

Virulence and phylogenetic groups of Escherichia coli cultured from raw sewage in Kuwait

In Kuwait, some untreated sewage is discharged into the sea which poses health risks. Therefore, we determined the virulence traits and the phylogenetic groups of E. coli cultured from raw sewage. Sewage was collected once every month for 12 months with culturing of a total of 140 E. coli isolates. E. coli was typed by the methods of Clermont. The five pathotypes of diarrheagenic E. coli (DEC), and extra-intestinal pathogenic E. coli (ExPEC) were detected by specific PCR assays. Four virulence genes which correlate with pathogenicity in animal models, were used for the first time for detection of ExPEC—vat (vacuolating auto-transporter toxin), fyuA (yersiniabactin receptor), chuA (heme-binding protein), and yfcV (major subunit of putative chaperon-usher fimbria). Most E. coli belonged to phylogenetic groups A (65[45%]) and B1 (34[24.3%]). Three (2.1%) isolates were DEC, while 14 (10%) isolates were ExPEC mostly in group B2 (57.1%). A relatively high prevalence of ExPEC in sewage has public health implications.

Whole-Genome Sequencing and Virulome Analysis of Escherichia coli Isolated from New Zealand Environments of Contrasting Observed Land Use

Generic Escherichia coli is commonly used as an indicator of fecal contamination to assess water quality and human health risk. Where measured E. coli exceedances occur, the presence of other pathogenic microorganisms, such as Shiga toxin-producing E. coli (STEC), is assumed, but confirmatory data are lacking. Putative E. coli isolates (n = 709) were isolated from water, sediment, soil, periphyton, and feces samples (n = 189) from five sites representing native forest and agricultural environments. Ten E. coli isolates (1.41%) were stx₂ positive, 19 (2.7%) were eae positive, and stx₁-positive isolates were absent. At the sample level, stx₂-positive E. coli (5 of 189, 2.6%) and eae-positive isolates (16 of 189, 8.5%) were rare. Using real-time PCR, these STEC-associated virulence factors were determined to be more prevalent in sample enrichments (stx₁, 23.9%; stx₂, 31.4%; eae, 53.7%) and positively correlated with generic E. coli isolate numbers (P < 0.05) determined using culture-based methods. Whole-genome sequencing (WGS) was undertaken on a subset of 238 isolates with assemblies representing seven E. coli phylogroups (A, B1, B2, C, D, E, and F), 22 Escherichia marmotae isolates, and 1 Escherichia ruysiae isolate. Virulence factors, including those from extraintestinal pathogenic E. coli, were extremely diverse in isolates from the different locations and were more common in phylogroup B2. Analysis of the virulome from WGS data permitted the identification of gene repertoires that may be involved in environmental fitness and broadly align with phylogroup. Although recovery of STEC isolates was low, our molecular data indicate that they are likely to be widely present in environmental samples containing diverse E. coli phylogroups. IMPORTANCE This study takes a systematic sampling approach to assess the public health risk of Escherichia coli recovered from freshwater sites within forest and farmland. The New Zealand landscape is dominated by livestock farming, and previous work has demonstrated that "recreational exposure to water" is a risk factor for human infection by Shiga toxin-producing Escherichia coli (STEC). Though STEC isolates were rarely isolated from water samples, STEC-associated virulence factors were identified more commonly from water sample culture enrichments and were associated with increased generic E. coli concentrations. Whole-genome sequencing data from both E. coli and newly described Escherichia spp. demonstrated the presence of virulence factors from E. coli pathotypes, including extraintestinal pathogenic E. coli. This has significance for understanding and interpreting the potential health risk from E. coli where water quality is poor and suggests a role of virulence factors in survival and persistence of E. coli and Escherichia spp.

Fecal Shedding of Multidrug Resistant Escherichia coli Isolates in Dogs Fed with Raw Meat-Based Diets in Brazil

The practice of feeding dogs raw meat-based diets (RMBDs) is growing in several countries, and the risks associated with the ingestion of pathogenic and antimicrobial-resistant Escherichia coli in dogs fed these diets are largely unknown. We characterized E. coli strains isolated from dogs fed either an RMBD or a conventional dry feed, according to the phylogroup, virulence genes, and antimicrobial susceptibility profiles of the bacteria. Two hundred and sixteen E. coli strains were isolated. Dogs fed RMBDs shed E. coli strains from the phylogroup E more frequently and were positive for the E. coli heat-stable enterotoxin 1-encoding gene. Isolates from RMBD-fed dogs were also frequently positive for multidrug-resistant E. coli isolates including extended-spectrum beta-lactamase (ESBL) producers. Whole-genome sequencing of seven ESBL-producing E. coli strains revealed that they predominantly harbored blaCTX-M-55, and two strains were also positive for the colistin-resistant gene mcr-1. These results suggest that feeding an RMBD can affect the dog’s microbiota, change the frequency of certain phylogroups, and increase the shedding of diarrheagenic E. coli. Also, feeding an RMBD seemed to be linked with the fecal shedding of multidrug-resistant E. coli, including the spread of strains harboring mobilizable colistin resistance and ESBL genes. This finding is of concern for both animal and human health.

Escherichia marmotae-a Human Pathogen Easily Misidentified as Escherichia coli

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

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

Prevalence of Escherichia coli strains in horticultural farms from Argentina: antibiotic resistance, biofilm formation, and phylogenetic affiliation

Escherichia coli is the bacteria most commonly used as an indicator of fecal contamination in agricultural environments. Moreover, E. coli is categorized as a priority pathogen due to its widespread antibiotic resistance. This study aimed to characterize E. coli strains isolated from 10 horticultural farms. Isolates were obtained from samples of vegetable crops (n = 62), the surrounding soil (n = 62), poultry litter (n = 8), and groundwater (n = 6). Phyllo-grouping assignment was performed on the total of E. coli isolates. Antibiograms and quantification of the minimal inhibitory concentration (MIC) were performed with antibiotics commonly used in humans. Biofilm formation capacity was studied by quantifying cells attached to culture tubes. Overall, 21 E. coli isolates were obtained. Three phylogenetic groups (A, B1, and C) and two Escherichia clade IV and IV-V were identified in the collection by polymerase chain reaction. Sixty-seven percent of the E. coli isolates were resistant to amoxicillin-clavulanic acid and/or ampicillin. Amoxicillin MIC values ranged from 11.9 to >190.5 µg/mL and ampicillin MIC values ranged from 3 to >190.5 µg/mL. All the E. coli isolates, resistant and non-resistant, had biofilm forming capacity. The presence of phenotypic resistance on fresh produce and environmental matrices could present significant opportunities for contamination that result in health risks for consumers. To the authors' best knowledge, this is the first environmental assessment of resistant E. coli occurrence in horticultural farms in South America.

Remarkable genomic diversity among Escherichia isolates recovered from healthy chickens

The genus Escherichia has been extensively studied and it is known to encompass a range of commensal and pathogenic bacteria that primarily inhabit the gastrointestinal tracts of warm-blooded vertebrates. However, the presence of E. coli as a model organism and potential pathogen has diverted attention away from commensal strains and other species in the genus. To investigate the diversity of Escherichia in healthy chickens, we collected fecal samples from antibiotic-free Lohmann Brown layer hens and determined the genome sequences of 100 isolates, 81 of which were indistinguishable at the HC0 level of the Hierarchical Clustering of Core Genome Multi-Locus Sequence Typing scheme. Despite initial selection on CHROMagar Orientation medium, which is considered selective for E. coli, in silico phylotyping and core genome single nucleotide polymorphism analysis revealed the presence of at least one representative of all major clades of Escherichia, except for E. albertii, Shigella, and E. coli phylogroup B2 and cryptic clade I. The most frequent phylogenomic groups were E. coli phylogroups A and B1 and E. ruysiae (clades III and IV). We compiled a collection of reference strains isolated from avian sources (predominantly chicken), representing every Escherichia phylogroup and species, and used it to confirm the phylogeny and diversity of our isolates. Overall, the isolates carried low numbers of the virulence and antibiotic resistance genes typically seen in avian pathogenic E. coli. Notably, the clades not recovered are ones that have been most strongly associated with virulence by other studies.

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.

Phylogeny and potential virulence of cryptic clade Escherichia coli species complex isolates derived from an arable field trial

Analysis of Escherichia coli taxonomy has expanded into a species-complex with the identification of divergent cryptic clades. A key question is the evolutionary trajectory of these clades and their relationship to isolates of clinical or veterinary importance. Since they have some environmental association, we screened a collection of E. coli isolated from a long-term spring barley field trial for their presence. While most isolates clustered into the enteric-clade, four of them clustered into Clade-V, and one in Clade-IV. The Clade -V isolates shared >96% intra-clade average nucleotide sequence identity but <91% with other clades. Although pan-genomics analysis confirmed their taxonomy as Clade -V (E. marmotae), retrospective phylogroup PCR did not discriminate them correctly. Differences in metabolic and adherence gene alleles occurred in the Clade -V isolates compared to E. coli sensu scricto. They also encoded the bacteriophage phage-associated cyto-lethal distending toxin (CDT) and antimicrobial resistance (AMR) genes, including an ESBL, blaOXA-453. Thus, the isolate collection encompassed a genetic diversity, and included cryptic clade isolates that encode potential virulence factors. The analysis has determined the phylogenetic relationship of cryptic clade isolates with E. coli sensu scricto and indicates a potential for horizontal transfer of virulence factors.

ECTyper: in silico Escherichia coli serotype and species prediction from raw and assembled whole-genome sequence data

Escherichia coli is a priority foodborne pathogen of public health concern and phenotypic serotyping provides critical information for surveillance and outbreak detection activities. Public health and food safety laboratories are increasingly adopting whole-genome sequencing (WGS) for characterizing pathogens, but it is imperative to maintain serotype designations in order to minimize disruptions to existing public health workflows. Multiple in silico tools have been developed for predicting serotypes from WGS data, including SRST2, SerotypeFinder and EToKi EBEis, but these tools were not designed with the specific requirements of diagnostic laboratories, which include: speciation, input data flexibility (fasta/fastq), quality control information and easily interpretable results. To address these specific requirements, we developed ECTyper (https://github.com/phac-nml/ecoli_serotyping) for performing both speciation within Escherichia and Shigella, and in silico serotype prediction. We compared the serotype prediction performance of each tool on a newly sequenced panel of 185 isolates with confirmed phenotypic serotype information. We found that all tools were highly concordant, with 92-97 % for O-antigens and 98-100 % for H-antigens, and ECTyper having the highest rate of concordance. We extended the benchmarking to a large panel of 6954 publicly available E. coli genomes to assess the performance of the tools on a more diverse dataset. On the public data, there was a considerable drop in concordance, with 75-91 % for O-antigens and 62-90 % for H-antigens, and ECTyper and SerotypeFinder being the most concordant. This study highlights that in silico predictions show high concordance with phenotypic serotyping results, but there are notable differences in tool performance. ECTyper provides highly accurate and sensitive in silico serotype predictions, in addition to speciation, and is designed to be easily incorporated into bioinformatic workflows.

Unraveling the ecological processes modulating the population structure of Escherichia coli in a highly polluted urban stream network

Escherichia coli dynamics in urban watersheds are affected by a complex balance among external inputs, niche modulation and genetic variability. To explore the ecological processes influencing E. coli spatial patterns, we analyzed its abundance and phylogenetic structure in water samples from a stream network with heterogeneous urban infrastructure and environmental conditions. Our results showed that environmental and infrastructure variables, such as macrophyte coverage, DIN and sewerage density, mostly explained E. coli abundance. Moreover, main generalist phylogroups A and B1 were found in high proportion, which, together with an observed negative relationship between E. coli abundance and phylogroup diversity, suggests that their dominance might be due to competitive exclusion. Lower frequency phylogroups were associated with sites of higher ecological disturbance, mainly involving simplified habitats, higher drainage infrastructure and septic tank density. In addition to the strong negative relationship between phylogroup diversity and dominance, the occurrence of these phylogroups would be associated with increased facilitated dispersal. Nutrients also contributed to explaining phylogroup distribution. Our study proposes the differential contribution of distinct ecological processes to the patterns of E. coli in an urban watershed, which is useful for the monitoring and management of fecal pollution.

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.

Extensive microbial diversity within the chicken gut microbiome revealed by metagenomics and culture

BACKGROUND

The chicken is the most abundant food animal in the world. However, despite its importance, the chicken gut microbiome remains largely undefined. Here, we exploit culture-independent and culture-dependent approaches to reveal extensive taxonomic diversity within this complex microbial community.

RESULTS

We performed metagenomic sequencing of fifty chicken faecal samples from two breeds and analysed these, alongside all (n = 582) relevant publicly available chicken metagenomes, to cluster over 20 million non-redundant genes and to construct over 5,500 metagenome-assembled bacterial genomes. In addition, we recovered nearly 600 bacteriophage genomes. This represents the most comprehensive view of taxonomic diversity within the chicken gut microbiome to date, encompassing hundreds of novel candidate bacterial genera and species. To provide a stable, clear and memorable nomenclature for novel species, we devised a scalable combinatorial system for the creation of hundreds of well-formed Latin binomials. We cultured and genome-sequenced bacterial isolates from chicken faeces, documenting over forty novel species, together with three species from the genus Escherichia, including the newly named species Escherichia whittamii.

CONCLUSIONS

Our metagenomic and culture-based analyses provide new insights into the bacterial, archaeal and bacteriophage components of the chicken gut microbiome. The resulting datasets expand the known diversity of the chicken gut microbiome and provide a key resource for future high-resolution taxonomic and functional studies on the chicken gut microbiome.

A review of the taxonomy, genetics, and biology of the genus Escherichia and the type species Escherichia coli

Historically, bacteriologists have relied heavily on biochemical and structural phenotypes for bacterial taxonomic classification. However, advances in comparative genomics have led to greater insights into the remarkable genetic diversity within the microbial world, and even within well-accepted species such as Escherichia coli. The extraordinary genetic diversity in E. coli recapitulates the evolutionary radiation of this species in exploiting a wide range of niches (i.e., ecotypes), including the gastrointestinal system of diverse vertebrate hosts as well as non-host natural environments (soil, natural waters, wastewater), which drives the adaptation, natural selection, and evolution of intragenotypic conspecific specialism as a strategy for survival. Over the last few years, there has been increasing evidence that many E. coli strains are very host (or niche)-specific. While biochemical and phylogenetic evidence support the classification of E. coli as a distinct species, the vast genomic (diverse pan-genome and intragenotypic variability), phenotypic (e.g., metabolic pathways), and ecotypic (host-/niche-specificity) diversity, comparable to the diversity observed in known species complexes, suggest that E. coli is better represented as a complex. Herein we review the taxonomic classification of the genus Escherichia and discuss how phenotype, genotype, and ecotype recapitulate our understanding of the biology of this remarkable bacterium.

Antimicrobial Resistance Profile and ExPEC Virulence Potential in Commensal Escherichia coli of Multiple Sources

We recently described the genetic antimicrobial resistance and virulence profile of a collection of 279 commensal E. coli of food-producing animal (FPA), pet, wildlife and human origin. Phenotypic antimicrobial resistance (AMR) and the role of commensal E. coli as reservoir of extra-intestinal pathogenic Escherichia coli (ExPEC) virulence-associated genes (VAGs) or as potential ExPEC pathogens were evaluated. The most common phenotypic resistance was to tetracycline (76/279, 27.24%), sulfamethoxazole/trimethoprim (73/279, 26.16%), streptomycin and sulfisoxazole (71/279, 25.45% both) among the overall collection. Poultry and rabbit were the sources mostly associated to AMR, with a significant resistance rate (p > 0.01) to quinolones, streptomycin, sulphonamides, tetracycline and, only for poultry, to ampicillin and chloramphenicol. Finally, rabbit was the source mostly associated to colistin resistance. Different pandemic (ST69/69*, ST95, ST131) and emerging (ST10/ST10*, ST23, ST58, ST117, ST405, ST648) ExPEC sequence types (STs) were identified among the collection, especially in poultry source. Both ST groups carried high number of ExPEC VAGs (pandemic ExPEC STs, mean = 8.92; emerging ExPEC STs, mean = 6.43) and showed phenotypic resistance to different antimicrobials (pandemic ExPEC STs, mean = 2.23; emerging ExPEC STs, mean = 2.43), suggesting their role as potential ExPEC pathogens. Variable phenotypic resistance and ExPEC VAG distribution was also observed in uncommon ExPEC lineages, suggesting commensal flora as a potential reservoir of virulence (mean = 3.80) and antimicrobial resistance (mean = 1.69) determinants.

Phylogenetic characterization of avian pathogenic Escherichia coli strains longitudinally isolated from broiler breeder flocks vaccinated with autogenous vaccine

Escherichia coli is the most common bacterial cause of infections in poultry farms. It is known for its genetic heterogenicity that complicates the protection of poultry health through immunoprophylaxis. In farms with continuous problems with colibacillosis, autogenous E. coli vaccine was implemented to the vaccination program instead of commercial vaccines. In this study, we investigated the effect of the autogenous vaccine on E. coli phylogroup diversity on 2 broiler breeder farms with 4 and 5 flocks, respectively. The first flocks on both farms were vaccinated with commercial vaccines, while application of autogenous vaccine was introduced in the second flock on both farms. In total, 113 strains were selected based on the target organs and age of chickens. Targeted organs were the peritoneum, liver, oviduct, and bone marrow, and analyzed strains were isolated from chickens older than 21 wk of age when problems with colibacillosis start emerging. The strains were phylotyped by PCR and allocated to phylogroups A, B1, B2, C, D, E, F or clades I–V. The results showed that autogenous vaccine could significantly affect the phylogroup shift of the strains. On farm A, application of the autogenous vaccine induced significantly lower prevalence (P = 0.01) of the phylogroups represented in the vaccine among the strains later isolated from the vaccinated flock, while on farm B, the results showed a decrease in the phylogenetic diversity with a dominant prevalence of group B2 despite the vaccine application. The results indicate that implementation of the autogenous vaccine can repress the majority of the strains, but also be unable to eliminate the presence of certain phylogroups, and thus lead to strain shift. Further detailed analyses of multilocus sequence typing and virulence genes will elucidate the pathogenic potential and selection of certain strains, with emphasis on B2 phylogroup.

Genetic Variation and Preliminary Indications of Divergent Niche Adaptation in Cryptic Clade II of Escherichia

The evolution, habitat, and lifestyle of the cryptic clade II of Escherichia, which were first recovered at low frequency from non-human hosts and later from external environments, were poorly understood. Here, the genomes of selected strains were analyzed for preliminary indications of ecological differentiation within their population. We adopted the delta bitscore metrics to detect functional divergence of their orthologous genes and trained a random forest classifier to differentiate the genomes according to habitats (gastrointestinal vs external environment). Model was built with inclusion of other Escherichia genomes previously demonstrated to have exhibited genomic traits of adaptation to one of the habitats. Overall, gene degradation was more prominent in the gastrointestinal strains. The trained model correctly classified the genomes, identifying a set of predictor genes that were informative of habitat association. Functional divergence in many of these genes were reflective of ecological divergence. Accuracy of the trained model was confirmed by its correct prediction of the habitats of an independent set of strains with known habitat association. In summary, the cryptic clade II of Escherichia displayed genomic signatures that are consistent with divergent adaptation to gastrointestinal and external environments.

Major role of iron uptake systems in the intrinsic extra-intestinal virulence of the genus Escherichia revealed by a genome-wide association study

The genus Escherichia is composed of several species and cryptic clades, including E. coli, which behaves as a vertebrate gut commensal, but also as an opportunistic pathogen involved in both diarrheic and extra-intestinal diseases. To characterize the genetic determinants of extra-intestinal virulence within the genus, we carried out an unbiased genome-wide association study (GWAS) on 370 commensal, pathogenic and environmental strains representative of the Escherichia genus phylogenetic diversity and including E. albertii (n = 7), E. fergusonii (n = 5), Escherichia clades (n = 32) and E. coli (n = 326), tested in a mouse model of sepsis. We found that the presence of the high-pathogenicity island (HPI), a ~35 kbp gene island encoding the yersiniabactin siderophore, is highly associated with death in mice, surpassing other associated genetic factors also related to iron uptake, such as the aerobactin and the sitABCD operons. We confirmed the association in vivo by deleting key genes of the HPI in E. coli strains in two phylogenetic backgrounds. We then searched for correlations between virulence, iron capture systems and in vitro growth in a subset of E. coli strains (N = 186) previously phenotyped across growth conditions, including antibiotics and other chemical and physical stressors. We found that virulence and iron capture systems are positively correlated with growth in the presence of numerous antibiotics, probably due to co-selection of virulence and resistance. We also found negative correlations between virulence, iron uptake systems and growth in the presence of specific antibiotics (i.e. cefsulodin and tobramycin), which hints at potential “collateral sensitivities” associated with intrinsic virulence. This study points to the major role of iron capture systems in the extra-intestinal virulence of the genus Escherichia.

Bacterial isolates belonging to the genus Escherichia can be human commensals but also opportunistic pathogens, with the ability to cause extra-intestinal infection. There is therefore the need to identify the genetic elements that favour extra-intestinal virulence, so that virulent bacterial isolates can be identified through genome analysis and potential treatment strategies be developed. To reduce the influence of host variability on virulence, we have used a mouse model of sepsis to characterize the virulence of 370 strains belonging to the genus Escherichia, for which whole genome sequences were also available. We have used a statistical approach called Genome-Wide Association Study (GWAS) to show how the presence of genes that encode for iron scavenging are significantly associated with the propensity of a bacterial isolate to cause extra-intestinal infections. Taking advantage of previously generated growth data on a subset of the strains and its correlation to virulence we generated hypothesis on the relationship between iron scavenging and growth in the presence of various antimicrobials, which could have implications for developing new treatment strategies.

Fecal indicator bacteria from environmental sources; strategies for identification to improve water quality monitoring

In tropical to temperate environments, fecal indicator bacteria (FIB), such as enterococci and Escherichia coli, can persist and potentially multiply, far removed from their natural reservoir of the animal gut. FIB isolated from environmental reservoirs such as stream sediments, beach sand and vegetation have been termed "naturalized" FIB. In addition, recent research suggests that the intestines of poikilothermic animals such as fish may be colonized by enterococci and E. coli, and therefore, these animals may contribute to FIB concentrations in the aquatic environment. Naturalized FIB that are derived from fecal inputs into the environment, and subsequently adapted to maintain their population within the non-host environment are termed "naturalized enteric FIB". In contrast, an additional theory suggests that some "naturalized" FIB diverged from enteric FIB many millions of years ago and are now normal inhabitants of the environment where they are referred to as "naturalized non-enteric FIB". In the case of the Escherichia genus, the naturalized non-enteric members are identified as E. coli during routine water quality monitoring. An over-estimation of the health risk could result when these naturalized, non-enteric FIB, (that is, not derived from avian or mammalian fecal contamination), contribute to water quality monitoring results. It has been postulated that these environmental FIB belonging to the genera Escherichia and Enterococcus can be differentiated from enteric FIB by genetic methods because they lack some of the genes required for colonization of the host intestine, and have acquired genes that aid survival in the environment. Advances in molecular tools such as next generation sequencing will aid the identification of genes peculiar or "enriched" in particular habitats to discriminate between enteric and environmental FIB. In this appraisal, we have reviewed the research studying "naturalized" FIB, and discussed the techniques for their differentiation from enteric FIB. This differentiation includes the important distinction between enteric FIB derived from fresh and non-recent fecal inputs, and those truly non-enteric environmental microbes, which are currently identified as FIB during routine water quality monitoring. The inclusion of tools for the identification of naturalized FIB (enteric or environmental) would be a valuable resource for future studies assessing water quality.

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.

Draft Genome Sequence of Escherichia marmotae E690, Isolated from Beef Cattle

We report here the draft genome sequence of an extended-spectrum β-lactamase (ESBL)-producing Escherichia species isolated from rectal feces collected from beef cattle in northern Spain. Analysis of the draft genome identified the strain as a member of the newly described species Escherichia marmotae.

We report here the draft genome sequence of an extended-spectrum β-lactamase (ESBL)-producing Escherichia species isolated from rectal feces collected from beef cattle in northern Spain. Analysis of the draft genome identified the strain as a member of the newly described species Escherichia marmotae.

Phylogenetic Characterization and Multilocus Sequence Typing of Extended-Spectrum Beta Lactamase-Producing Escherichia coli from Food-Producing Animals, Beef, and Humans in Southwest Nigeria

Multidrug-resistant extended-spectrum beta lactamase (ESBL)-producing Escherichia coli strains are emerging globally in both humans and animals. Antimicrobial susceptibility testing and ESBL screening were performed on pure cultures of 216 E. coli isolates from human and animal fecal samples as well as beef. Polymerase chain reaction was performed for the detection of resistance genes. Representative isolates of ESBL-producing E. coli were randomly selected for multilocus sequence typing and pulsed field gel electrophoresis (PFGE). Sixty of the isolates were identified as ESBL producers, and seven resistance genes were amplified in them: TEM (61.7%), blaCTX-M-15 (51.7%), AAC-6-LB (43.3%), blaCTX-M-1 (38.3%), blaCTX-M-9 (33.3%), blaCTX-M-2 (21.7%), and SHV (11.7%); they were classified into four phylogroups: A (25%), B1 (45%), B2 (20%), and D (10%). Thirty of these isolates were clustered into 10 sequence types with ST131 being mostly prevalent. Six PFGE types were discovered, each of which was shared by isolates from different subjects and had the same phylogroups and resistance gene profiles. There was a dissemination of PFGE types across various groups among humans, animals, and beef. This underlines the fact that the spread of ESBL E. coli could be from humans to animals, from animals to humans, as well as across animal species.

Plasmid Replicon Typing of Antibiotic-Resistant Escherichia coli From Clams and Marine Sediments

Unlike human isolates, environmental Escherichia coli isolates have not been thoroughly investigated for the diversity and transferability of antibiotic-resistant plasmids. In this study, antibiotic-resistant strains from marine sediment (n = 50) and clams (n = 53) were analyzed (i) for their plasmid content using a PCR-based plasmid replicon typing (PBRT) kit and (ii) for the transferability of plasmid-associated antibiotic resistance (AR) traits by mating experiments. Fifteen of the thirty replicons targeted by the PBRT kit were detected in the isolates; 8/15 were identified in both sediment and clam isolates, although at different frequencies. The most frequent replicons in sediment (74%) and in clam strains (66%) alike, were FIA, FIB, or FII, which are associated with the IncF group, followed by the I1α replicon, which was more frequent in clam (24.5%) than in sediment (10%) strains. More than 50% of the strains contained multiple replicons; although 15 were untypable, S1-PFGE analysis demonstrated that 14/15 carried no plasmids. All cryptic strains were successfully typed and were positive for IncF or IncI replicons. Antibiotic-resistant strains accounted for 63% of all isolates and were significantly (p < 0.05) more frequent in phylogroup A. Most (35%) multidrug-resistant (MDR) strains belonged to phylogroup A, too. Although 25/26 MDR strains were positive for IncF plasmids (the exception being a clam strain), the FII-FIB rep combination was predominant (63%) among the sediment isolates, whereas most clam isolates (40%) carried the FII replicon alone. In mating experiments, selected MDR strains carrying FIB, FII, and I1α replicons, used as the donors, transferred multiple ARs together with the IncF or IncI plasmids at high frequency. Since IncI plasmids are common in E. coli and Salmonella enterica isolates from poultry, our findings suggest an animal origin to the E. coli clam strains carrying IncI plasmids. They also suggest a role for IncI plasmids in the spread of ARs among environmental Enterobacteriaceae and, through the food chain, to human isolates. In conclusion, the PBRT kit proved to be a useful tool to identify plasmids carrying antibiotic-resistant genes and to shed light on the factors underpinning their diffusion.

Extended Phylogeny and Extraintestinal Virulence Potential of Commensal Escherichia coli from Piglets and Sows

Commensal Escherichia coli, naturally occurring in the intestinal tract, can be the origin of extraintestinal pathogenic E. coli (ExPEC) strains. ExPEC causes high mortality and significant economic losses in the swine industry in several countries and poses a serious threat to public health worldwide. The aim of this study was to analyze the extended phylogenetic structure and extraintestinal virulence potential in two groups of commensal E. coli isolates from post-weaning piglets and sows. The phylogenetic assignment to eight groups was determined using the revised Clermont phylogenetic typing method in quadruplex PCR. Identification of extraintestinal virulence genes (VGs) and adhesin operon genes was performed using multiplex or simplex PCR. The revised phylogenetic assignment allowed us to distinguish E. coli with significantly higher (groups C and F) or lower (group E) virulence potential in isolates from piglets. The majority of the tested VGs occurred more frequently in isolates from piglets than from sows, with statistically significant differences for seven genes: fimH, papAH, iutA, iroN, ompT, traT, and iss. Complete operons for type I and P fimbriae significantly prevailed among E. coli from piglets. This study provides insight into the extended phylogenetic structure of porcine commensal E. coli and showed that these strains, particularly from piglets, constitute a considerable reservoir of extraintestinal VGs and may increase the potential risk of extraintestinal infections.

Association of phylogenetic distribution and presence of integrons with multidrug resistance in Escherichia coli clinical isolates from children with diarrhoea

BACKGROUND

Escherichia coli strains include both commensal and virulent clones distributed in different phylogenetic groups. Antimicrobial resistance is an increasingly serious public health threat at the global level and integrons are important mobile genetic elements involved in resistance dissemination. This paper aims to determine the phylogenetic groups and presence of class 1 (intl1) and 2 (intl2) integrons in E. coli clinical isolates from children with diarrhoea, and to associate these characteristics with their antimicrobial resistance.

METHODS

Phylogeny and presence of integrons (intl1 and intl2) were analysed by PCR and amplicon sequencing in 70 E. coli isolates from children with and without diarrhoea (35 of each group) from Sinaloa, Mexico; these variables were analysed for correlation with the antimicrobial resistance profile of the isolates.

RESULTS

The most frequent phylogroups were A (42.9%) and B2 (15.7%). The E. coli isolates from children with diarrhoea were distributed in all phylogroups; while strains from children without diarrhoea were absent from phylogroups C, E, and clade I. The 17.1% of the isolates carried integrons (15.7% intI1 and 1.4% intI2); 28.6% of the isolates from children with diarrhoea showed the class 1 integron. Strains of phylogroup A showed the highest frequency of integrons (33.3%). The association of multidrug resistance and the presence of integrons was identified in 58.3% of strains isolated from children with diarrhoea included in phylogroups A and B2. The sequence analysis of intl1 and intl2 showed silent point mutations and similarities with plasmids of some APEC and AIEC strains.

CONCLUSION

Commensal E. coli strains are potential disseminators of antimicrobial resistance, and the improvement in the use of antimicrobials to treat childhood diarrhoea is essential for the control of such resistance.

Draft Genome Sequences of 16 Strains of Escherichia Cryptic Clade II Isolated from Intertidal Sediment in Hong Kong

The genus Escherichia includes several cryptic clades. Among them, the members of cryptic clade II have rarely been found, and their genome sequences remain largely uninvestigated. Here, we report the draft genome sequences of 16 strains of Escherichia cryptic clade II that were isolated from intertidal sediment in Hong Kong.

Coli Surface Antigen 26 Acts as an Adherence Determinant of Enterotoxigenic Escherichia coli and Is Cross-Recognized by Anti-CS20 Antibodies

The coli surface antigen 26 (CS26) of enterotoxigenic Escherichia coli (ETEC) had been described as a putative adhesive pilus based on the partial sequence of the crsH gene, detected in isolates from children with diarrhea in Egypt. However, its production and activity as adherence determinant has not been experimentally addressed. The crsH was identified as a homolog of genes encoding structural subunits of ETEC colonization factors (CFs) CS12, CS18, and CS20. These CFs, along with the recently discovered CS30, belong to the γ₂ family of pili assembled by the chaperone-usher pathway (CU pili). Further, the complete CS26 locus, crsHBCDEFG, was described in an O141 ETEC strain (ETEC 100664) obtained from a diarrhea case in The Gambia, during the Global Enterics Multicenter Study. Here, we report that CS26 is a pilus of ∼10 nm in diameter, with the capacity to increase the cell adherence of the non-pathogenic strain E. coli DH10B. As for other related pili, production of CS26 seems to be regulated by phase variation. Deletion of crsHBCDEFG in ETEC 100664 significantly decreased its adherence capacity, which was recovered by in trans complementation. Furthermore, CrsH was cross-recognized by polyclonal antibodies directed against the major structural subunit of CS20, CsnA, as determined by Western blotting and immunogold labeling. ETEC CS26+ strains were found to harbor the heat-labile enterotoxin only, within three different sequence types of phylogroups A and B1, the latter suggesting acquisition through independent events of horizontal transfer. Overall, our results demonstrate that CS26 is an adhesive pilus of human ETEC. In addition, cross-reactivity with anti-CsnA antibodies indicate presence of common epitopes in γ₂-CFs.

Fecal contamination in shallow temperate estuarine lagoon: Source of the pollution and environmental factors

In inner coastal waters such as lagoons, which are very turbid and rich in suspended particles, the survival of fecal pollution microorganisms may find favorable environments. In order to better characterize the sources and dynamics of fecal pollution in a strongly turbid environment, in situ observations were made in the Curonian Lagoon. A combination of traditional monitoring and molecular methods were used. To monitor the water quality, the indicator Escherichia coli (EC) was selected as a proxy for fecal contamination. E. coli concentration correlated with environmental parameters as pH, oxygen and turbidity. The main pollution sources are the sewage outlets in the lagoon area, while the pollution coming via rivers did not play a significant role. Still the human associated E. coli consisted only of 0 up to 20% of analyzed isolates, and did not correlate with the E. coli concentrations in the study sites. The role of birds, especially for potentially virulent E. coli may be underestimated in the lagoon.

ClermonTyping: an easy-to-use and accurate in silico method for Escherichia genus strain phylotyping

The genus Escherichia is composed of Escherichia albertii, E. fergusonii, five cryptic Escherichia clades and E. coli sensu stricto. Furthermore, the E. coli species can be divided into seven main phylogroups termed A, B1, B2, C, D, E and F. As specific lifestyles and/or hosts can be attributed to these species/phylogroups, their identification is meaningful for epidemiological studies. Classical phenotypic tests fail to identify non-sensu stricto E. coli as well as phylogroups. Clermont and colleagues have developed PCR assays that allow the identification of most of these species/phylogroups, the triplex/quadruplex PCR for E. coli phylogroup determination being the most popular. With the growing availability of whole genome sequences, we have developed the ClermonTyping method and its associated web-interface, the ClermonTyper, that allows a given strain sequence to be assigned to E. albertii, E. fergusonii, Escherichia clades I-V, E. coli sensu stricto as well as to the seven main E. coli phylogroups. The ClermonTyping is based on the concept of in vitro PCR assays and maintains the principles of ease of use and speed that prevailed during the development of the in vitro assays. This in silico approach shows 99.4 % concordance with the in vitro PCR assays and 98.8 % with the Mash genome-clustering tool. The very few discrepancies result from various errors occurring mainly from horizontal gene transfers or SNPs in the primers. We propose the ClermonTyper as a freely available resource to the scientific community at: http://clermontyping.iame-research.center/.

ClermonTyping: an easy-to-use and accurate in silico method for Escherichia genus strain phylotyping

The genus Escherichia is composed of Escherichia albertii, E. fergusonii, five cryptic Escherichia clades and E. coli sensu stricto. Furthermore, the E. coli species can be divided into seven main phylogroups termed A, B1, B2, C, D, E and F. As specific lifestyles and/or hosts can be attributed to these species/phylogroups, their identification is meaningful for epidemiological studies. Classical phenotypic tests fail to identify non-sensu stricto E. coli as well as phylogroups. Clermont and colleagues have developed PCR assays that allow the identification of most of these species/phylogroups, the triplex/quadruplex PCR for E. coli phylogroup determination being the most popular. With the growing availability of whole genome sequences, we have developed the ClermonTyping method and its associated web-interface, the ClermonTyper, that allows a given strain sequence to be assigned to E. albertii, E. fergusonii, Escherichia clades I-V, E. coli sensu stricto as well as to the seven main E. coli phylogroups. The ClermonTyping is based on the concept of in vitro PCR assays and maintains the principles of ease of use and speed that prevailed during the development of the in vitro assays. This in silico approach shows 99.4 % concordance with the in vitro PCR assays and 98.8 % with the Mash genome-clustering tool. The very few discrepancies result from various errors occurring mainly from horizontal gene transfers or SNPs in the primers. We propose the ClermonTyper as a freely available resource to the scientific community at: http://clermontyping.iame-research.center/.