Preprint
Article

The Effect of the Stringent Response and Oxidative Stress Response on Fitness Costs of De Novo Acquisition of Antibiotic Resistance

Altmetrics

Downloads

122

Views

77

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

09 November 2023

Posted:

09 November 2023

You are already at the latest version

Alerts
Abstract
Resistance evolution during exposure to non-lethal levels of antibiotics is influenced by various stress responses of bacteria which are known to affect growth rate. Here, we aim to disentangle how the interplay between resistance development and associated fitness costs is affected by stress responses. We performed de novo resistance evolution of wild-type strains and single-gene knockout strains in stress response pathways using four different antibiotics. Throughout resistance development, the increase in minimum inhibitory concentration (MIC) is accompanied by a gradual decrease in growth rate, most pronounced with amoxicillin or kanamycin. By measuring biomass yield on glucose and whole-genome sequences at intermediate and final timepoints, we identified two patterns of how the stress responses affect the correlation between MIC and growth rate. First, single-gene knockout E. coli strains associated with reactive oxygen species (ROS) acquire resistance faster and mutations related to antibiotic permeability and pumping out occur earlier. This increases the metabolic burden of resistant bacteria. Second, the ΔrelA knockout strain which has reduced (p)ppGpp synthesis, is restricted in its stringent response, leading to diminished growth rates. The ROS-related mutagenesis and the stringent response are vulnerabilities within the fitness of resistant strains that can possibly be targeted to prevent development of resistance.
Keywords: 
Subject: Biology and Life Sciences  -   Biochemistry and Molecular Biology

1. Introduction

Antibiotic resistance has become one of the major public health threats of the 21st century [1]. Overuse or misuse of antibiotics in hospitals, communities, and farms has led to the emergence, spread, and persistence of resistant strains, reducing the effectiveness of prevention and treatment of bacterial infections [2,3]. Understanding the mechanisms of de novo acquisition of antibiotic resistance is critical for the discovery of novel drug targets and the development of new antibiotics that will remain effective during their application. In addition, studying the physiological and metabolic costs of evolution of resistance may identify manners to reduce resistance by designing new therapies targeting physiological weaknesses associated with specific resistance mechanisms [4,5].
The development of antimicrobial resistance, both by de novo acquisition and horizontal transfer of resistance genes, is frequently accompanied by a decline in bacterial fitness [6,7,8]. Conversion of cellular homeostasis systems, such as ion- and pH-pumps, into antibiotic efflux pumps, causes a reduced capability to control the intracellular environment [9]. Mutations within genes associated with resistance, such as modification of drug target, degradation of antibiotics, reduced antibiotic uptake, and increased efflux, can result in diminished survival and reproductive capabilities among resistant populations [4]. Consequently, it has been suggested that limiting antibiotic usage, thus reducing the selection pressure favoring resistant populations, may allow more adaptable susceptible populations to outcompete and, ultimately, eradicate resistance [10]. However, several other responses of bacteria becoming resistant, such as co-selection of resistance genes alongside other functional traits can provide an adaptive advantage [11]. The acquisition of resistance plasmids did not reduce growth rates in E. coli [12]. Low-cost resistance mutations and compensatory evolution can compensate for the initial metabolic costs [10]. Therefore, a comprehensive understanding of the intricate relationship between mutation occurrence and fitness costs is essential for addressing the persistent challenge of antibiotic resistance.
Under antibiotic stress, bacteria initiate a variety of stress responses such as stringent and oxidative stress responses [13]. Bacterial stringent and oxidative stress responses are involved in the development of antibiotic resistance [14,15]. Bactericidal antibiotics interact with their targets in a way that increases the oxidation of NADH and thus produces respiratory chain byproducts such as reactive oxygen species (ROS) [16]. Although ROS at high concentrations kills bacteria, DNA damage caused by sub-lethal levels of ROS can activate the repair system induced by the SOS response, thereby increasing the mutation rate [14,17]. Ultimately, exposure to antibiotics can cause mutations in antimicrobial resistance-related genes. The single-gene knockout strain of (p)ppGpp synthase RelA has a reduced stringent response, resulting in lower growth rates and reduced ROS formation, reducing the development of drug resistance [15]. As these two stress responses also affect the growth rate, it is unclear to which extent the fitness costs during the development of resistance are related to the stress responses.
Important genes in the (p)ppGpp-related stringent response include hipA, hipB, relA, and rpoS. HipA and HipB are part of the type II toxin-antitoxin system and can affect glutamate-tRNA ligase [18]. HipA phosphorylates glutamate-tRNA ligase, causing uncharged tRNA (Glu) accumulation, thereby hindering the translation process [19]. The reduced translation rates increase (p)ppGpp levels, which, catalyzed by RelA/SpoT, activates the stringent response. As a result, RNA polymerase switches from transcribing growth and reproduction-related genes to stress-response-related genes, and (p)ppGpp accelerating amino acid biosynthesis [20]. HipB acts as an antitoxin component and can neutralize the toxicity of HipA [21]. Therefore, knocking out hipA or hipB may affect the stringent response induced by (p)ppGpp. RelA is a GTP pyrophosphokinase, that catalyzes the formation of (p)ppGpp, and a knockout of relA most likely directly affects the synthesis of (p)ppGpp [22]. (p)ppGpp binds to RNA polymerase and responds to stress through the sigma factor σS [23]. Therefore, directly knocking out the gene rpoS, encoding RNA polymerase sigma factor σS, has a potential impact on stress response. In the system dealing with ROS, sodA, sodB, soxR, and soxS play important roles. Superoxide dismutase SodA and SodB can destroy superoxide anion radicals [24]. Therefore, after knocking out these genes, the function of clearing excessively generated ROS in cells will be weakened to a certain extent. The gene soxR encodes a redox-sensitive transcriptional activator responsible for triggering the transcription of the superoxide response regulator SoxS, which plays a role in the removal of superoxide [25]. Consequently, the deletion of soxR or soxS can reduce the clearance of ROS.
In this study, we evaluate the correlation between changes in minimum inhibitory concentration (MIC) and growth rate as an indicator for fitness throughout the resistance acquisition process. For this assessment four single-gene knockout E. coli strains related to the stringent response were used and four related to the oxidative stress responses. The fitness costs associated with resistance was investigated for each strain by quantifying biomass yield on glucose. Evidence for the relationships between mutations and fitness was ultimately discovered through whole-genome sequencing. The overall picture that emerges is that initial mutations causing resistance come at a metabolic cost, that can subsequently be compensated by compensatory mutations. The complex cellular response to stress due to exposure to antibiotics, mediated by stringent response, oxidative stress response, and other cellular processes, allows the cell to overcome non-lethal concentrations of antimicrobials at minimal metabolic costs.

2. Results

2.1. Growth Rates Gradually Decreased during the Evolution of Resistance.

To investigate the effect of the stringent response and oxidative stress response on the fitness cost of resistant strains during de novo acquisition of antibiotic resistance, growth rates were measured at each MIC value reached. Four (p)ppGpp-related (ΔhipA, ΔhipB, ΔrelA, and ΔrpoS), four ROS-related (ΔsodA, ΔsodB, ΔsoxR, and ΔsoxS) single-gene knockout E. coli strains and the wild-type MG1655 strain were exposed to stepwise increasing levels of amoxicillin, enrofloxacin, kanamycin, and tetracycline. When these strains evolved higher resistance, the growth rate in antibiotic-free medium and the MIC were recorded, and it showed that generally the growth rate decreased with increasing MIC (Figure 1). Before the resistance evolution, the growth rates of all ancestor strains were tested, and no significant differences were found (Figure 1Q).
The growth rates during the amoxicillin resistance evolution of ΔhipA and ΔhipB were consistently similar to that of the wild type (WT) (Figure 1A). At the early stage of amoxicillin resistance evolution, when the MIC was 8 µg/mL, the growth rates of ΔrelA and ΔrpoS were clearly lower than WT, but when the MIC reached 2048 µg/mL, the growth rate of ΔrpoS was equal to that of WT (Figure 1E). For ROS-related knockout strains, the growth rate began to decrease at the middle stage of evolution when the MIC was 128 µg/mL, but it was similar to WT at the highest concentration (Figure 1I and M). When the MIC reached 128 µg/mL, the growth rates of ΔrelA, ΔsodB, ΔsoxR, and ΔsoxS were significantly reduced compared with WT (Figure 1R). The growth rate of all amoxicillin-resistant strains decreased gradually with the increase of MIC. Among them, ΔsoxR and ΔsoxS decreased the fastest and their slopes were 0.025 and 0.030, respectively, compared to 0.018 for the WT (Figure 1M).
The evolution process of enrofloxacin resistance was relatively complicated. Some strains such as ΔhipA, ΔrelA, and ΔrpoS cannot evolve to higher enrofloxacin concentrations, which makes them difficult to compare (Figure 1B and F). The other knockout strains had no significant difference compared with WT (Figure 1B, J, and N). However, enrofloxacin-resistant strains showed a lower slope of decline in growth rate than amoxicillin throughout the evolution.
During kanamycin resistance evolution, ΔhipA and ΔhipB had a higher growth rate than WT until reaching the middle phase, i.e., until the MIC reached 256 µg/mL (Figure 1C). Similarly, the growth rates of ΔsodA, ΔsodB, ΔsoxR, and ΔsoxS were slightly higher than those of WT at MIC’s lower than 128 µg/mL (Figure 1 K and O). However, their growth rates decreased significantly thereafter and reached the same level as WT at the highest concentration. Similar to the case of amoxicillin, the growth rate of ΔrelA during the kanamycin-resistance evolution was lower than that of the WT starting from MIC 16 µg/mL (Figure 1G). ΔrpoS showed a clear and gradual decreasing trend with an R-squared value of 0.957, indicating a high correlation between a decrease in growth rate and an increase in MIC (Figure 1G). Compared with WT, the growth rates of ΔrelA, ΔrpoS, ΔsodA, ΔsodB, ΔsoxR, and ΔsoxS were significantly lower than that of WT when the MIC was 256 µg/mL (Figure 1S). The growth rate reduction slopes of kanamycin-resistant evolution were steeper than those of amoxicillin, respectively, indicating that it has a higher fitness burden.
The growth rate during the evolution of tetracycline resistance was not significantly different from that of WT and either (p)ppGpp-associated or ROS-associated mutants (Figure 1D, H, L, and P). The highest MIC evolved by all the E. coli strains to the bacteriostatic antibiotic tetracycline was 128 µg/mL, which was lower than 2048 µg/mL or 4096 µg/mL to the bactericidal antibiotics. For all strains it is difficult to evolve very high resistance to tetracycline, limiting the comparison of growth rates between the strains.

2.2. Decrease in Biomass Yield on Glucose with Increased Resistance.

A clear correlation between the growth rate and MIC was observed during the amoxicillin or kanamycin resistance evolution, and both treatments growth rate showed a significant drop-down at the middle stage of resistance evolution in specific mutants compared to WT (Figure 1). To further understand the impact of stringent response and oxidative stress response on fitness cost during antimicrobial resistance evolution, we measured biomass yield on glucose during amoxicillin or kanamycin resistance evolution in three stages, initial, middle, and end stages (Figure 2).
Firstly, there was no significant difference between the ancestor WT MG1655 and the single-gene knockout mutants before resistance evolution (Figure 2). At the middle stage (MIC = 128 µg/mL) of amoxicillin resistance evolution, the biomass yield of ΔrelA, ΔsodB, ΔsoxR and ΔsoxS were significantly lower than that of the WT MG1655 (Figure 2A). Additionally, compared with the final resistant point of WT MG1655, the biomass yield of ΔrelA was significantly decreased (Figure 2A).
Except ΔhipA and ΔhipB in the middle stage of kanamycin resistant evolution, all the other mutants showed a lower biomass yield than that of the WT (Figure 2B). Similar to the amoxicillin resistance evolution, only the biomass yield of ΔrelA was significantly decreased compared to WT at the kanamycin evolution end point. The biomass yield of all final resistant strains was lower than their ancestor strains and middle stage of evolution resistant strains, respectively (Figure 2A and B). In summary, during the evolution of amoxicillin or kanamycin resistance, with the increase in resistance levels, the growth rate of specific (p)ppGpp- or ROS-related resistant strains gradually decreased, accompanied by decreased biomass yield on glucose.

2.3. Mutations in the Amoxicillin-Resistant Strains during the Resistance Evolution.

To detect the mutations and determine the correlation between fitness cost and mutations upon the development of antibiotic resistance, the whole genome of all amoxicillin or kanamycin resistant strains in the ancestor, middle, and end resistance points were sequenced. The shared mutations between the resistant strains at middle or end evolution points and their respective ancestor strains were excluded (Figure 3 and Figure 4).
All amoxicillin-resistant strains contained mutations in the ampC promoter region (Figure 3; Table 1). Most of these mutations appear in the middle stage and are preserved until the final stage. Some appear in the middle stage and disappear in the final stage, and some only appear in the end stage. The difference in these mutations may be a result of selection for different RNA polymerase recruitment abilities to retain the most efficient increase in transcription levels, ultimately increasing the expression of beta-lactamase AmpC [26]. Mutations occurring in the ampC promoter region are common among all strains, so their effects on fitness costs cannot be elucidated from this data.
Several other mutations emerged in different strains during the evolution of amoxicillin resistance (Figure 3). Resistance mechanisms involved in these mutations include reduced antibiotic permeability, increased antibiotic pumping, and altered antibiotic targets. Notably, in the middle of evolution, a cpxA mutation emerged in ΔhipA, ΔsodA, and ΔsoxR. In addition, other mutations in cpxA, but in different positions also occurred in a high population frequency at the ΔsodB and ΔsoxS middle point. Sensor histidine kinase CpxA initiates a kinase cascade leading to the activation of CpxR, which subsequently functions to enhance the expression of efflux complexes [27]. At the middle point of ΔsodB, other efflux pump-associated genes also mutated. Some other notable mutations at the middle point were envZ (ompB) in ΔsodB and ΔsoxS, and ompC in ΔsoxR. These mutations may be related to reduced antibiotic permeability [28,29]. To sum up, mutated genes related to antibiotic permeability and efflux pumps appeared simultaneously at high frequency in ΔsodB, ΔsoxR, and ΔsoxS. Compared to the wild type, these three strains had lower biomass yield at the middle stage of evolution (Figure 2A). This suggests that bacteria have evolved mutations that reduce antibiotic permeability and increase efflux, thereby increasing their energy expenditure burden.
At the end point of the evolution experiments, almost all the amoxicillin-resistant strains evolved mutations associated with reduced antibiotic permeability except ΔrelA (Figure 3). Surprisingly, only ΔrelA had significantly lower biomass yield on glucose at the endpoint than WT, while there were no obvious differences between the other strains (Figure 2A). Consistent with the biomass yield results observed at the midpoint, we hypothesize that in ΔrelA, the effect of glucose consumption due to the mutation is not as large as the effect of relA knockout itself. This phenomenon may be caused by the production of (p)ppGpp is inhibited to a certain extent after the deletion of relA. Along with the mutation's side effects, stringent responses are hindered. The consequence is that internal signals are in a configuration as if the cells are in a state of constant starvation. [30,31].

2.4. Mutations in the Kanamycin-Resistant Strains during the Resistance Evolution.

The common mutated gene of all kanamycin-resistant strains observed at the middle stage was fusA, which codes for the elongation factor G (EF-G) (Figure 4) (Table 2). FusA catalyzes the ribosomal translocation step during translation elongation, the mutations in fusA may avoid kanamycin binding EF-G and preventing translation [32]. Similar to ampC mutations in amoxicillin resistance, kanamycin-induced fusA mutations are drug target mutations that occur in every resistant strain, so their effects on fitness costs cannot be compared.
At the middle point of evolution, WT MG1655, ΔrpoS, and ΔsodA developed sbmA mutations (Figure 4). sbmA encodes a peptide antibiotic transporter, and mutations in this gene may reduce the permeability of antibiotics [33]. There were other mutated genes that occurred in ΔrpoS such as ubiA, oppD, and dxs, whose functions in antibiotic resistance were unclear, but these mutations may be responsible for the reduced growth rate and increased glucose consumption in the middle stage (Figure 1S and 2B). In addition to sbmA, another gene associated with reduced antibiotic uptake, trkH [34], was mutated in ΔsodA. In ΔsodB, trkH and other mutations emerged, potentially indicating a higher fitness burden than that of WT. In ΔsoxR and ΔsoxS, an antibiotic efflux associated gene kdpD was mutated earlier than in other strains [35], which may also be the reason why their fitness cost is higher than that of WT. Only ΔhipA, ΔhipB, and ΔrelA, did not have the antibiotic permeability and efflux related gene mutations at midpoint. However, a mutation in the spoT gene appeared in ΔrelA. spoT encodes another (p)ppGpp synthase [30]. Mutation of spoT potentially suggests that the bacterial stringent response to antibiotic stress is affected after knocking out relA. This in turn indicates that unlike other strains, which increase fitness cost due to genetic mutations, the relA knockout strain itself has a higher fitness cost in response to antibiotic stress, and it is not specific to a certain antibiotic. The absence of mutations related to antibiotic permeability and efflux at endpoint of ΔrelA kanamycin resistance further illustrates this point.

3. Discussion

In this study, we found that as resistance increases during amoxicillin or kanamycin evolution, the growth rate of resistant strains gradually decreases. This is to be expected because bacteria must spend the additional energy on either resistance mutations or on additional metabolic activity [9]. We tried to address the question of whether there is a common pattern for this increase in fitness cost, and whether other stress responses have synergistic or antagonistic effects. The answer we propose is that two types of growth inhibition patterns of resistant strains related to these two stress responses exist and that both the stringent response and oxidative stress response have an impact on the fitness cost during the evolution of resistance.
The reduced growth rates of ROS-related single-gene knockout strains ΔsodA, ΔsodB, ΔsoxR, and ΔsoxS (Figure 5) are due to mutations related to antibiotic permeability and efflux pumps that occur halfway through the evolution process. This conclusion concurs with the observation that strains with increased levels of ROS production have a faster rate of resistance development [14]. During the evolution of bactericidal antibiotic resistance, drug-target interactions stimulate NADH oxidation within the TCA cycle-dependent electron transport chain (ETC) [16]. The excessively produced superoxide compound in the ETC damages iron-sulfur clusters through the Fenton reaction, causing hydroxyl radical formation [36]. These superoxide and hydroxyl radicals, called ROS, can damage DNA [37]. In particular, guanine on the genome is oxidized by ROS to 8-hydroxy-2'-deoxyguanosine (8-HOdG), resulting in cytosine to thymine substitutions during DNA replication [38]. However, the non-lethal dose of ROS produced during non-lethal antibiotic exposure can still damage DNA and induce the formation of mutations [39]. This mutagenic effect is caused by activation of the cell's damage repair function through the SOS response [17]. The upregulation of the transcription level of error-prone DNA polymerase introduces errors during DNA damage repair [40]. These low-fidelity DNA polymerases can be considered as a tool used by bacteria to increase their environmental adaptability by increasing their mutation rate [41]. Mutations beneficial to antibiotic resistance are ultimately retained through selection. Single-gene knockout of superoxide dismutase SodA or SodB, or superoxide response regulon, SoxR or SoxS reduces the ability of cells to remove ROS, but increases the speed of resistance acquisition [14]. These earlier evolved mutations related to antibiotic permeability and efflux pumps all cost the cell more energy, ultimately leading to increased glucose consumption and reduced growth rate.
The growth rates of single-gene knockout strains related to the stringent response during the evolution process did not show any common features. This may be related to the different extents in which these strains trigger the alarmin (p)ppGpp under stringent responses. Only ΔrelA showed significant changes in both growth rate and biomass yield during the evolution of amoxicillin and kanamycin resistance. When relA is deleted, synthesis of (p)ppGpp is reduced, and amino acid starvation caused by the synthesis of resistance-related proteins cannot be offset in a timely manner [42]. The result is that deacylated tRNAs bind to ribosomes and hinder the translation elongation process [22]. Bacteria cannot promptly upregulate amino acid synthesis or proteolysis to increase the level of aminoacylated tRNA [43]. This will ultimately lead to reduced growth and increased energy requirements [44]. Moreover, the fitness burden caused by relA knockout is higher than that caused by specific resistance mutations. As a consequence, the growth rate of ΔrelA decreases at lower MIC levels than the growth rate of the WT (Figure 5A and B).
Two typical types of gene mutations were observed in antibiotic-resistant strains: Target-specific mutations, such as ampC upstream mutations in amoxicillin-resistant strains and fusA mutations in kanamycin-resistant strains; and non-target specific mutations, which mainly affects the intracellular antibiotic content, for example, reducing antibiotic uptake and increasing efflux. The mutations in the promoter region of ampC increase the beta-lactamase expression by upregulating its transcription level, which increases the energy consumption of bacteria and affects its growth rate [45]. Mutations of the EF-G encoding gene fusA also decrease the growth rate of bacteria by reducing the protein synthesis rate [46]. Similarly, mutations that affect intracellular antibiotic levels can also increase bacterial metabolic costs. Mutations related to outer membrane porins that reduce antibiotic permeability can also block the absorption of other beneficial compounds, including nutrients [47]. In addition to consuming extra energy, the activated efflux pumps can also export beneficial compounds out of the cell, ultimately increasing the metabolic burden of bacteria [48,49]. We speculate that this is why ΔrelA-resistant strains have not evolved mutations related to mechanisms such as antibiotic permeability and efflux pumps, because these mutations will cause the already nutrient-starved strains to become even more nutrient deficient.
During later stages of resistance evolution, almost all strains evolved mutations related to antibiotic permeability and efflux pumps. The final resistant strains all had significantly lower biomass yield on glucose compared to their ancestor. Moreover, the biomass yield of ΔrelA-resistant strains was significantly lower than that of other resistant strains. This further indicates that in the relA knockout, the ability to respond to stress is weakened due to limited synthesis of (p)ppGpp, resulting in a reduced growth rate and increased energy requirement to maintain reproduction. Furthermore, we found a slight increase in the growth rate of specific resistant strain during the later stages of amoxicillin, enrofloxacin, and kanamycin evolution. At this point we found mutations that have not been reported to be directly associated with antibiotic resistance, therefore, we hypothesize that this is the result of compensatory evolution, i.e., compensating for the metabolic burden caused by antibiotic-resistant adaptations [50,51]. Compensatory evolution allows bacteria to maintain antimicrobial-resistant properties, which is one of the reasons why resistant strains are difficult to completely eliminate [52,53,54].

4. Conclusion

During the evolution of resistance, the increase in MIC was positively correlated with a decrease in growth rate. Furthermore, in the middle stages of amoxicillin and kanamycin resistance evolution, resistant strains with reduced growth rates showed decreased biomass yield on glucose. In addition to the general decrease, we identified two distinct growth-limiting patterns that increase this effect through genetic mutation analysis during resistance development. First, the elevated mutation rate induced by oxidative stress in ROS scavenging-related gene knockout strains leads to the early emergence of non-specific resistance mutations, and these mutations are related to antibiotic permeability and pumping, which increase the physiological burden. Second, upon the knockout of the (p)ppGpp synthase gene relA, stringent response regulation became constrained, resulting in lower growth rates as a result of amino acid starvation. Genetic mutations not directly linked to antibiotic resistance are likely involved in compensatory evolution, ultimately contributing to the recovery of bacterial growth rates. This study shows how the stringent response and oxidative stress responses enhance the metabolic burdens, thereby reducing growth rates during development of resistance, indicating fitness costs associated with resistance development.

5. Materials and Methods

5.1. Bacterial Strains, Growth Media, and Culture Conditions

This study employed K-12 derivatives of E. coli, including the Keio Knockout strains ΔhipA, ΔhipB, ΔrelA, ΔrpoS, ΔsodA, ΔsodB, ΔsoxR, and ΔsoxS, as well as the almost identical wild-type MG1655 strain. Strains were cultured at 37 °C with shaking at 200 rpm in a phosphate-buffered minimal medium containing 100 mM NaH2PO4 and supplemented with 55 mM glucose [55].

5.2. MIC Determination

First, the MICs of the ancestral strains were measured. In a 96-well plate, strains with an initial OD600 of 0.05 were cultured in 150 μL of medium containing antibiotics at concentrations ranging from 0.5 to 4096 μg/mL in multiples of 2. The incubations took place in a spectrophotometer plate reader (Thermo Fisher Scientific) maintained at 37 °C with shaking. After 24 hours, the MIC was defined as the lowest antibiotic concentration at which the OD600 measurement value was below 0.2.
Subsequently, the ancestral strains of each bacterium were cultured in 5 mL tubes containing one of the four antibiotics amoxicillin, enrofloxacin, kanamycin, and tetracycline at one-quarter MIC. Through resistance evolution experiments, each strain gradually acquired resistance [14,15]. Specifically, after overnight incubation, if the OD600 of the culture containing antibiotics exceeded 75% of the OD600 of the culture without antibiotics, the bacterial strains from the antibiotic-containing culture were introduced into two fresh medium test tubes. One tube contained double the original antibiotic concentration, while the other retained the initial antibiotic concentration. After another day of incubation, if the OD600 of the culture with the higher antibiotic concentration surpassed 75% of the OD600 of the culture with the lower concentration, the bacterial culture with the higher antibiotic concentration was chosen. Otherwise, the culture with the lower antibiotic concentration was selected. This process was iteratively continued until the antibiotic concentration could no longer be increased. Throughout the resistance evolution, strains that developed increased antibiotic concentrations had their MICs measured using the 96-well plate method.

5.3. Growth Rate Measurements

During the evolution of resistance, strains measuring MIC were also tested for their growth rates in antibiotic-free conditions. The ancestor and resistant strains were cultured in 150 μL of medium with an initial OD600 of 0.05 in 96-well plates. Bacteria strains were grown overnight in a spectrophotometer plate reader and absorbance values were recorded every ten minutes. Growthrates-in-R (https://github.com/Pimutje/Growthrates-in-R/releases/tag/Growthrates) was used to plot growth curves and calculate growth rates. Each sample had at least three replicates. Data are presented as means ± SD, statistical significance was determined using a one-way ANOVA, *p < 0.05 **p < 0.01, ***p < 0.001.

5.4. Biomass Yield on Glucose Measurements

Glucose consumption was quantified using High-Performance Liquid Chromatography (HPLC) [56]. The ancestor and resistant strains were cultured in 5 mL of medium with an initial OD600 of 0.05 in test tubes, allowing growth until the OD600 reached 1.0. Subsequently, the supernatant was obtained by centrifuging at 12,000 rpm for 10 minutes and filtering it through 0.22 μm filters to prepare the samples. To make the calibration curve, glucose standards were prepared, ranging from 0 to 50 mM. These standards, along with the samples, were analyzed using an HPLC instrument (LC-20AT, Prominence, Shimadzu), which was equipped with a 300 × 7.8 mm Ion exclusion Rezex ROA Organic Acid H+ (8%) column (Phenomenex), a UV detector (SPD-20 A, 210 nm), and a refractive index detector (RID-20 A, 40 °C). The mobile phase consisted of 5 mM aqueous H2SO4, with a flow rate of 0.5 mL/min, and the analysis was conducted at 55 °C. The correlation of OD600 to cell dry weight (DW) was CDW = 0.56 × OD600 g/L [57]. The biomass yield on glucose was calculated by 0.56 g L-1/180.156 g mol-1 × Glucose consumption mol L-1 (gDW/gGlu). Each sample had at least three replicates. Data are presented as means ± SD, statistical significance was determined using a one-way ANOVA, *p < 0.05 **p < 0.01, ***p < 0.001.

5.5. Whole Genome Sequencing

Genomic DNA extracted from intermediate stages of amoxicillin (MIC = 128 µg/mL) and kanamycin (MIC = 256 µg/mL) resistance evolution, in addition to the final resistant strains and their respective ancestral strains, was processed using the DNeasy Blood and Tissue Kit (Qiagen). Subsequently, whole-genome sequencing was conducted employing the NextSeq 550 next-generation sequencing system (Illumina). Sequencing analysis adhered to established protocols [58], including read alignment to reference genomes via Bowtie2. Variant calling was executed employing Freebayes and Lofreq, with Snpeff facilitating variant annotation. Shared mutations between the resistant strains and their ancestral counterparts were excluded from the analysis. Specific single nucleotide polymorphisms (SNPs) and small insertions/deletions (indels) were documented for further investigation.

5.6. Quantification and Statistical Analysis

Statistical analysis employed IBM SPSS software. Specific statistical methods for each experiment are available in the figure legends and respective methods section.

Author Contributions

WQ and BtK conceived the project. MW assisted in the design of experiments. WQ performed experiments and analysis of the data. MJ, DK and WL performed the bioinformatic analysis. WQ and BtK wrote the manuscript. All authors critically reviewed the manuscript and approved the final version.

Acknowledgments

We thank Stanley Brul for stimulating discussions and Selina van Leeuwen for her assistance with the DNA and RNA sequencing. This study was financed by the Netherlands Food and Consumer Product Safety Authority (NVWA). The NVWA was not involved in design of the experiments, analysis of the data or writing the manuscript. Wenxi Qi acknowledges the China Scholarship Council for his PhD scholarship.

Conflicts of Interest

The authors declare no competing interests.

References

  1. Murray, C.J.; Ikuta, K.S.; Sharara, F.; Swetschinski, L.; Robles Aguilar, G.; Gray, A.; Han, C.; Bisignano, C.; Rao, P.; Wool, E.; et al. Global Burden of Bacterial Antimicrobial Resistance in 2019: A Systematic Analysis. Lancet 2022, 399, 629–655. [Google Scholar] [CrossRef] [PubMed]
  2. Llor, C.; Bjerrum, L. Antimicrobial Resistance: Risk Associated with Antibiotic Overuse and Initiatives to Reduce the Problem. Ther. Adv. drug Saf. 2014, 5, 229–241. [Google Scholar] [CrossRef] [PubMed]
  3. Mendelsohn, E.; Ross, N.; Zambrana-Torrelio, C.; Van Boeckel, T.P.; Laxminarayan, R.; Daszak, P. Global Patterns and Correlates in the Emergence of Antimicrobial Resistance in Humans. Proceedings. Biol. Sci. 2023, 290, 20231085. [Google Scholar] [CrossRef]
  4. Andersson, D.I.; Hughes, D. Antibiotic Resistance and Its Cost: Is It Possible to Reverse Resistance? Nat. Rev. Microbiol. 2010, 8, 260–271. [Google Scholar] [CrossRef]
  5. Pinheiro, F.; Warsi, O.; Andersson, D.I.; Lässig, M. Metabolic Fitness Landscapes Predict the Evolution of Antibiotic Resistance. Nat. Ecol. Evol. 2021, 5, 677–687. [Google Scholar] [CrossRef]
  6. Nishimoto, A.T.; Dao, T.H.; Jia, Q.; Ortiz-Marquez, J.C.; Echlin, H.; Vogel, P.; van Opijnen, T.; Rosch, J.W. Interspecies Recombination, Not de Novo Mutation, Maintains Virulence after β-Lactam Resistance Acquisition in Streptococcus Pneumoniae. Cell Rep. 2022, 41, 111835. [Google Scholar] [CrossRef]
  7. Ogunlana, L.; Kaur, D.; Shaw, L.P.; Jangir, P.; Walsh, T.; Uphoff, S.; MacLean, R.C. Regulatory Fine-Tuning of Mcr-1 Increases Bacterial Fitness and Stabilises Antibiotic Resistance in Agricultural Settings. ISME J. 2023, 17, 2058–2069. [Google Scholar] [CrossRef]
  8. Vanacker, M.; Lenuzza, N.; Rasigade, J.P. The Fitness Cost of Horizontally Transferred and Mutational Antimicrobial Resistance in Escherichia Coli. Front. Microbiol. 2023, 14, 1186920. [Google Scholar] [CrossRef] [PubMed]
  9. Händel, N.; Schuurmans, J.M.; Brul, S.; Ter Kuilea, B.H. Compensation of the Metabolic Costs of Antibiotic Resistance by Physiological Adaptation in Escherichia Coli. Antimicrob. Agents Chemother. 2013, 57, 3752–3762. [Google Scholar] [CrossRef]
  10. Andersson, D.I. The Biological Cost of Mutational Antibiotic Resistance: Any Practical Conclusions? Curr. Opin. Microbiol. 2006, 9, 461–465. [Google Scholar] [CrossRef]
  11. Baker-Austin, C.; Wright, M.S.; Stepanauskas, R.; McArthur, J. V. Co-Selection of Antibiotic and Metal Resistance. Trends Microbiol. 2006, 14, 176–182. [Google Scholar] [CrossRef] [PubMed]
  12. Darphorn, T.S.; Koenders-Van Sintanneland, B.B.; Grootemaat, A.E.; van der Wel, N.N.; Brul, S.; ter Kuile, B.H. Transfer Dynamics of Multi-Resistance Plasmids in Escherichia Coli Isolated from Meat. PLoS One 2022, 17, e0270205. [Google Scholar] [CrossRef] [PubMed]
  13. Dawan, J.; Ahn, J. Bacterial Stress Responses as Potential Targets in Overcoming Antibiotic Resistance. Microorganisms 2022, 10, 1385. [Google Scholar] [CrossRef] [PubMed]
  14. Qi, W.; Jonker, M.J.; de Leeuw, W.; Brul, S.; ter Kuile, B.H. Reactive Oxygen Species Accelerate de Novo Acquisition of Antibiotic Resistance in E. Coli. iScience 2023, 108373. [Google Scholar] [CrossRef]
  15. Qi, W.; Jonker, M.J.; Leeuw, W. de; Brul, S.; Kuile, B.H. ter Role of RelA-Synthesized (p)PpGpp in de Novo Acquisition of Antibiotic Resistance in E. Coli. Res. Sq. 2023. [Google Scholar] [CrossRef]
  16. Kohanski, M.A.; Dwyer, D.J.; Hayete, B.; Lawrence, C.A.; Collins, J.J. A Common Mechanism of Cellular Death Induced by Bactericidal Antibiotics. Cell 2007, 130, 797–810. [Google Scholar] [CrossRef]
  17. Pribis, J.P.; García-Villada, L.; Zhai, Y.; Lewin-Epstein, O.; Wang, A.Z.; Liu, J.; Xia, J.; Mei, Q.; Fitzgerald, D.M.; Bos, J.; et al. Gamblers: An Antibiotic-Induced Evolvable Cell Subpopulation Differentiated by Reactive-Oxygen-Induced General Stress Response. Mol. Cell 2019, 74, 785–800. [Google Scholar] [CrossRef] [PubMed]
  18. Kaspy, I.; Rotem, E.; Weiss, N.; Ronin, I.; Balaban, N.Q.; Glaser, G. HipA-Mediated Antibiotic Persistence via Phosphorylation of the Glutamyl-TRNA-Synthetase. Nat. Commun. 2013, 4. [Google Scholar] [CrossRef]
  19. Germain, E.; Castro-Roa, D.; Zenkin, N.; Gerdes, K. Molecular Mechanism of Bacterial Persistence by HipA. Mol. Cell 2013, 52, 248–254. [Google Scholar] [CrossRef]
  20. Magnusson, L.U.; Farewell, A.; Nyström, T. PpGpp: A Global Regulator in Escherichia Coli. Trends Microbiol. 2005, 13, 236–242. [Google Scholar] [CrossRef]
  21. Schumacher, M.A.; Piro, K.M.; Xu, W.; Hansen, S.; Lewis, K.; Brennan, R.G. Molecular Mechanisms of HipA-Mediated Multidrug Tolerance and Its Neutralization by HipB. Science (80-. ). 2009, 323, 396–401. [Google Scholar] [CrossRef]
  22. Arenz, S.; Abdelshahid, M.; Sohmen, D.; Payoe, R.; Starosta, A.L.; Berninghausen, O.; Hauryliuk, V.; Beckmann, R.; Wilson, D.N. The Stringent Factor RelA Adopts an Open Conformation on the Ribosome to Stimulate PpGpp Synthesis. Nucleic Acids Res. 2016, 44, 6471–6481. [Google Scholar] [CrossRef] [PubMed]
  23. Dalebroux, Z.D.; Swanson, M.S. PpGpp: Magic beyond RNA Polymerase. Nat. Rev. Microbiol. 2012, 10, 203–212. [Google Scholar] [CrossRef] [PubMed]
  24. Dwyer, D.J.; Kohanski, M.A.; Collins, J.J. Role of Reactive Oxygen Species in Antibiotic Action and Resistance. Curr. Opin. Microbiol. 2009, 12, 482–489. [Google Scholar] [CrossRef] [PubMed]
  25. Imlay, J.A. The Molecular Mechanisms and Physiological Consequences of Oxidative Stress: Lessons from a Model Bacterium. Nat. Rev. Microbiol. 2013, 11, 443–454. [Google Scholar] [CrossRef]
  26. Darphorn, T.S.; Hu, Y.; Sintanneland, B.B.K.; Brul, S.; Kuile, B.H. ter Multiplication of AmpC upon Exposure to a Beta-Lactam Antibiotic Results in a Transferable Transposon in Escherichia Coli. Int. J. Mol. Sci. 2021, 22, 9230. [Google Scholar] [CrossRef] [PubMed]
  27. Zhao, Z.; Xu, Y.; Jiang, B.; Qi, Q.; Tang, Y.-J.; Xian, M.; Wang, J.; Zhao, G. Systematic Identification of CpxRA-Regulated Genes and Their Roles in Escherichia Coli Stress Response. mSystems 2022, 7, e0041922. [Google Scholar] [CrossRef] [PubMed]
  28. Cai, S.J.; Inouye, M. EnvZ-OmpR Interaction and Osmoregulation in Escherichia Coli. J. Biol. Chem. 2002, 277, 24155–24161. [Google Scholar] [CrossRef] [PubMed]
  29. Choi, U.; Lee, C.R. Distinct Roles of Outer Membrane Porins in Antibiotic Resistance and Membrane Integrity in Escherichia Coli. Front. Microbiol. 2019, 10, 953. [Google Scholar] [CrossRef]
  30. Hobbs, J.K.; Boraston, A.B. (P)PpGpp and the Stringent Response: An Emerging Threat to Antibiotic Therapy. ACS Infect. Dis. 2019, 5, 1505–1517. [Google Scholar] [CrossRef]
  31. Wu, L.; Wang, Z.; Guan, Y.; Huang, X.; Shi, H.; Liu, Y.; Zhang, X. The (p)PpGpp-Mediated Stringent Response Regulatory System Globally Inhibits Primary Metabolism and Activates Secondary Metabolism in Pseudomonas Protegens H78. Appl. Microbiol. Biotechnol. 2020, 104, 3061–3079. [Google Scholar] [CrossRef] [PubMed]
  32. Feldman, M.B.; Terry, D.S.; Altman, R.B.; Blanchard, S.C. Aminoglycoside Activity Observed on Single Pre-Translocation Ribosome Complexes. Nat. Chem. Biol. 2010, 6, 54–62. [Google Scholar] [CrossRef] [PubMed]
  33. Mogre, A.; Veetil, R.T.; Seshasayee, A.S.N. Modulation of Global Transcriptional Regulatory Networks as a Strategy for Increasing Kanamycin Resistance of the Translational Elongation Factor-G Mutants in Escherichia Coli. G3 Genes, Genomes, Genet. 2017, 7, 3955–3966. [Google Scholar] [CrossRef]
  34. Oz, T.; Guvenek, A.; Yildiz, S.; Karaboga, E.; Tamer, Y.T.; Mumcuyan, N.; Ozan, V.B.; Senturk, G.H.; Cokol, M.; Yeh, P.; et al. Strength of Selection Pressure Is an Important Parameter Contributing to the Complexity of Antibiotic Resistance Evolution. Mol. Biol. Evol. 2014, 31, 2387. [Google Scholar] [CrossRef] [PubMed]
  35. Freeman, Z.N.; Dorus, S.; Waterfield, N.R. The KdpD/KdpE Two-Component System: Integrating K+ Homeostasis and Virulence. PLoS Pathog. 2013, 9, e1003201. [Google Scholar] [CrossRef] [PubMed]
  36. Dwyer, D.J.; Kohanski, M.A.; Hayete, B.; Collins, J.J. Gyrase Inhibitors Induce an Oxidative Damage Cellular Death Pathway in Escherichia Coli. Mol. Syst. Biol. 2007, 3, 91. [Google Scholar] [CrossRef] [PubMed]
  37. Dwyer, D.J.; Belenky, P.A.; Yang, J.H.; Cody MacDonald, I.; Martell, J.D.; Takahashi, N.; Chan, C.T.Y.; Lobritz, M.A.; Braff, D.; Schwarz, E.G.; et al. Antibiotics Induce Redox-Related Physiological Alterations as Part of Their Lethality. Proc. Natl. Acad. Sci. U. S. A. 2014, 111, E2100–E2109. [Google Scholar] [CrossRef] [PubMed]
  38. Foti, J.J.; Devadoss, B.; Winkler, J.A.; Collins, J.J.; Walker, G.C. Oxidation of the Guanine Nucleotide Pool Underlies Cell Death by Bactericidal Antibiotics. Science (80-. ). 2012, 336, 315–319. [Google Scholar] [CrossRef]
  39. Kohanski, M.A.; DePristo, M.A.; Collins, J.J. Sublethal Antibiotic Treatment Leads to Multidrug Resistance via Radical-Induced Mutagenesis. Mol. Cell 2010, 37, 311–320. [Google Scholar] [CrossRef] [PubMed]
  40. Goodman, M.F.; Mcdonald, J.P.; Jaszczur, M.M.; Woodgate, R. Insights into the Complex Levels of Regulation Imposed on Escherichia Coli DNA Polymerase V. DNA Repair (Amst). 2016, 44, 42–50. [Google Scholar] [CrossRef]
  41. Crane, J.K.; Alvarado, C.L.; Sutton, M.D. Role of the SOS Response in the Generation of Antibiotic Resistance in Vivo. Antimicrob. Agents Chemother. 2021, 65, e0001321. [Google Scholar] [CrossRef]
  42. Sinha, A.K.; Winther, K.S. The RelA Hydrolase Domain Acts as a Molecular Switch for (p)PpGpp Synthesis. Commun. Biol. 2021, 4, 434. [Google Scholar] [CrossRef]
  43. Irving, S.E.; Choudhury, N.R.; Corrigan, R.M. The Stringent Response and Physiological Roles of (Pp)PGpp in Bacteria. Nat. Rev. Microbiol. 2020, 19, 256–271. [Google Scholar] [CrossRef] [PubMed]
  44. Hauryliuk, V.; Atkinson, G.C.; Murakami, K.S.; Tenson, T.; Gerdes, K. Recent Functional Insights into the Role of (p)PpGpp in Bacterial Physiology. Nat. Rev. Microbiol. 2015, 13, 298–309. [Google Scholar] [CrossRef]
  45. Pérez-Gallego, M.; Torrens, G.; Castillo-Vera, J.; Moya, B.; Zamorano, L.; Cabot, G.; Hultenby, K.; Albertí, S.; Mellroth, P.; Henriques-Normark, B.; et al. Impact of AmpC Derepression on Fitness and Virulence: The Mechanism or the Pathway? MBio 2016, 7, e01783–16. [Google Scholar] [CrossRef] [PubMed]
  46. Macvanin, M.; Björkman, J.; Eriksson, S.; Rhen, M.; Andersson, D.I.; Hughes, D. Fusidic Acid-Resistant Mutants of Salmonella Enterica Serovar Typhimurium with Low Fitness in Vivo Are Defective in RpoS Induction. Antimicrob. Agents Chemother. 2003, 47, 3743–3749. [Google Scholar] [CrossRef]
  47. Knopp, M.; Andersson, D.I. Amelioration of the Fitness Costs of Antibiotic Resistance Due To Reduced Outer Membrane Permeability by Upregulation of Alternative Porins. Mol. Biol. Evol. 2015, 32, 3252–3263. [Google Scholar] [CrossRef]
  48. Nicoloff, H.; Perreten, V.; McMurry, L.M.; Levy, S.B. Role for Tandem Duplication and Lon Protease in AcrAB-TolC- Dependent Multiple Antibiotic Resistance (Mar) in an Escherichia Coli Mutant without Mutations in MarRAB or AcrRAB. J. Bacteriol. 2006, 188, 4413–4423. [Google Scholar] [CrossRef] [PubMed]
  49. Keeney, D.; Ruzin, A.; Mcaleese, F.; Murphy, E.; Bradford, P.A. MarA-Mediated Overexpression of the AcrAB Efflux Pump Results in Decreased Susceptibility to Tigecycline in Escherichia Coli. J. Antimicrob. Chemother. 2008, 61, 46–53. [Google Scholar] [CrossRef] [PubMed]
  50. Hernando-Amado, S.; Sanz-García, F.; Blanco, P.; Martínez, J.L. Fitness Costs Associated with the Acquisition of Antibiotic Resistance. Essays Biochem. 2017, 61, 37–48. [Google Scholar] [CrossRef]
  51. Lind, P.A.; Tobin, C.; Berg, O.G.; Kurland, C.G.; Andersson, D.I. Compensatory Gene Amplification Restores Fitness after Inter-Species Gene Replacements. Mol. Microbiol. 2010, 75, 1078–1089. [Google Scholar] [CrossRef]
  52. Andersson, D.I.; Hughes, D. Antibiotic Resistance and Its Cost: Is It Possible to Reverse Resistance? Nat. Rev. Microbiol. 2010, 8, 260–271. [Google Scholar] [CrossRef] [PubMed]
  53. Zampieri, M.; Enke, T.; Chubukov, V.; Ricci, V.; Piddock, L.; Sauer, U. Metabolic Constraints on the Evolution of Antibiotic Resistance. Mol. Syst. Biol. 2017, 13, 917. [Google Scholar] [CrossRef]
  54. Durão, P.; Balbontín, R.; Gordo, I. Evolutionary Mechanisms Shaping the Maintenance of Antibiotic Resistance. Trends Microbiol. 2018, 26, 677–691. [Google Scholar] [CrossRef] [PubMed]
  55. Evans, C.G.T.; Herbert, D.; Tempest, D.W. Chapter XIII The Continuous Cultivation of Micro-Organisms: 2. Construction of a Chemostat. Methods Microbiol. 1970, 2, 277–327. [Google Scholar] [CrossRef]
  56. Chen, T.; Brul, S.; Hugenholtz, J. Exploring the Potential of Bacillus Subtilis as Cell Factory for Food Ingredients and Special Chemicals. Microb. Cell Fact. 2023, 22, 200. [Google Scholar] [CrossRef] [PubMed]
  57. Link, H.; Anselment, B.; Weuster-Botz, D. Leakage of Adenylates during Cold Methanol/Glycerol Quenching of Escherichia Coli. Metabolomics 2008, 4, 240–247. [Google Scholar] [CrossRef]
  58. Qi, W.; Jonker, M.J.; Teichmann, L.; Wortel, M.; ter Kuile, B.H. The Influence of Oxygen and Oxidative Stress on de Novo Acquisition of Antibiotic Resistance in E. Coli and Lactobacillus Lactis. BMC Microbiol. 2023, 23, 279. [Google Scholar] [CrossRef]
Figure 1. Growth rates of ROS- or (p)ppGpp-related single-gene knockout strains during resistance evolution. (A-P) The growth rates of ROS- or (p)ppGpp-related amoxicillin- (A, E, I, and M), enrofloxacin- (B, F, J, and N), kanamycin- (C, G, K, and O), and tetracycline- (D, H, L, and P) resistant evolution at each minimum inhibitory concentration (MIC). The x-axis represents the MIC, while the y-axis represents the growth rate. The linear regression equation between the log of MIC and the growth rate and R squared (R2) of each mutant are shown in each figure. Data are presented as means ± SD, N ≥ 3. (Q-S) The growth rates comparison between the ancestor single-gene knockout strains (Q), the middle resistance evolution point (MIC = 128 µg/mL) of amoxicillin-resistant strains (R), and the middle resistance evolution point (MIC = 256 µg/mL) of kanamycin-resistant strains (R). Data are presented as means ± SD, statistical significance was determined using a one-way ANOVA, N ≥ 3, *p < 0.05, **p < 0.01, ***p < 0.001.
Figure 1. Growth rates of ROS- or (p)ppGpp-related single-gene knockout strains during resistance evolution. (A-P) The growth rates of ROS- or (p)ppGpp-related amoxicillin- (A, E, I, and M), enrofloxacin- (B, F, J, and N), kanamycin- (C, G, K, and O), and tetracycline- (D, H, L, and P) resistant evolution at each minimum inhibitory concentration (MIC). The x-axis represents the MIC, while the y-axis represents the growth rate. The linear regression equation between the log of MIC and the growth rate and R squared (R2) of each mutant are shown in each figure. Data are presented as means ± SD, N ≥ 3. (Q-S) The growth rates comparison between the ancestor single-gene knockout strains (Q), the middle resistance evolution point (MIC = 128 µg/mL) of amoxicillin-resistant strains (R), and the middle resistance evolution point (MIC = 256 µg/mL) of kanamycin-resistant strains (R). Data are presented as means ± SD, statistical significance was determined using a one-way ANOVA, N ≥ 3, *p < 0.05, **p < 0.01, ***p < 0.001.
Preprints 90087 g001
Figure 2. The biomass yield on glucose of the ancestor strains and amoxicillin- or kanamycin-resistant strains. The biomass yield on glucose of amoxicillin-resistant (A) and kanamycin-resistant (B) strains at the start, middle and final resistance evolution points in antibiotic-free medium. Data are presented as means ± SD, statistical significance was determined using a one-way ANOVA, N = 3, *p < 0.05 **p < 0.01, ***p < 0.001.
Figure 2. The biomass yield on glucose of the ancestor strains and amoxicillin- or kanamycin-resistant strains. The biomass yield on glucose of amoxicillin-resistant (A) and kanamycin-resistant (B) strains at the start, middle and final resistance evolution points in antibiotic-free medium. Data are presented as means ± SD, statistical significance was determined using a one-way ANOVA, N = 3, *p < 0.05 **p < 0.01, ***p < 0.001.
Preprints 90087 g002
Figure 3. Mutations in the amoxicillin-resistant strains during the resistance evolution. The mutated genes and their mutations at the middle (gray colors) resistance evolution point and the end (black colors) resistance evolution point during the resistance development. The x-axis represents the population frequency, while the y-axis indicates the mutations. The different colors of the mutated genes represent different functions according to the Comprehensive Antibiotic Resistance Database and UniProt. Orange color means genes associated with reduced antibiotic permeability, blue color signifies genes associated with antibiotic efflux pumps, red color indicates genes associated with antibiotic target alteration, purple color means genes associated with antibiotic inactivation, and black color means genes not directly associated with antibiotic resistance.
Figure 3. Mutations in the amoxicillin-resistant strains during the resistance evolution. The mutated genes and their mutations at the middle (gray colors) resistance evolution point and the end (black colors) resistance evolution point during the resistance development. The x-axis represents the population frequency, while the y-axis indicates the mutations. The different colors of the mutated genes represent different functions according to the Comprehensive Antibiotic Resistance Database and UniProt. Orange color means genes associated with reduced antibiotic permeability, blue color signifies genes associated with antibiotic efflux pumps, red color indicates genes associated with antibiotic target alteration, purple color means genes associated with antibiotic inactivation, and black color means genes not directly associated with antibiotic resistance.
Preprints 90087 g003
Figure 4. Mutations in the kanamycin-resistant strains during the resistance evolution. The mutated genes and their mutations at the middle (gray colors) resistance evolution point and the end (black colors) resistance evolution point during the resistance development. The x-axis represents the population frequency, while the y-axis indicates the mutations. The different colors of the mutated genes represent different functions according to the Comprehensive Antibiotic Resistance Database and UniProt. Orange color indicates genes associated with reduced antibiotic permeability, blue color means genes associated with antibiotic efflux pumps, red color denotes genes associated with antibiotic target alteration, and black color represents genes not directly associated with antibiotic resistance.
Figure 4. Mutations in the kanamycin-resistant strains during the resistance evolution. The mutated genes and their mutations at the middle (gray colors) resistance evolution point and the end (black colors) resistance evolution point during the resistance development. The x-axis represents the population frequency, while the y-axis indicates the mutations. The different colors of the mutated genes represent different functions according to the Comprehensive Antibiotic Resistance Database and UniProt. Orange color indicates genes associated with reduced antibiotic permeability, blue color means genes associated with antibiotic efflux pumps, red color denotes genes associated with antibiotic target alteration, and black color represents genes not directly associated with antibiotic resistance.
Preprints 90087 g004
Figure 5. Schematic diagram of two processes causing reduced growth rates as resistance increases. As the MICs for amoxicillin (A) and kanamycin (B) increase, the growth rates of all strains gradually decrease. This relates to two different processes. ROS-related knockout strains develop mutations related to antibiotic permeability and efflux pumps earlier (in mid-stage) than other strains. Because these mutations increase energy consumption, the fitness cost increases and the growth rate decreases. The other process is that in the early stages of resistance evolution, ΔrelA shows a decrease in growth rate, which may be related to it simulating starvation-induced growth arrest. Different colored dots represent the median of the growth rate of each strain at different MICs. Gray shading indicates the range of variation in WT growth rates. The green line connects the median of the growth rates of ΔrelA.
Figure 5. Schematic diagram of two processes causing reduced growth rates as resistance increases. As the MICs for amoxicillin (A) and kanamycin (B) increase, the growth rates of all strains gradually decrease. This relates to two different processes. ROS-related knockout strains develop mutations related to antibiotic permeability and efflux pumps earlier (in mid-stage) than other strains. Because these mutations increase energy consumption, the fitness cost increases and the growth rate decreases. The other process is that in the early stages of resistance evolution, ΔrelA shows a decrease in growth rate, which may be related to it simulating starvation-induced growth arrest. Different colored dots represent the median of the growth rate of each strain at different MICs. Gray shading indicates the range of variation in WT growth rates. The green line connects the median of the growth rates of ΔrelA.
Preprints 90087 g005
Table 1. Number of times that frequently mutated genes appear in 9 amoxicillin-resistant strains.
Table 1. Number of times that frequently mutated genes appear in 9 amoxicillin-resistant strains.
Times 2 3 4 5 6 7 8 9
Mid ompC
envZ
rpoD cpxA ampC
End ompX
mipA
rpoA
ftsI
prlF
rpoD cpxA envZ ompC ampC
Table 2. Number of times that frequently mutated genes appear in 9 kanamycin-resistant strains.
Table 2. Number of times that frequently mutated genes appear in 9 kanamycin-resistant strains.
Times 2 3 4 5 6 7 8 9
Mid trkH
kdpD
oppA
sbmA
oppD
fusA
End rpsL
oppB
pgsA
trkH
oppA
oppD
kdpD atpG sbmA fusA
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.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated