Complete Genome Sequences of Two Escherichia coli O145:H28 Outbreak Strains of Food Origin

Genomic Characterization of Two Shiga Toxin-Converting Bacteriophages Induced From Environmental Shiga Toxin-Producing Escherichia coli

Shiga toxin (Stx), encoded by stx genes located in prophage sequences, is the major agent responsible for the pathogenicity of Shiga toxin-producing Escherichia coli (STEC) and is closely associated with the development of hemolytic uremic syndrome (HUS). Although numerous Stx prophage sequences have been reported as part of STEC bacterial genomes, the information about the genomic characterization of Stx-converting bacteriophages induced from STEC strains is relatively scarce. The objectives of this study were to genomically characterize two Stx-converting phages induced from environmental STEC strains and to evaluate their correlations with published Stx-converting phages and STEC strains of different origins. The Stx1-converting phage Lys8385Vzw and the Stx2-converting phage Lys19259Vzw were induced from E. coli O103:H11 (RM8385) and E. coli O157:H7 (RM19259), respectively. Whole-genome sequencing of these phages was conducted on a MiSeq sequencer for genomic characterization. Phylogenetic analysis and comparative genomics were performed to determine the correlations between these two Stx-converting phages, 13 reference Stx-converting phages, and 10 reference STEC genomes carrying closely related Stx prophages. Both Stx-converting phages Lys8385Vzw and Lys19259Vzw had double-stranded DNA, with genome sizes of 50,953 and 61,072 bp, respectively. Approximately 40% of the annotated coding DNA sequences with the predicted functions were likely associated with the fitness for both phages and their bacterial hosts. The whole-genome–based phylogenetic analysis of these two Stx-converting phages and 13 reference Stx-converting phages revealed that the 15 Stx-converting phages were divided into three distinct clusters, and those from E. coli O157:H7, in particular, were distributed in each cluster, demonstrating the high genomic diversity of these Stx-converting phages. The genomes of Stx-converting phage Lys8385Vzw and Lys19259Vzw shared a high-nucleotide similarity with the prophage sequences of the selected STEC isolates from the clinical and environmental origin. The findings demonstrate the genomic diversity of Stx-converting phages induced from different STEC strains and provide valuable insights into the dissemination of stx genes among E. coli population via the lysogenization of Stx-converting phages.

Comparative genomic analysis of a Shiga toxin-producing Escherichia coli (STEC) O145:H25 associated with a severe pediatric case of hemolytic uremic syndrome in Davidson County, Tennessee, US

Background

Shiga toxin-producing E. coli (STECs) are foodborne pathogens associated with bloody diarrhea and hemolytic uremic syndrome (HUS). Although the STEC O157 serogroup accounts for the highest number of infections, HUS-related complications and deaths, the STEC non-O157, as a group, accounts for a larger proportion of STEC infections and lower HUS cases. There is limited information available on how to recognize non-O157 serotypes associated with severe disease. The objectives of this study were to describe a patient with STEC non-O157 infection complicated with HUS and to conduct a comparative whole genome sequence (WGS) analysis among the patient’s STEC clinical isolate and STEC O157 and non-O157 strains.

Results

The STEC O145:H25 strain EN1I-0044-2 was isolated from a pediatric patient with diarrhea, HUS and severe neurologic and cardiorespiratory complications, who was enrolled in a previously reported case-control study of acute gastroenteritis conducted in Davidson County, Tennessee in 2013. The strain EN1I-0044-2 genome sequence contained a chromosome and three plasmids. Two of the plasmids were similar to those present in O145:H25 strains whereas the third unique plasmid EN1I-0044-2_03 shared no similarity with other STEC plasmids, and it carried 23 genes of unknown function. Strain EN1I-0044-2, compared with O145:H25 and O157 serogroup strains shared chromosome- and plasmid-encoded virulence factors, including Shiga toxin, LEE type III secretion system, LEE effectors, SFP fimbriae, and additional toxins and colonization factors.

Conclusions

A STEC O145:H25 strain EN1I-0044-2 was isolated from a pediatric patient with severe disease, including HUS, in Davidson County, TN. Phylogenetic and comparison WGS analysis provided evidence that strain EN1I-0044-2 closely resembles O145:H25, and confirmed an independent evolutionary path of STEC O145:H25 and O145:H28 serotypes. The strain EN1I-0044-2 virulence make up was similar to other O145:H25 and O157 serogroups. It carried stx2 and the LEE pathogenicity island, and additional colonization factors and enterotoxin genes. A unique feature of strain EN1I-0044-2 was the presence of plasmid pEN1I-0044-2_03 carrying genes with functions to be determined. Further studies will be necessary to elucidate the role that newly acquired genes by O145:H25 strains play in pathogenesis, and to determine if they may serve as genetic markers of severe disease.

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.

Prediction, Diversity, and Genomic Analysis of Temperate Phages Induced From Shiga Toxin-Producing Escherichia coli Strains

Shiga toxin-producing Escherichia coli (STEC) is a notorious foodborne pathogen containing stx genes located in the sequence region of Shiga toxin (Stx) prophages. Stx prophages, as one of the mobile elements, are involved in the transfer of virulence genes to other strains. However, little is known about the diversity of prophages among STEC strains. The objectives of this study were to predict various prophages from different STEC genomes and to evaluate the effect of different stress factors on Stx prophage induction. Forty bacterial whole-genome sequences of STEC strains obtained from National Center for Biotechnology Information (NCBI) were used for the prophage prediction using PHASTER webserver. Eight of the STEC strains from different serotypes were subsequently selected to quantify the induction of Stx prophages by various treatments, including antibiotics, temperature, irradiation, and antimicrobial agents. After induction, Stx1-converting phage Lys8385Vzw and Stx2-converting phage Lys12581Vzw were isolated and further confirmed for the presence of stx genes using conventional PCR. Phage morphology was observed by transmission electron microscopy. The prediction results showed an average of 8–22 prophages, with one or more encoding stx, were predicted from each STEC genome obtained in this study. Additionally, the phylogenetic analysis revealed high genetic diversity of Stx prophages among the 40 STEC genomes. However, the sequences of Stx prophages in the genomes of STEC O45, O111, and O121 strains, in general, shared higher genetic homology than those in other serotypes. Interestingly, most STEC strains with two or more stx genes carried at least one each of Stx1 and Stx2 prophages. The induction results indicated EDTA and UV were the most effective inducers of Stx1 and Stx2 prophages of the 8 selected STECs, respectively. Additionally, both Stx-converting phages could infect non-pathogenic E. coli (WG5, DH5α, and MG1655) and form new lysogens. The findings of this study confirm that Stx prophages can be induced by environmental stress, such as exposure to solar radiation, and lysogenize other commensal E. coli strains.

Escherichia coli ST302: Genomic Analysis of Virulence Potential and Antimicrobial Resistance Mediated by Mobile Genetic Elements

aEPEC are associated with persistent diarrhea, and diarrheal outbreaks in both humans and animals worldwide. They are differentiated from typical EPEC by the lack of bundle-forming pili, and from EHEC by the lack of phage-mediated stx toxins. However, phylogenetic analyses often associate aEPEC with EHEC, promoting the hypothesis that aEPEC are the progenitors of EHEC, which is supported by aEPEC conversion to EHEC by stx-carrying phages. While aEPEC can cause disease outright, the potential to acquire stx, one of the most potent bacterial toxins known, merits close monitoring. Escherichia coli ST302 (O108:H9, O182:H9, O45:H9) are aEPEC that have been isolated from diarrheic human, pig and rabbit hosts, as well as in healthy pigs, however, no study to date has focused on E. coli ST302 strains. Through WGS and hybrid assembly we present the first closed chromosome, and two circularized plasmids of an ST302 strain - F2_18C, isolated from a healthy pig in Australia. A phylogenetic analysis placed E. coli ST302 strains in proximity to EHEC ST32 (O145:H28) strains. Public databases were interrogated for WGSs of E. coli ST302 strains and short-read gene screens were used to compare their virulence-associated gene (VAG) and antimicrobial resistance gene (ARG) cargo. E. coli ST302 strains carry diverse VAGs, including those that typically associated with extraintestinal pathogenic E. coli (ExPEC). Plasmid comparisons showed that pF2_18C_FIB shared homology with EHEC virulence plasmids such as pO103 while pF2_18C_HI2 is a large multidrug resistance IncHI2:ST3 plasmid. A comparison of 33 HI2:ST3 plasmids demonstrated that those of Australian origin have not acquired resistances to extended-spectrum beta-lactams, colistin, fosfomycin or rifampicin, unlike those originating from Asia. F2_18C was shown to carry two additional pathogenicity islands – ETT2, and the STEC-associated PAI_CL3, plasmid-associated heavy metal resistance genes, as well as several unoccupied stx-phage attachment sites. This study sheds light on the virulence and AMR potential of E. coli ST302 strains and informs AMR genomic surveillance.

Complete Genome Sequence of a Shiga Toxin-Converting Bacteriophage, Escherichia Phage Lys12581Vzw, Induced from an Outbreak Shiga Toxin-Producing Escherichia coli Strain

Although numerous Shiga toxin (Stx)-producing Escherichia coli (STEC) strains have been sequenced, genomic information on Stx-converting phages, highly related to the primary virulence factors of STEC, is scarce. Here, we report the complete genome sequence of a Stx-converting phage induced from an outbreak STEC O145 strain.

Although numerous Shiga toxin (Stx)-producing Escherichia coli (STEC) strains have been sequenced, genomic information on Stx-converting phages, highly related to the primary virulence factors of STEC, is scarce. Here, we report the complete genome sequence of a Stx-converting phage induced from an outbreak STEC O145 strain.

Evolution and typing of IncC plasmids contributing to antibiotic resistance in Gram-negative bacteria

The large, broad host range IncC plasmids are important contributors to the spread of key antibiotic resistance genes and over 200 complete sequences of IncC plasmids have been reported. To track the spread of these plasmids accurate typing to identify the closest relatives is needed. However, typing can be complicated by the high variability in resistance gene content and various typing methods that rely on features of the conserved backbone have been developed. Plasmids can be broadly typed into two groups, type 1 and type 2, using four features that differentiate the otherwise closely related backbones. These types are found in many different countries in bacteria from humans and animals. However, hybrids of type 1 and type 2 are also occasionally seen, and two further types, each represented by a single plasmid, were distinguished. Generally, the antibiotic resistance genes are located within a small number of resistance islands, only one of which, ARI-B, is found in both type 1 and type 2. The introduction of each resistance island generates a new lineage and, though they are continuously evolving via the loss of resistance genes or introduction of new ones, the island positions serve as valuable lineage-specific markers. A current type 2 lineage of plasmids is derived from an early type 2 plasmid but the sequences of early type 1 plasmids include features not seen in more recent type 1 plasmids, indicating a shared ancestor rather than a direct lineal relationship. Some features, including ones essential for maintenance or for conjugation, have been examined experimentally.

Complete Genome Sequences of Two Atypical Enteropathogenic Escherichia coli O145 Environmental Strains

Escherichia coli O145 strains RM14715 and RM14723 were isolated from wildlife feces near a leafy greens-growing region in Yuma, Arizona. Both strains carry a distinct genotype compared with the E. coli O145 strains isolated from Salinas Valley, California. Here we report complete genome sequences and annotations of RM14715 and RM14723.

Complete Genome Sequence of a Natural Escherichia coli O145:H11 Isolate That Belongs to Phylogroup A

Escherichia coli O145:H11 strain RM14721 was originally isolated from wildlife feces near a leafy greens-growing region in Yuma, AZ. This strain was initially positive for stx₁; however, in subsequent cultures, stx₁ was not detected by PCR. Here, we report the complete genome sequence and annotation of RM14721.

Conditional Function of Autoaggregative Protein Cah and Common cah Mutations in Shiga Toxin-Producing Escherichia coli

Cah is a calcium-binding autotransporter protein involved in autoaggregation and biofilm formation. Although cah is widespread in Shiga toxin-producing Escherichia coli (STEC), we detected mutations in cah at a frequency of 31.3% in this pathogen. In STEC O157:H7 supershedder strain SS17, a large deletion results in a smaller coding sequence, encoding a protein lacking the C-terminal 71 amino acids compared with Cah in STEC O157:H7 strain EDL933. We examined the function of Cah in biofilm formation and host colonization to better understand the selective pressures for cah mutations. EDL933-Cah played a conditional role in biofilm formation in vitro: it enhanced E. coli DH5α biofilm formation on glass surfaces under agitated culture conditions that prevented autoaggregation but inhibited biofilm formation under hydrostatic conditions that facilitated autoaggregation. This function appeared to be strain dependent since Cah-mediated biofilm formation was diminished when an EDL933 cah gene was expressed in SS17. Deletion of cah in EDL933 enhanced bacterial attachment to spinach leaves and altered the adherence pattern of EDL933 to bovine recto-anal junction squamous epithelial (RSE) cells. In contrast, in trans expression of EDL933 cah in SS17 increased its attachment to leaf surfaces, and in DH5α, it enhanced its adherence to RSE cells. Hence, the ecological function of Cah appears to be modulated by environmental conditions and other bacterial strain-specific properties. Considering the prevalence of cah in STEC and its role in attachment and biofilm formation, cah mutations might be selected in ecological niches in which inactivation of Cah would result in an increased fitness in STEC during colonization of plants or animal hosts.IMPORTANCE Shiga toxin-producing Escherichia coli (STEC) harbors genes encoding diverse adhesins, and many of these are known to play an important role in bacterial attachment and host colonization. We demonstrated here that the autotransporter protein Cah confers on E. coli DH5α cells a strong autoaggregative phenotype that is inversely correlated with its ability to form biofilms and plays a strain-specific role in plant and animal colonization by STEC. Although cah is widespread in the STEC population, we detected a mutation rate of 31.3% in cah, which is similar to that reported for rpoS and fimH The formation of cell aggregates due to increased bacterium-to-bacterium interactions may be disadvantageous to bacterial populations under conditions that favor a planktonic state in STEC. Therefore, a loss-of-function mutation in cah is likely a selective trait in STEC when autoaggregative properties become detrimental to bacterial cells and may contribute to the adaptability of STEC to fluctuating environments.

The Mobilome; A Major Contributor to Escherichia coli stx2-Positive O26:H11 Strains Intra-Serotype Diversity

Shiga toxin-producing Escherichia coli of serotype O26:H11/H- constitute a diverse group of strains and several clones with distinct genetic characteristics have been identified and characterized. Whole genome sequencing was performed using Illumina and PacBio technologies on eight stx2-positive O26:H11 strains circulating in France. Comparative analyses of the whole genome of the stx2-positive O26:H11 strains indicate that several clones of EHEC O26:H11 are co-circulating in France. Phylogenetic analysis of the French strains together with stx2-positive and stx-negative E. coli O26:H11 genomes obtained from Genbank indicates the existence of four clonal complexes (SNP-CCs) separated in two distinct lineages, one of which comprises the "new French clone" (SNP-CC1) that appears genetically closely related to stx-negative attaching and effacing E. coli (AEEC) strains. Interestingly, the whole genome SNP (wgSNP) phylogeny is summarized in the cas gene phylogeny, and a simple qPCR assay targeting the CRISPR array specific to SNP-CC1 (SP_O26-E) can distinguish between the two main lineages. The PacBio sequencing allowed a detailed analysis of the mobile genetic elements (MGEs) of the strains. Numerous MGEs were identified in each strain, including a large number of prophages and up to four large plasmids, representing overall 8.7-19.8% of the total genome size. Analysis of the prophage pool of the strains shows a considerable diversity with a complex history of recombination. Each clonal complex (SNP-CC) is characterized by a unique set of plasmids and phages, including stx-prophages, suggesting evolution through separate acquisition events. Overall, the MGEs appear to play a major role in O26:H11 intra-serotype clonal diversification.

Short-term evolution of Shiga toxin-producing Escherichia coli O157:H7 between two food-borne outbreaks

Shiga toxin-producing Escherichia coli (STEC) O157:H7 is a public health threat and outbreaks occur worldwide. Here, we investigate genomic differences between related STEC O157:H7 that caused two outbreaks, eight weeks apart, at the same restaurant. Short-read genome sequencing divided the outbreak strains into two sub-clusters separated by only three single-nucleotide polymorphisms in the core genome while traditional typing identified them as separate phage types, PT8 and PT54. Isolates did not cluster with local strains but with those associated with foreign travel to the Middle East/North Africa. Combined long-read sequencing approaches and optical mapping revealed that the two outbreak strains had undergone significant microevolution in the accessory genome with prophage gain, loss and recombination. In addition, the PT54 sub-type had acquired a 240 kbp multi-drug resistance (MDR) IncHI2 plasmid responsible for the phage type switch. A PT54 isolate had a general fitness advantage over a PT8 isolate in rich medium, including an increased capacity to use specific amino acids and dipeptides as a nitrogen source. The second outbreak was considerably larger and there were multiple secondary cases indicative of effective human-to-human transmission. We speculate that MDR plasmid acquisition and prophage changes have adapted the PT54 strain for human infection and transmission. Our study shows the added insights provided by combining whole-genome sequencing approaches for outbreak investigations.

An Environmental Shiga Toxin-Producing Escherichia coli O145 Clonal Population Exhibits High-Level Phenotypic Variation That Includes Virulence Traits

Shiga toxin-producing Escherichia coli (STEC) serotype O145 is one of the major non-O157 serotypes associated with severe human disease. Here we examined the genetic diversity, population structure, virulence potential, and antimicrobial resistance profiles of environmental O145 strains recovered from a major produce production region in California. Multilocus sequence typing analyses revealed that sequence type 78 (ST-78), a common ST in clinical strains, was the predominant genotype among the environmental strains. Similarly, all California environmental strains belonged to H28, a common H serotype in clinical strains. Although most environmental strains carried an intact fliC gene, only one strain retained swimming motility. Diverse stx subtypes were identified, including stx1a, stx2a, stx2c, and stx2e. Although no correlation was detected between the stx genotype and Stx1 production, high Stx2 production was detected mainly in strains carrying stx2a only and was correlated positively with the cytotoxicity of Shiga toxin. All environmental strains were capable of producing enterohemolysin, whereas only 10 strains were positive for anaerobic hemolytic activity. Multidrug resistance appeared to be common, as nearly half of the tested O145 strains displayed resistance to at least two different classes of antibiotics. The core virulence determinants of enterohemorrhagic E. coli were conserved in the environmental STEC O145 strains; however, there was large variation in the expression of virulence traits among the strains that were highly related genotypically, implying a trend of clonal divergence. Several cattle isolates exhibited key virulence traits comparable to those of the STEC O145 outbreak strains, emphasizing the emergence of hypervirulent strains in agricultural environments.

Sequence Variations in the Flagellar Antigen Genes fliCH25 and fliCH28 of Escherichia coli and Their Use in Identification and Characterization of Enterohemorrhagic E. coli (EHEC) O145:H25 and O145:H28

Enterohemorrhagic E. coli (EHEC) serogroup O145 is regarded as one of the major EHEC serogroups involved in severe infections in humans. EHEC O145 encompasses motile and non-motile strains of serotypes O145:H25 and O145:H28. Sequencing the fliC-genes associated with the flagellar antigens H25 and H28 revealed the genetic diversity of the fliCH25 and fliCH28 gene sequences in E. coli. Based on allele discrimination of these fliC-genes real-time PCR tests were designed for identification of EHEC O145:H25 and O145:H28. The fliCH25 genes present in O145:H25 were found to be very similar to those present in E. coli serogroups O2, O100, O165, O172 and O177 pointing to their common evolution but were different from fliCH25 genes of a multiple number of other E. coli serotypes. In a similar way, EHEC O145:H28 harbor a characteristic fliCH28 allele which, apart from EHEC O145:H28, was only found in enteropathogenic (EPEC) O28:H28 strains that shared some common traits with EHEC O145:H28. The real time PCR-assays targeting these fliCH25[O145] and fliCH28[O145] alleles allow better characterization of EHEC O145:H25 and EHEC O145:H28. Evaluation of these PCR assays in spiked ready-to eat salad samples resulted in specific detection of both types of EHEC O145 strains even when low spiking levels of 1-10 cfu/g were used. Furthermore these PCR assays allowed identification of non-motile E. coli strains which are serologically not typable for their H-antigens. The combined use of O-antigen genotyping (O145wzy) and detection of the respective fliCH25[O145] and fliCH28[O145] allele types contributes to improve identification and molecular serotyping of E. coli O145 isolates.

The master activator of IncA/C conjugative plasmids stimulates genomic islands and multidrug resistance dissemination

Dissemination of antibiotic resistance genes occurs mostly by conjugation, which mediates DNA transfer between cells in direct contact. Conjugative plasmids of the IncA/C incompatibility group have become a substantial threat due to their broad host-range, the extended spectrum of antimicrobial resistance they confer, their prevalence in enteric bacteria and their very efficient spread by conjugation. However, their biology remains largely unexplored. Using the IncA/C conjugative plasmid pVCR94ΔX as a prototype, we have investigated the regulatory circuitry that governs IncA/C plasmids dissemination and found that the transcriptional activator complex AcaCD is essential for the expression of plasmid transfer genes. Using chromatin immunoprecipitation coupled with exonuclease digestion (ChIP-exo) and RNA sequencing (RNA-seq) approaches, we have identified the sequences recognized by AcaCD and characterized the AcaCD regulon. Data mining using the DNA motif recognized by AcaCD revealed potential AcaCD-binding sites upstream of genes involved in the intracellular mobility functions (recombination directionality factor and mobilization genes) in two widespread classes of genomic islands (GIs) phylogenetically unrelated to IncA/C plasmids. The first class, SGI1, confers and propagates multidrug resistance in Salmonella enterica and Proteus mirabilis, whereas MGIVmi1 in Vibrio mimicus belongs to a previously uncharacterized class of GIs. We have demonstrated that through expression of AcaCD, IncA/C plasmids specifically trigger the excision and mobilization of the GIs at high frequencies. This study provides new evidence of the considerable impact of IncA/C plasmids on bacterial genome plasticity through their own mobility and the mobilization of genomic islands.

Multidrug resistance is a major health concern that complicates treatments of even the most common infections caused by bacteria. In recent years, IncA/C plasmids have emerged and spread in bacteria infecting humans, food-producing animals and food products, driving at the same time the dissemination of a broad spectrum of antibiotic resistance genes in environmental and in clinical settings. In this study, we have characterized the regulatory pathway that governs IncA/C plasmid dissemination. We have found that AcaCD, the master activator complex encoded by these plasmids, is not only essential for the dissemination of IncA/C plasmids but also activates unrelated mobile genetic elements in bacterial genomes, thereby further promoting the interspecies propagation of multidrug resistance and other adaptive traits at a very high frequency.