Campylobacter californiensis sp. nov., isolated from cattle and feral swine.

Nine Campylobacter strains were isolated from cattle and feral swine faeces: three were recovered during a 2007 Campylobacter-associated outbreak linked to a dairy, and the other six were isolated during a 2009-2010 survey of farms and ranches in Central California. The species identification of these strains could not be determined by 16S rRNA gene sequencing but were most similar to Campylobacter concisus and Campylobacter mucosalis. Additional atpA typing indicated that the nine strains composed a discrete novel clade related to C. concisus and C. mucosalis. A polyphasic study was undertaken here to clarify their taxonomic position. Phylogenetic analyses were performed based on 16S rRNA gene sequences and the concatenated sequences of 330 core genes. The core gene analysis placed the nine strains into a clade well separated from the other Campylobacter taxa, indicating that these strains represent a novel Campylobacter species. Pairwise digital DNA-DNA hybridization and average nucleotide identity values between these strains and other campylobacters are lower than 16 and 73%, respectively, further supporting their placement into a novel taxon. Standard phenotypic testing was performed. All strains are microaerobic or anaerobic, motile, Gram-negative, slightly-curved rods that are oxidase positive but catalase negative. Strains can be distinguished from the other catalase-negative Campylobacter species using phenotypic markers such as motility, oxidase activity, cephalothin resistance, hippuricase activity, growth at 30 °C, and α-haemolysis. The data presented here show that these strains represent a novel species within Campylobacter, for which the name Campylobacter californiensis sp. nov. (type strain RM6914T=LMG 32304T=CCUG 75329T) is proposed.

Rates of evolutionary change of resident Escherichia coli O157:H7 differ within the same ecological niche.

BACKGROUND: Shiga toxin-producing Escherichia coli (STEC) O157:H7 is a pathogen known to reside in cattle feedlots. This retrospective study examined 181 STEC O157:H7 strains collected over 23 years from a closed-system feedlot. All strains were subjected to short-read sequencing, with a subset of 36 also subjected to long-read sequencing. RESULTS: Over 96% of the strains fell into four phylogenetically distinct clades. Clade membership was associated with multiple factors including stx composition and the alleles of a well-characterized polymorphism (tir 255 T > A). Small plasmids (2.7 to 40 kb) were found to be primarily clade specific. Within each clade, chromosomal rearrangements were observed along with a core phageome and clade specific phages. Across both core and mobile elements of the genome, multiple SNP alleles were in complete linkage disequilibrium across all strains within specific clades. Clade evolutionary rates varied between 0.9 and 2.8 SNP/genome/year with two tir A allele clades having the lowest evolutionary rates. Investigation into possible causes of the differing rates was not conclusive but revealed a synonymous based mutation in the DNA polymerase III of the fastest evolving clade. Phylogenetic trees generated through our bioinformatic pipeline versus the NCBI's pathogen detection project were similar, with the two tir A allele clades matching individual NCBI SNP clusters, and the two tir T allele clades assigned to multiple closely-related SNP clusters. CONCLUSIONS: In one ecological niche, a diverse STEC O157:H7 population exhibited different rates of evolution that associated with SNP alleles in linkage disequilibrium in the core genome and mobile elements, including tir 255 T > A.

Whole genome sequencing of Moraxella bovis strains from North America reveals two genotypes with different genetic determinants.

BACKGROUND: Moraxella bovis and Moraxella bovoculi both associate with infectious bovine keratoconjunctivitis (IBK), an economically significant and painful ocular disease that affects cattle worldwide. There are two genotypes of M. bovoculi (genotypes 1 and 2) that differ in their gene content and potential virulence factors, although neither have been experimentally shown to cause IBK. M. bovis is a causative IBK agent, however, not all strains carry a complete assortment of known virulence factors. The goals of this study were to determine the population structure and depth of M. bovis genomic diversity, and to compare core and accessory genes and predicted outer membrane protein profiles both within and between M. bovis and M. bovoculi. RESULTS: Phylogenetic trees and bioinformatic analyses of 36 M. bovis chromosomes sequenced in this study and additional available chromosomes of M. bovis and both genotype 1 and 2 M. bovoculi, showed there are two genotypes (1 and 2) of M. bovis. The two M. bovis genotypes share a core of 2015 genes, with 121 and 186 genes specific to genotype 1 and 2, respectively. The two genotypes differ by their chromosome size and prophage content, encoded protein variants of the virulence factor hemolysin, and by their affiliation with different plasmids. Eight plasmid types were identified in this study, with types 1 and 6 observed in 88 and 56% of genotype 2 strains, respectively, and absent from genotype 1 strains. Only type 1 plasmids contained one or two gene copies encoding filamentous haemagglutinin-like proteins potentially involved with adhesion. A core of 1403 genes was shared between the genotype 1 and 2 strains of both M. bovis and M. bovoculi, which encoded a total of nine predicted outer membrane proteins. CONCLUSIONS: There are two genotypes of M. bovis that differ in both chromosome content and plasmid profiles and thus may not equally associate with IBK. Immunological reagents specifically targeting select genotypes of M. bovis, or all genotypes of M. bovis and M. bovoculi together could be designed from the outer membrane proteins identified in this study.

Shiga toxin sub-type 2a increases the efficiency of Escherichia coli O157 transmission between animals and restricts epithelial regeneration in bovine enteroids.

Specific Escherichia coli isolates lysogenised with prophages that express Shiga toxin (Stx) can be a threat to human health, with cattle being an important natural reservoir. In many countries the most severe pathology is associated with enterohaemorrhagic E. coli (EHEC) serogroups that express Stx subtype 2a. In the United Kingdom, phage type (PT) 21/28 O157 strains have emerged as the predominant cause of life-threatening EHEC infections and this phage type commonly encodes both Stx2a and Stx2c toxin types. PT21/28 is also epidemiologically linked to super-shedding (>103 cfu/g of faeces) which is significant for inter-animal transmission and human infection as demonstrated using modelling studies. We demonstrate that Stx2a is the main toxin produced by stx2a+/stx2c+ PT21/28 strains induced with mitomycin C and this is associated with more rapid induction of gene expression from the Stx2a-encoding prophage compared to that from the Stx2c-encoding prophage. Bacterial supernatants containing either Stx2a and/or Stx2c were demonstrated to restrict growth of bovine gastrointestinal organoids with no restriction when toxin production was not induced or prevented by mutation. Isogenic strains that differed in their capacity to produce Stx2a were selected for experimental oral colonisation of calves to assess the significance of Stx2a for both super-shedding and transmission between animals. Restoration of Stx2a expression in a PT21/28 background significantly increased animal-to-animal transmission and the number of sentinel animals that became super-shedders. We propose that while both Stx2a and Stx2c can restrict regeneration of the epithelium, it is the relatively rapid and higher levels of Stx2a induction, compared to Stx2c, that have contributed to the successful emergence of Stx2a+ E. coli isolates in cattle in the last 40 years. We propose a model in which Stx2a enhances E. coli O157 colonisation of in-contact animals by restricting regeneration and turnover of the colonised gastrointestinal epithelium.

Locus of Heat Resistance (LHR) in Meat-Borne Escherichia coli: Screening and Genetic Characterization.

Microbial resistance to processing treatments poses a food safety concern, as treatment tolerant pathogens can emerge. Occasional foodborne outbreaks caused by pathogenic Escherichia coli have led to human and economic losses. Therefore, this study screened for the extreme heat resistance (XHR) phenotype as well as one known genetic marker, the locus of heat resistance (LHR), in 4,123 E. coli isolates from diverse meat animals at different processing stages. The prevalences of XHR and LHR among the meat-borne E. coli were found to be 10.3% and 11.4%, respectively, with 19% agreement between the two. Finished meat products showed the highest LHR prevalence (24.3%) compared to other processing stages (0 to 0.6%). None of the LHR+E. coli in this study would be considered pathogens based on screening for virulence genes. Four high-quality genomes were generated by whole-genome sequencing of representative LHR+ isolates. Nine horizontally acquired LHRs were identified and characterized, four plasmid-borne and five chromosomal. Nine newly identified LHRs belong to ClpK1 LHR or ClpK2 LHR variants sharing 61 to 68% nucleotide sequence identity, while one LHR appears to be a hybrid. Our observations suggest positive correlation between the number of LHR regions present in isolates and the extent of heat resistance. The isolate exhibiting the highest degree of heat resistance possessed four LHRs belonging to three different variant groups. Maintenance of as many as four LHRs in a single genome emphasizes the benefits of the LHR in bacterial physiology and stress response.IMPORTANCE Currently, a "multiple-hurdle" approach based on a combination of different antimicrobial interventions, including heat, is being utilized during meat processing to control the burden of spoilage and pathogenic bacteria. Our recent study (M. Guragain, G. E. Smith, D. A. King, and J. M. Bosilevac, J Food Prot 83:1438-1443, 2020, https://doi.org/10.4315/JFP-20-103) suggests that U.S. beef cattle harbor Escherichia coli that possess the locus of heat resistance (LHR). LHR seemingly contributes to the global stress tolerance in bacteria and hence poses a food safety concern. Therefore, it is important to understand the distribution of the LHRs among meat-borne bacteria identified at different stages of different meat processing systems. Complete genome sequencing and comparative analysis of selected heat-resistant bacteria provide a clearer understanding of stress and heat resistance mechanisms. Further, sequencing data may offer a platform to gain further insights into the genetic background that provides optimal bacterial tolerance against heat and other processing treatments.

Rapid estimation of Salmonella enterica contamination level in ground beef - Application of the time-to-positivity method using a combination of molecular detection and direct plating.

Little progress has been made in decreasing the incidence rate of salmonellosis in the US over the past decade. Mitigating the contribution of contaminated raw meat to the salmonellosis incidence rate requires rapid methods for quantifying Salmonella, so that highly contaminated products can be removed before entering the food chain. Here we evaluated the use of Time-to-Positivity (TTP) as a rapid, semi-quantitative approach for estimating Salmonella contamination levels in ground beef. Growth rates of 14 Salmonella strains (inoculated at log 1 to -2 CFU/g) were characterized in lean ground beef mTSB enrichments and time-to-detection was determined using culture and molecular detection methods. Enrichments were sampled at five timepoints and results were used to construct a prediction model of estimated contamination level by TTP (superscript indicates time in hours) defined as TTP4: ≥5 CFU/g; TTP6: ≤5, ≥1 CFU/g; TTP8: ≤1, ≥0.01 CFU/g; with samples negative at 8 h estimated ≤0.01 CFU/g. Model performance measures showed high sensitivity (100%) and specificity (83% and 93% for two detection methods) for samples with a TTP4, with false negative rates of 0%.

Differences between predicted outer membrane proteins of genotype 1 and 2 Mannheimia haemolytica.

BACKGROUND: Mannheimia haemolytica strains isolated from North American cattle have been classified into two genotypes (1 and 2). Although members of both genotypes have been isolated from the upper and lower respiratory tracts of cattle with or without bovine respiratory disease (BRD), genotype 2 strains are much more frequently isolated from diseased lungs than genotype 1 strains. The mechanisms behind the increased association of genotype 2 M. haemolytica with BRD are not fully understood. To address that, and to search for interventions against genotype 2 M. haemolytica, complete, closed chromosome assemblies for 35 genotype 1 and 34 genotype 2 strains were generated and compared. Searches were conducted for the pan genome, core genes shared between the genotypes, and for genes specific to either genotype. Additionally, genes encoding outer membrane proteins (OMPs) specific to genotype 2 M. haemolytica were identified, and the diversity of their protein isoforms was characterized with predominantly unassembled, short-read genomic sequences for up to 1075 additional strains. RESULTS: The pan genome of the 69 sequenced M. haemolytica strains consisted of 3111 genes, of which 1880 comprised a shared core between the genotypes. A core of 112 and 179 genes or gene variants were specific to genotype 1 and 2, respectively. Seven genes encoding predicted OMPs; a peptidase S6, a ligand-gated channel, an autotransporter outer membrane beta-barrel domain-containing protein (AOMB-BD-CP), a porin, and three different trimeric autotransporter adhesins were specific to genotype 2 as their genotype 1 homologs were either pseudogenes, or not detected. The AOMB-BD-CP gene, however, appeared to be truncated across all examined genotype 2 strains and to likely encode dysfunctional protein. Homologous gene sequences from additional M. haemolytica strains confirmed the specificity of the remaining six genotype 2 OMP genes and revealed they encoded low isoform diversity at the population level. CONCLUSION: Genotype 2 M. haemolytica possess genes encoding conserved OMPs not found intact in more commensally prone genotype 1 strains. Some of the genotype 2 specific genes identified in this study are likely to have important biological roles in the pathogenicity of genotype 2 M. haemolytica, which is the primary bacterial cause of BRD.

Genomic-based identification of environmental and clinical Listeria monocytogenes strains associated with an abortion outbreak in beef heifers.

BACKGROUND: In a beef cattle facility an outbreak of abortions occurred over a 36-day period and included samples from two aborted (non-viable) fetuses and 21 post-abortion clinical cases. There are numerous etiologies, including clinical listeriosis. At the species level, Listeria monocytogenes is ubiquitous in cattle production environments, including soil, feed, and occasionally water sources, and is a common enteric resident of cattle and other mammals. There are four genetically distinct lineages of L. monocytogenes (I-IV), with most lineage III and IV isolates obtained from ruminants. Definitive diagnosis of L. monocytogenes as a causative agent in disease outbreaks relies upon case identification, appropriate sample collection, and laboratory confirmation. Furthermore, clearly establishing a relationship between a pathogen source and clinical disease is difficult. RESULTS: Of the two fetal and 21 clinical case submissions, 19 were positive for L. monocytogenes. Subsequent culture for L. monocytogenes from water and silage sources identified both as potential origins of infection. Using whole-genome sequencing and phylogenetic analyses, clinical, water and silage L. monocytogenes strains grouped into two of four lineages. All water and silage strains, plus 11 clinical strains placed in lineage III, with identical or nearly identical genomic sequences. The remaining eight clinical strains placed in lineage I, with seven having nearly identical sequences and one distinctly different. CONCLUSION: Three genetically distinct strains within two lineages of L. monocytogenes caused the abortion outbreak. The etiology of abortion in 11 cases was directly linked to water and silage contamination from a lineage III L. monocytogenes strain. The source of infection for the remaining abortion cases with two different strains from lineage I is unknown. This is the first report of L. monocytogenes genomics being used as part of an outbreak investigation of cattle abortion.

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.

Support vector machine applied to predict the zoonotic potential of E. coli O157 cattle isolates.

Sequence analyses of pathogen genomes facilitate the tracking of disease outbreaks and allow relationships between strains to be reconstructed and virulence factors to be identified. However, these methods are generally used after an outbreak has happened. Here, we show that support vector machine analysis of bovine E. coli O157 isolate sequences can be applied to predict their zoonotic potential, identifying cattle strains more likely to be a serious threat to human health. Notably, only a minor subset (less than 10%) of bovine E. coli O157 isolates analyzed in our datasets were predicted to have the potential to cause human disease; this is despite the fact that the majority are within previously defined pathogenic lineages I or I/II and encode key virulence factors. The predictive capacity was retained when tested across datasets. The major differences between human and bovine E. coli O157 isolates were due to the relative abundances of hundreds of predicted prophage proteins. This finding has profound implications for public health management of disease because interventions in cattle, such a vaccination, can be targeted at herds carrying strains of high zoonotic potential. Machine-learning approaches should be applied broadly to further our understanding of pathogen biology.

Genomic signatures of Mannheimia haemolytica that associate with the lungs of cattle with respiratory disease, an integrative conjugative element, and antibiotic resistance genes.

BACKGROUND: Mannheimia haemolytica typically resides in cattle as a commensal member of the upper respiratory tract microbiome. However, some strains can invade their lungs and cause respiratory disease and death, including those with multi-drug resistance. A nucleotide polymorphism typing system was developed for M. haemolytica from the genome sequences of 1133 North American isolates, and used to identify genetic differences between isolates from the lungs and upper respiratory tract of cattle with and without clinical signs of respiratory disease. RESULTS: A total of 26,081 nucleotide polymorphisms were characterized after quality control filtering of 48,403 putative polymorphisms. Phylogenetic analyses of nucleotide polymorphism genotypes split M. haemolytica into two major genotypes (1 and 2) that each were further divided into multiple subtypes. Multiple polymorphisms were identified with alleles that tagged genotypes 1 or 2, and their respective subtypes. Only genotype 2 M. haemolytica associated with the lungs of diseased cattle and the sequence of a particular integrative and conjugative element (ICE). Additionally, isolates belonging to one subtype of genotype 2 (2b), had the majority of antibiotic resistance genes detected in this study, which were assorted into seven combinations that ranged from 1 to 12 resistance genes. CONCLUSIONS: Typing of diverse M. haemolytica by nucleotide polymorphism genotypes successfully identified associations with diseased cattle lungs, ICE sequence, and antibiotic resistance genes. Management of cattle by their carriage of M. haemolytica could be an effective intervention strategy to reduce the prevalence of respiratory disease and supplemental needs for antibiotic treatments in North American herds.

Large genomic differences between Moraxella bovoculi isolates acquired from the eyes of cattle with infectious bovine keratoconjunctivitis versus the deep nasopharynx of asymptomatic cattle.

Moraxella bovoculi is a recently described bacterium that is associated with infectious bovine keratoconjunctivitis (IBK) or "pinkeye" in cattle. In this study, closed circularized genomes were generated for seven M. bovoculi isolates: three that originated from the eyes of clinical IBK bovine cases and four from the deep nasopharynx of asymptomatic cattle. Isolates that originated from the eyes of IBK cases profoundly differed from those that originated from the nasopharynx of asymptomatic cattle in genome structure, gene content and polymorphism diversity and consequently placed into two distinct phylogenetic groups. These results suggest that there are genetically distinct strains of M. bovoculi that may not associate with IBK.

Effect of proximity to a cattle feedlot on Escherichia coli O157:H7 contamination of leafy greens and evaluation of the potential for airborne transmission.

The impact of proximity to a beef cattle feedlot on Escherichia coli O157:H7 contamination of leafy greens was examined. In each of 2 years, leafy greens were planted in nine plots located 60, 120, and 180 m from a cattle feedlot (3 plots at each distance). Leafy greens (270) and feedlot manure samples (100) were collected six different times from June to September in each year. Both E. coli O157:H7 and total E. coli bacteria were recovered from leafy greens at all plot distances. E. coli O157:H7 was recovered from 3.5% of leafy green samples per plot at 60 m, which was higher (P < 0.05) than the 1.8% of positive samples per plot at 180 m, indicating a decrease in contamination as distance from the feedlot was increased. Although E. coli O157:H7 was not recovered from air samples at any distance, total E. coli was recovered from air samples at the feedlot edge and all plot distances, indicating that airborne transport of the pathogen can occur. Results suggest that risk for airborne transport of E. coli O157:H7 from cattle production is increased when cattle pen surfaces are very dry and when this situation is combined with cattle management or cattle behaviors that generate airborne dust. Current leafy green field distance guidelines of 120 m (400 feet) may not be adequate to limit the transmission of E. coli O157:H7 to produce crops planted near concentrated animal feeding operations. Additional research is needed to determine safe set-back distances between cattle feedlots and crop production that will reduce fresh produce contamination.

Genome sequencing and comparative genomics provides insights on the evolutionary dynamics and pathogenic potential of different H-serotypes of Shiga toxin-producing Escherichia coli O104.

BACKGROUND: Various H-serotypes of the Shiga toxin-producing Escherichia coli (STEC) O104, including H4, H7, H21, and H¯, have been associated with sporadic cases of illness and have caused food-borne outbreaks globally. In the U.S., STEC O104:H21 caused an outbreak associated with milk in 1994. However, there is little known on the evolutionary origins of STEC O104 strains, and how genotypic diversity contributes to pathogenic potential of various O104 H-antigen serotypes isolated from different ecological niches and/or geographical regions. RESULTS: Two STEC O104:H21 (milk outbreak strain) and O104:H7 (cattle isolate) strains were shot-gun sequenced, and the genomes were closed. The intimin (eae) gene, involved in the attaching-effacing phenotype of diarrheagenic E. coli, was not found in either strain. Examining various O104 genome sequences, we found that two "complete" left and right end portions of the locus of enterocyte effacement (LEE) pathogenicity island were present in 13 O104 strains; however, the central portion of LEE was missing, where the eae gene is located. In O104:H4 strains, the missing central portion of the LEE locus was replaced by a pathogenicity island carrying the aidA (adhesin involved in diffuse adherence) gene and antibiotic resistance genes commonly carried on plasmids. Enteroaggregative E. coli-specific virulence genes and European outbreak O104:H4-specific stx2-encoding Escherichia P13374 or Escherichia TL-2011c bacteriophages were missing in some of the O104:H4 genome sequences available from public databases. Most of the genomic variations in the strains examined were due to the presence of different mobile genetic elements, including prophages and genomic island regions. The presence of plasmids carrying virulence-associated genes may play a role in the pathogenic potential of O104 strains. CONCLUSIONS: The two strains sequenced in this study (O104:H21 and O104:H7) are genetically more similar to each other than to the O104:H4 strains that caused an outbreak in Germany in 2011 and strains found in Central Africa. A hypothesis on strain evolution and pathogenic potential of various H-serotypes of E. coli O104 strains is proposed.

Characterization of Escherichia coli O157:H7 strains from contaminated raw beef trim during "high event periods".

The development and implementation of effective antimicrobial interventions by the beef processing industry in the United States have dramatically reduced the incidence of beef trim contamination by Escherichia coli O157:H7. However, individual processing plants still experience sporadic peaks in contamination rates where multiple E. coli O157:H7-positive lots are clustered in a short time frame. These peaks have been referred to as "high event periods" (HEP) of contamination. The results reported here detail the characterization of E. coli O157:H7 isolates from 21 HEP across multiple companies and processing plants to gain insight regarding the mechanisms causing these incidents. Strain genotypes were determined by pulsed-field gel electrophoresis, and isolates were investigated for characteristics linking them to human illness. Through these analyses, it was determined that individual HEP show little to no diversity in strain genotypes. Hence, each HEP has one strain type that makes up most, if not all, of the contamination. This is shown to differ from the genotypic diversity of E. coli O157:H7 found on the hides of cattle entering processing plants. In addition, it was found that a large proportion (81%) of HEP are caused by strain types associated with human illness. These results pose a potential challenge to the current model for finished product contamination during beef processing.

The GEA pipeline for characterizing Escherichia coli and Salmonella genomes.

Salmonella enterica and Escherichia coli are major food-borne human pathogens, and their genomes are routinely sequenced for clinical surveillance. Computational pipelines designed for analyzing pathogen genomes should both utilize the most current information from annotation databases and increase the coverage of these databases over time. We report the development of the GEA pipeline to analyze large batches of E. coli and S. enterica genomes. The GEA pipeline takes as input paired Illumina raw reads files which are then assembled followed by annotation. Alternatively, assemblies can be provided as input and directly annotated. The pipeline provides predictive genome annotations for E. coli and S. enterica with a focus on the Center for Genomic Epidemiology tools. Annotation results are provided as a tab delimited text file. The GEA pipeline is designed for large-scale E. coli and S. enterica genome assembly and characterization using the Center for Genomic Epidemiology command-line tools and high-performance computing. Large scale annotation is demonstrated by an analysis of more than 14,000 Salmonella genome assemblies. Testing the GEA pipeline on E. coli raw reads demonstrates reproducibility across multiple compute environments and computational usage is optimized on high performance computers.

Pathogenomes and virulence profiles of representative big six non-O157 serogroup Shiga toxin-producing Escherichia coli.

Shiga toxin (Stx)-producing Escherichia coli (STEC) of non-O157:H7 serotypes are responsible for global and widespread human food-borne disease. Among these serogroups, O26, O45, O103, O111, O121, and O145 account for the majority of clinical infections and are colloquially referred to as the "Big Six." The "Big Six" strain panel we sequenced and analyzed in this study are reference type cultures comprised of six strains representing each of the non-O157 STEC serogroups curated and distributed by the American Type Culture Collection (ATCC) as a resource to the research community under panel number ATCC MP-9. The application of long- and short-read hybrid sequencing yielded closed chromosomes and a total of 14 plasmids of diverse functions. Through high-resolution comparative phylogenomics, we cataloged the shared and strain-specific virulence and resistance gene content and established the close relationship of serogroup O26 and O103 strains featuring flagellar H-type 11. Virulence phenotyping revealed statistically significant differences in the Stx-production capabilities that we found to be correlated to the strain's individual stx-status. Among the carried Stx1a, Stx2a, and Stx2d phages, the Stx2a phage is by far the most responsive upon RecA-mediated phage mobilization, and in consequence, stx2a + isolates produced the highest-level of toxin in this panel. The availability of high-quality closed genomes for this "Big Six" reference set, including carried plasmids, along with the recorded genomic virulence profiles and Stx-production phenotypes will provide a valuable foundation to further explore the plasticity in evolutionary trajectories in these emerging non-O157 STEC lineages, which are major culprits of human food-borne disease.

Pathogenomes of Shiga Toxin Positive and Negative Escherichia coli O157:H7 Strains TT12A and TT12B: Comprehensive Phylogenomic Analysis Using Closed Genomes.

Shiga toxin-producing Escherichia coli are zoonotic pathogens that cause food-borne human disease. Among these, the O157:H7 serotype has evolved from an enteropathogenic O55:H7 ancestor through the displacement of the somatic gene cluster and recurrent toxigenic conversion by Shiga toxin-converting bacteriophages. However, atypical strains that lack the Shiga toxin, the characteristic virulence hallmark, are circulating in this lineage. For this study, we analyzed the pathogenome and virulence inventories of the stx+ strain, TT12A, isolated from a patient with hemorrhagic colitis, and its respective co-isolated stx- strain, TT12B. Sequencing the genomes to closure proved critical to the cataloguing of subtle strain differentiating sequence and structural polymorphisms at a high-level of phylogenetic accuracy and resolution. Phylogenomic profiling revealed SNP and MLST profiles similar to the near clonal outbreak isolates. Their prophage inventories, however, were notably different. The attenuated atypical non-shigatoxigenic status of TT12B is explained by the absence of both the ΦStx1a- and ΦStx2a-prophages carried by TT12A, and we also recorded further alterations in the non-Stx prophage complement. Phenotypic characterization indicated that culture growth was directly impacted by the strains' distinct lytic phage complement. Altogether, our phylogenomic and phenotypic analyses show that these intimately related isogenic strains are on divergent Stx(+/stx-) evolutionary paths.

Phylogeny of Shiga toxin-producing Escherichia coli O157 isolated from cattle and clinically ill humans.

Cattle are a major reservoir for Shiga toxin-producing Escherichia coli O157 (STEC O157) and harbor multiple genetic subtypes that do not all associate with human disease. STEC O157 evolved from an E. coli O55:H7 progenitor; however, a lack of genome sequence has hindered investigations on the divergence of human- and/or cattle-associated subtypes. Our goals were to 1) identify nucleotide polymorphisms for STEC O157 genetic subtype detection, 2) determine the phylogeny of STEC O157 genetic subtypes using polymorphism-derived genotypes and a phage insertion typing system, and 3) compare polymorphism-derived genotypes identified in this study with pulsed field gel electrophoresis (PFGE), the current gold standard for evaluating STEC O157 diversity. Using 762 nucleotide polymorphisms that were originally identified through whole-genome sequencing of 189 STEC O157 human- and cattle-isolated strains, we genotyped a collection of 426 STEC O157 strains. Concatenated polymorphism alleles defined 175 genotypes that were tagged by a minimal set of 138 polymorphisms. Eight major lineages of STEC O157 were identified, of which cattle are a reservoir for seven. Two lineages regularly harbored by cattle accounted for the majority of human disease in this study, whereas another was rarely represented in humans and may have evolved toward reduced human virulence. Notably, cattle are not a known reservoir for E. coli O55:H7 or STEC O157:H(-) (the first lineage to diverge within the STEC O157 serogroup), which both cause human disease. This result calls into question how cattle may have originally acquired STEC O157. The polymorphism-derived genotypes identified in this study did not surpass PFGE diversity assessed by BlnI and XbaI digestions in a subset of 93 strains. However, our results show that they are highly effective in assessing the evolutionary relatedness of epidemiologically unrelated STEC O157 genetic subtypes, including those associated with the cattle reservoir and human disease.