Large-scale genome analysis of bovine commensal Escherichia coli reveals that bovine-adapted E. coli lineages are serving as evolutionary sources of the emergence of human intestinal pathogenic strains

Characterization of diarrheagenic Escherichia coli from different cattle production systems in Brazil

Diarrheagenic E. coli (DEC) can cause severe diarrhea and is a public health concern worldwide. Cattle are an important reservoir for this group of pathogens, and once introduced into the abattoir environment, these microorganisms can contaminate consumer products. This study aimed to characterize the distribution of DEC [Shiga toxin-producing E. coli (STEC), enteroinvasive E. coli (EIEC), enteropathogenic E. coli (EPEC), enterotoxigenic E. coli (ETEC), and enteroaggregative E. coli (EAEC)] from extensive and intensive cattle production systems in Brazil. Samples (n = 919) were collected from animal feces (n = 200), carcasses (n = 600), meat cuts (n = 90), employee feces (n = 9), and slaughterhouse water (n = 20). Virulence genes were detected by PCR in 10% of animal samples (94/919), with STEC (n = 81) as the higher prevalence, followed by EIEC (n = 8), and lastly EPEC (n = 5). Animals raised in an extensive system had a higher prevalence of STEC (average 48%, sd = 2.04) when compared to animals raised in an intensive system (23%, sd = 1.95) (Chi-square test, P < 0.001). From these animals, most STEC isolates only harbored stx2 (58%), and 7% were STEC LEE-positive isolates that were further identified as O157:H7. This study provides further evidence that cattle are potential sources of DEC, especially STEC, and that potentially pathogenic E. coli isolates are widely distributed in feces and carcasses during the slaughter process.

Characterization of Escherichia coli pathogenicity and drug resistance in yolk peritonitis

Yolk Peritonitis can lead to a rapid decline in egg production, which seriously affects the health of laying hens and the profitability of chicken farms. Escherichia coli (E. coli) is the most common cause of yolk peritonitis in laying hens. In this study, bacterial samples were collected from the ovaries and fallopian tubes of laying hens with suspected yolk peritonitis from a laying farm in Jiangsu Province, and their pathogenicity and drug resistance were investigated. Initially, morphological and biochemical detection methods were employed to isolate and identify the pathogenic bacteria. The results showed that a total of 16 strains of E. coli were isolated from laying hens with yolk peritonitis. Subsequently, the drug resistance and pathogenicity of a randomly selected E. coli strain were analyzed and predicted by genome sequencing technology, and the drug resistance of E. coli was verified by drug sensitivity test and PCR. Finally, the virulence was verified by infection experiment in mice. The study revealed that the egg-yolk peritonitis in laying hens was caused by E. coli infection, and the genome sequencing analysis revealed that the bacteria had multidrug resistance and high virulence. The drug susceptibility testing indicates that E. coli exhibited resistance to aminoglycosides, β-lactam, macrolides, fluoroquinolones, and sulfonamides. In this study, resistance genes including KdpE, aadA5, APH(3 ")-ID, APH(6)-ID, and TEM-1 were identified, and their expression levels varied across different stages of bacterial growth. The results of virulence analysis indicated a mortality rate of 50% in mice infected with E. coli at a concentration of 2.985 × 107 CFU/mL. E. coli infection resulted in damage to various tissues and organs in mice, with the intestinal tissue structure being the most severely affected. This study provides a reference for the study of drug resistance mechanisms in E. coli and provides valuable insights into the selection of drugs for the treatment of vitelline peritonitis.

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

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

Vaccines against extraintestinal pathogenic Escherichia coli (ExPEC): progress and challenges

The emergence of antimicrobial resistance (AMR) is a principal global health crisis projected to cause 10 million deaths annually worldwide by 2050. While the Gram-negative bacteria Escherichia coli is commonly found as a commensal microbe in the human gut, some strains are dangerously pathogenic, contributing to the highest AMR-associated mortality. Strains of E. coli that can translocate from the gastrointestinal tract to distal sites, called extraintestinal E. coli (ExPEC), are particularly problematic and predominantly afflict women, the elderly, and immunocompromised populations. Despite nearly 40 years of clinical trials, there is still no vaccine against ExPEC. One reason for this is the remarkable diversity in the ExPEC pangenome across pathotypes, clades, and strains, with hundreds of genes associated with pathogenesis including toxins, adhesins, and nutrient acquisition systems. Further, ExPEC is intimately associated with human mucosal surfaces and has evolved creative strategies to avoid the immune system. This review summarizes previous and ongoing preclinical and clinical ExPEC vaccine research efforts to help identify key gaps in knowledge and remaining challenges.

happi: a hierarchical approach to pangenomics inference

Recovering metagenome-assembled genomes (MAGs) from shotgun sequencing data is an increasingly common task in microbiome studies, as MAGs provide deeper insight into the functional potential of both culturable and non-culturable microorganisms. However, metagenome-assembled genomes vary in quality and may contain omissions and contamination. These errors present challenges for detecting genes and comparing gene enrichment across sample types. To address this, we propose happi, an approach to testing hypotheses about gene enrichment that accounts for genome quality. We illustrate the advantages of happi over existing approaches using published Saccharibacteria MAGs, Streptococcus thermophilus MAGs, and via simulation.

Comparative Genomic Analyses of Escherichia coli from a Meat Processing Environment in Relation to Their Biofilm Formation and Persistence

We investigated the phylogeny of biofilm forming (BF) and nonbiofilm forming (NBF) Escherichia coli (n = 114) from a beef processing environment as well as genetic elements in their BF and persistence via a comparative genomic analysis. Phylogroup B1 made up the largest proportion of both the BF (73.8%) and NBF (50.9%) groups. E. coli from all of the sources that were examined had mixed phylogroups, except for those that were recovered from equipment after cleaning, which were exclusively from phylogroup B1. Both the core genome and gene content trees showed a tree-wide spread of BF strains, with clusters, including both BF and NBF strains. Genome-wide association studies (GWAS) via Scoary or Pyseer did not find any genes or mutations that were overrepresented in the BF group. A retrospective analysis of phenotypes found a significant correlation (P < 0.05) between BF ability and curli production, cellulose synthesis, and/or mobility. However, the BF group also included strains that were negative for curli and cellulose and/or missing encoding genes for the two traits. All curli and cellulose encoding genes were present in most genomes, regardless of their BF status. The degree of motility was correlated with both curli and cellulose production, and 80 common genes were overrepresented in all three of the trait-positive groups. A PTS enzyme II, a subsidiary gluconate catabolism pathway, and an iron-dicitrate transport system were more abundant in the persisting E. coli group. These findings suggest gene function redundancy in E. coli for biofilm formation as well as additional substrate utilization and iron acquisition in its persistence. IMPORTANCE The persistence of potentially hazardous bacteria is a major challenge for meat processing environments, which are conducive for biofilm formation. Marker genes/phenotypes are commonly used to differentiate biofilm forming E. coli strains from their nonbiofilm forming counterparts. We took a comparative genomic analysis approach to analyze E. coli strains that were from the same environment but were differentiated by their biofilm forming ability. A diversification of the genes involved in the biofilm formation of E. coli was observed. Even though there is a correlation on the population level between biofilm formation and the expression of curli and cellulose, uncertainties exist on the individual strain level. Novel substrate utilization and iron acquisition could contribute to the persistence of E. coli. These findings not only advance our understanding of the ecology of E. coli with respect to its persistence but also show that a marker gene/phenotype driven approach for the biofilm control of E. coli may not be prudent.

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.

The Diversity of Escherichia coli Pathotypes and Vaccination Strategies against This Versatile Bacterial Pathogen

Escherichia coli (E. coli) is a gram-negative bacillus and resident of the normal intestinal microbiota. However, some E. coli strains can cause diseases in humans, other mammals and birds ranging from intestinal infections, for example, diarrhea and dysentery, to extraintestinal infections, such as urinary tract infections, respiratory tract infections, meningitis, and sepsis. In terms of morbidity and mortality, pathogenic E. coli has a great impact on public health, with an economic cost of several billion dollars annually worldwide. Antibiotics are not usually used as first-line treatment for diarrheal illness caused by E. coli and in the case of bloody diarrhea, antibiotics are avoided due to the increased risk of hemolytic uremic syndrome. On the other hand, extraintestinal infections are treated with various antibiotics depending on the site of infection and susceptibility testing. Several alarming papers concerning the rising antibiotic resistance rates in E. coli strains have been published. The silent pandemic of multidrug-resistant bacteria including pathogenic E. coli that have become more difficult to treat favored prophylactic approaches such as E. coli vaccines. This review provides an overview of the pathogenesis of different pathotypes of E. coli, the virulence factors involved and updates on the major aspects of vaccine development against different E. coli pathotypes.

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.

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

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

Enteropathogenic Escherichia coli (EPEC) expressing a non-functional bundle-forming pili (BFP) also leads to increased growth failure and intestinal inflammation in C57BL/6 mice

Bundle-forming pili (BFP) are implicated in the virulence of typical enteropathogenic E. coli (EPEC), resulting in enhanced colonization and mild to severe disease outcomes; hence, non-functional BFP may have a major influence on disease outcomes in vivo. Weaned antibiotic pre-treated C57BL/6 mice were orally infected with EPEC strain UMD901 (E2348/69 bfpA C129S); mice were monitored daily for body weight; stool specimens were collected daily; and intestinal tissues were collected at the termination of the experiment on day 3 post-infection. Real-time PCR was used to quantify fecal shedding and tissue burden. Intestinal inflammatory biomarkers lipocalin-2 (LCN-2) and myeloperoxidase (MPO) were also assessed. Infection caused substantial body weight loss, bloody diarrhea, and intestinal colonization with fecal and intestinal tissue inflammatory biomarkers that were comparable to those previously published with the wild-type typical EPEC strain. Here we further report on the evaluation of an EPEC infection model, showing how disruption of bfp function does not impair, and may even worsen diarrhea, colonization, and intestinal disruption and inflammation. More research is needed to understand the role of bfp in pathogenicity of EPEC infections in vivo.

Genomic analysis of Escherichia coli circulating in the Brazilian poultry sector

Escherichia coli are gut commensal bacteria and opportunistic pathogens, and the emergence of antimicrobial resistance threatens the safety of the food chain. To know the E. coli strains circulating in the Brazilian poultry sector is important since the country corresponds to a significant chicken meat production. Thus, we analyzed 90 publicly genomes available in a database using web-based tools. Genomic analysis revealed that sul alleles were the most detected resistance genes, followed by aadA, bla_CTX-M, and dfrA. Plasmids of the IncF family were important, followed by IncI1-Iα, Col-like, and p0111. Genes of specific metabolic pathways that contribute to virulence (terC and gad) were predominant, followed by sitA, traT, and iss. Additionally, pap, usp, vat, sfa/foc, ibeA, cnf1, eae, and sat were also predicted. In this regard, 11 E. coli were characterized as avian pathogenic E. coli and one as atypical enteropathogenic E. coli. Phylogenetic analysis confirmed the predominant occurrence of B1 but also A, D, B2, F, E, G, C, and Clade I phylogroups, whereas international clones ST38, ST73, ST117, ST155, and ST224 were predicted among 53 different sequence types identified. Serotypes O6:H1 and:H25 were prevalent, and fimH31 and fimH32 were the most representatives among the 36 FimH types detected. Finally, single nucleotide polymorphisms-based phylogenetic analysis confirmed high genomic diversity among E. coli strains. While international E. coli clones have adapted to the Brazilian poultry sector, the virulome background of these strains support a pathogenic potential to humans and animals, with lineages carrying resistance genes that can lead to hard-to-treat infections.

Genomic Diversity, Virulence Gene, and Prophage Arrays of Bovine and Human Shiga Toxigenic and Enteropathogenic Escherichia coli Strains Isolated in Hungary

Escherichia coli belonging to the enterohemorrhagic (EHEC), Shiga toxin-producing (STEC) and atypical enteropathogenic (aEPEC) pathotypes are significant foodborne zoonotic pathogens posing serious health risks, with healthy cattle as their main reservoir. A representative sampling of Hungarian cattle farms during 2017–2018 yielded a prevalence of 6.5 and 5.8% for STEC and aEPEC out of 309 samples. The draft genomes of twelve STEC (of them 9 EHEC) and four aEPEC of bovine origin were determined. For comparative purposes, we also included 3 EHEC and 2 aEPEC strains of human origin, as well four commensal isolates and one extraintestinal pathogenic E. coli (ExPEC) obtained from animals in a final set of 26 strains for a WGS-based analysis. Apart from key virulence genes, these isolates harbored several additional virulence genes with arrays characteristic for the site of isolation. The most frequent insertion site of Shiga toxin (stx) encoding prophages was yehV for the Stx1 prophage and wrbA and sbcB for Stx2. For O157:H7 strains, the locus of enterocyte effacement pathogenicity island was present at the selC site, with integration at pheV for other serotypes, and pheU in the case of O26:H11 strains. Several LEE-negative STEC and aEPEC as well as commensal isolates carried additional prophages, with an average of ten prophage regions per isolate. Comparative phylogenomic analysis showed no clear separation between bovine and human lineages among the isolates characterized in the current study. Similarities in virulence gene arrays and close phylogenetic relations of bovine and human isolates underline the zoonotic potential of bovine aEPEC and STEC and emphasize the need for frequent monitoring of these pathogens in livestock.

Comparative Genomics Applied to Systematically Assess Pathogenicity Potential in Shiga Toxin-Producing Escherichia coli O145:H28

Shiga toxin-producing Escherichia coli (STEC) O145:H28 can cause severe disease in humans and is a predominant serotype in STEC O145 environmental isolates. Here, comparative genomics was applied to a set of clinical and environmental strains to systematically evaluate the pathogenicity potential in environmental strains. While the core genes-based tree separated all O145:H28 strains from the non O145:H28 reference strains, it failed to segregate environmental strains from the clinical. In contrast, the accessory genes-based tree placed all clinical strains in the same clade regardless of their genotypes or serotypes, apart from the environmental strains. Loss-of-function mutations were common in the virulence genes examined, with a high frequency in genes related to adherence, autotransporters, and the type three secretion system. Distinct differences in pathogenicity islands LEE, OI-122, and OI-57, the acid fitness island, and the tellurite resistance island were detected between the O145:H28 and reference strains. A great amount of genetic variation was detected in O145:H28, which was mainly attributed to deletions, insertions, and gene acquisition at several chromosomal “hot spots”. Our study demonstrated a distinct virulence gene repertoire among the STEC O145:H28 strains originating from the same geographical region and revealed unforeseen contributions of loss-of-function mutations to virulence evolution and genetic diversification in STEC.

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.

Antibiotic resistance and phylogenetic profiling of Escherichia coli from dairy farm soils; organic versus conventional systems

•

First known comparison of antimicrobial resistance traits in E. coli strains from new zealand farms practicing organic and conventional husbandry.

•

Potential extended spectrum β-lactamase producing strains isolated from dairy farm environments.

•

Organic dairy farms tended to harbour fewer resistant isolates than those recovered from conventionally farmed counterparts.

•

Evidence for anthroponotic transmission of resistant strains of human origin to farm environments.

•

Implications for the spread of antimicrobial resistance traits from farm environments discussed.

The prevalence and spread of antimicrobial resistance (AMR) as a result of the persistent use and/or abuse of antimicrobials is a key health problem for health authorities and governments worldwide. A study of contrasting farming systems such as organic versus conventional dairy farming may help to authenticate some factors that may contribute to the prevalence and spread of AMR in their soils. A case study was conducted in organic and conventional dairy farms in the South Canterbury region of New Zealand.

A total of 814 dairy farm soil E. coli (DfSEC) isolates recovered over two years were studied. Isolates were recovered from each of two farms practicing organic, and another two practicing conventional husbandries. The E. coli isolates were examined for their antimicrobial resistance (AMR) against cefoxitin, cefpodoxime, chloramphenicol, ciprofloxacin, gentamicin, meropenem, nalidixic acid, and tetracycline. Phylogenetic relationships were assessed using an established multiplex PCR method. The AMR results indicated 3.7% of the DfSEC isolates were resistant to at least one of the eight selected antimicrobials. Of the resistant isolates, DfSEC from the organic dairy farms showed a lower prevalence of resistance to the antimicrobials tested, compared to their counterparts from the conventional farms. Phylogenetic analysis placed the majority (73.7%) of isolates recovered in group B1, itself dominated by isolates of bovine origin. The tendency for higher rates of resistance among strains from conventional farming may be important for future decision-making around farming practices Current husbandry practices may contribute to the prevalence and spread of AMR in the industry.

[Genome analysis-based studies on bacterial genetic diversity]

There are a huge number of bacterial species on earth, and a huge intra-species genomic diversity are also observed in many bacteria. The high ability of bacteria to acquire foreign DNA and the presence of various mobile genetic elements contribute the generation of such genomic diversity. During the biochemical and genetic analysis of a Pseudomonas aeruginosa toxin, called cytotoxin, and its converting phage, which I first engaged in my research carrier, I became very interested in the genetic diversity of bacteria and mobile genetic elements such as bacteriophages, and realized the usefulness and power of genome analysis. Since then, I have been involved in genome analyses of various pathogenic bacteria such as enterohemorrhagic Escherichia coli (EHEC), commensal bacteria of human and other animals, and bacteria or bacterial communities in natural environments. I was so lucky that I jumped in this research field at the very begging of genome analyses and experienced a very exciting time of surprisingly rapid advance in genome sequencing technologies which revolutionized a wide range of biology. In this article, I first review the main findings which our group obtained from the genome analyses on the P. aeruginosa cytotoxin converting phage and those on the evolution and genomic diversity of EHEC and related bacteria. The results of our analyses of Rickettsiaceae family genomes, which show surprisingly very low genomic diversity, and genome sequence-based analyses of an intrahospital bacterial outbreak and within-host genomic diversity are also summarized.

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

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

Development of a homogeneous time-resolved FRET (HTRF) assay for the quantification of Shiga toxin 2 produced by E. coli

Shiga toxin-producing Escherichia coli (STEC) is a major intestinal pathogen and causes serious gastrointestinal illness, which includes diarrhea, hemorrhagic colitis, and life-threatening hemolytic uremic syndrome. The major virulence factors of STEC are Shiga toxins (Stx1 and Stx2), which belong to the AB-type toxin family. Among several subtypes of Stx1 and Stx2, the production of Stx2a is thought to be a risk factor for severe STEC infections, but Stx2a production levels vary markedly between STEC strains, even strains with the same serotype. Therefore, quantitative analyses of Stx2 production by STEC strains are important to understand the virulence potential of specific lineages or sublineages. In this study, we developed a novel Stx2 quantification method by utilizing homogeneous time-resolved fluorescence resonance energy transfer (HTRF) technology. To determine suitable “sandwich” assay conditions, we tested 6 combinations of fluorescence-labeled monoclonal antibodies (mAbs) specific to Stx2 and compared the HTRF signal intensities obtained at various incubation times. Through this analysis, we selected the most suitable mAb pair, one recognizing the A subunit and the other recognizing the B subunit, thus together detecting Stx holotoxins. The optimal incubation time was also determined (18 h). Then, we optimized the concentrations of the two mAbs based on the range for linearity. The established HTRF assay detected 0.5 ng/ml of the highly purified recombinant Stx2a and Stx2e proteins and the working range was 1–64 ng/ml for both Stx2a and Stx2e. Through the quantification analysis of Stx proteins in STEC cell lysates, we confirmed that other Stx2 subtypes (Stx2b, Stx2c, Stx2d and Stx2g) can also be quantified at a certain level of accuracy, while this assay system does not detect Stx2f, which is highly divergent in sequence from other Stx2 subtypes, and Stx1. As the HTRF protocol we established is simple, this assay system should prove useful for the quantitative analysis of Stx2 production levels of a large number of STEC strains.

Comparative genomic analysis of Escherichia coli isolates from cases of bovine clinical mastitis identifies nine specific pathotype marker genes

Escherichia coli is a major causative agent of environmental bovine mastitis and this disease causes significant economic losses for the dairy industry. There is still debate in the literature as to whether mammary pathogenic E. coli (MPEC) is indeed a unique E. coli pathotype, or whether this infection is merely an opportunistic infection caused by any E. coli isolate being displaced from the bovine gastrointestinal tract to the environment and, then, into the udder. In this study, we conducted a thorough genomic analysis of 113 novel MPEC isolates from clinical mastitis cases and 100 bovine commensal E. coli isolates. A phylogenomic analysis indicated that MPEC and commensal E. coli isolates formed clades based on common sequence types and O antigens, but did not cluster based on mammary pathogenicity. A comparative genomic analysis of MPEC and commensal isolates led to the identification of nine genes that were part of either the core or the soft-core MPEC genome, but were not found in any bovine commensal isolates. These apparent MPEC marker genes were genes involved with nutrient intake and metabolism [adeQ, adenine permease; nifJ, pyruvate-flavodoxin oxidoreductase; and yhjX, putative major facilitator superfamily (MFS)-type transporter], included fitness and virulence factors commonly seen in uropathogenic E. coli (pqqL, zinc metallopeptidase, and fdeC, intimin-like adhesin, respectively), and putative proteins [yfiE, uncharacterized helix-turn-helix-type transcriptional activator; ygjI, putative inner membrane transporter; and ygjJ, putative periplasmic protein]. Further characterization of these highly conserved MPEC genes may be critical to understanding the pathobiology of MPEC.

Platanus_B: an accurate de novo assembler for bacterial genomes using an iterative error-removal process

De novo assembly of short DNA reads remains an essential technology, especially for large-scale projects and high-resolution variant analyses in epidemiology. However, the existing tools often lack sufficient accuracy required to compare closely related strains. To facilitate such studies on bacterial genomes, we developed Platanus_B, a de novo assembler that employs iterations of multiple error-removal algorithms. The benchmarks demonstrated the superior accuracy and high contiguity of Platanus_B, in addition to its ability to enhance the hybrid assembly of both short and nanopore long reads. Although the hybrid strategies for short and long reads were effective in achieving near full-length genomes, we found that short-read-only assemblies generated with Platanus_B were sufficient to obtain ≥90% of exact coding sequences in most cases. In addition, while nanopore long-read-only assemblies lacked fine-scale accuracies, inclusion of short reads was effective in improving the accuracies. Platanus_B can, therefore, be used for comprehensive genomic surveillances of bacterial pathogens and high-resolution phylogenomic analyses of a wide range of bacteria.

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

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

Prophages integrating into prophages: A mechanism to accumulate type III secretion effector genes and duplicate Shiga toxin-encoding prophages in Escherichia coli

Bacteriophages (or phages) play major roles in the evolution of bacterial pathogens via horizontal gene transfer. Multiple phages are often integrated in a host chromosome as prophages, not only carrying various novel virulence-related genetic determinants into host bacteria but also providing various possibilities for prophage-prophage interactions in bacterial cells. In particular, Escherichia coli strains such as Shiga toxin (Stx)-producing E. coli (STEC) and enteropathogenic E. coli (EPEC) strains have acquired more than 10 prophages (up to 21 prophages), many of which encode type III secretion system (T3SS) effector gene clusters. In these strains, some prophages are present at a single locus in tandem, which is usually interpreted as the integration of phages that use the same attachment (att) sequence. Here, we present phages integrating into T3SS effector gene cluster-associated loci in prophages, which are widely distributed in STEC and EPEC. Some of the phages integrated into prophages are Stx-encoding phages (Stx phages) and have induced the duplication of Stx phages in a single cell. The identified attB sequences in prophage genomes are apparently derived from host chromosomes. In addition, two or three different attB sequences are present in some prophages, which results in the generation of prophage clusters in various complex configurations. These phages integrating into prophages represent a medically and biologically important type of inter-phage interaction that promotes the accumulation of T3SS effector genes in STEC and EPEC, the duplication of Stx phages in STEC, and the conversion of EPEC to STEC and that may be distributed in other types of E. coli strains as well as other prophage-rich bacterial species.

Multiple prophages are often integrated in a bacterial host chromosome and some are present at a single locus in tandem. The most striking examples are Shiga toxin (Stx)-producing and enteropathogenic Escherichia coli (STEC and EPEC) strains, which usually contain more than 10 prophages (up to 21). Many of them encode a cluster of type III secretion system (T3SS) effector genes, contributing the acquisition of a large number of effectors (>30) by STEC and EPEC. Here, we describe prophages integrating into T3SS effector gene cluster-associated loci in prophages, which are widely distributed in STEC and EPEC. Two or three different attachment sequences derived from host chromosomes are present in some prophages, generating prophage clusters in various complex configurations. Of note, some of such phages integrating into prophages are Stx-encoding phages (Stx phages) and have induced the duplication of Stx phages. Thus, “prophage-in-prophage” represents an important inter-phage interaction as they can promote not only the accumulation of T3SS effectors in STEC and EPEC but also the duplication of Stx phages and the conversion of EPEC to STEC.

Combining Functional Genomics and Whole-Genome Sequencing to Detect Antibiotic Resistance Genes in Bacterial Strains Co-Occurring Simultaneously in a Brazilian Hospital

(1) Background: The rise of multi-antibiotic resistant bacteria represents an emergent threat to human health. Here, we investigate antibiotic resistance mechanisms in bacteria of several species isolated from an intensive care unit in Brazil. (2) Methods: We used whole-genome analysis to identify antibiotic resistance genes (ARGs) and plasmids in 34 strains of Gram-negative and Gram-positive bacteria, providing the first genomic description of Morganella morganii and Ralstonia mannitolilytica clinical isolates from South America. (3) Results: We identified a high abundance of beta-lactamase genes in resistant organisms, including seven extended-spectrum beta-lactamases (OXA-1, OXA-10, CTX-M-1, KPC, TEM, HYDRO, BLP) shared between organisms from different species. Additionally, we identified several ARG-carrying plasmids indicating the potential for a fast transmission of resistance mechanism between bacterial strains. Furthermore, we uncovered two pairs of (near) identical plasmids exhibiting multi-drug resistance. Finally, since many highly resistant strains carry several different ARGs, we used functional genomics to investigate which of them were indeed functional. In this sense, for three bacterial strains (Escherichia coli, Klebsiella pneumoniae, and M. morganii), we identified six beta-lactamase genes out of 15 predicted in silico as those mainly responsible for the resistance mechanisms observed, corroborating the existence of redundant resistance mechanisms in these organisms. (4) Conclusions: Systematic studies similar to the one presented here should help to prevent outbreaks of novel multidrug-resistant bacteria in healthcare facilities.

Genome evolution and the emergence of pathogenicity in avian Escherichia coli

Chickens are the most common birds on Earth and colibacillosis is among the most common diseases affecting them. This major threat to animal welfare and safe sustainable food production is difficult to combat because the etiological agent, avian pathogenic Escherichia coli (APEC), emerges from ubiquitous commensal gut bacteria, with no single virulence gene present in all disease-causing isolates. Here, we address the underlying evolutionary mechanisms of extraintestinal spread and systemic infection in poultry. Combining population scale comparative genomics and pangenome-wide association studies, we compare E. coli from commensal carriage and systemic infections. We identify phylogroup-specific and species-wide genetic elements that are enriched in APEC, including pathogenicity-associated variation in 143 genes that have diverse functions, including genes involved in metabolism, lipopolysaccharide synthesis, heat shock response, antimicrobial resistance and toxicity. We find that horizontal gene transfer spreads pathogenicity elements, allowing divergent clones to cause infection. Finally, a Random Forest model prediction of disease status (carriage vs. disease) identifies pathogenic strains in the emergent ST-117 poultry-associated lineage with 73% accuracy, demonstrating the potential for early identification of emergent APEC in healthy flocks.

Enteropathogenic Escherichia coli—A Summary of the Literature

Diarrheal disease is still a major public health concern, as it is still considered an important cause of death in children under five years of age. A few decades ago, the detection of enteropathogenic E. coli was made by detecting the O, H, and K antigens, mostly by agglutination. The recent protocols recommend the molecular methods for diagnosing EPEC, as they can distinguish between typical and atypical EPEC by identifying the presence/absence of specific virulence factors. EPEC are defined as diarrheagenic strains of E. coli that can produce attaching and effacing lesions on the intestinal epithelium while being incapable of producing Shiga toxins and heat-labile or heat-stable enterotoxins. The ability of these strains to produce attaching and effacing lesions enable them to cause localized lesions by attaching tightly to the surface of the intestinal epithelial cells, disrupting the surfaces of the cells, thus leading to the effacement of the microvilli. EPEC are classified on typical and atypical isolates, based on the presence or absence of E. coli adherence factor plasmids. All the EPEC strains are eae positive; typical EPEC strains are eae+, bfpA+, while atypical strains are eae+, bfpA−. No vaccines are currently available to prevent EPEC infections.

Genomic Epidemiology and Evolution of Escherichia coli in Wild Animals in Mexico

Escherichia coli is a clinically important bacterial species implicated in human- and livestock-associated infections worldwide. The bacterium is known to reside in the guts of humans, livestock, and wild animals.

Escherichia coli is a common bacterial species in the gastrointestinal tracts of warm-blooded animals and humans. Pathogenicity and antimicrobial resistance in E. coli may emerge via host switching from animal reservoirs. Despite its potential clinical importance, knowledge of the population structure of commensal E. coli within wild hosts and the epidemiological links between E. coli in nonhuman hosts and E. coli in humans is still scarce. In this study, we analyzed the whole-genome sequencing data of a collection of 119 commensal E. coli strains recovered from the guts of 55 mammal and bird species in Mexico and Venezuela in the 1990s. We observed low concordance between the population structures of E. coli isolates colonizing wild animals and the phylogeny, taxonomy, and ecological and physiological attributes of the host species, with distantly related E. coli strains often colonizing the same or similar host species and distantly related host species often hosting closely related E. coli strains. We found no evidence for recent transmission of E. coli genomes from wild animals to either domesticated animals or humans. However, multiple livestock- and human-related virulence factor genes were present in E. coli of wild animals, including virulence factors characteristic of Shiga toxin-producing E. coli (STEC) and atypical enteropathogenic E. coli (aEPEC), where several isolates from wild hosts harbored the locus of enterocyte effacement (LEE) pathogenicity island. Moreover, E. coli isolates from wild animal hosts often harbored known antibiotic resistance determinants, including those against ciprofloxacin, aminoglycosides, tetracyclines, and beta-lactams, with some determinants present in multiple, distantly related E. coli lineages colonizing very different host animals. We conclude that genome pools of E. coli colonizing the guts of wild animals and humans share virulence and antibiotic resistance genes, underscoring the idea that wild animals could serve as reservoirs for E. coli pathogenicity in human and livestock infections.

IMPORTANCEEscherichia coli is a clinically important bacterial species implicated in human- and livestock-associated infections worldwide. The bacterium is known to reside in the guts of humans, livestock, and wild animals. Although wild animals are recognized as potential reservoirs for pathogenic E. coli strains, the knowledge of the population structure of E. coli in wild hosts is still scarce. In this study, we used fine resolution of whole-genome sequencing to provide novel insights into the evolution of E. coli genomes from a small yet diverse collection of strains recovered within a broad range of wild animal species (including mammals and birds), the coevolution of E. coli strains with their hosts, and the genetics of pathogenicity of E. coli strains in wild hosts in Mexico. Our results provide evidence for the clinical importance of wild animals as reservoirs for pathogenic strains and highlight the need to include nonhuman hosts in the surveillance programs for E. coli infections.

Shiga Toxin-Producing E. coli in Animals: Detection, Characterization, and Virulence Assessment

Cattle and other ruminants are primary reservoirs for Shiga toxin-producing Escherichia coli (STEC) strains which have a highly variable, but unpredictable, pathogenic potential for humans. Domestic swine can carry and shed STEC, but only STEC strains producing the Shiga toxin (Stx) 2e variant and causing edema disease in piglets are considered pathogens of veterinary medical interest. In this chapter, we present general diagnostic workflows for sampling livestock animals to assess STEC prevalence, magnitude, and duration of host colonization. This is followed by detailed method protocols for STEC detection and typing at genetic and phenotypic levels to assess the relative virulence exerted by the strains.

Metabolic Traits of Bovine Shiga Toxin-Producing Escherichia Coli (STEC) Strains with Different Colonization Properties

Cattle harbor Shiga toxin-producing Escherichia coli (STEC) in their intestinal tract, thereby providing these microorganisms with an ecological niche, but without this colonization leading to any clinical signs. In a preceding study, genotypic characterization of bovine STEC isolates unveiled that their ability to colonize cattle persistently (STEC^per) or only sporadically (STEC^spo) is more closely associated with the overall composition of the accessory rather than the core genome. However, the colonization pattern could not be unequivocally linked to the possession of classical virulence genes. This study aimed at assessing, therefore, if the presence of certain phenotypic traits in the strains determines their colonization pattern and if these can be traced back to distinctive genetic features. STEC^spo strains produced significantly more biofilm than STEC^per when incubated at lower temperatures. Key substrates, the metabolism of which showed a significant association with colonization type, were glyoxylic acid and L-rhamnose, which were utilized by STEC^spo, but not or only by some STEC^per. Genomic sequences of the respective glc and rha operons contained mutations and frameshifts in uptake and/or regulatory genes, particularly in STEC^per. These findings suggest that STEC^spo conserved features leveraging survival in the environment, whereas the acquisition of a persistent colonization phenotype in the cattle reservoir was accompanied by the loss of metabolic properties and genomic mutations in the underlying genetic pathways.

Structural and Functional Characterization of Stx2k, a New Subtype of Shiga Toxin 2

Shiga toxin (Stx) is the major virulence factor of Shiga toxin-producing Escherichia coli (STEC). Stx evolves rapidly and, as such, new subtypes continue to emerge that challenge the efficacy of existing disease management and surveillance strategies. A new subtype, Stx2k, was recently identified in E. coli isolated from a wide range of sources including diarrheal patients, animals, and raw meats, and was poorly detected by existing immunoassays. In this study, the structure of Stx2kE167Q was determined at 2.29 Å resolution and the conservation of structure with Stx2a was revealed. A novel polyclonal antibody capable of neutralizing Stx2k and an immunoassay, with a 10-fold increase in sensitivity compared to assays using extant antibodies, were developed. Stx2k is less toxic than Stx2a in Vero cell assays but is similar to Stx2a in receptor-binding preference, thermostability, and acid tolerance. Although Stx2k does not appear to be as potent as Stx2a to Vero cells, the wide distribution and blended virulence profiles of the Stx2k-producing strains suggest that horizontal gene transfer through Stx2k-converting phages could result in the emergence of new and highly virulent pathogens. This study provides useful information and tools for early detection and control of Stx2k-producing E. coli, which could reduce public risk of infection by less-known STECs.