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.

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.

Escherichia coli O157:H7 tir 255 T > A allele strains differ in chromosomal and plasmid composition

Shiga toxin-producing Escherichia coli (STEC) O157:H7 strains with the T allele in the translocated intimin receptor polymorphism (tir) 255 A > T gene associate with human disease more than strains with an A allele; however, the allele is not thought to be the direct cause of this difference. We sequenced a diverse set of STEC O157:H7 strains (26% A allele, 74% T allele) to identify linked differences that might underlie disease association. The average chromosome and pO157 plasmid size and gene content were significantly greater within the tir 255 A allele strains. Eighteen coding sequences were unique to tir 255 A allele chromosomes, and three were unique to tir 255 T allele chromosomes. There also were non-pO157 plasmids that were unique to each tir 255 allele variant. The overall average number of prophages did not differ between tir 255 allele strains; however, there were different types between the strains. Genomic and mobile element variation linked to the tir 255 polymorphism may account for the increased frequency of the T allele isolates in human disease.

Genomic and Phenotypic Characterization of Shiga Toxin-Producing Escherichia albertii Strains Isolated from Wild Birds in a Major Agricultural Region in California

Escherichia albertii is an emerging foodborne pathogen. To better understand the pathogenesis and health risk of this pathogen, comparative genomics and phenotypic characterization were applied to assess the pathogenicity potential of E. albertii strains isolated from wild birds in a major agricultural region in California. Shiga toxin genes stx2f were present in all avian strains. Pangenome analyses of 20 complete genomes revealed a total of 11,249 genes, of which nearly 80% were accessory genes. Both core gene-based phylogenetic and accessory gene-based relatedness analyses consistently grouped the three stx2f-positive clinical strains with the five avian strains carrying ST7971. Among the three Stx2f-converting prophage integration sites identified, ssrA was the most common one. Besides the locus of enterocyte effacement and type three secretion system, the high pathogenicity island, OI-122, and type six secretion systems were identified. Substantial strain variation in virulence gene repertoire, Shiga toxin production, and cytotoxicity were revealed. Six avian strains exhibited significantly higher cytotoxicity than that of stx2f-positive E. coli, and three of them exhibited a comparable level of cytotoxicity with that of enterohemorrhagic E. coli outbreak strains, suggesting that wild birds could serve as a reservoir of E. albertii strains with great potential to cause severe diseases in humans.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08497-6.

Genome structural variation in Escherichia coli O157:H7

The human zoonotic pathogen Escherichia coli O157:H7 is defined by its extensive prophage repertoire including those that encode Shiga toxin, the factor responsible for inducing life-threatening pathology in humans. As well as introducing genes that can contribute to the virulence of a strain, prophage can enable the generation of large-chromosomal rearrangements (LCRs) by homologous recombination. This work examines the types and frequencies of LCRs across the major lineages of the O157:H7 serotype. We demonstrate that LCRs are a major source of genomic variation across all lineages of E. coli O157:H7 and by using both optical mapping and Oxford Nanopore long-read sequencing prove that LCRs are generated in laboratory cultures started from a single colony and that these variants can be recovered from colonized cattle. LCRs are biased towards the terminus region of the genome and are bounded by specific prophages that share large regions of sequence homology associated with the recombinational activity. RNA transcriptional profiling and phenotyping of specific structural variants indicated that important virulence phenotypes such as Shiga-toxin production, type-3 secretion and motility can be affected by LCRs. In summary, E. coli O157:H7 has acquired multiple prophage regions over time that act to continually produce structural variants of the genome. These findings raise important questions about the significance of this prophage-mediated genome contingency to enhance adaptability between environments.

Genomic analysis of shiga toxin-containing Escherichia coli O157:H7 isolated from Argentinean cattle

Cattle are the main reservoir of Enterohemorrhagic Escherichia coli (EHEC), with O157:H7 the distinctive serotype. EHEC is the main causative agent of a severe systemic disease, Hemolytic Uremic Syndrome (HUS). Argentina has the highest pediatric HUS incidence worldwide with 12–14 cases per 100,000 children. Herein, we assessed the genomes of EHEC O157:H7 isolates recovered from cattle in the humid Pampas of Argentina. According to phylogenetic studies, EHEC O157 can be divided into clades. Clade 8 strains that were classified as hypervirulent. Most of the strains of this clade have a Shiga toxin stx2a-stx2c genotype. To better understand the molecular bases related to virulence, pathogenicity and evolution of EHEC O157:H7, we performed a comparative genomic analysis of these isolates through whole genome sequencing. The isolates classified as clade 8 (four strains) and clade 6 (four strains) contained 13 to 16 lambdoid prophages per genome, and the observed variability of prophages was analysed. An inter strain comparison show that while some prophages are highly related and can be grouped into families, other are unique. Prophages encoding for stx2a were highly diverse, while those encoding for stx2c were conserved. A cluster of genes exclusively found in clade 8 contained 13 genes that mostly encoded for DNA binding proteins. In the studied strains, polymorphisms in Q antiterminator, the Q-stx2A intergenic region and the O and P γ alleles of prophage replication proteins are associated with different levels of Stx2a production. As expected, all strains had the pO157 plasmid that was highly conserved, although one strain displayed a transposon interruption in the protease EspP gene. This genomic analysis may contribute to the understanding of the genetic basis of the hypervirulence of EHEC O157:H7 strains circulating in Argentine cattle. This work aligns with other studies of O157 strain variation in other populations that shows key differences in Stx2a-encoding prophages.

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.

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.

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 TTP⁴: ≥5 CFU/g; TTP⁶: ≤5, ≥1 CFU/g; TTP⁸: ≤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 TTP⁴, with false negative rates of 0%.

Determination of gastrointestinal tract colonization sites from feedlot cattle transiently shedding or super-shedding Escherichia coli O157:H7 at harvest

AIMS

The objective of the study was to determine levels of Escherichia coli O157:H7 colonization in the gastrointestinal tract (GIT) of naturally shedding cattle shedding the pathogen at low- or super-shedder levels.

METHODS AND RESULTS

Over 2 years, feedlot cattle were sampled multiple times for faecal shedding of E. coli O157:H7. Just prior to harvest (1-2 days), animals that were super-shedders (≥10⁴  CFU per gram of faeces) were specifically identified, and based on the longer term screening data, pen cohorts that were low-shedders (years 1 and 2) or chronic-shedders (year 1) were also identified. At harvest, samples were collected from throughout the GIT, including the rectoanal junction (RAJ) for enumeration and enrichment of E. coli O157:H7. The mouth samples exhibited the greatest prevalence for the pathogen, and the abomasum and rumen exhibited the lowest prevalence (P < 0·05). Super-shedders had significantly greater prevalence for all GIT locations except the mouth and abomasum compared to low-shedders, but the super-shedders were the only animals with positive abomasum samples. Samples from the super-shedders were enumerable for most GIT locations, and the rectum and RAJ locations were the only locations that were significantly greater than other locations (P < 0·05).

CONCLUSIONS

Across all animals naturally exposed to E. coli O157:H7, the risk of ingestion is high, but rumen and abomasum are potential barriers to passage. In super-shedders, the passage through the GIT was greater, allowing colonization in the rectum and at the RAJ.

SIGNIFICANCE AND IMPACT OF THE STUDY

Escherichia coli O157:H7 low-shedding cattle had lower pathogen levels throughout the GIT, indicating intrinsic GIT factors to these cattle may reduce pathogen passage through the GIT, including the abomasum, and minimize risk of RAJ colonization.

Pathogenomes of Atypical Non-shigatoxigenic Escherichia coli NSF/SF O157:H7/NM: Comprehensive Phylogenomic Analysis Using Closed Genomes

The toxigenic conversion of Escherichia coli strains by Shiga toxin-converting (Stx) bacteriophages were prominent and recurring events in the stepwise evolution of enterohemorrhagic E. coli (EHEC) O157:H7 from an enteropathogenic (EPEC) O55:H7 ancestor. Atypical, attenuated isolates have been described for both non-sorbitol fermenting (NSF) O157:H7 and SF O157:NM serotypes, which are distinguished by the absence of Stx, the characteristic virulence hallmark of Stx-producing E. coli (STEC). Such atypical isolates either never acquired Stx-phages or may have secondarily lost stx during the course of infection, isolation, or routine subculture; the latter are commonly referred to as LST (Lost Shiga Toxin)-isolates. In this study we analyzed the genomes of 15 NSF O157:H7 and SF O157:NM strains from North America, Europe, and Asia that are characterized by the absence of stx, the virulence hallmark of STEC. The individual genomic basis of the Stx (-) phenotype has remained largely undetermined as the majority of STEC genomes in public genome repositories were generated using short read technology and are in draft stage, posing a major obstacle for the high-resolution whole genome sequence typing (WGST). The application of LRT (long-read technology) sequencing provided us with closed genomes, which proved critical to put the atypical non-shigatoxigenic NSF O157:H7 and SF O157:NM strains into the phylogenomic context of the stepwise evolutionary model. Availability of closed chromosomes for representative Stx (-) NSF O157:H7 and SF O157:NM strains allowed to describe the genomic basis and individual evolutionary trajectories underlying the absence of Stx at high accuracy and resolution. The ability of LRT to recover and accurately assemble plasmids revealed a strong correlation between the strains' featured plasmid genotype and chromosomally inferred clade, which suggests the coevolution of the chromosome and accessory plasmids. The identified ancestral traits in the pSFO157 plasmid of NSF O157:H7 strain LSU-61 provided additional evidence for its intermediate status. Taken together, these observations highlight the utility of LRTs for advancing our understanding of EHEC O157:H7/NM pathogenome evolution. Insights into the genomic and phenotypic plasticity of STEC on a lineage- and genome-wide scale are foundational to improve and inform risk assessment, biosurveillance, and prevention strategies.

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.

Complete Genome Sequencing of Four Arcobacter Species Reveals a Diverse Suite of Mobile Elements

Arcobacter species are recovered from a wide variety of sources, including animals, food, and both fresh and marine waters. Several Arcobacter species have also been recovered from human clinical samples and are thus associated tentatively with food- and water-borne human illnesses. Genome sequencing of the poultry isolate Arcobacter cibarius H743 and the Arcobacter acticola, Arcobacter pacificus, and Arcobacter porcinus type strains identified a large number and variety of insertion sequences. This study presents an analysis of these A. acticola, A. cibarius, A. pacificus, and A. porcinus IS elements. The four genomes sequenced here contain 276 complete and degenerate IS elements, representing 13 of the current 29 prokaryotic IS element families. Expansion of the analysis to include 15 other previously sequenced Arcobacter spp. added 73 complete and degenerate IS elements. Several of these IS elements were identified in two or more Arcobacter species, suggesting movement by horizontal gene transfer between the arcobacters. These IS elements are putatively associated with intragenomic deletions and inversions, and tentative movement of antimicrobial resistance genes. The A. cibarius strain H743 megaplasmid contains multiple IS elements common to the chromosome and, unusually, a complete ribosomal RNA locus, indicating that larger scale genomic rearrangements, potentially resulting from IS element-mediated megaplasmid cointegration and resolution may be occurring within A. cibarius and possibly other arcobacters. The presence of such a large and varied suite of mobile elements could have profound effects on Arcobacter biology and evolution.

A Computational Method to Quantify the Effects of Slipped Strand Mispairing on Bacterial Tetranucleotide Repeats

The virulence and pathogenicity of bacterial pathogens are related to their adaptability to changing environments. One process enabling adaptation is based on minor changes in genome sequence, as small as a few base pairs, within segments of genome called simple sequence repeats (SSRs) that consist of multiple copies of a short sequence (from one to several nucleotides), repeated in series. SSRs are found in eukaryotes as well as prokaryotes, and length variation in them occurs at frequencies up to a million-fold higher than bacterial point mutations through the process of slipped strand mispairing (SSM) by DNA polymerase during replication. The characterization of SSR length by standard sequencing methods is complicated by the appearance of length variation introduced during the sequencing process that obscures the lower abundance repeat number variants in a population. Here we report a computational approach to correct for sequencing process-induced artifacts, validated for tetranucleotide repeats by use of synthetic constructs of fixed, known length. We apply this method to a laboratory culture of Histophilus somni, prepared from a single colony, and demonstrate that the culture consists of populations of distinct sequence phase and length variants at individual tetranucleotide SSR loci.

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 (>10³ 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.

Enterohaemorrhagic E. coli (EHEC) O157 strains are found in cattle where they are asymptomatic, while human exposure can lead to severe symptoms including bloody diarrhoea and kidney damage due to the activity of Shiga toxin (Stx). The most serious symptoms in humans are associated with isolates that encode Stx subtype 2a. The advantage of these toxins in the animal reservoir is still not clear, however there is experimental evidence implicating Stx with increased bacterial adherence, immune modulation and suppression of predatory protozoa. In this study, the hypothesis that Stx2a is important for super-shedding and calf-to-calf transmission was tested by comparing excretion and transmission dynamics of E. coli O157 strains with and without Stx2a. While Stx2a did not alter excretion levels when calfs were orally challenge, it enabled colonisation of more in contact ‘sentinel’ animals in our transmission model. We show that Stx2a is generally induced more rapidly than Stx2c, resulting in increased levels of Stx2a expression. Both Stx2a and Stx2c were able to restrict cellular proliferation of epithelial cells in cultured bovine enteroids. Taken together, we propose that rapid production of Stx2a and its role in establishing E. coli O157 colonisation in the bovine gastrointestinal tract facilitate effective transmission and have led to its expansion in the cattle E. coli O157 population.

Occurrence of Escherichia coli O157:H7 in Pest Flies Captured in Leafy Greens Plots Grown Near a Beef Cattle Feedlot

Leafy greens are leading vehicles for Escherichia coli O157:H7 foodborne illness. Pest flies can harbor this pathogen and may disseminate it to produce. We determined the occurrence of E. coli O157:H7-positive flies in leafy greens planted up to 180 m from a cattle feedlot and assessed their relative risk to transmit this pathogen to leafy greens. The primary fly groups captured on sticky traps at the feedlot and leafy greens plots included house flies (Musca domestica L.), face flies (Musca autumnalis L.), stable flies (Stomoxys calcitrans L.), flesh flies (family Sarcophagidae), and blow flies (family Calliphoridae). E. coli O157:H7 carriage rates of house, face, flesh, and blow flies were similar (P > 0.05), ranging from 22.3 to 29.0 flies per 1,000 flies. In contrast, the carriage rate of stable flies was lower at 1.1 flies per 1,000 flies (P < 0.05). Differences in carriage rates are likely due to the uses of fresh bovine feces and manure by these different pest fly groups. E. coli O157:H7 carriage rates of total flies did not differ (P > 0.05) by distance (ranging from 0 to 180 m) from the feedlot. Most fly isolates were the same predominant pulsed-field gel electrophoresis types found in feedlot surface manure and leafy greens, suggesting a possible role for flies in transmitting E. coli O157:H7 to the leafy greens. However, further research is needed to clarify this role and to determine set-back distances between cattle production facilities and produce crops that will reduce the risk for pathogen contamination by challenging mechanisms like flies.

Complete Genome Sequences of the Arcobacter cryaerophilus Strains ATCC 43158T and ATCC 49615

Arcobacter cryaerophilus was originally recovered from aborted bovine and porcine fetuses, but it has been subsequently isolated from meat, water, and human clinical samples. This study describes the complete whole-genome sequences of two A. cryaerophilus strains, ATCC 43158T (=A 169/BT =LMG 24291T) and ATCC 49615 (=CDC D2610 =LMG 10829).

Complete Genome Sequence of the Arcobacter suis Type Strain LMG 26152

Arcobacter species are prevalent in pigs, and strains have been isolated from pig feces and pork meat; some Arcobacter strains may be porcine abortifacients. Arcobacter suis was recovered from pork meat in Spain. This study describes the whole-genome sequence of the A. suis type strain LMG 26152 (=F41T =CECT 7833T).

Complete Genome Sequence of the Arcobacter molluscorum Type Strain LMG 25693

As components of freshwater and marine microflora, Arcobacter spp. are often recovered from shellfish, such as mussels, clams, and oysters. Arcobacter molluscorum was isolated from mussels from the Ebro Delta in Catalonia, Spain. This article describes the whole-genome sequence of the A. molluscorum strain LMG 25693T (= F98-3T = CECT 7696T).

Complete Genome Sequence of the Arcobacter ellisii Type Strain LMG 26155

Arcobacter spp. are highly prevalent in contaminated environmental waters and have been recovered from both freshwater and seawater, with several species isolated from shellfish. Arcobacter ellisii was recovered from mussels collected in Catalonia, Spain. This study describes the whole-genome sequence of the A. ellisii type strain LMG 26155 (=F79-6T =CECT 7837T).

Draft Genome Sequences of Seven Strains of Shiga Toxin-Producing Escherichia coli O111 with Variation in Their Sensitivity to Novobiocin

Inclusion of novobiocin as a selective agent for enrichment media and selective agars inhibits the growth of some Shiga toxin-producing Escherichia coli (STEC) strains, particularly non-O157 STEC, which can yield false-negative detection results. Here, we report the draft genomic sequences of seven STEC O111 isolates with different sensitivities to novobiocin.

Closed Genome and Comparative Phylogenetic Analysis of the Clinical Multidrug Resistant Shigella sonnei Strain 866

Shigella sonnei is responsible for the majority of shigellosis infections in the US with over 500,000 cases reported annually. Here, we present the complete genome of the clinical multidrug resistant (MDR) strain 866, which is highly susceptible to bacteriophage infections. The strain has a circular chromosome of 4.85 Mb and carries a 113 kb MDR plasmid. This IncB/O/K/Z-type plasmid, termed p866, confers resistance to five different classes of antibiotics including ß-lactamase, sulfonamide, tetracycline, aminoglycoside, and trimethoprim. Comparative analysis of the plasmid architecture and gene inventory revealed that p866 shares its plasmid backbone with previously described IncB/O/K/Z-type Shigella spp. and Escherichia coli plasmids, but is differentiated by the insertion of antibiotic resistance cassettes, which we found associated with mobile genetic elements such as Tn3, Tn7, and Tn10. A whole genome-derived phylogenetic reconstruction showed the evolutionary relationships of S. sonnei strain 866 and the four established Shigella species, highlighting the clonal nature of S. sonnei.

Comparative genomics of Salmonella enterica serovar Montevideo reveals lineage-specific gene differences that may influence ecological niche association

Salmonella enterica serovar Montevideo has been linked to recent foodborne illness outbreaks resulting from contamination of products such as fruits, vegetables, seeds and spices. Studies have shown that Montevideo also is frequently associated with healthy cattle and can be isolated from ground beef, yet human salmonellosis outbreaks of Montevideo associated with ground beef contamination are rare. This disparity fuelled our interest in characterizing the genomic differences between Montevideo strains isolated from healthy cattle and beef products, and those isolated from human patients and outbreak sources. To that end, we sequenced 13 Montevideo strains to completion, producing high-quality genome assemblies of isolates from human patients (n=8) or from healthy cattle at slaughter (n=5). Comparative analysis of sequence data from this study and publicly available sequences (n=72) shows that Montevideo falls into four previously established clades, differentially occupied by cattle and human strains. The results of these analyses reveal differences in metabolic islands, environmental adhesion determinants and virulence factors within each clade, and suggest explanations for the infrequent association between bovine isolates and human illnesses.

Comparative Genomics of All Three Campylobacter sputorum Biovars and a Novel Cattle-Associated C. sputorum Clade

Campylobacter sputorum is a nonthermotolerant campylobacter that is primarily isolated from food animals such as cattle and sheep. C. sputorum is also infrequently associated with human illness. Based on catalase and urease activity, three biovars are currently recognized within C. sputorum: bv. sputorum (catalase negative, urease negative), bv. fecalis (catalase positive, urease negative), and bv. paraureolyticus (catalase negative, urease positive). A multi-locus sequence typing (MLST) method was recently constructed for C. sputorum. MLST typing of several cattle-associated C. sputorum isolates suggested that they are members of a divergent C. sputorum clade. Although catalase positive, and thus technically bv. fecalis, the taxonomic position of these strains could not be determined solely by MLST. To further characterize C. sputorum, the genomes of four strains, representing all three biovars and the divergent clade, were sequenced to completion. Here we present a comparative genomic analysis of the four C. sputorum genomes. This analysis indicates that the three biovars and the cattle-associated strains are highly related at the genome level with similarities in gene content. Furthermore, the four genomes are strongly syntenic with one or two minor inversions. However, substantial differences in gene content were observed among the three biovars. Finally, although the strain representing the cattle-associated isolates was shown to be C. sputorum, it is possible that this strain is a member of a novel C. sputorum subspecies; thus, these cattle-associated strains may form a second taxon within C. sputorum.

Comparative Genomic Analysis Identifies a Campylobacter Clade Deficient in Selenium Metabolism

The nonthermotolerant Campylobacter species C. fetus, C. hyointestinalis, C. iguaniorum, and C. lanienae form a distinct phylogenetic cluster within the genus. These species are primarily isolated from foraging (swine) or grazing (e.g., cattle, sheep) animals and cause sporadic and infrequent human illness. Previous typing studies identified three putative novel C. lanienae-related taxa, based on either MLST or atpA sequence data. To further characterize these putative novel taxa and the C. fetus group as a whole, 76 genomes were sequenced, either to completion or to draft level. These genomes represent 26 C. lanienae strains and 50 strains of the three novel taxa. C. fetus, C. hyointestinalis and C. iguaniorum genomes were previously sequenced to completion; therefore, a comparative genomic analysis across the entire C. fetus group was conducted (including average nucleotide identity analysis) that supports the initial identification of these three novel Campylobacter species. Furthermore, C. lanienae and the three putative novel species form a discrete clade within the C. fetus group, which we have termed the C. lanienae clade. This clade is distinguished from other members of the C. fetus group by a reduced genome size and distinct CRISPR/Cas systems. Moreover, there are two signature characteristics of the C. lanienae clade. C. lanienae clade genomes carry four to ten unlinked and similar, but nonidentical, flagellin genes. Additionally, all 76 C. lanienae clade genomes sequenced demonstrate a complete absence of genes related to selenium metabolism, including genes encoding the selenocysteine insertion machinery, selenoproteins, and the selenocysteinyl tRNA.

The impact of the bovine faecal microbiome on Escherichia coli O157:H7 prevalence and enumeration in naturally infected cattle

AIMS

The objective of this study was to determine if the faecal microbiome has an association with Escherichia coli O157:H7 prevalence and enumeration.

METHODS AND RESULTS

Pyrosequencing analysis of faecal microbiome was performed from feedlot cattle fed one of three diets: (i) 94 heifers fed low concentrate (LC) diet, (ii) 142 steers fed moderate concentrate (MC) diet, and (iii) 132 steers fed high concentrate (HC) diet. A total of 322 585 OTUs were calculated from 2,411,122 high-quality sequences obtained from 368 faecal samples. In the LC diet group, OTUs assigned to the orders Clostridiales and RF39 (placed within the class Mollicutes) were positively correlated with both E. coli O157:H7 prevalence and enumeration. In the MC diet group, OTUs assigned to Prevotella copri were positively correlated with both E. coli O157:H7 prevalence and enumeration, whereas OTUs assigned to Prevotella stercorea were negatively correlated with both E. coli O157:H7 prevalence and enumeration. In both the MC diet group and the HC diet group, OTUs assigned to taxa placed within Clostridiales were both positively and negatively correlated with both E. coli O157:H7 prevalence and enumeration. However, all correlations were weak. In both the MC diet group and the HC diet group, stepwise linear regression through backward elimination analyses indicated that these OTUs were significantly correlated (P < 0·001) with prevalence or enumeration, explaining as much as 50% of variability in E. coli O157:H7 prevalence or enumeration.

CONCLUSIONS

Individual colonic bacterial species have little impact on E. coli O157:H7 shedding but collectively groups of bacteria were strongly associated with pathogen shedding.

SIGNIFICANCE AND IMPACT OF THE STUDY

Bacterial groups in the bovine colon may impact faecal shedding of the zoonotic pathogen E. coli O157:H7, and manipulation of the intestinal microbiota to alter these bacteria may reduce shedding of this pathogen and foodborne illnesses.

The effects of monensin in diets fed to finishing beef steers and heifers on growth performance and fecal shedding of Escherichia coli O157:H7

Two experiments were conducted to study the effects of monensin dose on growth performance and O157:H7 shedding in finishing beef cattle. In Exp. 1, 198 heifers (298 ± 1.1 kg BW) were allocated to 1 of 2 treatments consisting of 1) 200 mg/heifer daily of monensin and 2) 400 mg/heifer daily of monensin and fed for 151 d. In Exp. 2, 199 steers (430 ± 1.9 kg BW) were stratified by BW and allocated to 1 of 2 treatments consisting of 1) 0 mg/steer daily of monensin and 2) 400 mg/steer daily of monensin and fed for 128 d. For both experiments, there were 4 pen replicates per treatment. For Exp. 1 and Exp. 2, the model included the fixed effect of treatment for growth performance measures and the fixed effects of treatment, time, and treatment × time interaction, respectively, for O157:H7 shedding. In Exp. 1, final BW was 1.9% greater for heifers fed 400 mg/d monensin than for heifers fed 200 mg/d monensin ( = 0.05). Furthermore, ADG was 4.9% greater ( = 0.05) and G:F was 3.1% greater ( = 0.04) when the heifers were fed 400 mg/d monensin vs. 200 mg/d monensin. Pen prevalence for O157:H7 ( = 0.96) and the percentage of animals in the pen shedding O157:H7 at enumerable levels ( = 0.82) did not differ between heifers fed 200 mg/d monensin and heifers fed 400 mg/d monensin over the 4 sampling periods. For Exp. 2, steers fed the supplement containing monensin had a 1.9% greater final BW ( = 0.04) and a 5.2% greater ADG ( = 0.02) than steers fed a control supplement without monensin. No differences in DMI or G:F were noted across the treatments ( ≥ 0.14). O157:H7 percentage of enumerable cattle within the pen was greater for the steers fed monensin than the control steers not fed monensin than the control steers not fed monensin ( = 0.02) over the 4 sampling periods. However, the percentage of animals in the pen shedding O157:H7 (prevalence positive) did not differ between treatments ( = 0.18), nor did the average fecal counts ( = 0.45). In conclusion, feeding a higher dose (400 mg/d) of monensin improved final BW and ADG compared with a low dose of monensin or a no-monensin control in steers and heifers across multiple years. The percentage of animals shedding O157:H7 at enumerable levels was greater for steers fed the monensin supplement than for steers fed the control supplement, yet the presence of monensin, irrespective of the dose, did not affect the percentage of animals in the pen shedding O157:H7.

Evolution of a zoonotic pathogen: investigating prophage diversity in enterohaemorrhagic Escherichia coli O157 by long-read sequencing

Enterohaemorrhagic Escherichia coli (EHEC) O157 is a zoonotic pathogen for which colonization of cattle and virulence in humans is associated with multiple horizontally acquired genes, the majority present in active or cryptic prophages. Our understanding of the evolution and phylogeny of EHEC O157 continues to develop primarily based on core genome analyses; however, such short-read sequences have limited value for the analysis of prophage content and its chromosomal location. In this study, we applied Single Molecule Real Time (SMRT) sequencing, using the Pacific Biosciences long-read sequencing platform, to isolates selected from the main sub-clusters of this clonal group. Prophage regions were extracted from these sequences and from published reference strains. Genome position and prophage diversity were analysed along with genetic content. Prophages could be assigned to clusters, with smaller prophages generally exhibiting less diversity and preferential loss of structural genes. Prophages encoding Shiga toxin (Stx) 2a and Stx1a were the most diverse, and more variable compared to prophages encoding Stx2c, further supporting the hypothesis that Stx2c-prophage integration was ancestral to acquisition of other Stx types. The concept that phage type (PT) 21/28 (Stx2a+, Stx2c+) strains evolved from PT32 (Stx2c+) was supported by analysis of strains with excised Stx-encoding prophages. Insertion sequence elements were over-represented in prophage sequences compared to the rest of the genome, showing integration in key genes such as stx and an excisionase, the latter potentially acting to capture the bacteriophage into the genome. Prophage profiling should allow more accurate prediction of the pathogenic potential of isolates.

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.

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

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

Disinfectant and Antimicrobial Susceptibility Profiles of the Big Six Non-O157 Shiga Toxin-Producing Escherichia coli Strains from Food Animals and Humans

The disinfectant and antimicrobial susceptibility profiles of 138 non-O157 Shiga toxin-producing Escherichia coli strains (STECs) from food animals and humans were determined. Antimicrobial resistance (AMR) was moderate (39.1% of strains) in response to 15 antimicrobial agents. Animal strains had a lower AMR prevalence (35.6%) than did human strains (43.9%) but a higher prevalence of the resistance profile GEN-KAN-TET. A decreasing prevalence of AMR was found among animal strains from serogroups O45 > O145 > O121 > O111 > O26 > O103 and among human strains from serogroups O145 > O103 > O26 > O111 > O121 > O45. One animal strain from serogroups O121 and O145 and one human strain from serogroup O26 had extensive drug resistance. A high prevalence of AMR in animal O45 and O121 strains and no resistance or a low prevalence of resistance in human strains from these serogroups suggests a source other than food animals for human exposure to these strains. Among the 24 disinfectants evaluated, all strains were susceptible to triclosan. Animal strains had a higher prevalence of resistance to chlorhexidine than did human strains. Both animal and human strains had a similar low prevalence of low-level benzalkonium chloride resistance, and animal and human strains had similar susceptibility profiles for most other disinfectants. Benzyldimethylammonium chlorides and C10AC were the primary active components in disinfectants DC&R and P-128, respectively, against non-O157 STECs. A disinfectant FS512 MIC ≥ 8 μg/ml was more prevalent among animal O121 strains (61.5%) than among human O121 strains (25%), which may also suggest a source of human exposure to STEC O121 other than food animals. Bacterial inhibition was not dependent solely on pH but was correlated with the presence of dissociated organic acid species and some undissociated acids.

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.

Sequence of Colonization Determines the Composition of Mixed Biofilms by Escherichia coli O157:H7 and O111:H8 Strains

Bacterial biofilms are one of the potential sources of cross-contamination in food processing environments. Shiga toxin-producing Escherichia coli (STEC) O157:H7 and O111:H8 are important foodborne pathogens capable of forming biofilms, and the coexistence of these two STEC serotypes has been detected in various food samples and in multiple commercial meat plants throughout the United States. Here, we investigated how the coexistence of these two STEC serotypes and their sequence of colonization could affect bacterial growth competition and mixed biofilm development. Our data showed that E. coli O157:H7 strains were able to maintain a higher cell percentage in mixed biofilms with the co-inoculated O111:H8 companion strains, even though the results of planktonic growth competition were strain dependent. On solid surfaces with preexisting biofilms, the sequence of colonization played a critical role in determining the composition of the mixed biofilms because early stage precolonization significantly affected the competition results between the E. coli O157:H7 and O111:H8 strains. The precolonizer of either serotype was able to outgrow the other serotype in both planktonic and biofilm phases. The competitive interactions among the various STEC serotypes would determine the composition and structure of the mixed biofilms as well as their potential risks to food safety and public health, which is largely influenced by the dominant strains in the mixtures. Thus, the analysis of mixed biofilms under various conditions would be of importance to determine the nature of mixed biofilms composed of multiple microorganisms and to help implement the most effective disinfection operations accordingly.

Canonical Single Nucleotide Polymorphisms (SNPs) for High-Resolution Subtyping of Shiga-Toxin Producing Escherichia coli (STEC) O157:H7

The objective of this study was to develop a canonical, parsimoniously-informative SNP panel for subtyping Shiga-toxin producing Escherichia coli (STEC) O157:H7 that would be consistent with epidemiological, PFGE, and MLVA clustering of human specimens. Our group had previously identified 906 putative discriminatory SNPs, which were pared down to 391 SNPs based on their prevalence in a test set. The 391 SNPs were screened using a high-throughput form of TaqMan PCR against a set of clinical isolates that represent the most diverse collection of O157:H7 isolates from outbreaks and sporadic cases examined to date. Another 30 SNPs identified by others were also screened using the same method. Two additional targets were tested using standard TaqMan PCR endpoint analysis. These 423 SNPs were reduced to a 32 SNP panel with the almost the same discriminatory value. While the panel partitioned our diverse set of isolates in a manner that was consistent with epidemiological data and PFGE and MLVA phylogenies, it resulted in fewer subtypes than either existing method and insufficient epidemiological resolution in 10 of 47 clusters. Therefore, another round of SNP discovery was undertaken using comparative genomic resequencing of pooled DNA from the 10 clusters with insufficient resolution. This process identified 4,040 potential SNPs and suggested one of the ten clusters was incorrectly grouped. After its removal, there were 2,878 SNPs, of which only 63 were previously identified and 438 occurred across multiple clusters. Among highly clonal bacteria like STEC O157:H7, linkage disequilibrium greatly limits the number of parsimoniously informative SNPs. Therefore, it is perhaps unsurprising that our panel accounted for the potential discriminatory value of numerous other SNPs reported in the literature. We concluded published O157:H7 SNPs are insufficient for effective epidemiological subtyping. However, the 438 multi-cluster SNPs we identified may provide the additional information required.

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.

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

Shiga toxin-producing Escherichia coli (STEC) O26 is the second leading E. coli serogroup responsible for human illness outbreaks behind E. coli O157:H7. Recent outbreaks have been linked to emerging pathogenic O26:H11 strains harboring stx 2 only. Cattle have been recognized as an important reservoir of O26 strains harboring stx 1; however the reservoir of these emerging stx 2 strains is unknown. The objective of this study was to identify nucleotide polymorphisms in human and cattle-derived strains in order to compare differences in polymorphism derived genotypes and virulence gene profiles between the two host species. Whole genome sequencing was performed on 182 epidemiologically unrelated O26 strains, including 109 human-derived strains and 73 non-human-derived strains. A panel of 289 O26 strains (241 STEC and 48 non-STEC) was subsequently genotyped using a set of 283 polymorphisms identified by whole genome sequencing, resulting in 64 unique genotypes. Phylogenetic analyses identified seven clusters within the O26 strains. The seven clusters did not distinguish between isolates originating from humans or cattle; however, clusters did correspond with particular virulence gene profiles. Human and non-human-derived strains harboring stx 1 clustered separately from strains harboring stx 2, strains harboring eae, and non-STEC strains. Strains harboring stx 2 were more closely related to non-STEC strains and strains harboring eae than to strains harboring stx 1. The finding of human and cattle-derived strains with the same polymorphism derived genotypes and similar virulence gene profiles, provides evidence that similar strains are found in cattle and humans and transmission between the two species may occur.

Reducing assembly complexity of microbial genomes with single-molecule sequencing

BACKGROUND

The short reads output by first- and second-generation DNA sequencing instruments cannot completely reconstruct microbial chromosomes. Therefore, most genomes have been left unfinished due to the significant resources required to manually close gaps in draft assemblies. Third-generation, single-molecule sequencing addresses this problem by greatly increasing sequencing read length, which simplifies the assembly problem.

RESULTS

To measure the benefit of single-molecule sequencing on microbial genome assembly, we sequenced and assembled the genomes of six bacteria and analyzed the repeat complexity of 2,267 complete bacteria and archaea. Our results indicate that the majority of known bacterial and archaeal genomes can be assembled without gaps, at finished-grade quality, using a single PacBio RS sequencing library. These single-library assemblies are also more accurate than typical short-read assemblies and hybrid assemblies of short and long reads.

CONCLUSIONS

Automated assembly of long, single-molecule sequencing data reduces the cost of microbial finishing to $1,000 for most genomes, and future advances in this technology are expected to drive the cost lower. This is expected to increase the number of completed genomes, improve the quality of microbial genome databases, and enable high-fidelity, population-scale studies of pan-genomes and chromosomal organization.

Comparative genomics of the Campylobacter lari group

The Campylobacter lari group is a phylogenetic clade within the epsilon subdivision of the Proteobacteria and is part of the thermotolerant Campylobacter spp., a division within the genus that includes the human pathogen Campylobacter jejuni. The C. lari group is currently composed of five species (C. lari, Campylobacter insulaenigrae, Campylobacter volucris, Campylobacter subantarcticus, and Campylobacter peloridis), as well as a group of strains termed the urease-positive thermophilic Campylobacter (UPTC) and other C. lari-like strains. Here we present the complete genome sequences of 11 C. lari group strains, including the five C. lari group species, four UPTC strains, and a lari-like strain isolated in this study. The genome of C. lari subsp. lari strain RM2100 was described previously. Analysis of the C. lari group genomes indicates that this group is highly related at the genome level. Furthermore, these genomes are strongly syntenic with minor rearrangements occurring only in 4 of the 12 genomes studied. The C. lari group can be bifurcated, based on the flagella and flagellar modification genes. Genomic analysis of the UPTC strains indicated that these organisms are variable but highly similar, closely related to but distinct from C. lari. Additionally, the C. lari group contains multiple genes encoding hemagglutination domain proteins, which are either contingency genes or linked to conserved contingency genes. Many of the features identified in strain RM2100, such as major deficiencies in amino acid biosynthesis and energy metabolism, are conserved across all 12 genomes, suggesting that these common features may play a role in the association of the C. lari group with coastal environments and watersheds.

Investigation of bacterial diversity in the feces of cattle fed different diets

The objective of this study is to investigate individual animal variation of bovine fecal microbiota including as affected by diets. Fecal samples were collected from 426 cattle fed 1 of 3 diets typically fed to feedlot cattle: 1) 143 steers fed finishing diet (83% dry-rolled corn, 13% corn silage, and 4% supplement), 2) 147 steers fed late growing diet (66% dry-rolled corn, 26% corn silage, and 8% supplement), and 3) 136 heifers fed early growing diet (70% corn silage and 30% alfalfa haylage). Bacterial 16S rRNA gene amplicons were determined from individual fecal samples using next-generation pyrosequencing technology. A total of 2,149,008 16S rRNA gene sequences from 333 cattle with at least 2,000 sequences were analyzed. Firmicutes and Bacteroidetes were dominant phyla in all fecal samples. At the genus level, Oscillibacter, Turicibacter, Roseburia, Fecalibacterium, Coprococcus, Clostridium, Prevotella, and Succinivibrio were represented by more than 1% of total sequences. However, numerous sequences could not be assigned to a known genus. Dominant unclassified groups were unclassified Ruminococcaceae and unclassified Lachnospiraceae that could be classified to a family but not to a genus. These dominant genera and unclassified groups differed (P < 0.001) with diets. A total of 176,692 operational taxonomic units (OTU) were identified in combination across all the 333 cattle. Only 2,359 OTU were shared across 3 diet groups. UniFrac analysis showed that bacterial communities in cattle feces were greatly affected by dietary differences. This study indicates that the community structure of fecal microbiota in cattle is greatly affected by diet, particularly between forage- and concentrate-based diets.

Disinfectant and antibiotic susceptibility profiles of Escherichia coli O157:H7 strains from cattle carcasses, feces, and hides and ground beef from the United States

The disinfectant and antibiotic susceptibility profiles of 344 Escherichia coli O157:H7 strains from cattle carcasses, feces, and hides and ground beef from the United States were determined. A low prevalence of antibiotic resistance was observed (14%). The highest prevalences of resistance were to sulfisoxazole (10.5%), tetracycline (9.9%), streptomycin (7%), and chloramphenicol (4.9%). Four strains were resistant to eight antibiotics (two strains from ground beef and one strain each from hide and preevisceration carcass swabs of cull cattle at harvest). Pulsed-field gel electrophoresis analysis of the E. coli O157:H7 strains revealed two major groups (designated 1 and 2) composed of 17 and 20 clusters, respectively. Clusters 1A, 1B, 1C, and 1G.1 were associated with multidrug-resistant strains. There was no observed correlation between disinfectant resistance and antibiotic resistance. Sixty-nine (20%) of the 344 strains were resistant to chlorhexidine or benzalkonium chloride or the MICs of benzyldimethyldodecylammonium chloride were elevated. Inducible resistance was observed at elevated concentrations of antibiotics (1.4%) and disinfectants (6.1%). The highest rate of disinfectant inducible resistance was to OdoBan, quaternary ammonium chlorides, and the surface disinfectants F25, FS512, and MG, which are used in dairies, restaurants, and food processing plants. High MICs (1,024 to 4,096 m g/ml) of acetic, lactic, and citric acids were found. The decreasing order of acid potency based on molar MICs (MICs(molar)) was acetic, citric, and lactic acid. The correlation of the concentration of dissociated organic acids and MICs(molar) strongly suggests that the observed inhibition of E. coli O157:H7 was primarily due to dissociated forms of the acids.

Biofilm formation by Shiga toxin-producing Escherichia coli O157:H7 and Non-O157 strains and their tolerance to sanitizers commonly used in the food processing environment

Shiga toxin-producing Escherichia coli (STEC) strains are important foodborne pathogens. Among these, E. coli O157:H7 is the most frequently isolated STEC serotype responsible for foodborne diseases. However, the non-O157 serotypes have been associated with serious outbreaks and sporadic diseases as well. It has been shown that various STEC serotypes are capable of forming biofilms on different food or food contact surfaces that, when detached, may lead to cross-contamination. Bacterial cells at biofilm stage also are more tolerant to sanitizers compared with their planktonic counterparts, which makes STEC biofilms a serious food safety concern. In the present study, we evaluated the potency of biofilm formation by a variety of STEC strains from serotypes O157:H7, O26:H11, and O111:H8; we also compared biofilm tolerance with two types of common sanitizers, a quaternary ammonium chloride-based sanitizer and chlorine. Our results demonstrated that biofilm formation by various STEC serotypes on a polystyrene surface was highly strain-dependent, whereas the two non-O157 serotypes showed a higher potency of pellicle formation at air-liquid interfaces on a glass surface compared with serotype O157:H7. Significant reductions of viable biofilm cells were achieved with sanitizer treatments. STEC biofilm tolerance to sanitization was strain-dependent regardless of the serotypes. Curli expression appeared to play a critical role in STEC biofilm formation and tolerance to sanitizers. Our data indicated that multiple factors, including bacterial serotype and strain, surface materials, and other environmental conditions, could significantly affect STEC biofilm formation. The high potential for biofilm formation by various STEC serotypes, especially the strong potency of pellicle formation by the curli-positive non-O157 strains with high sanitization tolerance, might contribute to bacterial colonization on food contact surfaces, which may result in downstream product contamination.

Comparison of bacterial communities in faeces of beef cattle fed diets containing corn and wet distillers' grain with solubles

AIM

The mammalian intestinal microflora has been shown to impact host physiology. In cattle, intestinal bacteria are also associated with faecal contamination of environmental sources and human illness via foodborne pathogens. Use of wet distillers' grains with solubles (WDGS) in cattle feed creates a gastrointestinal environment where some bacterial species are enriched. Here, we examine if a diet containing 40% WDGS results in fundamentally different microbial community structures.

METHODS AND RESULTS

The 20,002 16S r-RNA gene sequences from 20 beef cattle were analysed using Sanger sequencing methods. At the genus level, Prevotella (Gram negative) and Anaerobacter (Gram positive) were the most frequently occurring bacteria in our beef cattle faecal samples. Diet-associated differences in prevalence were noted for Prevotella but not Anaerobacter.

CONCLUSIONS

Diet affects community structure. Faecal communities of co-housed beef cattle are not identical.

SIGNIFICANCE AND IMPACT OF THE STUDY

It is known that a diet of 40% corn-based WDGS increases the generic Escherichia coli in the faeces and enriches E. coli O157:H7. The results from the current study suggest that in addition to previously observed changes in E. coli, the entire bacterial community structure is different for animals fed 40% corn-based WDGS compared to a traditional corn-finishing diet.

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.

Bacterial community analysis of beef cattle feedlots reveals that pen surface is distinct from feces

The surface of beef cattle feedlot pens is commonly conceptualized as being packed uncomposted manure. Despite the important role that the feedlot pen may play in the transmission of veterinary and zoonotic pathogens, the bacterial ecology of feedlot surface material is not well understood. Our present study characterized the bacterial communities of the beef cattle feedlot pen surface material using 3647 full-length 16S rDNA sequences, and we compared the community composition of feedlot pens to the fecal source material. The feedlot surface composite was represented by members of the phylum Actinobacteria (42%), followed by Firmicutes (24%), Bacteroidetes (24%), and Proteobacteria (9%). The feedlot pen surface material bacterial communities were clearly distinct from those of the feces from animals in the same pen. Comparisons with previously published results of feces from the animals in the same pen reveal that, of 139 genera identified, only 25 were present in both habitats. These results indicate that, microbiologically, the feedlot pen surface material is separate and distinct from the fecal source material, suggesting that bacteria that originate in cattle feces face different selection pressures and survival challenges during their tenure in the feedlot pen, as compared to their residence in the gastrointestinal tract.

Virulence-associated and antibiotic resistance genes of microbial populations in cattle feces analyzed using a metagenomic approach

The bovine fecal microbiota impacts human food safety as well as animal health. Although the bacteria of cattle feces have been well characterized using culture-based and culture-independent methods, techniques have been lacking to correlate total community composition with community function. We used high throughput sequencing of total DNA extracted from fecal material to characterize general community composition and examine the repertoire of microbial genes present in beef cattle feces, including genes associated with antibiotic resistance and bacterial virulence. Results suggest that traditional 16S sequencing using "universal" primers to generate full-length sequence may under represent Acitinobacteria and Proteobacteria. Over eight percent (8.4%) of the sequences from our beef cattle fecal pool sample could be categorized as virulence genes, including a suite of genes associated with resistance to antibiotic and toxic compounds (RATC). This is a higher proportion of virulence genes found in Sargasso sea, chicken cecum, and cow rumen samples, but comparable to the proportion found in Antarctic marine derived lake, human fecal, and farm soil samples. The quantitative nature of metagenomic data, combined with the large number of RATC classes represented in samples from widely different habitats indicates that metagenomic data can be used to track relative amounts of antibiotic resistance genes in individual animals over time. Consequently, these data can be used to generate sample-specific and temporal antibiotic resistance gene profiles to facilitate an understanding of the ecology of the microbial communities in each habitat as well as the epidemiology of antibiotic resistant gene transport between and among habitats.

Genotyping Escherichia coli O157:H7 for its ability to cause disease in humans

Escherichia coli are ubiquitous in the world, and for the most part are non-pathogenic and part of the normal lower gastrointestinal tract in mammals. However, some pathogenic isolates can cause severe disease that range from meningitis to hemorrhagic colitis (HC). In recent years, Shiga toxin-containing E. coli (STEC) have been a major cause of food borne and environmental cases of HC and hemolytic uremic syndrome. One STEC serotype, O157:H7, has been responsible for numerous food-associated outbreaks and recalls worldwide. The protocols in this unit will allow the reader to use real-time polymerase chain reaction genotyping to identify isolates that are more likely to cause disease in humans. The genotyping assay targets a single-nucleotide polymorphism (SNP) in the tir gene. The tir gene is located in a virulence operon called the locus for enterocyte effacement and functions as a receptor for the tight adherence of E. coli O157:H7 to epithelial cells. As more genomes are sequenced, informative SNPs that associate with phenotypes will be identified. Identifying isolates not only by their genus and species, but also by using other informative genomic traits will increase the general knowledge about their genetic diversity.

Prevalence and enumeration of Escherichia coli O157:H7 and Salmonella in U.S. abattoirs that process fewer than 1000 head of cattle per day

A significant portion (15 to 20%) of beef in the United States is produced in small beef processing plants that harvest fewer than 1000 cattle per day. However, there are little data on the prevalence and levels of Escherichia coli O157:H7 and Salmonella in these processing plants. To address this lack of data, hides (n=1995) and carcasses (n=1995) of cattle at seven small processing plants located across the United States were analyzed for E. coli O157:H7 and Salmonella. Across all plants, hide prevalence of E. coli O157:H7 and Salmonella was 71 and 91%, respectively. Twelve percent of hides had E. coli O157:H7 at enumerable levels (> or =40 CFU/100 cm2), while 36% of hides had Salmonella at enumerable levels. Across all plants, the prevalence of E. coli O157:H7 on preevisceration carcasses was 33%, with 2% at an enumerable level (> or = 0.8 CFU/ 100 cm2). Across all plants, Salmonella prevalence on preevisceration carcasses was 58%, with 8% at an enumerable level. Significant plant-to-plant variations in levels and prevalence of pathogens on carcasses were detected. Reduced levels of pathogens on carcasses were noted among small processors that had incorporated a hide-directed intervention. The results obtained are comparable to those observed previously for larger processors, showing that smaller beef processors face and address the same challenges as do larger beef processors.