Comparison of whole genome sequences from human and non-human Escherichia coli O26 strains

Predicting Host Association for Shiga Toxin-Producing E. coli Serogroups by Machine Learning

Escherichia coli is a species of bacteria that can be present in a wide variety of mammalian hosts and potentially soil environments. E. coli has an open genome and can show considerable diversity in gene content between isolates. It is a reasonable assumption that gene content reflects evolution of strains in particular host environments and therefore can be used to predict the host most likely to be the source of an isolate. An extrapolation of this argument is that strains may also have gene content that favors success in multiple hosts and so the possibility of successful transmission from one host to another, for example, from cattle to human, can also be predicted based on gene content. In this methods chapter, we consider the issue of Shiga toxin (Stx)-producing E. coli (STEC) strains that are present in ruminants as the main host reservoir and for which we know that a subset causes life-threatening infections in humans. We show how the genome sequences of E. coli isolated from both cattle and humans can be used to build a classifier to predict human and cattle host association and how this can be applied to score key STEC serotypes known to be associated with human infection. With the example dataset used, serogroups O157, O26, and O111 show the highest, and O103 and O145 the lowest, predictions for human association. The long-term ambition is to combine such machine learning predictions with phylogeny to predict the zoonotic threat of an isolate based on its whole genome sequence (WGS).

Multiplex PCR Assays for the Detection of One Hundred and Thirty Seven Serogroups of Shiga Toxin-Producing Escherichia coli Associated With Cattle

Escherichia coli carrying prophage with genes that encode for Shiga toxins are categorized as Shiga toxin-producing E. coli (STEC) pathotype. Illnesses caused by STEC in humans, which are often foodborne, range from mild to bloody diarrhea with life-threatening complications of renal failure and hemolytic uremic syndrome and even death, particularly in children. As many as 158 of the total 187 serogroups of E. coli are known to carry Shiga toxin genes, which makes STEC a major pathotype of E. coli. Seven STEC serogroups, called top-7, which include O26, O45, O103, O111, O121, O145, and O157, are responsible for the majority of the STEC-associated human illnesses. The STEC serogroups, other than the top-7, called “non-top-7” have also been associated with human illnesses, more often as sporadic infections. Ruminants, particularly cattle, are principal reservoirs of STEC and harbor the organisms in the hindgut and shed in the feces, which serves as a major source of food and water contaminations. A number of studies have reported on the fecal prevalence of top-7 STEC in cattle feces. However, there is paucity of data on the prevalence of non-top-7 STEC serogroups in cattle feces, generally because of lack of validated detection methods. The objective of our study was to develop and validate 14 sets of multiplex PCR (mPCR) assays targeting serogroup-specific genes to detect 137 non-top-7 STEC serogroups previously reported to be present in cattle feces. Each assay included 7–12 serogroups and primers were designed to amplify the target genes with distinct amplicon sizes for each serogroup that can be readily identified within each assay. The assays were validated with 460 strains of known serogroups. The multiplex PCR assays designed in our study can be readily adapted by most laboratories for rapid identification of strains belonging to the non-top-7 STEC serogroups associated with cattle.

Prevalence of Shiga toxin-producing Escherichia coli in pasture-based dairy herds

Shiga toxin-producing Escherichia coli strains (STEC) are food-borne pathogens. While E. coli O157:H7 is commonly associated with cattle, less is known about the prevalence of non-O157 STEC serogroups in bovines. This study evaluated the prevalence and virulence status of O157:H7 and six E. coli O-serogroups (O26, O103, O45, O145, O121, O111) in New Zealand dairy farms using molecular as well as culture-based methods. Fresh farm dairy effluent (FDE) (n = 36) and composite calf faeces (n = 12) were collected over three samplings from 12 dairy farms. All seven target serogroups were detected through molecular techniques. Of the 202 isolates which were serologically confirmed following traditional culturing and immunomagnetic separation (IMS), O103, O26, O45 and O121 were the most common serogroups, being found in 81, 47, 42 and 32% of the FDE and in 17, 33, 25 and 9% of the calf faeces respectively. The majority (157/202) of the isolates were negative for stx and eae virulence genes. The prevalence of the seven target STEC was low, and only nine O26 isolates (4%) were recovered from four of the farms. The study has highlighted the need for improving the isolation of Top 7 STEC from the stx-negative populations present in fresh dairy effluent and calf faeces. SIGNIFICANCE AND IMPACT OF THE STUDY: Shiga toxin-producing Escherichia coli (STEC) are important food-borne pathogens that can cause severe illness in humans. Cattle are asymptomatic reservoirs for STEC, and transmission to humans can be by consumption of food products or water contaminated with cattle faeces. Our study investigated the prevalence of O157:H7 and six E. coli serogroups of STEC (O26, O103, O45, O145, O121, O111) over time in the dairy reservoir and increases the knowledge and understanding of these pathogens on pasture-based farms. Such information is required to develop risk-assessment models aiming at limiting transmission of these STEC to human.

Molecular serogrouping ofEscherichia coli

O-antigens present on the surface ofEscherichia coliprovide antigenic specificity for the strain and are the main components for O-serogroup designation. Serotyping using O-group-specific antisera for the identification ofE. coliO-serogroups has been traditionally the gold-standard for distinguishingE. colistrains. Knowledge of the O-group is important for determining pathogenic lineage, classifyingE. colifor epidemiological studies, for determining virulence, and for tracing outbreaks of diseases and sources of infection. However, serotyping has limitations, as the antisera generated against each specific O-group may cross-react, many strains are non-typeable, and others can autoagglutinate or be rough (lacking an O-antigen). Currently, the nucleotide sequences are available for most of the 187 designatedE. coliO-groups. Public health and other laboratories are considering whole genome sequencing to develop genotypic methods to determine O-groups. These procedures require instrumentation and analysis that may not be accessible and may be cost-prohibitive at this time. In this review, we have identified unique gene sequences within the O-antigen gene clusters and have targeted these genes for identification of O-groups using the polymerase chain reaction. This information can be used to distinguish O-groups by developing other platforms forE. colidiagnostics in the future.

Comparative genomics reveals differences in mobile virulence genes of Escherichia coli O103 pathotypes of bovine fecal origin

Escherichia coli O103, harbored in the hindgut and shed in the feces of cattle, can be enterohemorrhagic (EHEC), enteropathogenic (EPEC), or putative non-pathotype. The genetic diversity particularly that of virulence gene profiles within O103 serogroup is likely to be broad, considering the wide range in severity of illness. However, virulence descriptions of the E. coli O103 strains isolated from cattle feces have been primarily limited to major genes, such as Shiga toxin and intimin genes. Less is known about the frequency at which other virulence genes exist or about genes associated with the mobile genetic elements of E. coli O103 pathotypes. Our objective was to utilize whole genome sequencing (WGS) to identify and compare major and putative virulence genes of EHEC O103 (positive for Shiga toxin gene, stx1, and intimin gene, eae; n = 43), EPEC O103 (negative for stx1 and positive for eae; n = 13) and putative non-pathotype O103 strains (negative for stx and eae; n = 13) isolated from cattle feces. Six strains of EHEC O103 from human clinical cases were also included. All bovine EHEC strains (43/43) and a majority of EPEC (12/13) and putative non-pathotype strains (12/13) were O103:H2 serotype. Both bovine and human EHEC strains had significantly larger average genome sizes (P < 0.0001) and were positive for a higher number of adherence and toxin-based virulence genes and genes on mobile elements (prophages, transposable elements, and plasmids) than EPEC or putative non-pathotype strains. The genome size of the three pathotypes positively correlated (R² = 0.7) with the number of genes carried on mobile genetic elements. Bovine strains clustered phylogenetically by pathotypes, which differed in several key virulence genes. The diversity of E. coli O103 pathotypes shed in cattle feces is likely reflective of the acquisition or loss of virulence genes carried on mobile genetic elements.

Population structure of Escherichia coli O26 : H11 with recent and repeated stx2 acquisition in multiple lineages

A key virulence factor of enterohaemorrhagic Escherichia coli (EHEC) is the bacteriophage-encoded Shiga toxin (Stx). Stxs are classified into two types, Stx1 and Stx2, and Stx2-producing strains are thought to cause more severe infections than strains producing only Stx1. Although O26 : H11 is the second most prevalent EHEC following O157 : H7, the majority of O26 : H11 strains produce Stx1 alone. However, Stx2-producing O26 strains have increasingly been detected worldwide. Through a large-scale genome analysis, we present a global phylogenetic overview and evolutionary timescale for E. coli O26 : H11. The origin of O26 has been estimated to be 415 years ago. Sequence type 21C1 (ST21C1), one of the two sublineages of ST21, the most predominant O26 : H11 lineage worldwide, emerged 213 years ago from one of the three ST29 sublineages (ST29C2). The other ST21 lineage (ST21C2) emerged 95 years ago from ST21C1. Increases in population size occurred in the late 20th century for all of the O26 lineages, but most remarkably for ST21C2. Analysis of the distribution of stx2-positive strains revealed the recent and repeated acquisition of the stx2 gene in multiple lineages of O26, both in ST21 and ST29. Other major EHEC virulence genes, such as type III secretion system effector genes and plasmid-encoded virulence genes, were well conserved in ST21 compared to ST29. In addition, more antimicrobial-resistance genes have accumulated in the ST21C1 lineage. Although current attention is focused on several highly virulent ST29 clones that have acquired the stx2 gene, there is also a considerable risk that the ST21 lineage could yield highly virulent clones.

Real-time genomic investigation underlying the public health response to a Shiga toxin-producing Escherichia coli O26:H11 outbreak in a nursery

Shiga toxin-producing Escherichia coli (STEC) is a significant cause of gastrointestinal infection and the haemolytic-uremic syndrome (HUS). STEC outbreaks are commonly associated with food but animal contact is increasingly being implicated in its transmission. We report an outbreak of STEC affecting young infants at a nursery in a rural community (three HUS cases, one definite case, one probable case, three possible cases and five carriers, based on the combination of clinical, epidemiological and laboratory data) identified using culture-based and molecular techniques. The investigation identified repeated animal contact (animal farming and petting) as a likely source of STEC introduction followed by horizontal transmission. Whole genome sequencing (WGS) was used for real-time investigation of the incident and revealed a unique strain of STEC O26:H11 carrying stx2a and intimin. Following a public health intervention, no additional cases have occurred. This is the first STEC outbreak reported from Israel. WGS proved as a useful tool for rapid laboratory characterization and typing of the outbreak strain and informed the public health response at an early stage of this unusual outbreak.

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.

Advances in Molecular Serotyping and Subtyping of Escherichia coli

Escherichia coli plays an important role as a member of the gut microbiota; however, pathogenic strains also exist, including various diarrheagenic E. coli pathotypes and extraintestinal pathogenic E. coli that cause illness outside of the GI-tract. E. coli have traditionally been serotyped using antisera against the ca. 186 O-antigens and 53 H-flagellar antigens. Phenotypic methods, including bacteriophage typing and O- and H- serotyping for differentiating and characterizing E. coli have been used for many years; however, these methods are generally time consuming and not always accurate. Advances in next generation sequencing technologies have made it possible to develop genetic-based subtyping and molecular serotyping methods for E. coli, which are more discriminatory compared to phenotypic typing methods. Furthermore, whole genome sequencing (WGS) of E. coli is replacing established subtyping methods such as pulsed-field gel electrophoresis, providing a major advancement in the ability to investigate food-borne disease outbreaks and for trace-back to sources. A variety of sequence analysis tools and bioinformatic pipelines are being developed to analyze the vast amount of data generated by WGS and to obtain specific information such as O- and H-group determination and the presence of virulence genes and other genetic markers.

Targeted Amplicon Sequencing for Single-Nucleotide-Polymorphism Genotyping of Attaching and Effacing Escherichia coli O26:H11 Cattle Strains via a High-Throughput Library Preparation Technique

Enterohemorrhagic Escherichia coli (EHEC) O26:H11, a serotype within Shiga toxin-producing E. coli (STEC) that causes severe human disease, has been considered to have evolved from attaching and effacing E. coli (AEEC) O26:H11 through the acquisition of a Shiga toxin-encoding gene. Targeted amplicon sequencing using next-generation sequencing technology of 48 phylogenetically informative single-nucleotide polymorphisms (SNPs) and three SNPs differentiating Shiga toxin-positive (stx-positive) strains from Shiga toxin-negative (stx-negative) strains were used to infer the phylogenetic relationships of 178 E. coli O26:H11 strains (6 stx-positive strains and 172 stx-negative AEEC strains) from cattle feces to 7 publically available genomes of human clinical strains. The AEEC cattle strains displayed synonymous SNP genotypes with stx2-positive sequence type 29 (ST29) human O26:H11 strains, while stx1 ST21 human and cattle strains clustered separately, demonstrating the close phylogenetic relatedness of these Shiga toxin-negative AEEC cattle strains and human clinical strains. With the exception of seven stx-negative strains, five of which contained espK, three stx-related SNPs differentiated the STEC strains from non-STEC strains, supporting the hypothesis that these AEEC cattle strains could serve as a potential reservoir for new or existing pathogenic human strains. Our results support the idea that targeted amplicon sequencing for SNP genotyping expedites strain identification and genetic characterization of E. coli O26:H11, which is important for food safety and public health.