Evidence of microevolution of Salmonella Typhimurium during a series of egg-associated outbreaks linked to a single chicken farm

Genomic insights of Salmonella isolated from dry fermented sausage production chains in Spain and France

The presence of Salmonella in dry fermented sausages is source of recalls and outbreaks. The genomic diversity of 173 Salmonella isolates from the dry fermented sausage production chains (pig carcasses, pork, and sausages) from France and Spain were investigated through their core phylogenomic relationships and accessory genome profiles. Ten different serovars and thirteen sequence type profiles were identified. The most frequent serovar from sausages was the monophasic variant of S. Typhimurium (1,4,[5],12:i:-, 72%) while S. Derby was in pig carcasses (51%). Phylogenomic clusters found in S. 1,4,[5],12:i:-, S. Derby, S. Rissen and S. Typhimurium serovars identified closely related isolates, with less than 10 alleles and 20 SNPs of difference, displaying Salmonella persistence along the pork production chain. Most of the S. 1,4,[5],12:i:- contained the Salmonella genomic island-4 (SGI-4), Tn21 and IncFIB plasmid. More than half of S. Derby strains contained the SGI-1 and Tn7. S. 1,4,[5],12:i:- genomes carried the most multidrug resistance genes (91% of the strains), whereas extended-spectrum β-lactamase genes were found in Typhimurium and Derby serovars. Salmonella monitoring and characterization in the pork production chains, specially S. 1,4,[5],12:i:- serovar, is of special importance due to its multidrug resistance capacity and persistence in dry fermented sausages.

Epidemiological links and antimicrobial resistance of clinical Salmonella enterica ST198 isolates: a nationwide microbial population genomic study in Switzerland

Salmonella is a leading cause of foodborne outbreaks and systemic infections worldwide. Emerging multi-drug resistant Salmonella lineages such as a ciprofloxacin-resistant subclade (CIP^R) within Salmonella enterica serovar Kentucky ST198 threaten the effective prevention and treatment of infections. To understand the genomic diversity and antimicrobial resistance gene content associated with S. Kentucky in Switzerland, we whole-genome sequenced 70 human clinical isolates obtained between 2010 and 2020. Most isolates belonged to ST198-CIP^R. High- and low-level ciprofloxacin resistance among CIP^R isolates was associated with variable mutations in ramR and acrB in combination with stable mutations in quinolone-resistance determining regions (QRDRs). Analysis of isolates from patients with prolonged ST198 colonization indicated subclonal adaptions with the ramR locus as a mutational hotspot. SNP analyses identified multiple clusters of near-identical isolates, which were often associated with travel but included spatiotemporally linked isolates from Switzerland. The largest SNP cluster was associated with travellers returning from Indonesia, and investigation of global data linked >60 additional ST198 salmonellosis isolates to this cluster. Our results emphasize the urgent need for implementing whole-genome sequencing as a routine tool for Salmonella surveillance and outbreak detection.

Salmonella Genomics in Public Health and Food Safety

The species Salmonella enterica comprises over 2,600 serovars, many of which are known to be intracellular pathogens of mammals, birds, and reptiles. It is now apparent that Salmonella is a highly adapted environmental microbe and can readily persist in a number of environmental niches, including water, soil, and various plant (including produce) species. Much of what is known about the evolution and diversity of nontyphoidal Salmonella serovars (NTS) in the environment is the result of the rise of the genomics era in enteric microbiology. There are over 340,000 Salmonella genomes available in public databases. This extraordinary breadth of genomic diversity now available for the species, coupled with widespread availability and affordability of whole-genome sequencing (WGS) instrumentation, has transformed the way in which we detect, differentiate, and characterize Salmonella enterica strains in a timely way. Not only have WGS data afforded a detailed and global examination of the molecular epidemiological movement of Salmonella from diverse environmental reservoirs into human and animal hosts, but they have also allowed considerable consolidation of the diagnostic effort required to test for various phenotypes important to the characterization of Salmonella. For example, drug resistance, serovar, virulence determinants, and other genome-based attributes can all be discerned using a genome sequence. Finally, genomic analysis, in conjunction with functional and phenotypic approaches, is beginning to provide new insights into the precise adaptive changes that permit persistence of NTS in so many diverse and challenging environmental niches.

Salmonella enterica from a soldier from the 1652 siege of Barcelona (Spain) supports historical transatlantic epidemic contacts

Ancient pathogen genomics is an emerging field allowing reconstruction of past epidemics. The demise of post-contact American populations may, at least in part, have been caused by paratyphoid fever brought by Europeans. We retrieved genome-wide data from two Spanish soldiers who were besieging the city of Barcelona in 1652, during the Reapers' War. Their ancestry derived from the Basque region and Sardinia, respectively, (at that time, this island belonged to the Spanish kingdom). Despite the proposed plague epidemic, we could not find solid evidence for the presence of the causative plague agent in these individuals. However, we retrieved from one individual a substantial fraction of the Salmonella enterica serovar Paratyphi C lineage linked to paratyphoid fever in colonial period Mexico. Our results support a growing body of evidence that Paratyphi C enteric fever was more prevalent in Europe and the Americas in the past than it is today.

Characterisation of IncI1 plasmids associated with change of phage type in isolates of Salmonella enterica serovar Typhimurium

Background

Acquisition of IncI1 plasmids by members of the Enterobacteriaceae sometimes leads to transfer of antimicrobial resistance and colicinogeny as well as change of phage type in Salmonella Typhimurium. Isolates of S. Typhimurium from a 2015 outbreak of food poisoning were found to contain an IncI1 plasmid implicated in change of phage type from PT135a to U307 not previously reported. The origin of the changes of phage type associated with this IncI1 plasmid was investigated. In addition, a comparison of its gene composition with that of IncI1 plasmids found in local isolates of S. Typhimurium typed as U307 from other times was undertaken. This comparison was extended to IncI1 plasmids in isolates of phage types PT6 and PT6 var. 1 which are thought to be associated with acquisition of IncI1 plasmids.

Results

Analysis of IncI1 plasmids from whole genome sequencing of isolates implicated a gene coding for a 1273 amino acid protein present only in U307 isolates as the likely source of change of phage type. The IncI1 plasmids from PT6 and PT6 var. 1 isolates all had the ibfA gene present in IncI1 plasmid R64. This gene inhibits growth of bacteriophage BF23 and was therefore the possible source of change of phage type. A fuller comparison of the genetic composition of IncI1 plasmids from U307 isolates and PT6 and PT6 var. 1 isolates along with two IncI1 plasmids from S. Typhimurium isolates not showing change of phage type was undertaken. Plasmids were classified as either ‘Delta’ or ‘Col’ IncI1 plasmids according to whether genes between repZ and the rfsF site showed high identity to genes in the same location in R64 or ColIb-P9 plasmids respectively. Comparison of the tra gene sets and the pil gene sets across the range of sequenced plasmids identified Delta and Col plasmids with almost identical sequences for both sets of genes. This indicated a genetic recombination event leading to a switch between Delta and Col gene sets at the rfsF site. Comparisons of other gene sets showing significant variation among the sequenced plasmids are reported. Searches of the NCBI GenBank database using DNA and protein sequences of interest from the sequenced plasmids identified global IncI1 plasmids with extensive regions showing 99 to 100% identity to some of the plasmids sequenced in this study indicating evidence for widespread distribution of these plasmids.

Conclusion

Two genes possibly associated with change of phage type were identified in IncI1 plasmids. IncI1 plasmids were classified as either ‘Delta’ or ‘Col’ plasmids and other sequences of significant variation among these plasmids were identified. This study offers a new perspective on the understanding of the gene composition of IncI1 plasmids. The sequences of newly sequenced IncI1 plasmids could be compared against the regions of significant sequence variation identified in this study to understand better their overall gene composition and relatedness to other IncI1 plasmids in the databases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02151-z.

Toward an Integrated Genome-Based Surveillance of Salmonella enterica in Germany

Despite extensive monitoring programs and preventative measures, Salmonella spp. continue to cause tens of thousands human infections per year, as well as many regional and international food-borne outbreaks, that are of great importance for public health and cause significant socio-economic costs. In Germany, salmonellosis is the second most common cause of bacterial diarrhea in humans and is associated with high hospitalization rates. Whole-genome sequencing (WGS) combined with data analysis is a high throughput technology with an unprecedented discriminatory power, which is particularly well suited for targeted pathogen monitoring, rapid cluster detection and assignment of possible infection sources. However, an effective implementation of WGS methods for large-scale microbial pathogen detection and surveillance has been hampered by the lack of standardized methods, uniform quality criteria and strategies for data sharing, all of which are essential for a successful interpretation of sequencing data from different sources. To overcome these challenges, the national GenoSalmSurv project aims to establish a working model for an integrated genome-based surveillance system of Salmonella spp. in Germany, based on a decentralized data analysis. Backbone of the model is the harmonization of laboratory procedures and sequencing protocols, the implementation of open-source bioinformatics tools for data analysis at each institution and the establishment of routine practices for cross-sectoral data sharing for a uniform result interpretation. With this model, we present a working solution for cross-sector interpretation of sequencing data from different sources (such as human, veterinarian, food, feed and environmental) and outline how a decentralized data analysis can contribute to a uniform cluster detection and facilitate outbreak investigations.

Escherichia coli ST8196 is a novel, locally evolved, and extensively drug resistant pathogenic lineage within the ST131 clonal complex

The H30Rx subclade of Escherichia coli ST131 is a clinically important, globally dispersed pathogenic lineage that typically displays resistance to fluoroquinolones and extended spectrum β-lactams. Isolates EC233 and EC234, variants of ST131-H30Rx with a novel sequence type (ST) 8196, isolated from unrelated patients presenting with bacteraemia at a Sydney Hospital in 2014 are characterised here. EC233 and EC234 are phylogroup B2, serotype O25:H4A, and resistant to ampicillin, amoxicillin, cefoxitin, ceftazidime, ceftriaxone, ciprofloxacin, norfloxacin and gentamicin and are likely clonal. Both harbour an IncFII_2 plasmid (pSPRC_Ec234-FII) that carries most of the resistance genes on an IS26 associated translocatable unit, two small plasmids and a novel IncI1 plasmid (pSPRC_Ec234-I). SNP-based phylogenetic analysis of the core genome of representatives within the ST131 clonal complex places both isolates in a subclade with three clinical Australian ST131-H30Rx clade-C isolates. A MrBayes phylogeny analysis of EC233 and EC234 indicates ST8196 share a most recent common ancestor with ST131-H30Rx strain EC70 isolated from the same hospital in 2013. Our study identified genomic hallmarks that define the ST131-H30Rx subclade in the ST8196 isolates and highlights a need for unbiased genomic surveillance approaches to identify novel high-risk MDR E. coli pathogens that impact healthcare facilities.

Evolution of Salmonella enterica serotype Typhimurium driven by anthropogenic selection and niche adaptation

Salmonella enterica serotype Typhimurium (S. Typhimurium) is a leading cause of gastroenteritis and bacteraemia worldwide, and a model organism for the study of host-pathogen interactions. Two S. Typhimurium strains (SL1344 and ATCC14028) are widely used to study host-pathogen interactions, yet genotypic variation results in strains with diverse host range, pathogenicity and risk to food safety. The population structure of diverse strains of S. Typhimurium revealed a major phylogroup of predominantly sequence type 19 (ST19) and a minor phylogroup of ST36. The major phylogroup had a population structure with two high order clades (α and β) and multiple subclades on extended internal branches, that exhibited distinct signatures of host adaptation and anthropogenic selection. Clade α contained a number of subclades composed of strains from well characterized epidemics in domesticated animals, while clade β contained multiple subclades associated with wild avian species. The contrasting epidemiology of strains in clade α and β was reflected by the distinct distribution of antimicrobial resistance (AMR) genes, accumulation of hypothetically disrupted coding sequences (HDCS), and signatures of functional diversification. These observations were consistent with elevated anthropogenic selection of clade α lineages from adaptation to circulation in populations of domesticated livestock, and the predisposition of clade β lineages to undergo adaptation to an invasive lifestyle by a process of convergent evolution with of host adapted Salmonella serotypes. Gene flux was predominantly driven by acquisition and recombination of prophage and associated cargo genes, with only occasional loss of these elements. The acquisition of large chromosomally-encoded genetic islands was limited, but notably, a feature of two recent pandemic clones (DT104 and monophasic S. Typhimurium ST34) of clade α (SGI-1 and SGI-4).

Phylogenetic analysis revealed that Salmonella Typhimurium ST313 isolated from humans and food in Brazil presented a high genomic similarity

Salmonella Typhimurium sequence type 313 (S. Typhimurium ST313) has caused invasive disease mainly in sub-Saharan Africa. In Brazil, ST313 strains have been recently described, and there is a lack of studies that assessed by whole genome sequencing (WGS)-the relationship of these strains. The aims of this work were to study the phylogenetic relationship of 70 S. Typhimurium genomes comparing strains of ST313 (n = 9) isolated from humans and food in Brazil among themselves, with other STs isolated in this country (n = 31) and in other parts of the globe (n = 30) by 16S rRNA sequences, the Gegenees software, whole genome multilocus sequence typing (wgMLST), and average nucleotide identity (ANI) for the genomes of ST313. Additionally, pangenome analysis was performed to verify the heterogeneity of these genomes. The phylogenetic analyses showed that the ST313 genomes were very similar among themselves. However, the ST313 genomes were usually clustered more distantly to other STs of strains isolated in Brazil and in other parts of the world. By pangenome calculation, the core genome was 2,880 CDSs and 4,171 CDSs singletons for all the 70 S. Typhimurium genomes studied. Considering the 10 ST313 genomes analyzed the core genome was 4,112 CDSs and 76 CDSs singletons. In conclusion, the ST313 genomes from Brazil showed a high similarity among them which information might eventually help in the development of vaccines and antibiotics. The pangenome analysis showed that the S. Typhimurium genomes studied presented an open pangenome, but specifically tending to become close for the ST313 strains.

A Simple and Robust Statistical Method to Define Genetic Relatedness of Samples Related to Outbreaks at the Genomic Scale - Application to Retrospective Salmonella Foodborne Outbreak Investigations

The investigation of foodborne outbreaks (FBOs) from genomic data typically relies on inspecting the relatedness of samples through a phylogenomic tree computed on either SNPs, genes, kmers, or alleles (i.e., cgMLST and wgMLST). The phylogenomic reconstruction is often time-consuming, computation-intensive and depends on hidden assumptions, pipelines implementation and their parameterization. In the context of FBO investigations, robust links between isolates are required in a timely manner to trigger appropriate management actions. Here, we propose a non-parametric statistical method to assert the relatedness of samples (i.e., outbreak cases) or whether to reject them (i.e., non-outbreak cases). With typical computation running within minutes on a desktop computer, we benchmarked the ability of three non-parametric statistical tests (i.e., Wilcoxon rank-sum, Kolmogorov-Smirnov and Kruskal-Wallis) on six different genomic features (i.e., SNPs, SNPs excluding recombination events, genes, kmers, cgMLST alleles, and wgMLST alleles) to discriminate outbreak cases (i.e., positive control: C+) from non-outbreak cases (i.e., negative control: C-). We leveraged four well-characterized and retrospectively investigated FBOs of Salmonella Typhimurium and its monophasic variant S. 1,4,[5],12:i:- from France, setting positive and negative controls in all the assays. We show that the approaches relying on pairwise SNP differences distinguished all four considered outbreaks in contrast to the other tested genomic features (i.e., genes, kmers, cgMLST alleles, and wgMLST alleles). The freely available non-parametric method written in R has been designed to be independent of both the phylogenomic reconstruction and the detection methods of genomic features (i.e., SNPs, genes, kmers, or alleles), making it widely and easily usable to anybody working on genomic data from suspected samples.

Differentiating Epidemic from Endemic or Sporadic Infectious Disease Occurrence

One important scope of work of epidemiology is the investigation of infectious diseases that cluster in time and place. Clusters of infectious disease may represent outbreaks or epidemics in which the cases share in common a point source exposure or an infectious agent in a chain of transmission pathways. Investigations of outbreaks of an illness can facilitate identification of a source, risk, or cause of the illness. However, most infectious disease episodes occur not as part of any apparent outbreaks but as sporadic infections. Multiple sporadic infections that occur steadily in time and place are referred to as endemic disease. How does one investigate sources and risk factors for sporadic or endemic infections? As part of the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology of Infectious Diseases, this review discusses limitations of traditional approaches and advantages of molecular epidemiology approaches to investigate sporadic and endemic infections. Using specific examples, the discussions show that most sporadic infections are actually part of unrecognized outbreaks and that what appears to be endemic disease occurrence is actually comprised of multiple small outbreaks. These molecular epidemiologic investigations have unmasked modes of transmission of infectious agents not known to cause outbreaks. They have also raised questions about the traditional ways to measure incidence and assess sources of drug-resistant infections in community settings. The discoveries made by the application of molecular microbiology methods in epidemiologic investigations have led to creation of new public health intervention strategies that have not been previously considered. *This article is part of a curated collection.

Serotype-specific evolutionary patterns of antimicrobial-resistant Salmonella enterica

Background

The emergence of antimicrobial-resistant (AMR) strains of the important human and animal pathogen Salmonella enterica poses a growing threat to public health. Here, we studied the genome-wide evolution of 90 S. enterica AMR isolates, representing one host adapted serotype (S. Dublin) and two broad host range serotypes (S. Newport and S. Typhimurium).

Results

AMR S. Typhimurium had a large effective population size, a large and diverse genome, AMR profiles with high diversity, and frequent positive selection and homologous recombination. AMR S. Newport showed a relatively low level of diversity and a relatively clonal population structure. AMR S. Dublin showed evidence for a recent population bottleneck, and the genomes were characterized by a larger number of genes and gene ontology terms specifically absent from this serotype and a significantly higher number of pseudogenes as compared to other two serotypes. Approximately 50% of accessory genes, including specific AMR and putative prophage genes, were significantly over- or under-represented in a given serotype. Approximately 65% of the core genes showed phylogenetic clustering by serotype, including the AMR gene aac (6′)-Iaa. While cell surface proteins were shown to be the main target of positive selection, some proteins with possible functions in AMR and virulence also showed evidence for positive selection. Homologous recombination mainly acted on prophage-associated proteins.

Conclusions

Our data indicates a strong association between genome content of S. enterica and serotype. Evolutionary patterns observed in S. Typhimurium are consistent with multiple emergence events of AMR strains and/or ecological success of this serotype in different hosts or habitats. Evolutionary patterns of S. Newport suggested that antimicrobial resistance emerged in one single lineage, Lineage IIC. A recent population bottleneck and genome decay observed in AMR S. Dublin are congruent with its narrow host range. Finally, our results suggest the potentially important role of positive selection in the evolution of antimicrobial resistance, host adaptation and serotype diversification in S. enterica.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1457-5) contains supplementary material, which is available to authorized users.

Microevolution of antimicrobial resistance and biofilm formation of Salmonella Typhimurium during persistence on pig farms

Salmonella Typhimurium and its monophasic variant S. 4,[5],12:i:- are the dominant serotypes associated with pigs in many countries. We investigated their population structure on nine farms using whole genome sequencing, and their genotypic and phenotypic variation. The population structure revealed the presence of phylogenetically distinct clades consisting of closely related clones of S. Typhimurium or S. 4,[5],12:i:- on each pig farm, that persisted between production cycles. All the S. 4,[5],12:i:- strains carried the Salmonella genomic island-4 (SGI-4), which confers resistance to heavy metals, and half of the strains contained the mTmV prophage, harbouring the sopE virulence gene. Most clonal groups were highly drug resistant due to the presence of multiple antimicrobial resistance (AMR) genes, and two clades exhibited evidence of recent on-farm plasmid-mediated acquisition of additional AMR genes, including an IncHI2 plasmid. Biofilm formation was highly variable but had a strong phylogenetic signature. Strains capable of forming biofilm with the greatest biomass were from the S. 4,[5],12:i:- and S. Typhimurium DT104 clades, the two dominant pandemic clones found over the last 25 years. On-farm microevolution resulted in enhanced biofilm formation in subsequent production cycle.

Detection of VIM-1-Producing Enterobacter cloacae and Salmonella enterica Serovars Infantis and Goldcoast at a Breeding Pig Farm in Germany in 2017 and Their Molecular Relationship to Former VIM-1-Producing S. Infantis Isolates in German Livestock Production

In 2011, VIM-1-producing Salmonella enterica serovar Infantis and Escherichia coli were isolated for the first time in four German livestock farms. In 2015/2016, highly related isolates were identified in German pig production. This raised the issue of potential reservoirs for these isolates, the relation of their mobile genetic elements, and potential links between the different affected farms/facilities. In a piglet-producing farm suspicious for being linked to some blaVIM-1 findings in Germany, fecal and environmental samples were examined for the presence of carbapenemase-producing Enterobacteriaceae and Salmonella spp. Newly discovered isolates were subjected to Illumina whole-genome sequencing (WGS) and S1 pulsed-field gel electrophoresis (PFGE) hybridization experiments. WGS data of these isolates were compared with those for the previously isolated VIM-1-producing Salmonella Infantis isolates from pigs and poultry. Among 103 samples, one Salmonella Goldcoast isolate, one Salmonella Infantis isolate, and one Enterobacter cloacae isolate carrying the blaVIM-1 gene were detected. Comparative WGS analysis revealed that the blaVIM-1 gene was part of a particular Tn21-like transposable element in all isolates. It was located on IncHI2 (ST1) plasmids of ∼290 to 300 kb with a backbone highly similar (98 to 100%) to that of reference pSE15-SA01028. SNP analysis revealed a close relationship of all VIM-1-positive S Infantis isolates described since 2011. The findings of this study demonstrate that the occurrence of the blaVIM-1 gene in German livestock is restricted neither to a certain bacterial species nor to a certain Salmonella serovar but is linked to a particular Tn21-like transposable element located on transferable pSE15-SA01028-like IncHI2 (ST1) plasmids, being present in all of the investigated isolates from 2011 to 2017.IMPORTANCE Carbapenems are considered one of few remaining treatment options against multidrug-resistant Gram-negative pathogens in human clinical settings. The occurrence of carbapenemase-producing Enterobacteriaceae in livestock and food is a major public health concern. Particularly the occurrence of VIM-1-producing Salmonella Infantis in livestock farms is worrisome, as this zoonotic pathogen is one of the main causes for human salmonellosis in Europe. Investigations on the epidemiology of those carbapenemase-producing isolates and associated mobile genetic elements through an in-depth molecular characterization are indispensable to understand the transmission of carbapenemase-producing Enterobacteriaceae along the food chain and between different populations to develop strategies to prevent their further spread.

Cross-Border Transmission of Salmonella Choleraesuis var. Kunzendorf in European Pigs and Wild Boar: Infection, Genetics, and Evolution

Salmonella enterica subspecies enterica serotype Choleraesuis is a swine adapted serovar. S. Choleraesuis variant Kunzendorf is responsible for the majority of outbreaks among pigs. S. Choleraesuis is rare in Europe, although there have been serious outbreaks in pigs including two outbreaks in Denmark in 1999–2000 and 2012–2013. Here, we elucidate the epidemiology, possible transmission routes and sources, and clonality of European S. Choleraesuis isolates including the Danish outbreak isolates. A total of 102 S. Choleraesuis isolates from different European countries and the United States, covering available isolates from the last two decades were selected for whole genome sequencing. We applied a temporally structured sequence analysis within a Bayesian framework to reconstruct a temporal and spatial phylogenetic tree. MLST type, resistance genes, plasmid replicons, and accessory genes were identified using bioinformatics tools. Fifty-eight isolates including 11 out of 12 strains from wild boars were pan-susceptible. The remaining isolates carried multiple resistance genes. Eleven different plasmid replicons in eight plasmids were determined among the isolates. Accessory genes were associated to the identified resistance genes and plasmids. The European S. Choleraesuis was estimated to have emerged in ∼1837 (95% credible interval, 1733–1983) with the mutation rate of 1.02 SNPs/genome/year. The isolates were clustered according to countries and neighbor countries. There were transmission events between strains from the United States and European countries. Wild boar and pig isolates were genetically linked suggesting cross-border transmission and transmission due to a wildlife reservoir. The phylogenetic tree shows that multiple introductions were responsible for the outbreak of 2012–2013 in Denmark, and suggests that poorly disinfected vehicles crossing the border into Denmark were potentially the source of the outbreak. Low levels of single nucleotide polymorphisms (SNPs) differences (0–4 SNPs) can be observed between clonal strains isolated from different organs of the same animal. Proper disinfection of livestock vehicles and improved quality control of livestock feed could help to prevent future spread of S. Choleraesuis or other more serious infectious diseases such as African swine fever (ASF) in the European pig production system.

Ω-Shaped Fiber-Optic Probe-Based Localized Surface Plasmon Resonance Biosensor for Real-Time Detection of Salmonella Typhimurium

A novel, Ω-shaped fiber-optic localized surface plasmon resonance (FOLSPR) biosensor was designed for sensitive real-time and label-free bacterial detection. The designed Ω-shaped fiber-optic probe exhibits an outstanding sensitivity, due to the effect of unique geometry on performance. The results show that refractive index (RI) sensitivity of the Ω-shaped fiber-optic probe is 14 times and 2.5 times higher than those of the straight-shaped and the U-shaped FOLSPR, respectively. In addition, the reason for the geometry and the bending radius effects on RI sensitivity was discussed by investigating the relationship between RI sensitivity and the bending area. The results show that RI sensitivity was enhanced with the increase of bending area, and the best RI sensitivity obtained by Ω-shaped FOLSPR was 64.582 (a.u.)/RIU. Combined with this newly designed Ω-shaped FOLSPR biosensor, a real-time, label-free, sensitive, and highly selective bacterial detection method was established. In this work, the aptamers immobilized on the surface of FOLSPR could specifically capture Salmonella Typhimurium, resulting in an intense change of the absorption peak. In line with this principle, the FOLSPR biosensor achieved high detection sensitivity for Salmonella Typhimurium down to 128 CFU/mL within a linear range from 5 × 10² to 1 × 10⁸ CFU/mL and showed good selectivity for Salmonella Typhimurium detection compared to other bacteria. Furthermore, the FOLSPR biosensor was successfully applied to the detection of Salmonella Typhimurium in a chicken sample with the recoveries of 85-123%. With these characteristics, the novel biosensor is a potential alternative tool in food analysis and environmental monitoring.

Genome Variation and Molecular Epidemiology of Salmonella enterica Serovar Typhimurium Pathovariants

Salmonella enterica serovar Typhimurium is one of approximately 2,500 distinct serovars of the genus Salmonella but is exceptional in its wide distribution in the environment, livestock, and wild animals. S Typhimurium causes a large proportion of nontyphoidal Salmonella (NTS) infections, accounting for a quarter of infections, second only to S. enterica serovar Enteritidis in incidence. S Typhimurium was once considered the archetypal broad-host-range Salmonella serovar due to its wide distribution in livestock and wild animals, and much of what we know of the interaction of Salmonella with the host comes from research using a small number of laboratory strains of the serovar (LT2, SL1344, and ATCC 14028). But it has become clear that these strains do not reflect the genotypic or phenotypic diversity of S Typhimurium. Here, we review the epidemiological record of S Typhimurium and studies of the host-pathogen interactions of diverse strains of S Typhimurium. We present the concept of distinct pathovariants of S Typhimurium that exhibit diversity of host range, distribution in the environment, pathogenicity, and risk to food safety. We review recent evidence from whole-genome sequencing that has revealed the extent of genomic diversity of S Typhimurium pathovariants, the genomic basis of differences in the level of risk to human and animal health, and the molecular epidemiology of prominent strains. An improved understanding of the impact of genome variation of bacterial pathogens on pathogen-host and pathogen-environment interactions has the potential to improve quantitative risk assessment and reveal how new pathogens evolve.

Evolution of Salmonella within Hosts

Within-host evolution has resulted in thousands of variants of Salmonella that exhibit remarkable diversity in host range and disease outcome, from broad host range to exquisite host restriction, causing gastroenteritis to disseminated disease such as typhoid fever. Within-host evolution is a continuing process driven by genomic variation that occurs during each infection, potentiating adaptation to a new niche resulting from changes in animal husbandry, the use of antimicrobials, and emergence of immune compromised populations. We discuss key advances in our understanding of the evolution of Salmonella within the host, inferred from (i) the process of host adaptation of Salmonella pathovars in the past, and (ii) direct observation of the generation of variation and selection of beneficial traits during single infections.

Salmonella is a bacterial pathogen with remarkable diversity in its host range and pathogenicity due to past within-host evolution in vertebrate species that modified ancestral mechanisms of pathogenesis.

Variation arising during infection includes point mutations, new genes acquired through horizontal gene transfer (HGT), deletions, and genomic rearrangements.

Beneficial mutations increase in frequency within the host and, if they retain the ability to be transmitted to subsequent hosts, may become fixed in the population.

Whole-genome sequencing of sequential isolates from clinical infections reveals within-host HGT and point mutations that impact therapy and clinical management.

HGT is the primary mechanism for evolution in prokaryotes and is synergised by complex networks of transfer involving the microbiome.

Within-host evolution of Salmonella, resulting in new pathovars, can proceed in the absence of HGT.

The distribution of bacterial doubling times in the wild

Generation time varies widely across organisms and is an important factor in the life cycle, life history and evolution of organisms. Although the doubling time (DT) has been estimated for many bacteria in the laboratory, it is nearly impossible to directly measure it in the natural environment. However, an estimate can be obtained by measuring the rate at which bacteria accumulate mutations per year in the wild and the rate at which they mutate per generation in the laboratory. If we assume the mutation rate per generation is the same in the wild and in the laboratory, and that all mutations in the wild are neutral, an assumption that we show is not very important, then an estimate of the DT can be obtained by dividing the latter by the former. We estimate the DT for five species of bacteria for which we have both an accumulation and a mutation rate estimate. We also infer the distribution of DTs across all bacteria from the distribution of the accumulation and mutation rates. Both analyses suggest that DTs for bacteria in the wild are substantially greater than those in the laboratory, that they vary by orders of magnitude between different species of bacteria and that a substantial fraction of bacteria double very slowly in the wild.

Retrospective genome-wide comparisons of Salmonella enterica serovar Enteritidis from suspected outbreaks in Singapore

The number of salmonellosis cases in Singapore has increased over the years. Salmonella enterica serovar Enteritidis has always been the most predominant serovar in the last five years. The National Public Health Laboratory assisted outbreak investigations by performing multilocus variable number tandem repeat analysis (MLVA) on isolates that were collected at the time of the investigations. Isolates were defined as belonging to a particular cluster if they had identical MLVA patterns. Whilst MLVA has been instrumental in outbreak investigations, it may not be useful when outbreaks are caused by an endemic MLVA type. In this study, we analysed 67 isolates from 12 suspected outbreaks with known epidemiological links to explore the use of next-generation sequencing (NGS) for defining outbreaks. We found that NGS can confidently group isolates into their respective outbreaks. The isolates from each suspected outbreak were closely related and differed by a maximum of 3 single nucleotide polymorphisms (SNPs). They were also clearly separated from isolates that belonged to different suspected outbreaks. This study provides an important insight and further evidence on the value of NGS for routine surveillance and outbreak detection of S. Enteritidis.

Salmonella enterica subsp. salamae serovar Sofia, a prevalent serovar in Australian broiler chickens, is also capable of transient colonisation in layers

1. Salmonella enterica subsp. salamae serovar sofia (S. sofia) is a prevalent strain of Salmonella in Australian broilers and has been isolated from broiler chickens, litter, dust, as well as pre- and post-processing carcasses, and retail chicken portions but has never been reported in commercial Australian layers or eggs. 2. To investigate whether a S. sofia isolate from a broiler could colonise layers, one-month-old Hyline brown layers were orally inoculated with S. sofia and colonisation was monitored for 2-4 weeks. 3. Overall, 30-40% of the chickens shed S. sofia from the cloaca between 6 and 14 d post-inoculation which then declined to 10% by d 21. Necropsy at 2 weeks post-inoculation revealed 80% of birds harboured S. sofia in the caecum, whilst, by 4 weeks post-infection, no chickens were colonised with S. sofia in the gastrointestinal tract, liver or spleen. Additionally, no aerosol 'bird to bird' transfer was evident. 4. This study demonstrated that laying hens can be colonised by broiler-derived S. sofia; however, this colonisation was transient, reaching a peak at 14 d post-inoculation, and was completely cleared by 28 d post-inoculation. The transience of colonisation of S. sofia in layers could be a factor explaining why S. sofia has never been detected when screening for Salmonella serotypes found in Australian laying hens or eggs.

Salmonella Enteritidis ST183: emerging and endemic biotypes affecting western European hedgehogs (Erinaceus europaeus) and people in Great Britain

The impacts of hedgehog (Erinaceus europaeus) Salmonella infection on public health and on animal welfare and conservation are unknown. We isolated Salmonella Enteritidis multi-locus sequence-type (ST)183 from 46/170 (27%) hedgehog carcasses (27 S. Enteritidis phage type (PT)11, 18 of a novel PT66 biotype and one with co-infection of these PTs) and from 6/208 (3%) hedgehog faecal samples (4 PT11, 2 PT66) from across Great Britain, 2012–2015. Whole genome phylogenetic analysis of the hedgehog isolates and ST183 from people in England and Wales found that PT11 and PT66 form two divergent clades. Hedgehog and human isolates were interspersed throughout the phylogeny indicating that infections in both species originate from a common population. PT11 was recovered from hedgehogs across England and Scotland, consistent with endemic infection. PT66 was isolated from Scotland only, possibly indicating a recent emergence event. People infected with ST183 were four times more likely to be aged 0–4 years than people infected by the more common ST11 S. Enteritidis. Evidence for human ST183 infection being non-foodborne included stronger correlation between geographic and genetic distance, and significantly increased likelihood of infection in rural areas, than for ST11. These results are consistent with hedgehogs acting as a source of zoonotic infection.

Comparison of advanced whole genome sequence-based methods to distinguish strains of Salmonella enterica serovar Heidelberg involved in foodborne outbreaks in Québec

Salmonella enterica serovar Heidelberg (S. Heidelberg) is one of the top serovars causing human salmonellosis. This serovar ranks second and third among serovars that cause human infections in Québec and Canada, respectively, and has been associated with severe infections. Traditional typing methods such as PFGE do not display adequate discrimination required to resolve outbreak investigations due to the low level of genetic diversity of isolates belonging to this serovar. This study evaluates the ability of four whole genome sequence (WGS)-based typing methods to differentiate among 145 S. Heidelberg strains involved in four distinct outbreak events and sporadic cases of salmonellosis that occurred in Québec between 2007 and 2016. Isolates from all outbreaks were indistinguishable by PFGE. The core genome single nucleotide variant (SNV), core genome multilocus sequence typing (MLST) and whole genome MLST approaches were highly discriminatory and separated outbreak strains into four distinct phylogenetic clusters that were concordant with the epidemiological data. The clustered regularly interspaced short palindromic repeats (CRISPR) typing method was less discriminatory. However, CRISPR typing may be used as a secondary method to differentiate isolates of S. Heidelberg that are genetically similar but epidemiologically unrelated to outbreak events. WGS-based typing methods provide a highly discriminatory alternative to PFGE for the laboratory investigation of foodborne outbreaks.

Comparative whole genome analysis of three consecutive Salmonella diarizonae isolates

Infections of very young children or immunocompromised people with Salmonella of higher subspecies are a well-known phenomenon often associated with contact to cold-blooded animals. We describe the molecular characterization of three S. enterica subsp. diarizonae strains, isolated consecutively over a period of several months from a hospital patient suffering from diarrhea and sepsis with fatal outcome. With the initial isolate the first complete genome sequence of a member of subsp. diarizonae is provided and based on this reference we revealed the genomic differences between the three isolates by use of next-generation sequencing and confirmed by phenotypical tests. Genome comparisons revealed mutations within gpt, hfq and purK in the first isolate as a sign of clonal variation rather than host-directed evolution. Furthermore, our work demonstrates that S. enterica subsp. diarizonae possess, besides a conserved set of known Salmonella Pathogenicity Islands, a variable portfolio of additional genomic islands of unknown function.

Evaluation of Whole Genome Mapping as a Fast and Automated Molecular Epidemiological Tool for the Study of Cronobacter spp. in Powdered Infant Formula Processing Facilities

Cronobacter has been identified as the causative agent of outbreaks or sporadic cases of meningitis, necrotizing enterocolitis, and septicemia associated with powdered infant formula. Food processing environments may provide a possible contamination route. The purpose of this study was to evaluate whole genome mapping (WGM) as a fast and automated molecular epidemiological method for characterizing Cronobacter spp. in the processing environment. This is the first study indicating the applicability of WGM to Cronobacter. WGM was compared with ribotyping, which is often used as an automated typing tool, and with pulsed-field gel electrophoresis, which is a well-known and highly discriminating tool that is also based on restriction site analysis. The comparison of the three tools was carried out on a subset of Cronobacter isolates collected from 2011 to 2014 through a monitoring program. The performance characteristics of WGM have not yet been described; therefore, in the current study its performance was evaluated based on five criteria: typeability, reproducibility, stability, epidemiological concordance, and the discrimination power. WGM was shown to produce typeable, reproducible, and stable results. With a similar cut-off of 98%, WGM was shown to have a discriminatory power equivalent to pulsed-field gel electrophoresis and higher than ribotyping. Future studies are needed to confirm the indicated cut-off level of 98%.

Comparative genomics of Australian and international isolates of Salmonella Typhimurium: correlation of core genome evolution with CRISPR and prophage profiles

Salmonella enterica subsp enterica serovar Typhimurium (S. Typhimurium) is a serovar with broad host range. To determine the genomic diversity of S. Typhimurium, we sequenced 39 isolates (37 Australian and 2 UK isolates) representing 14 Repeats Groups (RGs) determined primarily by clustered regularly interspaced short palindromic repeats (CRISPR). Analysis of single nucleotide polymorphisms (SNPs) among the 39 isolates yielded an average of 1,232 SNPs per isolate, ranging from 128 SNPs to 11,339 SNPs relative to the reference strain LT2. Phylogenetic analysis of the 39 isolates together with 66 publicly available genomes divided the 105 isolates into five clades and 19 lineages, with the majority of the isolates belonging to clades I and II. The composition of CRISPR profiles correlated well with the lineages, showing progressive deletion and occasional duplication of spacers. Prophage genes contributed nearly a quarter of the S. Typhimurium accessory genome. Prophage profiles were found to be correlated with lineages and CRISPR profiles. Three new variants of HP2-like P2 prophage, several new variants of P22 prophage and a plasmid-like genomic island StmGI_0323 were found. This study presents evidence of horizontal transfer from other serovars or species and provides a broader understanding of the global genomic diversity of S. Typhimurium.

Salmonella Typhimurium gastroenteritis leading to chronic prosthetic vascular graft infection

Introduction. It is estimated up to 6 % of prosthetic vascular grafts become infected. Staphylococcus aureus is predominant in early infection and coagulase-negative staphylococci are predominant in late infections. Enterobacteriaceae cause 14–40 % of prosthetic vascular graft infections. This is, to our knowledge the first reported case of Salmonella gastroenteritis causing chronic prosthetic vascular graft infection (PVGI).

Case presentation. A 57 years old lady presented with signs and symptoms of prosthetic vascular graft infection. Three years earlier, she had undergone a prosthetic axillo-femoral bypass graft for critical limb ischaemia. The infected prosthetic vascular graft was removed and Salmonella Typhimurium was isolated on culture. In the intervening period, Salmonella Typhimurium was isolated from a faecal specimen, collected during an episode of acute gastroenteritis. Whole-genome sequencing (WGS) showed that the respective Salmonella Typhimurium isolates differed by only a single nucleotide polymorphism (SNP). Salmonella Typhimurium was not isolated on culture of a faecal specimen collected five days following cessation of antimicrobial therapy. Six months after removal of the prosthetic graft, the patient remains under follow-up for her peripheral vascular disease, which currently requires no further surgical intervention.

Conclusion. This case has clear implications for the management of chronic PVGI. It is vital to collect high-quality surgical specimens for microbiological analysis and empirical choices of antibiotics are unlikely to cover all potential pathogens. It may also be prudent to enquire about a history of acute gastroenteritis when assessing patients presenting with chronic PVGI.

Evolution of Variable Number Tandem Repeats and Its Relationship with Genomic Diversity in Salmonella Typhimurium

Salmonella enterica serovar Typhimurium is the most common Salmonella serovar causing human infections in Australia and many other countries. A total of 12,112 S. Typhimurium isolates from New South Wales were analyzed by multi-locus variable number of tandem repeat (VNTR) analysis (MLVA) using five VNTRs from 2007 to 2014. We found that mid ranges of repeat units of 8–14 in VNTR locus STTR5, 6–13 in STTR6, and 9–12 in STTR10 were always predominant in the population (>50%). In vitro passaging experiments using MLVA type carrying extreme length alleles found that the majority of long length alleles mutated to short ones and short length alleles mutated to longer ones. Both data suggest directional mutability of VNTRs toward mid-range repeats. Sequencing of 28 isolates from a newly emerged MLVA type and its five single locus variants revealed that single nucleotide variation between isolates with up to two MLVA differences ranged from 0 to 12 single nucleotide polymorphisms (SNPs). However, there was no relationship between SNP and VNTR differences. A population genetic model of the joint distribution of VNTRs and SNPs variations was used to estimate the mutation rates of the two markers, yielding a ratio of 1 VNTR change to 6.9 SNP changes. When only one VNTR repeat difference was considered, the majority of pairwise SNP difference between isolates were 4 SNPs or fewer. Based on this observation and our previous findings of SNP differences of outbreak isolates, we suggest that investigation of S. Typhimurium community outbreaks should include cases of 1 repeat difference to increase sensitivity. This study offers new insights into the short-term VNTR evolution of S. Typhimurium and its application for epidemiological typing.

Whole genome sequencing of Salmonella Typhimurium illuminates distinct outbreaks caused by an endemic multi-locus variable number tandem repeat analysis type in Australia, 2014

Background

Salmonella Typhimurium (STM) is an important cause of foodborne outbreaks worldwide. Subtyping of STM remains critical to outbreak investigation, yet current techniques (e.g. multilocus variable number tandem repeat analysis, MLVA) may provide insufficient discrimination. Whole genome sequencing (WGS) offers potentially greater discriminatory power to support infectious disease surveillance.

Methods

We performed WGS on 62 STM isolates of a single, endemic MLVA type associated with two epidemiologically independent, food-borne outbreaks along with sporadic cases in New South Wales, Australia, during 2014. Genomes of case and environmental isolates were sequenced using HiSeq (Illumina) and the genetic distance between them was assessed by single nucleotide polymorphism (SNP) analysis. SNP analysis was compared to the epidemiological context.

Results

The WGS analysis supported epidemiological evidence and genomes of within-outbreak isolates were nearly identical. Sporadic cases differed from outbreak cases by a small number of SNPs, although their close relationship to outbreak cases may represent an unidentified common food source that may warrant further public health follow up. Previously unrecognised mini-clusters were detected.

Conclusions

WGS of STM can discriminate foodborne community outbreaks within a single endemic MLVA clone. Our findings support the translation of WGS into public health laboratory surveillance of salmonellosis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-016-0831-3) contains supplementary material, which is available to authorized users.

Phylogenetic structure of European Salmonella Enteritidis outbreak correlates with national and international egg distribution network

Outbreaks of Salmonella Enteritidis have long been associated with contaminated poultry and eggs. In the summer of 2014 a large multi-national outbreak of Salmonella Enteritidis phage type 14b occurred with over 350 cases reported in the United Kingdom, Germany, Austria, France and Luxembourg. Egg supply network investigation and microbiological sampling identified the source to be a Bavarian egg producer. As part of the international investigation into the outbreak, over 400 isolates were sequenced including isolates from cases, implicated UK premises and eggs from the suspected source producer. We were able to show a clear statistical correlation between the topology of the UK egg distribution network and the phylogenetic network of outbreak isolates. This correlation can most plausibly be explained by different parts of the egg distribution network being supplied by eggs solely from independent premises of the Bavarian egg producer (Company X). Microbiological sampling from the source premises, traceback information and information on the interventions carried out at the egg production premises all supported this conclusion. The level of insight into the outbreak epidemiology provided by whole-genome sequencing (WGS) would not have been possible using traditional microbial typing methods.

Tissue colonization and circulating T lymphocytes in laying hens upon oral challenge with Salmonella enterica serovars

Evaluating the potential of Salmonella serovars for tissue colonization and egg contamination in laying hens is critical due to widespread consumption of poultry and egg-containing products. The 2009 FDA Egg Rule was implemented to target the eradication of Salmonella enterica Enteritidis (SE) from layers; however, other Salmonella serovars, such as Heidelberg (SH) and Typhimurium (ST), have also been associated with poultry-related outbreaks. We conducted this study to see if serovars other than SE could colonize in laying hens, cause egg contamination, and modulate circulating T-cell populations. Laying hens were orally gavaged with 10⁷ colony forming units (CFU) of SE, SH, or ST and assessed for colonization in spleen, ovaries, and oviduct 10 d postchallenge. Splenic colonization was similar for all the serovars; however, colonization of ovaries and oviducts was significantly higher with SH compared to SE and ST. Furthermore, SH challenge resulted in egg contamination, while SE and ST did not result in contaminated eggs. Phenotypic evaluation of peripheral blood lymphocytes showed significant reduction in CD4 cells in SH-challenged birds and lower CD8α and CD8β cells in SE-challenged birds compared to controls. Our data showed that non-SE serovars have equal or higher potential to colonize reproductive tissues of laying hens and may be accompanied by altered lymphocyte populations.

Heatwaves differentially affect risk of Salmonella serotypes

OBJECTIVES

Given increasing frequency of heatwaves and growing public health concerns associated with foodborne disease, we examined the relationship between heatwaves and salmonellosis in Adelaide, Australia.

METHODS

Poisson regression analysis with Generalised Estimating Equations was used to estimate the effect of heatwaves and the impact of intensity, duration and timing on salmonellosis and specific serotypes notified from 1990 to 2012. Distributed lag non-linear models were applied to assess the non-linear and delayed effects of temperature during heatwaves on Salmonella cases.

RESULTS

Salmonella typhimurium PT135 notifications were sensitive to the effects of heatwaves with a twofold (IRR 2.08, 95% CI 1.14-3.79) increase in cases relative to non-heatwave days. Heatwave intensity had a significant effect on daily counts of overall salmonellosis with a 34% increase in risk of infection (IRR 1.34, 95% CI 1.01-1.78) at >41 °C. The effects of temperature during heatwaves on Salmonella cases and serotypes were found at lags of up to 14 days.

CONCLUSION

This study confirms heatwaves have a significant effect on Salmonella cases, and for the first time, identifies its impact on specific serotypes and phage types. These findings will contribute to the understanding of the impact of heatwaves on salmonellosis and provide insights that could mitigate their impact.

Pangenome Evolution in the Marine Bacterium Alteromonas

We have examined a collection of the free-living marine bacterium Alteromonas genomes with cores diverging in average nucleotide identities ranging from 99.98% to 73.35%, i.e., from microbes that can be considered members of a natural clone (like in a clinical epidemiological outbreak) to borderline genus level. The genomes were largely syntenic allowing a precise delimitation of the core and flexible regions in each. The core was 1.4 Mb (ca. 30% of the typical strain genome size). Recombination rates along the core were high among strains belonging to the same species (37.7-83.7% of all nucleotide polymorphisms) but they decreased sharply between species (18.9-5.1%). Regarding the flexible genome, its main expansion occurred within the boundaries of the species, i.e., strains of the same species already have a large and diverse flexible genome. Flexible regions occupy mostly fixed genomic locations. Four large genomic islands are involved in the synthesis of strain-specific glycosydic receptors that we have called glycotypes. These genomic regions are exchanged by homologous recombination within and between species and there is evidence for their import from distant taxonomic units (other genera within the family). In addition, several hotspots for integration of gene cassettes by illegitimate recombination are distributed throughout the genome. They code for features that give each clone specific properties to interact with their ecological niche and must flow fast throughout the whole genus as they are found, with nearly identical sequences, in different species. Models for the generation of this genomic diversity involving phage predation are discussed.

Genome Analysis of Listeria monocytogenes Sequence Type 8 Strains Persisting in Salmon and Poultry Processing Environments and Comparison with Related Strains

Listeria monocytogenes is an important foodborne pathogen responsible for the disease listeriosis, and can be found throughout the environment, in many foods and in food processing facilities. The main cause of listeriosis is consumption of food contaminated from sources in food processing environments. Persistence in food processing facilities has previously been shown for the L. monocytogenes sequence type (ST) 8 subtype. In the current study, five ST8 strains were subjected to whole-genome sequencing and compared with five additionally available ST8 genomes, allowing comparison of strains from salmon, poultry and cheese industry, in addition to a human clinical isolate. Genome-wide analysis of single-nucleotide polymorphisms (SNPs) confirmed that almost identical strains were detected in a Danish salmon processing plant in 1996 and in a Norwegian salmon processing plant in 2001 and 2011. Furthermore, we show that L. monocytogenes ST8 was likely to have been transferred between two poultry processing plants as a result of relocation of processing equipment. The SNP data were used to infer the phylogeny of the ST8 strains, separating them into two main genetic groups. Within each group, the plasmid and prophage content was almost entirely conserved, but between groups, these sequences showed strong divergence. The accessory genome of the ST8 strains harbored genetic elements which could be involved in rendering the ST8 strains resilient to incoming mobile genetic elements. These included two restriction-modification loci, one of which was predicted to show phase variable recognition sequence specificity through site-specific domain shuffling. Analysis indicated that the ST8 strains harbor all important known L. monocytogenes virulence factors, and ST8 strains are commonly identified as the causative agents of invasive listeriosis. Therefore, the persistence of this L. monocytogenes subtype in food processing facilities poses a significant concern for food safety.

Salmonella Typhimurium and Salmonella Sofia: Growth in and Persistence on Eggs under Production and Retail Conditions

Salmonellosis in Australia has been linked to eggs and egg products with specific serotypes associated with outbreaks. We compared attachment to and survival on egg shells and growth in eggs of two Salmonella serotypes, an egg outbreak associated Salmonella Typhimurium and a non-egg-associated Salmonella enterica ssp. II 1,4,12,27:b:[e,n,x] (S. Sofia). Experiments were conducted at combinations of 4, 15, 22, 37 and 42°C. No significant differences occurred between the serotypes in maximum growth rates, which were significantly greater (P < 0.001) in egg yolk (0.427 log₁₀ CFU/mL/h) compared to whole egg (0.312 log₁₀ CFU/mL/h) and egg white (0.029 log₁₀ CFU/mL/h). Attachment to egg shells varied by time (1 or 20 min) and temperature (4, 22 and 42°C), with S. Typhimurium isolates attaching at higher levels (P < 0.05) than S. Sofia after 1 min at 4°C and S. Typhimurium ATCC 14028 attaching at higher (P < 0.05) levels at 22°C. Survival on egg shells was not significantly different across isolates. Salmonella serotypes behaved similarly regarding growth in egg contents, attachment to egg shells and survival on eggs, indicating that other factors more likely contributed to reasons for S. Typhimurium being implicated in multiple egg-associated outbreaks.

It Is Not All about Single Nucleotide Polymorphisms: Comparison of Mobile Genetic Elements and Deletions in Listeria monocytogenes Genomes Links Cases of Hospital-Acquired Listeriosis to the Environmental Source

The control of food-borne outbreaks caused by Listeria monocytogenes in humans relies on the timely identification of food or environmental sources and the differentiation of outbreak-related isolates from unrelated ones. This study illustrates the utility of whole-genome sequencing for examining the link between clinical and environmental isolates of L. monocytogenes associated with an outbreak of hospital-acquired listeriosis in Sydney, Australia. Comparative genomic analysis confirmed an epidemiological link between the three clinical and two environmental isolates. Single nucleotide polymorphism (SNP) analysis showed that only two SNPs separated the three human outbreak isolates, which differed by 19 to 20 SNPs from the environmental isolates and 71 to >10,000 SNPs from sporadic L. monocytogenes isolates. The chromosomes of all human outbreak isolates and the two suspected environmental isolates were syntenic. In contrast to the genomes of background sporadic isolates, all epidemiologically linked isolates contained two novel prophages and a previously unreported clustered regularly interspaced short palindromic repeat (CRISPR)-associated (Cas) locus subtype sequence. The mobile genetic element (MGE) profile of these isolates was distinct from that of the other serotype 1/2b reference strains and sporadic isolates. The identification of SNPs and clonally distinctive MGEs strengthened evidence to distinguish outbreak-related isolates of L. monocytogenes from cocirculating endemic strains.

Genomic Variability of Serial Human Isolates of Salmonella enterica Serovar Typhimurium Associated with Prolonged Carriage

Salmonella enterica serovar Typhimurium is an important foodborne human pathogen that often causes self-limiting but severe gastroenteritis. Prolonged excretion of S. Typhimurium after the infection can lead to secondary transmissions. However, little is known about within-host genomic variation in bacteria associated with asymptomatic shedding. Genomes of 35 longitudinal isolates of S. Typhimurium recovered from 11 patients (children and adults) with culture-confirmed gastroenteritis were sequenced. There were three or four isolates obtained from each patient. Single nucleotide polymorphisms (SNPs) were analyzed in these isolates, which were recovered between 1 and 279 days after the initial diagnosis. Limited genomic variation (5 SNPs or fewer) was associated with short- and long-term carriage of S. Typhimurium. None of the isolates was shown to be due to reinfection. SNPs occurred randomly, and the majority of the SNPs were nonsynonymous. Two nonsense mutations were observed. A nonsense mutation in flhC rendered the isolate nonmotile, whereas the significance of a nonsense mutation in yihV is unknown. The estimated mutation rate is 1.49 × 10(-6) substitution per site per year. S. Typhimurium isolates excreted in stools following acute gastroenteritis in children and adults demonstrated limited genomic variability over time, regardless of the duration of carriage. These findings have important implications for the detection of possible transmission events suspected by public health genomic surveillance of S. Typhimurium infections.

Defining the Core Genome of Salmonella enterica Serovar Typhimurium for Genomic Surveillance and Epidemiological Typing

Salmonella enterica serovar Typhimurium is the most common Salmonella serovar causing foodborne infections in Australia and many other countries. Twenty-one S. Typhimurium strains from Salmonella reference collection A (SARA) were analyzed using Illumina high-throughput genome sequencing. Single nucleotide polymorphisms (SNPs) in 21 SARA strains ranged from 46 to 11,916 SNPs, with an average of 1,577 SNPs per strain. Together with 47 strains selected from publicly available S. Typhimurium genomes, the S. Typhimurium core genes (STCG) were determined. The STCG consist of 3,846 genes, a set that is much larger than that of the 2,882 Salmonella core genes (SCG) found previously. The STCG together with 1,576 core intergenic regions (IGRs) were defined as the S. Typhimurium core genome. Using 93 S. Typhimurium genomes from 13 epidemiologically confirmed community outbreaks, we demonstrated that typing based on the S. Typhimurium core genome (STCG plus core IGRs) provides superior resolution and higher discriminatory power than that based on SCG for outbreak investigation and molecular epidemiology of S. Typhimurium. STCG and STCG plus core IGR typing achieved 100% separation of all outbreaks compared to that of SCG typing, which failed to separate isolates from two outbreaks from background isolates. Defining the S. Typhimurium core genome allows standardization of genes/regions to be used for high-resolution epidemiological typing and genomic surveillance of S. Typhimurium.

Whole Genome Sequencing for the Retrospective Investigation of an Outbreak of Salmonella Typhimurium DT 8

Background: Salmonella enterica serovar Typhimurium DT8 is uncommon within the European Union. An increase in this phage type was reported in the summer of 2013 in the States of Jersey. Methods: A total of 21 human cases with this phage type were microbiologically confirmed. Salmonella isolates from mayonnaise made using raw eggs were also confirmed as being Salmonella Typhimurium DT8. The epidemiological investigations strongly supported a link between mayonnaise consumption and illness. Whole genome sequencing (WGS) was used to retrospectively investigate this outbreak with a view to assess the similarity between the suspect food and the human isolates and to characterise a known point source outbreak to assist in development of algorithms for outbreak detection. Results: Sequence data showed that the outbreak associated isolates, including the food isolates, formed a tightly clustered monophyletic group, with a maximum pairwise distance of 3 single nucleotide polymorphisms. Conclusions: WGS data is useful in confirming the causative agent of outbreaks where food and clinical isolates are available. This dataset, comprising a known outbreak, will be useful in the development of automatic algorithms for outbreak detection.

Delineating community outbreaks of Salmonella enterica serovar Typhimurium by use of whole-genome sequencing: insights into genomic variability within an outbreak

Whole-genome next-generation sequencing (NGS) was used to retrospectively examine 57 isolates from five epidemiologically confirmed community outbreaks (numbered 1 to 5) caused by Salmonella enterica serovar Typhimurium phage type DT170. Most of the human and environmental isolates confirmed epidemiologically to be involved in the outbreaks were either genomically identical or differed by one or two single nucleotide polymorphisms (SNPs), with the exception of those in outbreak 1. The isolates from outbreak 1 differed by up to 12 SNPs, which suggests that the food source of the outbreak was contaminated with more than one strain while each of the other four outbreaks was caused by a single strain. In addition, NGS analysis ruled in isolates that were initially not considered to be linked with the outbreak, which increased the total outbreak size by 107%. The mutation process was modeled by using known mutation rates to derive a cutoff value for the number of SNP difference to determine whether or not a case was part of an outbreak. For an outbreak with less than 1 month of ex vivo/in vivo evolution time, the maximum number of SNP differences between isolates is two or four using the lowest or highest mutation rate, respectively. NGS of S. Typhimurium significantly increases the resolution of investigations of community outbreaks. It can also inform a more targeted public health response by providing important supplementary evidence that cases of disease are or are not associated with food-borne outbreaks of S. Typhimurium.

Genomic signature of multidrug-resistant Salmonella enterica serovar typhi isolates related to a massive outbreak in Zambia between 2010 and 2012

Retrospectively, we investigated the epidemiology of a massive Salmonella enterica serovar Typhi outbreak in Zambia during 2010 to 2012. Ninety-four isolates were susceptibility tested by MIC determinations. Whole-genome sequence typing (WGST) of 33 isolates and bioinformatic analysis identified the multilocus sequence type (MLST), haplotype, plasmid replicon, antimicrobial resistance genes, and genetic relatedness by single nucleotide polymorphism (SNP) analysis and genomic deletions. The outbreak affected 2,040 patients, with a fatality rate of 0.5%. Most (83.0%) isolates were multidrug resistant (MDR). The isolates belonged to MLST ST1 and a new variant of the haplotype, H58B. Most isolates contained a chromosomally translocated region containing seven antimicrobial resistance genes, catA1, blaTEM-1, dfrA7, sul1, sul2, strA, and strB, and fragments of the incompatibility group Q1 (IncQ1) plasmid replicon, the class 1 integron, and the mer operon. The genomic analysis revealed 415 SNP differences overall and 35 deletions among 33 of the isolates subjected to whole-genome sequencing. In comparison with other genomes of H58, the Zambian isolates separated from genomes from Central Africa and India by 34 and 52 SNPs, respectively. The phylogenetic analysis indicates that 32 of the 33 isolates sequenced belonged to a tight clonal group distinct from other H58 genomes included in the study. The small numbers of SNPs identified within this group are consistent with the short-term transmission that can be expected over a period of 2 years. The phylogenetic analysis and deletions suggest that a single MDR clone was responsible for the outbreak, during which occasional other S. Typhi lineages, including sensitive ones, continued to cocirculate. The common view is that the emerging global S. Typhi haplotype, H58B, containing the MDR IncHI1 plasmid is responsible for the majority of typhoid infections in Asia and sub-Saharan Africa; we found that a new variant of the haplotype harboring a chromosomally translocated region containing the MDR islands of IncHI1 plasmid has emerged in Zambia. This could change the perception of the term "classical MDR typhoid" currently being solely associated with the IncHI1 plasmid. It might be more common than presently thought that S. Typhi haplotype H58B harbors the IncHI1 plasmid or a chromosomally translocated MDR region or both.

Development of a new molecular subtyping tool for Salmonella enterica serovar Enteritidis based on single nucleotide polymorphism genotyping using PCR

The lack of a sufficiently discriminatory molecular subtyping tool for Salmonella enterica serovar Enteritidis has hindered source attribution efforts and impeded regulatory actions required to disrupt its food-borne transmission. The underlying biological reason for the ineffectiveness of current molecular subtyping tools such as pulsed-field gel electrophoresis (PFGE) and phage typing appears to be related to the high degree of clonality of S. Enteritidis. By interrogating the organism's genome, we previously identified single nucleotide polymorphisms (SNP) distributed throughout the chromosome and have designed a highly discriminatory PCR-based SNP typing test based on 60 polymorphic loci. The application of the SNP-PCR method to DNA samples from S. Enteritidis strains (n = 55) obtained from a variety of sources has led to the differentiation and clustering of the S. Enteritidis isolates into 12 clades made up of 2 to 9 isolates per clade. Significantly, the SNP-PCR assay was able to further differentiate predominant PFGE types (e.g., XAI.0003) and phage types (e.g., phage type 8) into smaller subsets. The SNP-PCR subtyping test proved to be an accurate, precise, and quantitative tool for evaluating the relationships among the S. Enteritidis isolates tested in this study and should prove useful for clustering related S. Enteritidis isolates involved in outbreaks.

First Complete Genome Sequence of Salmonella enterica subsp. enterica Serovar Typhimurium Strain ATCC 13311 (NCTC 74), a Reference Strain of Multidrug Resistance, as Achieved by Use of PacBio Single-Molecule Real-Time Technology

We report the first complete genomic sequence of Salmonella enterica subsp. enterica serovar Typhimurium strain ATCC 13311, the leading food-borne pathogen and a reference strain used in drug resistance studies. De novo assembly with PacBio sequencing completed its chromosome and one plasmid. They will accelerate the investigation into multidrug resistance in Salmonella Typhimurium.

Analysis of Salmonella enterica serovar Typhimurium variable-number tandem-repeat data for public health investigation based on measured mutation rates and whole-genome sequence comparisons

Variable-number tandem repeats (VNTRs) mutate rapidly and can be useful markers for genotyping. While multilocus VNTR analysis (MLVA) is increasingly used in the detection and investigation of food-borne outbreaks caused by Salmonella enterica serovar Typhimurium (S. Typhimurium) and other bacterial pathogens, MLVA data analysis usually relies on simple clustering approaches that may lead to incorrect interpretations. Here, we estimated the rates of copy number change at each of the five loci commonly used for S. Typhimurium MLVA, during in vitro and in vivo passage. We found that loci STTR5, STTR6, and STTR10 changed during passage but STTR3 and STTR9 did not. Relative rates of change were consistent across in vitro and in vivo growth and could be accurately estimated from diversity measures of natural variation observed during large outbreaks. Using a set of 203 isolates from a series of linked outbreaks and whole-genome sequencing of 12 representative isolates, we assessed the accuracy and utility of several alternative methods for analyzing and interpreting S. Typhimurium MLVA data. We show that eBURST analysis was accurate and informative. For construction of MLVA-based trees, a novel distance metric, based on the geometric model of VNTR evolution coupled with locus-specific weights, performed better than the commonly used simple or categorical distance metrics. The data suggest that, for the purpose of identifying potential transmission clusters for further investigation, isolates whose profiles differ at one of the rapidly changing STTR5, STTR6, and STTR10 loci should be collapsed into the same cluster.