Cellular hysteresis as a principle to maximize the efficacy of antibiotic therapy
Edited by Shimon Bershtein, Ben-Gurion University, Beer-Sheva, Israel, and accepted by Editorial Board Member Daniel L. Hartl August 16, 2018 (received for review June 10, 2018)
Significance
Rapid evolution is central to the current antibiotic crisis. Sustainable treatments must thus take account of the bacteria’s potential for adaptation. We identified cellular hysteresis as a principle to constrain bacterial evolution. Cellular hysteresis is a persistent change in bacterial physiology, reminiscent of cellular memory, which is induced by one antibiotic and enhances susceptibility toward another antibiotic. Cellular hysteresis increases bacterial extinction in fast sequential treatments and reduces selection of resistance by favoring responses specific to the induced physiological effects. Fast changes between antibiotics are key, because they create the continuously high selection conditions that are difficult to counter by bacteria. Our study highlights how an understanding of evolutionary processes can help to outsmart human pathogens.
Abstract
Antibiotic resistance has become one of the most dramatic threats to global health. While novel treatment options are urgently required, most attempts focus on finding new antibiotic substances. However, their development is costly, and their efficacy is often compromised within short time periods due to the enormous potential of microorganisms for rapid adaptation. Here, we developed a strategy that uses the currently available antibiotics. Our strategy exploits cellular hysteresis, which is the long-lasting, transgenerational change in cellular physiology that is induced by one antibiotic and sensitizes bacteria to another subsequently administered antibiotic. Using evolution experiments, mathematical modeling, genomics, and functional genetic analysis, we demonstrate that sequential treatment protocols with high levels of cellular hysteresis constrain the evolving bacteria by (i) increasing extinction frequencies, (ii) reducing adaptation rates, and (iii) limiting emergence of multidrug resistance. Cellular hysteresis is most effective in fast sequential protocols, in which antibiotics are changed within 12 h or 24 h, in contrast to the less frequent changes in cycling protocols commonly implemented in hospitals. We found that cellular hysteresis imposes specific selective pressure on the bacteria that disfavors resistance mutations. Instead, if bacterial populations survive, hysteresis is countered in two distinct ways, either through a process related to antibiotic tolerance or a mechanism controlled by the previously uncharacterized two-component regulator CpxS. We conclude that cellular hysteresis can be harnessed to optimize antibiotic therapy, to achieve both enhanced bacterial elimination and reduced resistance evolution.
Sign up for PNAS alerts.
Get alerts for new articles, or get an alert when an article is cited.
Natural environments are often temporally dynamic. They produce continuously changing selective constraints that are a particular challenge for organisms to adapt to (1). Similar dynamic conditions may be used in human therapy to limit the ability of pathogens for resistance evolution. Antibiotic resistance is a global threat (2), enhanced by the ongoing emergence of new resistance mechanisms (3, 4). Sequential treatments may be one option to counter emerging resistance: The alternation of antibiotics may produce the challenging selection conditions that are known from many fluctuating natural environments. Indeed, previous evolution experiments consistently identified reduced resistance evolution and increased extinction at sublethal drug concentrations in sequential treatment protocols (5–8). The efficacy of sequential treatments has generally been attributed to the evolution of collateral sensitivity (9, 10), i.e., genetic trade-offs. Collateral sensitivity describes the case where resistance mutations pleiotropically increase sensitivity to other antibiotics. Intriguingly, collateral sensitivity only occasionally explains the observed dynamics of resistance evolution (5), indicating the activity of unknown selective pressures.
An alternative selective pressure may be caused by antibiotic-induced physiological changes in bacterial cells. Such inducible changes can affect motility (11), biofilm formation (12), and heat shock response (13). Moreover, certain bacterial responses can be stabilized across generations. For example, the lactose utilization response in Escherichia coli is maintained for >10 generations after removal of lactose, by the inheritance of stable proteins (14). Positive autoregulation, a common feature of two-component sensors, can similarly stabilize responses for several days, as observed with PhoQ/PhoP (15). Upon specific environmental change, these inherited responses can have deleterious fitness effects. For example, salt-induced expression of the PhoE porin increases sensitivity to acid (16). Furthermore, structural damage from past stress encounters may be amplified in new environments. We here use the term “cellular hysteresis” (17) to summarize the inducible and long-lasting physiological effects on cellular integrity and function, spanning both discrete and quantitative changes. The term is more general than phenotypic memory, as it comprises positive and negative fitness effects induced by previous exposures. Although cellular hysteresis likely determines pathogen survival during sequential antibiotic treatments, it is not part of current treatment concepts. Consideration of negative hysteresis may enhance treatment efficacy when exposure to one antibiotic temporarily increases susceptibility to a subsequently administered second antibiotic.
The objectives of our study were to test the potential of negative hysteresis to enhance efficacy of sequential therapy. We used the pathogen Pseudomonas aeruginosa as a model—the second-most critical threat of a multidrug-resistant pathogen (18). We first characterized the hysteresis landscape for P. aeruginosa for three clinically relevant, bactericidal antibiotics with distinct cellular targets: ciprofloxacin (CIP), gentamicin (GEN), and carbenicillin (CAR). We then investigated how different levels of cellular hysteresis modulate the evolutionary adaptive response to sequential treatment using high-throughput experimental evolution, mathematical modeling, whole-genome sequencing, and functional genetic analysis of potential targets of selection. Finally, we validated the potential of antibiotic hysteresis for the inhibition and prediction of resistance evolution by second-order experimental evolution.
Results and Discussion
Cellular Hysteresis Depends on the Order of Drug Switches.
We determined how short exposures to nonlethal concentrations of the three considered antibiotics CIP, GEN, and CAR affected later antibiotic treatment, by performing time-kill experiments (Fig. 1A). The inferred hysteresis landscape included negative, positive, and neutral effects and showed strong directionality. For the drug pair CAR/GEN, the sign of hysteresis effects was dependent on drug order (Fig. 1B): Preexposure to GEN, which is known to cause translational stress, protected cells from killing by CAR, while preexposure to CAR, which causes cell envelope stress, increased bactericidal activity of GEN (Fig. 1 A–C and Movie S1). Similar directionality was observed for the drug pair CIP/GEN (Fig. 1B). Neutral or mild hysteresis effects were observed for the pair CIP/CAR. Altogether, the hysteresis landscape indicated that population survival during the antibiotic switch is strongly history-dependent, thus emphasizing the importance of drug order for the design of effective treatment. Order dependence and the observation that negative hysteresis can be achieved with drug pairs known to interact either synergistically (GEN+CAR) or antagonistically (CIP+GEN) (19) demonstrate that hysteresis and drug interaction are not necessarily linked. Hysteresis is also distinct from collateral sensitivity, as it occurs immediately, i.e., without acquiring resistance mutations.
Movie S1.
Fig. 1.
Further analysis of negative hysteresis caused by the CAR → GEN switch revealed that pretreatments as short as 15 min (equivalent to 1/3 generation time) and concentrations as low as 3 µg/mL [1/32 of the minimal inhibitory concentration (MIC)] were sufficient for the induction of accelerated killing (Fig. 1 D and E). Bactericidal activity increased 400-fold when the 15-min pretreatments were performed with higher concentrations, i.e., 2× MIC. These findings indicate a specific and robust physiological effect that may be exploited to increase the efficacy of antibiotic therapy, for example, by using sequential treatment. The increased killing is likely explained by β-lactam−induced acceleration of cellular influx of the aminoglycoside (20). The reduced killing in the reversed direction may be explained by aminoglycoside-induced up-regulation of efflux pumps (21, 22).
Negative Hysteresis Increases Treatment Efficacy and Constrains Evolution of Resistance.
By modulating bacterial killing, we hypothesized that cellular hysteresis influences the rate of resistance evolution during sequential treatment dependent on drug order and the resulting cumulative level of hysteresis. To test these predictions, we conducted a high-throughput evolution experiment with 190 replicate populations over a total of 96 transfers (each with 12-h growth intervals; total of ∼500 generations; Fig. 2A and SI Appendix, Fig. S2). We included three main types of sequential protocols to disentangle the influence of hysteresis from the frequency and regularity of drug switches (Fig. 2B, columns 1–3). Evolutionary dynamics were tightly monitored, with continuous measurements of population growth (every 15 min). These measurements revealed that the evolutionary dynamics were separated into three main phases across the 96 transfers (Fig. 3A): (i) an initial phase up to roughly transfer 12, during which populations adapted rapidly and treatments varied strongly in evolutionary dynamics; (ii) an optimization phase from transfer 12 up to approximately transfer 48, during which two treatments (monotherapy, slow regular protocol) were almost fully adapted, while the other two (fast regular and random protocols) still produced increases in growth, yet at lower rates as during the first phase; and (iii) the long-term dynamics from transfer 48 onward, during which only small growth increases and little variation among treatments were observed. For a more detailed analysis, we focused on the early dynamics up to transfer 12, because these encompass the strongest differences in the tested variables and cover a clinically relevant time period of 6 d.
Fig. 2.
Fig. 3.
For the early dynamics, fast and random sequences led to significant improvements in three independently characterized, complementary measures for treatment efficacy (Fig. 2B; see SI Appendix, Tables S1–S6 for statistics), including (i) higher extinction frequencies (inferred from absence of growth during experimental evolution), (ii) lower adaptation rates (calculated from growth characteristics measured during the evolution experiment), and (iii) lower levels of evolved multidrug resistance (MDR, determined from antibiotic dose–response curves for individual bacterial clones isolated from the evolving populations; SI Appendix, Fig. S3). Importantly, the cumulative levels of hysteresis were significantly correlated to both adaptation rates and evolved MDR (Spearman rank correlation, rho ≥ 0.73, P ≤ 0.01; Fig. 2C). Evolved MDR was also significantly associated with switching rate, but to a lesser degree (Spearman rank correlation, rho = 0.66, P = 0.019), while there was no significant relationship between switching and adaptation rates (SI Appendix, Fig. S5). Increases in both negative hysteresis level and switching rate led to higher extinction frequencies, even though the effect was not statistically significant. We conclude that, even though switching rate is important, the consideration of negative hysteresis is sufficient to predict treatment efficacy.
Surprisingly, the fast sequential protocols resulted in significantly fewer resistance types and thus less genotypic diversity than the slow regular treatments (Figs. 2B and 3B and SI Appendix, Fig. S4 and Table S3). These findings contrast with expectations from population genetic theory, because fast switching should have prevented competitive exclusion and, instead, caused coexistence of multiple types that continuously varied in relative frequency parallel to antibiotic exposure. To assess these dynamics, we developed and analyzed a mathematical model tailored to the design of the evolution experiment. Under standard conditions (without hysteresis), the model indeed predicted coexistence of several types under fast sequential conditions (Fig. 2D and SI Appendix, Fig. S6). Importantly, when we added hysteresis effects to the model, we found increased selection pressure and a reduction of diversity, especially for the fast and random sequential protocols. These observations suggest that the inducible physiological effects act as a strong selective constraint during sequential treatment and influence diversity within the evolving populations.
The mathematical model indicated that negative hysteresis increases selection intensity (Fig. 2D), yet the observed outcome was not MDR—as would be expected from competitive release (23)—but rather a constrained ability to evolve MDR (Fig. 2B). Thus, we hypothesized that negative hysteresis selects for traits specifically directed against the inducible physiological effects rather than resistance. This idea was supported by our additional analysis of growth rate. Almost all drug protocols resulted in significantly reduced growth rates under drug-free conditions (Fig. 2B, last column), but the three sequences (#5, #12, #15) with high levels of negative hysteresis and almost no evolved MDR showed the strongest growth reductions of up to 41%. The combination of reduced growth and no MDR is indicative of antibiotic tolerance (24), an evolutionary strategy, in which bacteria evade killing by slow growth and which could have been favored through selection by negative hysteresis.
Negative Hysteresis Favors Genetic Changes Mediating Tolerance and an Unknown Response.
To further assess the selective impact of negative hysteresis, we characterized the genes that have likely been the targets of selection using whole-genome sequencing and functional genetic analysis. The genomic characterization identified different sets of mutations to be favored by the main treatment types (Fig. 4 and Dataset S1). Intriguingly, the single nonresistant isolate from protocol #12 harbored a mutation that mediated a phenotypic response related to antibiotic tolerance. In detail, we followed concepts and methodology described previously (25), to test for tolerance with the help of time-kill experiments. These experiments revealed absence of resistance but reduced killing on all three antibiotics (Figs. 3B and 5B and SI Appendix, Table S7), consistent with antibiotic tolerance. This isolate had two mutations, a point mutation in the ispA gene (leading to amino acid change Y249D) and a frame shift in the gcvT2 gene. Because the gcvT2 mutation occurred across treatment groups (Fig. 4A), the ispA mutation is likely the adaptive change that caused reduced growth, possibly due to the toxic accumulation of isoprenyl diphosphates, as previously recorded for a ΔispA E. coli mutant (26). Sequence #12 was enriched for CAR-induced sensitization toward GEN (Fig. 2A). A reassessment of the CAR → GEN transition showed that bacterial cells of this isolate could no longer be sensitized (Fig. 5C). We conclude that, in this single case, selection by negative hysteresis in sequence #12 has likely been countered by the emergence of a process related to antibiotic tolerance, mediated through a mutation in ispA.
Fig. 4.
Fig. 5.
In several other fast sequential protocols, negative hysteresis was countered by mutations in cpxS (SI Appendix, Fig. S7). This gene is related to the E. coli envelope stress response system CpxA−CpxR, which is activated by misfolded proteins, as caused by aminoglycosides (27), and involved in intrinsic resistance to these drugs in E. coli (28). Mutations in cpxS were significantly enriched in fast sequential treatments (SI Appendix, Table S6), including those with little indication of antibiotic tolerance (e.g., normal growth under drug-free conditions; protocol #7; Fig. 2B). To explore its function, we reintroduced one prevalent cpxS mutation (leading to T163P) into the ancestral background and compared it to the ancestor and the similarly generated mexR (T130P) mutant as a control. MexR regulates the multidrug efflux pump MexAB-OprM, which can extrude different drug classes (29), potentially conferring MDR. MexR mutations are unlikely favored by negative hysteresis, as this gene was frequently and exclusively mutated in slow sequential treatments (Fig. 4A). Our analysis revealed that resistance against CIP and CAR was moderate for the cpxS mutant, but strongly increased for the mexR mutant, while neither mutation altered resistance to GEN (Fig. 5A and SI Appendix, Fig. S8). Importantly, using two independent methods, we consistently found that the CAR-induced sensitization toward GEN was abolished in the cpxS mutant, but still present in the mexR mutant (Fig. 5C), which is also consistent with growth rate dynamics in the evolution experiment (SI Appendix, Fig. S9). In conclusion, mutations in cpxS were favored in several fast sequential protocols and are thus likely to represent a general response to selection by negative hysteresis, apparently independent of antibiotic tolerance.
The above result for the mexR mutant additionally suggests that canonical resistance mechanisms do not abolish negative hysteresis. To further test this point, we characterized CAR → GEN hysteresis for defined mutants in mexA, nalC, and nalD, all affecting the MexAB-OprM multidrug efflux system. Despite their increased resistance to the antibiotic of the pretreatment (i.e., CAR), negative hysteresis and thus enhanced killing by GEN could be readily induced, with high (up to 600 µg/mL) and also subinhibitory doses (SI Appendix, Fig. S10). We conclude that bacteria can suffer from negative hysteresis even if they are resistant against the drug of the sensitizing pretreatment.
An Independent Experimental Test Validates the Importance of Negative Hysteresis.
To specifically test the consequences of negative hysteresis, we took inspiration from Lewontin’s theoretical work on evolutionary historicity (30), which highlighted that the order of events influences the evolutionary outcome. We thus repeated evolution experiments with the reversed order of drugs for the most effective sequence, #12, and the least effective sequence, #13 (inferred from evolved MDR; Fig. 2B). The reverse sequences had the same drug proportions and number of switches as the original sequences, but the direction of transitions was opposite. As a consequence, all was equal except that the cumulative level of negative hysteresis was decreased by 10% in the first case and increased by 11% in the second case (Fig. 6A). As expected, reversing #12 decreased extinction frequency (Fig. 6B) and significantly increased resistance (Fig. 6C and SI Appendix, Table S7). Conversely, reversing #13 increased extinction, although resistance was not affected, most likely because only few populations survived and could be used for resistance analysis. These results demonstrate that cellular hysteresis can determine the efficacy of sequential therapy.
Fig. 6.
Conclusions
Our results show that antibiotics can induce long-lasting changes in bacterial physiology that enhance or inhibit the bactericidal activity of other antibiotics, thereby revealing a principle, cellular hysteresis, that can be exploited to optimize antibiotic therapy. Cellular hysteresis can still act in bacteria that are resistant against the pretreatment drug, highlighting its clinical potential, where antibiotic resistance is widespread. Fast switching between antibiotics is key, because it can increase the cumulative effect of negative hysteresis, leading, in our experiments, to improvements in three complementary measures of treatment efficacy. Negative hysteresis exerted specific selection that did not favor resistance mutations, but rather mutations that counter the inducible physiological effects, such as those identified here for ispA and cpxS. We also showed that slow drug changes enhance resistance evolution and target different sets of genes. Our findings may thus explain the limited success of antibiotic cycling in the clinic, where antibiotics are usually changed once per month or less often (31). Negative hysteresis may further explain the particular success of one of the very few clinical studies with one antibiotic switch in less than a day. This study from 1988 demonstrated that the staggered application of first a β-lactam and then, 4 h later, an aminoglycoside—exactly the switch required for negative hysteresis—causes a significant reduction and often full clearance of P. aeruginosa from the lungs of a small cohort of cystic fibrosis patients (32). In this particular study, clearance could not be achieved by the corresponding combination treatment [i.e., simultaneous dosing (32)]. In our own previous work, the CIP+GEN combination did not produce any clearance, while the CAR+GEN combination at equivalent effective dose led to clearance rates comparable to those observed here for sequential treatments with high levels of negative hysteresis (19), generally supporting the possible power of temporal changes for effective treatment. Novel approaches such as multilayer liposomes (33) may then become useful for sequential drug delivery. A further exploration of inducible physiological effects may thus help to find new ways for improving antibiotic therapy—using the available drugs in a rational and refined way.
Materials and Methods
Material.
P. aeruginosa UCBPP-PA14 (34) was grown at 37 °C in M9 medium supplemented with glucose (2 g/L), citrate (0.5 g/L), casamino acids (1 g/L), and an antibiotic, as indicated.
Time-Kill Experiments.
Exponential phase cells (5 × 107 cfu/mL) were pretreated with sublethal antibiotic concentrations for the indicated time. The medium was exchanged to expose cells to a second antibiotic at IC75, followed by cfu counting for 6 h. Hysteresis effects were quantified as the average log10 difference in cfu counts of pretreated and control cultures. Negative values indicate sensitization, and positive values indicate protective effects. Bacterial killing was confirmed by time-lapse microscopy, monitoring bacteria on agarose pads (35) using a Zeiss LSM 700.
Dose–Response Curves.
A standardized inoculum (5 × 105 cfu) was incubated with defined antibiotic concentrations in 96-well plates for 12 h at 37 °C, followed by optical density measurements at 600 nm (OD600).
Main Evolution Experiment.
We performed serial dilution evolution experiments (6, 23, 36) and selection with 16 different antibiotic sequences (see also SI Appendix). Sequences #1 to #4 had constant environments. Sequences #5 to #16 contained equal frequencies of CIP, GEN, and CAR, but differed in hysteresis levels, due to drug order and switching rate. Each treatment had 12 replicates (founded with 5 × 105 cells from six independent PA14 starting cultures) and 96 serial transfers (2% transfer volume), each separated by 12 h. Antibiotic selection was applied at IC75 in 96-well plates, and growth was monitored by OD600 every 15 min (EON; BioTek Instruments; 180-rpm double orbital shaking). Evolved material was conserved at −80 °C in 10% (vol/vol) DMSO. Resistance evolution was assessed using the integral of the growth curve divided by the integral for the untreated reference evolving in parallel (relative area under curve or relative biomass, Fig. 3). Low values denote sensitivity to treatment, a value of 1 uninhibited growth (dynamics for all populations are shown in SI Appendix, Fig. S11). Adaptation rate was calculated with a sliding window approach as X−1, where X is the transfer at which the mean relative biomass of a sliding window of 12 transfers reaches 0.75 for the first time. This measure is comparable with the previously described rate of adaptation (36), defined for constant environments, yet not applicable to fluctuating environments, in which growth often oscillates. Extinction frequencies were determined at the end of the evolution experiment by counting cases unable to grow in drug-free media.
Characterization of Evolved Isolates.
We measured antibiotic dose–response curves for 880 evolved isolates from evolved populations after transfers 12 and 48. Resistance profiles were obtained as in Dose–Response Curves, with concentrations from 1/8 MIC to 16× MIC, assessed as the relative area of dose–response curves for isolates and corresponding ancestors measured on the same plate (SI Appendix, Fig. S3). For treatment comparisons, we defined MDR scores as the sum of resistance values on the three antibiotics. Subpopulations were identified by hierarchical clustering of resistance profiles. We characterized growth in drug-free medium. See SI Appendix.
Mathematical Model.
We developed a deterministic model to explore the ability of different antibiotic protocols to limit population growth by the evolution of resistant types. We modified a logistic growth model (competition for space) that included mutation, . Different genotypes were included with density xi and growth rate ri. Each genotype had three growth rates, for each of the possible treatments, ri = {riCIP, riGEN, riCAR}. The mutation rate qji determined the change of genotype j to another genotype i. The carrying capacity was defined by K. To simulate serial transfers, the mixture of types was diluted by a dilution factor DF at the end of each season. If the density of a genotype fell below the cutoff κ during dilution, it was lost and could only reappear via mutation. Following dilution, treatments could either switch or be repeated. The model was parameterized according to the evolution experiment: K = 108 cells, DF = 50 applied every 12 h, κ = 10. Population size was K/4 (IC75) directly before the first transfer. Growth dynamics were generated for a simple system with four competing genotypes, the nonresistant wt and three mutants, individually resistant to CIP, GEN, or CAR, parameterized according to the results of the monotreatments #1 to #3 (growth rate table R, SI Appendix, Fig. S6). Some mutant growth rates were lower than those of the wt on particular antibiotics, denoting collateral sensitivity, consistent with previous results (37). Switches between antibiotics allowed for hysteresis effects, which we included by multiplying the respective growth rates from table R with the corresponding entry from table S, showing the antibiotic-induced physiological effects, experimentally inferred for the four genotypes (SI Appendix, Fig. S6). Using this model, we generated growth dynamics for mixed populations for the different sequential treatments, either with or without hysteresis. From the modeled dynamics, we inferred the selective pressure, as defined by K/xwt, and the within-population diversity, as calculated from Shannon entropy.
Genomics and Functional Genetic Analysis.
Whole-genome sequencing was performed for 30 evolved isolates from different subpopulations at the early time point, and the three subpopulations of #8 at the late time point (SI Appendix). For defined mutants (SI Appendix), we assessed the change in growth upon pretreatment by OD600 every 15 min. We additionally used flow cytometry (Guava EasyCyte HT Blue-Green; Merck KGaA) with hourly samples, Live/Dead staining, and three technical replicates. For staining, cells were incubated for 10 min with propidium iodide (P4170-25MG; Sigma-Aldrich) and thiazole orange (390062-250MG; Sigma-Aldrich). We assessed antibiotic tolerance for isolate 12-1a-E2-4 (isolated after transfer 12 from sequence #12) via minimal duration of killing (25).
Replay Evolution Experiment.
The experiment was performed as the main evolution experiment, using sequences #12, #12rev, #13, #13rev, and an untreated control. Sequences #12 and #13 were the same as the first 12 transfers in the main experiment, and #12rev and #13rev were their respective reverse sequences. Resistance was quantified by the fold changes in IC75.
Statistical Analyses.
Data analysis was performed with R (38). Statistics, P values, and explanatory notes are provided in SI Appendix, Tables S1–S6.
Data Availability
The data are supplied as Datasets S1–S5. Sequence data are available from NCBI, BioProject PRJNA484297.
Data Availability
Data deposition: The sequences reported in this paper have been deposited in the NCBI Sequence Read Archive (SRA) database, https://www.ncbi.nlm.nih.gov/sra (Bioproject no. PRJNA484297).
Acknowledgments
We thank T. Bollenbach, T. Dagan, C. Eschenbrenner, D. Falush, M. Habig, A. Read, J. Rolff, M. Sixt, and the H.S. lab for advice; and G. Hemmrich-Stanisak, T. Naujoks, C. Noack, and M. Vollstedt from the Institute of Clinical Molecular Biology for DNA sequencing, as supported by the German Science Foundation Cluster of Excellence EXC 608. This research was funded by the German Science Foundation Grant SCHU 1415/12 (to H.S.), a Young Scientist Grant from the Zentrum Molekulare Biowissenschaften Kiel (to R.R.), the International Max Planck Research School for Evolutionary Biology (R.R.), the Swedish Research Council (D.I.A.), the Leibniz Science Campus EvoLUNG (H.S.), and the Max-Planck Society (C.S.G., A.T., and H.S.).
Supporting Information
Appendix (PDF)
- Download
- 1.59 MB
Dataset_S01 (XLSX)
- Download
- 17.55 KB
Dataset_S02 (TXT)
- Download
- 1.03 MB
Dataset_S03 (TXT)
- Download
- 226.97 KB
Dataset_S04 (TXT)
- Download
- 3.92 KB
Dataset_S05 (TXT)
- Download
- 1.63 KB
Movie S1.
Illustration of negative hysteresis by time-lapse microscopy. Time-lapse of a time-kill experiment over 6 h using pretreatment with 100 μg/ml carbenicillin followed by exposure to 2 μg/ml gentamicin. Differential interference contrast (DIC) and corresponding fluorescence images were taken every 10 min, and composite images were assembled for the video. Dead cells are shown in red due to propidium iodide staining.
- Download
- 957.73 KB
References
1
R Bürger, M Lynch, Evolution and extinction in a changing environment: A quantitative-genetic analysis. Evolution 49, 151–163 (1995).
2
SP Carroll, et al., Applying evolutionary biology to address global challenges. Science 346, 1245993 (2014).
3
Y-Y Liu, et al., Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: A microbiological and molecular biological study. Lancet Infect Dis 16, 161–168 (2016).
4
GV Bloemberg, et al., Acquired resistance to bedaquiline and delamanid in therapy for tuberculosis. N Engl J Med 373, 1986–1988, and erratum (2015) 375:e29 (2015).
5
S Kim, TD Lieberman, R Kishony, Alternating antibiotic treatments constrain evolutionary paths to multidrug resistance. Proc Natl Acad Sci USA 111, 14494–14499 (2014).
6
R Roemhild, C Barbosa, RE Beardmore, G Jansen, H Schulenburg, Temporal variation in antibiotic environments slows down resistance evolution in pathogenic Pseudomonas aeruginosa. Evol Appl 8, 945–955 (2015).
7
A Fuentes-Hernandez, et al., Using a sequential regimen to eliminate bacteria at sublethal antibiotic dosages. PLoS Biol 13, e1002104 (2015).
8
M Yoshida, et al., Time-programmable drug dosing allows the manipulation, suppression and reversal of antibiotic drug resistance in vitro. Nat Commun 8, 15589 (2017).
9
W Szybalski, V Bryson, Genetic studies on microbial cross resistance to toxic agents. I. Cross resistance of Escherichia coli to fifteen antibiotics. J Bacteriol 64, 489–499 (1952).
10
L Imamovic, MOA Sommer, Use of collateral sensitivity networks to design drug cycling protocols that avoid resistance development. Sci Transl Med 5–204ra132 (2013).
11
JF Linares, I Gustafsson, F Baquero, JL Martinez, Antibiotics as intermicrobial signaling agents instead of weapons. Proc Natl Acad Sci USA 103, 19484–19489 (2006).
12
LR Hoffman, et al., Aminoglycoside antibiotics induce bacterial biofilm formation. Nature 436, 1171–1175 (2005).
13
RA VanBogelen, FC Neidhardt, Ribosomes as sensors of heat and cold shock in Escherichia coli. Proc Natl Acad Sci USA 87, 5589–5593 (1990).
14
G Lambert, E Kussell, Memory and fitness optimization of bacteria under fluctuating environments. PLoS Genet 10, e1004556, and erratum (2014) 10:e1004556 (2014).
15
S Ram, M Goulian, The architecture of a prototypical bacterial signaling circuit enables a single point mutation to confer novel network properties. PLoS Genet 9, e1003706 (2013).
16
Z Lazim, TJ Humphrey, RJ Rowbury, Induction of the PhoE porin by NaCl as the basis for salt-induced acid sensitivity in Escherichia coli. Lett Appl Microbiol 23, 269–272 (1996).
17
J-W Veening, WK Smits, OP Kuipers, Bistability, epigenetics, and bet-hedging in bacteria. Annu Rev Microbiol 62, 193–210 (2008).
18
; World Health Organization, Global priority list of antibiotic-resistant bacteria to guide research, discovery, and development of new antibiotics. Available at www.who.int/en/news-room/detail/27-02-2017-who-publishes-list-of-bacteria-for-which-new-antibiotics-are-urgently-needed. Accessed May 16, 2017. (2017).
19
C Barbosa, R Beardmore, H Schulenburg, G Jansen, Antibiotic combination efficacy (ACE) networks for a Pseudomonas aeruginosa model. PLoS Biol 16, e2004356 (2018).
20
PH Plotz, BD Davis, Synergism between streptomycin and penicillin: A proposed mechanism. Science 135, 1067–1068 (1962).
21
T Hay, S Fraud, CH-F Lau, C Gilmour, K Poole, Antibiotic inducibility of the mexXY multidrug efflux operon of Pseudomonas aeruginosa: Involvement of the MexZ anti-repressor ArmZ. PLoS One 8, e56858 (2013).
22
Y Morita, ML Sobel, K Poole, Antibiotic inducibility of the MexXY multidrug efflux system of Pseudomonas aeruginosa: Involvement of the antibiotic-inducible PA5471 gene product. J Bacteriol 188, 1847–1855 (2006).
23
R Pena-Miller, et al., When the most potent combination of antibiotics selects for the greatest bacterial load: The smile-frown transition. PLoS Biol 11, e1001540 (2013).
24
O Fridman, A Goldberg, I Ronin, N Shoresh, NQ Balaban, Optimization of lag time underlies antibiotic tolerance in evolved bacterial populations. Nature 513, 418–421 (2014).
25
A Brauner, O Fridman, O Gefen, NQ Balaban, Distinguishing between resistance, tolerance and persistence to antibiotic treatment. Nat Rev Microbiol 14, 320–330 (2016).
26
S Fujisaki, et al., Disruption of the structural gene for farnesyl diphosphate synthase in Escherichia coli. J Biochem 137, 395–400 (2005).
27
BD Davis, LL Chen, PC Tai, Misread protein creates membrane channels: An essential step in the bactericidal action of aminoglycosides. Proc Natl Acad Sci USA 83, 6164–6168 (1986).
28
RL Guest, TL Raivio, Role of the Gram-negative envelope stress response in the presence of antimicrobial agents. Trends Microbiol 24, 377–390 (2016).
29
N Masuda, et al., Substrate specificities of MexAB-OprM, MexCD-OprJ, and MexXY-oprM efflux pumps in Pseudomonas aeruginosa. Antimicrob Agents Chemother 44, 3322–3327 (2000).
30
RC Lewontin, The principle of historicity in evolution. Wistar Inst Symp Monogr 5, 81–94 (1967).
31
P Abel zur Wiesch, R Kouyos, S Abel, W Viechtbauer, S Bonhoeffer, Cycling empirical antibiotic therapy in hospitals: Meta-analysis and models. PLoS Pathog 10, e1004225 (2014).
32
JP Guggenbichler, F Allerberger, MP Dierich, R Schmitzberger, E Semenitz, Spaced administration of antibiotic combinations to eliminate Pseudomonas from sputum in cystic fibrosis. Lancet 2, 749–750 (1988).
33
SW Morton, et al., A nanoparticle-based combination chemotherapy delivery system for enhanced tumor killing by dynamic rewiring of signaling pathways. Sci Signal 7, ra44 (2014).
34
LG Rahme, et al., Common virulence factors for bacterial pathogenicity in plants and animals. Science 268, 1899–1902 (1995).
35
JW Young, et al., Measuring single-cell gene expression dynamics in bacteria using fluorescence time-lapse microscopy. Nat Protoc 7, 80–88 (2011).
36
M Hegreness, N Shoresh, D Damian, D Hartl, R Kishony, Accelerated evolution of resistance in multidrug environments. Proc Natl Acad Sci USA 105, 13977–13981 (2008).
37
C Barbosa, et al., Alternative evolutionary paths to bacterial antibiotic resistance cause distinct collateral effects. Mol Biol Evol 34, 2229–2244 (2017).
38
; R Core Team, R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna), R Version 3.3.1. Available at www.R-project.org/. Accessed June 21, 2016. (2016).
Information & Authors
Information
Published in
Classifications
Copyright
Copyright © 2018 the Author(s). Published by PNAS. This open access article is distributed under Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND).
Data Availability
Data deposition: The sequences reported in this paper have been deposited in the NCBI Sequence Read Archive (SRA) database, https://www.ncbi.nlm.nih.gov/sra (Bioproject no. PRJNA484297).
Submission history
Published online: September 12, 2018
Published in issue: September 25, 2018
Keywords
Acknowledgments
We thank T. Bollenbach, T. Dagan, C. Eschenbrenner, D. Falush, M. Habig, A. Read, J. Rolff, M. Sixt, and the H.S. lab for advice; and G. Hemmrich-Stanisak, T. Naujoks, C. Noack, and M. Vollstedt from the Institute of Clinical Molecular Biology for DNA sequencing, as supported by the German Science Foundation Cluster of Excellence EXC 608. This research was funded by the German Science Foundation Grant SCHU 1415/12 (to H.S.), a Young Scientist Grant from the Zentrum Molekulare Biowissenschaften Kiel (to R.R.), the International Max Planck Research School for Evolutionary Biology (R.R.), the Swedish Research Council (D.I.A.), the Leibniz Science Campus EvoLUNG (H.S.), and the Max-Planck Society (C.S.G., A.T., and H.S.).
Notes
This article is a PNAS Direct Submission. S.B. is a guest editor invited by the Editorial Board.
Authors
Competing Interests
The authors declare no conflict of interest.
Metrics & Citations
Metrics
Citation statements
Altmetrics
Citations
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.
Cited by
Loading...
View Options
View options
PDF format
Download this article as a PDF file
DOWNLOAD PDFGet Access
Login options
Check if you have access through your login credentials or your institution to get full access on this article.
Personal login Institutional LoginRecommend to a librarian
Recommend PNAS to a LibrarianPurchase options
Purchase this article to access the full text.