Genomic evidence of environmental and resident Salmonella Senftenberg and Montevideo contamination in the pistachio supply-chain

Genomic characteristics of Salmonella Montevideo and Pomona: impact of isolation source on antibiotic resistance, virulence and metabolic capacity

Salmonella enterica is known for its disease-causing serotypes, including Montevideo and Pomona. These serotypes have been found in various environments, including river water, sediments, food, and animals. However, the global spread of these serotypes has increased, leading to many reported infections and outbreaks. The goal of this study was the genomic analysis of 48 strains of S. Montevideo and S. Pomona isolated from different sources, including clinical. Results showed that environmental strains carried more antibiotic resistance genes than the clinical strains, such as genes for resistance to aminoglycosides, chloramphenicol, and sulfonamides. Additionally, the type 4 secretion system, was only found in environmental strains. .Also many phosphotransferase transport systems were identified and the presence of genes for the alternative pathway Entner-Doudoroff. The origin of isolation may have a significant impact on the ability of Salmonella isolates to adapt and survive in different environments, leading to genomic flexibility and a selection advantage.

Phenotypic Characteristics That May Contribute to Persistence of Salmonella Strains in the Pistachio Supply Chain

Salmonella enterica subsp. enterica strain diversity in California pistachios is limited; some strains have persisted in the pistachio supply chain for ≥10 years. Representative isolates of six persistent strains and three sporadic strains isolated from California pistachios were selected to evaluate copper resistance, growth in pistachio hull slurry, biofilm formation, desiccation tolerance, and survival during subsequent storage. The presence of a copper homeostasis and silver-resistance island sequence in three of the persistent strains was associated with an increase in tolerance to CuSO4 from 7.5 mM to 15 mM under anaerobic but not aerobic conditions; all isolates were resistant to ≥120 mM Cu-EDTA under both anerobic and aerobic conditions. When inoculated into pistachio hull slurry at 2.75 ± 0.04 log CFU/mL and incubated at 30 °C, the populations of Salmonella Enteritidis strain A (sporadic) increased to significantly lower levels than the other strains at 16, 20, 24, and 28 h but not at 40 and 48 h. Maximum populations of 8.70-8.85 log CFU/mL were observed for all strains at ≥40 h of incubation. All nine Salmonella strains produced weak to strong biofilms after 4 days at 25 °C; seven strains, including two sporadic strains, produced moderate biofilms, and Salmonella Liverpool strain A (persistent) produced a strong biofilm. The rdar+ and rdar- morphotypes were observed in both persistent and sporadic Salmonella strains. Population declines of 5.03 log were observed for Salmonella Enteritidis strain A within 18 h of drying on filter paper whereas reductions of 0.50-1.25 log were observed for the other eight Salmonella strains. Population reductions (3.98-5.12 log) of these eight strains were not significantly different after storage at 25 ± 1 °C and 35% relative humidity for 50 days. The phenotypic characteristics evaluated here do not independently account for the persistence of a small number of Salmonella strains associated with the California pistachio production chain.

Glycogen plays a key role in survival of Salmonella Typhimurium on dry surfaces and in low-moisture foods

Salmonella enterica is known to survive in desiccate environments and is often associated with low-moisture foods (LMFs). In this work, S. Typhimurium ATCC 14028 was found to survive better by achieving the least reductions (3.17 ± 0.20 Log CFU reduction) compared to S. Tennessee ATCC 10722 (3.82 ± 0.13 Log CFU reduction) and S. Newport ATCC 6962 (6.03 ± 0.36 Log CFU reduction) after 30 days on surfaces with a relative humidity of 49% at ambient temperature. A metabolomic analysis revealed that S. Typhimurium was still active in energy metabolism after 24 h in the desiccate environment and glycogen, an energy reserve, was drastically reduced. We followed up on the glycogen levels over 30 days and found indeed a sharp decline on the first day. However, the glycogens detected on day 7 were significantly higher (P < 0.05) and thereafter remained stable above the original levels until day 30. The expression levels of both glycogen anabolism- and catabolism-related genes (csrA, glgA, glgC, glgX) were significantly up-regulated at all tested points (P < 0.05). The glgA and glgC insertion mutants displayed weaker survivability on both dry surfaces and in representative LMFs (flour and milk powder) compared to the wild-type strain. This work highlights the role of glycogen during different periods of desiccation, which may bring novel insight into mitigating Salmonella by disrupting glycogen metabolism.

Characterizing the Genetic Diversity of Salmonella Isolated from U.S. Raw Inshell Pistachios Using Whole Genome Sequencing

The genetic diversity of 169 Salmonella isolates from pistachios collected from California storage silos during the 2010, 2011, and 2012 harvests (silo survey isolates) was determined by analyzing the whole genome sequence data using the CFSAN SNP pipeline developed by the U.S. Food and Drug Administration's Center for Food Safety and Applied Nutrition. Salmonella isolates clustered by serovars Agona, Enteritidis, Montevideo, Sandiego, Senftenberg, Liverpool, Tennessee, and Worthington in the phylogenetic tree. Within each serovar, isolates grouped into one or two clusters (≤14 SNPs). Two distinct clusters (>14 SNPs; A and B) were identified for Salmonella Enteritidis, Montevideo, and Liverpool for a total of 11 unique strains. Sequences of representative silo survey isolates clustered with sequences of Salmonella strains isolated from U.S. pistachio-associated samples collected between 2008 and 2018 available on the National Center for Biotechnology Information database, and, in all but two cases, not with sequences of Salmonella strains recovered from raw California almonds from 2001 through 2013. The genomic evidence suggests that strains of Salmonella Agona, Liverpool Cluster A, Montevideo Clusters A and B, Senftenberg, and Worthington have persisted in the California pistachio environment for ≥3 years and some of these strains have been reported exclusively in association with pistachios.

Genetic diversity of Salmonella enterica isolated over 13 years from raw California almonds and from an almond orchard

A comparative genomic analysis was conducted for 171 Salmonella isolates recovered from raw inshell almonds and raw almond kernels between 2001 and 2013 and for 30 Salmonella Enteritidis phage type (PT) 30 isolates recovered between 2001 and 2006 from a 2001 salmonellosis outbreak-associated almond orchard. Whole genome sequencing was used to measure the genetic distance among isolates by single nucleotide polymorphism (SNP) analyses and to predict the presence of plasmid DNA and of antimicrobial resistance (AMR) and virulence genes. Isolates were classified by serovars with Parsnp, a fast core-genome multi aligner, before being analyzed with the CFSAN SNP Pipeline (U.S. Food and Drug Administration Center for Food Safety and Applied Nutrition). Genetically similar (≤18 SNPs) Salmonella isolates were identified among several serovars isolated years apart. Almond isolates of Salmonella Montevideo (2001 to 2013) and Salmonella Newport (2003 to 2010) differed by ≤9 SNPs. Salmonella Enteritidis PT 30 isolated between 2001 and 2013 from survey, orchard, outbreak, and clinical samples differed by ≤18 SNPs. One to seven plasmids were found in 106 (62%) of the Salmonella isolates. Of the 27 plasmid families that were identified, IncFII and IncFIB plasmids were the most predominant. AMR genes were identified in 16 (9%) of the survey isolates and were plasmid encoded in 11 of 16 cases; 12 isolates (7%) had putative resistance to at least one antibiotic in three or more drug classes. A total of 303 virulence genes were detected among the assembled genomes; a plasmid that harbored a combination of pef, rck, and spv virulence genes was identified in 23% of the isolates. These data provide evidence of long-term survival (years) of Salmonella in agricultural environments.

Evaluation of multiplex nanopore sequencing for Salmonella serotype prediction and antimicrobial resistance gene and virulence gene detection

In a previous study, Multiplex-nanopore-sequencing based whole genome sequencing (WGS) allowed for accurate in silico serotype prediction of Salmonella within one day for five multiplexed isolates, using both SISTR and SeqSero2. Since only ten serotypes were tested in our previous study, the conclusions above were yet to be evaluated in a larger scale test. In the current study we evaluated this workflow with 69 Salmonella serotypes and also explored the feasibility of using multiplex-nanopore-sequencing based WGS for antimicrobial resistance gene (AMR) and virulence gene detection. We found that accurate in silico serotype prediction with nanopore-WGS data was achieved within about five hours of sequencing at a minimum of 30× Salmonella genome coverage, with SeqSero2 as the serotype prediction tool. For each tested isolate, small variations were observed between the AMR/virulence gene profiles from the Illumina and Nanopore sequencing platforms. Taking results generated using Illumina data as the benchmark, the average precision value per isolate was 0.99 for both AMR and virulence gene detection. We found that the resistance gene identifier - RGI identified AMR genes with nanopore data at a much lower accuracy compared to Abricate, possibly due to RGI's less stringent minimum similarity and coverage by default for database matching. This study is an evaluation of multiplex-nanopore-sequencing based WGS as a cost-efficient and rapid Salmonella classification method, and a starting point for future validation and verification of using it as a AMR/virulence gene profiling tool for the food industry. This study paves the way for the application of nanopore sequencing in surveillance, tracking, and risk assessment of Salmonella across the food supply chain.

Metabolomics-based response of Salmonella to desiccation stress and skimmed milk powder storage

The strong survival ability of Salmonella in low-moisture foods (LMFs) has been of public concern, and is considered a threat to people's health. Recently, the development of omics technology has promoted research on the molecular mechanisms of the desiccation stress response of pathogenic bacteria. However, multiple analytical aspects related to their physiological characteristics remain unclear. We explored the physiological metabolism changes of S. enterica Enteritidis exposed to a 24 h-desiccation treatment and a subsequent 3-month desiccation storage in skimmed milk powder (SMP) with an approach of gas chromatography-mass spectrometry (GC-MS) and ultra-performance liquid chromatography-Q Exactive-mass spectrometry (UPLC-QE-MS). A total of 8,292 peaks were extracted, of which 381 were detected by GC-MS and 7,911 peaks were identified by LC-MS/MS, respectively. Through analyses of differentially expressed metabolites (DEMs) and key pathways, a total of 58 DEMs emerged from the 24 h-desiccation treatment, which exhibited the highest relevance for five metabolic pathways, involving glycine, serine, and threonine metabolism, pyrimidine metabolism, purine metabolism, vitamin B6 metabolism, and pentose phosphate pathway. After 3-month SMP storage, 120 DEMs were identified, which were related to several regulatory pathways including arginine and proline metabolism, serine and threonine metabolism, β-alanine metabolism, glycerolipid metabolism, and glycolysis. The analyses of key enzyme activities of XOD, PK, and G6PDH and ATP content provided further evidence that supported the metabolic responses such as nucleic acid degradation, glycolysis, and ATP production played an important role in Salmonella's adaptation to desiccation stress. This study enables a better understanding of metabolomics-based responses of Salmonella at the initial stage of desiccation stress and the following long-term adaptive stage. Meanwhile, the identified discriminative metabolic pathways may serve as potentially useful targets in developing strategies for the control and prevention of desiccation-adapted Salmonella in LMFs.

Antimicrobial resistance and genomic characterization of Salmonella enterica serovar Senftenberg isolates in production animals from the United States

In the USA, Salmonella enterica subspecies enterica serovar Senftenberg is among the top five serovars isolated from food and the top 11 serovars isolated from clinically ill animals. Human infections are associated with exposure to farm environments or contaminated food. The objective of this study was to characterize S. Senftenberg isolates from production animals by analyzing phenotypic antimicrobial resistance profiles, genomic features and phylogeny. Salmonella Senftenberg isolates (n = 94) from 20 US states were selected from NVSL submissions (2014-2017), tested against 14 antimicrobial drugs, and resistance phenotypes determined. Resistance genotypes were determined using whole genome sequencing analysis with AMRFinder and the NCBI and ResFinder databases with ABRicate. Plasmids were detected using PlasmidFinder. Integrons were detected using IntFinder and manual alignment with reference genes. Multilocus-sequence-typing (MLST) was determined using ABRicate with PubMLST database, and phylogeny was determined using vSNP. Among 94 isolates, 60.6% were resistant to at least one antimicrobial and 39.4% showed multidrug resistance. The most prevalent resistance findings were for streptomycin (44.7%), tetracycline (42.6%), ampicillin (36.2%) and sulfisoxazole (32.9%). The most commonly found antimicrobial resistance genes were aac(6')-Iaa (100%), aph(3″)-Ib and aph(6)-Id (29.8%) for aminoglycosides, followed by bla _TEM-1 (26.6%) for penicillins, sul1 (25.5%) and sul2 (23.4%) for sulfonamides and tetA (23.4%) for tetracyclines. Quinolone-resistant isolates presented mutations in gyrA and/or parC genes. Class 1 integrons were found in 37 isolates. Thirty-six plasmid types were identified among 77.7% of the isolates. Phylogenetic analysis identified two distinct lineages of S. Senftenberg that correlated with the MLST results. Isolates were classified into two distinct sequence types (ST): ST14 (97.9%) and ST 185 (2.1%). The diversity of this serotype suggests multiple introductions into animal populations from outside sources. This study provided antimicrobial susceptibility and genomic characteristics of S. Senftenberg clinical isolates from production animals in the USA during 2014 to 2017. This study will serve as a base for future studies focused on the phenotypic and molecular antimicrobial characterization of S. Senftenberg isolates in animals. Monitoring of antimicrobial resistance to detect emergence of multidrug-resistant strains is critical.

Changes in the genomes and methylomes of three Salmonella enterica serovars after long-term storage in ground black pepper

Low moisture foods (LMFs) have traditionally been recognized as safe for consumption, as most bacteria require higher water content to grow. However, outbreaks due to LMF foods are increasing, and the microbial pathogen Salmonella enterica is frequently implicated. S. enterica can survive in LMFs for years, but few serovars have been studied, and the mechanisms which underlie this longevity are not well understood. Here, we determine that S. enterica serovars S. Tennessee, S. Anatum, and S. Reading but not S. Oranienburg can survive in the ground black pepper for 6 years. S. Reading was not previously associated with any LMF. Using both Illumina and Pacific Biosciences sequencing technologies, we also document changes in the genomes and methylomes of the surviving serovars over this 6-year period. The three serovars acquired a small number of single nucleotide polymorphisms (SNPs) including seven substitutions (four synonymous, two non-synonymous, and one substitution in a non-coding region), and two insertion-deletions. Nine distinct N6-methyladenine (m6A) methylated motifs across the three serovars were identified including five which were previously known, G^m6ATC, CAG^m6AG, BATGC^m6AT, CRT^m6AYN⁶CTC, and CC^m6AN⁷TGAG, and four novel serovar-specific motifs, GRT^m6AN⁸TTYG, GA^m6ACN⁷GTA, GAA m^6ACY, and CAA^m6ANCC. Interestingly, the BATGCAT motif was incompletely methylated (35–64% sites across the genome methylated), suggesting a possible role in gene regulation. Furthermore, the number of methylated BATGC^m6AT motifs increased after storage in ground black pepper for 6 years from 475 to 657 (S. Tennessee), 366 to 608 (S. Anatum), and 525 to 570 (S. Reading), thus warranting further study as an adaptive mechanism. This is the first long-term assessment of genomic changes in S. enterica in a low moisture environment, and the first study to examine the methylome of any bacteria over a period of years, to our knowledge. These data contribute to our understanding of S. enterica survival in LMFs, and coupled with further studies, will provide the information necessary to design effective interventions which reduce S. enterica in LMFs and maintain a healthy, safe food supply.