Use of Genomics to Investigate Historical Importation of Shiga Toxin-Producing Escherichia coli Serogroup O26 and Nontoxigenic Variants into New Zealand

High-resolution genomic analysis to investigate the impact of the invasive brushtail possum (Trichosurus vulpecula) and other wildlife on microbial water quality assessments

Escherichia coli are routine indicators of fecal contamination in water quality assessments. Contrary to livestock and human activities, brushtail possums (Trichosurus vulpecula), common invasive marsupials in Aotearoa/New Zealand, have not been thoroughly studied as a source of fecal contamination in freshwater. To investigate their potential role, Escherichia spp. isolates (n = 420) were recovered from possum gut contents and feces and were compared to those from water, soil, sediment, and periphyton samples, and from birds and other introduced mammals collected within the Mākirikiri Reserve, Dannevirke. Isolates were characterized using E. coli-specific real-time PCR targeting the uidA gene, Sanger sequencing of a partial gnd PCR product to generate a gnd sequence type (gST), and for 101 isolates, whole genome sequencing. Escherichia populations from 106 animal and environmental sample enrichments were analyzed using gnd metabarcoding. The alpha diversity of Escherichia gSTs was significantly lower in possums and animals compared with aquatic environmental samples, and some gSTs were shared between sample types, e.g., gST535 (in 85% of samples) and gST258 (71%). Forty percent of isolates gnd-typed and 75% of reads obtained by metabarcoding had gSTs shared between possums, other animals, and the environment. Core-genome single nucleotide polymorphism (SNP) analysis showed limited variation between several animal and environmental isolates (<10 SNPs). Our data show at an unprecedented scale that Escherichia clones are shared between possums, other wildlife, water, and the wider environment. These findings support the potential role of possums as contributors to fecal contamination in Aotearoa/New Zealand freshwater. Our study deepens the current knowledge of Escherichia populations in under-sampled wildlife. It presents a successful application of high-resolution genomic methods for fecal source tracking, thereby broadening the analytical toolbox available to water quality managers. Phylogenetic analysis of isolates and profiling of Escherichia populations provided useful information on the source(s) of fecal contamination and suggest that comprehensive invasive species management strategies may assist in restoring not only ecosystem health but also water health where microbial water quality is compromised.

Phylogenetic relationship and virulence composition of Escherichia coli O26:H11 cattle and human strain collections in Scotland; 2002-2020

O26 is the commonest non-O157 Shiga toxin (stx)-producing Escherichia coli serogroup reported in human infections worldwide. Ruminants, particularly cattle, are the primary reservoir source for human infection. In this study, we compared the whole genomes and virulence profiles of O26:H11 strains (n = 99) isolated from Scottish cattle with strains from human infections (n = 96) held by the Scottish Escherichia coli O157/STEC Reference Laboratory, isolated between 2002 and 2020. Bovine strains were from two national cross-sectional cattle surveys conducted between 2002-2004 and 2014-2015. A maximum likelihood phylogeny was constructed from a core-genome alignment with the O26:H11 strain 11368 reference genome. Genomes were screened against a panel of 2,710 virulence genes using the Virulence Finder Database. All stx-positive bovine O26:H11 strains belonged to the ST21 lineage and were grouped into three main clades. Bovine and human source strains were interspersed, and the stx subtype was relatively clade-specific. Highly pathogenic stx2a-only ST21 strains were identified in two herds sampled in the second cattle survey and in human clinical infections from 2010 onwards. The closest pairwise distance was 9 single-nucleotide polymorphisms (SNPs) between Scottish bovine and human strains and 69 SNPs between the two cattle surveys. Bovine O26:H11 was compared to public EnteroBase ST29 complex genomes and found to have the greatest commonality with O26:H11 strains from the rest of the UK, followed by France, Italy, and Belgium. Virulence profiles of stx-positive bovine and human strains were similar but more conserved for the stx2a subtype. O26:H11 stx-negative ST29 (n = 17) and ST396 strains (n = 5) were isolated from 19 cattle herds; all were eae-positive, and 10 of these herds yielded strains positive for ehxA, espK, and Z2098, gene markers suggestive of enterohaemorrhagic potential. There was a significant association (p < 0.001) between nucleotide sequence percent identity and stx status for the bacteriophage insertion site genes yecE for stx2 and yehV for stx1. Acquired antimicrobial resistance genes were identified in silico in 12.1% of bovine and 17.7% of human O26:H11 strains, with sul2, tet, aph(3″), and aph(6″) being most common. This study describes the diversity among Scottish bovine O26:H11 strains and investigates their relationship to human STEC infections.

Occurrence and genetic characterization of Shiga toxin-producing Escherichia coli on bovine and pork carcasses and the environment from transport trucks

Shiga toxin-producing Escherichia coli (STEC) are foodborne pathogens causing severe diseases. The ability of STEC to produce disease is associated with Shiga toxin (Stx) production. We investigated the occurrence of STEC on bovine and pork carcasses and walls of trucks where they were transported, and we characterized virulence genes and serotypes of STEC strains. We compared the whole genomic sequencing of a STEC O157:H7 strain isolated from a bovine carcass in this work and a STEC O157:H7 strain isolated from a child with HUS, both isolated in 2019. We studied the relationship between these isolates and others collected in the database. The results show a 40% of STEC and two different serogroups were identified (O130 and O157). STEC O157:H7 were isolated from bovine carcasses and harbored stx2, eae, ehxA, katP, espP, stcE, ECSP_0242/1773/2687/2870/2872/3286/3620 and were classified as lineage I/II. In STEC non-O157 isolates, three isolates were isolated from bovine carcasses and harbored the serogroup O130 and one strain isolated from pork carcasses was O-non-typeable. All STEC non-O157 harbored sxt1 gene. The analysis from the whole genome showed that both STEC O157:H7 strains belonged to the hypervirulent clade 8, ST11, phylogroup E, carried the allele tir 255 T > A T, and they were not clonal. The analysis of information allows us to conclude that the STEC strains circulate in pork and bovine carcasses arriving in transport. This situation represents a risk for the consumers and the need to implement an integrated STEC control in the food chain.

Genetic and phylogenetic characterization of Shiga toxin-producing Escherichia coli and enteropathogenic E. coli from livestock in Jiangsu by using whole-genome sequencing

AIMS

There are knowledge gaps regarding STEC and EPEC strains in livestock in Jiangsu, China. This study aimed to evaluate the potential public health significance of STEC and EPEC strains isolated from livestock by determining the serotypes, virulence profiles, and genetic relationship with international STEC strains.

METHODS AND RESULTS

A total of 68 STEC and 37 EPEC strains were obtained from 231 fecal sheep samples and 70 fecal cattle samples. By using whole-genome sequencing (WGS) analysis, all STEC belonged to 15 O:H serotypes and the most prevalent serotypes were O6:H10 (19.1%), O155:H21 (14.7%), and O21:H25 (10.3%). The main Shiga toxin gene subtypes detected were stx1c (41.2%), stx1a (26.5%), stx2b (14.7%) and stx2k (14.7%). Only the STEC from cattle carried eae gene. Other adherence-associated or toxin-related genes, including lpfA (70.6%), iha (48.5%), subA (54.4%), and ehxA (33.8%), were found in STEC. All EPEC strains were bfpA-negative, and the predominant eae variants were eae-β1 (62.2%), eae-ζ (21.6%), and eae-θ (8.1%). The core-genome multi-locus sequence typing (cgMLST) analysis revealed nine scattered clusters in STEC and one dominant cluster in EPEC. The strains with the same serotypes, including O22:H8 and O43:H2 in the two towns, possessed a closely genomic distance. The core genome single nucleotide polymorphism (cgSNP) showed that part of STEC strains in this study were clustered with isolates possessing the same serotypes from the Netherlands, Sweden, and Xinjiang of China. Five serotypes of STEC isolates were associated with the clinical STEC strains from databases.

CONCLUSION

This study provided the diverse serotypes and the virulence genes profiles in STEC and EPEC strains. Local strains possessed widely diverse and scattered clusters by cgMLST. Closely genomic correlation with clinical isolates displayed that part of the STEC strains may threaten to public health.

SIGNIFICANCE AND IMPACT OF THE STUDY

Non-O157 STEC strains act as important pathogens for human infections. This study supports the increased surveillance work of non-O157 STEC rather than just O157 STEC in this region.

Whole-Genome Phylogenetic Analysis Reveals a Wide Diversity of Non-O157 STEC Isolated From Ground Beef and Cattle Feces

Shiga toxin-producing Escherichia coli (STEC) causes foodborne outbreaks that can lead to complications such as hemolytic uremic syndrome. Their main reservoir is cattle, and ground beef has been frequently associated with disease and outbreaks. In this study, we attempted to understand the genetic relationship among STEC isolated in Chile from different sources, their relationship to STEC from the rest of the world, and to identify molecular markers of Chilean STEC. We sequenced 62 STEC isolated in Chile using MiSeq Illumina. In silico typing was determined using tools of the Center Genomic Epidemiology, Denmark University (CGE/DTU). Genomes of our local STEC collection were compared with 113 STEC isolated worldwide through a core genome MLST (cgMLST) approach, and we also searched for distinct genes to be used as molecular markers of Chilean isolates. Genomes in our local collection were grouped based on serogroup and sequence type, and clusters were formed within local STEC. In the worldwide STEC analysis, Chilean STEC did not cluster with genomes of the rest of the world suggesting that they are not phylogenetically related to previously described STEC. The pangenome of our STEC collection was 11,650 genes, but we did not identify distinct molecular markers of local STEC. Our results showed that there may be local emerging STEC with unique features, nevertheless, no molecular markers were detected. Therefore, there might be elements such as a syntenic organization that might explain differential clustering detected between local and worldwide STEC.

Genomic epidemiology and carbon metabolism of Escherichia coli serogroup O145 reflect contrasting phylogenies

Shiga toxin-producing Escherichia coli (STEC) are a leading cause of foodborne outbreaks of human disease, but they reside harmlessly as an asymptomatic commensal in the ruminant gut. STEC serogroup O145 are difficult to isolate as routine diagnostic methods are unable to distinguish non-O157 serogroups due to their heterogeneous metabolic characteristics, resulting in under-reporting which is likely to conceal their true prevalence. In light of these deficiencies, the purpose of this study was a twofold approach to investigate enhanced STEC O145 diagnostic culture-based methods: firstly, to use a genomic epidemiology approach to understand the genetic diversity and population structure of serogroup O145 at both a local (New Zealand) (n = 47) and global scale (n = 75) and, secondly, to identify metabolic characteristics that will help the development of a differential media for this serogroup. Analysis of a subset of E. coli serogroup O145 strains demonstrated considerable diversity in carbon utilisation, which varied in association with eae subtype and sequence type. Several carbon substrates, such as D-serine and D-malic acid, were utilised by the majority of serogroup O145 strains, which, when coupled with current molecular and culture-based methods, could aid in the identification of presumptive E. coli serogroup O145 isolates. These carbon substrates warrant subsequent testing with additional serogroup O145 strains and non-O145 strains. Serogroup O145 strains displayed extensive genetic heterogeneity that was correlated with sequence type and eae subtype, suggesting these genetic markers are good indicators for distinct E. coli phylogenetic lineages. Pangenome analysis identified a core of 3,036 genes and an open pangenome of >14,000 genes, which is consistent with the identification of distinct phylogenetic lineages. Overall, this study highlighted the phenotypic and genotypic heterogeneity within E. coli serogroup O145, suggesting that the development of a differential media targeting this serogroup will be challenging.

Whole Genome Sequencing Differentiates Presumptive Extended Spectrum Beta-Lactamase Producing Escherichia coli along Segments of the One Health Continuum

Antimicrobial resistance (AMR) has important implications for the continued use of antibiotics to control infectious diseases in both beef cattle and humans. AMR along the One Health continuum of the beef production system is largely unknown. Here, whole genomes of presumptive extended-spectrum β-lactamase E. coli (ESBL-EC) from cattle feces (n = 40), feedlot catch basins (n = 42), surrounding streams (n = 21), a beef processing plant (n = 4), municipal sewage (n = 30), and clinical patients (n = 25) are described. ESBL-EC were isolated from ceftriaxone selective plates and subcultured on ampicillin selective plates. Agreement of genotype-phenotype prediction of AMR ranged from 93.2% for ampicillin to 100% for neomycin, trimethoprim/sulfamethoxazole, and enrofloxacin resistance. Overall, β-lactam (100%; bla_EC, bla_TEM-1, bla_SHV, bla_OXA, bla_CTX-M-), tetracycline (90.1%; tet(A), tet(B)) and folate synthesis (sul2) antimicrobial resistance genes (ARGs) were most prevalent. The ARGs tet(C), tet(M), tet(32),bla_CTX-M-1, bla_CTX-M-14, bla_OXA-1, dfrA18, dfrA19, catB3, and catB4 were exclusive to human sources, while bla_TEM-150, bla_SHV-11–12,dfrA12, cmlA1, and cmlA5 were exclusive to beef cattle sources. Frequently encountered virulence factors across all sources included adhesion and type II and III secretion systems, while IncFIB(AP001918) and IncFII plasmids were also common. Specificity and prevalence of ARGs between cattle-sourced and human-sourced presumptive ESBL-EC likely reflect differences in antimicrobial use in cattle and humans. Comparative genomics revealed phylogenetically distinct clusters for isolates from human vs. cattle sources, implying that human infections caused by ESBL-EC in this region might not originate from beef production sources.