Whole Genome Sequencing demonstrates that Geographic Variation of Escherichia coli O157 Genotypes Dominates Host Association

Commensal Escherichia coli Strains of Bovine Origin Competitively Mitigated Escherichia coli O157:H7 in a Gnotobiotic Murine Intestinal Colonization Model with or without Physiological Stress

Cattle are a primary reservoir of enterohemorrhagic Escherichia coli (EHEC) O157:H7. Currently, there are no effective methods of eliminating this important zoonotic pathogen from cattle, and colonization resistance in relation to EHEC O157:H7 in cattle is poorly understood. We developed a gnotobiotic EHEC O157:H7 murine model to examine aspects of the cattle pathogen-microbiota interaction, and to investigate competitive suppression of EHEC O157:H7 by 18 phylogenetically distinct commensal E. coli strains of bovine origin. As stress has been suggested to influence enteric colonization by EHEC O157:H7 in cattle, corticosterone administration (±) to incite a physiological stress response was included as an experimental variable. Colonization of the intestinal tract (IT) of mice by the bovine EHEC O157:H7 strain, FRIK-2001, mimicked characteristics of bovine IT colonization. In this regard, FRIK-2001 successfully colonized the IT and temporally incited minimal impacts on the host relative to other EHEC O157:H7 strains, including on the renal metabolome. The presence of the commensal E. coli strains decreased EHEC O157:H7 densities in the cecum, proximal colon, and distal colon. Moreover, histopathologic changes and inflammation markers were reduced in the distal colon of mice inoculated with commensal E. coli strains (both propagated separately and communally). Although stress induction affected the behavior of mice, it did not influence EHEC O157:H7 densities or disease. These findings support the use of a gnotobiotic murine model of enteric bovine EHEC O157:H7 colonization to better understand pathogen-host-microbiota interactions toward the development of effective on-farm mitigations for EHEC O157:H7 in cattle, including the identification of bacteria capable of competitively colonizing the IT.

Host Association and Spatial Proximity Shape but Do Not Constrain Population Structure in the Mutualistic Symbiont Xenorhabdus bovienii

To what extent are generalist species cohesive evolutionary units rather than a compilation of recently diverged lineages? We examine this question in the context of host specificity and geographic structure in the insect pathogen and nematode mutualist Xenorhabdus bovienii. This bacterial species partners with multiple nematode species across two clades in the genus Steinernema. We sequenced the genomes of 42 X. bovienii strains isolated from four different nematode species and three field sites within a 240-km² region and compared them to globally available reference genomes. We hypothesized that X. bovienii would comprise several host-specific lineages, such that bacterial and nematode phylogenies would be largely congruent. Alternatively, we hypothesized that spatial proximity might be a dominant signal, as increasing geographic distance might lower shared selective pressures and opportunities for gene flow. We found partial support for both hypotheses. Isolates clustered largely by nematode host species but did not strictly match the nematode phylogeny, indicating that shifts in symbiont associations across nematode species and clades have occurred. Furthermore, both genetic similarity and gene flow decreased with geographic distance across nematode species, suggesting differentiation and constraints on gene flow across both factors, although no absolute barriers to gene flow were observed across the regional isolates. Several genes associated with biotic interactions were found to be undergoing selective sweeps within this regional population. The interactions included several insect toxins and genes implicated in microbial competition. Thus, gene flow maintains cohesiveness across host associations in this symbiont and may facilitate adaptive responses to a multipartite selective environment. IMPORTANCE Microbial populations and species are notoriously hard to delineate. We used a population genomics approach to examine the population structure and the spatial scale of gene flow in Xenorhabdus bovienii, an intriguing species that is both a specialized mutualistic symbiont of nematodes and a broadly virulent insect pathogen. We found a strong signature of nematode host association, as well as evidence for gene flow connecting isolates associated with different nematode host species and collected from distinct study sites. Furthermore, we saw signatures of selective sweeps for genes involved with nematode host associations, insect pathogenicity, and microbial competition. Thus, X. bovienii exemplifies the growing consensus that recombination not only maintains cohesion but can also allow the spread of niche-beneficial alleles.

Assessment of Prevalence and Diversity of Antimicrobial Resistant Escherichia coli from Retail Meats in Southern California

Retail meat products may serve as reservoirs and conduits for antimicrobial resistance, which is frequently monitored using Escherichia coli as indicator bacteria. In this study, E. coli isolation was conducted on 221 retail meat samples (56 chicken, 54 ground turkey, 55 ground beef, and 56 pork chops) collected over a one-year period from grocery stores in southern California. The overall prevalence of E. coli in retail meat samples was 47.51% (105/221), with E. coli contamination found to be significantly associated with meat type and season of sampling. From antimicrobial susceptibility testing, 51 isolates (48.57%) were susceptible to all antimicrobials tested, 54 (51.34%) were resistant to at least 1 drug, 39 (37.14%) to 2 or more drugs, and 21 (20.00%) to 3 or more drugs. Resistance to ampicillin, gentamicin, streptomycin, and tetracycline were significantly associated with meat type, with poultry counterparts (chicken or ground turkey) exhibiting higher odds for resistance to these drugs compared to non-poultry meats (beef and pork). From the 52 E. coli isolates selected to undergo whole-genome sequencing (WGS), 27 antimicrobial resistance genes (ARGs) were identified and predicted phenotypic AMR profiles with an overall sensitivity and specificity of 93.33% and 99.84%, respectively. Clustering assessment and co-occurrence networks revealed that the genomic AMR determinants of E. coli from retail meat were highly heterogeneous, with a sparsity of shared gene networks.

Prevalence and Implications of Shiga Toxin-Producing E. coli in Farm and Wild Ruminants

Shiga-toxin-producing Escherichia coli (STEC) is a food-borne pathogen that causes human gastrointestinal infections across the globe, leading to kidney failure or even death in severe cases. E. coli are commensal members of humans and animals' (cattle, bison, and pigs) guts, however, may acquire Shiga-toxin-encoded phages. This acquisition or colonization by STEC may lead to dysbiosis in the intestinal microbial community of the host. Wildlife and livestock animals can be asymptomatically colonized by STEC, leading to pathogen shedding and transmission. Furthermore, there has been a steady uptick in new STEC variants representing various serotypes. These, along with hybrids of other pathogenic E. coli (UPEC and ExPEC), are of serious concern, especially when they possess enhanced antimicrobial resistance, biofilm formation, etc. Recent studies have reported these in the livestock and food industry with minimal focus on wildlife. Disturbed natural habitats and changing climates are increasingly creating wildlife reservoirs of these pathogens, leading to a rise in zoonotic infections. Therefore, this review comprehensively surveyed studies on STEC prevalence in livestock and wildlife hosts. We further present important microbial and environmental factors contributing to STEC spread as well as infections. Finally, we delve into potential strategies for limiting STEC shedding and transmission.

Clarification of relationship between single-nucleotide polymorphism panels of Shiga toxin-producing Escherichia coli O157:H7/H- strains

Eighty strains of enterohemorrhagic Escherichia coli O157:H7/H- were analyzed by three single-nucleotide polymorphism (SNP) panels using whole-genome sequencing data. The partial concordance of SNP types among the different SNP panels was observed on minimum spanning trees reconstructed with SNP data. As for lineage I/II strains, some of the clade 7 strains belonged to one unique SNP type as determined by three panels, suggesting that clade 7 should be divided into at least two genotypes, namely, the unique type and the rest. In addition, clade 8 contained two unique genotypes, which was consistent with the previous prediction. Similarly, for lineage II, clade 12 should be divided into three genotype strains. In contrast, many strains of several clades belonging to lineage I were clustered into the same node on each minimum spanning tree upon testing with the three SNP panels. Previous studies reported that lineage I diverged more recently than lineages I/II and II. Such low diversity in lineage I in this study may have arisen because this lineage has not accumulated SNPs because of its relatively recent divergence. Based on the concordance observed in this study, some of the previously published O157 genotype distribution data were successfully interpreted to clarify the clade distribution, which was well supported by previous literature.

Pathogenomes and variations in Shiga toxin production among geographically distinct clones of Escherichia coli O113:H21

Infections with globally disseminated Shiga toxin-producing Escherichia coli (STEC) of the O113:H21 serotype can progress to severe clinical complications, such as hemolytic uremic syndrome (HUS). Two phylogeographically distinct clonal complexes have been established by multi locus sequence typing (MLST). Infections with ST-820 isolates circulating exclusively in Australia have caused severe human disease, such as HUS. Conversely, ST-223 isolates prevalent in the US and outside Australia seem to rarely cause severe human disease but are frequent contaminants. Following a genomic epidemiology approach, we wanted to gain insights into the underlying cause for this disparity. We examined the plasticity in the genome make-up and Shiga toxin production in a collection of 20 ST-820 and ST-223 strains isolated from produce, the bovine reservoir, and clinical cases. STEC are notorious for assembly into fragmented draft sequences when using short-read sequencing technologies due to the extensive and partly homologous phage complement. The application of long-read technology (LRT) sequencing yielded closed reference chromosomes and plasmids for two representative ST-820 and ST-223 strains. The established high-resolution framework, based on whole genome alignments, single nucleotide polymorphism (SNP)-typing and MLST, includes the chromosomes and plasmids of other publicly available O113:H21 sequences and allowed us to refine the phylogeographical boundaries of ST-820 and ST-223 complex isolates and to further identify a historic non-shigatoxigenic strain from Mexico as a quasi-intermediate. Plasmid comparison revealed strong correlations between the strains' featured pO113 plasmid genotypes and chromosomally inferred ST, which suggests coevolution of the chromosome and virulence plasmids. Our pathogenicity assessment revealed statistically significant differences in the Stx_2a-production capabilities of ST-820 as compared to ST-223 strains under RecA-induced Stx phage mobilization, a condition that mimics Stx-phage induction. These observations suggest that ST-820 strains may confer an increased pathogenic potential in line with the strain-associated epidemiological metadata. Still, some of the tested ST-223 cultures sourced from contaminated produce or the bovine reservoir also produced Stx at levels comparable to those of ST-820 isolates, which calls for awareness and for continued surveillance of this lineage.

Development of Single Nucleotide Polymorphism (SNP)-Based Triplex PCR Marker for Serotype-Specific Escherichia coli Detection

Single-nucleotide polymorphisms (SNPs) are one of the most common forms of genetic variation and as such are powerful tools for the identification of bacterial strains, their genetic diversity, phylogenetic analysis, and outbreak surveillance. In this study, we used 15 sets of SNP-containing primers to amplify and sequence the target Escherichia coli. Based on the combination of the 15-sequence primer sets, each SNP site encompassing forward and reverse primer sequences (620–919 bp) were aligned and an SNP-based marker was designed. Each SNP marker exists in at least two SNP sites at the 3′ end of each primer; one natural and the other artificially created by transition or transversion mutation. Thus, 12 sets of SNP primers (225–488 bp) were developed for validation by amplifying the target E. coli. Finally, a temperature gradient triplex PCR kit was designed to detect target E. coli strains. The selected primers were amplified in three genes (ileS, thrB, and polB), with fragment sizes of 401, 337, and 232 bp for E. coli O157:H7, E. coli, and E. coli O145:H28, respectively. This allele-specific SNP-based triplex primer assay provides serotype-specific detection of E. coli strains in one reaction tube. The developed marker would be used to diagnose, investigate, and control food-borne E. coli outbreaks.

Genome Sequence Analysis and Characterization of Shiga Toxin 2 Production by Escherichia coli O157:H7 Strains Associated With a Laboratory Infection

A laboratory-acquired E. coli O157:H7 infection with associated severe sequelae including hemolytic uremic syndrome occurred in an individual working in the laboratory with a mixture of nalidixic acid-resistant (Nal^R) O157:H7 mutant strains in a soil-biochar blend. The patient was hospitalized and treated with an intravenous combination of metronidazole and levofloxacin. The present study investigated the source of this severe laboratory acquired infection and further examined the influence of the antibiotics used during treatment on the expression and production of Shiga toxin. Genomes of two Stx_2a-and eae-positive O157:H7 strains isolated from the patient's stool were sequenced along with two pairs of the wt strains and their derived Nal^R mutants used in the laboratory experiments. High-resolution SNP typing determined the strains' individual genetic relatedness and unambiguously identified the two laboratory-derived Nal^R mutant strains as the source of the researcher's life-threatening disease, rather than a conceivable ingestion of unrelated O157:H7 isolates circulating at the same time. It was further confirmed that in sublethal doses, the antibiotics increased toxin expression and production. Our results support a simultaneous co-infection with clinical strains in the laboratory, which were the causative agents of previous O157:H7 outbreaks, and further that the administration of antibiotics may have impacted the outcome of the infection.

Genomic Characterization of Endemic and Ecdemic Non-typhoidal Salmonella enterica Lineages Circulating Among Animals and Animal Products in South Africa

In Africa, the burden of illness caused by non-typhoidal Salmonella enterica is disproportionally high; however, whole-genome sequencing (WGS) efforts are overwhelmingly concentrated in world regions with lower burdens. While WGS is being increasingly employed in South Africa to characterize Salmonella enterica, the bulk of these efforts have centered on characterizing human clinical strains. Thus, very little is known about lineages circulating among animals in the country on a genomic scale. Here, we used WGS to characterize 63 Salmonella enterica strains isolated from livestock, companion animals, wildlife, and animal products in South Africa over a 60-year period. Genomes were assigned to serotypes Dublin, Hadar, Enteritidis, and Typhimurium (n = 18, 8, 13, and 24 strains, respectively) and sequence types (STs) ST10 (all S. Dublin), ST33 (all S. Hadar), ST11/ST366 (n = 12 and 1 S. Enteritidis, respectively), and ST19/ST34 (n = 23 and 1 S. Typhimurium, respectively; via seven-gene multi-locus sequence typing). Within-ST phylogenies were constructed using genomes sequenced in this study, plus publicly available genomes representative of each ST's (i) global (n = 2,802 and 1,569 S. Dublin and Hadar genomes, respectively) and (ii) African (n = 716 and 343 S. Enteritidis and Typhimurium genomes, respectively) population. For S. Dublin ST10, a largely antimicrobial-susceptible, endemic lineage circulating among humans, animals, and food in South Africa was identified, as well as a lineage that was likely recently introduced from the United States. For S. Hadar ST33, multiple South African lineages harboring streptomycin and tetracycline resistance-conferring genes were identified. African S. Enteritidis ST11 could be primarily partitioned into one largely antimicrobial-susceptible and one largely multidrug-resistant (MDR) clade, with South African isolates confined to the largely antimicrobial-susceptible clade. S. Typhimurium ST19/ST34 strains sequenced here were distributed across the African S. Typhimurium ST19/ST34 phylogeny, representing a diverse range of lineages, including numerous MDR lineages. Overall, this study provides critical insights into endemic and ecdemic non-typhoidal Salmonella enterica lineages circulating among animals, foods, and humans in South Africa and showcases the utility of WGS in characterizing animal-associated strains from a world region with a high salmonellosis burden.

Monitoring the Microevolution of Salmonella enterica in Healthy Dairy Cattle Populations at the Individual Farm Level Using Whole-Genome Sequencing

Livestock represent a possible reservoir for facilitating the transmission of the zoonotic foodborne pathogen Salmonella enterica to humans; there is also concern that strains can acquire resistance to antimicrobials in the farm environment. Here, whole-genome sequencing (WGS) was used to characterize Salmonella strains (n = 128) isolated from healthy dairy cattle and their associated environments on 13 New York State farms to assess the diversity and microevolution of this important pathogen at the level of the individual herd. Additionally, the accuracy and concordance of multiple in silico tools are assessed, including: (i) two in silico serotyping tools, (ii) combinations of five antimicrobial resistance (AMR) determinant detection tools and one to five AMR determinant databases, and (iii) one antimicrobial minimum inhibitory concentration (MIC) prediction tool. For the isolates sequenced here, in silico serotyping methods outperformed traditional serotyping and resolved all un-typable and/or ambiguous serotype assignments. Serotypes assigned in silico showed greater congruency with the Salmonella whole-genome phylogeny than traditional serotype assignments, and in silico methods showed high concordance (99% agreement). In silico AMR determinant detection methods additionally showed a high degree of concordance, regardless of the pipeline or database used (≥98% agreement among susceptible/resistant assignments for all pipeline/database combinations). For AMR detection methods that relied exclusively on nucleotide BLAST, accuracy could be maximized by using a range of minimum nucleotide identity and coverage thresholds, with thresholds of 75% nucleotide identity and 50-60% coverage adequate for most pipeline/database combinations. In silico characterization of the microevolution and AMR dynamics of each of six serotype groups (S. Anatum, Cerro, Kentucky, Meleagridis, Newport, Typhimurium/Typhimurium variant Copenhagen) revealed that some lineages were strongly associated with individual farms, while others were distributed across multiple farms. Numerous AMR determinant acquisition and loss events were identified, including the recent acquisition of cephalosporin resistance-conferring bla CMY- and bla CTX-M-type beta-lactamases. The results presented here provide high-resolution insight into the temporal dynamics of AMR Salmonella at the scale of the individual farm and highlight both the strengths and limitations of WGS in tracking zoonotic pathogens and their associated AMR determinants at the livestock-human interface.

Genomic analysis of Shiga toxin-producing Escherichia coli O157:H7 from cattle and pork-production related environments

Three E. coli O157:H7 outbreaks have been attributed to contaminated pork in Alberta, Canada, recently. This study investigates the phylogenetic relatedness of E. coli O157:H7 from pigs, cattle, and pork-production environments for source attribution. Limited strain diversity was observed using five conventional subtyping methods, with most or all strains being in one subgroup. Whole-genome single nucleotide polymorphism analysis confirmed the recent ancestry of the isolates from all three sources. Most environmental isolates clustered closer with pig isolates than cattle isolates. Also, a direct link was observed between 2018-outbreak environmental isolates and isolates collected from a pig farm in 2018. The majority of pig isolates harbor only one Shiga toxin gene, stx_2a, while 70% (35/50) of the cattle isolates have both stx_1a and stx_2a. The results show some E. coli O157:H7 strains could establish persistence on pig farms and as such, pigs can be a significant source of the organism.

Genomic Epidemiology of Shiga Toxin-Producing Escherichia coli Isolated from the Livestock-Food-Human Interface in South America

Simple Summary

Shiga toxin-producing Escherichia coli (STEC) are zoonotic pathogens that cause food-borne diseases in humans, where cattle and derived products play a key role as reservoirs and vehicles. We analyzed the genomic data of STEC strains circulating at the livestock-food-human interface in South America, extracting clinically and epidemiologically relevant information (serotypes, virulome, resistance genes, sequence types, and phylogenomics). This study included 130 STEC genomes obtained from cattle (n = 51), beef (n = 48), and human (n = 31) samples. The successful expansion of O157:H7 (ST11) and non-O157 (ST16, ST21, ST223, ST443, ST677, ST679, ST2388) clones is highlighted, suggesting common activities, such as multilateral trade and travel. Circulating STEC strains analyzed exhibit high genomic diversity and harbor several genetic determinants associated with severe illness in humans, highlighting the need to establish official surveillance of this pathogen that should be focused on detecting molecular determinants of virulence and clonal relatedness, in the whole beef production chain.

Abstract

Shiga toxin-producing Escherichia coli (STEC) are zoonotic pathogens responsible for causing food-borne diseases in humans. While South America has the highest incidence of human STEC infections, information about the genomic characteristics of the circulating strains is scarce. The aim of this study was to analyze genomic data of STEC strains isolated in South America from cattle, beef, and humans; predicting the antibiotic resistome, serotypes, sequence types (STs), clonal complexes (CCs) and phylogenomic backgrounds. A total of 130 whole genome sequences of STEC strains were analyzed, where 39.2% were isolated from cattle, 36.9% from beef, and 23.8% from humans. The ST11 was the most predicted (20.8%) and included O-:H7 (10.8%) and O157:H7 (10%) serotypes. The successful expansion of non-O157 clones such as ST16/CC29-O111:H8 and ST21/CC29-O26:H11 is highlighted, suggesting multilateral trade and travel. Virulome analyses showed that the predominant stx subtype was stx2a (54.6%); most strains carried ehaA (96.2%), iha (91.5%) and lpfA (77.7%) genes. We present genomic data that can be used to support the surveillance of STEC strains circulating at the livestock-food-human interface in South America, in order to control the spread of critical clones “from farm to table”.

Phylogeny, Prevalence, and Shiga Toxin (Stx) Production of Clinical Escherichia coli O157 Clade 2 Strains Isolated in Shimane Prefecture, Japan

This study investigated the genetic and pathogenic variation of the subgroups of clade 2 strains of Shiga toxin (Stx)-producing Escherichia coli (STEC) O157. A total of 111 strains of STEC O157 isolated in Shimane prefecture, Japan, were classified in clade 2 (n = 39), clade 3 (n = 16), clade 4/5 (n = 3), clade 7 (n = 14), clade 8 (n = 17), and clade 12 (n = 22) by single-nucleotide polymorphism analysis and lineage-specific polymorphism assay-6. These results showed a distinct difference from our previous study in which clade 3 strains were the most prevalent strains in three other prefectures in Japan, indicating that the clade distribution of O157 strains was different in different geographic areas in Japan. Phylogenetic analysis using insertion sequence (IS) 629 distribution data showed that clade 2 strains formed two clusters, designated 2a and 2b. Stx2 production by cluster 2b strains was significantly higher than by cluster 2a strains (P < 0.01). In addition, population genetic analysis of the clade 2 strains showed significant linkage disequilibrium in the IS629 distribution of the strains in clusters 2a and 2b (P < 0.05). The Φ_PT values calculated using the IS629 distribution data indicated that strains in clusters 2a and 2b were genetically different (P < 0.001). Cluster 2b strains are a highly pathogenic phylogenetic group and their geographic spread may be a serious public health concern.

The population genetics of pathogenic Escherichia coli

Escherichia coli is a commensal of the vertebrate gut that is increasingly involved in various intestinal and extra-intestinal infections as an opportunistic pathogen. Numerous pathotypes that represent groups of strains with specific pathogenic characteristics have been described based on heterogeneous and complex criteria. The democratization of whole-genome sequencing has led to an accumulation of genomic data that render possible a population phylogenomic approach to the emergence of virulence. Few lineages are responsible for the pathologies compared with the diversity of commensal strains. These lineages emerged multiple times during E. coli evolution, mainly by acquiring virulence genes located on mobile elements, but in a specific chromosomal phylogenetic background. This repeated emergence of stable and cosmopolitan lineages argues for an optimization of strain fitness through epistatic interactions between the virulence determinants and the remaining genome.

Phylogeographic Analysis Reveals Multiple International transmission Events Have Driven the Global Emergence of Escherichia coli O157:H7

BACKGROUND

Shiga toxin-producing Escherchia coli (STEC) O157:H7 is a zoonotic pathogen that causes numerous food and waterborne disease outbreaks. It is globally distributed, but its origin and the temporal sequence of its geographical spread are unknown.

METHODS

We analyzed whole-genome sequencing data of 757 isolates from 4 continents, and performed a pan-genome analysis to identify the core genome and, from this, extracted single-nucleotide polymorphisms. A timed phylogeographic analysis was performed on a subset of the isolates to investigate its worldwide spread.

RESULTS

The common ancestor of this set of isolates occurred around 1890 (1845-1925) and originated from the Netherlands. Phylogeographic analysis identified 34 major transmission events. The earliest were predominantly intercontinental, moving from Europe to Australia around 1937 (1909-1958), to the United States in 1941 (1921-1962), to Canada in 1960 (1943-1979), and from Australia to New Zealand in 1966 (1943-1982). This pre-dates the first reported human case of E. coli O157:H7, which was in 1975 from the United States.

CONCLUSIONS

Inter- and intra-continental transmission events have resulted in the current international distribution of E. coli O157:H7, and it is likely that these events were facilitated by animal movements (eg, Holstein Friesian cattle). These findings will inform policy on action that is crucial to reduce the further spread of E. coli O157:H7 and other (emerging) STEC strains globally.

Impact of whole genome sequencing on the investigation of food-borne outbreaks of Shiga toxin-producing Escherichia coli serogroup O157:H7, England, 2013 to 2017

We aim to provide insight and guidance on the utility of whole genome sequencing (WGS) data for investigating food-borne outbreaks of Shiga toxin-producing Escherichia coli (STEC) O157:H7 in England between 2013 and 2017. Analysis of WGS data delivered an unprecedented level of strain discrimination when compared with multilocus variable number tandem repeat analysis. The robustness of the WGS method ensured confidence in the microbiological identification of linked cases, even when epidemiological links were obscured. There was evidence that phylogeny derived from WGS data can be used to trace the geographical origin of an isolate. Further analysis of the phylogenetic data provided insight on the evolutionary context of emerging pathogenic strains. Publically available WGS data linked to the clinical, epidemiological and environmental context of the sequenced strain has improved trace back investigations during outbreaks. Expanding the use of WGS-based typing analysis globally will ensure the rapid implementation of interventions to protect public health, inform risk assessment and facilitate the management of national and international food-borne outbreaks of STEC O157:H7.

Mining whole genome sequence data to efficiently attribute individuals to source populations

Whole genome sequence (WGS) data could transform our ability to attribute individuals to source populations. However, methods that efficiently mine these data are yet to be developed. We present a minimal multilocus distance (MMD) method which rapidly deals with these large data sets as well as methods for optimally selecting loci. This was applied on WGS data to determine the source of human campylobacteriosis, the geographical origin of diverse biological species including humans and proteomic data to classify breast cancer tumours. The MMD method provides a highly accurate attribution which is computationally efficient for extended genotypes. These methods are generic, easy to implement for WGS and proteomic data and have wide application.

Comparative genomics reveals different population structures associated with host and geographic origin in antimicrobial-resistant Salmonella enterica

Genetic variation in a pathogen, including the causative agent of salmonellosis, Salmonella enterica, can occur as a result of eco-evolutionary forces triggered by dissimilarities of ecological niches. Here, we applied comparative genomics to study 90 antimicrobial resistant (AMR) S. enterica isolates from bovine and human hosts in New York and Washington states to understand host- and geographic-associated population structure. Results revealed distinct presence/absence profiles of functional genes and pseudogenes (e.g., virulence genes) associated with bovine and human isolates. Notably, bovine isolates contained significantly more transposase genes but fewer transposase pseudogenes than human isolates, suggesting the occurrence of large-scale transposition in genomes of bovine and human isolates at different times. The high correlation between transposase genes and AMR genes, as well as plasmid replicons, highlights the potential role of horizontally transferred transposons in promoting adaptation to antibiotics. By contrast, a number of potentially geographic-associated single-nucleotide polymorphisms (SNPs), rather than geographic-associated genes, were identified. Interestingly, 38% of these SNPs were in genes annotated as cell surface protein-encoding genes, including some essential for antibiotic resistance and host colonization. Overall, different evolutionary forces and limited recent inter-population transmission appear to shape AMR S. enterica population structure in different hosts and geographic origins.

High Genomic Diversity and Heterogenous Origins of Pathogenic and Antibiotic-Resistant Escherichia coli in Household Settings Represent a Challenge to Reducing Transmission in Low-Income Settings

Escherichia coli is present in multiple hosts and environmental compartments as a normal inhabitant, temporary or persistent colonizer, and as a pathogen. Transmission of E. coli between hosts and with the environment is considered to occur more often in areas with poor sanitation. We performed whole-genome comparative analyses on 60 E. coli isolates from soils and fecal sources (cattle, chickens, and humans) in households in rural Bangladesh. Isolates from household soils were in multiple branches of the reconstructed phylogeny, intermixed with isolates from fecal sources. Pairwise differences between all strain pairs were large (minimum, 189 single nucleotide polymorphisms [SNPs]), suggesting high diversity and heterogeneous origins of the isolates. The presence of multiple virulence and antibiotic resistance genes is indicative of the risk that E. coli from soil and feces represent for the transmission of variants that pose potential harm to people. Analysis of the accessory genomes of the Bangladeshi E. coli relative to E. coli genomes available in NCBI identified a common pool of accessory genes shared among E. coli isolates in this geographic area. Together, these findings indicate that in rural Bangladesh, a high level of E. coli in soil is likely driven by contributions from multiple and diverse E. coli sources (human and animal) that share an accessory gene pool relatively unique to previously published E. coli genomes. Thus, interventions to reduce environmental pathogen or antimicrobial resistance transmission should adopt integrated One Health approaches that consider heterogeneous origins and high diversity to improve effectiveness and reduce prevalence and transmission.IMPORTANCEEscherichia coli is reported in high levels in household soil in low-income settings. When E. coli reaches a soil environment, different mechanisms, including survival, clonal expansion, and genetic exchange, have the potential to either maintain or generate E. coli variants with capabilities of causing harm to people. In this study, we used whole-genome sequencing to identify that E. coli isolates collected from rural Bangladeshi household soils, including pathogenic and antibiotic-resistant variants, are diverse and likely originated from multiple diverse sources. In addition, we observed specialization of the accessory genome of this Bangladeshi E. coli compared to E. coli genomes available in current sequence databases. Thus, to address the high level of pathogenic and antibiotic-resistant E. coli transmission in low-income settings, interventions should focus on addressing the heterogeneous origins and high diversity.

Comparable stx2a expression and phage production levels between Shiga toxin-producing Escherichia coli strains from human and bovine origin

Shiga toxin-producing Escherichia coli (STEC) can cause diarrhoea and severe diseases in humans, such as haemolytic uraemic syndrome. STEC virulence is considered to correlate with the amount of Shiga toxins (Stx) produced, especially Stx2, whose subtype Stx2a is most frequently associated with high virulence. Stx are encoded in prophages, which play an important role in STEC pathogenesis. The aim of this study was to evaluate stx_2a expression levels and Stx2a phage production using qPCR and the double-agar-layer method in 29 STEC strains, corresponding to serotypes O26:H11 (6), O91:H21 (1), O145:H- (11) and O157:H7 (11), isolated from cattle and humans. Results were then tested for possible associations with serotype, origin or some genetic features. We observed heterogeneous levels of stx_2a expression and Stx2a phage production. However, statistical comparisons identified a higher stx_2a expression in response to mitomycin C in strains isolated from cattle than in those from humans. At the same time, compared to stx_2a /stx_2c strains, stx_2a strains showed a higher increase in phage production under induced conditions. Notably, most of the strains studied, regardless of serotype and origin, carried inducible Stx2a phages and evidenced expression of stx_2a that increased along with phage production levels under induced conditions.

Application of kernel smoothing to estimate the spatio-temporal variation in risk of STEC O157 in England

Identifying geographical areas with significantly higher or lower rates of infectious diseases can provide important aetiological clues to inform the development of public health policy and interventions designed to reduce morbidity. We applied kernel smoothing to estimate the spatial and spatio-temporal variation in risk of STEC O157 infection in England between 2009 and 2015, and to explore differences between the residential locations of cases reporting travel and those not reporting travel. We provide evidence that the distribution of STEC O157 infection in England is non-uniform with respect to the distribution of the at-risk population; that the spatial distribution of the three main genetic lineages infecting humans (I, II and I/II) differs significantly and that the spatio-temporal risk is highly dynamic. Our results also indicate that cases of STEC O157 reporting travel within or outside the UK are more likely to live in the south/south-east of the country, meaning that their residential location may not reflect the location of exposure that led to their infection. We suggest that the observed variation in risk reflects exposure to sources of STEC O157 that are geographically prescribed. These differences may be related to a combination of changes in the strains circulating in the ruminant reservoir, animal movements (livestock, birds or wildlife) or the behavior of individuals prior to infection. Further work to identify the importance of behaviours and exposures reported by cases relative to residential location is needed.

Risk factors and dynamics of verotoxigenic Escherichia coli O157:H7 on cattle farms: An observational study combining information from questionnaires, spatial data and molecular analyses

The increasing number of human cases infected with a highly virulent type of verotoxigenic Escherichia coli (VTEC) O157:H7 in Sweden is the result of domestic transmission originating in regional clusters of infected cattle farms. To control the spread of the bacteria a comprehensive picture of infection dynamics, routes of transmission between farms and risk factors for persistence is urgently needed. The aim of the study was to investigate different aspects of the epidemiology of VTEC O157:H7 on the Swedish island of Öland by combining information from environmental sampling of VTEC O157:H7 from 80 farms with information from farmer questionnaires, spatial and molecular analyses. The farms were sampled in the spring and fall of 2014 and on four of them additional samples were collected during summer and winter. The results show a high prevalence of VTEC O157:H7 and a high proportion of strains belonging to the virulent clade 8. Farms that became infected between samplings were all located in an area with high cattle density. The most important risk factors identified are generally associated with biosecurity and indicate that visitors travelling between farms may be important for transmission. In addition, whole genome sequencing of a subset of isolates from the four farms where additional sampling was performed revealed ongoing local transmission that cannot be observed with a lower resolution typing method. Our observations also show that VTEC O157:H7 may persist in the farm environment for extended periods of time, suggesting that specific on-farm measures to reduce environmental prevalence and spread between groups of animals may be required in these cases.

Genome-wide identification of geographical segregated genetic markers in Salmonella enterica serovar Typhimurium variant 4,[5],12:i:

Salmonella enterica ser. Typhimurium monophasic variant 4,[5],12:i:- has been associated with food-borne epidemics worldwide and swine appeared to be the main reservoir in most of the countries of isolation. However, the monomorphic nature of this serovar has, so far, hindered identification of the source due to expansion of clonal lineages in multiple hosts and food producing systems. Since geographically structured genetic signals can shape bacterial populations, identification of biogeographical markers in S. 1,4,[5],12:i:- genomes can contribute to improving source attribution. In this study, the phylogeographical structure of 148 geographically and temporally related Italian S. 1,4,[5],12:i:- has been investigated. The Italian isolates belong to a large population of clonal S. Typhimurium/1,4,[5],12:i:- isolates collected worldwide in two decades showing up to 2.5% of allele differences. Phylogenetic reconstruction revealed that isolates from the same geographical origin form highly supported monophyletic groups, suggesting discrete geographical segregation. These monophyletic groups are characterized by the gene content of a large sopE-containing prophage. Within this prophage, genome-wide comparison identified several genes overrepresented in strains of Italian origin. This suggests that certain lineages may be characterized by the acquisition of specific accessory genetic markers useful for improving identification of the source in ongoing epidemics.

Comparative Genomic Studies of Salmonella Heidelberg Isolated From Chicken- and Turkey-Associated Farm Environmental Samples

Salmonella is one of the leading causes of human foodborne gastroenteritis in the United States. In addition, Salmonella contributes to morbidity and mortality in livestock. The control of Salmonella is an increasing problematic issue in livestock production due to lack of effective control methods and the constant adaptation of Salmonella to new management practices, which is often related to horizontal acquisition of virulence or antibiotic resistance genes. Salmonella enterica serotype Heidelberg is one of the most commonly isolated serotypes in all poultry production systems in North America. Emergence and persistence of multi-drug resistant Salmonella Heidelberg isolates further impact the poultry production and public health. We hypothesized that distinct poultry production environments affect Salmonella genomic content, and by consequence its survival and virulence abilities. This study compared the genomic composition of S. Heidelberg isolated from environmental samples (19 chicken and 12 turkey isolates) of different breeder farms (16 chicken and 8 turkey farms) in the Midwest, United States. Whole genome comparison of 31 genomes using RAST and SEED identified differences in specific sub-systems in isolates between the chicken- and turkey-associated farm environmental samples. Genes associated with the type IV secretion system (n = 12) and conjugative transfer (n = 3) were absent in turkey farm isolates compared to the chicken ones (p-value < 0.01); Further, turkey farm isolates were enriched in prophage proteins (n = 53; p-value < 0.01). Complementary studies using PHASTER showed that prophages were all Caudovirales phages and were more represented in turkey environmental isolates than the chicken isolates. This study corroborates that isolates from distinct farm environment show differences in S. Heidelberg genome content related to horizontal transfer between bacteria or through viral infections. Complementary microbiome studies of these samples would provide critical insights on sources of these variations. Overall, our findings enhance the understanding of Salmonella genome plasticity and may aid in the development of future effective management practices to control Salmonella.

Geogenomic Segregation and Temporal Trends of Human Pathogenic Escherichia coli O157:H7, Washington, USA, 2005-20141

The often-noted and persistent increased incidence of Escherichia coli O157:H7 infections in rural areas is not well understood. We used a cohort of E. coli O157:H7 cases reported in Washington, USA, during 2005–2014, along with phylogenomic characterization of the infecting isolates, to identify geographic segregation of and temporal trends in specific phylogenetic lineages of E. coli O157:H7. Kernel estimation and generalized additive models demonstrated that pathogen lineages were spatially segregated during the period of analysis and identified a focus of segregation spanning multiple, predominantly rural, counties for each of the main clinical lineages, Ib, IIa, and IIb. These results suggest the existence of local reservoirs from which humans are infected. We also noted a secular increase in the proportion of lineage IIa and IIb isolates. Spatial segregation by phylogenetic lineage offers the potential to identify local reservoirs and intervene to prevent continued transmission.

Variability in Characterizing Escherichia coli from Cattle Feces: A Cautionary Tale

Shiga toxin-producing Escherichia coli (STEC) are diverse bacteria, with seven serogroups (O26, O45, O103, O111, O121, O145, O157; "Top 7") of interest due to their predominance in human disease. Confirmation of STEC relies on a combination of culturing, immunological and molecular assays, but no single gold standard for identification exists. In this study, we compared analysis of STEC between three independent laboratories (LAB) using different methodologies. In LAB A, colonies of Top 7 were picked after serogroup-specific immunomagnetic separation of feces from western-Canadian slaughter cattle. A fraction of each colony was tested by PCR (stx1, stx2, eae, O group), and Top 7 isolates were saved as glycerol stocks (n = 689). In LAB B, a subsample of isolates (n = 171) were evaluated for stx1 and stx2 using different primer sets. For this, approximately half of the PCR were performed using original DNA template provided by LAB A and half using DNA extracted from sub-cultured isolates. All Top 7 isolates were sub-cultured by LAB A and shipped to LAB C for traditional serotyping (TS) to determine O and H groups, with PCR-confirmation of virulence genes using a third set of primers. By TS, 76% of O groups (525/689) matched PCR-determined O groups. Lowest proportions (p < 0.05) of O group matches between PCR and TS (62.6% and 69.8%) occurred for O26 and O45 serogroups, respectively. PCR-detection of stx differed most between LAB A and LAB C. Excluding isolates where O groups by PCR and TS did not match, detection of stx1 was most consistent (p < 0.01) for O111 and O157:H7/NM. In contrast, for O45 and O103, stx1 was detected in >65% of isolates by LAB A and <5% by LAB C. Stx2 was only detected by LAB C in isolates of serogroups O121, O145, and O157:H7/NM. LAB B also detected stx2 in O26 and O157:H12/H29, while LAB A detected stx2 in all serogroups. Excluding O111 and O157:H7/NM, marked changes in stx detection were observed between initial isolation and sub-cultures of the same isolate. While multiple explanations exist for discordant O-typing between PCR and TS and for differences in stx detection across labs, these data suggest that assays for STEC classification may require re-evaluation and/or standardization.

Genotypic Features of Clinical and Bovine Escherichia coli O157 Strains Isolated in Countries with Different Associated-Disease Incidences

There is great geographical variation in the frequency of Escherichia coli O157 infections that correlates with important differences in the bovine reservoir of each country. Our group carried out a broad molecular characterization of human and bovine E. coli O157 strains circulating in Argentina using different methodologies. Our data allows us to conclude that in Argentina, a high homogeneity is observed in both cattle and human strains, with almost exclusive circulation of strains belonging to the hypervirulent clade 8 described by Manning. The aim of this review was to compare the genetic background of E. coli O157 strains isolated in countries that have conducted similar studies, to try to correlate specific O157 genotypes with the incidence and severity of E. coli O157 associated diseases. The characteristics of the strains that cause disease in humans reflect the predominant genotypes in cattle in each of the countries analyzed. The main features clearly linked to high incidence or severity of E. coli O157 infections are lineage-specific polymorphism assay-6 lineage I/II, clade 8 strains and probably, clade 6 strains, the stx_2a/stx_2c genotype, the presence of q933 and q21 simultaneously, and putative virulence factor EC_3286. In countries with an absence of these features in O157 strains, the overall incidence of O157 disease is low. Argentina, where these characteristics are detected in most strains, shows the highest incidence of hemolytic uremic syndrome (HUS) worldwide.

The role of meat in foodborne disease: Is there a coming revolution in risk assessment and management?

Meat has featured prominently as a source of foodborne disease and a public health concern. For about the past 20 years the risk management paradigm has dominated international thinking about food safety. Control through the supply chain is supported by risk management concepts, as the public health risk at the point of consumption becomes the accepted outcome based measure. Foodborne pathogens can be detected at several points in the supply chain and determining the source of where these pathogens arise and how they behave throughout meat production and processing are important parts of risk based approaches. Recent improvements in molecular and genetic based technologies and data analysis for investigating source attribution and pathogen behaviour have enabled greater insights into how foodborne outbreaks occur and where controls can be implemented. These new approaches will improve our understanding of the role of meat in foodborne disease and are expected to have a significant impact on our understanding in the coming years.

Virulence profile of Escherichia coli O157 strains isolated from surface water in cattle breeding areas

Shiga toxin-producing Escherichia coli (STEC) O157:H7 is a worldwide concern. Cattle are their main reservoir and may contaminate watercourses through manure. We characterized a collection of 38 STEC O157:H7 strains isolated from surface water in feedlots areas (puddles inside pens formed after the rainfall or by spill around drinking troughs, and small water courses and lagoons, formed by runoff). Nineteen (50·0%) strains harboured stx_2a /stx_2c genes, 18 (47·4%) stx_2c and one stx_1a /stx_2c . All strains harboured eae, ehxA, rfb_O157 and fliC_H7 genes, and the putative virulence determinants ECSP_0242, ECSP_2687 and ECSP_3620. All isolates tested as Lineage I/II by lineage-specific polymorphism assay-6. Nineteen (50%) belonged to the high virulent clade 8. The q₂₁ allele was found in all strains and q₉₃₃ /q₂₁ alleles in 17 (44·7%). By XbaI-pulsed-field gel electrophoresis, 29 strains were grouped into seven clusters. Four clusters grouped isolates from distant places separated by 150-250 km. This may be related to vectors, like birds, involved in their spread. Otherwise, three clusters contained isolates recovered at same places with intervals of 1-9 months. This could be explained by the high environmental persistence of STEC O157:H7. These strains recovered from surface water showed similar genotypes to those found in the bovine reservoir and in human diseases, and could be linked to the high incidence of haemolytic uremic syndrome in Argentina.

SIGNIFICANCE AND IMPACT OF THE STUDY

The challenge for the growing global demand for food is to find sustained production strategies without collateral effects. Intensive livestock operations generate large volumes of manure that can contaminate a finite resource, the water. This study shows how water contaminated by confined feeding operations can transport dangerous pathogens and warns to pay more attention to control and sanitation systems to prevent this type of pollution.

Shiga Toxin-Producing Escherichia coli O157 Shedding Dynamics in an Australian Beef Herd

Shiga toxin-producing Escherichia coli (STEC) O157 is an important foodborne pathogen that can be transmitted to humans both directly and indirectly from the feces of beef cattle, its primary reservoir. Numerous studies have investigated the shedding dynamics of E. coli O157 by beef cattle; however, the spatiotemporal trends of shedding are still not well understood. Molecular tools can increase the resolution through the use of strain typing to explore transmission dynamics within and between herds and identify strain-specific characteristics that may influence pathogenicity and spread. Previously, the shedding dynamics and molecular diversity, through the use of multilocus variable number of tandem repeat analysis (MLVA) of STEC O157, were separately investigated in an Australian beef herd over a 9-month study period. Variation in shedding was observed over time, and 33 MLVA types were identified. The study presented here combines the two datasets previously published with an aim to clarify the relationship between epidemiological variables and strain types. Three major genetic clusters (GCs) were identified that were significantly associated with the location of the cattle in different paddocks. No significant association between GCs and individual cow was observed. Results from this molecular epidemiological study provide evidence for herd-level clonal replacement over time that may have been triggered by movement to a new paddock. In conclusion, this study has provided further insight into STEC O157 shedding dynamics and pathogen transmission. Knowledge gaps remain regarding the relationship of strain types and the shedding dynamics of STEC O157 by beef cattle that could be further clarified through the use of whole-genome sequencing.

Application of whole genome sequence data in analyzing the molecular epidemiology of Shiga toxin-producing Escherichia coli O157:H7/H

Seventeen clusters of Shiga toxin-producing Escherichia coli O157:H7/- (O157) strains, determined by cluster analysis of pulsed-field gel electrophoresis patterns, were analyzed using whole genome sequence (WGS) data to investigate this pathogen's molecular epidemiology. The 17 clusters included 136 strains containing strains from nine outbreaks, with each outbreak caused by a single source contaminated with the organism, as shown by epidemiological contact surveys. WGS data of these strains were used to identify single nucleotide polymorphisms (SNPs) by two methods: short read data were directly mapped to a reference genome (mapping derived SNPs) and common SNPs between the mapping derived SNPs and SNPs in assembled data of short read data (common SNPs). Among both SNPs, those that were detected in genes with a gap were excluded to remove ambiguous SNPs from further analysis. The effectiveness of both SNPs was investigated among all the concatenated SNPs that were detected (whole SNP set); SNPs were divided into three categories based on the genes in which they were located (i.e., backbone SNP set, O-island SNP set, and mobile element SNP set); and SNPs in non-coding regions (intergenic region SNP set). When SNPs from strains isolated from the nine single source derived outbreaks were analyzed using an unweighted pair group method with arithmetic mean tree (UPGMA) and a minimum spanning tree (MST), the maximum pair-wise distances of the backbone SNP set of the mapping derived SNPs were significantly smaller than those of the whole and intergenic region SNP set on both UPGMAs and MSTs. This significant difference was also observed when the backbone SNP set of the common SNPs were examined (Steel-Dwass test, P≤0.01). When the maximum pair-wise distances were compared between the mapping derived and common SNPs, significant differences were observed in those of the whole, mobile element, and intergenic region SNP set (Wilcoxon signed rank test, P≤0.01). When all the strains included in one complex on an MST or one cluster on a UPGMA were designated as the same genotype, the values of the Hunter-Gaston Discriminatory Power Index for the backbone SNP set of the mapping derived and common SNPs were higher than those of other SNP sets. In contrast, the mobile element SNP set could not robustly subdivide lineage I strains of tested O157 strains using both the mapping derived and common SNPs. These results suggested that the backbone SNP set were the most effective for analysis of WGS data for O157 in enabling an appropriation of its molecular epidemiology.

Whole-Genome Sequencing of Drug-Resistant Salmonella enterica Isolates from Dairy Cattle and Humans in New York and Washington States Reveals Source and Geographic Associations

Multidrug-resistant (MDR) Salmonella enterica can be spread from cattle to humans through direct contact with animals shedding Salmonella as well as through the food chain, making MDR Salmonella a serious threat to human health. The objective of this study was to use whole-genome sequencing to compare antimicrobial-resistant (AMR) Salmonella enterica serovars Typhimurium, Newport, and Dublin isolated from dairy cattle and humans in Washington State and New York State at the genotypic and phenotypic levels. A total of 90 isolates were selected for the study (37 S Typhimurium, 32 S Newport, and 21 S Dublin isolates). All isolates were tested for phenotypic antibiotic resistance to 12 drugs using Kirby-Bauer disk diffusion. AMR genes were detected in the assembled genome of each isolate using nucleotide BLAST and ARG-ANNOT. Genotypic prediction of phenotypic resistance resulted in a mean sensitivity of 97.2 and specificity of 85.2. Sulfamethoxazole-trimethoprim resistance was observed only in human isolates (P < 0.05), while resistance to quinolones and fluoroquinolones was observed only in 6 S Typhimurium isolates from humans in Washington State. S Newport isolates showed a high degree of AMR profile similarity, regardless of source. S Dublin isolates from New York State differed from those from Washington State based on the presence/absence of plasmid replicons, as well as phenotypic AMR susceptibility/nonsusceptibility (P < 0.05). The results of this study suggest that distinct factors may contribute to the emergence and dispersal of AMR S. enterica in humans and farm animals in different regions.IMPORTANCE The use of antibiotics in food-producing animals has been hypothesized to select for AMR Salmonella enterica and associated AMR determinants, which can be transferred to humans through different routes. Previous studies have sought to assess the degree to which AMR livestock- and human-associated Salmonella strains overlap, as well as the spatial distribution of Salmonella's associated AMR determinants, but have often been limited by the degree of resolution at which isolates can be compared. Here, a comparative genomics study of livestock- and human-associated Salmonella strains from different regions of the United States shows that while many AMR genes and phenotypes were confined to human isolates, overlaps between the resistomes of bovine and human-associated Salmonella isolates were observed on numerous occasions, particularly for S Newport. We have also shown that whole-genome sequencing can be used to reliably predict phenotypic resistance across Salmonella isolated from bovine sources.

Whole Genome Sequencing for Genomics-Guided Investigations of Escherichia coli O157:H7 Outbreaks

Multi isolate whole genome sequencing (WGS) and typing for outbreak investigations has become a reality in the post-genomics era. We applied this technology to strains from Escherichia coli O157:H7 outbreaks. These include isolates from seven North America outbreaks, as well as multiple isolates from the same patient and from different infected individuals in the same household. Customized high-resolution bioinformatics sequence typing strategies were developed to assess the core genome and mobilome plasticity. Sequence typing was performed using an in-house single nucleotide polymorphism (SNP) discovery and validation pipeline. Discriminatory power becomes of particular importance for the investigation of isolates from outbreaks in which macrogenomic techniques such as pulse-field gel electrophoresis or multiple locus variable number tandem repeat analysis do not differentiate closely related organisms. We also characterized differences in the phage inventory, allowing us to identify plasticity among outbreak strains that is not detectable at the core genome level. Our comprehensive analysis of the mobilome identified multiple plasmids that have not previously been associated with this lineage. Applied phylogenomics approaches provide strong molecular evidence for exceptionally little heterogeneity of strains within outbreaks and demonstrate the value of intra-cluster comparisons, rather than basing the analysis on archetypal reference strains. Next generation sequencing and whole genome typing strategies provide the technological foundation for genomic epidemiology outbreak investigation utilizing its significantly higher sample throughput, cost efficiency, and phylogenetic relatedness accuracy. These phylogenomics approaches have major public health relevance in translating information from the sequence-based survey to support timely and informed countermeasures. Polymorphisms identified in this work offer robust phylogenetic signals that index both short- and long-term evolution and can complement currently employed typing schemes for outbreak ex- and inclusion, diagnostics, surveillance, and forensic studies.