Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 16 November 2023
Sec. Antimicrobials, Resistance and Chemotherapy
This article is part of the Research Topic Salmonella spp.- Transmission, Pathogenesis, Host-pathogen interaction, Prevention and Treatment View all 20 articles

Adaptive laboratory evolution of Salmonella enterica in acid stress

  • 1Department of Microbiology, University of Massachusetts, Amherst, MA, United States
  • 2Department of Food Science, University of Massachusetts, Amherst, MA, United States

Introduction: Adaptive laboratory evolution (ALE) studies play a crucial role in understanding the adaptation and evolution of different bacterial species. In this study, we have investigated the adaptation and evolution of Salmonella enterica serovar Enteritidis to acetic acid using ALE.

Materials and methods: Acetic acid concentrations below the minimum inhibitory concentration (sub-MIC) were used. Four evolutionary lineages (EL), namely, EL1, EL2, EL3, and EL4, of S. Enteritidis were developed, each demonstrating varying levels of resistance to acetic acid.

Results: The acetic acid MIC of EL1 remained constant at 27 mM throughout 70 days, while the MIC of EL2, EL3, and EL4 increased throughout the 70 days. EL4 was adapted to the highest concentration of acetic acid (30 mM) and demonstrated the highest increase in its MIC against acetic acid throughout the study, reaching an MIC of 35 mM on day 70. The growth rates of the evolved lineages increased over time and were dependent on the concentration of acetic acid used during the evolutionary process. EL4 had the greatest increase in growth rate, reaching 0.33 (h−1) after 70 days in the presence of 30 mM acetic acid as compared to EL1, which had a growth rate of 0.2 (h−1) after 70 days with no exposure to acetic acid. Long-term exposure to acetic acid led to an increased MIC of human antibiotics such as ciprofloxacin and meropenem against the S. enterica evolutionary lineages. The MIC of ciprofloxacin for EL1 stayed constant at 0.016 throughout the 70 days while that of EL4 increased to 0.047. Bacterial whole genome sequencing revealed single-nucleotide polymorphisms in the ELs in various genes known to be involved in S. enterica virulence, pathogenesis, and stress response including phoP, phoQ, and fhuA. We also observed genome deletions in some of the ELs as compared to the wild-type S. Enteritidis which may have contributed to the bacterial acid adaptation.

Discussion: This study highlights the potential for bacterial adaptation and evolution under environmental stress and underscores the importance of understanding the development of cross resistance to antibiotics in S. enterica populations. This study serves to enhance our understanding of the pathogenicity and survival strategies of S. enterica under acetic acid stress.

1. Introduction

Adaptive laboratory evolution (ALE) is a technique used to study the evolution of organisms in response to selective pressures over time. In laboratory settings, ALE is typically performed with organisms possessing a short generation time such as bacteria. Through ALE, bacteria are subjected to increasingly harsh conditions over several generations, allowing them to adapt and evolve in response to stress (Dragosits and Mattanovich, 2013). The technique of ALE has proven to be highly effective in providing a deeper understanding of the genetic mechanisms and processes involved in bacterial adaptation, as well as in facilitating research on the evolution of bacterial populations over time (Gresham and Dunham, 2014; LaCroix et al., 2017; Sandberg et al., 2019).

Salmonella enterica is a pathogenic bacterium that can cause serious foodborne illnesses in humans (Spector and Kenyon, 2012; Ren et al., 2015). It frequently encounters acidic conditions in the external environment and inside hosts, and has evolved strategies to survive and adapt to acid stress (Álvarez-Ordóñez et al., 2012). The use of organic acids for fermentation and flavoring in the food industry provides opportunities for contaminating S. enterica to become acid-adapted (Coban, 2020). Additionally, the use of acidified solutions for cleaning and sanitization in households and industries creates environments conducive to the development of acid adaptation and acid tolerance response in surviving S. enterica (Wang et al., 2019). In this study, we have used an evolutionary approach to understand the adaptation of S. enterica to acid stress by subjecting it to sub-inhibitory concentrations of acetic acid over more than 1,000 bacterial generations. Subsequently, we were able to analyze the bacterial genome and identify the genetic mutations resulting from growth under acidic conditions.

Adaptive laboratory evolution of S. enterica in environmental stresses is a relatively unexplored field, making our study pioneering and essential. The results of the study have implications for food safety since understanding the evolution of Salmonella populations with reduced acid susceptibility can impact the efficacy of food preservation methods that rely on acidification. Moreover, prior adaptation to acid stress has been reported to provide cross-protection against multiple other forms of stresses that are frequently encountered by the pathogen in industrial settings (Xu et al., 2008). In this study, we evaluate whether continued exposure to acetic acid provides cross-resistance in S. enterica populations against human antibiotics. Understanding the relationship between acetic acid exposure and antibiotic resistance is crucial in combating the global issue of antibiotic resistance. Analyzing the acid-adaptive capabilities of S. enterica will contribute to our knowledge of the dynamics of bacterial evolution in response to stress. Overall, this study emphasizes the importance of understanding the genetic basis of bacterial adaptation to acid stress which will help in the development of effective preventative strategies.

2. Materials and methods

2.1. Bacterial isolates and growth media

The bacterial isolate used in this study is Salmonella enterica subsp. enterica serovar Enteritidis (ATCC BAA-1045, phage type 30). This bacterial strain was originally isolated from a salmonellosis outbreak in raw almonds (Brandl et al., 2008). It has been reported to survive long-term desiccation and nutrient starvation stress, which makes it an interesting target to study its adaptation to acid stress (Deng et al., 2012). The bacterial strain was maintained as frozen stocks in a −80°C freezer in trypticase soy broth (TSB; Sigma-Aldrich) supplemented with 15% glycerol. Frozen cultures of S. enterica were revived by streaking on tryptic soy agar (TSA; Sigma-Aldrich) and incubated at 37°C for 18 to 20 h to create a bacterial lawn. The bacteria were further streaked on trypticase soy agar (TSA, Sigma-Aldrich). Single colonies of bacteria were then inoculated in tryptic soy broth (TSB, BD diagnostic systems) overnight for 18 h.

2.2. Quantification of MIC (minimum inhibitory concentration) of acetic acid

The MIC of acetic acid against the evolutionary lineages of S. Enteritidis was determined using the broth dilution method (Wiegand et al., 2008) with minor modifications. In brief, in 50 mL tubes, a bacterial inoculum of 107–108 CFU/ml was added in 20 mL TSB containing acetic acid concentrations from 24 to 34 mM. The acetic acid concentration in each tube was incremented by 1 mM. The tubes were incubated under shaking conditions at 37°C for ∼18 h. TSB with 200 mM acetic acid and TSB without any acetic acid were also inoculated with the bacteria and served as negative and positive controls, respectively. The OD600 was quantified after 20 h, and OD600 ≤ 0.1 was considered an inhibition of growth. To determine the MIC or minimum bactericidal concentration (MBC) of acetic acid, inoculum was taken from the tubes with OD600 ≤ 0.1 and plated on TSA plates. Recovery of growth was an indication that the acetic acid concentration was MIC. This experiment was repeated in triplicates for each of the replicates of all the evolutionary lineages. The minimum inhibitory concentration of acetic acid of the evolutionary lineages was determined every 5 days after the initiation of the ALE study.

2.3. Adaptive laboratory evolution of S. Enteritidis in acetic acid

At the beginning of the adaptive evolutionary study, Salmonella enterica subsp. enterica serovar Enteritidis (referred to as WT) was grown in TSB overnight. The overnight S. Enteritidis culture was used to quantify the MIC of acetic acid [explained in section 2.2. Quantification of MIC (minimum inhibitory concentration) of acetic acid]. The MIC of WT S. Enteritidis was quantified to be 27 mM. The WT culture was grown overnight and after 20 h, the bacterial optical density of the overnight culture was adjusted, and inoculum of 107–108CFU/mL was added to three 50 mL conical tubes containing 20 mL of TSB. These tubes were labeled as EL1a, EL1b, and EL1c. Here, EL stands for evolutionary lineage. Inoculum from the WT was also added to three tubes containing 26 mM acetic acid in TSB that were labeled EL2a, EL2b, and EL2c. EL1 served as our control evolutionary lineage, and EL2 served as our first evolutionary lineage that was grown and adapted to sub-MIC of acetic acid. The six tubes were incubated at 37°C for 18–20 h. Following that, the OD600 of the six tubes was measured and recorded. An inoculum of 106 CFU/ml from each tube was collected and added to six new tubes for day 2. This process was repeated throughout the ALE study, and a similar inoculum was transferred to fresh tubes for each of the EL daily. The MIC of the six tubes was quantified every 5 days, and on day 20, the MIC of EL2 (all three replicates) increased to 29 mM. Hence, on day 20, inoculum from EL2a-EL2c was also added to three tubes containing 28 mM acetic acid (sub-MIC acetic acid) in TSB and labeled EL3a-EL3c. On day 30, the MIC of EL3 increased to 31 mM, and inoculum from EL3a-EL3c was added to three tubes containing 30 mM acetic acid in TSB and labeled EL4a-EL4c. The OD600 for each of the tubes was measured after growing them for 20 h at 37°C under shaking conditions. Samples from the previous day were stored at −80°C over the course of the evolutionary process. The number of bacterial generations ‘n’ was calculated using the following equation (Lee et al., 2011):

n = l o g ( N / N ) 0 / l o g ( 2 ) .

where N is the final number of cells in the conical tube after 20 h at the time of passage to the next day’s tubes. N0 is the initial number of cells that are transferred to each conical tube at the beginning of ALE for that day. The initial and final numbers of cells were estimated daily by measuring the OD600 using a spectrophotometer and using plate counts. The generation number calculation assumes that each cell is viable, the death rate is negligible, the cells are growing exponentially throughout the ALE experiment, and the cells are dividing by binary fission.

2.4. Quantification of growth rate of the S. Enteritidis evolutionary lineages

The growth rate of the bacterial lineages was quantified using previously established methods (Herring et al., 2006; Hall et al., 2013) with some modifications. After 20 h of growth of the evolutionary lineages at 37°C, the inoculum was adjusted to approximately 104CFU/ml and added inside the wells of a 96-well microtiter plate (Thermo Fisher Scientific, Waltham, MA, USA, model number 266120). Bacterial cultures collected at the end of day 2 and day 70 from the ELs were used to generate the growth curves inside the oCelloScope. For WT S. Enteritidis, EL1 and EL2, along with the bacterial inoculum, the wells also contained TSB with 26 mM acetic acid. For EL3 and EL4, the wells in the microtiter plate contained 28 mM and 30 mM acetic acid in TSB, respectively. The concentration of acetic acid added to the wells for each EL was 1 mM less than their MIC (sub-MIC acetic acid). The microtiter plates were placed inside the oCelloScope for 24 h to monitor the progression of bacterial growth. The exponential growth rate of all three replicates of each of the evolutionary lineages was measured. The background corrected absorption (BCA) algorithm and the BCA normalized algorithm of the software UniExplorer (version 10.1) were used to quantify the kinetics of bacterial growth. The kinetics of growth of the WT S. Enteritidis and all the evolutionary lineages in TSB without acetic acid stress were also monitored using the oCelloScope.

2.5. Quantification of MIC of antibiotics against the S. Enteritidis evolutionary lineages

The MIC of the following antibiotics: vancomycin, meropenem, ciprofloxacin, gentamycin, and streptomycin, were determined using MTS strips (MIC test strips, Liofilchem) (Matuschek et al., 2018; van den Bijllaardt et al., 2018). Antibiotics with different mechanisms of action were chosen to gain a thorough understanding of their effects on the evolutionary lineages. Vancomycin and meropenem inhibit bacterial cell wall biosynthesis (Yarlagadda et al., 2018), ciprofloxacin inhibits bacterial DNA synthesis, and gentamycin and streptomycin are known to inhibit bacterial protein synthesis (Athamneh et al., 2014).

The test strips for each antibiotic had a concentration range from 0.016 to 256 μg/ml. Overnight cultures of the evolutionary lineages were grown in their corresponding acetic acid concentrations, and the next day, the OD600 was adjusted to 107CFU/mL in TSB. Cotton swab applicators were dipped in the bacterial cultures and streaked on TSA to create bacterial lawns. The MTS test strips were carefully positioned at the center of the TSA plates using sterile forceps. The plates were incubated for 18–20 h at 37°C, following which the zones of inhibition around the test strips were recorded. The minimum concentration of the antibiotic, which caused inhibition of bacterial growth, was denoted as the MIC. This experiment was repeated in triplicates for each of the replicates of all the evolutionary lineages. MTS test strips have been used successfully for determining the MIC of antibiotics against different pathogens, and the strips have been reported to provide MIC values similar to those obtained using broth microdilution methods (Hakvoort et al., 2020; Koeth et al., 2022).

2.6. Bacterial whole-genome sequencing

The WT strain of S. Enteritidis and evolved lineages after day 70 of the ALE study were streaked on TSA and incubated overnight at 37°C. Genomic DNA extraction and sequencing were performed at SeqCenter (Pittsburgh, PA, USA). For the evolved lineages, after 70 days, the liquid cultures were plated on TSA, and colonies were randomly chosen for sequencing. DNA extraction was performed using Zymo DNA Miniprep (bead-beating lysis). For the WT strain, sample libraries were generated using the Illumina DNA Prep kit and IDT 10bp UDI indices. Sequencing was performed on a NextSeq 2000 sequencer, producing 151-bp paired-end reads. Low-quality reads were trimmed, and adapter sequences were removed from Illumina sequences using bcl-convert version 3.9.3. Oxford Nanopore Technologies (ONT) PCR-free ligation library preparation was also used to generate long sequence reads. ONT reads were adapted and quality-trimmed using porechop33 version 0.2.3_seqan2.1.1 (RRID:SCR_016967). Unicycler version 0.4.8 was used to produce a hybrid assembly from the Illumina and ONT reads. Assembly quality was assessed using QUAST version 5.0.2 (Gurevich et al., 2013). The NCBI Nucleotide BLAST database was then used to characterize the complete genome assembly and identify plasmid sequences. Gene models were predicted with functional annotations using Prokka version 1.14.5 with default parameters + ‘–rfam’ (Tatusova et al., 2016; Haft et al., 2018). Culturing and DNA extraction of the evolved lineages were performed using the same procedure as the WT, but 151-bp paired-end libraries were generated using only Illumina sequencing as described above.

2.7. Genome assembly and identification of polymorphisms

All software mentioned below was used with default parameters unless specified otherwise. BWA version 0.7.15 (RRID:SCR_010910) was used to map the evolved lineages resequencing data against the ancestral reference genome (Li, 2013). Samtools (v 1.14) (RRID:SCR_002105) was used to index the sorted BAM files. Bamaddrg was used to add read groups to the indexed and sorted BAM files (Danecek et al., 2021). Joint genotyping of the evolved strains was then performed using Freebayes (v 1.3.1) (RRID:SCR_010761) (Garrison and Marth, 2012). GATK (v. 4.0.6) (RRID:SCR_001876) was used to convert the resulting VCF files into a table format (McKenna et al., 2010). A SnpEff database was built for the S. Enteritidis BAA-1045 WT reference genome, and SNP annotation prediction was performed in the evolved strains using SnpEff (v. 4.1) (Cingolani et al., 2012). A second variant calling software, Breseq (v. 0.35.4) (RRID:SCR_010810), was used by SeqCenter to align and compare evolved lineage sequence reads to the ancestor (Deatherage and Barrick, 2014). Assemblies were generated from the Illumina reads of each evolved lineage using SPAdes (v. 3.13.1) (RRID:SCR_000131) with the parameters, “-k 21,33,55 –careful” (Anton Bankevich et al., 2012). Large-Scale Blast Score Ratio (LS-BSR) was used to investigate gene presence/absence mutations using default parameters, along with additional parameters “-b blastn -c cdhit” (Sahl et al., 2014). Read-depth analysis was performed using Samtools (v. 1.14).

2.8. Statistical analyses

All bacterial growth measurements were biologically triplicated, and their differences were examined by a two-tailed t- test assuming unequal variance (Welch’s t- test). GraphPad Prism was used to generate all the graphs in this project. Statistical significance was calculated with GraphPad Prism using two-way ANOVA. A p-value of ≤ 0.05 was considered statistically significant.

3. Results

3.1. Experimental design of the ALE study

Salmonella Enteritidis was adapted in different concentrations of acetic acid in TSB over 70 days. Figure 1 outlines the evolution of four lines of S. enterica in acetic acid, which we are calling “evolutionary lineages” (EL). The MIC of acetic acid for wild-type (WT) S. Enteritidis was 27 mM at the beginning of the experiment (Figure 2). On day 0, two ELs were developed from WT S. Enteritidis. In the first EL, EL1, S. Enteritidis was grown in TSB without any acid stress, and in the second EL, EL2, S. Enteritidis was grown in 26 mM acetic acid (sub-MIC) in TSB. Each EL was grown in triplicates at 37°C for 18 h daily before they were transferred into fresh media as described in the section “2. Materials and methods”. The MIC of acetic acid of EL1 and EL2 was measured every 5 days, and the OD600 was monitored daily before transfer. After 20 ALE days, the MIC of acetic acid for EL2 rose to 29 mM, and a new lineage, EL3, was created from EL2 (Figure 1). An inoculum of 108CFU/ml was transferred from the three replicates of EL2 to generate three replicates of EL3. EL3 was grown in 28 mM acetic acid for the remainder of the experiment. Here, ALE days refer to the number of days that have passed since the initiation of the ALE study. After 30 ALE days, the MIC of EL3 increased to 31 mM, and a new EL, EL4, was initiated (Figure 1). EL4 was grown in 30 mM acetic acid in TSB with daily transfers in fresh media. The number of generations for each EL was quantified from OD600 values and is shown in Figure 1. After 70 days, the adaptive evolution process was halted, and all the ELs were frozen at −80°C. The MIC of acetic acid of EL2–EL4 was found to increase over time throughout the duration of the ALE study (Figure 2). The MIC of acetic acid in EL1 remained constant throughout the study, and the highest increase in the MIC of acetic acid was observed in EL4.

FIGURE 1
www.frontiersin.org

Figure 1. Adaptive evolution of S. Enteritidis evolutionary lineages (EL) with daily transfers in acetic acid. The MIC of acetic acid of wild-type S. Enteritidis was quantified to be 27 mM, and hence, EL2 was grown in 26 mM acetic acid (sub-MIC) while EL1 was grown in no acid stress (control group). The MIC of EL1 and EL2 was quantified every 5 days. On day 20, the MIC of EL2 increased to 29 mM. So, we started EL3 from EL2, and it was grown in 28 mM acetic acid with daily transfers for the remainder of the ALE study. On day 30, the MIC of EL3 went up to 31 mM; hence, we started EL4 from EL3, and EL4 was grown in 30 mM acetic acid with daily transfers until day 70 of the ALE study. The OD600 for each of the EL was quantified every day before and after daily transfers, and the values were used to quantify the number of generations. The number of generations of each of the ELs has been shown in green.

FIGURE 2
www.frontiersin.org

Figure 2. Change in acetic acid MIC of the ELs of S. Enteritidis. The MIC of acetic acid of the evolutionary lineages was measured every 5 days until day 70 of the ALE study. Bacterial cultures (1 mL) were collected from the three replicates of EL1–EL4 every 5 days to quantify the MIC of acetic acid. The arrows indicate the initiation of EL3 and EL4 from EL2 and EL3, respectively. The evolutionary lineage EL3 was initiated from EL2 (on ALE day 20) when the acetic acid MIC of EL2 increased for the first time and a similar process was used for the initiation of EL4 (on ALE day 30).

3.2. Effect of ALE on the growth rate of the adapted evolutionary lineages

The exponential growth rate of the adapted ELs was evaluated in the presence and absence of acetic acid using the oCelloScope (Figure 3 and Table 1). The oCelloScope has previously been used to study changes in growth rates after treatment with antimicrobial compounds (Fredborg et al., 2013; Canali et al., 2018; Ghoshal et al., 2022b). The growth rates were quantified on ALE days 2 for EL1 and EL2 and for ALE days 22 and 32 for EL3 and EL4, respectively. ALE days 22 and 32 were 2 days after the initial transfer of EL3 and EL4, respectively. Growth rates of all ELs were calculated from the exponential growth phase of the growth curves. The growth curves of EL1 and EL4 generated using the oCelloScope are shown in Supplementary Figure 1.

FIGURE 3
www.frontiersin.org

Figure 3. Exponential growth rate of EL1 (A) quantified on ALE days 2 and 70 and EL4 (B) quantified on ALE days 32 and 70 of S. Enteritidis evolutionary lines. ALE day 32 was the second day after the initiation of EL4 from EL3. p-value of 0.05 was considered to be statistically significant using two-tailed t-test. (***p < 0.001, ns, non-significant, = p > 0.05).

TABLE 1
www.frontiersin.org

Table 1. Growth rates of adapted evolutionary lines.

No significant differences were observed in the growth rate after days 2 and 70 for EL1 (p > 0.05) (Figure 3A), while for EL4, the growth rate in 30 mM acetic acid on day 70 was significantly higher than that on day 32 (p < 0.001) (Figure 3B). However, no significant differences in growth rates in the absence of acetic acid were observed in either EL1 after days 2 and 70 or in EL4 after days 32 and 70 (p > 0.05). Additionally, the exponential growth rates of EL2 and EL3 have been represented graphically in Supplementary Figure 2. The exponential growth rate of all the ELs 2 days after their initiation and after 70 ALE days in the presence and absence of acetic acid is shown in Table 1. The results indicate that the growth rates on ALE day 70 were significantly higher than those on days 2, 22, and 32 for EL2 (p < 0.01), EL3 (p < 0.01), and EL4 (p < 0.001), respectively. On day 70, EL4 had the highest growth rate in 30 mM acetic acid followed by EL3 in 28 mM acetic acid and EL2 in 26 mM acetic acid. The lowest growth rate was observed for EL1 in the presence of 26 mM acetic acid after day 2, and no change in its growth rate was observed after day 70 (p > 0.05).

3.3. Effect of ALE on the MIC of human antibiotics in the evolutionary lineages

The MIC of ciprofloxacin, gentamycin, meropenem, streptomycin, and vancomycin was quantified in all the ELs using MTS (MIC test strips). MIC values of the antibiotics were quantified every 10 days starting from 2 days after the initiation of the ELs and ending after ALE day 70. The MIC of the antibiotics tested against EL1 and EL4 is shown in Figure 4. There were no significant differences in the MIC values of the antibiotics against EL1 between days 2 and 70. However, there was a significant difference in the MIC of ciprofloxacin (p < 0.05), gentamycin (p < 0.01), and meropenem (p < 0.001) against EL4 between days 32 and 70. Significant differences were also observed in EL3 against ciprofloxacin (p < 0.001) between ALE day 22 (2 days after the initiation of the evolutionary lineage) and day 70 of the ALE study. The MIC values of antibiotics against all the ELs are reported in Supplementary Table 1.

FIGURE 4
www.frontiersin.org

Figure 4. Change in MIC of human antibiotics against the S. Enteritidis evolutionary lineages EL1 (quantified on ALE days 2 and 70) and EL4 (quantified on ALE days 32 and 70). The antibiotics tested are meropenem and ciprofloxacin (A), gentamycin and vancomycin (B), and streptomycin (C). EL1 was grown and adapted in TSB without stress throughout the ALE study, and EL4 was adapted in 30 mM acetic acid with initial transfer from EL3 on ALE day 30. p-value of 0.05 was considered to be statistically significant using two-tailed t-test. (***p < 0.001, ns, non-significant, = p > 0.05).

3.4. Comparative genomic analysis of acetic acid-evolved lineages

Genomic analysis of the evolutionary lineages was conducted after 70 ALE days to investigate the genetic determinants of the increased MIC against acetic acid and the human antibiotics tested. We generated a complete genome assembly of the WT S. Enteritidis using Illumina and Oxford Nanopore Technologies (ONT) sequencing. WT S. Enteritidis strain BAA-1045 whole-genome Illumina and Oxford Nanopore Technologies sequencing data are available for wild-type and evolved lineages (Illumina only) under the NCBI BioProject PRJNA1007881. The genome was assembled into four contigs consisting of the circular genome (∼4.67 Mb) and three plasmids with approximate sizes of 60 kb (S. Enteritidis strain RM2968 plasmid pRM2968-2), 56 kb (S. Enteritidis strain 56-3991 plasmid pSE56-3991), and 6.2 kb (E. coli strain RHB18-C20 plasmid pRHB18-C20_4) based on the best hits (e-value < 0.01, query coverage >99%, and percent identity >99.9%) against the NCBI Nucleotide Collection (nt) database (Johnson et al., 2008).

This high-quality genome assembly and annotation was used as a reference to map Illumina whole-genome sequencing data for each EL. The mutations observed within each evolved lineage are detailed in Table 2. We identified 8, 10, 3, and 5 base substitutions in EL1, EL2, EL3, and EL4, respectively. Base substitutions were observed in a total of 15 different genes (10 missense; 3 synonymous; 2 frameshift). Notably, missense mutations were observed in the histidine kinase gene, phoQ, among lineages EL2, EL3, and EL4. Moreover, missense mutations were observed in EL2 within phoP, which is the regulatory protein of phoQ. For the PhoP protein, mutations were observed in the DNA-binding domain, while in the PhoQ protein, mutations occurred in the HAMP domain and the histidine kinase domain (Figure 5).

TABLE 2
www.frontiersin.org

Table 2. Mutations in evolutionary lines after adaptation in acetic acid for 70 ALE days identified using genomic analysis.

FIGURE 5
www.frontiersin.org

Figure 5. Schematics of PhoQ and PhoP proteins. Ellipses show protein domains with their respective Interpro identifiers for PhoP (A) and PhoQ (B). Yellow symbols represent locations of missense mutations. The EL (evolutionary lineage), replicate (each EL had three replicates namely R1-R3), and amino acid substitution are provided for each site.

We used the Large-Scale Blast Ratio (LS-BSR) pipeline (Sahl et al., 2014) to compare the genetic changes in the acid-adapted evolutionary lineages after 70 days. This helped us investigate the potential gain or loss of entire coding sequences in the evolved lineages (Figure 7). Additionally, a 2,207 bp deletion was observed in replicate 1 and replicate 2 of EL3 and all three replicates of EL4 (Figure 6). Read-depth analysis was used to map the deletion of genes between ackA and hxpA in EL4 and has been compared to WT S. Enteritidis as shown in Figure 6. This deleted region in the evolved lineage is relative to the ancestral genome as the genomic region is present in the ancestral genome, but no reads map to that region in the evolved lines, although they do map to the flanking regions.

FIGURE 6
www.frontiersin.org

Figure 6. Read-depth analysis of the deletion observed between the genes ackA and hxp2 in EL4. The read depth of EL4 has been compared to the wild-type (WT) strain of S. enterica serovar Enteritidis BAA-1045. Genes ackA and hxpA encode for acetate kinase and hexitol phosphatase in S. Enteritidis, respectively, and hypo1 and hypo2 stand for hypothetical protein 1 and hypothetical protein 2, respectively.

FIGURE 7
www.frontiersin.org

Figure 7. Large-Scale Blast Score Ratio (LS-BSR) heatmap of the three replicates of the evolutionary lineages EL1-EL4 after 70 ALE days.

The BLASTX and InterPro analysis suggests that there were two hypothetical proteins within the deleted range that are homologous to genes yfbV and yfbU. The read-depth analysis shows that a fraction of ackA and hxpA has been deleted, while the two hypothetical proteins yfbV and yfbU have been deleted entirely in EL4 as compared to the WT. The gene yfbV is a cytoplasmic protein that is predicted to be involved in chromosome segregation in E. coli (Miyakoshi et al., 2019), and the yfbU gene is known to play a role in DNA damage response in E. coli (Amitai et al., 2009). Studying the boundaries of the deleted region in EL4 indicates that by examining the extent of the reduction in read depth, it is also possible to estimate the size or length of the deleted region. This information is crucial for understanding the impact of the deletion on the affected genes. In our study, since only a fraction of ackA and hxpA is deleted, it is possible that a truncated version of the proteins is being expressed.

4. Discussion

This study describes the adaptation and evolution of Salmonella Enteritidis to different concentrations of acetic acid in TSB using adaptive laboratory evolution (ALE). The increase in MIC values of EL2, EL3, and EL4 over time suggests that S. Enteritidis can adapt to sub-inhibitory concentrations of acetic acid during prolonged acid exposure.

The oCelloScope is a valuable tool for studying the effects of antimicrobial compounds on bacterial growth kinetics and can be used to investigate the adaptation of bacteria to various stress conditions (Fredborg et al., 2013, 2015; McLaughlin et al., 2017). We used it to study the effects of acetic acid stress on the exponential hourly growth rate of the evolved lineages and compared them to the growth rate of WT S. Enteritidis. Our results indicate that the growth rates of the evolved lineages in the presence of acetic acid were dependent on the concentration of acid used during the evolutionary process (Figure 3, Table 1 and Supplementary Figure 2). Since EL1 was adapted without any acid stress throughout the ALE study, there were no significant differences in the exponential growth rate of its three replicates between ALE days 2 and 70 after exposure to 26 mM acetic acid in the oCelloScope. The exponential growth rates of EL1 recovered from days 2 and 70 in the presence and absence of 26 mM acetic acid were 0.2 (h–1) and 0.4 (h–1), respectively. These values were similar to the exponential growth rate of WT S. Enteritidis in the presence and absence of 26 mM acetic acid as monitored using the oCelloScope. The growth rates of EL2–EL4 on ALE day 70 were significantly higher than the values 2 days after the initiation of the corresponding evolutionary lineages (p < 0.01). Other researchers have demonstrated that adaptation in stress environments caused an increase in the growth rate of E. coli during ALE studies (Fong et al., 2005; Cheng et al., 2014; LaCroix et al., 2015).

Next, we aimed to investigate whether long-term exposure to acetic acid in the evolutionary lineages led to cross-resistance against human antibiotics (Figure 4). The different classes of antibiotics used in this study against S. enterica have been tested in earlier studies (Geornaras and von Holy, 2001; Jones et al., 2002; Muhammad et al., 2010; Chang et al., 2021). Studies in E. coli have found that resistance against antibiotics changes over the course of ALE studies (Maeda et al., 2020; Marciano et al., 2022). E. coli strains sensitive to carbapenem were found to develop resistance against carbapenems by the 15th generation in another ALE study (Geng et al., 2022). Ciprofloxacin has historically been used as the first line of defense against S. enterica infections, especially against many multidrug-resistant (MDR) strains (Askoura and Hegazy, 2020). Carbapenems (meropenem) have also been used successfully for the treatment of S. enterica infections, especially the ones that are resistant to cefotaxime and ciprofloxacin (Jean et al., 2005). Using antibiotics with distinct mechanisms of action provides us with a more comprehensive understanding of the evolution of resistance in evolved populations.

We did not observe significant changes in the MIC of the antibiotics tested against EL1 between days 2 and 70, likely due to the lack of any selective pressure on EL1 in the ALE study. In contrast, significant increases in the MIC values of ciprofloxacin, gentamycin, and meropenem against EL 4 (p < 0.05) (Figure 4 and Supplementary Table 1) were observed after day 70. This highlights that exposure to acid stress provides cross-adaptation against antibiotic stress in S. enterica acid-evolved populations. This result is in line with previous studies demonstrating that exposure to low concentrations of stress over time in S. enterica leads to increased MIC against antibiotics (Gullberg et al., 2011; Wistrand-Yuen et al., 2018).

Even though the MIC against the various antibiotics, especially, ciprofloxacin, gentamycin, and meropenem, was found to significantly increase after 70 ALE days, their resultant MIC values still fell within the susceptible range as defined by established clinical breakpoints (FDA, 2023). Despite this, there remains a pertinent cause for concern regarding the increase in MIC values as this suggests that the ELs exposed to acid stress displayed a reduced susceptibility to the antibiotics over time. The rise in MIC values, albeit within the susceptible range, potentially signifies the presence of underlying mechanisms that contribute to reduced antibiotic effectiveness. Such alterations in MIC values hint at adaptive responses in S. enterica under acid stress, which could subsequently impact the overall efficacy of antibiotic treatment strategies. This finding also raises concerns about the potential for acidic conditions, such as those encountered during food processing, which may inadvertently contribute to the development of antibiotic-resistant bacterial populations.

The observed differences in growth rate and increase in MIC of acetic acid and antibiotics over time suggest the presence of genetic changes enabling the ELs to adapt to the selective pressure of acid. Long-term exposures to acid stress have been reported to be associated with mutations at the genomic level in E. coli (Du et al., 2020; Zeng et al., 2021), which can alter the bacterial cell structure, physiology, and metabolism leading to multidrug resistance (MDR). However, when the ELs from day 70 were desiccated on the surface of stainless-steel coupons for 18 h and subsequently exposed to 100 mM acetic acid using previously established methods (Ghoshal et al., 2022a), no significant differences between the log reduction of the ELs and WT S. Enteritidis were observed (Supplementary Figure 3). This could be due to the fact that bacterial adaptation to a specific stressor often involves the optimization of cellular mechanisms and regulatory pathways to withstand the stress within a limited concentration range (Mongold et al., 1999). In the case of acetic acid, the acid-adapted bacteria may have undergone genotypic and phenotypic changes such as alterations in membrane permeability, efflux pump activity, intracellular pH regulation, and activation of stress response systems that enable them to survive and proliferate in the presence of lower concentrations of acetic acid. However, when exposed to significantly higher concentrations, the stress exceeds the threshold of their adaptive response, leading to a loss of viability and growth inhibition.

Mutations in the outer membrane lipoproteins and membrane transporters are known to alter the binding properties of bacterial membranes which increase resistance against various antibiotics such as vancomycin and meropenem (Langsrud et al., 2004; Maldonado et al., 2016). We observed mutations in the outer membrane lipoprotein ybaY as well as mutations in the outer membrane transporter fhuA in the acid-adapted S. Enteritidis lineages, and it is possible that these mutations rendered the outer membrane less susceptible to antibiotics. Additionally, periplasmic binding proteins help gram-negative pathogens sense their environment, regulate the uptake of small molecules, and even provide antibiotic resistance in S. enterica (Blair et al., 2009). Mutations observed in the thiamine-binding periplasmic protein in EL3 and EL4 as well as in the metal-binding protein zinT in EL2 could have possibly played a role in providing antibiotic cross-resistance after acid adaptation.

In E. coli, exposure to antimicrobial compounds such as benzalkonium chloride (BAC) has been reported to induce cross-resistance against other antimicrobial compounds and antibiotics (Chapman, 2003; Langsrud et al., 2004). The acquisition of antibiotic resistance in these studies was mainly attributed to the differential expression of bacterial efflux pumps and several outer membrane proteins. Changes in the expression of efflux pumps promote resistance against several classes of antibiotics (Jibril et al., 2021; Dawan et al., 2022). A variety of Salmonella strains also differentially express enzymes that can modify and inactivate antibiotics such as streptomycin (Mengistu et al., 2020). Upregulation in the expression of several transcription factors, which are global regulators of stress, is also known to trigger multidrug resistance in E. coli (Duval and Lister, 2013). Future RNA sequencing studies on the acid-adapted lineages of S. Enteritidis will help to determine the influence of efflux pumps, membrane proteins, and transcription factors in providing cross-resistance against antibiotics. The two-component signal transduction system PhoP/Q is reported to provide general cellular defense against stressful conditions and antibiotic resistance in E. coli and S. Enteritidis (Lázár et al., 2014; Hu et al., 2023). A hypothesis is that the mutations observed in phoP and phoQ observed in our study played an important role in triggering cross-adaptation against antibiotics in the acid-evolved S. enterica lineages.

Additionally, the phenomenon of persistence is triggered in bacteria after exposure to various stressful environments which renders them resistant to future stressors such as exposure to antibiotics (Pacios et al., 2020). Persistence is known to provide E. coli with tolerance and resistance against several antibiotics (Li and Zhang, 2007), and short-term exposure to acid stress is reported to trigger persistence in S. enterica (Leyer and Johnson, 1992). We hypothesize that reduced susceptibility to the antibiotics could be in part due to the presence of persister cells in the subpopulation of the acid-evolved S. enterica lineages. It is also possible that cross-resistance to antibiotics in S. Enteritidis acid-evolved lineages is due to a combination of the factors described above.

Genomic analysis of the acetic acid-evolved ELs revealed base substitutions in 15 different genes (Table 2), including missense mutations in multiple ELs in phoQ and phoP. The PhoQ protein is known to control the transcription of a large class of genes (Prost et al., 2007), and the PhoP/Q system is triggered under conditions of low pH and in the presence of antimicrobial and cationic peptides in S. enterica (Bearson et al., 1998; Prost et al., 2007; Yuan et al., 2017). Activation of this system is known to regulate S. enterica virulence and pathogenesis (Prost et al., 2007). In addition, the PhoP/Q protein system has been reported to be involved in stress resistance and antibiotic resistance and in the regulation of envelope composition in S. enterica (Yuan et al., 2017; Shprung et al., 2021; Dawan and Ahn, 2022).

The presence of mutations in the phoP/Q genes in all three replicates of EL2, EL3, and EL4 and its absence in EL1 indicate that these mutations are a result of the adaptation to acetic acid stress. One of the mutations observed in phoQ in EL2 was in the HAMP domain of the protein (Figure 5), which is thought to play an important role in transmitting signal to the catalytic domain of the protein (Lemmin et al., 2013). Base substitutions in the protein PhoQ Asn255Ileu have been reported where an amino acid with a neutral polar R group was replaced by a neutral non-polar R group under stress conditions (Matamouros et al., 2015) which is similar to our case (Thr271Pro). Researchers found that the Asn255Ileu substitution significantly increased the activity of the PhoQ protein (over 10-fold increase in activity) and that the mutant protein was found to spend more time in the activation state than the wild-type protein (Matamouros et al., 2015). It is possible that the mutation in our system also led to a similar increase in the activity of the PhoP protein. We also observed base substitutions in the sensor histidine kinase domain of phoQ in EL3 and EL4, which serves as the catalytic domain and regulates/activates the phosphoryl-state of the protein (Lemmin et al., 2013). Additionally, we observed a mutation in EL2 in the OmpR/PhoB-type DNA-binding domain of phoP (Figure 5). In this case, arginine, which is a basic amino acid, was replaced by glutamine, which is polar in nature, and at low pH, such amino acids acquire an overall positive charge. This positively charged amino acid in DNA-binding domains of proteins aids in more effective binding with negatively charged DNA (Cherstvy, 2009), and such modifications in response regulators are often associated with differential changes in gene expression of downstream genes (Rhee et al., 2008). It is possible that the mutations in phoP/Q among the ELs contributed to the reduced susceptibility against the antibiotics causing a differential expression of stress and acid adaptation genes. Mutations in the phoPQ system are associated with increased resistance against antibiotics in other gram-negative pathogens such as Pseudomonas aeruginosa (Miller et al., 2011). In our study, we have chosen to prioritize the investigation of missense mutations within the phoP/Q genes as a primary focus, while avoiding gene deletion studies. This strategic decision is grounded in several considerations. First, the existing body of literature provides substantial evidence suggesting that missense mutations in phoP/Q genes play a pivotal role in bacterial adaptations to stress conditions, including acid stress (Prost et al., 2007; Matamouros et al., 2015; Pitt et al., 2018). Moreover, planned RNA sequencing studies in future will offer valuable insights into the dynamic changes in the expression of phoP/Q genes in response to acid stress. These findings will complement our observations of missense mutations and help establish a multifaceted understanding of the molecular mechanisms underlying the adaptation process. By focusing on missense mutations in conjunction with transcriptomic analysis, we aim to provide a comprehensive perspective on the significance of phoP/Q genes in stress adaptation without resorting to gene deletions. This will help minimize potential disruptions to the overall cellular machinery and preserve the ecological relevance of our study.

In EL2, we found a mutation in the fhuA gene which encodes for the ferrichrome outer membrane transporter. This mutation was present in all three replicates of EL2, suggesting that it was not random and probably associated with the adaptive evolution in acetic acid. The gene fhuA is involved in the transport of antibiotics such as albomycin and rifamycin, and bacteriophages and toxins such as colicin M inside E. coli cells (Braun et al., 2002). Mutations in this gene have been reported in E. coli and S. enterica under stress conditions, which enable the bacteria to escape from future attacks of antibiotics while maintaining its siderophore uptake activity that is essential for the survival of the bacteria (Wang et al., 2018). Additionally, two of the three replicates of EL2 showed mutations in the rpoB gene, which is usually targeted by antibiotics such as rifampicin in E. coli and S. enterica (Björkman et al., 1998). Several studies have shown that mutations in the rpoB gene lead to resistance against rifampicin and other similar antibiotics in E. coli and S. enterica (Björkman et al., 1998; Garibyan et al., 2003; Rodríguez-Verdugo et al., 2013). Mutations in rpoB are also known to be triggered during stress conditions in antibiotic-free environments (Rodríguez-Verdugo et al., 2013). This is similar to our study, where mutations in rpoB gene were triggered during adaptation to acid stress in the absence of antibiotics. Another gene that was mutated in all three replicates of EL3 and EL4 is rfbE, which is involved in the synthesis of the O-antigen in E. coli and the rbf gene cluster is known to be involved in the synthesis of the O-antigen of cell wall LPS in S. enterica and other gram-negative Enterobacteriaceae (Iguchi et al., 2011; Cota et al., 2015). The O-antigen is involved in the pathogenesis and virulence of S. enterica (Murray et al., 2003). Mutations in the O-antigen in the cell wall LPS are likely to affect the conformation and function of various surface proteins. The composition of the O-antigen affects bacterial surface charge, which is known to alter membrane permeability under changing pH conditions (Lerouge and Vanderleyden, 2002). We measured the surface charge of the ELs after 70 ALE days in acid stress and compared them to the surface charge of WT S. Enteritidis to study whether acid adaptation caused any changes in the surface charge. No significant differences in the surface charge of the ELs were observed as compared to the WT S. Enteritidis after 70 ALE days (p > 0.05) (Supplementary Table 2). Furthermore, the gene thiB that encodes for the thiamine-binding periplasmic protein was found to be mutated in the three replicates of EL3 and EL4. The presence of environmental stresses has been reported to cause mutations in genes involved in thiamine biosynthesis and transport in yeast cells (Kartal et al., 2018).

Finally, the deletion of 2207 bp of the genome in two replicates of EL3 and all three replicates of EL4 was observed. There were two major genes in this region, namely, ackA, which encodes for acetate kinase, and hxpA, which encodes for hexitol phosphatase. The gene ackA encodes for acetate kinase that is involved in the production of acetate from acetyl phosphate (Kumari et al., 2000). Deletion of this gene could be attributed to the presence of large amounts of acetate available during acetic acid stress, and hence, the need for this enzyme was significantly reduced in the ELs. Even though the ackA gene is also involved in the conversion of acetate to acetyl phosphate, which is a precursor to acetyl-CoA, there is an alternative pathway for the synthesis of acetyl-CoA from acetate (Kumari et al., 1995). This pathway is catalyzed by the enzyme acetyl-CoA synthase, and it is commonly used by cells during stressful conditions such as nutrient limitation, low oxygen pressure (De Mets et al., 2019), and, possibly, acid stress.

Since EL2 is the primary ancestor of the evolutionary lineages EL3 and EL4, we studied the genomes of EL2–EL4 to look for any common mutations that may have occurred before ALE day 20 and subsequently passed on to EL3 and EL4. We did not observe any common mutations in EL2–EL4 suggesting that most of the mutations observed in EL2 happened after ALE day 20 when EL3 diverged from EL2. Even though replicate 2 of EL2 and all three replicates of EL3 and EL4 had mutations in the phoQ gene, the positions of this mutation in EL3 and EL4 were different from EL2, indicating these occurred after ALE day 20. Missense mutations in rfbE, phoQ, and thiB were identical between EL3 and EL4, which suggests that these mutations occurred between ALE day 20 (after the split of EL3 from EL2) and ALE day 30 (before the initiation of EL4 from EL3). SNPs in certain genes such as tuf1 and oadB were observed in EL4 alone, indicating these mutations occurred after ALE day 30 when EL4 was split from EL3. EL3 and EL4 demonstrated a 2207 base pair mutation between genes ackA and hxpA. In EL3 and EL4, this mutation was observed in replicates 1 and 2, which possibly indicates that they were carried from EL3 to EL4.

Although our study is one of the first of its kind to study acid adaptation in S. enterica, prior studies have revealed the effects of acid adaptation in E. coli, which is a related gram-negative pathogen. Harden et al. (2015) documented that E. coli K12 lineages exposed to acid conditions exhibit enhanced growth rates and improved fitness. The same phenomenon was replicated in our investigation, where we observed an increase in growth rates among acid-evolved S. enterica lineages when subjected to acidic conditions. Other researchers have also independently corroborated these findings by demonstrating increased growth rates and enhanced survival in acid-evolved E. coli lineages exposed to acidic environments (He et al., 2017; Du et al., 2020). These researchers also reported that the acid-evolved E. coli lineages had mutations in the RNA polymerase subunits. Our study supports these findings as we also observed mutations in the RNA polymerase subunit beta in EL2. These findings suggest that RNA polymerases play an important role in acid stress adaptation in many gram-negative pathogens. Mutations in the RNA polymerase can lead to changes in gene expression, which could affect the production of acid-resistant proteins, the transport of protons out of the cell, or other cellular processes that are important for acid survival.

The results of this study highlight that adaptation to acetic acid stress is associated with genotypic and phenotypic changes in S. enterica. Similar genetic and phenotypic changes in Salmonella are expected if other organic acids with similar chain lengths and mechanisms of action against S. enterica were used during the adaptive process. The increased MIC of the ELs against both acetic acid and human antibiotics raises concerns regarding the development of cross-resistance and the potential emergence of multidrug-resistant S. enterica. Although the specific mechanisms associated with the cross-resistance remain to be elucidated, it is plausible that the mutations identified in the phoPQ system and other stress-related genes along with the gene deletions contribute to a broader adaptive response, leading to increased resistance against multiple stressors, including antibiotics.

5. Conclusion

Understanding the genetic mechanisms associated with acid stress adaptation can lead to targeted interventions and therapies to combat antimicrobial resistance. The uniqueness of our study lies in the fact that it is the first of its kind (to the best of our knowledge) to apply adaptive laboratory evolution for studying the adaptation of S. enterica under acid stress. The findings of this study have implications for food safety and demonstrate the utility of the ALE approach in studying the evolution and adaptation of foodborne pathogens to environmental stresses, particularly acid stress. Organic acids, especially acetic acid, are often used as additives for flavoring and preserving different foods and for cleaning and sanitization (Mani-López et al., 2012; Coban, 2020). It is possible that S. enterica, present as a contaminant in the food industry, is exposed to sub-lethal levels of acetic acid. The results of this study highlight that under such conditions, S. enterica will be able to adapt and evolve in acid stress. This might lead to enhanced persistence and survival of the bacteria, leading to a higher risk of contamination of food products during processing. The increase in MIC against human antibiotics in the acid-adapted S. enterica strains can also pose a significant public health risk. The results of this ALE study emphasize the need for continuous evaluation and optimization of sanitization protocols to ensure their efficacy against emerging antibiotic-resistant pathogens. Finally, this study underscores the need for comprehensive approaches to combat antimicrobial resistance, including the prudent use of acid-based antimicrobials and the development of alternative antimicrobial strategies.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Author contributions

MG: Conceptualization, Data curation, Investigation, Methodology, Writing – original draft, Writing – review and editing. TB: Formal analysis, Methodology, Writing – original draft. JG: Conceptualization, Software, Supervision, Visualization, Writing – review and editing. LM: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – review and editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was by a Rapid Research Grant through the University of Massachusetts, Amherst Graduate School (to MG), supported by the National Institute of Food and Agriculture (NIFA), the U.S. Department of Agriculture (USDA), and the Center for Agriculture, Food, and the Environment and the Department of Food Science at the University of Massachusetts-Amherst, under project number MAS00567, with support from the Foundational and Applied Science Program (grant number 2020-67017-30786) (to LM).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1285421/full#supplementary-material

References

Álvarez-Ordóñez, A., Prieto, M., Bernardo, A., Hill, C., and López, M. (2012). The acid tolerance response of Salmonella spp.: An adaptive strategy to survive in stressful environments prevailing in foods and the host. Food Res. Int. 45, 482–492.

Google Scholar

Amitai, S., Kolodkin-Gal, I., Hananya-Meltabashi, M., Sacher, A., and Engelberg-Kulka, H. (2009). Escherichia coli MazF leads to the simultaneous selective synthesis of both “death proteins” and “survival proteins”. PLoS Genet. 5:e1000390. doi: 10.1371/journal.pgen.1000390

PubMed Abstract | CrossRef Full Text | Google Scholar

Anton Bankevich, S., Dmitry Antipov, N., Alexey, A., Gurevich, Mikhail Dvorkin, Alexander, S., et al. (2012). SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477.

Google Scholar

Askoura, M., and Hegazy, W. A. H. (2020). Ciprofloxacin interferes with Salmonella Typhimurium intracellular survival and host virulence through repression of Salmonella pathogenicity island-2 (SPI-2) genes expression. Pathog. Dis. 78:ftaa011. doi: 10.1093/femspd/ftaa011

PubMed Abstract | CrossRef Full Text | Google Scholar

Athamneh, A. I. M., Alajlouni, R. A., Wallace, R. S., Seleem, M. N., and Senger, R. S. (2014). Phenotypic profiling of antibiotic response signatures in Escherichia coli using raman spectroscopy. Antimicrob. Agents Chemother. 58, 1302–1314. doi: 10.1128/AAC.02098-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Bearson, B. L., Wilson, L., and Foster, J. W. (1998). A low pH-inducible, PhoPQ-dependent acid tolerance response protects Salmonella typhimurium against inorganic acid stress. J. Bacteriol. 180, 2409–2417. doi: 10.1128/JB.180.9.2409-2417.1998

PubMed Abstract | CrossRef Full Text | Google Scholar

Björkman, J., Hughes, D., and Andersson, D. I. (1998). Virulence of antibiotic-resistant Salmonella typhimurium. Proc. Natl. Acad. Sci. U.S.A. 95, 3949–3953.

Google Scholar

Blair, J. M. A., La Ragione, R. M., Woodward, M. J., and Piddock, L. J. V. (2009). Periplasmic adaptor protein AcrA has a distinct role in the antibiotic resistance and virulence of Salmonella enterica serovar Typhimurium. J. Antimicrob. Chemother. 64, 965–972. doi: 10.1093/jac/dkp311

PubMed Abstract | CrossRef Full Text | Google Scholar

Brandl, M. T., Pan, Z., Huynh, S., Zhu, Y., and McHugh, T. H. (2008). Reduction of Salmonella enteritidis population sizes on almond kernels with infrared heat. J. Food Protect. 71, 897–902. doi: 10.4315/0362-028x-71.5.897

PubMed Abstract | CrossRef Full Text | Google Scholar

Braun, M., Killmann, H., Maier, E., Benz, R., and Braun, V. (2002). Diffusion through channel derivatives of the Escherichia coli FhuA transport protein. Eur. J. Biochem. 269, 4948–4959. doi: 10.1046/j.1432-1033.2002.03195.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Canali, C., Spillum, E., Valvik, M., Agersnap, N., and Olesen, T. (2018). Real-time digital bright field technology for rapid antibiotic susceptibility testing. Methods Mol. Biol. 1736, 75–84. doi: 10.1007/978-1-4939-7638-6_7

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, M.-X., Zhang, J.-F., Sun, Y.-H., Li, R.-S., Lin, X.-L., Yang, L., et al. (2021). Contribution of different mechanisms to ciprofloxacin resistance in Salmonella spp. Front. Microbiol. 12:663731. doi: 10.3389/fmicb.2021.663731663731

CrossRef Full Text | Google Scholar

Chapman, J. S. (2003). Disinfectant resistance mechanisms, cross-resistance, and co-resistance. Int. Biodeteriorat. Biodegradat. 51, 271–276. doi: 10.1093/jac/49.4.631

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, K.-K., Lee, B.-S., Masuda, T., Ito, T., Ikeda, K., Hirayama, A., et al. (2014). Global metabolic network reorganization by adaptive mutations allows fast growth of Escherichia coli on glycerol. Nat. Commun. 5:3233. doi: 10.1038/ncomms4233

PubMed Abstract | CrossRef Full Text | Google Scholar

Cherstvy, A. G. (2009). Positively charged residues in DNA-binding domains of structural proteins follow sequence-specific positions of DNA phosphate groups. J. Phys. Chem. B 113, 4242–4247. doi: 10.1021/jp810009s

PubMed Abstract | CrossRef Full Text | Google Scholar

Cingolani, P., Platts, A., Wang, L. L., Coon, M., Nguyen, T., Wang, L., et al. (2012). A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff. Fly 6, 80–92. doi: 10.4161/fly.19695

PubMed Abstract | CrossRef Full Text | Google Scholar

Coban, H. B. (2020). Organic acids as antimicrobial food agents: Applications and microbial productions. Bioprocess Biosyst. Eng. 43, 569–591.

Google Scholar

Cota, I., Sánchez-Romero, M. A., Hernández, S. B., Pucciarelli, M. G., García-del Portillo, F., and Casadesús, J. (2015). Epigenetic control of Salmonella enterica O-Antigen chain length: A tradeoff between virulence and bacteriophage resistance. PLoS Genet. 11:e1005667. doi: 10.1371/journal.pgen.1005667

PubMed Abstract | CrossRef Full Text | Google Scholar

Danecek, P., Bonfield, J. K., Liddle, J., Marshall, J., Ohan, V., Pollard, M. O., et al. (2021). Twelve years of SAMtools and BCFtools. GigaScience 10:giab008. doi: 10.1093/gigascience/giab008

PubMed Abstract | CrossRef Full Text | Google Scholar

Dawan, J., and Ahn, J. (2022). Bacterial stress responses as potential targets in overcoming antibiotic resistance. Microorganisms 10:1385.

Google Scholar

Dawan, J., Li, Y., Lu, F., He, X., and Ahn, J. (2022). Role of efflux pump-mediated antibiotic resistance in quorum sensing-regulated biofilm formation by Salmonella typhimurium. Pathogens 11:147. doi: 10.3390/pathogens11020147

PubMed Abstract | CrossRef Full Text | Google Scholar

De Mets, F., Van Melderen, L., and Gottesman, S. (2019). Regulation of acetate metabolism and coordination with the TCA cycle via a processed small RNA. Proc. Natl. Acad. Sci. U.S.A. 116, 1043–1052. doi: 10.1073/pnas.1815288116

PubMed Abstract | CrossRef Full Text | Google Scholar

Deatherage, D. E., and Barrick, J. E. (2014). Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using breseq. Methods Mol. Biol. 1151, 165–188. doi: 10.1007/978-1-4939-0554-6_12

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, X., Li, Z., and Zhang, W. (2012). Transcriptome sequencing of Salmonella enterica serovar Enteritidis under desiccation and starvation stress in peanut oil. Food Microbiol. 30, 311–315. doi: 10.1016/j.fm.2011.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Dragosits, M., and Mattanovich, D. (2013). Adaptive laboratory evolution – principles and applications for biotechnology. Microb. Cell Fact. 12:64.

Google Scholar

Du, B., Olson, C. A., Sastry, A. V., Fang, X., Phaneuf, P. V., Chen, K., et al. (2020). Adaptive laboratory evolution of Escherichia coli under acid stress. Microbiology 166, 141–148. doi: 10.1099/mic.0.000867

PubMed Abstract | CrossRef Full Text | Google Scholar

Duval, V., and Lister, I. M. (2013). MarA, SoxS and Rob of Escherichia coli - Global regulators of multidrug resistance, virulence and stress response. Int. J. Biotechnol. Wellness Ind. 2, 101–124. doi: 10.6000/1927-3037.2013.02.03.2

PubMed Abstract | CrossRef Full Text | Google Scholar

FDA (2023). Federal NARMS partners. Available online at: https://www.fda.gov/animal-veterinary/national-antimicrobial-resistance-monitoring-system/resources (accessed August 7, 2023).

Google Scholar

Fong, S., Joyce, A., and Palsson, B. (2005). Parallel adaptive evolution cultures of Escherichia coli lead to convergent growth phenotypes with different gene expression states. Genome Res. 15, 1365–1372. doi: 10.1101/gr.3832305

PubMed Abstract | CrossRef Full Text | Google Scholar

Fredborg, M., Andersen, K. R., Jørgensen, E., Droce, A., Olesen, T., Jensen, B. B., et al. (2013). Real-time optical antimicrobial susceptibility testing. J. Clin. Microbiol. 51, 2047–2053.

Google Scholar

Fredborg, M., Rosenvinge, F. S., Spillum, E., Kroghsbo, S., Wang, M., and Sondergaard, T. E. (2015). Rapid antimicrobial susceptibility testing of clinical isolates by digital time-lapse microscopy. Eur. J. Clin. Microbiol. Infect. Dis. 34, 2385–2394.

Google Scholar

Garibyan, L., Huang, T., Kim, M., Wolff, E., Nguyen, A., Nguyen, T., et al. (2003). Use of the rpoB gene to determine the specificity of base substitution mutations on the Escherichia coli chromosome. DNA Repair 2, 593–608. doi: 10.1016/s1568-7864(03)00024-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Garrison, E., and Marth, G. (2012). Haplotype-based variant detection from short-read sequencing. BioArxiv. [Preprint]. doi: 10.48550/arXiv.1207.3907

CrossRef Full Text | Google Scholar

Geng, J., Liu, H., Chen, S., Long, J., Jin, Y., Yang, H., et al. (2022). Comparative genomic analysis of Escherichia coli strains obtained from continuous imipenem stress evolution. FEMS Microbiol. Lett. 369:fnac015. doi: 10.1093/femsle/fnac015

PubMed Abstract | CrossRef Full Text | Google Scholar

Geornaras, I., and von Holy, A. (2001). Antimicrobial susceptibilities of isolates of Staphylococcus aureus, Listeria species and Salmonella serotypes associated with poultry processing. Int. J. Food Microbiol. 70, 29–35. doi: 10.1016/s0168-1605(01)00517-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghoshal, M., Ryu, V., and McLandsborough, L. (2022b). Evaluation of the efficacy of antimicrobials against pathogens on food contact surfaces using a rapid microbial log reduction detection method. Int. J. Food Microbiol. 373:109699. doi: 10.1016/j.ijfoodmicro.2022.109699

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghoshal, M., Chuang, S., Zhang, Y., and McLandsborough, L. (2022a). Efficacy of acidified oils against Salmonella in low-moisture environments. Appl. Environ. Microbiol. 88:e00935-22. doi: 10.1128/aem.00935-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Gresham, D., and Dunham, M. J. (2014). The enduring utility of continuous culturing in experimental evolution. Genomics 104(6, Part A), 399–405. doi: 10.1016/j.ygeno.2014.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Gullberg, E., Cao, S., Berg, O. G., Ilbäck, C., Sandegren, L., Hughes, D., et al. (2011). Selection of resistant bacteria at very low antibiotic concentrations. PLoS Pathog. 7:e1002158. doi: 10.1371/journal.ppat.1002158

PubMed Abstract | CrossRef Full Text | Google Scholar

Gurevich, A., Saveliev, V., Vyahhi, N., and Tesler, G. (2013). QUAST: Quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075.

Google Scholar

Haft, D. H., DiCuccio, M., Badretdin, A., Brover, V., Chetvernin, V., O’Neill, K., et al. (2018). RefSeq: An update on prokaryotic genome annotation and curation. Nucleic Acids Res. 46, D851–D860. doi: 10.1093/nar/gkx1068

PubMed Abstract | CrossRef Full Text | Google Scholar

Hakvoort, H., Bovenkamp, E., Greenwood-Quaintance, K. E., Schmidt-Malan, S. M., Mandrekar, J. N., Schuetz, A. N., et al. (2020). Imipenem-relebactam susceptibility testing of gram-negative bacilli by agar dilution, disk diffusion, and gradient strip methods compared with broth Microdilution. J. Clin. Microbiol. 58:e00695-20. doi: 10.1128/JCM.00695-20

PubMed Abstract | CrossRef Full Text | Google Scholar

Hall, B. G., Acar, H., Nandipati, A., and Barlow, M. (2013). Growth rates made easy. Mol. Biol. Evol. 31, 232–238.

Google Scholar

Harden, M. M., He, A., Creamer, K., Clark, M. W., Hamdallah, I., and Martinez, K. A. II, et al. (2015). Acid-adapted strains of Escherichia coli K-12 obtained by experimental evolution. Appl. Environ. Microbiol. 81, 1932–1941. doi: 10.1128/AEM.03494-14

PubMed Abstract | CrossRef Full Text | Google Scholar

He, A., Penix, S. R., Basting, P. J., Griffith, J. M., Creamer, K. E., Camperchioli, D., et al. (2017). Acid evolution of Escherichia coli K-12 eliminates amino acid decarboxylases and reregulates catabolism. Appl. Environ. Microbiol. 83:e00442-17. doi: 10.1128/AEM.00442-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Herring, C. D., Raghunathan, A., Honisch, C., Patel, T., Applebee, M. K., Joyce, A. R., et al. (2006). Comparative genome sequencing of Escherichia coli allows observation of bacterial evolution on a laboratory timescale. Nat. Genet. 38, 1406–1412. doi: 10.1038/ng1906

PubMed Abstract | CrossRef Full Text | Google Scholar

Higginson, E. E., Ramachandran, G., Hazen, T. H., Kania, D. A., Rasko, D. A., Pasetti, M. F., et al. (2018). Improving our understanding of Salmonella enterica serovar paratyphi B through the engineering and testing of a live attenuated vaccine strain. mSphere 3:e00474-00418. doi: 10.1128/msphere.00474-00418

CrossRef Full Text | Google Scholar

Hu, M., Zhang, Y., Huang, X., He, M., Zhu, J., Zhang, Z., et al. (2023). PhoPQ regulates quinolone and cephalosporin resistance formation in Salmonella enteritidis at the transcriptional level. mBio 14:e0339522. doi: 10.1128/mbio.03395-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Iguchi, A., Shirai, H., Seto, K., Ooka, T., Ogura, Y., Hayashi, T., et al. (2011). Wide distribution of O157-antigen biosynthesis gene clusters in Escherichia coli. PLoS One 6:e23250. doi: 10.1371/journal.pone.0023250

PubMed Abstract | CrossRef Full Text | Google Scholar

Jean, S.-S., Lee, Y.-T., Guo, S.-M., and Hsueh, P.-R. (2005). Recurrent infections caused by cefotaxime- and ciprofloxacin-resistant Salmonella enterica serotype choleraesuis treated successfully with imipenem. J. Infect. 51, e163–e165. doi: 10.1016/j.jinf.2004.12.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Jibril, A. H., Okeke, I. N., Dalsgaard, A., Menéndez, V. G., and Olsen, J. E. (2021). Genomic analysis of antimicrobial resistance and resistance plasmids in Salmonella Serovars from poultry in nigeria. Antibiotics 10:99. doi: 10.3390/antibiotics10020099

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, M., Zaretskaya, I., Raytselis, Y., Merezhuk, Y., McGinnis, S., and Madden, T. L. (2008). NCBI BLAST: A better web interface. Nucleic Acids Res. 36(Suppl._2), W5–W9. doi: 10.1093/nar/gkn201

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, R. N., Rhomberg, P. R., Varnam, D. J., and Mathai, D. (2002). A comparison of the antimicrobial activity of meropenem and selected broad-spectrum antimicrobials tested against multi-drug resistant gram-negative bacilli including bacteraemic Salmonella spp.: Initial studies for the MYSTIC programme in India. Int. J. Antimicrob. Agents 20, 426–431. doi: 10.1016/s0924-8579(02)00210-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Kartal, B., Akçay, A., and Palabiyik, B. (2018). Oxidative stress upregulates the transcription of genes involved in thiamine metabolism. Turk J. Biol. 42, 447–452.

Google Scholar

Koeth, L. K., DiFranco-Fisher, J. M., Hardy, D. J., Palavecino, E. L., Carretto, E., and Windau, A. (2022). Multilaboratory comparison of Omadacycline MIC test strip to broth microdilution MIC against gram-negative, gram-positive, and fastidious bacteria. J. Clin. Microbiol. 60:e0141021. doi: 10.1128/JCM.01410-21

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumari, S., Beatty, C. M., Browning, D. F., Busby, S. J. W., Simel, E. J., Hovel-Miner, G., et al. (2000). Regulation of acetyl coenzyme A Synthetase in Escherichia coli. J. Bacteriol. 182, 4173–4179. doi: 10.1128/JB.182.15.4173-4179.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumari, S., Tishel, R., Eisenbach, M., and Wolfe, A. J. (1995). Cloning, characterization, and functional expression of acs, the gene which encodes acetyl coenzyme A synthetase in Escherichia coli. J. Bacteriol. 177, 2878–2886. doi: 10.1128/jb.177.10.2878-2886.1995

PubMed Abstract | CrossRef Full Text | Google Scholar

LaCroix, R. A., Palsson, B. O., and Feist, A. M. (2017). A model for designing adaptive laboratory evolution experiments. Appl. Environ. Microbiol. 83:e03115-16. doi: 10.1128/AEM.03115-16

PubMed Abstract | CrossRef Full Text | Google Scholar

LaCroix, R. A., Sandberg, T. E., O’Brien, E. J., Utrilla, J., Ebrahim, A., Guzman, G. I., et al. (2015). Use of adaptive laboratory evolution to discover key mutations enabling rapid growth of Escherichia coli K-12 MG1655 on glucose minimal medium. Appl. Environ. Microbiol. 81, 17–30. doi: 10.1128/AEM.02246-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Langsrud, S., Sundheim, G., and Holck, A. L. (2004). Cross-resistance to antibiotics of Escherichia coli adapted to benzalkonium chloride or exposed to stress-inducers. J. Appl. Microbiol. 96, 201–208.

Google Scholar

Lázár, V., Nagy, I., Spohn, R., Csörgõ, B., Györkei, Á, Nyerges, Á, et al. (2014). Genome-wide analysis captures the determinants of the antibiotic cross-resistance interaction network. Nat. Commun. 5:4352. doi: 10.1038/ncomms5352

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, D. H., Feist, A. M., Barrett, C. L., and Palsson, B. (2011). Cumulative number of cell divisions as a meaningful timescale for adaptive laboratory evolution of Escherichia coli. PLoS One 6:e26172. doi: 10.1371/journal.pone.0026172

PubMed Abstract | CrossRef Full Text | Google Scholar

Lemmin, T., Soto, C. S., Clinthorne, G., DeGrado, W. F., and Dal Peraro, M. (2013). Assembly of the Transmembrane domain of E. coli PhoQ Histidine Kinase: Implications for signal transduction from molecular simulations. PLoS Comput. Biol. 9:e1002878. doi: 10.1371/journal.pcbi.1002878

PubMed Abstract | CrossRef Full Text | Google Scholar

Lerouge, I., and Vanderleyden, J. (2002). O-antigen structural variation: Mechanisms and possible roles in animal/plant–microbe interactions. FEMS Microbiol. Rev. 26, 17–47. doi: 10.1111/j.1574-6976.2002.tb00597.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Leyer, G. J., and Johnson, E. A. (1992). Acid adaptation promotes survival of Salmonella spp. in cheese. Appl Environ. Microbiol. 58, 2075–2080. doi: 10.1128/aem.58.6.2075-2080.1992

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H. (2013). Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. BioArxiv. [Preprint]. doi: 10.48550/arXiv.1303.3997.

CrossRef Full Text | Google Scholar

Li, Y., and Zhang, Y. (2007). PhoU is a persistence switch involved in persister formation and tolerance to multiple antibiotics and stresses in Escherichia coli. Antimicrob. Agents Chemother. 51, 2092–2099. doi: 10.1128/AAC.00052-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Maeda, T., Iwasawa, J., Kotani, H., Sakata, N., Kawada, M., Horinouchi, T., et al. (2020). High-throughput laboratory evolution reveals evolutionary constraints in Escherichia coli. Nat. Commun. 11:5970. doi: 10.1038/s41467-020-19713-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Maldonado, R. F., Sá-Correia, I., and Valvano, M. A. (2016). Lipopolysaccharide modification in Gram-negative bacteria during chronic infection. FEMS Microbiol. Rev. 40, 480–493.

Google Scholar

Mani-López, E., García, H. S., and ópez-Malo, A. L. (2012). Organic acids as antimicrobials to control Salmonella in meat and poultry products. Food Res. Int. 45, 713–721.

Google Scholar

Marciano, D. C., Wang, C., Hsu, T.-K., Bourquard, T., Atri, B., Nehring, R. B., et al. (2022). Evolutionary action of mutations reveals antimicrobial resistance genes in Escherichia coli. Nat. Commun. 13:3189. doi: 10.1038/s41467-022-30889-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Matamouros, S., Hager, K., and Miller, S. (2015). HAMP domain rotation and tilting movements associated with signal transduction in the PhoQ Sensor Kinase. mBio 6:e00616-15. doi: 10.1128/mBio.00616-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Matuschek, E., Åhman, J., Webster, C., and Kahlmeter, G. (2018). Antimicrobial susceptibility testing of colistin – evaluation of seven commercial MIC products against standard broth microdilution for Escherichia coli, Klebsiella pneumoniae, Pseudomonas aeruginosa, and Acinetobacter spp. Clin. Microbiol. Infect. 24, 865–870. doi: 10.1016/j.cmi.2017.11.020

PubMed Abstract | CrossRef Full Text | Google Scholar

McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The genome analysis toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303. doi: 10.1101/gr.107524.110

PubMed Abstract | CrossRef Full Text | Google Scholar

McLaughlin, H. P., Gargis, A. S., Michel, P., Sue, D., and Weigel, L. M. (2017). Optical screening for rapid antimicrobial susceptibility testing and for observation of phenotypic diversity among strains of the genetically clonal species Bacillus anthracis. J. Clin. Microbiol. 55, 959–970. doi: 10.1128/JCM.02209-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Mengistu, G., Dejenu, G., Tesema, C., Arega, B., Awoke, T., Alemu, K., et al. (2020). Epidemiology of streptomycin resistant Salmonella from humans and animals in Ethiopia: A systematic review and meta-analysis. PLoS One 15:e0244057. doi: 10.1371/journal.pone.0244057

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, A. K., Brannon, M. K., Stevens, L., Johansen, H. K., Selgrade, S. E., Miller, S. I., et al. (2011). PhoQ mutations promote lipid A modification and polymyxin resistance of Pseudomonas aeruginosa found in colistin-treated cystic fibrosis patients. Antimicrob. Agents Chemother. 55, 5761–5769. doi: 10.1128/AAC.05391-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyakoshi, M., Matera, G., Maki, K., Sone, Y., and Vogel, J. (2019). Functional expansion of a TCA cycle operon mRNA by a 3’ end-derived small RNA. Nucleic Acids Res. 47, 2075–2088. doi: 10.1093/nar/gky1243

PubMed Abstract | CrossRef Full Text | Google Scholar

Mongold, J. A., Bennett, A. F., and Lenski, R. E. (1999). Evolutionary adaptation to temperature: Extension of the upper thermal limit of Escherichia coli. Evolution 53, 386–394. doi: 10.1111/j.1558-5646.1999.tb03774.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Muhammad, M., Muhammad, L. U., Ambali, A.-G., Mani, A. U., Azard, S., and Barco, L. (2010). Prevalence of Salmonella associated with chick mortality at hatching and their susceptibility to antimicrobial agents. Vet. Microbiol. 140, 131–135. doi: 10.1016/j.vetmic.2009.07.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Murray, G. L., Attridge, S. R., and Morona, R. (2003). Regulation of Salmonella typhimurium lipopolysaccharide O antigen chain length is required for virulence; identification of FepE as a second Wzz. Mol. Microbiol. 47, 1395–1406. doi: 10.1046/j.1365-2958.2003.03383.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pacios, O., Blasco, L., Bleriot, I., Fernandez-Garcia, L., Ambroa, A., López, M., et al. (2020). (p)ppGpp and Its role in bacterial persistence: New challenges. Antimicrob. Agents Chemother. 64:e01283-20.

Google Scholar

Pitt, M. E., Elliott, A. G., Cao, M. D., Ganesamoorthy, D., Karaiskos, I., Giamarellou, H., et al. (2018). Multifactorial chromosomal variants regulate polymyxin resistance in extensively drug-resistant Klebsiella pneumoniae. Microb. Genom. 4:e000158. doi: 10.1099/mgen.0.000158

PubMed Abstract | CrossRef Full Text | Google Scholar

Prost, L. R., Daley, M. E., Le Sage, V., Bader, M. W., Le Moual, H., Klevit, R. E., et al. (2007). Activation of the bacterial sensor kinase PhoQ by Acidic pH. Mol. Cell 26, 165–174.

Google Scholar

Ren, J., Sang, Y., Ni, J., Tao, J., Lu, J., Zhao, M., et al. (2015). Acetylation regulates survival of Salmonella enterica Serovar typhimurium under acid stress. Appl. Environ. Microbiol. 81, 5675–5682. doi: 10.1128/AEM.01009-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Rhee, J. E., Sheng, W., Morgan, L. K., Nolet, R., Liao, X., and Kenney, L. J. (2008). Amino acids important for DNA recognition by the response regulator OmpR. J. Biol. Chem. 283, 8664–8677. doi: 10.1074/jbc.M705550200

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodríguez-Verdugo, A., Gaut, B. S., and Tenaillon, O. (2013). Evolution of Escherichia coli rifampicin resistance in an antibiotic-free environment during thermal stress. BMC Evol. Biol. 13:50. doi: 10.1186/1471-2148-13-50

PubMed Abstract | CrossRef Full Text | Google Scholar

Sahl, J. W., Caporaso, J. G., Rasko, D. A., and Keim, P. (2014). The large-scale blast score ratio (LS-BSR) pipeline: A method to rapidly compare genetic content between bacterial genomes. PeerJ 2:e332. doi: 10.7717/peerj.332

PubMed Abstract | CrossRef Full Text | Google Scholar

Sandberg, T. E., Salazar, M. J., Weng, L. L., Palsson, B. O., and Feist, A. M. (2019). The emergence of adaptive laboratory evolution as an efficient tool for biological discovery and industrial biotechnology. Metab. Eng. 56, 1–16. doi: 10.1016/j.ymben.2019.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Shprung, T., Wani, N. A., Wilmes, M., Mangoni, M. L., Bitler, A., Shimoni, E., et al. (2021). Opposing effects of PhoPQ and PmrAB on the properties of Salmonella enterica serovar typhimurium: Implications on resistance to antimicrobial peptides. Biochemistry 60, 2943–2955. doi: 10.1021/acs.biochem.1c00287

PubMed Abstract | CrossRef Full Text | Google Scholar

Spector, M. P., and Kenyon, W. J. (2012). Resistance and survival strategies of Salmonella enterica to environmental stresses. Food Res. Int. 45, 455–481.

Google Scholar

Tatusova, T., DiCuccio, M., Badretdin, A., Chetvernin, V., Nawrocki, E. P., Zaslavsky, L., et al. (2016). NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 44, 6614–6624.

Google Scholar

van den Bijllaardt, W., Schijffelen, M. J., Bosboom, R. W., Cohen Stuart, J., Diederen, B., and Kampinga, G. (2018). Susceptibility of ESBL Escherichia coli and Klebsiella pneumoniae to fosfomycin in the Netherlands and comparison of several testing methods including Etest, MIC test strip, Vitek2, Phoenix and disc diffusion. J. Antimicrob. Chemother. 73, 2380–2387. doi: 10.1093/jac/dky214

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Sheng, H., Xu, W., Huang, J., Meng, L., Cao, C., et al. (2019). Diversity of serotype, genotype, and antibiotic susceptibility of Salmonella prevalent in pickled ready-to-eat meat. Front. Microbiol. 10:2577. doi: 10.3389/fmicb.2019.02577

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Chen, X., Hu, Y., Zhu, G., White, A. P., and Köster, W. (2018). Evolution and sequence diversity of FhuA in Salmonella and Escherichia. Infect. Immun. 86:e00573-18. doi: 10.1128/IAI.00573-18

PubMed Abstract | CrossRef Full Text | Google Scholar

Wiegand, I., Hilpert, K., and Hancock, R. E. W. (2008). Agar and broth dilution methods to determine the minimal inhibitory concentration (MIC) of antimicrobial substances. Nat. Protoc. 3, 163–175.

Google Scholar

Wistrand-Yuen, E., Knopp, M., Hjort, K., Koskiniemi, S., Berg, O. G., and Andersson, D. I. (2018). Evolution of high-level resistance during low-level antibiotic exposure. Nat. Commun. 9:1599.

Google Scholar

Xu, H., Lee, H. Y., and Ahn, J. (2008). Cross-protective effect of acid-adapted Salmonella enterica on resistance to lethal acid and cold stress conditions. Lett. Appl. Microbiol. 47, 290–297. doi: 10.1111/j.1472-765x.2008.02429.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yarlagadda, V., Sarkar, P., Samaddar, S., Manjunath, G. B., Mitra, S. D., Paramanandham, K., et al. (2018). Vancomycin analogue restores Meropenem activity against NDM-1 gram-negative pathogens. ACS Infect. Dis. 4, 1093–1101. doi: 10.1021/acsinfecdis.8b00011

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, J., Jin, F., Glatter, T., and Sourjik, V. (2017). Osmosensing by the bacterial PhoQ/PhoP two-component system. Proc. Natl. Acad. Sci. U.S.A. 114, E10792–E10798. doi: 10.1073/pnas.1717272114

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, J., Wu, L., Liu, Z., Lv, Y., Feng, J., Wang, W., et al. (2021). Gain-of-function mutations in acid stress response (evgS) protect Escherichia coli from killing by gallium nitrate, an antimicrobial Candidate. Antimicrob. Agents Chemother. 65:e01595-20. doi: 10.1128/AAC.01595-20

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Salmonella enterica, adaptive laboratory evolution (ALE), acetic acid stress adaptation, minimum inhibitory concentration (MIC), bacterial whole genome sequencing, missense mutations

Citation: Ghoshal M, Bechtel TD, Gibbons JG and McLandsborough L (2023) Adaptive laboratory evolution of Salmonella enterica in acid stress. Front. Microbiol. 14:1285421. doi: 10.3389/fmicb.2023.1285421

Received: 29 August 2023; Accepted: 16 October 2023;
Published: 16 November 2023.

Edited by:

Sébastien Holbert, INRA Centre Val de Loire, France

Reviewed by:

Arnaud Bridier, Agence Nationale de Sécurité Sanitaire de l’Alimentation, de l’Environnement et du Travail (ANSES), France
Ahmed Ghamry Abdelhamid, The Ohio State University, United States

Copyright © 2023 Ghoshal, Bechtel, Gibbons and McLandsborough. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Lynne McLandsborough, bG1AZm9vZHNjaS51bWFzcy5lZHU=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.