Comparative genomics of two super-shedder isolates of Escherichia coli O157:H7

Investigating the adherence factors of Escherichia coli at the bovine recto-anal junction

Shiga toxin-producing Escherichia coli (STEC) are major foodborne pathogens that result in thousands of hospitalizations each year in the United States. Cattle, the natural reservoir, harbor STEC asymptomatically at the recto-anal junction (RAJ). The molecular mechanisms that allow STEC and non-STEC E. coli to adhere to the RAJ are not fully understood, in part because most adherence studies utilize human cell culture models. To identify a set of bovine-specific E. coli adherence factors, we used the primary RAJ squamous epithelial (RSE) cell-adherence assay to coculture RSE cells from healthy Holstein cattle with diverse E. coli strains from bovine and nonbovine sources. We hypothesized that a comparative genomic analysis of the strains would reveal factors associated with RSE adherence. After performing adherence assays with historical strains from the E. coli Reference Center (n = 62) and strains newly isolated from the RAJ (n = 15), we used the bioinformatic tool Roary to create a pangenome of this collection. We classified strains as either low or high adherence and using the Scoary program compiled a list of accessory genes correlated with the "high adherence" strains. While none of the correlations were statistically significant, several gene clusters were associated with the high-adherence phenotype, including two that encode uncharacterized proteins. We also demonstrated that non-STEC E. coli strains from the RAJ are more adherent than other isolates and can outcompete STEC in coculture with RSEs. Further analysis of adherence-associated gene clusters may lead to an improved understanding of the molecular mechanisms of RSE adherence and may help develop probiotics targeting STEC in cattle.

IMPORTANCE

E. coli strains that produce Shiga toxin cause foodborne illness in humans but colonize cattle asymptomatically. The molecular mechanisms that E. coli uses to adhere to cattle cells are largely unknown. Various strategies are used to control E. coli in livestock and limit the risk of outbreaks. These include vaccinating animals against common E. coli strains and supplementing their feed with probiotics to reduce the carriage of pathogens. No strategy is completely effective, and probiotics often fail to colonize the animals. We sought to clarify the genes required for E. coli adherence in cattle by quantifying the attachment to bovine cells in a diverse set of bacteria. We also isolated nonpathogenic E. coli from healthy cows and showed that a representative isolate could outcompete pathogenic strains in cocultures. We propose that the focused study of these strains and their adherence factors will better inform the design of probiotics and vaccines for livestock.

In Silico Analysis of Shiga Toxin-Producing Escherichia coli O157:H7 Strains from Presumptive Super- and Low-Shedder Cattle

Cattle are the primary reservoir for STEC O157, with some shedding >104 CFU/g in feces, a phenomenon known as super-shedding (SS). The mechanism(s) responsible for SS are not understood but have been attributed to the environment, host, and pathogen. This study aimed to compare genetic characteristics of STEC O157 strains from cattle in the same commercial feedlot pens with SS or low-shedding (LS) status. Strains from SS (n = 35) and LS (n = 28) collected from 11 pens in three feedlots were analyzed for virulence genes, Shiga toxin-carrying bacteriophage insertion sites, and phylogenetic relationships. In silico analysis showed limited variation regarding virulence gene profiles. Stx-encoding prophage insertion sites mrlA and wrbA for stx1a and stx2a, respectively, were all occupied, but two isolates had fragments of the stx-carrying phage in mrlA and wrbA loci without stx1a and stx2a. All strains screened for lineage-specific polymorphism assay (LSPA-6) were 111111, lineage I. Of the isolates, 61 and 2 were clades 1 and 8, respectively. Phylogenetic analysis revealed that pens with more than one SS had multiple distantly related clusters of SS and LS isolates. Although virulence genes and lineage were largely similar within and across feedlots, multiple genetic origins of strains within a single feedlot pen illustrate challenges for on-farm control of STEC.

Comparative genomics analysis and characterization of Shiga toxin-producing Escherichia coli O157:H7 strains reveal virulence genes, resistance genes, prophages and plasmids

Escherichia coli O157:H7 is a foodborne pathogen that has been linked to global disease outbreaks. These diseases include hemorrhagic colitis and hemolytic uremic syndrome. It is vital to know the features that make this strain pathogenic to understand the development of disease outbreaks. In the current study, a comparative genomic analysis was carried out to determine the presence of structural and functional features of O157:H7 strains obtained from 115 National Center for Biotechnology Information database. These strains of interest were analysed in the following programs: BLAST Ring Image Generator, PlasmidFinder, ResFinder, VirulenceFinder, IslandViewer 4 and PHASTER. Five strains (ECP19-198, ECP19-798, F7508, F8952, H2495) demonstrated a great homology with Sakai because of a few regions missing. Five resistant genes were identified, however, Macrolide-associated resistance gene mdf(A) was commonly found in all genomes. Majority of the strains (97%) were positive for 15 of the virulent genes (espA, espB, espF, espJ, gad, chuA, eae, iss, nleA, nleB, nleC, ompT, tccP, terC and tir). The plasmid analysis demonstrated that the IncF group was the most prevalent in the strains analysed. The prophage and genomic island analysis showed a distribution of bacteriophages and genomic islands respectively. The results indicated that structural and functional features of the many O157:H7 strains differ and may be a result of obtaining mobile genetic elements via horizontal gene transfer. Understanding the evolution of O157:H7 strains pathogenicity in terms of their structural and functional features will enable the development of detection and control of transmission strategies.

Alteration of a Shiga toxin-encoding phage associated with a change in toxin production level and disease severity in Escherichia coli

Among the nine clades of Shiga toxin (Stx)-producing Escherichia coli O157:H7, clade 8 is thought to be highly pathogenic, as it causes severe disease more often than other clades. Two subclades have been proposed, but there are conflicting reports on intersubclade differences in Stx2 levels, although Stx2 production is a risk factor for severe disease development. The global population structure of clade 8 has also yet to be fully elucidated. Here, we present genome analyses of a global clade 8 strain set (n =510), including 147 Japanese strains sequenced in this study. The complete genome sequences of 18 of the 147 strains were determined to perform detailed clade-wide genome analyses together with 17 publicly available closed genomes. Intraclade variations in Stx2 production level and disease severity were also re-evaluated within the phylogenetic context. Based on phylogenomic analysis, clade 8 was divided into four lineages corresponding to the previously proposed SNP genotypes (SGs): SG8_30, SG8_31A, SG8_31B and SG8_32. SG8_30 and the common ancestor of the other SGs were first separated, with SG8_31A and SG8_31B emerging from the latter and SG8_32 emerging from SG8_31B. Comparison of 35 closed genomes revealed the overall structure of chromosomes and pO157 virulence plasmids and the prophage contents to be well conserved. However, Stx2a phages exhibit notable genomic diversity, even though all are integrated into the argW locus, indicating that subtype changes in Stx2a phage occurred from the γ subtype to its variant (γ_v1) in SG8_31A and from γ to δ in SG8_31B and SG8_32 via replacement of parts or almost entire phage genomes, respectively. We further show that SG8_30 strains (all carrying γ Stx2a phages) produce significantly higher levels of Stx2 and cause severe disease more frequently than SG8_32 strains (all carrying δ Stx2a phages). Clear conclusions on SG8_31A and SG8_31B cannot be made due to the small number of strains available, but as SG8_31A (carrying γ_v1 Stx2a phages) contains strains that produce much more Stx2 than SG8_30 strains, attention should also be paid to this SG.

Super Shedding in Enteric Pathogens: A Review

Super shedding occurs when a small number of individuals from a given host population shed high levels of a pathogen. Beyond this general definition, various interpretations of the shedding patterns have been proposed to identify super shedders, leading to the description of the super shedding phenomenon in a wide range of pathogens, in particular enteric pathogens, which are of considerable interest. Several underlying mechanisms may explain this observation, including factors related to the environment, the gut microbiota, the pathogen itself (i.e., genetic polymorphism), and the host (including immune factors). Moreover, data suggest that the interplay of these parameters, in particular at the host-pathogen-gut microbiota interface, is of crucial importance for the determination of the super shedding phenotype in enteric pathogens. As a phenomenon playing an important role in the epidemics of enteric diseases, the evidence of super shedding has highlighted the need to develop various control strategies.

Comparative Analysis of Three Bradyrhizobium diazoefficiens Genomes Show Specific Mutations Acquired during Selection for a Higher Motility Phenotype and Adaption to Laboratory Conditions

Microbial genomes are being extensively studied using next-generation sequencing technologies in order to understand the changes that occur under different selection regimes. In this work, the number and type of mutations that have occurred in three Bradyrhizobium diazoefficiens USDA 110^T strains under laboratory conditions and during selection for a more motile phenotypic variant were analyzed. Most of the mutations found in both processes consisted of single nucleotide polymorphisms, single nucleotide deletions or insertions. In the case of adaptation to laboratory conditions, half of the changes occurred within intergenic regions, and around 80% were insertions. When the more motile phenotypic variant was evaluated, eight single nucleotide polymorphisms and an 11-bp deletion were found, although none of them was directly related to known motility or chemotaxis genes. Two mutants were constructed to evaluate the 11-bp deletion affecting the alpha subunit of 2-oxoacid:acceptor oxidoreductase (AAV28_RS30705-blr6743). The results showed that this single deletion was not responsible for the enhanced motility phenotype. IMPORTANCE The genetic and genomic changes that occur under laboratory conditions in Bradyrhizobium diazoefficiens genomes remain poorly studied. Only a few genome sequences of this important nitrogen-fixing species are available, and there are no genome-wide comparative analyses of related strains. In the present work, we sequenced and compared the genomes of strains derived from a parent strain, B. diazoefficiens USDA 110, that has undergone processes of repeated culture in the laboratory environment, or phenotypic selection toward antibiotic resistance and enhanced motility. Our results represent the first analysis in B. diazoefficiens that provides insights into the specific mutations that are acquired during these processes.

EDL933 Strains of Escherichia coli O157 can Demonstrate Genetic Diversity and Differential Adherence to Bovine Recto-Anal Junction Squamous Epithelial Cells

Background:

Differences between Escherichia coli O157 (O157) strains are well-established with some of these strains being associated with major outbreaks in the US. EDL933 is one such O157 strain that caused a multistate outbreak in 1982 and has since been used as a prototype in various O157-related experiments.

Objective:

As O157 can readily acquire genetic mutations, we sought to determine if the genetic and phenotypic profiles of EDL933 strains from different sources would be consistent.

Methods:

We evaluated wild-type O157 strains stocked as EDL933 from three different laboratories, in the strain typing Polymorphic Amplified Typing Sequence (PATS) and the bovine rectal-anal junction squamous epithelial (RSE) cell- and HEp-2 cell- adherence assays. In addition, we also verified if Shiga toxins (Stx), the Locus of Enterocyte Effacement (LEE) or curli fimbriae contributed to the adherence phenotypes observed using mutant and wild-type EDL933 isolates.

Results:

Our results showed differences in PATS profiles and RSE cell-adherence phenotype, with no influence from the Stx or LEE genes, between EDL933 from different sources. Interestingly, the EDL933 strain that demonstrated the most contrasting diffuse adherence phenotype on RSE cells, EDL933-T, had decreased curli production that may have contributed to this phenotype.

Conclusion:

Our observations suggest that a comprehensive characterization of bacterial isolates, even if assigned to the same strain type prior to use in experiments, is warranted to ensure consistency and reproducibility of results.

Strain and host-cell dependent role of type-1 fimbriae in the adherence phenotype of super-shed Escherichia coli O157:H7

Super-shed (SS) Escherichia coli O157 (E. coli O157) demonstrate a strong, aggregative, locus of enterocyte effacement (LEE)-independent adherence phenotype on bovine recto-anal junction squamous epithelial (RSE) cells, and harbor polymorphisms in non-LEE-adherence-related loci, including in the type 1 fimbriae operon. To elucidate the role of type 1 fimbriae in strain- and host-specific adherence, we evaluated the entire Fim operon (FimB-H) and its adhesion (FimH) deletion mutants in four E. coli O157 strains, SS17, SS52, SS77 and EDL933, and evaluated the adherence phenotype in bovine RSE and human HEp-2 adherence assays. Consistent with the prevailing dogma that fimH expression is genetically switched off in E. coli O157, the ΔfimHSS52, ΔfimB-HSS52, ΔfimB-HSS17, and ΔfimHSS77 mutants remained unchanged in adherence phenotype to RSE cells. In contrast, the ΔfimHSS17 and ΔfimB-HSS77 mutants changed from a wild-type strong and aggregative, to a moderate and diffuse adherence phenotype, while both ΔfimHEDL933 and ΔfimB-HEDL933 mutants demonstrated enhanced binding to RSE cells (p < 0.05). Additionally, both ΔfimHSS17 and ΔfimHEDL933 were non-adherent to HEp-2 cells (p < 0.05). Complementation of the mutant strains with their respective wild-type genes restored parental phenotypes. Microscopy revealed that the SS17 and EDL933 strains indeed carry type 1 fimbriae-like structures shorter than those seen in uropathogenic E. coli. Taken together, these results provide compelling evidence for a strain and host cell type-dependent role of fimH and the fim operon in E. coli O157 adherence that needs to be further evaluated.

Farm animal contact is associated with progression to Hemolytic uremic syndrome in patients with Shiga toxin-producing Escherichia coli - Indiana, 2012-2018

Background

Hemolytic uremic syndrome (HUS) is a life-threatening complication of Shiga toxin-producing Escherichia coli (STEC) infection. The relationship between STEC exposure and severity of clinical outcomes is not well documented. We examined whether direct contact with farm animals increased the likelihood of HUS among Indiana residents diagnosed with STEC.

Methods

Exposure data for laboratory-confirmed STEC cases among Indiana residents during 2012–2018 were retrieved. Logistic regression and mediation analysis were performed to determine the extent to which a history of direct contact with farm animals was associated with post-diarrheal HUS independent of age and mediated by stx2 gene presence.

Results

A total of 784 STEC cases were retrieved. Of these, 46 (6%) developed HUS. Complete exposure data were available for 600 (77%) cases. A total of 24 (52%) HUS patients reported direct contact with farm animals, while 114 (21%) STEC patients who did not develop HUS reported this exposure. Among all STEC cases, HUS was associated with direct farm animal contact after adjusting for age (OR = 3.40, 95% CI: 1.81, 6.40). Detection of stx2 genes mediated 12% of the association between farm animal contact and HUS.

Conclusions

Direct farm animal contact was a risk factor for development of HUS among laboratory-confirmed STEC cases, independent of stx2 presence. Direct farm animal contact should be considered a potential predictor of progression to HUS when patients present for care and the mechanism for its effect on virulence investigated.

•

Independent of stx2 presence, contact with farm animals is a risk factor for the development of HUS.

•

Stx2 gene detection mediated 12.2% of the association between farm animal contact and HUS.

•

Exposure source may impact virulence of STEC and thus the severity of clinical outcomes.

Genetic and Functional Analyses of Virulence Potential of an Escherichia coli O157:H7 Strain Isolated From Super-Shedder Cattle

Shiga toxin (Stx)-producing Escherichia coli (STEC) O157:H7 is an enteric pathogen that causes life-threatening disease in humans, with cattle being major natural reservoirs. A group of STEC O157:H7 with a dramatic combination of high virulence potentials and super-shedder bovine origin have been isolated. Here, an STEC O157:H7 isolate, JEONG-1266, was analyzed by comparative genomics, stx genotyping, and phenotypic analyses. The phylogenetic typing and whole-genome comparison consistently showed that JEONG-1266 is genetically close to EC4115 (one of 2006 Spinach outbreak isolates) and SS17 (an isolate from super-shedder cattle) strains, all of which belong to lineage I/II and Clade 8. Both lineage I/II and Clade 8 are known to be mostly associated with clinical strains with high virulence and severe clinical symptoms. Further, JEONG-1266, like EC4115 and SS17, harbors stx2a/stx2c genes, and carries Stx-encoding prophages, specifically the φstx2a-γ subtype. Possession of the φstx2a-γ subtype of Stx-encoding prophages and production of Stx2a have been shown to be a key signature associated with hypervirulent STEC O157:H7 strains. In silico virulence typing elucidated JEONG-1266, EC4115, and SS17 shared a highly conserved profile of key virulence genes at the nucleotide sequence level. Consistently, phenotypic data showed that JEONG-1266 expressed a high level of Stx2 toxins and had the full capacity of adhesion in vitro. Taken together, our study suggests that JEONG-1266 may represent an emerging STEC O157:H7 group, which are hypervirulent strains that originate from super-shedders, that can be a threat to food safety and public health.

Is Shiga Toxin-Producing Escherichia coli O45 No Longer a Food Safety Threat? The Danger is Still Out There

Many Shiga toxin-producing Escherichia coli (STEC) strains, including the serogroups of O157 and most of the top six non-O157 serotypes, are frequently associated with foodborne outbreaks. Therefore, they have been extensively studied using next-generation sequencing technology. However, related information regarding STEC O45 strains is scarce. In this study, three environmental E. coli O45:H16 strains (RM11911, RM13745, and RM13752) and one clinical E. coli O45:H2 strain (SJ7) were sequenced and used to characterize virulence factors using two reference E. coli O45:H2 strains of clinical origin. Subsequently, whole-genome-based phylogenetic analysis was conducted for the six STEC O45 strains and nine other reference STEC genomes, in order to evaluate their evolutionary relationship. The results show that one locus of enterocyte effacement pathogenicity island was found in all three STEC O45:H2 strains, but not in the STEC O45:H16 strains. Additionally, E. coli O45:H2 strains were evolutionarily close to E. coli O103:H2 strains, sharing high homology in terms of virulence factors, such as Stx prophages, but were distinct from E. coli O45:H16 strains. The findings show that E. coli O45:H2 may be as virulent as E. coli O103:H2, which is frequently associated with severe illness and can provide genomic evidence to facilitate STEC surveillance.

Supershed Escherichia coli O157:H7 Has Potential for Increased Persistence on the Rectoanal Junction Squamous Epithelial Cells and Antibiotic Resistance

Supershedding cattle shed Escherichia coli O157:H7 (O157) at ≥ 104 colony-forming units/g feces. We recently demonstrated that a supershed O157 (SS-O157) strain, SS-17, hyperadheres to the rectoanal junction (RAJ) squamous epithelial (RSE) cells which may contribute to SS-O157 persistence at this site in greater numbers, thereby increasing the fecal O157 load characterizing the supershedding phenomenon. In order to verify if this would be the signature adherence profile of any SS-O157, we tested additional SS-O157 isolates (n = 101; each from a different animal) in the RSE cell adherence assay. Similar to SS-17, all 101 SS-O157 exhibited aggregative adherence on RSE cells, with 56% attaching strongly (>10 bacteria/cell; hyperadherent) and 44% attaching moderately (1-10 bacteria/cells). Strain typing using Polymorphic Amplified Typing Sequences (PATS) analysis assigned the 101 SS-O157 into 5 major clades but not to any predominant genotype. Interestingly, 69% of SS-O157 isolates were identical to human O157 outbreak strains based on pulsed field gel electrophoresis profiles (CDC PulseNet Database), grouped into two clades by PATS distinguishing them from remaining SS-O157, and were hyperadherent on RSE cells. A subset of SS-O157 isolates (n = 53) representing different PATS and RSE cell adherence profiles were analyzed for antibiotic resistance (AR). Several SS-O157 (30/53) showed resistance to sulfisoxazole, and one isolate was resistant to both sulfisoxazole and tetracycline. Minimum inhibitory concentration (MIC) tests confirmed some of the resistance observed using the Kirby-Bauer disk diffusion test. Each SS-O157 isolate carried at least 10 genes associated with AR. However, genes directly associated with AR were rarely amplified: aac (3)-IV in 2 isolates, sul2 in 3 isolates, and tetB in one isolate. The integrase gene, int, linked with integron-based AR acquisition/transmission, was identified in 92% of SS-O157 isolates. Our results indicate that SS-O157 isolates could potentially persist longer at the bovine RAJ but exhibit limited resistance towards clinical antibiotics.

Construction of a Lectin-Glycan Interaction Network from Enterohemorrhagic Escherichia coli Strains by Multi-omics Analysis

Enterohemorrhagic Escherichia coli (EHEC) causes hemorrhagic colitis and hemolytic uremic syndrome. EHEC infection begins with bacterial adherence to the host intestine via lectin-like adhesins that bind to the intestinal wall. However, EHEC-related lectin–glycan interactions (LGIs) remain unknown. Here, we conducted a genome-wide investigation of putative adhesins to construct an LGI network. We performed microarray-based transcriptomic and proteomic analyses with E. coli EDL933. Using PSORTb-based analysis, potential outer-membrane-embedded adhesins were predicted from the annotated genes of 318 strains. Predicted proteins were classified using TMHMM v2.0, SignalP v5.0, and LipoP v1.0. Functional and protein–protein interaction analyses were performed using InterProScan and String databases, respectively. Structural information of lectin candidate proteins was predicted using Iterative Threading ASSEmbly Refinement (I-TASSER) and Spatial Epitope Prediction of Protein Antigens (SEPPA) tools based on 3D structure and B-cell epitopes. Pathway analysis returned 42,227 Gene Ontology terms; we then selected 2585 lectin candidate proteins by multi-omics analysis and performed homology modeling and B-cell epitope analysis. We predicted a total of 24,400 outer-membrane-embedded proteins from the genome of 318 strains and integrated multi-omics information into the genomic information of the proteins. Our integrated multi-omics data will provide a useful resource for the construction of LGI networks of E. coli.

Precision public health to inhibit the contagion of disease and move toward a future in which microbes spread health

Antimicrobial resistance continues to outpace the development of new chemotherapeutics. Novel pathogens continue to evolve and emerge. Public health innovation has the potential to open a new front in the war of "our wits against their genes" (Joshua Lederberg). Dense sampling coupled to next generation sequencing can increase the spatial and temporal resolution of microbial characterization while sensor technologies precisely map physical parameters relevant to microbial survival and spread. Microbial, physical, and epidemiological big data could be combined to improve prospective risk identification. However, applied in the wrong way, these approaches may not realize their maximum potential benefits and could even do harm. Minimizing microbial-human interactions would be a mistake. There is evidence that microbes previously thought of at best "benign" may actually enhance human health. Benign and health-promoting microbiomes may, or may not, spread via mechanisms similar to pathogens. Infectious vaccines are approaching readiness to make enhanced contributions to herd immunity. The rigorously defined nature of infectious vaccines contrasts with indigenous "benign or health-promoting microbiomes" but they may converge. A "microbial Neolithic revolution" is a possible future in which human microbial-associations are understood and managed analogously to the macro-agriculture of plants and animals. Tradeoffs need to be framed in order to understand health-promoting potentials of benign, and/or health-promoting microbiomes and infectious vaccines while also discouraging pathogens. Super-spreaders are currently defined as individuals who play an outsized role in the contagion of infectious disease. A key unanswered question is whether the super-spreader concept may apply similarly to health-promoting microbes. The complex interactions of individual rights, community health, pathogen contagion, the spread of benign, and of health-promoting microbiomes including infectious vaccines require study. Advancing the detailed understanding of heterogeneity in microbial spread is very likely to yield important insights relevant to public health.