Genomic analysis of an outbreak of Shiga toxin-producing Escherichia coli O183:H18 in the United Kingdom, 2023.

In June 2023, UKHSA surveillance systems detected an outbreak of severe gastrointestinal symptoms caused by a rare serotype of Shiga toxin-producing Escherichia coli, STEC O183:H18. There were 26 cases aged 6 months to 74 years (42 % cases were aged 0-9 years), distributed across the UK with onset dates range between 22 May 2023 and 4 July 2023. The epidemiological and food chain investigations were inconclusive, although meat products made from beef mince were implicated as a potential vehicle. The outbreak strain belonged to sequence type (ST) 657 and harboured a Shiga toxin (stx) subtype stx2a located on a prophage that was unique in the UKHSA stx-encoding bacteriophage database. Plasmid encoded, putative virulence genes subA, ehxA, saa, iha, lpfA and iss were detected, however, the established STEC virulence genes involved in attachment to the gut mucosa (eae and aggR) were absent. The acquisition of stx across the global population structure of ST657 appeared to correspond with the presence of subA, ehxA, saa, iha, lpfA and iss. During the outbreak investigation, we used long read sequencing to characterise the plasmid and prophage content of this atypical STEC, to look for evidence to explain its recent emergence. Although we were unable to determine source and transmission route of the outbreak strain, the genomic analysis revealed potential clues as to how novel strains for STEC evolve. With the implementation of PCR capable of detecting all STEC, and genome sequencing for typing and virulence profiling, we have the tools to enable us to monitor the changing landscape of STEC. Improvements in the standardised collection of epidemiological data and trace-back strategies within the food industry, will ensure we have a surveillance system capable of alerting us to emerging threats to public health.

Surveillance of antimicrobial resistant Shiga toxin-producing E. coli O157:H7 in England, 2016-2020.

OBJECTIVES: Shiga toxin-producing Escherichia coli (STEC) O157:H7 are zoonotic pathogens and transmission to humans occurs via contaminated food or contact with infected animals. The aim of this study was to describe the frequency, and distribution across the phylogeny, of antimicrobial resistance (AMR) determinants in STEC O157:H7 isolated from human cases in England. METHODS: Short-read whole-genome sequencing data from 1473 isolates of STEC O157:H7 from all seven sub-lineages (Ia-Ic, IIa-IIc and I/II) were mapped to genes known to confer phenotypic resistance to 10 different classes of antibiotic. Long-read sequencing was used to determine the location and genomic architecture of the AMR determinants within phylogenetic clusters exhibiting multidrug resistance. RESULTS: Overall, 216/1473 (14.7%) isolates had at least one AMR determinant, although the proportion of isolates exhibiting AMR varied by sub-lineage. The highest proportion of AMR determinants were detected in sub-lineages Ib (28/64, 43.7%), I/II (18/51, 35.3%) and IIc (122/440, 27.7%). In all sub-lineages, the most commonly detected AMR determinants conferred resistance to the aminoglycosides, tetracyclines and sulphonamides, while AMR determinants conferring resistance to fluroquinolones, macrolides and third-generation cephalosporins were rarely detected. Long-read sequencing analysis showed that the AMR determinants were co-located on the chromosome in sub-lineages Ib and lineage I/II, whereas those associated with sub-lineage IIc were encoded on the chromosome and/or large plasmids. CONCLUSIONS: AMR genes were unevenly distributed across the different sub-lineages of STEC O157:H7 and between different clades within the same sub-lineage. Long-read sequencing facilitates tracking the transmission of AMR at the pathogen and mobile genetic element level.

Escherichia coli encoding Shiga toxin subtype Stx2f causing human infections in England, 2015-2022.

Introduction. Shiga toxin-producing Escherichia coli (STEC) belong to a diverse group of gastrointestinal pathogens defined by the presence of Shiga toxin genes (stx) of which there are at least ten subtypes (Stx1a-Stx1d and Stx2a-Stx2g).Gap Statement. Initially thought to be associated with mild symptoms, more recently STEC encoding stx2f have been isolated from cases of haemolytic uraemic syndrome (HUS) and the clinical significance and public health burden require further investigation.Aim. We analysed clinical outcomes and genome-sequencing data linked to patients infected with STEC encoding-stx2f in England to assess the risk to public health.Methodology. One hundred and twelve E. coli (n=58 isolates encoded stx2f; n=54 isolates E. coli belonging to CC122 or CC722 that had eae but were negative for stx) isolated from patients' faecal specimens between 2015 and 2022 were genome sequenced and linked to epidemiological and clinical outcome data. All isolates were investigated for the presence of virulence genes and a maximum-likelihood phylogeny of isolates belonging to CC122 and CC722 was constructed.Results. There were 52 cases infected with STEC harbouring stx2f between 2015 and 2022, with the majority identified in 2022. Most cases resided in the North of England (n=39/52, 75 %), were female (n=31, 59.6 %) and/or aged five and under (n=29, 55.8 %). Clinical outcome data were available for 40/52 cases (76.9 %) and 7/40(17.5 %) were diagnosed with STEC-HUS. In the two most common clonal complexes, CC122 and CC722, the presence of the stx2f-encoding prophage correlated with the presence of additional virulence genes, astA, bfpA and cdt, located on an 85kbp IncFIB plasmid.Conclusions. Certain serotypes of E. coli harbouring stx2f cause severe clinical outcomes, including STEC-HUS. Public health advice and possible interventions are limited, as little is known about the animal and environmental reservoirs and transmission routes. We recommend more comprehensive and standardized collection of microbiological and epidemiological data, and routine sharing of sequencing data between public health agencies worldwide.

Microevolution of Campylobacter jejuni during long-term infection in an immunocompromised host.

Campylobacteriosis typically manifests as a short-lived, self-limiting gastrointestinal infection in humans, however prolonged infection can be seen in cases with underlying immunodeficiency. Public Health England received 25 isolates of Campylobacter jejuni from an individual with combined variable immunodeficiency over a period of 15 years. All isolates were typed and archived at the time of receipt. Whole genome sequencing (WGS) and antimicrobial susceptibility testing were performed to examine the relatedness of the isolates and to investigate the changes in the genome that had taken place over the course of the infection. Genomic typing methods were compared to conventional phenotypic methods, and revealed that the infection was caused by a single, persistent strain of C. jejuni belonging to clonal complex ST-45, with evidence of adaptation and selection in the genome over the course of the infection. Genomic analysis of sequence variants associated with antimicrobial resistance identified the genetic background behind rRNA gene mutations causing variable levels of resistance to erythromycin. This application of WGS to examine a persistent case of campylobacteriosis provides insight into the mutations and selective pressures occurring over the course of an infection, some of which have important clinical relevance.

Comparison of phenotypic and WGS-derived antimicrobial resistance profiles of Campylobacter jejuni and Campylobacter coli isolated from cases of diarrhoeal disease in England and Wales, 2015-16.

OBJECTIVES: To compare and evaluate phenotypic and genotypic methods for the detection of antimicrobial resistance (AMR) in Campylobacter jejuni and Campylobacter coli in England and Wales. METHODS: WGS data from 528 isolates of Campylobacter spp. (452 C. jejuni and 76 C. coli) from human (494), food (21) and environmental (2) sources, collected between January 2015 and December 2016, and from the PHE culture collection (11) were mapped to genes known to be associated with phenotypic resistance to antimicrobials in the genus. Phenotypic antibiotic susceptibility (erythromycin, ciprofloxacin, tetracycline, gentamicin and streptomycin) testing using an in-agar dilution method was performed on all isolates. RESULTS: Concordance between phenotypic resistance and the presence of corresponding AMR determinants was 97.5% (515/528 isolates). Only 13 out of 528 isolates (10 C. jejuni and 3 C. coli) had discordant interpretations for at least one of the five antibiotics tested, equating to a total of 15 (0.6%) discrepancies out of 2640 isolate/antimicrobial combinations. Seven discrepant results were genotypically resistant but phenotypically susceptible (major errors) and eight discrepant results were genotypically susceptible but phenotypically resistant (very major errors). CONCLUSIONS: The use of this bioinformatics approach for predicting AMR from WGS data for routine public health surveillance is a reliable method for real-time monitoring of changing AMR patterns in isolates of C. jejuni and C. coli.

Phylogenomic analysis of gastroenteritis-associated Clostridium perfringens in England and Wales over a 7-year period indicates distribution of clonal toxigenic strains in multiple outbreaks and extensive involvement of enterotoxin-encoding (CPE) plasmids.

Clostridium perfringens is a major enteric pathogen known to cause gastroenteritis in human adults. Although major outbreak cases are frequently reported, only limited whole-genome sequencing (WGS) based studies have been performed to understand the genomic epidemiology and virulence gene content of outbreak-associated C. perfringens strains. We performed phylogenomic analysis on 109 C. perfringens isolates (human and food) obtained from disease cases in England and Wales between 2011 and 2017. Initial findings highlighted the enhanced discriminatory power of WGS in profiling outbreak C. perfringens strains, when compared to the current Public Health England referencing laboratory technique of fluorescent amplified fragment length polymorphism analysis. Further analysis identified that isogenic C. perfringens strains were associated with nine distinct care-home-associated outbreaks over the course of a 5-year interval, indicating a potential common source linked to these outbreaks or transmission over time and space. As expected, the enterotoxin cpe gene was encoded in all but 4 isolates (96.3 %; 105/109), with virulence plasmids encoding cpe (particularly pCPF5603 and pCPF4969 plasmids) extensively distributed (82.6 %; 90/109). Genes encoding accessory virulence factors, such as beta-2 toxin, were commonly detected (46.7 %; 51/109), and genes encoding phage proteins were also frequently identified. Overall, this large-scale genomic study of gastroenteritis-associated C. perfringens suggested that three major cpe-encoding (toxinotype F) genotypes underlie these outbreaks: strains carrying (1) pCPF5603 plasmid, (2) pCPF4969 plasmid and (3) chromosomal-cpe strains. Our findings substantially expanded our knowledge on type F C. perfringens involved in human-associated gastroenteritis, with further studies required to fully probe the dissemination and regional reservoirs of this enteric pathogen, which may help devise effective prevention strategies to reduce the food-poisoning disease burden in vulnerable patients, such as the elderly.

Utility of Whole Genome Sequencing To Describe the Persistence and Evolution of Listeria monocytogenes Strains within Crabmeat Processing Environments Linked to Two Outbreaks of Listeriosis.

This article describes the identification and investigation of two extended outbreaks of listeriosis in which crabmeat was identified as the vehicle of infection. Comparing contemporary and retrospective typing data of Listeria monocytogenes isolates from clinical cases and from food and food processing environments allowed the detection of cases going back several years. This information, combined with the analysis of routinely collected enhanced surveillance data, helped to direct the investigation and identify the vehicle of infection. Retrospective whole genome sequencing (WGS) analysis of isolates provided robust microbiological evidence of links between cases, foods, and the environments in which they were produced and demonstrated that for some cases and foods, identified by fluorescent amplified fragment length polymorphism, the molecular typing method in routine use at the time, were not part of the outbreak. WGS analysis also showed that the strains causing illness had persisted in two food production environments for many years and in one producer had evolved into two strains over a period of around 8 years. This article demonstrates the value of reviewing L. monocytogenes typing data from clinical cases together with that from foods as a means of identifying potential vehicles and sources of infection in outbreaks of listeriosis. It illustrates the importance of reviewing retrospective L. monocytogenes typing alongside enhanced surveillance data to characterize extended outbreaks and inform control measures. Also, this article highlights the advantages of WGS analysis for strain discrimination and clarification of evolutionary relationships that refine outbreak investigations and improve our understanding of L. monocytogenes in the food chain.

Tracking of Listeria monocytogenes in meat establishment using Whole Genome Sequencing as a food safety management tool: A proof of concept.

Repeated Listeria outbreaks particularly associated with Ready-To-Eat (RTE) delicatessen meat products have been reported annually at global level. The most frequent scenario that led to foodborne outbreaks was the post-thermal treatment cross-contamination of deli meat products during slicing and modified atmosphere packaging (MAP). The precondition for such cross contamination is the previous introduction of Listeria into meat processing facilities and subsequent colonization of the production environment, associated with formation of biofilms resilient to common sanitation procedures regularly applied in meat establishments. The use of Whole Genome Sequencing (WGS) can facilitate the understanding of contamination and colonization routes of pathogens within the food production environment and enable efficient pathogen tracking among different departments. This study aimed to: a) provide a proof of concept on practical use of WGS in a meat establishment to define the entry routes and spread pattern of L. monocytogenes, and b) to consider the regular use of WGS in meat processing establishments as a strong support of food safety management system. The results revealed that Listeria spp. was present in slaughter line, chilling chambers, deboning, slicing, MAP, as well as in corridors and dispatch (53 positive samples, out of 240). Eight L. monocytogenes isolates (out of 53) were identified from the slaughterhouse, chilling chambers, deboning, MAP and dispatch. L. monocytogenes isolates were of three different serotypes (1/2a, 1/2c, 4b) and correspondingly of three MLST sequence types. Overall, two pairs of L. monocytogenes isolates were genetically identical, i.e. two serotype 4b isolates (ST1), isolated from water drain at dispatch unit and two isolates obtained from slaughterhouse (floorwall junction at the carcass wash point) and MAP (water drain). These findings indicated that L. monocytogenes isolates identified in meat processing units (MAP, chilling chamber and dispatch) originated from the slaughter line. Further, all eight L. monocytogenes isolates were confirmed to be biofilm producers on glass and stainless steel surfaces. The identification of the main entry routes of L. monocytogenes into meat establishments and tracking the routes for spread of the pathogen are of essential importance to define appropriate risk mitigation strategies for L. monocytogenes in meat production environment. The routine use of WGS for bacterial characterization, as a strong support of food safety management system in meat establishments, will require the cost-effective approach. It may encompass in-house sequencing when sequencing equipment is used for multiple applications (e.g. WGS of pathogens, starter cultures and spoilage organisms).

Exploiting Bacterial Whole-Genome Sequencing Data for Evaluation of Diagnostic Assays: Campylobacter Species Identification as a Case Study.

The application of whole-genome sequencing (WGS) to problems in clinical microbiology has had a major impact on the field. Clinical laboratories are now using WGS for pathogen identification, antimicrobial susceptibility testing, and epidemiological typing. WGS data also represent a valuable resource for the development and evaluation of molecular diagnostic assays, which continue to play an important role in clinical microbiology. To demonstrate this application of WGS, this study used publicly available genomic data to evaluate a duplex real-time PCR (RT-PCR) assay that targets mapA and ceuE for the detection of Campylobacter jejuni and Campylobacter coli, leading global causes of bacterial gastroenteritis. In silico analyses of mapA and ceuE primer and probe sequences from 1,713 genetically diverse C. jejuni and C. coli genomes, supported by RT-PCR testing, indicated that the assay was robust, with 1,707 (99.7%) isolates correctly identified. The high specificity of the mapA-ceuE assay was the result of interspecies diversity and intraspecies conservation of the target genes in C. jejuni and C. coli Rare instances of a lack of specificity among C. coli isolates were due to introgression in mapA or sequence diversity in ceuE The results of this study illustrate how WGS can be exploited to evaluate molecular diagnostic assays by using publicly available data, online databases, and open-source software.

MOST: a modified MLST typing tool based on short read sequencing.

Multilocus sequence typing (MLST) is an effective method to describe bacterial populations. Conventionally, MLST involves Polymerase Chain Reaction (PCR) amplification of housekeeping genes followed by Sanger DNA sequencing. Public Health England (PHE) is in the process of replacing the conventional MLST methodology with a method based on short read sequence data derived from Whole Genome Sequencing (WGS). This paper reports the comparison of the reliability of MLST results derived from WGS data, comparing mapping and assembly-based approaches to conventional methods using 323 bacterial genomes of diverse species. The sensitivity of the two WGS based methods were further investigated with 26 mixed and 29 low coverage genomic data sets from Salmonella enteridis and Streptococcus pneumoniae. Of the 323 samples, 92.9% (n = 300), 97.5% (n = 315) and 99.7% (n = 322) full MLST profiles were derived by the conventional method, assembly- and mapping-based approaches, respectively. The concordance between samples that were typed by conventional (92.9%) and both WGS methods was 100%. From the 55 mixed and low coverage genomes, 89.1% (n = 49) and 67.3% (n = 37) full MLST profiles were derived from the mapping and assembly based approaches, respectively. In conclusion, deriving MLST from WGS data is more sensitive than the conventional method. When comparing WGS based methods, the mapping based approach was the most sensitive. In addition, the mapping based approach described here derives quality metrics, which are difficult to determine quantitatively using conventional and WGS-assembly based approaches.

Utility of real-time amplification of selected 16S rRNA gene sequences as a tool for detection and identification of microbial signatures directly from clinical samples.

The potential of incorporating a real-time PCR for amplification and detection of 16S rRNA gene signatures directly from clinical samples was assessed as a tool for diagnostics. Universal PCR primers spanning short variable regions (~500 bp) were optimized for real-time PCR and tested in comparison with a longer fragment (~1400 bp) generated from block-based amplification. Real-time PCR had improved sensitivity of detection (8% increase), decreased amplification time and simplified downstream processing. The real-time PCR primers also offered an improvement in detection of bacteria from samples that demonstrated inhibition with the block-based primers and in the resolution of mixed-sequence traces. In addition to testing primer sensitivity, the effect of amplifying and sequencing two different variable regions of the 16S rRNA gene on organism identification was compared. By amplifying and sequencing a shorter variable region, the number of species that were identified to the species level was reduced by 18%.

Human papillomavirus genotype detection and viral load in paired genital and urine samples from both females and males.

The ability to detect type-specific high risk HPV (HR-HPV) infections in samples from females and males is important for monitoring the epidemiology of HPV and the impact of vaccination. Type-specific detection concordance between paired urine and genital samples from females (n = 264) undergoing routine colposcopy and males (n = 88) attending a genito-urinary medicine clinic was evaluated using an in-house genotyping assay. The overall inter-rater agreement (κ) was 0.781 for female pairs and 0.346 for male pairs. Female urine had sensitivity for detection of HPV16/18 and HR-HPV of 75% and 84%, respectively, while male urine had sensitivities of 13% and 28%, respectively. Genital samples had a higher HPV DNA copy number than urine although a small proportion (10%) of urine samples had a higher copy number than the corresponding genital sample. The proportion of females with normal cytology positive for HPV16/18 was 19%, increasing to 57% in moderate or severely dyskaryotic samples. The same trend was seen in the corresponding urine (19-43%) compounded by the reduced sensitivity of this sample type. The HPV16 viral load in female genital samples, but not in urine, was weakly associated with cervical disease stage. Despite reduced sensitivity, urine appears to be an appropriate surrogate sample for type-specific HPV detection in females for epidemiological objectives. The lower sensitivity and lack of association between viral load and disease stage in urine suggest that urine may not be useful for clinical management of HPV infection. The utility of urine for type-specific detection in males is less certain.

Amoxicillin therapy of poultry flocks: effect upon the selection of amoxicillin-resistant commensal Campylobacter spp.

BACKGROUND: The aim of this study was to investigate the effect of amoxicillin therapy of poultry flocks upon the persistence of commensal Campylobacter spp. and the incidence of antibiotic resistance. METHODS: Four poultry flocks naturally colonized with Campylobacter were treated with amoxicillin and monitored before, during and up to 4 weeks post-treatment. The numbers of Campylobacter were determined and the isolates speciated and typed by flaA short variable region (SVR) sequence analysis and PFGE. The susceptibility of the isolates to antibiotics, presence of the Cj0299 gene encoding a beta-lactamase and beta-lactamase production (nitrocefin hydrolysis) were also determined. RESULTS: Amoxicillin-resistant Campylobacter were isolated from Flock 1 before and during treatment, but Campylobacter were not detected afterwards. Flock 2 was colonized by amoxicillin-susceptible strains throughout sampling. No amoxicillin-resistant isolates arose during or after treatment. Flock 3 contained amoxicillin-susceptible and -resistant types pre-treatment. Resistant isolates were detected during treatment, while antibiotic-susceptible isolates re-emerged at 3 weeks post-treatment. All Campylobacter isolates from Flock 4 were amoxicillin resistant, irrespective of sampling time. All but one of the 82 amoxicillin-resistant (MICs 16 to >128 mg/L) Campylobacter jejuni and Campylobacter coli tested for the presence of Cj0299 carried the gene and all of these produced beta-lactamase. Co-amoxiclav remained active against amoxicillin-resistant isolates. CONCLUSIONS: Amoxicillin therapy had little effect on the numbers of amoxicillin-resistant commensal Campylobacter except for one flock where amoxicillin-resistant Campylobacter temporarily dominated. Amoxicillin therapy did not select amoxicillin-resistant isolates from a previous susceptible strain. Co-amoxiclav remained active against amoxicillin-resistant isolates.

Phase-variable surface structures are required for infection of Campylobacter jejuni by bacteriophages.

This study characterizes the interaction between Campylobacter jejuni and the 16 phages used in the United Kingdom typing scheme by screening spontaneous mutants of the phage-type strains and transposon mutants of the sequenced strain NCTC 11168. We show that the 16 typing phages fall into four groups based on their patterns of activity against spontaneous mutants. Screens of transposon and defined mutants indicate that the phage-bacterium interaction for one of these groups appears to involve the capsular polysaccharide (CPS), while two of the other three groups consist of flagellatropic phages. The expression of CPS and flagella is potentially phase variable in C. jejuni, and the implications of these findings for typing and intervention strategies are discussed.

Applicability of a rapid duplex real-time PCR assay for speciation of Campylobacter jejuni and Campylobacter coli directly from culture plates.

A rapid duplex real-time polymerase chain reaction (PCR) assay for speciation of Campylobacter jejuni and Campylobacter coli using the ABI Prism 7700 sequence detection system (Applied Biosystems) was developed based on two of the genes used in a conventional multiplex PCR. A rapid turnaround time of 3 h was achieved with the use of boiled cell lysates. Applicability of the assay was tested with 6015 random campylobacter strains referred to the Campylobacter Reference Unit, with 97.6% being identified as either C. jejuni or C. coli by this technique. Rapidity, combined with specificity and sensitivity, makes this method for routine campylobacter speciation attractive to any laboratory with a Taqman system.