Utilizing whole genome sequencing to characterize Listeria spp. persistence and transmission patterns in a farmstead dairy processing facility and its associated farm environment.

Farmstead dairy processing facilities may be particularly susceptible to Listeria spp. contamination due to the close physical proximity of their processing environments (PE) to associated dairy farm environments (FE). In this case study, we supported the implementation of interventions focused on improving (i) cleaning and sanitation efficacy, (ii) hygienic zoning, and (iii) sanitary equipment/facility design and maintenance in a farmstead dairy processing facility, and evaluated their impact on Listeria spp. detection in the farmstead's PE over 1 year. Detection of Listeria spp. in the farmstead's PE was numerically reduced from 50% to 7.5% after 1 year of intervention implementation, suggesting that these interventions were effective at improving Listeria spp. control. In addition, environmental samples were also collected from the farmstead's FE to evaluate the risk of the FE as a potential source of Listeria spp. in the PE. Overall, detection of Listeria spp. was higher in samples collected from the FE (75%, 27/36) compared with samples collected from the PE (24%, 29/120). Whole genome sequencing (WGS) performed on select isolates collected from the PE and FE supported the identification of 6 clusters (range of 3 to 15 isolates per cluster) that showed ≤ 50 high quality single nucleotide polymorphism (hqSNP) differences. Of these 6 clusters, 3 (i.e., clusters 2, 4, and 5) contained isolates that were collected from both the PE and FE, suggesting that transmission between these 2 environments was likely. Moreover, all cluster 2 isolates represented a clonal complex (CC) of L. monocytogenes commonly associated with dairy farm environmental reservoirs (i.e., CC666), which may support that the farmstead's FE represented an upstream source of the cluster 2 isolates that were found in the PE. Overall, our data underscore that, while the FE can represent a potential upstream source of Listeria spp. contamination in a farmstead dairy processing facility, implementation of targeted interventions can help effectively minimize Listeria spp. contamination in the PE.

A longitudinal study on the bacterial quality of baby spinach cultivated in Arizona and California.

In the U.S., baby spinach is mostly produced in Arizona (AZ) and California (CA). Characterizing the impact of growing region on the bacterial quality of baby spinach can inform quality management practices in industry. Between December 2021 and December 2022, baby spinach was sampled after harvest and packaging for microbiological testing, including shelf-life testing of packaged samples that were stored at 4°C. Samples were tested to (i) determine bacterial concentration, and (ii) obtain and identify bacterial isolates. Packaged samples from the Salinas, CA, area (n = 13), compared to those from the Yuma, AZ, area (n = 9), had a significantly higher bacterial concentration, on average, by 0.78 log10 CFU/g (P < 0.01, based on aerobic, mesophilic plate count data) or 0.67 log10 CFU/g (P < 0.01, based on psychrotolerant plate count data); the bacterial concentrations of harvest samples from the Yuma and Salinas areas were not significantly different. Our data also support that an increase in preharvest temperature is significantly associated with an increase in the bacterial concentration on harvested and packaged spinach. A Fisher's exact test and linear discriminant analysis (effect size), respectively, demonstrated that (i) the genera of 2,186 bacterial isolates were associated (P < 0.01) with growing region and (ii) Pseudomonas spp. and Exiguobacterium spp. were enriched in spinach from the Yuma and Salinas areas, respectively. Our findings provide preliminary evidence that growing region and preharvest temperature may impact the bacterial quality of spinach and thus could inform more targeted strategies to manage produce quality. IMPORTANCE: In the U.S., most spinach is produced in Arizona (AZ) and California (CA) seasonally; typically, spinach is cultivated in the Yuma, AZ, area during the winter and in the Salinas, CA, area during the summer. As the bacterial quality of baby spinach can influence consumer acceptance of the product, it is important to assess whether the bacterial quality of baby spinach can vary between spinach-growing regions. The findings of this study provide insights that could be used to support region-specific quality management strategies for baby spinach. Our results also highlight the value of further evaluating the impact of growing region and preharvest temperature on the bacterial quality of different produce commodities.

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 were also performed and assessed for effectiveness. Listeria prevalence among environmental sponges collected from the four packinghouses was 20% (range of 5-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.

Assessment of the exposure to cytotoxic Bacillus cereus group genotypes through HTST milk consumption.

This study addresses the limited tools available for assessing food safety risks from cytotoxic Bacillus cereus group strains in contaminated food. We quantified the growth, in skim milk broth, of 17 cytotoxic B. cereus strains across 6 phylogenetic groups with various virulence gene profiles. The strains did not grow in HTST milk at 4 or 6°C. At 10°C, 15 strains exhibited growth; at 8°C, one strain grew; and all strains grew at temperatures ≥ 14°C. Using growth data from 16 strains, we developed linear secondary growth models and an exposure assessment model. This model, simulating a 5-stage HTST milk supply chain and up to 35 d of consumer storage with an initial contamination of 100 cfu/mL, estimated that 2.81 ± 0.66% and 4.13 ± 2.53% of milk containers would surpass 105 cfu/mL of B. cereus by d 21 and 35, respectively. A sensitivity analysis identified the initial physiological state of cells (Q0) as the most influential variable affecting predictions for specific isolates. What-if scenarios indicated that increases in mean and variability of consumer storage temperatures significantly affected the predicted B. cereus concentrations in milk. This model serves as an initial tool for risk-based food safety decision making regarding low-level B. cereus contamination.

Risk Assessment Predicts Most of the Salmonellosis Risk in Raw Chicken Parts is Concentrated in Those Few Products with High Levels of High-Virulence Serotypes of Salmonella.

Salmonella prevalence declined in U.S. raw poultry products since adopting prevalence-based Salmonella performance standards, but human illnesses did not reduce proportionally. We used Quantitative Microbial Risk Assessment (QMRA) to evaluate public health risks of raw chicken parts contaminated with different levels of all Salmonella and specific high- and low-virulence serotypes. Lognormal Salmonella level distributions were fitted to 2012 USDA-FSIS Baseline parts survey and 2023 USDA-FSIS HACCP verification sampling data. Three different Dose-Response (DR) approaches included (i) a single DR for all serotypes, (ii) DR that reduces Salmonella Kentucky ST152 virulence, and (iii) multiple serotype-specific DR models. All scenarios found risk concentrated in the few products with high Salmonella levels. Using a single DR model with Baseline data (µ = -3.19, σ = 1.29 Log CFU/g), 68% and 37% of illnesses were attributed to the 0.7% and 0.06% of products with >1 and >10 CFU/g Salmonella, respectively. Using distributions from 2023 HACCP data (µ = -5.53, σ = 2.45), 99.8% and 99.0% of illnesses were attributed to the 1.3% and 0.4% of products with >1 and >10 CFU/g Salmonella, respectively. Scenarios with serotype-specific DR models showed more concentrated risk at higher levels. Baseline data showed 92% and 67% and HACCP data showed >99.99% and 99.96% of illnesses attributed to products with >1 and >10 CFU/g Salmonella, respectively. Regarding serotypes using Baseline or HACCP input data, 0.002% and 0.1% of illnesses were attributed to the 0.2% and 0.4% of products with >1 CFU/g of Kentucky ST152, respectively, while 69% and 83% of illnesses were attributed to the 0.3% and 0.6% of products with >1 CFU/g of Enteritidis, Infantis, or Typhimurium, respectively. Therefore, public health risk in chicken parts is concentrated in finished products with high levels and specifically high levels of high-virulence serotypes. Low-virulence serotypes like Kentucky contribute few human cases.

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 nonsporeformers as well as sporeformers.

Laboratory pasteurization count (LPC) enumerates thermoduric bacteria and is one parameter used to assess raw milk quality. No regulatory limit has presently been set for LPC, but 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. Although it is occasionally used as a proxy for levels of bacterial spores in raw milk, limited knowledge is available on the types of organisms that are enumerated by LPC in contemporary raw milk supplies. Although 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 nonsporeformers, and gram-negatives, respectively; 0.5% 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 nonsporeformers were predominantly 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 MALDI-TOF mass spectrometry (MALDI-TOF MS), a rapid microbial identification tool that is more readily available to dairy producers than 16S rDNA PCR and sequencing. Although our data indicated agreement between 16S rDNA sequencing 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 MS, 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.