Salmonella nomenclature in the genomic era: a time for change

Salmonella enterica nomenclature has evolved over the past one hundred years into a highly sophisticated naming convention based on the recognition of antigens by specific antibodies. This serotyping scheme has led to the definition of over 2500 serovars which are well understood, have standing in nomenclature and, for the majority, biological relevance. Therefore, it is highly desirable for any change in naming convention to maintain backwards compatibility with the information linked to these serovars. The routine use of whole genome sequencing and the well-established link between sequence types and serovars presents an opportunity to update the scheme by incorporating the phylogenetically relevant sequence data whilst preserving the best of serotyping nomenclature. Advantages include: overcoming the variability in antibody preparations; removing the need to use laboratory animals and implementing a truly universal system. However, the issue of trying to reproduce the phenotyping gold standard needs to be relaxed if we are to fully embrace the genomic era. We have used whole genome sequence data from over 46,000 isolates of Salmonella enterica subspecies enterica to define clusters in two stages: Multi Locus Sequence Typing followed by antigen prediction. Sequence type—serotype discrepancies were resolved using core SNP clustering to determine the phylogenetic groups and this was confirmed by overlaying the antigenic prediction onto the core SNP clusters and testing the separation of clusters using cgMLST Hierarchical Clustering. This allowed us to define any major antigenic clusters within an ST—here called the MAC type and written as ST-serovar. Using this method, 99.96% of Salmonella isolates reported in the UK were assigned a MAC type and linked to a serovar name taken from the Kauffmann and White scheme. We propose a change for reporting of Salmonella enterica sub-types using the ST followed by serovar

Identification and characterisation of enteroaggregative Escherichia coli subtypes associated with human disease

Enteroaggregative E. coli (EAEC) are a major cause of diarrhoea worldwide. Due to their heterogeneity and carriage in healthy individuals, identification of diagnostic virulence markers for pathogenic strains has been difficult. In this study, we have determined phenotypic and genotypic differences between EAEC strains of sequence types (STs) epidemiologically associated with asymptomatic carriage (ST31) and diarrhoeal disease (ST40). ST40 strains demonstrated significantly enhanced intestinal adherence, biofilm formation, and pro-inflammatory interleukin-8 secretion compared with ST31 isolates. This was independent of whether strains were derived from diarrhoea patients or healthy controls. Whole genome sequencing revealed differences in putative virulence genes encoding aggregative adherence fimbriae, E. coli common pilus, flagellin and EAEC heat-stable enterotoxin 1. Our results indicate that ST40 strains have a higher intrinsic potential of human pathogenesis due to a specific combination of virulence-related factors which promote host cell colonization and inflammation. These findings may contribute to the development of genotypic and/or phenotypic markers for EAEC strains of high virulence

Use of whole-genus genome sequence data to develop a multilocus sequence typing tool that accurately identifies Yersinia isolates to the species and subspecies levels

The genus Yersinia is a large and diverse bacterial genus consisting of human-pathogenic species, a fish-pathogenic species, and a large number of environmental species. Recently, the phylogenetic and population structure of the entire genus was elucidated through the genome sequence data of 241 strains encompassing every known species in the genus. Here we report the mining of this enormous data set to create a multilocus sequence typing-based scheme that can identify Yersinia strains to the species level to a level of resolution equal to that for whole-genome sequencing. Our assay is designed to be able to accurately subtype the important human-pathogenic species Yersinia enterocolitica to whole-genome resolution levels. We also report the validation of the scheme on 386 strains from reference laboratory collections across Europe. We propose that the scheme is an important molecular typing system to allow accurate and reproducible identification of Yersinia isolates to the species level, a process often inconsistent in nonspecialist laboratories. Additionally, our assay is the most phylogenetically informative typing scheme available for Y. enterocolitica

Whole Genome Sequencing for Public Health Surveillance of Shiga Toxin-Producing Escherichia coli Other than Serogroup O157

Shiga toxin-producing Escherichia coli (STEC) are considered to be a significant threat to public health due to the severity of gastrointestinal symptoms associated with human infection. In England STEC O157 is the most commonly detected STEC serogroup, however, the implementation of PCR at local hospital laboartories has resulted in an increase in the detection of STEC other than serogroup O157 (non-O157 STEC). The aim of this study was to evaluate the use of whole genome sequencing (WGS) for routine public health surveillance of non-O157 STEC by comparing this approach to phenotypic serotyping and PCR for subtyping the stx-encoding genes. Of the 102 isolates where phenotypic and genotypic serotyping could be compared, 98 gave fully concordant results. The most common non-O157 STEC serogroups detected were O146 (22) and O26 (18). All but one of the 38 isolates that could not be phenotypically serotyped (designated O unidentifiable or O rough) were serotyped using the WGS data. Of the 73 isolates where a flagella type was available by traditional phenotypic typing, all results matched the H-type derived from the WGS data. Of the 140 sequenced non-O157 isolates, 52 (37.1%) harboured stx1 only, 42 (30.0%) had stx2 only, 46 (32.9%) carried stx1 and stx2. Of these, stx subtyping PCR results were available for 131 isolates and 121 of these had concordant results with the stx subtype derived from the WGS data. Non-specific primer binding during PCR amplification, due to the similarity of the stx2 subtype gene sequences was the most likely cause. The results of this study showed WGS provided a reliable and robust one-step process for characterisation of STEC. Deriving the full serotype from WGS data in real time has enabled us to report a higher level of strain discrimination while stx subtyping provides data on the pathogenic potential of each isolate, enabling us to predict clinical outcome of each case and to monitor the emergence of hyper-virulent strains

The origins of haplotype 58 (H58) Salmonella enterica serovar Typhi

Antimicrobial resistance (AMR) poses a serious threat to the clinical management of typhoid fever. AMR in Salmonella Typhi (S. Typhi) is commonly associated with the H58 lineage, a lineage that arose comparatively recently before becoming globally disseminated. To better understand when and how H58 emerged and became dominant, we performed detailed phylogenetic analyses on contemporary genome sequences from S. Typhi isolated in the period spanning the emergence. Our dataset, which contains the earliest described H58 S. Typhi organism, indicates that ancestral H58 organisms were already multi-drug resistant (MDR). These organisms emerged spontaneously in India in 1987 and became radially distributed throughout South Asia and then globally in the ensuing years. These early organisms were associated with a single long branch, possessing mutations associated with increased bile tolerance, suggesting that the first H58 organism was generated during chronic carriage. The subsequent use of fluoroquinolones led to several independent mutations in gyrA. The ability of H58 to acquire and maintain AMR genes continues to pose a threat, as extensively drug-resistant (XDR; MDR plus resistance to ciprofloxacin and third generation cephalosporins) variants, have emerged recently in this lineage. Understanding where and how H58 S. Typhi originated and became successful is key to understand how AMR drives successful lineages of bacterial pathogens. Additionally, these data can inform optimal targeting of typhoid conjugate vaccines (TCVs) for reducing the potential for emergence and the impact of new drug-resistant variants. Emphasis should also be placed upon the prospective identification and treatment of chronic carriers to prevent the emergence of new drug resistant variants with the ability to spread efficiently

Rapid geographical source attribution of Salmonella enterica serovar Enteritidis genomes using hierarchical machine learning.

Salmonella enterica serovar Enteritidis is one of the most frequent causes of Salmonellosis globally and is commonly transmitted from animals to humans by the consumption of contaminated foodstuffs. In the UK and many other countries in the Global North, a significant proportion of cases are caused by consumption of imported food products or contracted during foreign travel, therefore making the rapid identification of the geographical source of new infections a requirement for robust public health outbreak investigations. Herein, we detail the development and application of a hierarchical machine learning model to rapidly identify and trace the geographical source of S. Enteritidis infections from whole genome sequencing data. 2,313 S. Enteritidis genomes, collected by the UKHSA between 2014-2019, were used to train a 'local classifier per node' hierarchical classifier to attribute isolates to 4 continents, 11 sub-regions and 38 countries (53 classes). The highest classification accuracy was achieved at the continental level followed by the sub-regional and country levels (macro F1: 0.954, 0.718, 0.661 respectively). A number of countries commonly visited by UK travellers were predicted with high accuracy (hF1: >0.9). Longitudinal analysis and validation with publicly accessible international samples indicated that predictions were robust to prospective external datasets. The hierarchical machine learning framework provided granular geographical source prediction directly from sequencing reads in <4 minutes per sample, facilitating rapid outbreak resolution and real-time genomic epidemiology. The results suggest additional application to a broader range of pathogens and other geographically structured problems, such as antimicrobial resistance prediction, is warranted

Increased phage resistance through lysogenic conversion accompanying emergence of monophasic Salmonella Typhimurium ST34 pandemic strain

Salmonella enterica serovar Typhimurium (S. Typhimurium) comprises a group of closely related human and animal pathogens that account for a large proportion of all Salmonella infections globally. The epidemiological record of S. Typhimurium in Europe is characterized by successive waves of dominant clones, each prevailing for approximately 10–15 years before replacement. Succession of epidemic clones may represent a moving target for interventions aimed at controlling the spread and impact of this pathogen on human and animal health. Here, we investigate the relationship of phage sensitivity and population structure of S. Typhimurium using data from the Anderson phage typing scheme. We observed greater resistance to phage predation of epidemic clones circulating in livestock over the past decades compared to variants with a restricted host range implicating increased resistance to phage in the emergence of epidemic clones of particular importance to human health. Emergence of monophasic S. Typhimurium ST34, the most recent dominant multidrug-resistant clone, was accompanied by increased resistance to phage predation during clonal expansion, in part by the acquisition of the mTmII prophage that may have contributed to the fitness of the strains that replaced ancestors lacking this prophage

Increased phage resistance through lysogenic conversion accompanying emergence of monophasic Salmonella Typhimurium ST34 pandemic strain

Salmonella enterica serovar Typhimurium ( S. Typhimurium) comprises a group of closely related human and animal pathogens that account for a large proportion of all Salmonella infections globally. The epidemiological record of S. Typhimurium in Europe is characterized by successive waves of dominant clones, each prevailing for approximately 10-15 years before replacement. Succession of epidemic clones may represent a moving target for interventions aimed at controlling the spread and impact of this pathogen on human and animal health. Here, we investigate the relationship of phage sensitivity and population structure of S. Typhimurium using data from the Anderson phage typing scheme. We observed greater resistance to phage predation of epidemic clones circulating in livestock over the past decades compared to variants with a restricted host range implicating increased resistance to phage in the emergence of epidemic clones of particular importance to human health. Emergence of monophasic S. Typhimurium ST34, the most recent dominant multidrug-resistant clone, was accompanied by increased resistance to phage predation during clonal expansion, in part by the acquisition of the mTmII prophage that may have contributed to the fitness of the strains that replaced ancestors lacking this prophage

Characterization of a pESI-like plasmid and analysis of multidrug-resistant Salmonella enterica Infantis isolates in England and Wales

Salmonella enterica serovar Infantis is the fifth most common Salmonella serovar isolated in England and Wales. Epidemiological, genotyping and antimicrobial-resistance data for S. enterica Infantis isolates were used to analyse English and Welsh demographics over a 5 year period. Travel cases associated with S. enterica Infantis were mainly from Asia, followed by cases from Europe and North America. Since 2000, increasing numbers of S. enterica Infantis had multidrug resistance determinants harboured on a large plasmid termed ‘plasmid of emerging S. enterica Infantis’ (pESI). Between 2013 and 2018, 42 S. enterica Infantis isolates were isolated from humans and food that harboured resistance determinants to multiple antimicrobial classes present on a pESI-like plasmid, including extended-spectrum β-lactamases (ESBLs; blaCTX-M-65). Nanopore sequencing of an ESBL-producing human S. enterica Infantis isolate indicated the presence of two regions on an IncFIB pESI-like plasmid harbouring multiple resistance genes. Phylogenetic analysis of the English and Welsh S. enterica Infantis population indicated that the majority of multidrug-resistant isolates harbouring the pESI-like plasmid belonged to a single clade maintained within the population. The blaCTX-M-65 ESBL isolates first isolated in 2013 comprise a lineage within this clade, which was mainly associated with South America. Our data, therefore, show the emergence of a stable resistant clone that has been in circulation for some time in the human population in England and Wales, highlighting the necessity of monitoring resistance in this serovar

Enteroaggregative escherichia coli have evolved independently as distinct complexes within the E. Coli population with varying ability to cause disease

Enteroaggregative E. Coli (EAEC) is an established diarrhoeagenic pathotype. The association with virulence gene content and ability to cause disease has been studied but little is known about the population structure of EAEC and how this pathotype evolved. Analysis by Multi Locus Sequence Typing of 564 EAEC isolates from cases and controls in Bangladesh, Nigeria and the UK spanning the past 29 years, revealed multiple successful lineages of EAEC. The population structure of EAEC indicates some clusters are statistically associated with disease or carriage, further highlighting the heterogeneous nature of this group of organisms. Different clusters have evolved independently as a result of both mutational and recombination events; the EAEC phenotype is distributed throughout the population of E. coli