Microbial food spoilage: impact, causative agents and control strategies.

Microbial food spoilage is a major contributor to food waste and, hence, to the negative environmental sustainability impacts of food production and processing. Globally, it is estimated that 15-20% of food is wasted, with waste, by definition, occurring after primary production and harvesting (for example, in households and food service establishments). Although the causative agents of food spoilage are diverse, many microorganisms are major contributors across different types of foods. For example, the genus Pseudomonas causes spoilage in various raw and ready-to-eat foods. Aerobic sporeformers (for example, members of the genera Bacillus, Paenibacillus and Alicyclobacillus) cause spoilage across various foods and beverages, whereas anaerobic sporeformers (for example, Clostridiales) cause spoilage in a range of products that present low-oxygen environments. Fungi are also important spoilage microorganisms, including in products that are not susceptible to bacterial spoilage due to their low water activity or low pH. Strategies that can reduce spoilage include improved control of spoilage microorganisms in raw material and environmental sources as well as application of microbicidal or microbiostatic strategies (for example, to products and packaging). Emerging tools (for example, systems models and improved genomic tools) represent an opportunity for rational design of systems, processes and products that minimize microbial food spoilage.

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.

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.

Troubleshooting high laboratory pasteurization counts in organic raw milk requires characterization of dominant thermoduric bacteria, which includes non-sporeformers as well as sporeformers.

Laboratory Pasteurization Count (LPC) enumerates thermoduric bacteria and is one parameter used to assess raw milk quality. While there is currently no regulatory limit for LPC, LPC data are used by some dairy processors and cooperatives to designate raw milk quality premiums paid to farmers and may also be used for troubleshooting bacterial contamination issues. Despite occasionally being used as a proxy for levels of bacterial spores in raw milk, there is limited knowledge of the types of organisms that are enumerated by LPC in contemporary raw milk supplies. While historical studies have reported that thermoduric bacteria quantified by LPC may predominantly represent Gram-positive cocci, updated knowledge on microbial populations enumerated by LPC in contemporary organic raw milk supplies is needed. To address this gap, organic raw milk samples from across the United States (n = 94) were assessed using LPC, and bacterial isolates were characterized. LPC ranged from below detection (<0.70 log cfu/mL) to 4.07 log cfu/mL, with a geometric mean of 1.48 log cfu/mL. Among 380 isolates characterized by 16S rDNA sequencing, 52.6%, 44.5%, and 2.4% were identified as Gram-positive sporeformers, Gram-positive non-sporeformers, and Gram-negatives, respectively, and 0.5% that could not be categorized into those groups because they could only be assigned a higher level of taxonomy. Isolates identified as Gram-positive sporeformers were predominantly Bacillus (168/200) and Gram-positive non-sporeformers were predominately Brachybacterium (56/169) and Kocuria (47/169). To elucidate if the LPC level can be an indicator of the type of thermoduric (e.g., sporeforming bacteria) present in raw milk, we evaluated the proportion of sporeformers in raw milk samples with LPC of ≤100 cfu/mL, 100 to 200 cfu/mL, and ≥200 cfu/mL (51%, 67%, and 35%), showing a trend for sporeformers to represent a smaller proportion of the total thermoduric population when LPC increases, although overall linear regression showed no significant association between the proportion of sporeformers and the LPC concentration. Hence, LPC level alone provides no insight into the makeup of the thermoduric population in raw milk and further characterization is needed to elucidate the bacterial drivers of elevated LPC in raw milk. We therefore further characterized the isolates from this study using matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI-TOF MS), a rapid microbial identification tool that is more readily available to dairy producers than 16S rDNA PCR and sequencing. While our data indicated agreement between 16S and MALDI-TOF MS for 66.6% of isolates at the genus level, 24.2% and 9.2% could not be reliably identified or were mischaracterized using MALDI-TOF, respectively. This suggests that further optimization of this method is needed to allow for accurate characterization of thermoduric organisms commonly found in raw milk. Ultimately, our study provides a contemporary perspective on thermoduric bacteria selected by the LPC method and establishes that the LPC alone is not sufficient for identifying the bacterial drivers of LPC levels. Further development of rapid characterization methods that are accessible to producers, cooperatives, and processors will support milk quality troubleshooting efforts and ultimately improve outcomes for dairy industry community members.

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.

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.

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.

Population dynamics of Listeria spp., Salmonella spp., and Escherichia coli on fresh produce: A scoping review.

Collation of the current scope of literature related to population dynamics (i.e., growth, die-off, survival) of foodborne pathogens on fresh produce can aid in informing future research directions and help stakeholders identify relevant research literature. A scoping review was conducted to gather and synthesize literature that investigates population dynamics of pathogenic and non-pathogenic Listeria spp., Salmonella spp., and Escherichia coli on whole unprocessed fresh produce (defined as produce not having undergone chopping, cutting, homogenization, irradiation, or pasteurization). Literature sources were identified using an exhaustive search of research and industry reports published prior to September 23, 2021, followed by screening for relevance based on strict, a priori eligibility criteria. A total of 277 studies that met all eligibility criteria were subjected to an in-depth qualitative review of various factors (e.g., produce commodities, study settings, inoculation methodologies) that affect population dynamics. Included studies represent investigations of population dynamics on produce before (i.e., pre-harvest; n = 143) and after (i.e., post-harvest; n = 144) harvest. Several knowledge gaps were identified, including the limited representation of (i) pre-harvest studies that investigated population dynamics of Listeria spp. on produce (n = 13, 9% of pre-harvest studies), (ii) pre-harvest studies that were carried out on non-sprouts produce types grown using hydroponic cultivation practices (n = 7, 5% of pre-harvest studies), and (iii) post-harvest studies that reported the relative humidity conditions under which experiments were carried out (n = 56, 39% of post-harvest studies). These and other knowledge gaps summarized in this scoping review represent areas of research that can be investigated in future studies.

Selection of mutant Listeria phages under food-relevant conditions can enhance application potential.

Bacteriophages are viruses that infect and kill bacteria. Currently, phage products are available for the control of the pathogen Listeria monocytogenes in food products in the United States. In this study, we explore whether experimental evolution can be used to generate phages with improved abilities to function under specific food-relevant conditions. Ultra-pasteurized oat and whole milk were chosen as test matrices as they represent different food groups, yet have similar physical traits and macronutrient composition. We showed that (i) wild-type phage LP-125 infection kinetics are different in the two matrices and (ii) LP-125 has a significantly higher burst size in oat milk. From this, we attempted to evolve LP-125 to have improved infection kinetics in whole milk. Ancestral LP-125 was passaged through 10 rounds of amplification in milk conditions. Plaque-purified DNA samples from milk-selected phages were isolated and sequenced, and mutations present in the isolated phages were identified. We found two nonsynonymous substitutions in LP125_108 and LP125_112 genes, which encode putative baseplate-associated glycerophosphoryl diester phosphodiesterase and baseplate protein, respectively. Protein structural modeling showed that the substituted amino acids in the mutant phages are predicted to localize to surface-exposed helices on the corresponding structures, which might affect the surface charge of proteins and their interaction with the bacterial cell. The phage containing the LP125_112 mutation adsorbed significantly faster than the ancestral phage in both oat and whole milk. Follow-up experiments suggest that fat content may be a key factor for the expression of the phenotype of this mutation. IMPORTANCE Bacteriophages are one of the tools available to control the foodborne pathogen, Listeria monocytogenes. Phage products must work under a broad range of food conditions to be an effective control for L. monocytogenes. Here, we show that the experimental evolution of phages can be used to generate new phages with phenotypes useful under specific conditions. We used this approach to select for a mutant phage that more efficiently binds to L. monocytogenes that is grown in whole milk and oat milk. We show that the fat content of these milks is necessary for the expression of this phenotype. Our findings show that experimental evolution can be used to select for improved phages with better performance under specific conditions. This approach has the potential to support the development of condition-specific phage-based biocontrols in the food industry.

Comparative genomics unveils extensive genomic variation between populations of Listeria species in natural and food-associated environments.

Comprehending bacterial genomic variation linked to distinct environments can yield novel insights into mechanisms underlying differential adaptation and transmission of microbes across environments. Gaining such insights is particularly crucial for pathogens as it benefits public health surveillance. However, the understanding of bacterial genomic variation is limited by a scarcity of investigations in genomic variation coupled with different ecological contexts. To address this limitation, we focused on Listeria, an important bacterial genus for food safety that includes the human pathogen L. monocytogenes, and analyzed a large-scale genomic dataset collected by us from natural and food-associated environments across the United States. Through comparative genomics analyses on 449 isolates from the soil and 390 isolates from agricultural water and produce processing facilities representing L. monocytogenes, L. seeligeri, L. innocua, and L. welshimeri, we find that the genomic profiles strongly differ by environments within each species. This is supported by the environment-associated subclades and differential presence of plasmids, stress islands, and accessory genes involved in cell envelope biogenesis and carbohydrate transport and metabolism. Core genomes of Listeria species are also strongly associated with environments and can accurately predict isolation sources at the lineage level in L. monocytogenes using machine learning. We find that the large environment-associated genomic variation in Listeria appears to be jointly driven by soil property, climate, land use, and accompanying bacterial species, chiefly representing Actinobacteria and Proteobacteria. Collectively, our data suggest that populations of Listeria species have genetically adapted to different environments, which may limit their transmission from natural to food-associated environments.

More than mcr: canonical plasmid- and transposon-encoded mobilized colistin resistance genes represent a subset of phosphoethanolamine transferases.

Mobilized colistin resistance genes (mcr) may confer resistance to the last-resort antimicrobial colistin and can often be transmitted horizontally. mcr encode phosphoethanolamine transferases (PET), which are closely related to chromosomally encoded, intrinsic lipid modification PET (i-PET; e.g., EptA, EptB, CptA). To gain insight into the evolution of mcr within the context of i-PET, we identified 69,814 MCR-like proteins present across 256 bacterial genera (obtained by querying known MCR family representatives against the National Center for Biotechnology Information [NCBI] non-redundant protein database via protein BLAST). We subsequently identified 125 putative novel mcr-like genes, which were located on the same contig as (i) ≥1 plasmid replicon and (ii) ≥1 additional antimicrobial resistance gene (obtained by querying the PlasmidFinder database and NCBI's National Database of Antibiotic Resistant Organisms, respectively, via nucleotide BLAST). At 80% amino acid identity, these putative novel MCR-like proteins formed 13 clusters, five of which represented putative novel MCR families. Sequence similarity and a maximum likelihood phylogeny of mcr, putative novel mcr-like, and ipet genes indicated that sequence similarity was insufficient to discriminate mcr from ipet genes. A mixed-effect model of evolution (MEME) indicated that site- and branch-specific positive selection played a role in the evolution of alleles within the mcr-2 and mcr-9 families. MEME suggested that positive selection played a role in the diversification of several residues in structurally important regions, including (i) a bridging region that connects the membrane-bound and catalytic periplasmic domains, and (ii) a periplasmic loop juxtaposing the substrate entry tunnel. Moreover, eptA and mcr were localized within different genomic contexts. Canonical eptA genes were typically chromosomally encoded in an operon with a two-component regulatory system or adjacent to a TetR-type regulator. Conversely, mcr were represented by single-gene operons or adjacent to pap2 and dgkA, which encode a PAP2 family lipid A phosphatase and diacylglycerol kinase, respectively. Our data suggest that eptA can give rise to "colistin resistance genes" through various mechanisms, including mobilization, selection, and diversification of genomic context and regulatory pathways. These mechanisms likely altered gene expression levels and enzyme activity, allowing bona fide eptA to evolve to function in colistin resistance.

Escherichia coli B-Strains Are Intrinsically Resistant to Colistin and Not Suitable for Characterization and Identification of mcr Genes.

Antimicrobial resistance is an increasing threat to human and animal health. Due to the rise of multi-, extensive, and pandrug resistance, last resort antibiotics, such as colistin, are extremely important in human medicine. While the distribution of colistin resistance genes can be tracked through sequencing methods, phenotypic characterization of putative antimicrobial resistance (AMR) genes is still important to confirm the phenotype conferred by different genes. While heterologous expression of AMR genes (e.g., in Escherichia coli) is a common approach, so far, no standard methods for heterologous expression and characterization of mcr genes exist. E. coli B-strains, designed for optimum protein expression, are frequently utilized. Here, we report that four E. coli B-strains are intrinsically resistant to colistin (MIC 8-16 µg/mL). The three tested B-strains that encode T7 RNA polymerase show growth defects when transformed with empty or mcr-expressing pET17b plasmids and grown in the presence of IPTG; K-12 or B-strains without T7 RNA polymerase do not show these growth defects. E. coli SHuffle T7 express carrying empty pET17b also skips wells in colistin MIC assays in the presence of IPTG. These phenotypes could explain why B-strains were erroneously reported as colistin susceptible. Analysis of existing genome data identified one nonsynonymous change in each pmrA and pmrB in all four E. coli B-strains; the E121K change in PmrB has previously been linked to intrinsic colistin resistance. We conclude that E. coli B-strains are not appropriate heterologous expression hosts for identification and characterization of mcr genes. IMPORTANCE Given the rise in multidrug, extensive drug, and pandrug resistance in bacteria and the increasing use of colistin to treat human infections, occurrence of mcr genes threatens human health, and characterization of these resistance genes becomes more important. We show that three commonly used heterologous expression strains are intrinsically resistant to colistin. This is important because these strains have previously been used to characterize and identify new mobile colistin resistance (mcr) genes. We also show that expression plasmids (i.e., pET17b) without inserts cause cell viability defects when carried by B-strains with T7 RNA polymerase and grown in the presence of IPTG. Our findings are important as they will facilitate improved selection of heterologous strains and plasmid combinations for characterizing AMR genes, which will be particularly important with a shift to Culture-independent diagnostic tests where bacterial isolates become increasingly less available for characterization.

Culture-independent bacterial cell extraction from fluid milk and oat-based beverage for basic qualitative microscopy.

Butterfat and protein complicate attempts to extract bacterial cells from milk by centrifugation for use in basic microscopy. Some types of bacteria preferentially separate into the butterfat layer upon centrifugation and are lost when this layer is discarded, and the action of bacterial protease enzymes can cause milk proteins to precipitate and partition into the centrifugal pellet. Butterfat and precipitated protein remaining in the centrifugal pellet along with the desired bacterial cells can confound the results of differential staining and microscopy. Oat- and other plant-based beverages, which are often manufactured by dairy processors on shared equipment, present similar hurdles to bacterial extraction and microscopic visualization because of the presence of oils, starch granules, and dietary fiber particles in these products. Herein we describe methods for centrifugal separation of bacterial cells for microscopy from unflavored milk, chocolate milk, and oat-based beverage. Cell suspensions prepared through these methods were used for phase-contrast microscopy, Gram staining, and viability staining. These techniques can be used to provide rapid, culture-independent diagnostic information when bacterial cells are expected to be present in high concentrations, as in the event of sporadic product spoilage or mass product spoilage incidents.

Restriction endonuclease cleavage of phage DNA enables resuscitation from Cas13-induced bacterial dormancy.

Type VI CRISPR systems protect against phage infection using the RNA-guided nuclease Cas13 to recognize viral messenger RNA. Upon target recognition, Cas13 cleaves phage and host transcripts non-specifically, leading to cell dormancy that is incompatible with phage propagation. However, whether and how infected cells recover from dormancy is unclear. Here we show that type VI CRISPR and DNA-cleaving restriction-modification (RM) systems frequently co-occur and synergize to clear phage infections and resuscitate cells. In the natural type VI CRISPR host Listeria seeligeri, we show that RM cleaves the phage genome, thus removing the source of phage transcripts and enabling cells to recover from Cas13-induced cellular dormancy. We find that phage infections are neutralized more effectively when Cas13 and RM systems operate together. Our work reveals that type VI CRISPR immunity is cell-autonomous and non-abortive when paired with RM, and hints at other synergistic roles for the diverse host-directed immune systems in bacteria.

Low-cost, on-farm intervention to reduce spores in bulk tank raw milk benefits producers, processors, and consumers.

Bacterial spores, which are found in raw milk, can survive harsh processing conditions encountered in dairy manufacturing, including pasteurization and drying. Low-spore raw milk is desirable for dairy industry stakeholders, especially those who want to extend the shelf life of their product, expand their distribution channels, or reduce product spoilage. A recent previous study showed that an on-farm intervention that included washing towels with chlorine bleach and drying them completely, as well as training milking parlor employees to focus on teat end cleaning, significantly reduced spore levels in bulk tank raw milk. As a follow up to that previous study, here we calculate the costs associated with that previously described intervention as ranging from $9.49 to $13.35 per cow per year, depending on farm size. A Monte Carlo model was used to predict the shelf life of high temperature, short time fluid milk processed from raw milk before and after this low-cost intervention was applied, based on experimental data collected in a previous study. The model predicted that 18.24% of half-gallon containers of fluid milk processed from raw milk receiving no spore intervention would exceed the pasteurized milk ordinance limit of 20,000 cfu/mL by 17 d after pasteurization, while only 16.99% of containers processed from raw milk receiving the spore intervention would reach this level 17 d after pasteurization (a reduction of 1.25 percentage points and a 6.85% reduction). Finally, a survey of consumer milk use was conducted to determine how many consumers regularly consume fluid milk near or past the date printed on the package (i.e., code date), which revealed that over 50% of fluid milk consumers surveyed continue to consume fluid milk after this date, indicating that a considerable proportion of consumers are exposed to fluid milk that is likely to have high levels spore-forming bacterial growth and possibly associated quality defects (e.g., flavor or odor defects). This further highlights the importance of reducing spore levels in raw milk to extend pasteurized fluid milk shelf life and thereby reducing the risk of adverse consumer experiences. Processors who are interested in extending fluid milk shelf life by controlling the levels of spores in the raw milk supply should consider incentivizing low-spore raw milk through premium payments to producers.

Salmonella enterica serovar Cerro displays a phylogenetic structure and genomic features consistent with virulence attenuation and adaptation to cattle.

Salmonella enterica subsp. enterica (S.) serovar Cerro is rarely isolated from human clinical cases of salmonellosis but represents the most common serovar isolated from cattle without clinical signs of illness in the United States. In this study, using a large, diverse set of 316 isolates, we utilized genomic methods to further elucidate the evolutionary history of S. Cerro and to identify genomic features associated with its apparent virulence attenuation in humans. Phylogenetic analyses showed that within this polyphyletic serovar, 98.4% of isolates (311/316) represent a monophyletic clade within section Typhi and the remaining 1.6% of isolates (5/316) form a monophyletic clade within subspecies enterica Clade A1. Of the section Typhi S. Cerro isolates, 93.2% of isolates (290/311) clustered into a large clonal clade comprised of predominantly sequence type (ST) 367 cattle and environmental isolates, while the remaining 6.8% of isolates (21/311), primarily from human clinical sources, clustered outside of this clonal clade. A tip-dated phylogeny of S. Cerro ST367 identified two major clades (I and II), one of which overwhelmingly consisted of cattle isolates that share a most recent common ancestor that existed circa 1975. Gene presence/absence and rarefaction curve analyses suggested that the pangenome of section Typhi S. Cerro is open, potentially reflecting the gain/loss of prophage; human isolates contained the most open pangenome, while cattle isolates had the least open pangenome. Hypothetically disrupted coding sequences (HDCs) displayed clade-specific losses of intact speC and sopA virulence genes within the large clonal S. Cerro clade, while loss of intact vgrG, araH, and vapC occurred in all section Typhi S. Cerro isolates. Further phenotypic analysis suggested that the presence of a premature stop codon in speC does not abolish ornithine decarboxylase activity in S. Cerro, likely due to the activity of the second ornithine decarboxylase encoded by speF, which remained intact in all isolates. Overall, our study identifies specific genomic features associated with S. Cerro's infrequent isolation from humans and its apparent adaptation to cattle, which has broader implications for informing our understanding of the evolutionary events facilitating host adaptation in Salmonella.

Structural Equation Models Suggest That On-Farm Noncrop Vegetation Removal Is Not Associated with Improved Food Safety Outcomes but Is Linked to Impaired Water Quality.

While growers have reported pressures to minimize wildlife intrusion into produce fields through noncrop vegetation (NCV) removal, NCV provides key ecosystem services. To model food safety and environmental tradeoffs associated with NCV removal, published and publicly available food safety and water quality data from the Northeastern United States were obtained. Because data on NCV removal are not widely available, forest-wetland cover was used as a proxy, consistent with previous studies. Structural equation models (SEMs) were used to quantify the effect of forest-wetland cover on (i) food safety outcomes (e.g., detecting pathogens in soil) and (ii) water quality (e.g., nutrient levels). Based on the SEMs, NCV was not associated with or had a protective effect on food safety outcomes (more NCV was associated with a reduced likelihood of pathogen detection). The probabilities of detecting Listeria spp. in soil (effect estimate [EE] = -0.17; P = 0.005) and enterohemorrhagic Escherichia coli in stream samples (EE = -0.27; P < 0.001) were negatively associated with the amount of NCV surrounding the sampling site. Larger amounts of NCV were also associated with lower nutrient, salinity, and sediment levels, and higher dissolved oxygen levels. Total phosphorous levels were negatively associated with the amount of NCV in the upstream watershed (EE = -0.27; P < 0.001). Similar negative associations (P < 0.05) were observed for other physicochemical parameters, such as nitrate (EE = -0.38). Our findings suggest that NCV should not be considered an inherent produce safety risk or result in farm audit demerits. This study also provides a framework for evaluating environmental tradeoffs associated with using specific preharvest food safety strategies. IMPORTANCE Currently, on-farm food safety decisions are typically made independently of conservation considerations, often with detrimental impacts on agroecosystems. Comanaging agricultural environments to simultaneously meet conservation and food safety aims is complicated because farms are closely linked to surrounding environments, and management decisions can have unexpected environmental, economic, and food safety consequences. Thus, there is a need for research on the conservation and food safety tradeoffs associated with implementing specific preharvest food safety practices. Understanding these tradeoffs is critical for developing adaptive comanagement strategies and ensuring the short- and long-term safety, sustainability, and profitability of agricultural systems. This study quantifies tradeoffs and synergies between food safety and environmental aims, and outlines a framework for modeling tradeoffs and synergies between management aims that can be used to support future comanagement research.

A Validated Preharvest Sampling Simulation Shows that Sampling Plans with a Larger Number of Randomly Located Samples Perform Better than Typical Sampling Plans in Detecting Representative Point-Source and Widespread Hazards in Leafy Green Fields.

Commercial leafy greens customers often require a negative preharvest pathogen test, typically by compositing 60 produce sample grabs of 150 to 375 g total mass from lots of various acreages. This study developed a preharvest sampling Monte Carlo simulation, validated it against literature and experimental trials, and used it to suggest improvements to sampling plans. The simulation was validated by outputting six simulated ranges of positive samples that contained the experimental number of positive samples (range, 2 to 139 positives) recovered from six field trials with point source, systematic, and sporadic contamination. We then evaluated the relative performance between simple random, stratified random, or systematic sampling in a 1-acre field to detect point sources of contamination present at 0.3% to 1.7% prevalence. Randomized sampling was optimal because of lower variability in probability of acceptance. Optimized sampling was applied to detect an industry-relevant point source [3 log(CFU/g) over 0.3% of the field] and widespread contamination [-1 to -4 log(CFU/g) over the whole field] by taking 60 to 1,200 sample grabs of 3 g. More samples increased the power of detecting point source contamination, as the median probability of acceptance decreased from 85% with 60 samples to 5% with 1,200 samples. Sampling plans with larger total composite sample mass increased power to detect low-level, widespread contamination, as the median probability of acceptance with -3 log(CFU/g) contamination decreased from 85% with a 150-g total mass to 30% with a 1,200-g total mass. Therefore, preharvest sampling power increases by taking more, smaller samples with randomization, up to the constraints of total grabs and mass feasible or required for a food safety objective. IMPORTANCE This study addresses a need for improved preharvest sampling plans for pathogen detection in leafy green fields by developing and validating a preharvest sampling simulation model, avoiding the expensive task of physical sampling in many fields. Validated preharvest sampling simulations were used to develop guidance for preharvest sampling protocols. Sampling simulations predicted that sampling plans with randomization are less variable in their power to detect low-prevalence point source contamination in a 1-acre field. Collecting larger total sample masses improved the power of sampling plans in detecting widespread contamination in 1-acre fields. Hence, the power of typical sampling plans that collect 150 to 375 g per composite sample can be improved by taking more, randomized smaller samples for larger total sample mass. The improved sampling plans are subject to feasibility constraints or to meet a particular food safety objective.

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.