Next Article in Journal
The Probiotic Bacillus subtilis MB40 Improves Immunity in a Porcine Model of Listeriosis
Previous Article in Journal
Gut Microbiota and Respiratory Infections: Insights from Mendelian Randomization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Phenotype Nutritional Profiles and Phosphate Metabolism Genes in Four Serovars of Salmonella enterica from Water Sources

Produce Safety and Microbiology Research Unit, Agricultural Research Service, United States Department of Agriculture, Albany, CA 94710, USA
*
Author to whom correspondence should be addressed.
Microorganisms 2023, 11(8), 2109; https://doi.org/10.3390/microorganisms11082109
Submission received: 7 June 2023 / Revised: 3 August 2023 / Accepted: 14 August 2023 / Published: 18 August 2023
(This article belongs to the Section Food Microbiology)

Abstract

:
The surveillance of foods for Salmonella is hindered by bias in common enrichment media where serovars implicated in human illness are outgrown by less virulent serovars. We examined four Salmonella serovars, two common in human illness (Enteritidis and Typhimurium) and two that often dominate enrichments (Give and Kentucky), for factors that might influence culture bias. The four serovars had similar growth kinetics in Tryptic Soy Broth and Buffered Peptone Water. Phenotype microarray analysis with 950 chemical substrates to assess nutrient utilization and stress resistance revealed phenotype differences between serovars. Strains of S. Enteritidis had better utilization of plant-derived sugars such as xylose, mannitol, rhamnose, and fructose, while S. Typhimurium strains were able to metabolize tagatose. Strains of S. Kentucky used more compounds as phosphorus sources and grew better with inorganic phosphate as the sole phosphorus source. The sequences of nine genes involved in phosphate metabolism were compared, and there were differences between serovars in the catalytic ATP-binding domain of the histidine kinase phoR. Analysis of the predicted PhoR amino acid sequences from additional Salmonella genomes indicated a conservation of sequences each within the Typhimurium, Give, and Enteritidis serovars. However, three different PhoR versions were observed in S. Kentucky.

1. Introduction

The most common agent of foodborne illness is non-typhoidal Salmonella enterica [1,2]. Recency-weighted statistical modeling attributes most salmonellosis as resulting from contaminated seeded vegetables, eggs, chicken, produce, and beef [3]. There are over 2650 different serovars of Salmonella, but most human illness is caused by approximately 100 serovars, with the U.S. Centers for Disease Control and Prevention recognizing the top 20 serovars implicated in human illness in annual reports. The serovars Typhimurium and Enteritidis are among the top 5 causes of salmonellosis in the U.S. and the top 2 in the European Union and China [4,5,6]. Contamination of foods and water with multiple serovars is not uncommon. Using CRISPR-SeroSeq, Cox et al. reported an average of five different Salmonella serovars present on broiler carcasses [7]. A recent outbreak of salmonellosis in France associated with dried pork sausages was traced back to two different serovars, Bovismorbificans and monophasic Typhimurium [8]. Similarly, several surveys of surface water report the isolation of multiple serovars of Salmonella from the same sample [9,10].
Laboratory methods for Salmonella surveillance and traceback investigations of Salmonella outbreaks involve suspension of the food sample in a primary enrichment or pre-enrichment medium and incubation for 18–24 h at 35 °C–37 °C, followed by secondary enrichments and plating on various selective and differential media. Standard methods used in various surveillance laboratories prescribe several different complex, liquid media for pre-enrichments and primary enrichments [11,12,13]. Standard pre-enrichment media are dilute, complex media that are based on the protein hydrolysates of casein or peptone that allow for the recovery of stressed and injured cells that might be present in contaminated foods [11,12,14]. Buffered Peptone Water (BPW) is recommended in the methods published by the USDA’s Food Safety and Inspection Service (USDA-FSIS) and the International Organization for Standardization (ISO) to pre-enrich Salmonella from chicken and other foods [11,12,13]. Tryptic Soy Broth (TSB) is a common peptide-based, laboratory growth medium that is routinely used for Salmonella growth and is also used in non-selective enrichment of Salmonella from eggs and spices according to the Food and Drug Administration (FDA) method [11].
Culture bias in laboratory media influences the serovars of Salmonella that arise from samples containing multiple subtypes [7,15,16]. The Salmonella serovar Kentucky is prevalent in poultry in the U.S. and is frequently isolated in Salmonella enrichments, where it outgrows strains of serovars such as Typhimurium and Enteritidis [7,17,18]. Similarly, S. Give is frequently isolated from surface water and from ruminants [9,19,20,21,22]. The bias can be problematic when conducting food surveillance or outbreak tracebacks because the serovars Kentucky and Give are rarely associated with human illness in the U.S. [4]. In the U.S., S. Kentucky and S. Give are more fit in commonly used Salmonella enrichment protocols and can outgrow and mask the presence of strains of S. Typhimurium and S. Enteritidis [7,15,22,23]. In Europe, the Middle East, and several African countries, an antibiotic-resistant S. Kentucky Sequence Type 198 has arisen to cause multiple cases of human illness. While ST198 has caused human illness in the U.S., the number of cases is markedly lower [24,25].
In BPW and TSB, peptides serve as the sources of carbon, nitrogen, sulfur, phosphorus, and energy for the microbiota. Differences between strains in the utilization of media components may affect the final recovery of strains in the mixture. These differences may be intensified when nutrients affect many roles in cellular metabolism. For example, proteomic analysis of the S. Typhimurium strain 14028s revealed over 300 proteins that are regulated by phosphorus-limiting conditions [26].
Phosphorus metabolism via the uptake of inorganic phosphate (Pi) has effects on magnesium and calcium homeostasis [27,28]. In Salmonella, phosphate homeostasis is controlled by protein synthesis by responding to reductions in ATP, destabilized ribosomes, or low cytoplasmic magnesium [29]. Phosphate levels are sensed and controlled through the actions of PhoB/PhoR, a two-component regulatory system in which PhoR is an integral membrane-bound sensor histidine kinase, which autophosphorylates when Pi levels are low, and transfers that phosphoryl group to PhoB. PhoB is a transcriptional regulator responsible for the moderation of gene expression of Pho genes. One of the operons affected by PhoB is the pstSCAB operon, which encodes a high-affinity ABC-type transporter responsible for the transport of Pi into the cell [30,31,32]. When intracellular Pi concentrations are high, the transporter is turned off through the interaction of the PhoU protein with the PstB component of the transporter. This PhoU–PstB interaction is mediated through the action of the sensor kinase PhoR [29,31]. Also, during periods of low Pi, the Pho regulon plays a role in maintaining the stability of the stationary phase and stress response sigma factor RpoS through the actions of the small protein IraP. IraP assures the continued action of RpoS by binding to and inhibiting the RssB adaptor protein so that RssB is unable to bind to RpoS and target it for degradation [30,32,33].
In addition to phosphate metabolism, low nutrient conditions induce ppGpp, which slows protein synthesis and may induce a stress response [34]. Differences between strains in the utilization or efficiency of utilization of basic nutrient compounds result in fitness variances based on the competition for nutrients and lead to differences in recovery of subtypes [35,36]. It is possible that minor differences in the sequence of proteins used for nutrient utilization may affect efficiency of action, which in turn might influence recovery in mixed cultures.
A systematic study of the nutritional capacity of Salmonella serovars is needed to assess potential nutrient utilization differences between serovars. This information would be helpful in the design of culture media to bolster poor performers such as Typhimurium and Enteritidis and/or lower the competitiveness of good performers such as Kentucky and Give during routine surveillance of food and water. Similarly, the genetic differences resulting in such nutritional differences also need investigation.

2. Materials and Methods

2.1. Strains and Growth Conditions

The strains used in this study are listed in Table 1. Most of the strains were isolated in our laboratory as part of a large surface-water survey [10]. Strains were chosen based on the serovar, and strains within serovars were of different Pulsed Field Gel Electrophoresis (PFGE) pulsotypes isolated at different times and from different locations. Serovar identity was previously determined using sero-agglutination [10]. Because S. Kentucky strains were not well represented in that survey, additional S. Kentucky strains from the strain collection in the Produce Safety and Microbiology research unit were added. In addition to the strains used for growth and nutrition experiments in Table 1, genome information for additional analysis was downloaded from Enterobase [37,38].
For general use, strains were routinely grown in Trypticase Soy Broth (TSB, Difco, Becton Dickinson, Franklin Lakes, NJ, USA) and Trypticase Soy Agar via incubation overnight at 37 °C. Liquid cultures were incubated in a shaking incubator (Infors, Annapolis Junction, MD, USA) set to 150 rpm. The Phosphate Buffered Saline (PBS) contained 150 mM of NaCl and 10 mM of sodium phosphate, pH 7.2.

2.2. Growth Kinetics

Growth kinetics of individual strains were measured in TSB and BPW (Difco). Single colonies were inoculated into 4 mL of TSB or BPW and incubated with shaking for 16 h–18 h at 37 °C. A portion (1.0–1.5 mL) of the resulting culture was centrifuged at 13,000× g in a microfuge for 3–5 min, and the cell pellets were resuspended in PBS to an A600 of 0.05–0.1. This cell suspension (35 µL) was used to inoculate 35 mL of fresh, room-temperature medium, resulting in a starting cell concentration of approximately 105 CFU/mL. The TSB and BPW cultures were incubated with shaking at 37 °C. Turbidity was monitored by measuring A600 at regular intervals. Cell density was monitored by making serial dilutions of cultures in PBS and plating onto TSA. Plates were incubated for 18 h–24 h at 37 °C, and the colonies were counted to calculate the CFU/mL. These values were log-transformed, and the data were plotted as log CFU/mL.
Growth kinetics were modeled with a modified Gompertz equation using Prism 9.5 (GraphPad, San Diego, CA, USA) [39]. Growth curves were generated, and non-linear regression was performed using the Gompertz option in Prism. Best-fit values were derived from the curves in Prism according to Equation (1):
Y = YM × (Y0/YM)^exp(−K × X)
where Y0 is the starting population, YM is the maximum population level, K determines the lag time, and 1/K is the X value of the curve inflection point. Therefore, the 1/K value represents the duration of lag phase. Generation time was calculated during lag phase using the following equation:
g = t/n
where g is generation time, t represents the time interval, and n is the number of generations within the time interval. The number of generations within a time interval calculated from the growth curves using the following equation:
Nt = No × 2n
where Nt is the cfu/mL at time t, No is the cfu/mL at time 0, and n is the number of generations during the time interval.

2.3. Phenotype Microarrays

Metabolic characteristics of selected strains were determined using phenotype microarrays (Omnilog, Biolog, Hayward, CA, USA). Bacteria were grown on Biolog Universal Growth Agar with 5% sheep blood, and bacteria were harvested and inoculated into Phenotype Microarray (PM) plates according to the manufacturer’s instructions [40]. PM plates 1–10 were used to assess the phenotype characteristics of strains RM15596, RM17366, RM18387, RM18411, RM18433, RM18465, RM19068, and RM20875 using pre-configured 96-well plates (Biolog). PM plates 1–10 measured growth resulting from 950 different compounds to test for metabolism of carbon sources (PM-1 and PM-2), nitrogen and peptide nitrogen sources (PM-3, PM-6, PM-7, and PM-8), phosphorus and sulfur sources (PM-4), nutritional supplements (PM-5), osmotic and ionic effects (PM-9), and pH effects (PM-10) [40]. Phenotypes for each strain were determined in three replicate experiments and analyzed using Data Analysis v. 1.7 software (Biolog). Data were expressed as Omnilog Units (OUs), which are a spectrophotometric-based measure of metabolism based on the intensity of a redox dye in each well of the PM plates. The range of OUs is 0–500. Each plate contained a negative control well to serve as a background value. Results for strain replicate experiments were compared to ensure reproducibility, and then strains and serovars were compared with each other using the Data Analysis v. 1.7 software. The background values were subtracted from the experimental OU values, and averages and standard deviations were calculated in the Data Analysis v. 1.7 software. The OU values of the wells for each compound were compared between serovars. For comparison purposes, the average OUs for each well were differentiated into three levels: <20 OUs, indicating little to no utilization; 20–100 OUs, indicating moderate levels of utilization; and >100 OUs, indicating high levels of utilization. Levels such as these have been used for phenotype microarray analysis by other investigators [41,42].

2.4. Genome Sequencing and CRISPR Sequence Comparisons

The genomes of RM15596, RM17366, RM18387, RM18411, RM18433, RM18465, RM19068, and RM20875 were sequenced using the Illumina Mi-Seq platform with 500-cycle v2 paired-end sequencing kits (Illumina, San Diego, CA, USA). DNA libraries were prepared and quantified using the KAPA Low-Throughput Library Preparation Kit and KAPA Library Quantification Kit (Roche, Kapa Biosystems, Wilmington, MA, USA) with modifications as described in [43] and validated on Agilent 2100 Bioanalyzer chips using the Agilent 12,000 kit (Agilent, Santa Clara, CA, USA). Sequences were trimmed and draft sequences were uploaded to Enterobase, where sequence quality was assessed, de novo assemblies were constructed, and the seven-gene Multi-Locus Sequence Type (MLST) was determined [37,38]. Serovars were confirmed with SeqSero2 [44]. Sequence coverage ranged from 67X to 300X. Genomes were annotated in RAST v2.0 (https://rast.nmpdr.org/rast.cgi, accessed on 2 May 2023) [45].
Genome sequence data were deposited with the NCBI under the BioProject ID PRJNA875079. The raw sequencing data were sent to the Sequence Read Archive and are identified by the BioSample accession numbers SAMN30597303–SAMN30597310. Strains were identified by their ID numbers: RM15596, RM17366, RM18387, RM18411, RM18433, RM18465, RM19068, and RM20875.

2.5. Gene and Protein Sequence Alignments

The DNA and predicted protein sequences of several genes involved in inorganic phosphate metabolism were compared. The identity and function of the genes of interest are shown in Table 2. The sequences of phoB, phoR, phoU, phoP, phoQ, iraP, pstA, pstB, pstC, pstS, rpoS, and rssB from the sequenced strains (indicated with an asterisk (*) in Table 1) were imported into Geneious Prime. DNA and predicted protein sequences were aligned using Geneious and/or MegAlign Pro v. 15.0 (DNAStar, Inc., Madison, WI, USA).
Thirty-eight additional Salmonella Enteritidis, Give, Kentucky, and Typhimurium genomes were used for added analysis of phoR sequences. Some of the additional genomes were from Salmonella strains in our collection that were previously sequenced [10]. Additional Salmonella genomes were downloaded from Enterobase [37,38]. The additional genomes used are indicated in Supplementary Table S1. Genome sequences were annotated in RAST, and the phoR gene sequences were determined. Alignments of predicted PhoR amino acid sequences were conducted in Geneious Prime.

3. Results

3.1. The Serovars Had Similar Growth Kinetics in TSB and BPW

Growth in TSB and BPW was measured in the strains indicated in Table 1. Growth curves over 24 h are shown in Figure 1A. Figure 1B shows the first 10 h of the total 24 h growth to highlight the lag and log phases. The data for the individual strains were grouped according to serovar. The growth curves had similar shapes to each other, and the maximum CFU/mL numbers reached were similar for the serovars in pure culture in each medium. Non-linear regression analysis was performed on the growth data, and the model revealed that the serovars grew at similar rates in both media (Table 3). The lag phase in the Enteritidis strains was longer than in the other serovars (3.0 h vs. 2.3 h). Generation times were similar in TSB and BPW at approximately 0.5 h, with serovars Enteritidis and Typhimurium growing approximately 10% faster in pure cultures. The CFU/mL counts were approximately 1 log lower in BPW, consistent with BPW being a more dilute medium than TSB.

3.2. The Serovars Had Phenotype Differences in Nutrient Utilization

A phenotype microarray analysis with pre-configured 96-well plates was used to assess the effect of and/or utilization of 950 compounds for carbon, nitrogen, phosphorus, and sulfur sources, as well as nutritional supplements and osmotic and pH tolerance. The phenotypes demonstrated in the wells of the plates was reflected by redox dye intensity as explained in the Section 2 and were measured in Omnilog Units (OUs), which range from 0–500. For analysis purposes, OUs were differentiated into three categories: little to no growth (<20 OUs), intermediate growth (20–100 OUs), and high growth (>100 OUs). The data were examined according to the individual strains and average of each serovar. Under these criteria, there were 157 wells among the 960 analyzed that showed differential OU categories between the serovars (Supplementary Table S2). These could be further differentiated by of the number of OUs between the different levels. There were 125 wells where there was >100 OUs between the highest utilizing serovar for a compound and one of the other serovars. A total of 84 of these wells measured nitrogen utilization of dipeptides, 70 of which had the highest levels of utilization (OUs ranging from 100–130) in serovar Enteritidis strains. The remaining 41 phenotype wells that showed differential utilization levels according to serovar are shown in Figure 2. The data for the individual strains are given in Supplementary Table S2. Of the 41 traits, 21 were carbon sources, 6 were nitrogen sources, 9 were phosphorus sources, 2 were sulfur sources, and 3 measured osmotic sensitivity. S. Enteritidis strains had the highest OU numbers among the carbon sources, showing utilization of several compounds that were not used by at least one of the other serovars. These included plant-derived sugars such as xylose, mannitol, rhamnose, fructose, and melibiose. S. Typhimurium strains showed little utilization of 8 of the carbon sources among this set of 21. However, S. Typhimurium were the only strains positive for utilization of D-tagatose, which is found in dairy products and can be derived from lactose. There was variability between the serovars regarding the non-peptide sources of nitrogen in the compounds shown in Figure 2. The average values of the S. Kentucky strains were the most diverse among the four serovars in the use of phosphorus sources, showing more or faster growth with inorganic phosphate, thiophosphate, phosphoenol pyruvate, and several nucleotide phosphates as the sole phosphorus sources. The S. Kentucky strain RM16068 had higher utilization of the phosphorus compounds than did the S. Kentucky strain RM17366, but both strains performed better than strains of the other serovars for these compounds (Supplementary Table S2). Due to the better performance of the S. Kentucky strains in the phosphate wells and the importance of the Pi utilization to several aspects of metabolism, we examined the phosphorus assimilation genes in the strains of different serovars.

3.3. PhoR Gene and Amino Acid Sequences Differed between Serovars

The predicted protein sequences of the phoB, phoR, phoU, pstA, pstB, pstC, iraP, rpoS, and rrsB genes among strains of the different serovars were compared, and the lengths of the nucleic acid and predicted protein sequences were conserved. For the sequenced strains in this study (those with an asterisk (*) in Table 1), the predicted amino acid sequences of PhoU, PstB, PstC, IraP, RpoS, and RrsB were identical (Supplementary Figure S1). Five of these six genes had SNP differences between the eight strains and four serovars (Supplementary Figure S2), but the changes resulted in the same predicted protein sequence. The DNA sequences for the iraP gene were identical in all the strains (Supplementary Figure S2). Among these six genes, the largest number of SNP differences were in phoU, which ranged from 0–24 SNPs between the strains (Supplementary Figure S2).
The remaining three genes involved in Pi metabolism (phoB, phoR, and pstA) had differences in the nucleic acid and in the predicted protein sequences between the eight strains and four serovars (Figure 3; Supplementary Figures S1 and S2). The phoBR and pstA genes were the same size, with phoB at 690 bp, phoR at 1296 bp, and pstA at 891 bp. The ranges of SNP differences were 0–11 for phoB, 0–13 for phoR, and 0–24 for pstA. The predicted protein sequences of the transcription regulator PhoB were identical in seven of the eight strains; the sole difference was S. Kentucky RM17366, which had a Met73Leu substitution. The predicted amino acid sequence of pstA varied only for the S. Give strains RM18465 and RM20875, which each had an Ile120Val substitution. This variant would result in a conserved substitution in this inner membrane protein responsible for Pi translocation across the membrane [46]. The gene with the most predicted amino acid substitutions was the histidine kinase phoR (Figure 3). There were four different PhoR protein sequences among the eight strains. The two S. Typhimurium strains and one of the S. Kentucky strains (RM17366) shared identical PhoR amino acid sequences. Both S. Enteritidis strains had an Ala347Val substitution of this sequence. The two S. Give strains and S. Kentucky strain RM19068 had a Val346Ala substitution, and RM19068 had an additional Thr340Lys substitution. The valine–alanine substitutions were conservative amino acid substitutions. Threonine–lysine resulted in R groups with a charge difference (neutral—basic, respectively). All of the amino acid changes in the predicted PhoR sequence were in the catalytic ATP-binding (CA) domain of the protein responsible for autophosphorylation of this sensor histidine kinase [47].
To determine if these PhoR amino acid differences between serovars were present in additional Salmonella genomes, additional PhoR sequences were compared from the genomes of serovars Enteritidis, Give, Kentucky, and Typhimurium. Including the 8 strains highlighted in Figure 3, the PhoR sequences from 10 Enteritidis strains, 10 Give strains, 15 Kentucky strains, and 11 Typhimurium strains were compared (Figure 4). All the S. Typhimurium PhoR sequences were identical to each other. All the S. Enteritidis had the same Ala347Val substation indicated above (Figure 3), and 9 of the 10 were identical. S. Enteritidis strain 1729765 contained a Leu22Phe, which is a conserved substitution. Nine of the S. Give PhoR sequences were identical and contained the Val346Ala substitution indicated above (Figure 3). The final S. Give strain (MS170160) also contained the Val346Ala in addition to a Phe261Leu substitution. While the PhoR sequences were mostly conserved within the serovars Enteritidis, Give, and Typhimurium, there were three different PhoR sequences among the S. Kentucky strains. Seven S. Kentucky strains (RM17366, FSIS1500990, NY-NY14650, CFSAN011780, NY-FSL-R8-2848, FSIS1502818, and FSIS1503823) had PhoR sequences that were identical to the S. Typhimurium strains (except for a deletion of Pro64 in NY-FSL-R8-2848). Interestingly, these seven strains belong to Sequence Type 152. Three S. Kentucky strains (45587, FDA778036 C2 1-1, and FDA744661) had the Val346Ala substitution observed in the S. Give strains. These three strains belonged to ST696, ST314, and ST314, respectively. Finally, five S. Kentucky strains (RM19068, FDA354917, 37002, 1663660, and SAL-21-VL-ON-BC-0032) contained the Val346Ala and Thr340Lys substitutions. These four strains were all ST198.

4. Discussion

While nutritional differences have long been used to distinguish species of bacteria from each other, such differences are also described within species. Strains of S. serovar Paratyphi are distinguished by isolates that can and cannot ferment D-tartrate [48]. Enrichment culture was developed to favor growth of a target organism and limit growth of potential competitors for the same nutrient sources in a culture medium, but it will always be a competition for nutrients between the microbiota present in the inoculating sample. One reason for the adoption of culture-independent methods for microbial community analysis is because of issues of bias in various enrichment media. Direct comparisons of culture-dependent and culture-independent methods reveal large differences in community breakdowns between the two methods in vastly different ecosystems, including enumeration of antibiotic-resistant bacteria in surface waters, characterizations of plant endophytic microbiota, and dental water systems [49,50,51]. However, in public health laboratories that conduct surveillance for foodborne pathogens, the isolation of individual strains is still important. Metagenomics methods for surveillance are beginning to be discussed and developed [52,53], but there is still a need for bacterial isolates for subtyping and inclusion in national databases to link foodborne illness and outbreaks [54].
In assessing factors that might affect fitness in culture media, we studied strains of four different Salmonella serovars. In pure cultures, growth kinetics were similar for the serovars in TSB and BPW. While the Enteritidis and Typhimurium strains had generation times about 10% faster than Kentucky and Give, this was within standard error. However, such slight differences might affect fitness when in mixed cultures, and this is a subject of further study. Higher levels of growth in TSB reflect the fact that TSB is a richer medium than BPW. BPW contains peptone (10 g/L), NaCl, and sodium and potassium phosphates for a phosphate concentration of approximately 36 mM [14]. Peptone is a complex medium made from the hydrolysis of animal proteins. TSB comprises peptides sourced from digests of casein (15 g/L) and soybean (5 g/L), NaCl, 14 mM of glucose, and 14 mM of potassium phosphate [14]. The protein hydrolysates serve as the sources of carbon, nitrogen, sulfur, and energy [14].
We used phenotype microarrays to determine potential nutritional differences, and there were over 150 phenotype differences noted. It is interesting that Enteritidis strains were most adept at using dipeptides as nitrogen sources. S. Enteritidis strains are often outgrown by other Salmonella serovars when in mixed serovar cultures. Deaven et al. showed that strains of serovars Enteritidis and Typhimurium were masked in 71 and 78% of samples, respectively, when they were detected in fresh water samples [9]. In samples of naturally contaminated chicken carcasses, Cox et al. revealed that serovar Enteritidis could be detected after enrichments in Buffered Peptone Water from three of the eight carcasses using the molecular method CRISPR-SeroSeq but was found through culture-dependent methods only once [7]. Enteritidis and Typhimurium strains both displayed the potential to grow on various sugars such as mannitol, xylose, rhamnose, while Typhimurium strains utilized tagatose. Growth on these sugars has been documented for Salmonella Typhimurium previously, and xylose is used in Salmonella selective plating media, including Xylose Lysine Desoxycholate Agar and Xylose Lysine Tergitol 4 Agar [55,56,57,58,59]. Addition of any of these sugars to BPW or other complex pre-enrichment media might aid in improving the growth of Enteritidis and Typhimurium strains when competing with serovar Kentucky and/or Give strains in enrichments; this will be addressed in future studies.
The number of different compounds that may be used as phosphorus sources and the better growth with the individual phosphorus sources by S. Kentucky strains is curious. Phosphorus is usually assimilated by bacteria in the form of inorganic phosphates (Pi) [60]. Phosphorus is not predominant in nature, which is why excess phosphates (as well as nitrates and sulfates) in agricultural and sewage runoff lead to eutrophication in surface waters and grasslands [61]. The concentration of phosphates used in laboratory culture media is artificial, and many bacteria are stressed or inhibited by high phosphate concentrations [62,63,64]. S. Kentucky might be better at scavenging for phosphorus in various conditions. Therefore, we thought a closer look at phosphate metabolism in these strains might be useful.
Of the nine genes assessed, the gene product that showed differences based on serovar was phoR. The PhoB/PhoR two-component regulatory signal transduction system regulates the expression of genes involved in the transport of Pi and other P-containing compounds [65]. These are not the only genes comprising the Pho regulon, which is estimated to contain at least 137 gene products in Escherichia coli [65,66]. Proteomic analysis of the S. Typhimurium strain 14028s revealed over 300 proteins that are regulated by Pi-limiting conditions, and at least 27 of these are regulated by PhoB [26]. PhoR autophosphorylates upon reduced Pi levels and transfers that phosphoryl group to PhoB, thereby activating PhoB. When Pi levels are sufficient, PhoR removes the phosphoryl group from PhoB [29]. It was reported recently that Salmonella limits Pi uptake to maintain cellular processes involving Mg2+ homeostasis. Excess Pi uptake leads to increased ATP, which in turn chelates intracellular Mg2+ [28]. Therefore, regulation of aspects of cell metabolism that control nutrient assimilation, cellular processes, and survival of stresses and transition into the stationary phase are complex. Sequence variability in a central two-component regulatory system might result in slight differences in functionality of the proteins. The increased variability in phoR compared to the other genes led to examination of more phoR gene sequences from additional Salmonella strains.
The highest phoR sequence variability was among the S. Kentucky strains. S. Kentucky is a polyphyletic serovar with at least two different subgroups [67], and our comparison of 15 S. Kentucky strains found three different predicted protein sequences for PhoR. Interestingly, S. Kentucky ST152 and ST198 each had unique PhoR substitutions. In a study of another polyphyletic serovar (S. Mississippi), Cheng et al. concluded that the O and H antigenic elements that determine serovar may be inherited from different ancestors, leading to different sequence types with the same serovar [68]. More study of the phylogenetically distinct S. Kentucky strains is needed to determine if similar independent acquisitions of antigenic elements occurred.
Nearly all of the variability in PhoR sequences was in the CA domain, i.e., the part of the protein involved in ATP binding and catalysis. This region is located between amino acid residues 267 and 431. Immediately upstream of the CA domain is the DHp domain (located between amino acid residues 193 and 267), which is responsible for dimerization and histidine phosphorylation. The autophosphorylation of PhoR involves the transfer of a phosphoryl group from the CA domain to an essential His-213 residue in the DHp domain [47]. The PhoR protein has opposing kinase and phosphatase functions in the CA and DHp domains, and it was predicted that these two activities are controlled by physical constraint and a swing movement in the protein to prevent access of the CA domain to the DHp domain [65]. Such constraint between the domains has been described in studies of the crystal structure of similar sensor histidine kinases [69,70]. The mechanism by which PhoR receives the environmental phosphate signal is unknown [46], but in Salmonella, the phoBR operon is partly responsible for maintaining intracellular phosphate levels between 1 and 10 mM [65,71,72]. Phosphate homeostasis, in turn, involves multiple functions in the cell [26,28,65]. While the different PhoR proteins in these Salmonella are obviously functional (as the operon is essential), it might be possible that sequence differences in this CA domain have variable kinetics. We could not find any studies in the literature that involved direct testing of amino acid variations in the CA domain of PhoR, but it is tempting to speculate that changes in this region might affect the efficiency of function. Because of the complexity of the interactions of the proteins involved in phosphate metabolism and transport as well as the central role phosphate homeostasis plays in cellular metabolism, small differences in efficiency might affect competition between strains when in complex culture backgrounds.
In addition to the nutritional and genomic differences reported here, there are other factors beyond culture media that could contribute to unequal enrichment of Salmonella serovars. Salmonella strains were shown to enter into a viable but non-culturable (VBNC) state when incubated on dried fruit [73]. While a cocktail of Salmonella strains was not detectable via enrichment culture from dried apple slices after 46 days of incubation, live/dead staining showed 25% and 32% of the total population were viable after storage for 50 days at 25 °C and 110 days at 4 °C, respectively [73]. Bacteria can enter the VBNC state after exposure to stresses such as temperature, sanitizers, osmotic pressure, pH changes, and nutrient starvation [74,75,76]. Similar conditions may simply lead to stressed cells that are less likely to survive enrichment [77]. Several of these conditions may be encountered when Salmonella is present as a contaminant in and on foods. Additionally, Salmonella can survive inside of protozoa, which might mask its presence and render it more resistant to disinfectants [74,78]. In enrichment culture, the target organism must compete with other microbes present on or in the foods. This natural microbiota can vary depending on several factors, including the type of food, location, time of year, or natural variation [79]. Similarly, other microorganisms present in food samples may produce compounds such as bacteriocins that kill Salmonella [80]. Concentration effects of different Salmonella serovars that may contaminate the same food product may also affect the serovars that result after the enrichment process. S. Kentucky and S. Give may dominate enrichment cultures simply because there are many more of those types in the samples being tested. Surveillance for pathogens in foods is a complex process with many variables that affect the outcome. The comparison of different Salmonella serovars in growth characteristics is the beginning of a process to determine if the methods might be altered to achieve a fuller picture of the subtypes that might be present in contaminated foods.

5. Conclusions

To assess differential retrieval of Salmonella serovars in enrichment cultures, this study tested growth of four different serovars of Salmonella and revealed similar growth kinetics for the strain in pure cultures in TSB and BPW. Nutrient utilization differences between the serovars provided a basis of experimentation for supplements to enrichment culture media to enhance recovery of serovars of clinical interest that might be missed in routine food surveillance. In-depth analysis of predicted protein sequences of genes involved in phosphorus metabolism indicated differences in phoR sequence among the serovars. Further experimentation regarding the efficiency of various PhoR protein sequences may shed light on the efficiency of nutrient metabolism. Variations in the outcomes of enrichment culture may be affected by many factors, and small populations of specific subtypes can be missed by culture-based methods. Molecular and/or immunogenic methods that detect living cells may be better suited to finding small numbers of contaminants in situations where isolated colonies are not necessary.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms11082109/s1, Table S1: Salmonella genomes used in this study; Table S2: Heat map of all 157 phenotype differences between the strains and serovars; Figure S1: Comparisons of nine protein sequences between sequenced strains in this study; Figure S2: Comparisons of nine gene sequences between sequenced strains in this study.

Author Contributions

Conceptualization, L.G.; methodology, L.G.; validation, A.A.N.; formal analysis, L.G. and A.A.N.; investigation, L.G. and A.A.N.; resources, L.G.; data curation, L.G. and A.A.N.; writing—original draft preparation, L.G.; writing—review and editing, L.G.; visualization, L.G. and A.A.N.; supervision, L.G.; project administration, L.G.; funding acquisition, L.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by funds from U.S. Department of Agriculture under Agricultural Research Service CRIS project 2030-42000-052-00D.

Data Availability Statement

Genome sequence data were deposited with the NCBI under the BioProject ID PRJNA875079. The raw sequencing data were sent to the Sequence Read Archive and are identified by BioSample accession numbers SAMN30597303–SAMN30597310. Strains were identified by their ID numbers: RM15596, RM17366, RM18387, RM18411, RM18433, RM18465, RM19068, and RM20875.

Acknowledgments

The authors thank S. Huynh, W. Miller, T. Lie, J. Leadbetter, and W. McCleary for advice and helpful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Walter, E.J.S.; Griffin, P.M.; Bruce, B.B.; Hoekstra, R.M. Estimating the number of illnesses caused by agents transmitted commonly through food: A scoping review. Foodborne Path. Dis. 2021, 18, 841–858. [Google Scholar] [CrossRef] [PubMed]
  2. Dewey-Mattia, D.; Manikonda, K.; Hall, A.J.; Wise, M.E.; Crowe, S.J. Surveillance for foodborne disease outbreaks—United States, 2009–2015. MMWR 2018, 67, 1–11. [Google Scholar] [CrossRef] [PubMed]
  3. Batz, M.B.; Richardson, L.C.; Bazaco, M.C.; Parker, C.C.; Chirtel, S.J.; Cole, D.; Golden, N.J.; Griffin, P.M.; Gu, W.; Schmitt, S.K.; et al. Recency-weighted statistical modeling approach to attribute illnesses caused by 4 pathogens to food sources using outbreak data, United States. Emerg. Infect. Dis. 2021, 27, 214–222. [Google Scholar] [CrossRef]
  4. U.S. Centers for Disease Control and Prevention. National Enteric Disease Surveillance: Salmonella Annual Report. 2016. Available online: https://www.cdc.gov/nationalsurveillance/pdfs/2016-Salmonella-report-508.pdf (accessed on 15 August 2022).
  5. European Food Safety Authority. The European Union One Health 2021 Zoonoses Report. EFSA J. 2022, 20, 7666. [Google Scholar] [CrossRef]
  6. Wang, Y.; Liu, Y.; Lyu, N.; Li, Z.; Ma, S.; Cao, D.; Pan, Y.; Hu, Y.; Huang, H.; Gao, G.F.; et al. The temporal dynamics of antimicrobial-resistant Salmonella enterica and predominant serovars in China. Natl. Sci. Rev. 2023, 10, nwac269. [Google Scholar] [CrossRef] [PubMed]
  7. Cox, N.A.; Berrang, M.E.; House, S.L.; Medina, D.; Cook, K.L.; Shariat, N.W. Population analyses reveal preenrichment method and selective enrichment media affect Salmonella serovars detected on broiler carcasses. J. Food Prot. 2019, 82, 1688–1696. [Google Scholar] [CrossRef]
  8. Pardol de la Gandara, M.; Fournet, N.; Bonifalt, L.; Lefèvre, S.; Chemaly, M.; Grastilleur, C.; Cadel-Six, S.; Fach, P.; Pignault, A.; Brisabols, A.; et al. Countrywide multi-serotype outbreak of Salmonella Bovismorbificans ST142 and monophasic Salmonella Typhimurium ST34 associated with dried pork sausages in France, September 2020 to January 2021. Euro Surveill. 2023, 28, 2200123. [Google Scholar] [CrossRef]
  9. Deaven, A.M.; Ferreira, C.M.; Reed, E.A.; See, J.R.C.; Lee, N.A.; Almaraz, E.; Rios, P.C.; Marogi, J.G.; Lamendella, R.; Zheng, J.; et al. Salmonella genomics and population analyses reveal high inter- and intraserovar diversity in freshwater. Appl. Environ. Microbiol. 2021, 87, e02594-20. [Google Scholar] [CrossRef]
  10. Gorski, L.; Liang, A.S.; Walker, S.; Carychao, D.; Aviles Noriega, A.; Mandrell, R.E.; Cooley, M.B. Salmonella enterica serovar diversity, distribution, and prevalence in public-access waters from a Central California coastal leafy green-growing region from 2011 to 2016. Appl. Environ. Microbiol. 2022, 88, e01834-21. [Google Scholar] [CrossRef]
  11. Andrews, W.H.; Wang, H.; Jacobson, A.; Ge, B.; Zhang, G.; Hammack, T. BAM Chapter 5: Salmonella. Available online: https://www.fda.gov/food/laboratory-methods-food/bam-chapter-5-salmonella (accessed on 1 June 2023).
  12. U. S. Department of Agriculture Food Safety Inpection Service. Isolation and Identification of Salmonella from Meat, Poultry, Pasteurized Egg, Carcass, and Environmental Sponges. Available online: https://www.fsis.usda.gov/sites/default/files/media_file/documents/MLG-4.13.pdf (accessed on 1 June 2023).
  13. ISO 6579-1:2017; Microbiology of the Food Chain—Horizontal Method for the Detection, Enumeration and Serotyping of Salmonella—Part1: Detection of Salmonella spp. International Organization for Standardization: Geneva, Switzerland, 2022. Available online: https://www.iso.org/standard/56712.html (accessed on 24 July 2023).
  14. Becton Dickinson. Difco & BBL Manual; Becton Dickinson and Company: Sparks, MD, USA, 2003. [Google Scholar]
  15. Gorski, L. Selective enrichment media bias the types of Salmonella enterica strains isolated from mixed strain cultures and complex enrichment broths. Pathogenesis 2012, 7, e34722. [Google Scholar] [CrossRef]
  16. Larsen, B.R.; Richardson, K.E.; Obe, T.; Schaeffer, C.; Shariat, N.W. Mixed Salmonella cultures reveal competitive advantages between strains during pre-enrichment and selective enrichment. J. Food Saf. 2021, 41, e12934. [Google Scholar] [CrossRef]
  17. Shah, D.H.; Paul, N.C.; Sischo, W.C.; Crespo, R.; Guard, J. Population dynamics and antimicrobial resistance of the most prevalent poultry-associated Salmonella serotypes. Poultry Sci. 2017, 96, 687–702. [Google Scholar] [CrossRef]
  18. Parveen, S.; Taabodi, M.; Schwarz, J.G.; Oscar, T.P.; Harter-Dennis, J.; White, D.G. Prevalence and antimicrobial resistance of Salmonella recovered from processed poultry. J. Food Prot. 2007, 70, 2466–2472. [Google Scholar] [CrossRef]
  19. McEgan, R.; Chandler, J.C.; Goodridge, L.D.; Danyluk, M.D. Diversity of Salmonella isolates from Central Florida surface waters. Appl. Environ. Microbiol. 2014, 80, 6819–6827. [Google Scholar] [CrossRef]
  20. Maurer, J.J.; Martin, G.; Hernandez, S.; Cheng, Y.; Gerner-Smidt, P.; Hise, K.; D’Angelo, M.T.; Cole, D.; Sanchez, S.; Madden, M.; et al. Diversity and persistence of Salmonella enterica strains in rural landscapes in the Southeastern United States. Pathogenesis 2015, 10, e0128937. [Google Scholar] [CrossRef] [PubMed]
  21. Callahan, M.T.; Van Kessel, J.S.; Micallef, S.A. Salmonella enterica recovery from river waters of the Maryland Eastern Shore reveals high serotype diversity and some multidrug resistance. Environ. Res. 2019, 168, 7–13. [Google Scholar] [CrossRef] [PubMed]
  22. Higgins, R.; Désilets, A.; Cantin, M.; Messier, S.; Khakhria, R.; Ismaïl, J.; Mulvey, M.R.; Daignault, D.; Caron, H. Outbreak of Salmonella Give in the province of Quebec. Can. Vet. J. 1997, 38, 780–781. [Google Scholar]
  23. Giardin, F.; Mezger, N.; Hächler, H.; Bovier, P.A. Salmonella serovar Give: An unusual pathogen causing splenic abscess. Eur. J. Clin. Microbiol. Infect. Dis. 2006, 25, 272–274. [Google Scholar] [CrossRef]
  24. Le Hello, S.; Harrois, D.; Bouchrif, B.; Sontag, L.; Elhani, D.; Guibert, V.; Zerouali, K.; Weill, F.-X. Highly drug-resistant Salmonella enterica serotype Kentucky ST198-X1: A microbiological study. Lancet Infect. Dis. 2013, 13, 672–679. [Google Scholar] [CrossRef]
  25. Haley, B.J.; Kim, S.W.; Haendiges, J.; Keller, E.; Torpey, D.; Kim, A.; Crocker, K.; Myers, R.A.; Van Kessel, J.S. Salmonella enterica serovar Kentucky recovered from human clinical cases in Maryland, USA (2011–2015). Zoonoses Pub. Health 2019, 66, 382–392. [Google Scholar] [CrossRef]
  26. Jiang, J.; Yu, K.; Qi, L.; Liu, Y.; Cheng, S.; Wu, M.; Wang, Z.; Fu, J.; Liu, X. A proteomic view of Salmonella Typhimurium in response to phosphate limitation. Proteomes 2018, 6, 19. [Google Scholar] [CrossRef] [PubMed]
  27. Groisman, E.A. The pleiotropic two-component regulatory system PhoP-PhoQ. J. Bacteriol. 2001, 183, 1835–1842. [Google Scholar] [CrossRef]
  28. Bruna, R.E.; Kendra, C.G.; Groisman, E.A.; Pontes, M.H. Limitation of phosphate assimilation maintains cytoplasmic magnesium homeostasis. Proc. Natl. Acad. Sci. USA 2021, 118, e2021370118. [Google Scholar] [CrossRef] [PubMed]
  29. Pontes, M.H.; Groisman, E.A. Protein synthesis controls phosphate homeostasis. Genes Dev. 2018, 32, 79–92. [Google Scholar] [CrossRef] [PubMed]
  30. Lamarche, M.G.; Wanner, B.L.; Crépin, S.; Harel, J. The phosphate regulon and bacterial virulence: A regulatory network connecting phosphate homeostasis and pathogenesis. FEMS Microbiol Rev. 2008, 32, 461–473. [Google Scholar] [CrossRef]
  31. Choi, S.; Jeong, G.; Choi, E.; Lee, E.-J. A dual regulatory role of the PhoU protein in Salmonella Typhimurium. mBio 2022, 13, 00811–00822. [Google Scholar] [CrossRef]
  32. Santos-Beneit, F. The Pho regulon: A huge regulatory network in bacteria. Front. Microbiol. 2015, 6, 402. [Google Scholar] [CrossRef]
  33. Bougdour, A.; Gottesman, S. ppGpp regulation of RpoS degradation via an anti-adaptor protein IraP. Proc. Natl. Acad. Sci. USA 2007, 104, 12896–12901. [Google Scholar] [CrossRef]
  34. Steinchen, W.; Zegarra, V.; Bange, G. (p)ppGpp: Magic modulators of bacterial physiology and metabolism. Front. Microbiol. 2020, 11, 2072. [Google Scholar] [CrossRef]
  35. Lanzl, M.I.; van Mastrigt, O.; Zwietering, M.H.; Abee, T.; den Besten, H.M.W. Role of substrate availability in the growth of Campylobacter co-cultured with extended spectrum beta-lactamase-producing Escherichia coli in Bolton broth. Intl. J. Food Microbiol. 2022, 363, 109518. [Google Scholar] [CrossRef]
  36. Yang, D.-D.; Alexander, A.; Kinnersley, M.; Cook, E.; Caudy, A.; Rosebrock, A.; Rosenzweig, F. Fitness and productivity increase with ecotypic diversity among Escherichia coli strains that coevolved in a simple, constant environment. Appl. Environ. Microbiol. 2020, 86, e00051-20. [Google Scholar] [CrossRef] [PubMed]
  37. Alikhan, N.F.; Zhou, Z.; Sergeant, M.J.; Achtman, M. A genomic overview of the population structure of Salmonella. PLoS Genet. 2018, 14, e1007261. [Google Scholar] [CrossRef] [PubMed]
  38. Zhou, Z.; Alikhan, N.F.; Mohamed, K.; Fan, Y.; the Agama Study Group; Achtman, M. The EnteroBase user’s guide, with case studies on Salmonella transmissions, Yersinia pestis phylogeny, and Escherichia core genomic diversity. Genome Res. 2020, 32, 138–152. [Google Scholar] [CrossRef] [PubMed]
  39. Zwietering, M.H.; Jongenburger, I.; Rombouts, F.M.; van’t Riet, K. Modeling of the bacterial growth curve. Appl. Environ. Microbiol. 1990, 56, 1875–1881. [Google Scholar] [CrossRef] [PubMed]
  40. Zhou, L.; Lei, X.-H.; Bochner, B.R.; Wanner, B.L. Phenotype microarray analysis of Escherichia coli K-12 mutants with deletions of all two-component systems. J. Bacteriol. 2003, 185, 4956–4972. [Google Scholar] [CrossRef]
  41. Tohsato, Y.; Mori, H. Phenotype profiling of single gene deletion mutants of E. coli using Biolog technology. Genome Inform. 2008, 21, 42–52. [Google Scholar] [CrossRef]
  42. Ceapa, C.; Lambert, J.; van Limpt, K.; Wels, M.; Smokvina, T.; Knol, J.; Kleerebezem, M. Correlation of Lactobacillus rhamnosus genotypes and carbohydrate utilization signatures determined by phenotype profiling. Appl. Environ. Microbiol. 2015, 81, 5458–5470. [Google Scholar] [CrossRef]
  43. Parker, C.T.; Huynh, S.; Alexander, A.; Oliver, A.S.; Cooper, K.K. Genomic characterization of Salmonella typhimurium DT104 strains associated with cattle and beef products. Pathogens 2021, 10, 529. [Google Scholar] [CrossRef]
  44. Zhang, S.; den Bakker, H.C.; Li, S.; Chen, J.; Dinsmore, B.A.; Lane, C.; Lauer, A.C.; Fields, P.I.; Deng, X. SeqSero2: Rapid and Improved Salmonella Serotype Determination Using Whole-Genome Sequencing Data. Appl. Environ. Microbiol. 2019, 85, e01746-19. [Google Scholar] [CrossRef]
  45. Aziz, R.K.; Bartels, D.; Best, A.A.; DeJongh, M.; Disz, T.; Edwards, R.A.; Formsma, K.; Gerdes, S.; Glass, E.M.; Kubal, M.; et al. The RAST server: Rapid annotations using subsystems technology. BMC Genom. 2008, 9, 75. [Google Scholar] [CrossRef]
  46. Gardner, S.G.; Miller, J.B.; Dean, T.; Robinson, T.; Erickson, M.; Ridge, P.G.; McCleary, W.R. Genetic analysis, structural modeling, and direct coupling analysis suggest a mechanism for phosphate signaling in Escherichia coli. BMC Genet. 2015, 16, S2. [Google Scholar] [CrossRef] [PubMed]
  47. Carmany, D.O.; Hollingsworth, K.; McCleary, W.R. Genetic and biochemical studies of phosphatase activity of PhoR. J. Bacteriol. 2003, 185, 1112–1115. [Google Scholar] [CrossRef]
  48. Malorny, B.; Bunge, C.; Helmuth, R. Discrimination of d-tartrate-fermenting and -nonfermenting Salmonella enterica subsp. enterica isolates by genotypic and phenotypic methods. J. Clin. Microbiol. 2003, 41, 4292–4297. [Google Scholar] [CrossRef] [PubMed]
  49. Garcia-Armisen, T.; Aznil, A.; Cornelis, P.; Chevreuil, M.; Servais, P. Identification of antimicrobial resistant bacteria in rivers: Insights into the cultivation bias. Water Res. 2013, 47, 4938–4947. [Google Scholar] [CrossRef] [PubMed]
  50. Singh, R.; Stine, O.C.; Smith, D.L.; Spitznagel, J.K., Jr.; Labib, M.E.; Williams, H.N. Microbial diversity of biofilms in dental unit water systems. Appl. Environ. Microbiol. 2003, 69, 3412–3420. [Google Scholar] [CrossRef] [PubMed]
  51. Anguita-Maeso, M.; Olivares-Garcia, C.; Haro, C.; Imperial, J.; Navas-Cortés, J.A.; Landa, B.B. Culture-dependent and culture-independent characterization of the olive xylem microbiota: Effect of sap extraction methods. Front. Plant Sci. 2020, 10, 1708. [Google Scholar] [CrossRef]
  52. Uelze, L.; Becker, N.; Borowiak, M.; Busch, U.; Dangel, A.; Deneke, C.; Fischer, J.; Flieger, A.; Hepner, S.; Huber, I.; et al. Toward an integreated genome-based surveillance of Salmonella enterica in Germany. Front. Microbiol. 2021, 12, 626941. [Google Scholar] [CrossRef]
  53. Forbes, J.D.; Knox, N.C.; Ronholm, J.; Pagotto, F.; Reimer, A. Metagenomics: The next culture-independent game changer. Front. Microbiol. 2017, 8, 1069. [Google Scholar] [CrossRef]
  54. Carleton, H.A.; Besser, J.; Williams-Newkirk, A.J.; Huang, A.; Trees, E.; Gerner-Smidt, P. Metagenomic approaches for public health surveillance of foodborne infections: Opportunities and challenges. Foodborne Path. Dis. 2019, 16, 474–479. [Google Scholar] [CrossRef]
  55. Miller, R.G.; Tate, C.R. Xylose-lysine-tergitol 4: An improved selective agar medium for the isolation of Salmonella. Poult. Sci. 1991, 70, 2429–2432. [Google Scholar] [CrossRef]
  56. Mortlock, R.P.; Old, D.C. Utilization of D-xylose by wild-type strains of Salmonella typhimurium. J. Bacteriol. 1979, 137, 173–178. [Google Scholar] [CrossRef]
  57. Berkowitz, D. D-mannitol utilization in Salmonella typhimurium. J. Bacteriol. 1971, 105, 232–240. [Google Scholar] [CrossRef] [PubMed]
  58. Taylor, W.I.; Schelhart, D. Isolation of Shigellae. Appl. Microbiol. 1968, 16, 1383–1386. [Google Scholar] [CrossRef] [PubMed]
  59. Akhy, M.T.; Brown, C.M.; Old, D.C. L-rhamnose utilisation in Salmonella typhimurium. J. Appl. Microbiol. 1984, 56, 269–274. [Google Scholar] [CrossRef]
  60. Wanner, B.L. Phosphorus assimilation and its control of gene expression in Escherichia coli. In The Molecular Basis of Bacterial Metabolism; Hauska, G., Thauer, R.K., Eds.; Springer: Berlin/Heidelberg, Germany, 1990. [Google Scholar]
  61. Correll, D.L. The role of phosphorus in the eutrophication of receiving waters: A review. J. Environ. Qual. 1998, 27, 261–266. [Google Scholar] [CrossRef]
  62. Conrad, R.; Klose, M.; Claus, P. Phosphate inhibits acetotrophic methanogenesis on rice roots. Appl. Environ. Microbiol. 2000, 66, 828–831. [Google Scholar] [CrossRef]
  63. Zhang, Z.-Z.; Hu, H.-Y.; Xu, J.-J.; Shi, Z.-J.; Deng, R.; Ji, Z.-Q.; Shi, M.-L.; Jin, R.-C. Effects of inorganic phosphate on a high-rate anammox system: Performance and microbial community. Ecol. Engin. 2017, 101, 201–210. [Google Scholar] [CrossRef]
  64. Korkeala, H.; Pekkanen, T.J. The effect of pH and postassium phosphate buffer on the toxicity of cadmium for bacteria. Acta Vet. Scand. 1978, 19, 93–101. [Google Scholar] [CrossRef]
  65. Gardner, S.G.; McCleary, W.R. Control of the phoBR regulon in Escherichia coli. EcoSal Plus 2019, 8, ESP-0006-2019. [Google Scholar] [CrossRef]
  66. VanBogelen, R.A.; Olson, E.R.; Wanner, B.L.; Neidhardt, F.C. Global analysis of proteins synthesized during phophorus restriction in Escherichia coli. J. Bacteriol. 1996, 178, 4344–4366. [Google Scholar] [CrossRef]
  67. Vosik, D.; Tewari, D.; Dettinger, L.; M’ikanatha, N.M.; Shariat, N.W. CRISPR typing and antibiotic resistance correlates with polyphyletic distribution in human isolates of Salmonella Kentucky. Foodborne Path. Dis. 2018, 15, 101–108. [Google Scholar] [CrossRef] [PubMed]
  68. Cheng, R.A.; Orsi, R.H.; Wiedmann, M. Phylogenetic clustering suggests that distinct clades of Salmonella enterica serovar Mississippi are endemic in Australia, the United Kingdom, and the United States. mSphere 2021, 6, e00485-21. [Google Scholar] [CrossRef] [PubMed]
  69. Marina, A.; Waldburger, C.D.; Hendrickson, W.A. Structure of the entire cytoplasmid portion of a sensor histidine-kinase protein. EMBO J. 2005, 24, 4247–4259. [Google Scholar] [CrossRef] [PubMed]
  70. Wang, C.; Sang, J.; Wang, J.; Su, M.; Downey, J.S.; Wu, Q.; Wang, S.; Cai, Y.; Xu, X.; Wu, J.; et al. Mechanistic insights revealed by the crystal structure of a histidine kinase with signal transducer and sensor domains. PLoS Biol. 2013, 11, e1001493. [Google Scholar] [CrossRef]
  71. Xavier, K.B.; Kossman, M.; Santos, H.; Boos, W. Kinetic analysis by in vivo 31P nuclear magnetic resonance of internal Pi during the uptake of sn-glycerol-3-phosphate by the pho regulon-dependent Ugp system and the glp regulon-dependent GlpT system. J. Bacteriol. 1995, 177, 699–704. [Google Scholar] [CrossRef]
  72. Shulman, R.G.; Brown, T.R.; Ugubil, K.; Ogawa, S.; Cohen, S.M.; den Hollander, J.A. Cellular applications of 31P and 13C nuclear magnetic resonance. Science 1979, 205, 160–166. [Google Scholar] [CrossRef]
  73. Jayeola, V.; Farber, J.M.; Kathariou, S. Induction of the viable-but-nonculturable state in Salmonella contaminating dried fruit. Appl. Environ. Microbiol. 2022, 88, e01733-21. [Google Scholar] [CrossRef]
  74. Liu, H.; Whitehouse, C.A.; Li, B. Presence and persistence of Salmonella in water: The impact on microbial quality of water and food safety. Front. Public Health 2018, 6, 159. [Google Scholar] [CrossRef]
  75. Purevdorj-Gage, L.; Nixon, B.; Bodine, K.; Xu, Q.; Doerrler, W.T. Differential effect of food sanitizers on formation of viable but nonculturable Salmonella enterica in poultry. J. Food Prot. 2018, 81, 386–393. [Google Scholar] [CrossRef]
  76. Zhang, Y.; Liao, X.; Feng, J.; Liu, D.; Chen, C.; Ding, T. Inducation of viable but nonculturable Salmonella spp. in liquid eggs by mild heat and subsequent resuscitation. Food Microbiol. 2023, 109, 104127. [Google Scholar] [CrossRef]
  77. Wu, V.C.H. A review of microbial injury and recovery methods in food. Food Microbiol. 2008, 25, 735–744. [Google Scholar] [CrossRef] [PubMed]
  78. Brandl, M.T.; Rosenthal, B.M.; Haxo, A.F.; Berk, S.G. Enhanced survival of Salmonella enterica in vesicles released by a soilborne Tetrahymena species. Appl. Environ. Microbiol. 2005, 71, 1562–1569. [Google Scholar] [CrossRef] [PubMed]
  79. Brandl, M.T.; Mammel, M.K.; Simko, I.; Richter, T.K.S.; Gebru, S.T.; Leonard, S.R. Weather factors, soil microbiome, and bacteria-fungi interactions as drivers of epiphytic phyllosphere communities of romaine lettuce. Food Microbiol. 2023, 113, 104260. [Google Scholar] [CrossRef] [PubMed]
  80. Ottesen, A.R.; Gonzalez, A.; Bell, R.; Arce, C.; Rideout, S.; Allard, M.; Evans, P.; Strain, E.; Musser, S.; Knight, R.; et al. Co-enriching microflora associated with culture based methods to detect Salmonella from tomato phyllosphere. Pathogenesis 2013, 8, e73079. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Growth curves of the Salmonella serovars in TSB and BPW. Growth data from the Salmonella strains were averaged and are displayed according to serovar. (A) Log CFU/mL over 24 h; (B) the same graph as in (A) but showing only the first 10 h to highlight the lag and log phases. Enteritidis is shown with black circles, Give is shown with orange squares, Kentucky is shown with blue triangles, and Typhimurium is shown with pink diamonds. TSB curves are solid lines, and BPW curves are dashed lines.
Figure 1. Growth curves of the Salmonella serovars in TSB and BPW. Growth data from the Salmonella strains were averaged and are displayed according to serovar. (A) Log CFU/mL over 24 h; (B) the same graph as in (A) but showing only the first 10 h to highlight the lag and log phases. Enteritidis is shown with black circles, Give is shown with orange squares, Kentucky is shown with blue triangles, and Typhimurium is shown with pink diamonds. TSB curves are solid lines, and BPW curves are dashed lines.
Microorganisms 11 02109 g001
Figure 2. Heat map of 41 phenotype differences between the serovars. Strains were analyzed via phenotype microarray at least three times and background-subtracted, and the averages for strains were grouped into serovars. The numbers (OU range 0–500) are a measure of intensity of redox dye that reflects usage of the compound tested. The type of pre-configured microplates tested, the phenotype tested, and the chemical compound in the well are indicated. Cells are colored to indicate ranges of OU values (<20 OUs: lavender; 20–100 OUs: yellow; >100 OUs: blue).
Figure 2. Heat map of 41 phenotype differences between the serovars. Strains were analyzed via phenotype microarray at least three times and background-subtracted, and the averages for strains were grouped into serovars. The numbers (OU range 0–500) are a measure of intensity of redox dye that reflects usage of the compound tested. The type of pre-configured microplates tested, the phenotype tested, and the chemical compound in the well are indicated. Cells are colored to indicate ranges of OU values (<20 OUs: lavender; 20–100 OUs: yellow; >100 OUs: blue).
Microorganisms 11 02109 g002
Figure 3. Alignments of predicted protein sequences of PhoB, PhoR, and PstA. Protein sequences were predicted from DNA sequences of the strains indicated. Strains are labeled with the serovar and their RM number. Disagreements in protein sequence are indicated by a gap in the green consensus bar and are shown by a line in the corresponding site in the individual sequences. Green line indicates consensus between the strains. Different color vertical lines in the individual strain sequences indicate amino acid identities between strains.
Figure 3. Alignments of predicted protein sequences of PhoB, PhoR, and PstA. Protein sequences were predicted from DNA sequences of the strains indicated. Strains are labeled with the serovar and their RM number. Disagreements in protein sequence are indicated by a gap in the green consensus bar and are shown by a line in the corresponding site in the individual sequences. Green line indicates consensus between the strains. Different color vertical lines in the individual strain sequences indicate amino acid identities between strains.
Microorganisms 11 02109 g003
Figure 4. Alignment of PhoR amino acid sequences from 46 Salmonella strains. Green line indicates consensus between the strains. Different color vertical lines in the individual strain sequences indicate amino acid identities between strains.
Figure 4. Alignment of PhoR amino acid sequences from 46 Salmonella strains. Green line indicates consensus between the strains. Different color vertical lines in the individual strain sequences indicate amino acid identities between strains.
Microorganisms 11 02109 g004
Table 1. Salmonella strains used in this study.
Table 1. Salmonella strains used in this study.
Strain IDSerovarSequence TypeSource, Year
RM18465 *,1Give654Water/sediment, 2014
RM20875 *Give158Water/sediment, 2013
RM15596 *Enteritidis814Water/sediment, 2011
RM15890EnteritidisND 2Water/sediment, 2011
RM16950EnteritidisNDWater/sediment, 2011
RM18433 *Enteritidis11Water/sediment, 2013
RM20787EnteritidisNDWater/sediment, 2015
RM1316KentuckyNDGround beef, 1997
RM7890KentuckyNDGround chicken, 2008
RM17366 *Kentucky152Water/sediment, 2013
RM19068 *Kentucky198Water/sediment, 2015
RM17930TyphimuriumNDWater/sediment, 2013
RM16783TyphimuriumNDWater/sediment, 2012
RM18411 *Typhimurium19 3Water/sediment, 2014
RM18387 *Typhimurium19 3Water/sediment, 2014
1 Strains marked with an asterisk (*) were subjected to whole genome sequencing. 2 ND: not determined. 3 While these are the same sequence type, their whole genome MLST profile is different, indicating they are different strains.
Table 2. Genes involved in phosphate metabolism in Salmonella.
Table 2. Genes involved in phosphate metabolism in Salmonella.
Gene/OperonIdentityFunction
phoBRTwo-component histidine kinase response regulatory proteinsRespond to environmental Pi levels and transcriptional regulation of genes involved in phosphate transfer and utilization
pstSCABPhosphate ABC transporterHigh-affinity Pi transporter
phoUPhosphate-specific transport system accessoryBinds to and controls function of PstSCAB transfer system Binds Mg2+ and Mn2+
rpoSSigma factor SCellular responses to stresses
iraPRpoS adaptor proteinInhibits RpoS proteolysis during Pi starvation
rrsBAdaptor proteinRegulator of RpoS degradation
Table 3. Growth kinetics of Salmonella serovars in TSB and BPW at 37 °C.
Table 3. Growth kinetics of Salmonella serovars in TSB and BPW at 37 °C.
MediumSerovarLag Phase Length (h)Generation Time (h)Maximum CFU/mL
TSB 1Enteritidis3.00.511.32 × 1010
Give2.30.582.83 × 1010
Kentucky2.20.563.23 × 1010
Typhimurium2.20.493.11 × 1010
BPW 1Enteritidis2.30.506.69 × 109
Give2.20.576.93 × 109
Kentucky2.20.578.07 × 109
Typhimurium2.30.517.95 × 109
1 TSB: Trypticase Soy Broth; BPW: Buffered Peptone Water.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gorski, L.; Noriega, A.A. Comparison of Phenotype Nutritional Profiles and Phosphate Metabolism Genes in Four Serovars of Salmonella enterica from Water Sources. Microorganisms 2023, 11, 2109. https://doi.org/10.3390/microorganisms11082109

AMA Style

Gorski L, Noriega AA. Comparison of Phenotype Nutritional Profiles and Phosphate Metabolism Genes in Four Serovars of Salmonella enterica from Water Sources. Microorganisms. 2023; 11(8):2109. https://doi.org/10.3390/microorganisms11082109

Chicago/Turabian Style

Gorski, Lisa, and Ashley Aviles Noriega. 2023. "Comparison of Phenotype Nutritional Profiles and Phosphate Metabolism Genes in Four Serovars of Salmonella enterica from Water Sources" Microorganisms 11, no. 8: 2109. https://doi.org/10.3390/microorganisms11082109

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop