Genomic diversity of Salmonella enterica -The UoWUCC 10K genomes project

A new Salmonella enterica serovar that was isolated from a wild sparrow presents a distinct genetic, metabolic and virulence profile

Salmonella enterica is a ubiquitous and clinically-important bacterial pathogen, able to infect and cause different diseases in a wide range of hosts. Here, we report the isolation and characterization of a new S. enterica serovar (13,23:i:-; S. Tirat-Zvi), belonging to the Havana supper-lineage that was isolated from a wild house sparrow (Passer domesticus) in Israel. Whole genome sequencing and complete assembly of its genome indicated a plasmid-free, 4.7 Mb genome that carries the Salmonella pathogenicity islands 1-6, 9, 19 and an integrative and conjugative element (ICE), encoding arsenic resistance genes. Phenotypically, S. Tirat-Zvi isolate TZ282 was motile, readily formed biofilm, more versatile in carbon source utilization than S. Typhimurium and highly tolerant to arsenic, but impaired in host cell invasion. In-vivo infection studies indicated that while S. Tirat-Zvi was able to infect and cause an acute inflammatory enterocolitis in young chicks, it was compromised in mice colonization and did not cause an inflammatory colitis in mice compared to S. Typhimurium. We suggest that these phenotypes reflect the distinctive ecological niche of this new serovar and its evolutionary adaptation to passerine birds, as a permissive host. Moreover, these results further illuminate the genetic, phenotypic and ecological diversity of S. enterica pathovars.

Genomic Features and Phylogenetic Analysis of Antimicrobial-Resistant Salmonella Mbandaka ST413 Strains

In recent years, Salmonella enterica subsp. enterica serovar Mbandaka (S. Mbandaka) has been increasingly isolated from laying hens and shell eggs around the world. Moreover, this serovar has been identified as the causative agent of several salmonellosis outbreaks in humans. Surprisingly, little is known about the characteristics of this emerging serovar, and therefore, we investigated antimicrobial resistance, virulence, and prophage genes of six selected Brazilian strains of Salmonella Mbandaka using Whole Genome Sequencing (WGS). Multi-locus sequence typing revealed that the tested strains belong to Sequence Type 413 (ST413), which has been linked to recent multi-country salmonellosis outbreaks in Europe. A total of nine resistance genes were detected, and the most frequent ones were aac(6')-Iaa, sul1, qacE, blaOXA-129, tet(B), and aadA1. A point mutation in ParC at the 57th position (threonine → serine) associated with quinolone resistance was present in all investigated genomes. A 112,960 bp IncHI2A plasmid was mapped in 4/6 strains. This plasmid harboured tetracycline (tetACDR) and mercury (mer) resistance genes, genes contributing to conjugative transfer, and genes involved in plasmid maintenance. Most strains (four/six) carried Salmonella genomic island 1 (SGI1). All S. Mbandaka genomes carried seven pathogenicity islands (SPIs) involved in intracellular survival and virulence: SPIs 1-5, 9, and C63PI. The virulence genes csgC, fimY, tcfA, sscA, (two/six), and ssaS (one/six) were absent in some of the genomes; conversely, fimA, prgH, and mgtC were present in all of them. Five Salmonella bacteriophage sequences (with homology to Escherichia phage phiV10, Enterobacteria phage Fels-2, Enterobacteria phage HK542, Enterobacteria phage ST64T, Salmonella phage SW9) were identified, with protein counts between 31 and 54, genome lengths of 24.7 bp and 47.7 bp, and average GC content of 51.25%. In the phylogenetic analysis, the genomes of strains isolated from poultry in Brazil clustered into well-supported clades with a heterogeneous distribution, primarily associated with strains isolated from humans and food. The phylogenetic relationship of Brazilian S. Mbandaka suggests the presence of strains with high epidemiological significance and the potential to be linked to foodborne outbreaks. Overall, our results show that isolated strains of S. Mbandaka are multidrug-resistant and encode a rather conserved virulence machinery, which is an epidemiological hallmark of Salmonella strains that have successfully disseminated both regionally and globally.

Themisto: a scalable colored k-mer index for sensitive pseudoalignment against hundreds of thousands of bacterial genomes

MOTIVATION

Huge datasets containing whole-genome sequences of bacterial strains are now commonplace and represent a rich and important resource for modern genomic epidemiology and metagenomics. In order to efficiently make use of these datasets, efficient indexing data structures-that are both scalable and provide rapid query throughput-are paramount.

RESULTS

Here, we present Themisto, a scalable colored k-mer index designed for large collections of microbial reference genomes, that works for both short and long read data. Themisto indexes 179 thousand Salmonella enterica genomes in 9 h. The resulting index takes 142 gigabytes. In comparison, the best competing tools Metagraph and Bifrost were only able to index 11 000 genomes in the same time. In pseudoalignment, these other tools were either an order of magnitude slower than Themisto, or used an order of magnitude more memory. Themisto also offers superior pseudoalignment quality, achieving a higher recall than previous methods on Nanopore read sets.

AVAILABILITY AND IMPLEMENTATION

Themisto is available and documented as a C++ package at https://github.com/algbio/themisto available under the GPLv2 license.

Functionally distinct mutations within AcrB underpin antibiotic resistance in different lifestyles

Antibiotic resistance is a pressing healthcare challenge and is mediated by various mechanisms, including the active export of drugs via multidrug efflux systems, which prevent drug accumulation within the cell. Here, we studied how Salmonella evolved resistance to two key antibiotics, cefotaxime and azithromycin, when grown planktonically or as a biofilm. Resistance to both drugs emerged in both conditions and was associated with different substitutions within the efflux-associated transporter, AcrB. Azithromycin exposure selected for an R717L substitution, while cefotaxime for Q176K. Additional mutations in ramR or envZ accumulated concurrently with the R717L or Q176K substitutions respectively, resulting in clinical resistance to the selective antibiotics and cross-resistance to other drugs. Structural, genetic, and phenotypic analysis showed the two AcrB substitutions confer their benefits in profoundly different ways. R717L reduces steric barriers associated with transit through the substrate channel 2 of AcrB. Q176K increases binding energy for cefotaxime, improving recognition in the distal binding pocket, resulting in increased efflux efficiency. Finally, we show the R717 substitution is present in isolates recovered around the world.

Heightened variability observed in resistance and virulence genes across salmonella Kentucky isolates from poultry environments in British Columbia, Canada

Many niche-dependent barriers along the poultry production continuum favour the survival of certain Salmonella serovars over others. Historically, the presence of particular serovars has been determined by niche-specific genes which encode resistance to selective pressures such as host defenses and industrial antimicrobial practices. Over the past decade, Canada has witnessed unexplained shifts in the Salmonella landscape in the poultry sector. Several formerly minor Salmonella serovars, including S. Kentucky and S. Reading, have recently increased in prevalence in live chickens and turkeys, respectively, in British Columbia (BC). The purpose of this research was to investigate the genomic features of the top poultry-associated Salmonella spp. in BC, to probe for serovar-specific characteristics that could address the recently shifting balance of serovars along the poultry continuum. By examining the quantity and diversity of antimicrobial resistance (AMR) genes, virulence factors (VFs), Salmonella Pathogenicity Islands (SPIs), and plasmids across 50 poultry-associated S. enterica isolates using whole genome sequencing and antimicrobial resistance profiling, we have identified serovar-specific differences that likely influence niche survival. Specifically, isolates in our collection from predominantly human pathogenic serovars (S. I 4, [5], 12:i: , S. Typhimurium, and S. Enteritidis) were found to share the IncFIB(S) and IncFII(S) plasmids which carry important VFs known to aid in human host infection. Additionally, these strains held the lowest number of AMR genes, and the highest number of unique SPIs which also facilitate virulence. However, other serovars containing a greater diversity and abundance of resistance genes have been increasing across the poultry sector. S. Kentucky was found to carry unique AMR genes, VFs, SPIs, and plasmids that could bolster persistence in farm and processing environments. Overall, S. Kentucky also had comparatively high levels of intra-serovar genetic variability when compared to other prominent serovars from our collection. In addition, one of our two S. Reading isolates had high carriage of both AMR genes and VFs relative to other isolates in our collection. As the poultry-associated Salmonella landscape continues to evolve in Canada, future studies should monitor the genetic composition of prominent serovars across poultry production to maintain up-to-date risk assessments of these foodborne pathogens to consumers.

Genomic Epidemiology and Multilevel Genome Typing of Australian Salmonella enterica Serovar Enteritidis

Salmonella enterica serovar Enteritidis is one of the leading causes of salmonellosis in Australia. In this study, a total of 568 S. Enteritidis isolates from two Australian states across two consecutive years were analyzed and compared to international strains, using the S. Enteritidis multilevel genome typing (MGT) database, which contained 40,390 publicly available genomes from 99 countries. The Australian S. Enteritidis isolates were divided into three phylogenetic clades (A, B, and C). Clades A and C represented 16.4% and 3.5% of the total isolates, respectively, and were of local origin. Clade B accounted for 80.1% of the isolates which belonged to seven previously defined lineages but was dominated by the global epidemic lineage. At the MGT5 level, three out of five top sequence types (STs) in Australia were also top STs in Asia, suggesting that a fair proportion of Australian S. Enteritidis cases may be epidemiologically linked with Asian strains. In 2018, a large egg-associated local outbreak was caused by a recently defined clade B lineage prevalent in Europe and was closely related, but not directly linked, to three European isolates. Additionally, over half (54.8%) of predicted multidrug resistance (MDR) isolates belonged to 10 MDR-associated MGT-STs, which were also frequent in Asian S. Enteritidis . Overall, this study investigated the genomic epidemiology of S. Enteritidis in Australia, including the first large local outbreak, using MGT. The open MGT platform enables a standardized and sharable nomenclature that can be effectively applied to public health for unified surveillance of S. Enteritidis nationally and globally. IMPORTANCE Salmonella enterica serovar Enteritidis is a leading cause of foodborne infections. We previously developed a genomic typing database (MGTdb) for S. Enteritidis to facilitate global surveillance of this pathogen. In this study, we examined the genomic features of Australian S. Enteritidis using the MGTdb and found that Australian S. Enteritidis is mainly epidemiologically linked with Asian strains (especially strains carrying antimicrobial resistance genes), followed by European strains. The first large-scale egg-associated local outbreak in Australia was caused by a recently defined lineage prevalent in Europe, and three European isolates in the MGTdb were closely related but not directly linked to this outbreak. In summary, the S. Enteritidis MGTdb open platform is shown to be a potentially powerful tool for national and global public health surveillance of this pathogen.

EnteroBase: hierarchical clustering of 100 000s of bacterial genomes into species/subspecies and populations

The definition of bacterial species is traditionally a taxonomic issue while bacterial populations are identified by population genetics. These assignments are species specific, and depend on the practitioner. Legacy multilocus sequence typing is commonly used to identify sequence types (STs) and clusters (ST Complexes). However, these approaches are not adequate for the millions of genomic sequences from bacterial pathogens that have been generated since 2012. EnteroBase (http://enterobase.warwick.ac.uk) automatically clusters core genome MLST allelic profiles into hierarchical clusters (HierCC) after assembling annotated draft genomes from short-read sequences. HierCC clusters span core sequence diversity from the species level down to individual transmission chains. Here we evaluate HierCC's ability to correctly assign 100 000s of genomes to the species/subspecies and population levels for Salmonella, Escherichia, Clostridoides, Yersinia, Vibrio and Streptococcus. HierCC assignments were more consistent with maximum-likelihood super-trees of core SNPs or presence/absence of accessory genes than classical taxonomic assignments or 95% ANI. However, neither HierCC nor ANI were uniformly consistent with classical taxonomy of Streptococcus. HierCC was also consistent with legacy eBGs/ST Complexes in Salmonella or Escherichia and with O serogroups in Salmonella. Thus, EnteroBase HierCC supports the automated identification of and assignment to species/subspecies and populations for multiple genera. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Polyphyly in widespread Salmonella enterica serovars and using genomic proximity to choose the best reference genome for bioinformatics analyses

Salmonella is the most common cause of gastroenteritis in the world. Over the past 5 years, whole-genome analysis has led to the high-resolution characterization of clinical and foodborne Salmonella responsible for typhoid fever, foodborne illness or contamination of the agro-food chain. Whole-genome analyses are simplified by the availability of high-quality, complete genomes for mapping analysis and for calculating the pairwise distance between genomes, but unfortunately some difficulties may still remain. For some serovars, the complete genome is not available, or some serovars are polyphyletic and knowing the serovar alone is not sufficient for choosing the most appropriate reference genome. For these serovars, it is essential to identify the genetically closest complete genome to be able to carry out precise genome analyses. In this study, we explored the genomic proximity of 650 genomes of the 58 Salmonella enterica subsp. enterica serovars most frequently isolated in humans and from the food chain in the United States (US) and in Europe (EU), with a special focus on France. For each serovar, to take into account their genomic diversity, we included all the multilocus sequence type (MLST) profiles represented in EnteroBase with 10 or more genomes (on 19 July 2021). A phylogenetic analysis using both core- and pan-genome approaches was carried out to identify the genomic proximity of all the Salmonella studied and 20 polyphyletic serovars that have not yet been described in the literature. This study determined the genetic proximity between all 58 serovars studied and revealed polyphyletic serovars, their genomic lineages and MLST profiles. Finally, we enhanced the open-access databases with 73 new genomes and produced a list of high-quality complete reference genomes for 48 S. enterica subsp. enterica serovars among the most isolated in the US, EU, and France.

Genome-based analysis of infrequent Salmonella serotypes through the Thai pork production chain

Salmonella is a prevalent zoonotic foodborne pathogen. Swine and pork are implicated as important sources of salmonellosis in humans. In Chiang Mai and Lamphun Provinces in northern Thailand, there has been a high prevalence of Salmonella persistence for over a decade. Infection is usually with dominant S. enterica serotypes, including serotypes Rissen and 1,4,[5],12:i:-. However, other serotypes also contribute to disease but are less well characterized. The whole genome sequencing data of 43 S. enterica serotypes isolated from pork production chain through 2011-2014, were used to evaluate genetic diversity and ascertain the possible source of Salmonella contamination based on Core Genome Multilocus Sequence Typing (cgMLST) approach. The Salmonella serotypes recovered from farms and slaughterhouses were re-circulating by swine environmental contamination. Conversely, the Salmonella contamination in the retail market represents cross-contamination from multiple sources, including contaminated foodstuffs. Salmonella contamination in the pork production chain has the competency for host cell adhesion, host cell invasion, and intracellular survival, which is enough for the pathogenicity of salmonellosis. In addition, all of these isolates were multi-drug resistant Salmonella, which contained at least 10 antimicrobial resistance genes. This result indicated that these S. enterica serotypes also pose a significant public health risk. Our findings support the need for appropriate surveillance of food-animal products going to market to reduce public exposure to highly pathogenic, multi-drug resistant Salmonella. Acquiring information would motivate all stakeholders to reinforce sanitation standards throughout the pork production chain in order to eradicate Salmonella contamination and reduce the risk of salmonellosis in humans.

Integration of the Salmonella Typhimurium Methylome and Transcriptome Reveals That DNA Methylation and Transcriptional Regulation Are Largely Decoupled under Virulence-Related Conditions

Despite being in a golden age of bacterial epigenomics, little work has systematically examined the plasticity and functional impacts of the bacterial DNA methylome. Here, we leveraged single-molecule, real-time sequencing (SMRT-seq) to examine the m⁶A DNA methylome of two Salmonella enterica serovar Typhimurium strains: 14028s and a ΔmetJ mutant with derepressed methionine metabolism, grown in Luria broth or medium that simulates the intracellular environment. We found that the methylome is remarkably static: >95% of adenosine bases retain their methylation status across conditions. Integration of methylation with transcriptomic data revealed limited correlation between changes in methylation and gene expression. Further, examination of the transcriptome in ΔyhdJ bacteria lacking the m⁶A methylase with the most dynamic methylation pattern in our data set revealed little evidence of YhdJ-mediated gene regulation. Curiously, despite G(m⁶A)TC motifs being particularly resistant to change across conditions, incorporating dam mutants into our analyses revealed two examples where changes in methylation and transcription may be linked across conditions. This includes the novel finding that the ΔmetJ motility defect may be partially driven by hypermethylation of the chemotaxis gene tsr. Together, these data redefine the S. Typhimurium epigenome as a highly stable system that has rare but important roles in transcriptional regulation. Incorporating these lessons into future studies will be critical as we progress through the epigenomic era.

Phylogenomic Analysis of Salmonella enterica subsp. enterica Serovar Bovismorbificans from Clinical and Food Samples Using Whole Genome Wide Core Genes and kmer Binning Methods to Identify Two Distinct Polyphyletic Genome Pathotypes

Salmonella enterica subsp. enterica serovar Bovismorbificans has caused multiple outbreaks involving the consumption of produce, hummus, and processed meat products worldwide. To elucidate the intra-serovar genomic structure of S. Bovismorbificans, a core-genome analysis with 2690 loci (based on 150 complete genomes representing Salmonella enterica serovars developed as part of this study) and a k-mer-binning based strategy were carried out on 95 whole genome sequencing (WGS) assemblies from Swiss, Canadian, and USA collections of S. Bovismorbificans strains from foodborne infections. Data mining of a digital DNA tiling array of legacy SARA and SARB strains was conducted to identify near-neighbors of S. Bovismorbificans. The core genome analysis and the k-mer-binning methods identified two polyphyletic clusters, each with emerging evolutionary properties. Four STs (2640, 142, 1499, and 377), which constituted the majority of the publicly available WGS datasets from >260 strains analyzed by k-mer-binning based strategy, contained a conserved core genome backbone with a different evolutionary lineage as compared to strains comprising the other cluster (ST150). In addition, the assortment of genotypic features contributing to pathogenesis and persistence, such as antimicrobial resistance, prophage, plasmid, and virulence factor genes, were assessed to understand the emerging characteristics of this serovar that are relevant clinically and for food safety concerns. The phylogenomic profiling of polyphyletic S. Bovismorbificans in this study corresponds to intra-serovar variations observed in S. Napoli and S. Newport serovars using similar high-resolution genomic profiling approaches and contributes to the understanding of the evolution and sequence divergence of foodborne Salmonellae. These intra-serovar differences may have to be thoroughly understood for the accurate classification of foodborne Salmonella strains needed for the uniform development of future food safety mitigation strategies.

Through the Looking Glass: Genome, Phenome, and Interactome of Salmonella enterica

This review revisits previous concepts on biological phenomenon contributing to the success of the Salmonella enterica subspecies I as a pathogen and expands upon them to include progress in epidemiology based on whole genome sequencing (WGS). Discussion goes beyond epidemiological uses of WGS to consider how phenotype, which is the biological character of an organism, can be correlated with its genotype to develop a knowledge of the interactome. Deciphering genome interactions with proteins, the impact of metabolic flux, epigenetic modifications, and other complex biochemical processes will lead to new therapeutics, control measures, environmental remediations, and improved design of vaccines.

Salmonella enterica subsp. enterica Welikade: guideline for phylogenetic analysis of serovars rarely involved in foodborne outbreaks

BACKGROUND

Salmonella spp. is a major foodborne pathogen with a wide variety of serovars associated with human cases and food sources. Nevertheless, in Europe a panel of ten serovars is responsible for up to 80% of confirmed human cases. Clustering studies by single nucleotide polymorphism (SNP) core-genome phylogenetic analysis of outbreaks due to these major serovars are simplified by the availability of many complete genomes in the free access databases. This is not the case for outbreaks due to less common serovars, such as Welikade, for which no reference genomes are available. In this study, we propose a method to solve this problem. We propose to perform a core genome MLST (cgMLST) analysis based on hierarchical clustering using the free-access EnteroBase to select the most suitable genome to use as a reference for SNP phylogenetic analysis. In this study, we applied this protocol to a retrospective analysis of a Salmonella enterica serovar Welikade (S. Welikade) foodborne outbreak that occurred in France in 2016. Finally, we compared the cgMLST and SNP analyses. SNP phylogenetic reconstruction was carried out considering the effect of recombination events identified by the ClonalFrameML tool. The accessory genome was also explored by phage content and virulome analyses.

RESULTS

Our findings revealed high clustering concordance using cgMLST and SNP analyses. Nevertheless, SNP analysis allowed for better assessment of the genetic distance among strains. The results revealed epidemic clones of S. Welikade circulating within the poultry and dairy sectors in France, responsible for sporadic and non-sporadic human cases between 2012 and 2019.

CONCLUSIONS

This study increases knowledge on this poorly described serovar and enriches public genome databases with 42 genomes from human and non-human S. Welikade strains, including the isolate collected in 1956 in Sri Lanka, which gave the name to this serovar. This is the first genomic analysis of an outbreak due to S. Welikade described to date.

Investigating the dynamics of Salmonella contamination in integrated poultry companies using a whole genome sequencing approach

The study of non-typhoid Salmonella in broiler integrations has been limited by the resolution of typing techniques. Although serotyping of Salmonella isolates is used as a traditional approach, it is not of enough resolution to clearly understand the dynamics of this pathogen within poultry companies. The aim of this research was to investigate the epidemiology and population dynamics of Salmonella serotypes in 2 poultry integrations using a whole genome sequencing approach. Two hundred and forty-three Salmonella isolates recovered from the broiler production chain of 2 integrated poultry companies were whole genome sequenced and analyzed with dedicated databases and bioinformatic software. The analyses of sequences revealed that S. Infantis was the most frequent serotype (82.3%). Most isolates showed a potential for resistance against medically important antibiotics and disinfectants. Furthermore, 97.5% of isolates harbored the pESI-like mega plasmid, that plays an important role in the global dissemination of AMR. SNP tree analysis showed that there were clones that are niche-specific while other ones were distributed throughout the broiler production chains. In this study, we demonstrated the potential of whole genome sequencing analysis for a comprehensive understanding of Salmonella distribution in integrated poultry companies. Data obtained with these techniques allow determination of the presence of genetic factors that play an important role in the environmental fitness and pathogenicity of Salmonella.

Identification of a New Serovar of Salmonella enterica in Mediterranean Buffalo Calves (Bubalus bubalis)

Simple Summary

The presence of Salmonella enterica in buffalo (Bubalus bubalis) farms is a cause for great concern, as it is responsible for serious economic losses and is a zoonotic agent capable of contaminating the mozzarella production chain. Salmonella induced disease mainly affects buffalo calves and is characterized by severe gastrointestinal lesions, profuse diarrhea, severe dehydration, septicemia and death; in adults it can also induce termination of pregnancy. The main source of infection in the herd is asymptomatic elderly animals that release high bacterial charges through the feces. Other sources of infection are contaminated fodder and water; the infection can also be transmitted by rodents, wild birds, insects and humans. There are several serotypes of Salmonella reported in the literature capable of infecting the buffalo. This work describes an important outbreak of Salmonellosis that occurred in young subjects of a buffalo farm located in the province of Salerno. It is supported by a new serotype of Salmonella so far never described and included in the Kauffmann–White classification.

Abstract

This case report describes for the first-time cases of severe gastroenteritis in water buffalo calves due to a new serovar of Salmonella enterica. The study was carried out on fecal matrix collected from live water buffalo calves that showed profuse diarrhea, severe dehydration and fever, exhibiting a systemic course. Culture and molecular investigations identified the pathogens isolated from intestinal contents as two Salmonella serovars, Salmonella enterica enterica O:35 and a new serovar of Salmonella enterica. The isolates showed multi-drug resistance. Timely diagnosis associated with a targeted antimicrobial treatment were found to be sufficient for the survival and recovery of the infected animals. Herd vaccines prepared from isolated pathogens were used to prevent further deaths of the calves.

An accessible, efficient and global approach for the large-scale sequencing of bacterial genomes

We have developed an efficient and inexpensive pipeline for streamlining large-scale collection and genome sequencing of bacterial isolates. Evaluation of this method involved a worldwide research collaboration focused on the model organism Salmonella enterica, the 10KSG consortium. Following the optimization of a logistics pipeline that involved shipping isolates as thermolysates in ambient conditions, the project assembled a diverse collection of 10,419 isolates from low- and middle-income countries. The genomes were sequenced using the LITE pipeline for library construction, with a total reagent cost of less than USD$10 per genome. Our method can be applied to other large bacterial collections to underpin global collaborations.

A Review of Significant European Foodborne Outbreaks in the Last Decade

ABSTRACT

Foodborne diseases remain a global public health challenge worldwide. The European surveillance system of multistate foodborne outbreaks integrates elements from public and animal health and the food chain for early detection, assessment, and control. This review includes descriptions of the significant outbreaks that occurred in Europe in the last decade. Their significance and relevance to public health is derived from the changes, improvements, and novelties that pushed toward building a safer food system in the European Union, certainly driven by the One Health approach. In 2011, a point source monoclonal outbreak of infections caused by Escherichia coli serotype O104:H4 in sprouted seeds resulted in hundreds of cases of hemolytic uremic syndrome and several fatalities. In 2015, a prolonged outbreak of Listeria monocytogenes infections caused by contamination of frozen corn in Europe resulted in 47 cases and nine deaths. In 2016, a persistent polyclonal outbreak of Salmonella Enteritidis was linked to the consumption of eggs and was associated with hundreds of cases. The outbreak evaluations highlight the importance of rapid sharing of data (e.g., sequencing and tracing data) and the need for harmonizing bioinformatics outputs and computational approaches to facilitate detection and investigation of foodborne illnesses. These outbreaks led to development of a legal framework for a European collaboration platform for sharing whole genome sequence data and enabled the enforcement of existing hygiene and food safety provisions and the development of new hygiene guidelines and best practices. This review also briefly touches on the new trends in information technologies that are being explored for food traceability and safety. These technologies could enhance the traceability of food throughout the supply chain and redirect the conventional tracing system toward a digitized supply chain.

HIGHLIGHTS

The Invasin and Complement-Resistance Protein Rck of Salmonella is More Widely Distributed than Previously Expected

The rck open reading frame (ORF) on the pefI-srgC operon encodes an outer membrane protein responsible for invasion of nonphagocytic cell lines and resistance to complement-mediated killing. Until now, the rck ORF was only detected on the virulence plasmids of three serovars of Salmonella subsp. enterica (i.e., Bovismorbificans, Enteritidis, and Typhimurium). The increasing number of Salmonella genome sequences allowed us to use a combination of reference sequences and whole-genome multilocus sequence typing (wgMLST) data analysis to probe the presence of the operon and of rck in a wide array of isolates belonging to all Salmonella species and subspecies. We established the presence of partial or complete operons in 61 subsp. enterica serovars as well as in 4 other subspecies with various syntenies and frequencies. The rck ORF itself was retrieved in 36 subsp. enterica serovars and in two subspecies with either chromosomal or plasmid-borne localization. It displays high conservation of its sequence within the genus, and we demonstrated that most of the allelic variations identified did not alter the virulence properties of the protein. However, we demonstrated the importance of the residue at position 38 (at the level of the first extracellular loop of the protein) in the invasin function of Rck. Altogether, our results highlight that rck is not restricted to the three formerly identified serovars and could therefore have a more important role in virulence than previously expected. Moreover, this work raises questions about the mechanisms involved in rck acquisition and about virulence plasmid distribution and evolution.

IMPORTANCE The foodborne pathogen Salmonella is responsible for a wide variety of pathologies depending on the infected host, the infecting serovars, and its set of virulence factors. However, the implication of each of these virulence factors and their role in the specific host-pathogen interplay are not fully understood. The significance of our research is in determining the distribution of one of these factors, the virulence plasmid-encoded invasin and resistance to complement killing protein Rck. In addition to providing elements of reflection concerning the mechanisms of acquisition of specific virulence genes in certain serotypes, this work will help to understand the role of Rck in the pathogenesis of Salmonella.

AT Homopolymer Strings in Salmonella enterica Subspecies I Contribute to Speciation and Serovar Diversity

Adenine and thymine homopolymer strings of at least 8 nucleotides (AT 8+mers) were characterized in Salmonella enterica subspecies I. The motif differed between other taxonomic classes but not between Salmonella enterica serovars. The motif in plasmids was possibly associated with serovar. Approximately 12.3% of the S. enterica motif loci had mutations. Mutability of AT 8+mers suggests that genomes undergo frequent repair to maintain optimal gene content, and that the motif facilitates self-recognition; in addition, serovar diversity is associated with plasmid content. A theory that genome regeneration accounts for both persistence of predominant Salmonella serovars and serovar diversity provides a new framework for investigating root causes of foodborne illness.

Genomic population structure associated with repeated escape of Salmonella enterica ATCC14028s from the laboratory into nature

Salmonella enterica serovar Typhimurium strain ATCC14028s is commercially available from multiple national type culture collections, and has been widely used since 1960 for quality control of growth media and experiments on fitness (“laboratory evolution”). ATCC14028s has been implicated in multiple cross-contaminations in the laboratory, and has also caused multiple laboratory infections and one known attempt at bioterrorism. According to hierarchical clustering of 3002 core gene sequences, ATCC14028s belongs to HierCC cluster HC20_373 in which most internal branch lengths are only one to three SNPs long. Many natural Typhimurium isolates from humans, domesticated animals and the environment also belong to HC20_373, and their core genomes are almost indistinguishable from those of laboratory strains. These natural isolates have infected humans in Ireland and Taiwan for decades, and are common in the British Isles as well as the Americas. The isolation history of some of the natural isolates confirms the conclusion that they do not represent recent contamination by the laboratory strain, and 10% carry plasmids or bacteriophages which have been acquired in nature by HGT from unrelated bacteria. We propose that ATCC14028s has repeatedly escaped from the laboratory environment into nature via laboratory accidents or infections, but the escaped micro-lineages have only a limited life span. As a result, there is a genetic gap separating HC20_373 from its closest natural relatives due to a divergence between them in the late 19^th century followed by repeated extinction events of escaped HC20_373.

Clades of closely related bacteria exist in nature. Individual isolates from such clades are often distinguishable by genomic sequencing because genomic sequence differences can be acquired over a few years due to neutral drift and natural selection. The evolution of laboratory strains is often largely frozen, physically due to storage conditions and genetically due to long periods of storage. Thus, laboratory strains can normally be readily distinguished from natural isolates because they show much less diversity. However, laboratory strain ATCC14028s shows modest levels of sequence diversity because it has been shipped around the world to multiple laboratories and is routinely used for analyses of laboratory evolution. Closely related natural isolates also exist, but their genetic diversity is not dramatically greater at the core genome level. Indeed, many scientists doubt that such isolates are natural, and interpret them as undetected contamination by the laboratory strain. We present data indicating that ATCC14028s has repeatedly escaped from the laboratory through inadvertent contamination of the environment, infection of technical staff and deliberate bioterrorism. The escapees survive in nature long enough that some acquire mobile genomic elements by horizontal gene transfer, but eventually they go extinct. As a result, even extensive global databases of natural isolates lack closely related isolates whose ancestors diverged from ATCC14028s within the last 100 years.

Prevalence, Serotypes, and Antimicrobial Resistance of Salmonella from Mink Feces and Feed in the United States

Nontyphoidal Salmonella is a leading cause of foodborne illnesses, and concentrated animal production such as commercial mink farming can be a reservoir. The objective of this study was to determine the prevalence, serotypes, antimicrobial resistance, virulence genes and mobile genetic elements, and phylogenetic characterization of Salmonella from mink fecal (n = 42) and feed (n = 8) samples at a commercial mink farm in the United States. Salmonella was detected from 11 (26.2%) fecal and 1 (12.5%) feed samples. Twelve fecal isolates obtained from the positive samples were identified into four serotypes: Salmonella Uganda (eight isolates; two isolates were obtained from one sample), Salmonella Reading (two isolates), Salmonella Heidelberg (one isolate), and Salmonella Orion (one isolate). Two isolates from the positive feed sample were identified as Salmonella Cerro. Twelve isolates, one isolate per positive sample, were whole genome sequenced for further characterization. The two Salmonella Reading isolates were multidrug resistant (MDR) with the classical ASSuT (ampicillin, streptomycin, sulfisoxazole, and tetracycline) phenotype. Genotypically, the isolates were correspondingly positive, both by polymerase chain reaction and by whole-genome sequencing (WGS), for bla_TEM-1, aph(3″)-Ib and aph(6)-Id, sul2, and tet(A). WGS additionally identified trimethoprim resistance gene dfrA8 in Salmonella Reading, and fosfomycin resistance gene fosA7 in Salmonella Heidelberg. All isolates carried Salmonella-specific multidrug resistance (mdsA/B), gold resistance, and intimin-like virulence genes. In addition, Salmonella Orion carried copper and silver resistance gene clusters. The two Salmonella Reading isolates also carried a cytolethal distension toxin (cdtB) gene. Salmonella Reading isolates belonged to ST412, a predominant sequence type among turkey and human isolates obtained in relation to recent salmonellosis outbreaks in North America. Phylogenetically, Salmonella Uganda, Salmonella Heidelberg, and Salmonella Reading were mostly associated with historic human/clinical, and poultry and swine source isolates. This study indicates that mink production can be a reservoir for foodborne pathogens such as Salmonella with MDR serovars commonly associated with human illnesses.

Investigation of Salmonella Phage-Bacteria Infection Profiles: Network Structure Reveals a Gradient of Target-Range from Generalist to Specialist Phage Clones in Nested Subsets

Bacteriophages that lyse Salmonella enterica are potential tools to target and control Salmonella infections. Investigating the host range of Salmonella phages is a key to understand their impact on bacterial ecology, coevolution and inform their use in intervention strategies. Virus-host infection networks have been used to characterize the "predator-prey" interactions between phages and bacteria and provide insights into host range and specificity. Here, we characterize the target-range and infection profiles of 13 Salmonella phage clones against a diverse set of 141 Salmonella strains. The environmental source and taxonomy contributed to the observed infection profiles, and genetically proximal phages shared similar infection profiles. Using in vitro infection data, we analyzed the structure of the Salmonella phage-bacteria infection network. The network has a non-random nested organization and weak modularity suggesting a gradient of target-range from generalist to specialist species with nested subsets, which are also observed within and across the different phage infection profile groups. Our results have implications for our understanding of the coevolutionary mechanisms shaping the ecological interactions between Salmonella phages and their bacterial hosts and can inform strategies for targeting Salmonella enterica with specific phage preparations.

High-Resolution Genomic Comparisons within Salmonella enterica Serotypes Derived from Beef Feedlot Cattle: Parsing the Roles of Cattle Source, Pen, Animal, Sample Type, and Production Period

Salmonella enterica is a major foodborne pathogen, and contaminated beef products have been identified as one of the primary sources of Salmonella-related outbreaks. Pathogenicity and antibiotic resistance of Salmonella are highly serotype and subpopulation specific, which makes it essential to understand high-resolution Salmonella population dynamics in cattle. Time of year, source of cattle, pen, and sample type (i.e., feces, hide, or lymph nodes) have previously been identified as important factors influencing the serotype distribution of Salmonella (e.g., Anatum, Lubbock, Cerro, Montevideo, Kentucky, Newport, and Norwich) that were isolated from a longitudinal sampling design in a research feedlot. In this study, we performed high-resolution genomic comparisons of Salmonella isolates within each serotype using both single-nucleotide polymorphism-based maximum-likelihood phylogeny and hierarchical clustering of core-genome multilocus sequence typing. The importance of the aforementioned features in clonal Salmonella expansion was further explored using a supervised machine learning algorithm. In addition, we identified and compared the resistance genes, plasmids, and pathogenicity island profiles of the isolates within each subpopulation. Our findings indicate that clonal expansion of Salmonella strains in cattle was mainly influenced by the randomization of block and pen, as well as the origin/source of the cattle, i.e., regardless of sampling time and sample type (i.e., feces, lymph node, or hide). Further research is needed concerning the role of the feedlot pen environment prior to cattle placement to better understand carryover contributions of existing strains of Salmonella and their bacteriophages. IMPORTANCE Salmonella serotypes isolated from outbreaks in humans can also be found in beef cattle and feedlots. Virulence factors and antibiotic resistance are among the primary defense mechanisms of Salmonella, and are often associated with clonal expansion. This makes understanding the subpopulation dynamics of Salmonella in cattle critical for effective mitigation. There remains a gap in the literature concerning subpopulation dynamics within Salmonella serotypes in feedlot cattle from the beginning of feeding up until slaughter. Here, we explore Salmonella population dynamics within each serotype using core-genome phylogeny and hierarchical classifications. We used machine learning to quantitatively parse the relative importance of both hierarchical and longitudinal clustering among cattle host samples. Our results reveal that Salmonella populations in cattle are highly clonal over a 6-month study period and that clonal dissemination of Salmonella in cattle is mainly influenced spatially by experimental block and pen, as well by the geographical origin of the cattle.

Genomic Analysis of Antibiotic-Resistant and -Susceptible Escherichia coli Isolated from Bovine Sources in Maputo, Mozambique

This study reports a genomic analysis of Escherichia coli isolates recovered from 25 bovine fecal composite samples collected from four different production units in Maputo city and around Maputo Province, Mozambique. The genomes were analyzed to determine the presence of antibiotic resistance genes (ARGs), genetic relatedness, and virulence factors known to cause diseases in humans. Whole-genome sequencing was conducted on 28 isolates using an Illumina NextSeq 500 sequencing platform. The genomes were analyzed using BLASTN for the presence of resistance genes and virulence factors, as well as to determine their phylogenetic groups, sequence types (ST), and ST complexes (ST Cplxs). The majority of the isolates (85%) were identified as members of phylogenetic groups B1, with fewer isolates identified as members of group A, and a single isolate identified as group "E/Clade I." The ST analysis demonstrated a higher level of diversity than the phylogenetic group analysis. Sixteen different STs, five ST Cplxs, and seven singleton complexes were identified. A strain identified as a novel ST (ST9215) showed a high level of similarity with an isolate recovered from a wild animal in the Gambia. Seven different ARGs were identified, with tet(B) being the most frequently detected, followed by aph(3″)-Ib, aph(6)-Id, sul2, bla_TEM-1B, and dfrA1. Three isolates encoded β-lactam-conferring point mutations in the ampC promoter (-42C>T). In total, 51 different virulence factors were identified among the genomes. This study demonstrates that E. coli from bovine sources in Mozambique encoded multiple antibiotic resistance elements, plasmids, and virulence factors. To the best of our knowledge, this is the first genomic description of antibiotic-resistant E. coli isolated from bovine sources in Mozambique.