Phenotypic and genomic characterizations of Klebsiella pneumoniae ssp. pneumoniae and Rahnella inusitata strains reveal no clear association between genetic content and ropy phenotype.

Ropy defect of pasteurized fluid milk is a type of spoilage which manifests itself by an increased viscosity, slimy body, and string-like flow during pouring. This defect has, among other causes, been attributed to the growth, proliferation and exopolysaccharide production by coliform bacteria, which are most commonly introduced in milk as post-pasteurization contaminants. As we identified both Klebsiella pneumoniae ssp. pneumoniae and Rahnella inusitata that were linked to a ropy defect, the goal of this study was to characterize 3 K. pneumoniae ssp. pneumoniae strains and 2 R. inusitata for (1) their ability to grow and cause ropy defect in milk at 6°C and 21°C and to (2) probe the genetic basis for observed ropy phenotype. Although all K. pneumoniae ssp. pneumoniae and R. inusitata strains showed net growth of >4 log10 over 48 h in UHT milk at 21°C, only R. inusitata strains displayed growth during 28-d incubation period at 6°C (>6 log10). Two out of 3 K. pneumoniae ssp. pneumoniae strains were capable of causing the ropy defect in milk at 21°C, as supported by an increase in the viscosity of milk and string-like flow during pouring; these 2 strains were originally isolated from raw milk. Only one R. inusitata strains was able to cause the ropy defect in milk; this strain was able to cause the defect at both 6°C and 21°C, and was originally isolated from a pasteurized milk. These findings suggest that the potential of K. pneumoniae ssp. pneumoniae and R. inusitata to cause ropy defect in milk is a strain-dependent characteristic. Comparative genomics provided no definitive answer on genetic basis for the ropy phenotype. However, for K. pneumoniae ssp. pneumoniae, genes rffG, rffH, rfbD, and rfbC involved in biosynthesis and secretion of enterobacterial common antigen (ECA) could only be found in the 2 strains that produced ropy defect, and for R. inusitata a set of 2 glycosyltransferase- and flippase genes involved in nucleotide sugar biosynthesis and export could only be identified in the ropy strain. Although these results provide some initial information for potential markers for strains that can cause ropy milk, the relationship between genetic content and ropiness in milk remains poorly understood and merits further investigation.

Whole-Genome Sequencing-Based Characterization of Listeria Isolates from Produce Packinghouses and Fresh-Cut Facilities Suggests Both Persistence and Reintroduction of Fully Virulent L. monocytogenes.

The contamination of ready-to-eat produce with Listeria monocytogenes (LM) can often be traced back to environmental sources in processing facilities and packinghouses. To provide an improved understanding of Listeria sources and transmission in produce operations, we performed whole-genome sequencing (WGS) of LM (n = 169) and other Listeria spp. (n = 107) obtained from 13 produce packinghouses and three fresh-cut produce facilities. Overall, a low proportion of LM isolates (9/169) had inlA premature stop codons, and a large proportion (83/169) had either or both of the LIPI-3 or LIPI-4 operons, which have been associated with hypervirulence. The further analysis of the WGS data by operation showed a reisolation (at least 2 months apart) of highly related isolates (<10 hqSNP differences) in 7/16 operations. Two operations had highly related strains reisolated from samples that were collected at least 1 year apart. The identification of isolates collected during preproduction (i.e., following sanitation but before the start of production) that were highly related to isolates collected during production (i.e., after people or products have entered and begun moving through the operation) provided evidence that some strains were able to survive standard sanitation practices. The identification of closely related isolates (<20 hqSNPs differences) in different operations suggests that cross-contamination between facilities or introductions from common suppliers may also contribute to Listeria transmission. Overall, our data suggest that the majority of LM isolates collected from produce operations are fully virulent and that both persistence and reintroduction may lead to the repeat isolation of closely related Listeria in produce operations. IMPORTANCE Listeria monocytogenes is of particular concern to the produce industry due to its frequent presence in natural environments as well as its ability to survive in packinghouses and fresh-cut processing facilities over time. The use of whole-genome sequencing, which provides high discriminatory power for the characterization of Listeria isolates, along with detailed source data (isolation date and sample location) shows that the presence of Listeria in produce operations appears to be due to random and continued reintroduction as well as to the persistence of highly related strains in both packinghouses and fresh-cut facilities. These findings indicate the importance of using high-resolution characterization approaches for root cause analyses of Listeria contamination issues. In cases of repeat isolation of closely related Listeria in a given facility, both persistence and reintroduction need to be considered as possible root causes.

Recent Evolution and Genomic Profile of Salmonella enterica Serovar Heidelberg Isolates from Poultry Flocks in Brazil.

Salmonella enterica serovar Heidelberg is isolated from poultry-producing regions around the world. In Brazil, S. Heidelberg has been frequently detected in poultry flocks, slaughterhouses, and chicken meat. The goal of the present study was to assess the population structure, recent temporal evolution, and some important genetic characteristics of S. Heidelberg isolated from Brazilian poultry farms. Phylogenetic analysis of 68 S. Heidelberg genomes sequenced here and additional whole-genome data from NCBI demonstrated that all isolates from the Brazilian poultry production chain clustered into a monophyletic group, here called S. Heidelberg Brazilian poultry lineage (SH-BPL). Bayesian analysis defined the time of the most recent common ancestor (tMRCA) as 2004, and the overall population size (Ne) was constant until 2008, when an ∼10-fold Ne increase was observed until circa 2013. SH-BPL presented at least two plasmids with replicons ColpVC (n = 68; 100%), IncX1 (n = 66; 97%), IncA/C2 (n = 65; 95.5%), ColRNAI (n = 43; 63.2%), IncI1 (n = 32; 47%), ColMG828, Col156, IncHI2A, IncHI2, IncQ1, IncX4, IncY, and TrfA (each with n < 4; <4% each). Antibiotic resistance genes were found, with high frequencies of fosA7 (n = 68; 100%), mdf(A) (n = 68; 100%), tet(34) (n = 68; 100%), sul2 (n = 64; 94.1%), and blaCMY-2 (n = 56; 82.3%), along with an overall multidrug resistance (MDR) profile. Ten Salmonella pathogenicity islands (SPI1 to SPI5, SPI9, and SPI11 to SPI14) and 139 virulence genes were also detected. The SH-BPL profile was like those of other previous S. Heidelberg isolates from poultry around the world in the 1990s. In conclusion, the present study demonstrates the recent introduction (2004) and high level of dissemination of an MDR S. Heidelberg lineage in Brazilian poultry operations. IMPORTANCES. Heidelberg is the most frequent serovar in several broiler farms from the main Brazilian poultry-producing regions. Therefore, avian-source foods (mainly chicken carcasses) commercialized in the country and exported to other continents are contaminated with this foodborne pathogen, generating several national and international economic losses. In addition, isolates of this serovar are usually resistant to antibiotics and can cause human invasive and septicemic infection, representing a public health concern. This study demonstrates the use of whole-genome sequencing (WGS) to obtain epidemiological information for one S. Heidelberg lineage highly spread among Brazilian poultry farms. This information will help to define biosecurity measures to control this important Salmonella serovar in Brazilian and worldwide poultry operations.

Evolution of Listeria monocytogenes in a Food Processing Plant Involves Limited Single-Nucleotide Substitutions but Considerable Diversification by Gain and Loss of Prophages.

Whole-genome sequencing (WGS) is becoming the standard method for subtyping Listeria monocytogenes Interpretation of WGS data for isolates from foods and associated environments is, however, challenging due to a lack of detailed data on Listeria evolution in processing facilities. Here, we used previously collected WGS data for 40 L. monocytogenes isolates obtained from a cold-smoked salmon processing facility between 1998 and 2015 to probe the L. monocytogenes molecular evolution in this facility, combined with phenotypic assessment of selected isolates. Isolates represented three clusters (1, 2, and 3); cluster 3 isolates (n = 32) were obtained over 18 years. The average mutation rate for cluster 3 was estimated as 1.15 × 10-7 changes per nucleotide per year (∼0.35 changes per genome per year); the most recent common ancestors (MRCAs) of subclusters 3a and 3b were estimated to have occurred around 1958 and 1974, respectively, within the age of the facility, suggesting long-term persistence in this facility. Extensive prophage diversity was observed within subclusters 3a and 3b, which have one shared and six unique prophage profiles for each subcluster (with 16 prophage profiles found among all 40 isolates). The plasmid-borne sanitizer tolerance operon bcrABC was found in all cluster 2 and 3 isolates, while the transposon-borne sanitizer tolerance gene qacH was found in one cluster 1 isolate; presence of these genes was correlated with the ability to survive increased concentrations of sanitizers. Selected isolates showed significant variation in the ability to attach to surfaces, with persistent isolates attaching better than transient isolates at 21°C.IMPORTANCE Knowledge about the genetic evolution of L. monocytogenes in food processing facilities over multiple years is generally lacking. This information is critical to interpret WGS findings involving food or food-associated isolates. This study suggests that L. monocytogenes that persists in processing facilities may evolve with a low single-nucleotide mutation rate mostly driven by negative (i.e., purifying) selection but with rapid diversification of prophages. Hence, isolation of L. monocytogenes with few single-nucleotide polymorphism (SNP) differences in different locations (e.g., supplier plants and receiving plants) is possible, highlighting the importance of epidemiological and detailed isolate metadata for interpreting WGS data in traceback investigation. Our study also shows how advanced WGS data analyses can be used to support root cause analysis efforts and may, for example, pinpoint the time when a persistence event started (which then potentially could be linked to facility changes, introduction of new equipment, etc.).

A century of gray: A genomic locus found in 2 distinct Pseudomonas spp. is associated with historical and contemporary color defects in dairy products worldwide.

Some gram-negative bacteria, including Pseudomonas spp., can grow at refrigeration temperatures and cause flavor, odor, and texture defects in fluid milk. Historical and modern cases exist of gray and blue color defects in fluid milk due to Pseudomonas, and several recent reports have detailed fresh cheese spoilage associated with blue-pigment-forming Pseudomonas. Our goal was to investigate the genomes of pigmented Pseudomonas isolates responsible for historical and modern pigmented spoilage of dairy products in the United States to determine the genetic basis of pigment-forming phenotypes. We performed whole genome sequencing of 9 Pseudomonas isolates: 3 from recent incidents of gray-pigmented fluid milk (Pseudomonas fluorescens group), 1 from blue-pigmented cheese (P. fluorescens group), 2 from a historical blue milk spoilage incident (Pseudomonas putida group), and 3 with no evidence for blue or gray pigment formation (2 from P. fluorescens group and 1 from Pseudomonas chlororaphis group). All 6 isolates collected from products with a gray or blue pigment defect were confirmed to produce pigment using potato dextrose agar or pasteurized milk. A subset of 2 isolates was selected for inoculation into milk and onto the surface of a model cheese for subsequent color measurement. These isolates produced different colors on potato dextrose agar, but produced nearly identical color defects in milk and on model cheese. For the same subset of 2 isolates, the gray color defect in milk was produced only in containers with ample headspace and not in full containers, suggesting that oxygen is vital for pigment formation. This work also demonstrated that a Pseudomonas isolate from cheese can produce a pigment defect in milk, and vice versa. Comparative genomics identified an accessory locus encoding tryptophan biosynthesis genes that was present in all isolates that produced gray or blue pigment under laboratory conditions and was only previously reported in 2 P. fluorescens isolates responsible for blue mozzarella in Italy. Because this locus was found in genetically distant isolates belonging to different Pseudomonas species groups, it may have been acquired via horizontal gene transfer. These data suggest that several past and present gray- or blue-pigmented dairy spoilage events share a common genetic etiology that transcends species-level identification and merits further investigation to determine mechanistic details and modes of prevention.

Characteristics and distribution of Listeria spp., including Listeria species newly described since 2009.

The genus Listeria is currently comprised of 17 species, including 9 Listeria species newly described since 2009. Genomic and phenotypic data clearly define a distinct group of six species (Listeria sensu strictu) that share common phenotypic characteristics (e.g., ability to grow at low temperature, flagellar motility); this group includes the pathogen Listeria monocytogenes. The other 11 species (Listeria sensu lato) represent three distinct monophyletic groups, which may warrant recognition as separate genera. These three proposed genera do not contain pathogens, are non-motile (except for Listeria grayi), are able to reduce nitrate (except for Listeria floridensis), and are negative for the Voges-Proskauer test (except for L. grayi). Unlike all other Listeria species, species in the proposed new genus Mesolisteria are not able to grow below 7 °C. While most new Listeria species have only been identified in a few countries, the availability of molecular tools for rapid characterization of putative Listeria isolates will likely lead to future identification of isolates representing these new species from different sources. Identification of Listeria sensu lato isolates has not only allowed for a better understanding of the evolution of Listeria and virulence characteristics in Listeria but also has practical implications as detection of Listeria species is often used by the food industry as a marker to detect conditions that allow for presence, growth, and persistence of L. monocytogenes. This review will provide a comprehensive critical summary of our current understanding of the characteristics and distribution of the new Listeria species with a focus on Listeria sensu lato.

Prevalence and distribution of Listeria monocytogenes inlA alleles prone to phase variation and inlA alleles with premature stop codon mutations among human, food, animal, and environmental isolates.

In Listeria monocytogenes, 18 mutations leading to premature stop codons (PMSCs) in the virulence gene inlA have been identified to date. While most of these mutations represent nucleotide substitutions, a frameshift deletion in a 5' seven-adenine homopolymeric tract (HT) in inlA has also been reported. This HT may play a role in phase variation and was first identified among L. monocytogenes lineage II ribotype DUP-1039C isolates. In order to better understand the distribution of different inlA mutations in this ribotype, a newly developed multiplex real-time PCR assay was used to screen 368 DUP-1039C isolates from human, animal, and food-associated sources for three known 5' inlA HT alleles: (i) wild-type (WT) (A7), (ii) frameshift (FS) (A6), and (iii) guanine interruption (A2GA4) alleles. Additionally, 228 DUP-1039C isolates were screened for all inlA PMSCs; data on the presence of all inlA PMSCs for the other 140 isolates were obtained from previous studies. The statistical analysis based on 191 epidemiologically unrelated strains showed that strains with inlA PMSC mutations (n = 41) were overrepresented among food-associated isolates, while strains encoding full-length InlA (n = 150) were overrepresented among isolates from farm animals and their environments. Furthermore, the A6 allele was overrepresented and the A7 allele was underrepresented among food isolates, while the A6 allele was underrepresented among farm and animal isolates. Our results indicate that genetic variation in inlA contributes to niche adaptation within the lineage II subtype DUP-1039C.

Transcriptomic Analysis of the Adaptation of Listeria monocytogenes to Growth on Vacuum-Packed Cold Smoked Salmon.

The foodborne pathogen Listeria monocytogenes is able to survive and grow in ready-to-eat foods, in which it is likely to experience a number of environmental stresses due to refrigerated storage and the physicochemical properties of the food. Little is known about the specific molecular mechanisms underlying survival and growth of L. monocytogenes under different complex conditions on/in specific food matrices. Transcriptome sequencing (RNA-seq) was used to understand the transcriptional landscape of L. monocytogenes strain H7858 grown on cold smoked salmon (CSS; water phase salt, 4.65%; pH 6.1) relative to that in modified brain heart infusion broth (MBHIB; water phase salt, 4.65%; pH 6.1) at 7°C. Significant differential transcription of 149 genes was observed (false-discovery rate [FDR], <0.05; fold change, ≥2.5), and 88 and 61 genes were up- and downregulated, respectively, in H7858 grown on CSS relative to the genes in H7858 grown in MBHIB. In spite of these differences in transcriptomes under these two conditions, growth parameters for L. monocytogenes were not significantly different between CSS and MBHIB, indicating that the transcriptomic differences reflect how L. monocytogenes is able to facilitate growth under these different conditions. Differential expression analysis and Gene Ontology enrichment analysis indicated that genes encoding proteins involved in cobalamin biosynthesis as well as ethanolamine and 1,2-propanediol utilization have significantly higher transcript levels in H7858 grown on CSS than in that grown in MBHIB. Our data identify specific transcriptional profiles of L. monocytogenes growing on vacuum-packaged CSS, which may provide targets for the development of novel and improved strategies to control L. monocytogenes growth on this ready-to-eat food.

Genomic comparison of sporeforming bacilli isolated from milk.

BACKGROUND: Sporeformers in the order Bacillales are important contributors to spoilage of pasteurized milk. While only a few Bacillus and Viridibacillus strains can grow in milk at 6°C, the majority of Paenibacillus isolated from pasteurized fluid milk can grow under these conditions. To gain a better understanding of genomic features of these important spoilage organisms and to identify candidate genomic features that may facilitate cold growth in milk, we performed a comparative genomic analysis of selected dairy associated sporeformers representing isolates that can and cannot grow in milk at 6°C. RESULTS: The genomes for seven Paenibacillus spp., two Bacillus spp., and one Viridibacillus sp. isolates were sequenced. Across the genomes sequenced, we identified numerous genes encoding antimicrobial resistance mechanisms, bacteriocins, and pathways for synthesis of non-ribosomal peptide antibiotics. Phylogenetic analysis placed genomes representing Bacillus, Paenibacillus and Viridibacillus into three distinct well supported clades and further classified the Paenibacillus strains characterized here into three distinct clades, including (i) clade I, which contains one strain able to grow at 6°C in skim milk broth and one strain not able to grow under these conditions, (ii) clade II, which contains three strains able to grow at 6°C in skim milk broth, and (iii) clade III, which contains two strains unable to grow under these conditions. While all Paenibacillus genomes were found to include multiple copies of genes encoding ß-galactosidases, clade II strains showed significantly higher numbers of genes encoding these enzymes as compared to clade III strains. Genome comparison of strains able to grow at 6°C and strains unable to grow at this temperature identified numerous genes encoding features that might facilitate the growth of Paenibacillus in milk at 6°C, including peptidases with cold-adapted features (flexibility and disorder regions in the protein structure) and cold-adaptation related proteins (DEAD-box helicases, chaperone DnaJ). CONCLUSIONS: Through a comparative genomics approach we identified a number of genomic features that may relate to the ability of selected Paenibacillus strains to cause spoilage of refrigerated fluid milk. With additional experimental evidence, these data will facilitate identification of targets to detect and control Gram positive spore formers in fluid milk.

Phosphotransferase system-dependent extracellular growth of listeria monocytogenes is regulated by alternative sigma factors σL and σH.

Alternative sigma (σ) factors and phosphotransferase systems (PTSs) play pivotal roles in the environmental adaptation and virulence of Listeria monocytogenes. The growth of the L. monocytogenes parent strain 10403S and 15 isogenic alternative σ factor mutants was assessed in defined minimal medium (DM) with PTS-dependent or non-PTS-dependent carbon sources at 25°C or 37°C. Overall, our results suggested that the regulatory effect of alternative σ factors on the growth of L. monocytogenes is dependent on the temperature and the carbon source. One-way analysis of variance (one-way ANOVA) showed that the factor "strain" had a significant effect on the maximum growth rate (µmax), lag phase duration (λ), and maximum optical density (ODmax) in PTS-dependent carbon sources (P < 0.05) but not in a non-PTS-dependent carbon source. Also, the ODmax was not affected by strain for any of the three PTS-dependent carbon sources at 25°C but was affected by strain at 37°C. Monitoring by quantitative real-time PCR (qRT-PCR) showed that transcript levels for lmo0027, a glucose-glucoside PTS permease (PTS(Glc)-1)-encoding gene, were higher in the absence of σ(L), and lower in the absence of σ(H), than in the parent strain. Our data thus indicate that σ(L) negatively regulates lmo0027 and that the increased µmax observed for the ΔsigL strain in DM with glucose may be associated with increased expression of PTS(Glc)-1 encoded by lmo0027. Our findings suggest that σ(H) and σ(L) mediate the PTS-dependent growth of L. monocytogenes through complex transcriptional regulations and fine-tuning of the expression of specific pts genes, including lmo0027. Our findings also reveal a more important and complex role of alternative σ factors in the regulation of growth in different sugar sources than previously assumed.

Application of environmental monitoring programs and root cause analysis to identify and implement interventions to reduce or eliminate Listeria populations in apple packinghouses.

Controlling Listeria in produce packinghouses can be challenging due to the large number of potential contamination routes. For example, repeated isolation of the same Listeria subtype in a packinghouse could indicate persistence in the packinghouse or reintroduction of the same Listeria from an upstream source. To improve understanding of Listeria transmission patterns in packinghouses, we performed a longitudinal study in four apple packinghouses, including testing of 1,339 environmental sponges and whole genome sequencing (WGS)-based characterization of 280 isolates. Root cause analysis and subsequent intervention implementation was also performed and assessed for effectiveness. Listeria prevalence among environmental sponges collected from the four packinghouses was 20% (range of 5 to 31% for individual packinghouses). Sites that showed high Listeria prevalence included drains, forklift tires and forks, forklift stops, and waxing area equipment frames. A total of 240/280 WGS-characterized isolates were represented in 41 clusters, each containing two or more isolates that differed by ≤50 high-quality single nucleotide polymorphisms (hqSNPs); 21 clusters were isolated from one packinghouse over ≥2 samplings (suggesting persistence or possibly reintroduction), while 11 clusters included isolates from >2 packinghouses, suggesting common upstream sources. Some interventions successfully (i) reduced Listeria detection on forklift tires and forks (across packinghouses) and (ii) mitigated packinghouse-specific Listeria issues (e.g., in catch pans). However, interventions that lacked enhanced equipment disassembly when persistence was suspected typically appeared to be unsuccessful. Overall, while our data suggest a combination of intensive environmental sampling with subtyping and root cause analysis can help identify effective interventions, implementation of effective interventions continues to be a challenge in packinghouses.

Strain-specific Growth Parameters are Important to Accurately Model Bacterial Growth on Baby Spinach in Simulation Models.

Digital tools to predict produce shelf life have the potential to reduce food waste and improve consumer satisfaction. To address this need, we (i) performed an observational study on the microbial quality of baby spinach, (ii) completed growth experiments of bacteria that are representative of the baby spinach microbiota, and (iii) developed an initial simulation model of bacterial growth on baby spinach. Our observational data showed that the predominant genera found on baby spinach were Pseudomonas, Pantoea and Exiguobacterium. Rifampicin-resistant mutants (rifR mutants) of representative bacterial subtypes were subsequently generated to obtain strain-specific growth parameters on baby spinach. These experiments showed that: (i) it is difficult to select rifR mutants that do not have fitness costs affecting growth (9 of 15 rifR mutants showed substantial differences in growth, compared to their corresponding wild-type strain) and (ii) based on estimates from primary growth models, the mean (geometric) maximum population of rifR mutants on baby spinach (7.6 log10 CFU/g, at 6°C) appears lower than that of the spinach microbiota (9.6 log10 CFU/g, at 6°C), even if rifR mutants did not have substantial growth-related fitness costs. Thus, a simulation model, parameterized with the data obtained here as well as literature data on home refrigeration temperatures, underestimated bacterial growth on baby spinach. The root mean square error of the simulation's output, compared against data from the observational study, was 1.11 log10 CFU/g. Sensitivity analysis was used to identify key parameters (e.g., strain maximum population) that impact the simulation model's output, allowing for prioritization of future data collection to improve the simulation model. Overall, this study provides a roadmap for the development of models to predict bacterial growth on leafy vegetables with strain-specific parameters and suggests that additional data are required to improve these models.

Taxonomy, ecology, and relevance to food safety of the genus Listeria with a particular consideration of new Listeria species described between 2010 and 2022.

Since 2010, the genus Listeria has had the addition of 22 new species that more than tripled the number of species identified until 2010. Sixteen of these 22 new species are distantly related to the type species, Listeria monocytogenes, and several of these present phenotypes that distinguish them from classical Listeria species (L. monocytogenes, Listeria innocua, Listeria ivanovii, Listeria seeligeri, Listeria welshimeri, and Listeria grayi). These 22 newly described species also show that Listeria is more genetically diverse than previously estimated. While future studies and surveys are needed to clarify the distribution of these species, at least some of these species may not be widely spread, while other species may be frequently found spread to human-related settings (e.g., farms and processing facilities), and others may be adapted to specific environmental habitats. Here, we review the taxonomic, phylogenetic, and ecological characteristics of these new Listeria species identified since 2010 and re-iterate the suggestion of re-classification of some species into three new genera: Murraya, Mesolisteria, and Paenilisteria. We also provide a review of current detection issues and the relevance to food safety related to the identification of these new species. For example, several new non-pathogenic species could be misidentified as the pathogen L. monocytogenes, based on methods that do not target L. monocytogenes-specific virulence genes/factors, leading to unnecessary product recalls. Moreover, eight species in the proposed new genus Mesolisteria are not good indicators of environmental conditions that could allow L. monocytogenes to grow since Mesolisteria species are unable to grow at low temperatures.

Role of Whole Genome Sequencing in Assessing Resident and Transient Listeria monocytogenes in a Produce Packinghouse.

Whole genome sequencing (WGS) is a powerful tool that may be used to assist in identifying Listeria contamination sources and movement within environments, and to assess persistence. This study investigated sites in a produce packinghouse where Listeria had been historically isolated; and aimed to characterize dispersal patterns and identify cases of transient and resident Listeria. Environmental swab samples (n = 402) were collected from 67 sites at two time-points on three separate visits. Each sample was tested for Listeria, and Listeria isolates were characterized by partial sigB sequencing to determine species and allelic type (AT). Representative isolates from the three most common L. monocytogenes ATs (n = 79) were further characterized by WGS. Of the 144 Listeria species positive samples (35.8%), L. monocytogenes was the most prevalent species. L. monocytogenes was often coisolated with another species of Listeria. WGS identified cases of sporadic and continued reintroduction of L. monocytogenes from the cold storages into the packinghouse and demonstrated cases of L. monocytogenes persistence over 2 years in cold storages, drains, and on a forklift. Nine distinct clusters were found in this study. Two clusters showed evidence of persistence. Isolates in these two clusters (N = 11, with one historical isolate) were obtained predominantly and over multiple samplings from cold storages, with sporadic movement to sites in the packing area, suggesting residence in cold storages with opportunistic dispersal within the packinghouse. The other seven clusters demonstrated evidence of transient Listeria, as isolation was sporadic over time and space during the packing season. Our data provide important insights into likely L. monocytogenes harborage points and transfer in a packinghouse, which is key to root cause analysis. While results support Listeria spp. as a suitable indicator organism for environmental monitoring surveys, findings were unable to establish a specific species as an index organism for L. monocytogenes. Findings also suggest long-term persistence with substantial SNP diversification, which may assist in identifying potential contamination sources and implementing control measures.

Intensive Environmental Sampling and Whole Genome Sequence-based Characterization of Listeria in Small- and Medium-sized Dairy Facilities Reveal Opportunities for Simplified and Size-appropriate Environmental Monitoring Strategies.

Small- and medium-sized dairy processing facilities (SMDFs) may face unique challenges with respect to controlling Listeria in their processing environments, e.g., due to limited resources. The aim of this study was to implement and evaluate environmental monitoring programs (EMPs) for Listeria control in eight SMDFs in a ∼1-year longitudinal study; this included a comparison of pre-operation (i.e., after cleaning and sanitation and prior to production) and mid-operation (i.e., at least 4 h into production) sampling strategies. Among 2,072 environmental sponge samples collected across all facilities, 272 (13%) were positive for Listeria. Listeria prevalence among pre- and mid-operation samples (15% and 17%, respectively), was not significantly different. Whole genome sequencing (WGS) performed on select isolates to characterize Listeria persistence patterns revealed repeated isolation of closely related Listeria isolates (i.e., ≤20 high-quality single nucleotide polymorphism [hqSNP] differences) in 5/8 facilities over >6 months, suggesting Listeria persistence and/or reintroduction was relatively common among the SMDFs evaluated here. WGS furthermore showed that for 41 sites where samples collected pre- and mid-operation were positive for Listeria, Listeria isolates obtained were highly related (i.e., ≤10 hqSNP differences), suggesting that pre-operation sampling alone may be sufficient and more effective for detecting sites of Listeria persistence. Importantly, our data also showed that only 1/8 of facilities showed a significant decrease in Listeria prevalence over 1 year, indicating continued challenges with Listeria control in at least some SMDFs. We conclude that options for simplified Listeria EMPs (e.g., with a focus on pre-operation sampling, which allows for more rapid identification of likely persistence sites) may be valuable for improved Listeria control in SMDFs.

Genomic characterization provides new insight into Salmonella phage diversity.

BACKGROUND: Salmonella is a widely distributed foodborne pathogen that causes tens of millions of salmonellosis cases globally every year. While the genomic diversity of Salmonella is increasingly well studied, our knowledge of Salmonella phage genomic diversity is still rather limited, despite the contributions of both lysogenic and lytic phages to Salmonella virulence, diversity and ecology (e.g., through horizontal gene transfer and Salmonella lysis). To gain a better understanding of phage diversity in a specific ecological niche, we sequenced 22 Salmonella phages isolated from a number of dairy farms from New York State (United States) and analyzed them using a comparative genomics approach. RESULTS: Classification of the 22 phages according to the presence/absence of orthologous genes allowed for classification into 8 well supported clusters. In addition to two phage clusters that represent novel virulent Salmonella phages, we also identified four phage clusters that each contained previously characterized phages from multiple continents. Our analyses also identified two clusters of phages that carry putative virulence (e.g., adhesins) and antimicrobial resistance (tellurite and bicyclomycin) genes as well as virulent and temperate transducing phages. Insights into phage evolution from our analyses include (i) identification of DNA metabolism genes that may facilitate nucleotide synthesis in phages with a G+C % distinct from Salmonella, and (ii) evidence of Salmonella phage tailspike and fiber diversity due to both single nucleotide polymorphisms and major re-arrangements, which may affect the host specificity of Salmonella phages. CONCLUSIONS: Genomics-based characterization of 22 Salmonella phages isolated from dairy farms allowed for identification of a number of novel Salmonella phages. While the comparative genomics analyses of these phages provide a number of new insights in the evolution and diversity of Salmonella phages, they only represent a first glimpse into the diversity of Salmonella phages that is likely to be discovered when phages from different environments are characterized.

Protein level identification of the Listeria monocytogenes sigma H, sigma L, and sigma C regulons.

BACKGROUND: Transcriptional regulation by alternative sigma (σ) factors represents an important mechanism that allows bacteria to rapidly regulate transcript and protein levels in response to changing environmental conditions. While the role of the alternative σ factor σB has been comparatively well characterized in L. monocytogenes, our understanding of the roles of the three other L. monocytogenes alternative σ factors is still limited. In this study, we employed a quantitative proteomics approach using Isobaric Tags for Relative and Absolute Quantitation (iTRAQ) to characterize the L. monocytogenes σL, σH, and σC protein regulons. Proteomic comparisons used a quadruple alternative σ factor mutant strain (ΔBCHL) and strains expressing a single alternative σ factor (i.e., σL, σH, and σC; strains ΔBCH, ΔBCL, and ΔBHL) to eliminate potential redundancies between σ factors. RESULTS: Among the three alternative σ factors studied here, σH provides positive regulation for the largest number of proteins, consistent with previous transcriptomic studies, while σL appears to contribute to negative regulation of a number of proteins. σC was found to regulate a small number of proteins in L. monocytogenes grown to stationary phase at 37°C. Proteins identified as being regulated by multiple alternative σ factors include MptA, which is a component of a PTS system with a potential role in regulation of PrfA activity. CONCLUSIONS: This study provides initial insights into global regulation of protein production by the L. monocytogenes alternative σ factors σL, σH, and σC. While, among these σ factors, σH appears to positively regulate the largest number of proteins, we also identified PTS systems that appear to be co-regulated by multiple alternative σ factors. Future studies should not only explore potential roles of alternative σ factors in activating a "cascade" of PTS systems that potentially regulate PrfA, but also may want to explore the σL and σC regulons under different environmental conditions to identify conditions where these σ factors may regulate larger numbers of proteins or genes.

The Number and Type of Chaperone-Usher Fimbriae Reflect Phylogenetic Clade Rather than Host Range in Salmonella.

Salmonella is one of the most successful foodborne pathogens worldwide, owing in part to its ability to colonize or infect a wide range of hosts. Salmonella serovars are known to encode a variety of different fimbriae (hairlike organelles that facilitate binding to surfaces); however, the distribution, number, and sequence diversity of fimbriae encoded across different lineages of Salmonella were unknown. We queried whole-genome sequence (WGS) data for 242 Salmonella enterica subsp. enterica (subspecies enterica) isolates from the top 217 serovars associated with isolation from humans and agricultural animals; this effort identified 2,894 chaperone-usher (CU)-type fimbrial usher sequences, representing the most conserved component of CU fimbriae. On average, isolates encoded 12 different CU fimbrial ushers (6 to 18 per genome), although the distribution varied significantly (P = 1.328E-08) by phylogenetic clade, with isolates in section Typhi having significantly fewer fimbrial ushers than isolates in clade A2 (medians = 10 and 12 ushers, respectively). Characterization of fimbriae in additional non-enterica subspecies genomes suggested that 8 fimbrial ushers were classified as being unique to subspecies enterica isolates, suggesting that the majority of fimbriae were most likely acquired prior to the divergence of subspecies enterica. Characterization of mobile elements suggested that plasmids represent an important vehicle facilitating the acquisition of a wide range of fimbrial ushers, particularly for the acquisition of fimbriae from other Gram-negative genera. Overall, our results suggest that differences in the number and type of fimbriae encoded most likely reflect differences in phylogenetic clade rather than differences in host range. IMPORTANCE Fimbriae of the CU assembly pathway represent important organelles that mediate Salmonella's interactions with host tissues and abiotic surfaces. Our analyses provide a comprehensive overview of the diversity of CU fimbriae in Salmonella spp., highlighting that the majority of CU fimbriae are distributed broadly across multiple subspecies and suggesting that acquisition most likely occurred prior to the divergence of subspecies enterica. Our data also suggest that plasmids represent the primary vehicles facilitating the horizontal transfer of diverse CU fimbriae in Salmonella. Finally, the observed high sequence similarity between some ushers suggests that different names may have been assigned to closely related fimbrial ushers that likely should be represented by a single designation. This highlights the need to establish standard criteria for fimbria classification and nomenclature, which will also facilitate future studies seeking to associate virulence factors with adaptation to or differences in the likelihood of causing disease in a given host.

Development of a Modeling Tool To Assess and Reduce Regulatory and Recall Risks for Cold-Smoked Salmon Due to Listeria monocytogenes Contamination.

ABSTRACT: Although public health risk assessments for Listeria monocytogenes (Lm) have been published for various foods, firm-level decision making on interventions targeting Lm involves considerations of both public health and enterprise risks. Smoked seafood is a ready-to-eat product with a high incidence of Lm contamination and has been associated with several recalls. We used cold-smoked salmon as a model product to develop a decision support tool (the regulatory and recall risk [3R] model) to estimate (i) baseline regulatory and recall (RR) risks (i.e., overall risks of a lot sampled and found positive for Lm, e.g., by food regulatory agencies) due to Lm contamination and (ii) the RR risk reduction that can be achieved through interventions with underlying mechanisms such as reducing the prevalence and/or level of Lm and retarding or preventing Lm growth. Given that a set number of samples (e.g., 10) are tested for a given lot, the RR risk equals the likelihood of detecting Lm in at least one sample. Under the baseline scenario, which assumes a 4% Lm prevalence and no interventions, the median predicted RR risk for a given production lot was 0.333 (95% credible interval: 0.288, 0.384) when 10 25-g samples were tested. Nisin treatments, which reduce both the prevalence and initial level of Lm, reduced RR risks in a concentration-dependent manner to 0.109 (0.074, 0.146) with 5 ppm, 0.049 (0.024, 0.083) with 10 ppm, and 0.017 (0.007, 0.033) with 20 ppm. In general, more effective reduction in RR risks can be achieved by reducing Lm prevalence than by retarding Lm growth; the RR risk was reduced to 0.182 (0.153, 0.213) by a 50% prevalence reduction but to only 0.313 (0.268, 0.367) by bacteriostatic growth inhibitors. Sensitivity analysis indicated that prevalence and initial level of Lm and storage temperature have the greatest impact on predicting RR risks, suggesting that reliable data for these parameters will improve model performance.

Development of a Genomics-Based Approach To Identify Putative Hypervirulent Nontyphoidal Salmonella Isolates: Salmonella enterica Serovar Saintpaul as a Model.

While differences in human virulence have been reported across nontyphoidal Salmonella (NTS) serovars and associated subtypes, a rational and scalable approach to identify Salmonella subtypes with differential ability to cause human diseases is not available. Here, we used NTS serovar Saintpaul (S. Saintpaul) as a model to determine if metadata and associated whole-genome sequence (WGS) data in the NCBI Pathogen Detection (PD) database can be used to identify (i) subtypes with differential likelihoods of causing human diseases and (ii) genes and single nucleotide polymorphisms (SNPs) potentially responsible for such differences. S. Saintpaul SNP clusters (n = 211) were assigned different epidemiology types (epi-types) based on statistically significant over- or underrepresentation of human clinical isolates, including human associated (HA; n = 29), non-human associated (NHA; n = 23), and other (n = 159). Comparative genomic analyses identified 384 and 619 genes overrepresented among isolates in 5 HA and 4 NHA SNP clusters most significantly associated with the respective isolation source. These genes included 5 HA-associated virulence genes previously reported to be present on Gifsy-1/Gifsy-2 prophages. Additionally, premature stop codons in 3 and 7 genes were overrepresented among the selected HA and NHA SNP clusters, respectively. Tissue culture experiments with strains representing 4 HA and 3 NHA SNP clusters did not reveal evidence for enhanced invasion or intracellular survival for HA strains. However, the presence of sodCI (encoding a superoxide dismutase), found in 4 HA and 1 NHA SNP clusters, was positively correlated with intracellular survival in macrophage-like cells. Post hoc analyses also suggested a possible difference in intracellular survival among S. Saintpaul lineages. IMPORTANCE Not all Salmonella isolates are equally likely to cause human disease, and Salmonella control strategies may unintentionally focus on serovars and subtypes with high prevalence in source populations but are rarely associated with human clinical illness. We describe a framework leveraging WGS data in the NCBI PD database to identify Salmonella subtypes over- and underrepresented among human clinical cases. While we identified genomic signatures associated with HA/NHA SNP clusters, tissue culture experiments failed to identify consistent phenotypic characteristics indicative of enhanced human virulence of HA strains. Our findings illustrate the challenges of defining hypo- and hypervirulent S. Saintpaul and potential limitations of phenotypic assays when evaluating human virulence, for which in vivo experiments are essential. Identification of sodCI, an HA-associated virulence gene associated with enhanced intracellular survival, however, illustrates the potential of the framework and is consistent with prior work identifying specific genomic features responsible for enhanced or reduced virulence of nontyphoidal Salmonella.