Phylogenetic structure of Shiga toxin-producing Escherichia coli O157:H7 from sub-lineage to SNPs

Re-analysis of an outbreak of Shiga toxin-producing Escherichia coli O157:H7 associated with raw drinking milk using Nanopore sequencing

The aim of this study was to compare Illumina and Oxford Nanopore Technology (ONT) sequencing data to quantify genetic variation to assess within-outbreak strain relatedness and characterise microevolutionary events in the accessory genomes of a cluster of 23 genetically and epidemiologically linked isolates related to an outbreak of Shiga toxin-producing Escherichia coli O157:H7 caused by the consumption of raw drinking milk. There were seven discrepant variants called between the two technologies, five were false-negative or false-positive variants in the Illumina data and two were false-negative calls in ONT data. After masking horizontally acquired sequences such as prophages, analysis of both short and long-read sequences revealed the 20 isolates linked to the outbreak in 2017 had a maximum SNP distance of one SNP between each other, and a maximum of five SNPs when including three additional strains identified in 2019. Analysis of the ONT data revealed a 47 kbp deletion event in a terminal compound prophage within one sample relative to the remaining samples, and a 0.65 Mbp large chromosomal rearrangement (inversion), within one sample relative to the remaining samples. Furthermore, we detected two bacteriophages encoding the highly pathogenic Shiga toxin (Stx) subtype, Stx2a. One was typical of Stx2a-phage in this sub-lineage (Ic), the other was atypical and inserted into a site usually occupied by Stx2c-encoding phage. Finally, we observed an increase in the size of the pO157 IncFIB plasmid (1.6 kbp) in isolates from 2019 compared to those from 2017, due to the duplication of insertion elements within the plasmids from the more recently isolated strains. The ability to characterize the accessory genome in this way is the first step to understanding the significance of these microevolutionary events and their impact on the genome plasticity and virulence between strains of this zoonotic, foodborne pathogen.

Clinical and public health implications of increasing notifications of LEE-negative Shiga toxin-producing Escherichia coli in England, 2014-2022

Introduction. Shiga toxin-producing Escherichia coli (STEC) belong to a diverse group of gastrointestinal pathogens. The pathogenic potential of STEC is enhanced by the presence of the pathogenicity island called the Locus of Enterocyte Effacement (LEE), including the intimin encoding gene eae.Gap statement. STEC serotypes O128:H2 (Clonal Complex [CC]25), O91:H14 (CC33), and O146:H21 (CC442) are consistently in the top five STEC serotypes isolated from patients reporting gastrointestinal symptoms in England. However, they are eae/LEE-negative and perceived to be a low risk to public health, and we know little about their microbiology and epidemiology.Aim. We analysed clinical outcomes and genome sequencing data linked to patients infected with LEE-negative STEC belonging to CC25 (O128:H2, O21:H2), CC33 (O91:H14) and, and CC442 (O146:H21, O174:H21) in England to assess the risk to public health.Results. There was an almost ten-fold increase between 2014 and 2022 in the detection of all STEC belonging to CC25, CC33 and CC442 (2014 n=38, 2022 n=336), and a total of 1417 cases. There was a higher proportion of female cases (55-70 %) and more adults than children, with patients aged between 20-40 and >70 most at risk across the different serotypes. Symptoms were consistent across the three dominant serotypes O91:H14 (CC33), O146:H21 (CC442) and O128:H2 (CC25) (diarrhoea >75 %; bloody diarrhoea 25-32 %; abdominal pain 64-72 %; nausea 37-45 %; vomiting 10-24 %; and fever 27-30 %). Phylogenetic analyses revealed multiple events of acquisition and loss of different stx-encoding prophage. Additional putative virulence genes were identified including iha, agn43 and subA.Conclusions. Continued monitoring and surveillance of LEE-negative STEC infections is essential due to the increasing burden of infectious intestinal disease, and the risk that highly pathogenic strains may emerge following acquisition of the Shiga toxin subtypes associated with the most severe clinical outcomes.

Persistence of microbiological hazards in food and feed production and processing environments

Listeria monocytogenes (in the meat, fish and seafood, dairy and fruit and vegetable sectors), Salmonella enterica (in the feed, meat, egg and low moisture food sectors) and Cronobacter sakazakii (in the low moisture food sector) were identified as the bacterial food safety hazards most relevant to public health that are associated with persistence in the food and feed processing environment (FFPE). There is a wide range of subtypes of these hazards involved in persistence in the FFPE. While some specific subtypes are more commonly reported as persistent, it is currently not possible to identify universal markers (i.e. genetic determinants) for this trait. Common risk factors for persistence in the FFPE are inadequate zoning and hygiene barriers; lack of hygienic design of equipment and machines; and inadequate cleaning and disinfection. A well-designed environmental sampling and testing programme is the most effective strategy to identify contamination sources and detect potentially persistent hazards. The establishment of hygienic barriers and measures within the food safety management system, during implementation of hazard analysis and critical control points, is key to prevent and/or control bacterial persistence in the FFPE. Once persistence is suspected in a plant, a 'seek-and-destroy' approach is frequently recommended, including intensified monitoring, the introduction of control measures and the continuation of the intensified monitoring. Successful actions triggered by persistence of L. monocytogenes are described, as well as interventions with direct bactericidal activity. These interventions could be efficient if properly validated, correctly applied and verified under industrial conditions. Perspectives are provided for performing a risk assessment for relevant combinations of hazard and food sector to assess the relative public health risk that can be associated with persistence, based on bottom-up and top-down approaches. Knowledge gaps related to bacterial food safety hazards associated with persistence in the FFPE and priorities for future research are provided.

Genomic epidemiology of the clinically dominant clonal complex 1 in the Listeria monocytogenes population in the UK

Listeria monocytogenes is a food-borne pathogen, typically affecting the elderly, immunocompromised patients and pregnant women. The aim of this study was to determine the population structure of L. monocytogenes clonal complex 1 (CC1) in the UK and describe the genomic epidemiology of this clinically significant CC. We interrogated a working dataset of 4073 sequences of L. monocytogenes isolated between January 2015 and December 2020 from human clinical specimens, food and/or food-production environments. A minimum spanning tree was reconstructed to determine the population structure of L. monocytogenes in the UK. Subsequent analysis focused on L. monocytogenes CC1, as the cause of the highest proportion of invasive listeriosis in humans. Sequencing data was integrated with metadata on food and environmental isolates, and information from patient questionnaires, including age, sex and clinical outcomes. All isolates either belonged to lineage I (n=1299/4073, 32%) or lineage II (n=2774/4073, 68%), with clinical isolates from human cases more likely to belong to lineage I (n=546/928, 59%) and food isolates more likely to belong to lineage II (n=2352/3067, 77%). Of the four largest CCs, CC1 (n=237) had the highest proportion of isolates from human cases of disease (CC1 n=160/237, 67.5 %; CC121 n=13/843, 2 %; CC9 n=53/360, 15 %; CC2 n=69/339, 20%). Within CC1, most cases were female (n=95/160, 59%, P=0.01771) and the highest proportion of cases were in people >60 years old (39/95, 41%, P=1.314×10-6) with a high number of them aged 20-39 years old (n=35/95, 37%) most linked to pregnancy-related listeriosis (n=29/35, 83%). Most of the male cases were in men aged over 60 years old (40/65, 62%), and most of the fatal cases in both males and females were identified in this age group (42/55, 76%). Phylogenetic analysis revealed 23 5 SNP single linkage clusters comprising 80/237 (34 %) isolates with cluster sizes ranging from 2 to 19. Five 5 SNP clusters comprised isolates from human cases and an implicated food item. Expanding the analysis to 25 SNP single linkage clusters resolved an additional two clusters linking human cases to a potential food vehicle. Analysis of demographic and clinical outcome data identified CC1 as a clinically significant cause of invasive listeriosis in the elderly population and in women of child-bearing age. Phylogenetic analysis revealed the population structure of CC1 in the UK comprised small, sparsely populated genomic clusters. Only clusters containing isolates from an implicated food vehicle, or food processing or farming environments, were resolved, emphasizing the need for clinical, food and animal-health agencies to share sequencing data in real time, and the importance of a One Health approach to public-health surveillance of listeriosis.

Analysis of Escherichia coli O157 strains in cattle and humans between Scotland and England & Wales: implications for human health

For the last two decades, the human infection frequency of Escherichia coli O157 (O157) in Scotland has been 2.5-fold higher than in England and Wales. Results from national cattle surveys conducted in Scotland and England and Wales in 2014/2015 were combined with data on reported human clinical cases from the same time frame to determine if strain differences in national populations of O157 in cattle could be associated with higher human infection rates in Scotland. Shiga toxin subtype (Stx) and phage type (PT) were examined within and between host (cattle vs human) and nation (Scotland vs England and Wales). For a subset of the strains, whole genome sequencing (WGS) provided further insights into geographical and host association. All three major O157 lineages (I, II, I/II) and most sub-lineages (Ia, Ib, Ic, IIa, IIb, IIc) were represented in cattle and humans in both nations. While the relative contribution of different reservoir hosts to human infection is unknown, WGS analysis indicated that the majority of O157 diversity in human cases was captured by isolates from cattle. Despite comparable cattle O157 prevalence between nations, strain types were localized. PT21/28 (sub-lineage Ic, Stx2a+) was significantly more prevalent in Scottish cattle [odds ratio (OR) 8.7 (2.3-33.7; P<0.001] and humans [OR 2.2 (1.5-3.2); P<0.001]. In England and Wales, cattle had a significantly higher association with sub-lineage IIa strains [PT54, Stx2c; OR 5.6 (1.27-33.3); P=0.011] while humans were significantly more closely associated with sub-lineage IIb [PT8, Stx1 and Stx2c; OR 29 (4.9-1161); P<0.001]. Therefore, cattle farms in Scotland were more likely to harbour Stx2a+O157 strains compared to farms in E and W (P<0.001). There was evidence of limited cattle strain migration between nations and clinical isolates from one nation were more similar to cattle isolates from the same nation, with sub-lineage Ic (mainly PT21/28) exhibiting clear national association and evidence of local transmission in Scotland. While we propose the higher rate of O157 clinical cases in Scotland, compared to England and Wales, is a consequence of the nationally higher level of Stx2a+O157 strains in Scottish cattle, we discuss the multiple additional factors that may also contribute to the different infection rates between these nations.

Whole-Genome Analysis of Bordetella pertussis MT27 Isolates from School-Associated Outbreaks: Single-Nucleotide Polymorphism Diversity and Threshold of the Outbreak Strains

Bordetella pertussis, the causative agent of whooping cough, can cause pertussis outbreaks in humans, especially in school-aged children. Here, we performed whole-genome sequencing of 51 B. pertussis isolates (epidemic strain MT27) collected from patients infected during 6 school-associated outbreaks lasting less than 4 months. We compared their genetic diversity with that of 28 sporadic isolates (non-outbreak MT27 isolates) based on single-nucleotide polymorphisms (SNPs). Our temporal SNP diversity analysis revealed a mean SNP accumulation rate (time-weighted average) of 0.21 SNPs/genome/year during the outbreaks. The outbreak isolates showed a mean of 0.74 SNP differences (median, 0; range, 0 to 5) between 238 isolate pairs, whereas the sporadic isolates had a mean of 16.12 SNP differences (median, 17; range 0 to 36) between 378 isolate pairs. A low SNP diversity was observed in the outbreak isolates. Receiver operating characteristic analysis demonstrated that the optimal cutoff value to distinguish between the outbreak and sporadic isolates was 3 SNPs (Youden's index of 0.90 with a true-positive rate of 0.97 and a false-positive rate of 0.07). Based on these results, we propose an epidemiological threshold of ≤3 SNPs per genome as a reliable marker of B. pertussis strain identity during pertussis outbreaks that span less than 4 months. IMPORTANCE Bordetella pertussis is a highly infectious bacterium that easily causes pertussis outbreaks in humans, especially in school-aged children. In detection and investigation of outbreaks, excluding non-outbreak isolates is important for understanding the bacterial transmission routes. Currently, whole-genome sequencing is widely used for outbreak investigations, and the genetic relatedness of outbreak isolates is assessed based on differences in the number of single-nucleotide polymorphisms (SNPs) in the genomes of different isolates. The optimal SNP threshold defining strain identity has been proposed for many bacterial pathogens, but not for B. pertussis. In this study, we performed whole-genome sequencing of 51 B. pertussis outbreak isolates and identified a genetic threshold of ≤3 SNPs per genome as a marker defining the strain identity during pertussis outbreaks. This study provides a useful marker for identifying and analyzing pertussis outbreaks and can serve as a basis for future epidemiological studies on pertussis.

Two Outbreaks of Foodborne Gastrointestinal Infection Linked to Consumption of Imported Melons, United Kingdom, March to August 2021

The aim of this study was to describe two foodborne outbreaks caused by contaminated imported melon and make recommendations for future practice. Between March and July 2021, there was an outbreak of 113 cases of Salmonella Braenderup in the UK (62% female, median age 61 years, 33% hospitalized). Analytical epidemiological studies identified Galia melons as the vehicle of infection (OR 671.9, 95% CI 39.0-58,074.0, p < 0.001). Subsequently, the outbreak strain was isolated from two samples of Galia melon imported from Latin America. In July and August 2021, there was an outbreak of 17 cases of Shiga toxin-producing Escherichia coli (STEC) O157:H7 in the UK (53% female, median age 21 years, 35% were hospitalized). Review of the STEC surveillance questionnaire data, followed by the analysis of responses from a modified hypothesis-generating questionnaire, implicated eating precut watermelon from retailer B sourced from Europe as the vehicle of infection. Outbreaks of gastrointestinal pathogens caused by contaminated food of nonanimal origin are a global public health concern. Given the difficulty in removing pathogens from the flesh of ready-to-eat fruit and vegetables, public health interventions should target all steps of the food chain prior to consumption, from cultivation on the farm to processing/packing and distribution.

Distinct evolutionary trajectories in the Escherichia coli pangenome occur within sequence types

The Escherichia coli species contains a diverse set of sequence types and there remain important questions regarding differences in genetic content within this population that need to be addressed. Pangenomes are useful vehicles for studying gene content within sequence types. Here, we analyse 21 E. coli sequence type pangenomes using comparative pangenomics to identify variance in both pangenome structure and content. We present functional breakdowns of sequence type core genomes and identify sequence types that are enriched in metabolism, transcription and cell membrane biogenesis genes. We also uncover metabolism genes that have variable core classification, depending on which allele is present. Our comparative pangenomics approach allows for detailed exploration of sequence type pangenomes within the context of the species. We show that ongoing gene gain and loss in the E. coli pangenome is sequence type-specific, which may be a consequence of distinct sequence type-specific evolutionary drivers.

Whole genome sequencing and characteristics of extended-spectrum beta-lactamase producing Escherichia coli isolated from poultry farms in Banaskantha, India

Worldwide dissemination of extended-spectrum -lactamase (ESBL)-producing Escherichia coli constitutes an emerging global health issue, with animal food products contributing as potential reservoirs. ESBL E. coli infection is associated with the high mortality and mobility rate in developing countries due to less susceptibility to antibiotics. The present study aimed to elucidate the molecular characteristics and sequence-based analysis of ESBL E. coli in the Gujarat state of India. This study included 108 E. coli strains were isolated from different poultry farms (broiler and layer) in the Banaskantha District. PCR was employed to identify genotypic ESBL-producing antimicrobial resistance genes. Overall, a high occurrence of ESBL genes was found in poultry farms due to the high usage of antimicrobials. The PCR analysis revealed that 79.62% of isolates were detected positive with one or more ESBL genes. Among them, bla_TEM (63.88%) was found to be the predominant genotype, followed by bla_SHV (30.55%) and bla_OXA (28.70%). In the bla_CTX-M group, a higher occurrence was observed in bla_CTX-M-9 (23.14%), followed by bla_CTX-M-2 (24.07%) and bla_CTX-M-1 (22.22%). We used the whole-genome sequencing (WGS) method to evaluate the antimicrobial resistance genes, virulence factors, single nucleotide polymorphisms (SNPs), plasmid replicons, and plasmid-mediated AMR genes of one ESBL E. coli isolated. We examined the genetic relatedness of a human pathogenic E. coli strain by comparing its sequence with the broad geographical reference E. coli sequences. Escherichia coli ST 681 was determined using multi-locus sequence typing. We compared our findings to the reference sequence of Escherichia coli str. K- 12 substr. MG1655. We found 24,937 SNPs with 21,792 in the genic region, 3,145 in the intergenic region, and six InDels across the genome. The WGS analysis revealed 46 antimicrobial resistance genes and seven plasmid-mediated AMR genes viz., tetA, qnrS1, dfrA14, sul2, aph(3”)-lb, aph(6)-ld, and Aph(3’)-la. The ST 681 was found to have Cib, traT, and terC virulence factors and two plasmid replicons, IncFII(pHN7A8) and IncI1-I(Alpha). This study revealed a higher occurrence of ESBL E. coli detected in poultry.

Outbreak of STEC O157:H7 linked to a milk pasteurisation failure at a dairy farm in England, 2019

In November 2019, an outbreak of Shiga toxin-producing Escherichia coli O157:H7 was detected in South Yorkshire, England. Initial investigations established consumption of milk from a local dairy as a common exposure. A sample of pasteurised milk tested the next day failed the phosphatase test, indicating contamination of the pasteurised milk by unpasteurised (raw) milk. The dairy owner agreed to immediately cease production and initiate a recall. Inspection of the pasteuriser revealed a damaged seal on the flow divert valve. Ultimately, there were 21 confirmed cases linked to the outbreak, of which 11 (52%) were female, and 12/21 (57%) were either <15 or >65 years of age. Twelve (57%) patients were treated in hospital, and three cases developed haemolytic uraemic syndrome. Although the outbreak strain was not detected in the milk samples, it was detected in faecal samples from the cattle on the farm. Outbreaks of gastrointestinal disease caused by milk pasteurisation failures are rare in the UK. However, such outbreaks are a major public health concern as, unlike unpasteurised milk, pasteurised milk is marketed as ‘safe to drink’ and sold to a larger, and more dispersed, population. The rapid, co-ordinated multi-agency investigation initiated in response to this outbreak undoubtedly prevented further cases.

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

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

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

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

Outbreak of Shiga toxin-producing Escherichia coli O157 linked with consumption of a fast-food product containing imported cucumbers, United Kingdom, August 2020

BACKGROUND

In August 2020, an outbreak of Shiga toxin-producing Escherichia coli (STEC) O157:H7 occurred in the United Kingdom. Whole genome sequencing revealed that these cases formed a genetically distinct cluster.

METHODS

Hypotheses generated from case interviews were tested in analytical studies, and results informed environmental sampling and food chain analysis. A case-case study used non-outbreak 'comparison' STEC cases; a case-control study used a market research panel to recruit controls.

RESULTS

A total of 36 cases were identified; all cases reported symptom onset between August 3 and August 16, 2020. The majority of cases (83%) resided in the Midlands region of England and in Wales. A high proportion of cases reported eating out, with one fast-food restaurant chain mentioned by 64% (n = 23) of cases. Both the case-case study (adjusted odds ratio (aOR) 31.8, 95% confidence interval (CI) 1.6-624.9) and the case-control study (aOR 9.19, 95% CI 1.0-82.8) revealed statistically significant results, showing that the consumption of a specific fast-food product was independently associated with infection.

CONCLUSIONS

Consumption of a specific fast-food product was a likely cause of this outbreak. The only ingredient specific to the product was cucumbers. The supply of cucumbers was immediately halted, and no further cases have been identified.