
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39261449
52166
10.1038/s41467-024-52166-z
Article
The evolution of heteroresistance via small colony variants in Escherichia coli following long term exposure to bacteriostatic antibiotics
Gil-Gil Teresa 1
http://orcid.org/0000-0003-3430-0901
Berryhill Brandon A. 12
Manuel Joshua A. 1
Smith Andrew P. 1
McCall Ingrid C. 1
Baquero Fernando 3
http://orcid.org/0000-0003-3600-2780
Levin Bruce R. blevin@emory.edu

1
1 https://ror.org/03czfpz43 grid.189967.8 0000 0004 1936 7398 Department of Biology, Emory University, Atlanta, GA 30322 USA
2 https://ror.org/03czfpz43 grid.189967.8 0000 0004 1936 7398 Program in Microbiology and Molecular Genetics, Graduate Division of Biological and Biomedical Sciences, Laney Graduate School, Emory University, Atlanta, GA 30322 USA
3 grid.466571.7 0000 0004 1756 6246 Hospital Universitario Ramón y Cajal, Instituto Ramón y Cajal de Investigación Sanitaria, and Centro de Investigación Médica en Red, Epidemiología y Salud Pública (CIBERESP), Madrid, Spain
11 9 2024
11 9 2024
2024
15 793630 10 2023
27 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Traditionally, bacteriostatic antibiotics are agents able to arrest bacterial growth. Despite being traditionally viewed as unable to kill bacterial cells, when they are used clinically the outcome of these drugs is frequently as effective as when a bactericidal drug is used. We explore the dynamics of Escherichia coli after exposure to two ribosome-targeting bacteriostatic antibiotics, chloramphenicol and azithromycin, for thirty days. The results of our experiments provide evidence that bacteria exposed to these drugs replicate, evolve, and generate a sub-population of small colony variants (SCVs) which are resistant to multiple drugs. These SCVs contribute to the evolution of heteroresistance and rapidly revert to a susceptible state once the antibiotic is removed. Stated another way, exposure to bacteriostatic drugs selects for the evolution of heteroresistance in populations previously lacking this trait. More generally, our results question the definition of bacteriostasis as populations exposed to bacteriostatic drugs are replicating despite the lack of net growth.

Traditionally, bacteriostatic antibiotics were thought to simply arrest bacterial growth. In this work, the authors explore the dynamics of Escherichia coli when exposed to bacteriostatic antibiotics, finding that long-term exposure to bacteriostatic antibiotics favors the evolution of heteroresistant small colony variants.

Subject terms

Antibiotics
Experimental evolution
Antimicrobial resistance
https://doi.org/10.13039/100000057 U.S. Department of Health & Human Services | NIH | National Institute of General Medical Sciences (NIGMS) R35GM136407 Levin Bruce R. https://doi.org/10.13039/100000060 U.S. Department of Health & Human Services | NIH | National Institute of Allergy and Infectious Diseases (NIAID) U19AI158080 Levin Bruce R. issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Antibiotics can be broadly classified as being bactericidal or bacteriostatic based on whether they kill bacteria or simply arrest their growth1. However, it is known that at high drug concentrations bacteriostatic drugs can and do kill bacteria. Intuitively, it would make sense to treat an infection with drugs that readily kill the infecting bacteria, the bactericidal drugs, and thereby eliminate the reliance on the host’s immune system to clear the infection, as would be the case with bacteriostatic antibiotics. For this reason, bacteriostatic drugs have been considered “weaker” than bactericidal drugs and are not recommended for the treatment of severe infections or infections in immunodeficient patients2,3.

This distinction between bactericidal and bacteriostatic drugs is manifest in quantitative experimental studies of the pharmacodynamics (PD) of antibiotics and bacteria. These studies focus on the rates and dynamics of the drug’s ability to kill exposed populations of bacteria4–7. Many of the studies concerning why antibiotics fail to control bacterial infections have focused on phenomena solely belonging to bactericidal antibiotics such as persistence and tolerance, that is the ability to survive exposure to high concentrations of an antibiotic of a subpopulation or the whole bacterial population, respectively6,8–10. These traits are qualities only associated with bactericidal drugs: the bacteriostatic drugs do not kill, and consequently the minority surviving or slower-dying populations are not revealed. With much of the clinical application of antibiotics focusing on bactericidal drugs and the majority of the research on the PD of antibiotics and bacteria also focusing on bactericidal drugs, PD research on bacteriostatic drugs has been relatively neglected.

However, in recent years, clinicians have given less importance to the antibiotic’s ability to kill bacteria in vitro and instead have focused on the outcome of treatment with these drugs. By this criterion, in many cases bacteriostatic antibiotics are as effective as bactericidal even in severe infections, with the possible exception of immunosuppressed patients2,11. The increase in use of bacteriostatic antibiotics could reduce the selection pressure for resistance to bactericidal agents used in these critical infections. However, the shift to using bacteriostatic antibiotics requires the development of quantitative measures of the PD of antibiotics that arrest the growth of, rather than kill, bacteria12. One cannot solely characterize the PD of bacteriostatic antibiotics by the minimum concentration required to prevent the replication of exposed bacteria, the MIC13. Furthermore, the evolution of genomic resistance for these agents is rare and tolerance or persistence that occurs in bactericidal antibiotics are difficult to define and perhaps impossible to detect with these bacteriostatic agents14,15. This raises questions about the population and evolutionary dynamics of bacteria confronted with these drugs. If the bacteria exposed to bacteriostatic antibiotics are not replicating, one would not expect them to evolve. Therefore, more considerations of the pharmaco-, population, and evolutionary dynamics of bacteriostatic antibiotics are needed both academically and clinically.

In this study, we explore the pharmaco-, population, and evolutionary dynamics of Escherichia coli exposed to two bacteriostatic antibiotics of different classes, chloramphenicol (CHL) and azithromycin (AZM), over 30 days. While resistance did evolve over the 30 days, we did not detect any previously described genotypic changes that could account for antibiotic resistance, such as antibiotic efflux, decrease in membrane permeability, or antibiotic inactivation for CHL16; nor the overproduction of efflux pumps or mutations in genes encoding the 23S rRNA subunit for AZM17. The results of our experiments provide evidence that: (i) long term exposure to these ribosome-targeting bacteriostatic antibiotics does not change the net density of exposed populations, (ii) despite the fact that the population’s net density does not change, bacteria exposed to these drugs replicate, evolve, and generate small colony variants, and (iii) the selective pressure mediated by these drugs favors the evolution of heteroresistant populations, i.e. the emergence of resistant minority populations, in a strain previously lacking this trait.

Results

Long term exposure to bacteriostatic antibiotics

We begin our investigation into the effects of long term exposure to ribosome-targeting bacteriostatic drugs by evaluating the impact that these agents have on bacterial survival over 30 days. We exposed four independent cultures of ~105 CFU/mL of E. coli MG1655 in glucose-limited minimal media to super-MIC (Minimum Inhibitory Concentration) concentrations of CHL and AZM for 30 days without transferring or diluting (Fig. 1). One thousand μg/mL of glucose was used since this concentration allows cultures to grow to maximum stationary densities while avoiding nutrient limitation. The MIC for CHL and AZM of MG1655 were estimated by broth microdution in this glucose-limited minimal media and found to be 6.25 ug/mL for both drugs13. Super-MIC concentrations of each drug that were shown to be bacteriostatic with minimal killing and/or growth were used in the experiments (Supplementary Fig. 1). Cultures exposed to these drugs were sampled every five days. We observed that over the course of the experiment, the control cultures containing no drugs reached their maximum stationary phase density of approximately 109 CFU/mL and then went down by approximately two logs over the course of 30 days, while the densities in the cultures containing the drugs remained stable, with at most a half-log change in density in the drug-treated populations. Notably, in drug-treated cultures where the bacteria were not lost, two distinct colony morphologies emerged. The colonies were either similar to the ancestral wild-type E. coli or were much smaller bacterial colonies, small colony variants (SCVs). This evolution occurred while under strong selective pressure from these drugs. There was no change in colony size in the drug-free controls, with a limit of detection of 10−3 SCVs per colony of the ancestral population. These SCVs started to emerge at day 10 of the long term experiment in cultures exposed to CHL and at day 15 in cultures exposed to AZM. Their frequency increased from 2–5% when they were first observed in the presence of both antibiotics to 10–44% in CHL and 10–50% in AZM by day 30. To assess the maintenance of the activity of the antibiotics after 30 days, bacteria resistant to the antibiotic in each culture was added at approximately 106 CFU/mL and over the course of 24 h each culture grew approximately three orders of magnitude; while when the susceptible, wild-type E. coli strain was added to these 30-day supernatants it was not able to increase in density (Supplementary Fig. 2). This residual growth indicates that at 30 days the antibiotic is still at a super-MIC concentration and thereby is the limiting factor for growth in the long term experiments. Moreover, 0- and 30-day supernatants were filtered, and the amount of antibiotic present tested through a disk diffusion assay. The diameter of both time points was 1 cm for the CHL and 2.8 cm for the AZM supernatants confirming the stability of both antibiotics during the long term experiments.Fig. 1 Long term exposure of E. coli to bacteriostatic drugs.

Density in CFU/mL of E. coli MG1655 measured every five days for 30 days of four independent biological replicas (I-IV). A E. coli exposed to 4x MIC CHL (chloramphenicol); B E. coli exposed to 3x MIC AZM (azithromycin); C Drug-free control.

To determine how evolution is occurring in the apparent absence of net growth in the presence of the antibiotic, we performed a long term experiment using a conditionally non-replicative plasmid 18 to identify if growth is occurring at the same rate as death (Fig. 2). After 10 days the plasmid frequency decreased by 100-fold, and after 20 days the plasmid frequency was near the limit of detection. That means that the plasmid containing cells are progressively diluted (as each cell division gives rise to a plasmid-free descendant) and indicates that the population is growing, replicating at least once a day and dying at the same rate. To explain how these results are consistent with the hypothesis that growth is occurring at the same rate as death, we constructed a mathematical and computer-simulation model of the plasmid being lost both with and without bacteriostatic drugs (Supplementary Text, Supplementary Equations 1–4, and Supplementary Fig. 3).Fig. 2 Long term experiment with a non-replicative plasmid.

Ratio of the plasmid-containing cells to the total number of cells in CFU/mL of E. coli MG1655 with the non-replicative plasmid pAM34 which bears an ampicillin resistance cassette. The total cell density and density of cells bearing the plasmid were measured every 10 days for 30 days of four independent biological replicas for four conditions: i) minimal media with CHL (chloramphenicol), ii) minimal media with AZM (azithromycin), iii) minimal media with ampicillin– maximal plasmid retention, iv) minimal media with no antibiotic– maximal plasmid loss. In the absence of CHL or AZM the initial cultures can only grow two logs until arriving to a maximum stationary phase, explaining the maximum loss of two log in the absence of ampicillin.

Small colony variants characterization

To determine what these SCVs are, we isolated six independently generated SCVs of MG1655 (CHL A, CHL B, CHL C, AZM A, AZM B, and AZM C) and characterized them phenotypically and genotypically. Firstly, to determine if the SCVs are a form of resistance that has emerged over the long term experiment, we estimated the MIC to the drugs they were previously exposed to (Fig. 3). Each SCV showed at least an 8-fold change over the ancestral MG1655’s MIC. On the other hand, normal-sized colonies from each culture were isolated and changes in MIC respecting the ancestral MG1655 were not observed in any of them. Each SCV has a distinct antibiotic susceptibility profile in terms of collateral susceptibilities and cross-resistances (Supplementary Fig. 4).Fig. 3 MIC of the SCVs to their respective drugs.

Three SCVs were isolated from each condition from day 30 of the long term experiment, grown up in 1.5x E. coli MG1655’s MIC for their respective drug, and then E-tested. Shown is the fold change respecting the MIC of the final evolved control populations, whose MIC was 6.25 μg/mL as it was at the start of the long term experiment.

Notably, although these SCVs are resistant to the bacteriostatic antibiotic to which they were exposed, there is very little net growth in the long term culture over the course of 30 days. To determine why a marked increase in density does not occur despite the evolution of resistance, we performed OD-based growth experiments of the SCVs with different concentrations of antibiotics. We found that even though these mutants are resistant, their growth rates and the maximum optical densities decreased proportionally to the drug concentration and their lag time was substantially increased (Supplementary Fig. 5). This result is consistent with previous observations13.

Antibiotic resistance in the small colonies obtained from both CHL and AZM cultures was unstable, that is, when streaked on LB plates without the drug the MIC of the SCVs decreased. After the genotypic characterization of these SCVs we do find genetic differences in most of them (Supplementary Table 1), but we could not find a clear mechanism that would explain this resistance. While we could not identify causal genetic differences between the SCVs and their ancestor MG1655, since both CHL and AZM prevent translation, we sought to determine if there were differences in the rate or total level of translation. Testing the level of a constitutively expressed enzyme, beta-galactosidase, we found that the six SCVs had higher levels of translation than the background MG1655 (Fig. 4). Moreover, we found that the rate of translation decreases in a concentration-dependent manner to both CHL and AZM (Supplementary Fig. 6). This is consistent with recent observations indicating partially compromised translation as a result of exposure to tetracycline, a translational inhibitor, which potentially accounts for the collateral susceptibility to this drug observed in Supplementary Fig. 419. Therefore, at the drug concentrations used in the long term experiments, these SCVs are readily selected for over the ancestor, not only for their level of antibiotic resistance but also for their level of translation, explaining their ability to increase in frequency when rare.Fig. 4 An assay of beta-galactosidase as a proxy for translation level.

The beta-galactosidase levels as determined by absorbance for the ancestral E. coli MG1655 and the six SCVs. Shown are the means and standard deviations of three independent biological replicates.

Heteroresistance

Antibiotic heteroresistance (HR) is defined as, “a phenotype in which a bacterial isolate contains subpopulations of cells that show a substantial reduction in antibiotic susceptibility compared with the main population”, and is detected via a population analysis profile (PAP) test20. The revertant populations obtained from these SCVs have a lower MIC than those of their small colony ancestor. These revertant populations are capable of rapidly regenerating the SCVs, which have a higher MIC – meeting the definition of antibiotic HR. The ancestral E. coli MG1655 is not capable of generating resistant subpopulations, as shown via PAP test (Supplementary Fig. 7 Panels A and B). In Fig. 5, we show PAP tests of CHL B (Panel A) and AZM C (Panel B). Both SCVs are shown to be heteroresistant based on the above criteria. Moreover, highly resistant colonies isolated from these PAP tests were found to be unstable and revert to the initial SCV state in 15 days, as shown by a decrease in the MIC (Fig. 3) as well as a decrease in the frequency of antibiotic-resistant subpopulations (Fig. 5). As shown in Supplementary Fig. 7 Panels C and D, all SCVs obtained from the long term experiment meet the criteria for HR. Pursuing the genetic basis of this reversion, we performed a genotypic characterization of these revertants. We did find extra mutations in CHL B and AZM B (Supplementary Table 2), but once again we did not find mutations that can explain reversion.Fig. 5 PAP tests of a CHL and an AZM SCV.

A PAP tests of a chloramphenicol SCV (blue line), the most resistant isolate of this clone (green line), and the most resistant isolate after being grown without antibiotic pressure for 15 days (orange line). B PAP tests of an azithromycin SCV (blue line), the most resistant isolate of this clone (green line), and the most resistant isolate after being grown without antibiotic pressure for 15 days (orange line).

Mathematical model and computer simulations

To explore the dynamics of the emergence of the SCVs and subsequently the evolution of HR described in the above experimental results, we constructed a mathematical and computer-simulation model. In particular, we are interested in explaining the timeframe and conditions in which the SCVs arise, why the SCVs fail to ascend to a nutrient-limited density, and why the evolution of HR is stable. In Fig. 6 we depict the model graphically and in the Supplementary Text and Supplementary Equations 5–12 we describe the model and its equations.Fig. 6 Diagram of a semi-stochastic model of the evolution of HR.

The variables N, S, and H are, respectively, the wildtype E. coli MG1655, SCVs, and the heteroresistant bacteria in cells/mL. The parameters µns, µsn, µsh and µhs are the transition rates, per cell per hour, between the different states.

In Fig. 7A we illustrate how the presence of an antibiotic selects for the emergence and ascent of resistant SCVs and a heteroresistant population from an initially wild-type population. In Fig. 7B, we show that an initial population of SCVs when grown without antibiotics will rapidly transition and give rise to a heteroresistant population which ascends and becomes limited by the resource. In Fig. 7C, we show that exposing the heteroresistant population to the drug more rapidly selects for the emergence and dominance of resistant SCVs than when exposing the wild type population to the same concentraton of the drug. In Fig. 7D, we show the changes in average MIC for the three scenarios depicted in panels A, B, and C. In Supplementary Fig. 8 we show the predicted dynamics of the above model with differing concentrations of the treating antibiotic. Notably, the emergence of the SCVs and HR is dependent on this concentration—critically, the drug concentration cannot exceed the MIC of the most resistant population. For a more detailed consideration of this result and the impact that having multiple resistant states has on the emergence of HR, see ref. 21. Additionally, in Supplementary Fig. 9, we explore the sensitivity of our model to the governing parameters: the transitions rates and relative fitness costs. Variations in the transition rates from the wild-type to the SCVs (Supplementary Fig. 9A, B) simply alter the rate at which SCVs emerge and dominate, without changing the overall outcome. Similarly, adjusting the transition rates from SCVs to the heteroresistant population (Supplementary Fig. 9C, D) only affects the frequency and timing of SCV dominance. Lastly, differing fitness costs of the SCVs (Supplementary Fig. 9E, F) merely impact the time required for SCV dominance. Together, these mathematical simulations illustrate that our observations hold true across a wide range of biologically meaningful parameter values.Fig. 7 Computer simulations of the proposed model of evolved HR.

Parameters used for these simulations are e = 5 × 10−7 μg/cell; vmaxN = 1.0, vmaxS = 0.5, vmaxH = 1.0 per cell per hour; vmin = −0.01 per cell per hour; K = 1; k = 1; MICN = 1.0, MICS = 8.0, MICH = 2.0 µg/ml; da = 0 µg/hour; µns = 1×10-8, µsn = 1 × 10-8, µsh = 1 × 10−3, µhs = 1 × 10−3 per cell per hour. The variables N, S, and H are, respectively, the wildtype E. coli MG1655, SCVs, and the heteroresistant bacteria in cells/mL. Note, points which do not appear to have error bars have standard deviations which are too small to be graphed. One hundred independent runs of the models were carried out in each condition. A H and S selection from an initial sensitive population. Parameters are N = 1 × 106 cells/mL; A = 3 μg/mL. B H selection from an initial SCV population in the absence of the drug. Parameters are S = 1 × 106 cells/mL; A = 0 μg/mL. C S selection from an initial H population in the presence of the drug. H = 1 × 106 cells/mL; A = 3 μg/mL. D Changes in the average MIC of the system over the simulations shown (A, B, and C).

Discussion

The canonical distinction between bacteriostatic and bactericidal antibiotics has deeply influenced their clinical usage. Traditionally, bacteriostatic antibiotics were considered “weaker drugs”, but this traditional view is questionable. Drugs which are classified as bacteriostatic can and do kill bacteria in a concentration dependent manner22. Moreover, meta-analysis studies do not demonstrate differences in the clinical success of therapy with either type of drugs even in severe infections11. Indeed, bactericidal agents could be reserved for life-threatening infections, particularly in immunocompromised patients or those suffering from chronic infections. A more extended use of bacteriostatic drugs could be beneficial to spare the use and overuse of bactericidal antibiotics which fosters resistance. A limitation to progress in the extended use of bacteriostatic drugs is the shortage of pharmacodynamic (PD) data with these drugs.

The distinction between bacteriostatic and bactericidal antibiotics is confounded by another factor: The primary cellular and molecular targets do not necessarily differ between these classes of drugs. Several bacteriostatic antibiotics have mechanisms of action that one would anticipate being bactericidal: mecillinam and cefsulodin inhibit cell wall synthesis23; novobiocin inhibits DNA gyrase24; and rifampin inhibits the DNA-dependent RNA polymerase which is bacteriostatic for E. coli and bactericidal for Mycobacterium25,26. Most notably, drugs which target the ribosome and inhibit protein synthesis can be either bactericidal (such as gentamicin or tobramycin) or bacteriostatic (such as the macrolides, phenicols, tetracyclines, oxazolidinones, and spectinomycin)27. Interestingly, the potential bactericidal effect of one of the ribosome-targeting drugs, chloramphenicol, is prevented by the production of (p)ppGpp in the exposed cells28. This suggests that the difference between a bacteriostatic drug and a bactericidal one is a property of the treated cell rather than the antibiotic29. All the above examples illustrate the need for a better understanding of the PD, population biology, and evolutionary biology of treatment with bacteriostatic antibiotics.

In this study, we present evidence that bacteriostatic drugs of two different classes (the phenicols and macrolides) inhibit the growth of E. coli for extended periods, i.e., 30 days, and moreover maintain the culture in a kind of stationary phase where the density of viable bacterial cells is stable. Although, unlike stationary phase we found there to be an abundance of the limiting resource, implying that the cultures remain drug-limited even after a month. Most interestingly, despite the fact that the bacteria in the population appear to not be replicating due to the lack of net growth, evolution still occurred. A population of small colony variants (SCVs) emerged and ascended to become a significant population of bacteria. We attribute this evolution to the fact that even though the population at large was neither increasing nor decreasing in the cultures with bacteriostatic drugs, the population was replicating at a rate roughly equal to that at which it was dying. This finding is unanticipated and inconsistent with the common perception that bacteriostatic drugs simply arrest bacterial growth. This result questions the definition of bacteriostasis.

Our mathematical and computer-simulation model demonstrates that in the presence of an antibiotic, even when starting with an antibiotic susceptible population that generates resistant sub-populations at an extremely low rate, the resistant sub-population will ascend to dominate. When the selection pressure is removed this resistant sub-population will increase in susceptibility, however the population which then emerges to dominance is not the initial antibiotic susceptible population but is instead a population which readily generates resistant sub-populations. Stated another way, our model illustrates and our results confirm how HR can evolve in a population otherwise lacking this trait.

These SCVs were found to be highly resistant not only to the challenging agents, but even to some types of bactericidal agents, such as aminoglycosides and rifampin. Curiously, the resistance of the SCVs to chloramphenicol (CHL) and azithromycin (AZM) are not due to the canonical resistance mechanisms16,30. SCVs have been implicated in treatment failure, primarily in Staphylococcus aureus, but there are limited reports of SCVs being associated with treatment failure in E. coli31,32. Although SCVs of E. coli have been observed for a long time, they are under-described both in their genetic basis and their role in bacterial communities. Two clinical isolates of E. coli SCVs have been studied, both of which have distinct metabolic disruptions by mutations in heme and thymidine biosynthesis pathways31,32. Lab-generated SCVs have been described with mutations in Q8, hemin, and lipoic acid biosynthesis33–35. However, none of these pathways appear to be disrupted in our SCVs.

We have yet to determine the genetic and molecular basis of the observed SCVs, but they appear to be distinct from the previously described mechanisms33–37. Certain mutations observed here (Supplementary Table 1) might account for changes in MICs in the SCVs. For instance, in the case of AZM-induced SCVs, missense variants of the citG gene, encoding the 2-(5”-triphosphoribosyl)-3’-dephosphocoenzyme-A synthase, were consistently found. This mutation might alter members of the GntR family of transcriptional regulators which could influence the DNA-binding properties of the regulator resulting in repression or activation of transcription, or could directly impact ATP synthesis38, leading to the generation of the SCVs. We did find extra mutations in CHL C, AZM A, and AZM C involving the AcrAB efflux pump, 50S ribosomal proteins L4 and L22, and the Lon protease, all of which can modify the antibiotic susceptibility. The mutations in AZM C found in acrAB could interfere with AZM membrane transport39. The mutations found in CHL C and AZM A in the gene encoding the 50S ribosomal protein L4 (a known chloramphenicol binding site) and L22, respectively, are a known mechanism of macrolide resistance40,41. The mutation in the Lon protease in AZM A is frequently associated with changes in antibiotic susceptibility as well42. Taken together, these changes could account for the increase in the MIC to AZM and CHL but cannot account for the emergence of the SCVs. We found a mutation in a CHL-induced SCV in the gene encoding the Tyrosine-tRNA ligase. This mutation could account for the generation of SCVs, since inhibitors of this ligase strongly decrease bacterial growth, but this mutation was only found in one of the six SCVs43. Most interestingly, we were unable to find any SNPs, duplications, or INDELs in one of the isolated CHL SCVs. However, all mutations found in the antibiotic-resistant SCVs were also found in the antibiotic-suspectable revertants (Supplementary Table 2), thus even if the mutations presented could account for the transition to resistance, they cannot account for the transition to susceptibility. This is also true of the transition to the SCV phenotype. We cannot exclude that the transition from normal cells to SCVs could be due to a kind of “structural epistasis” (not a genetic one) due to the altered interactions between intracellular machinery and molecules when the shape of the cell is altered by antibiotics44.

Unexpectedly, the antibiotic resistance observed here is transient, as would be anticipated for heteroresistance (HR), suggesting a high fitness cost of the mutations detected. In support of this HR hypothesis we found that these SCVs meet all the criteria set forth for HR: there are subpopulations present at a frequency greater than 10-7, with an MIC higher than 8x that of the main population, and reversion of the resistant subpopulation occurs in short order21. To our knowledge, this is the first report of both the spontaneous evolution of HR as well as HR to bacteriostatic drugs.

Methods

Bacterial strains

E. coli MG1655 was obtained from the Levin Lab bacterial collection13. pAM34 with the origin of replication under the control of an IPTG promoter and an ampicillin resistance cassette was obtained from Calin Guet from IST Austria18.

Growth media

LB (Lysogeny Broth) (244620) was obtained from BD. The DM (Davis Minimal) minimal base without dextrose (15758-500G-F) was obtained from Sigma Aldrich. MHII plates were made from MHII broth (212322) obtained from BD. Glucose (41095–5000) was obtained from Acros. LB agar (244510) for plates was obtained from BD.

Growth conditions

Unless otherwise stated, all experiments were conducted at 37 °C.

Sampling bacterial densities

The densities of bacteria were estimated by serial dilution in 0.85% saline and plating. The total density of bacteria was estimated on LB agar plates.

Antibiotics

Chloramphenicol (23660) was obtained from United States Pharmacopeia. Azithromycin (3771) was obtained from TOCRIS. Ampicillin (A9518-25G) was obtained from Sigma Aldrich. Isopropyl β-D-1-thiogalactopyranoside (IPTG; I56000-5.0) was obtained from Research Products International. All E-tests were obtained from bioMérieux.

Estimating minimal inhibitory concentrations

Antibiotic MICs were estimated using E-tests on MHII plates or via broth microdilution45,46.

Long term experiments

Flasks containing 10 mL of DM with 1000 μg/mL of glucose and an initial density of 105 CFU/mL cells were grown at 37 °C with shaking for 30 days without transferring. For each biological replicate, initial inoculums were obtained from independent overnights of E. coli cultures. E. coli MG1655 was grown with either no drug, 4x MIC of CHL or 3x MIC of AZM. Samples were taken every five days and plated on LB agar plates. Twenty SCVs from each condition were isolated from the last time point and grown in the absence of antibiotic. Only those whose MIC was higher than the ancestral MG1655 after two overnight passages in the absence of the drug were selected.

Long term experiments with non-replicative plasmid

Flasks containing 10 mL of DM with 1000 μg/mL of glucose and an initial density of 105 CFU/mL cells were grown at 37 °C with shaking for 30 days. E. coli MG1655 pAM34 was grown with either 4x MIC of CHL, 3x MIC of AZM, ampicillin, or no selection. All experiments were performed in the absence of IPTG. pAM34 is a single copy number plasmid. Samples were taken every five days and plated on LB agar plates as well as 100 μg/mL ampicillin and 0.5 mM IPTG LB agar plates.

Sequencing

Complete genomes were obtained with hybrid Illumina/Nanopore sequencing by SeqCenter. Samples were extracted from single colonies using Zymo Quick-DNA HMW MagBead Kit. Oxford Nanopore Sequencing library prep was performed with PCR-free ligation library prep using ONT’s V14 chemistry. Long read sequencing was performd using R10.4.1 flowcells on a GridION with basecalling performed by Guppy in Super High Accuracy mode. Illumina libraries were prepared and sequenced per SeqCenter’s standards. Quality control and adapter trimming was performed with bcl-convert and porechop for Illumina and ONT sequencing respectively. Hybrid assembly with Illumina and ONT reads was performed with Unicycler47. Assembly statistics were recorded with QUAST48. Assembly annotation was performed with Prokka49.

Beta-galactosidase assay

The BetaRedTM β-Galactosidase Assay Kit (Sigma Aldrich) was used to experimentally measure the translation level. Flasks containing 10 mL of DM with 1000 μg/mL of glucose and an initial density of 105 CFU/mL cells were grown at 37 °C with shaking overnight in the absence or presence of different concentrations of CHL or AZM. The manufacturer protocol was followed, and OD 570 nm was measured. The experiments were performed in triplicate.

Growth rate estimation

Exponential growth rates were estimated from changes in optical density (OD600) in a Bioscreen C. For this, 24 h stationary phase cultures were diluted in LB or glucose-limited liquid media to an initial density of approximately 105cells per mL. Five replicas were made for each estimate by adding 300 µL of the suspensions to the wells of the Bioscreen plates. The plates were incubated at 37 °C and shaken continuously. Estimates of the OD (600 nm) were made every five minutes for 24 h in LB and 48 h in glucose-limited medium. Normalization, replicate means and error, growth rate, lag, and maximum OD were found using a novel R Bioscreen C analysis tool accessible at https://josheclf.shinyapps.io/bioscreen_app.

Residual growth

After 30 days the cultures were centrifuged and filtered through a 0.22 μm filter. A strain resistant to CHL (Strain 1012 from the US Center for Disease Control’s MuGSI Isolate Bank which is cmlA5 positive), AZM (Strain 1007 from the US Center for Disease Control’s MuGSI Isolate Bank which is mph(A) positive), or just MG1655 were added to the supernatants and allowed to grow for 24 h.

Antibiotic stability

Initial (day 0) and final (day 30) cultures from the long term experiment in the presence of CHL or AZM were centrifugated and the supernatants filtrated. Twenty μL of these supernatants were spotted on a 7 mm blank disk and MG1655 was used as the lawn. After 24 h of incubation, the zone of inhibition was recorded.

Population analysis profile tests

Single colonies of the corresponding strain were inoculated into 10 mL of LB and grown overnight at 37 °C with shaking. Cultures were serially diluted in saline and all dilutions (100 to 10−7) plated on LB agar plates containing 0, 1, 2, 4, 8, and 16 times the MIC of CHL or AZM. Plates were grown at 37 °C for 48 h before the density of surviving colonies was estimated.

Numerical solutions (simulations)

For our numerical analysis of the coupled, ordered differential equations presented (Supplementary Equations 1–12) we used Berkeley Madonna with the parameters presented in the captions of Fig. 7 and Supplementary Figs. 8 and 9. Copies of the Berkeley Madonna programs used for these simulations are available at www.eclf.net. Parameters were obtained from refs. 21,50,51. To calculate the average MIC in our simulations, we take a weighted average of the MIC of each population.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Supplementary Information

Peer Review File

Reporting summary

Source data

Source data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52166-z.

Acknowledgements

We would like to thank Danielle Steed for her discussion of the clinical implications and relevance of this work. We would also like to thank Jason Chen and LM Bradley for their unwavering assistance and support. We would also like to thank the other members of the Levin Lab. BRL would like to thank the US National Institute of General Medical Sciences for their funding support via R35 GM 136407 and the US National Institute of Allergy and Infectious Diseases for their funding support via U19 AI 158080-02.

Author contributions

Conceptualization: T.G.G., B.A.B., J.A.M., F.B., B.R.L. Methodology: T.G.G., B.A.B., J.A.M., I.C.M., B.R.L. Investigation: T.G.G., B.A.B., J.A.M., A.P.S., I.C.M. Visualization: T.G.G., J.A.M. Funding Acquisition: B.R.L. Project Administration: B.R.L. Supervision: B.R.L. Writing– Initial Draft: T.G.G., B.A.B., J.A.M., A.P.S., F.B., B.R.L. Writing– Review & Editing: T.G.G., B.A.B., J.A.M., A.P.S., I.C.M., F.B., B.R.L.

Peer review

Peer review information

Nature Communications thanks Quentin Leclerc, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Data availability

All the data generated are available in this manuscript and its supporting supplementary material. The genomic data generated in this study have been deposited in the NCBI database under accession code PRJNA1032893. Source data are provided with this paper.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Teresa Gil-Gil, Brandon A. Berryhill.
==== Refs
References

1. Bernatová S Following the mechanisms of bacteriostatic versus bactericidal action using Raman spectroscopy Molecules 2013 18 13188 13199 10.3390/molecules181113188 24284484
Bernatová, S. et al. Following the mechanisms of bacteriostatic versus bactericidal action using Raman spectroscopy. Molecules 18, 13188–13199 (2013).24284484 10.3390/molecules181113188
2. Pankey GA Sabath LD Clinical Relevance of Bacteriostatic versus Bactericidal Mechanisms of Action in the Treatment of Gram-Positive Bacterial Infections Clin. Infect. Dis. 2004 38 864 870 10.1086/381972 14999632
Pankey, G. A. & Sabath, L. D. Clinical Relevance of Bacteriostatic versus Bactericidal Mechanisms of Action in the Treatment of Gram-Positive Bacterial Infections. Clin. Infect. Dis. 38, 864–870 (2004).14999632 10.1086/381972
3. Nemeth J Oesch G Kuster SP Bacteriostatic versus bactericidal antibiotics for patients with serious bacterial infections: systematic review and meta-analysis J. Antimicrob. Chemother. 2015 70 382 395 10.1093/jac/dku379 25266070
Nemeth, J., Oesch, G. & Kuster, S. P. Bacteriostatic versus bactericidal antibiotics for patients with serious bacterial infections: systematic review and meta-analysis. J. Antimicrob. Chemother. 70, 382–395 (2015).25266070 10.1093/jac/dku379
4. Regoes RR Pharmacodynamic functions: a multiparameter approach to the design of antibiotic treatment regimens Antimicrob. Agents Chemother. 2004 48 3670 3676 10.1128/AAC.48.10.3670-3676.2004 15388418
Regoes, R. R. et al. Pharmacodynamic functions: a multiparameter approach to the design of antibiotic treatment regimens. Antimicrob. Agents Chemother. 48, 3670–3676 (2004).15388418 10.1128/AAC.48.10.3670-3676.2004
5. William A Killing and regrowth of bacteria in vitro Scand. J. Infect. Dis. Suppl. 1991 74 63 70
William, A. Killing and regrowth of bacteria in vitro. Scand. J. Infect. Dis. Suppl. 74, 63–70 (1991).
6. Bonapace CR Friedrich LV Bosso JA White RL Determination of antibiotic effect in an in vitro pharmacodynamic model: comparison with an established animal model of infection Antimicrobial agents Chemother. 2002 46 3574 3579 10.1128/AAC.46.11.3574-3579.2002
Bonapace, C. R., Friedrich, L. V., Bosso, J. A. & White, R. L. Determination of antibiotic effect in an in vitro pharmacodynamic model: comparison with an established animal model of infection. Antimicrobial agents Chemother. 46, 3574–3579 (2002).10.1128/AAC.46.11.3574-3579.2002
7. Craig WA Ebert SC Killing and regrowth of bacteria in vitro: a review Scand. J. Infect. Dis. Suppl. 1990 74 63 70 2097720
Craig, W. A. & Ebert, S. C. Killing and regrowth of bacteria in vitro: a review. Scand. J. Infect. Dis. Suppl. 74, 63–70 (1990).2097720
8. Huemer M Mairpady Shambat S Brugger SD Zinkernagel AS Antibiotic resistance and persistence—Implications for human health and treatment perspectives EMBO Rep. 2020 21 e51034 10.15252/embr.202051034 33400359
Huemer, M., Mairpady Shambat, S., Brugger, S. D. & Zinkernagel, A. S. Antibiotic resistance and persistence—Implications for human health and treatment perspectives. EMBO Rep. 21, e51034 (2020).33400359 10.15252/embr.202051034
9. Levin-Reisman I Brauner A Ronin I Balaban NQ Epistasis between antibiotic tolerance, persistence, and resistance mutations Proc. Natl. Acad. Sci. USA 2019 116 14734 14739 10.1073/pnas.1906169116 31262806
Levin-Reisman, I., Brauner, A., Ronin, I. & Balaban, N. Q. Epistasis between antibiotic tolerance, persistence, and resistance mutations. Proc. Natl. Acad. Sci. USA 116, 14734–14739 (2019).31262806 10.1073/pnas.1906169116
10. Balaban NQ Publisher Correction: Definitions and guidelines for research on antibiotic persistence Nat. Rev. Microbiol 2019 17 460 10.1038/s41579-019-0207-4 31036919
Balaban, N. Q. et al. Publisher Correction: Definitions and guidelines for research on antibiotic persistence. Nat. Rev. Microbiol 17, 460 (2019).31036919 10.1038/s41579-019-0207-4
11. Wald-Dickler N Holtom P Spellberg B Busting the Myth of “Static vs Cidal”: A Systemic Literature Review Clin. Infect. Dis. 2018 66 1470 1474 10.1093/cid/cix1127 29293890
Wald-Dickler, N., Holtom, P. & Spellberg, B. Busting the Myth of “Static vs Cidal”: A Systemic Literature Review. Clin. Infect. Dis. 66, 1470–1474 (2018).29293890 10.1093/cid/cix1127
12. Rizk, M. L. et al. Considerations for Dose Selection and Clinical Pharmacokinetics/Pharmacodynamics for the Development of Antibacterial Agents. Antimicrob Agents Chemother 63 (2019). 10.1128/aac.02309-18.
13. Berryhill BA What’s the Matter with MICs: Bacterial Nutrition, Limiting Resources, and Antibiotic Pharmacodynamics Microbiol. Spectr. 2023 11 e04091 04022 10.1128/spectrum.04091-22 37130356
Berryhill, B. A. et al. What’s the Matter with MICs: Bacterial Nutrition, Limiting Resources, and Antibiotic Pharmacodynamics. Microbiol. Spectr. 11, e04091–04022 (2023).37130356 10.1128/spectrum.04091-22
14. Lopatkin, A. J. et al. Clinically relevant mutations in core metabolic genes confer antibiotic resistance. Science 371 (2021). 10.1126/science.aba0862.
15. Martela K Kuźniewicz R Pluskiewicz W Tabor E Zagórski P Calcium supplementation from dairy products in the diet of women over the age of 55 from Zabrze (the Silesian Osteoactive Study - nutritional part) Reumatologia 2018 56 382 387 10.5114/reum.2018.80716 30647485
Martela, K., Kuźniewicz, R., Pluskiewicz, W., Tabor, E. & Zagórski, P. Calcium supplementation from dairy products in the diet of women over the age of 55 from Zabrze (the Silesian Osteoactive Study - nutritional part). Reumatologia 56, 382–387 (2018).30647485 10.5114/reum.2018.80716
16. Fernández M Mechanisms of resistance to chloramphenicol in Pseudomonas putida KT2440 Antimicrob. Agents Chemother. 2012 56 1001 1009 10.1128/AAC.05398-11 22143519
Fernández, M. et al. Mechanisms of resistance to chloramphenicol in Pseudomonas putida KT2440. Antimicrob. Agents Chemother. 56, 1001–1009 (2012).22143519 10.1128/AAC.05398-11
17. Heidary M Mechanism of action, resistance, synergism, and clinical implications of azithromycin J. Clin. Lab Anal. 2022 36 e24427 10.1002/jcla.24427 35447019
Heidary, M. et al. Mechanism of action, resistance, synergism, and clinical implications of azithromycin. J. Clin. Lab Anal. 36, e24427 (2022).35447019 10.1002/jcla.24427
18. Cronan JE pBR322 vectors having tetracycline-dependent replication Plasmid 2016 84-85 20 26 10.1016/j.plasmid.2016.02.004 26876942
Cronan, J. E. pBR322 vectors having tetracycline-dependent replication. Plasmid 84-85, 20–26 (2016).26876942 10.1016/j.plasmid.2016.02.004
19. Khusainov I Bactericidal effect of tetracycline in E. coli strain ED1a may be associated with ribosome dysfunction Nat. Commun. 2024 15 4783 10.1038/s41467-024-49084-5 38839776
Khusainov, I. et al. Bactericidal effect of tetracycline in E. coli strain ED1a may be associated with ribosome dysfunction. Nat. Commun. 15, 4783 (2024).38839776 10.1038/s41467-024-49084-5
20. Andersson DI Nicoloff H Hjort K Mechanisms and clinical relevance of bacterial heteroresistance Nat. Rev. Microbiol. 2019 17 479 496 10.1038/s41579-019-0218-1 31235888
Andersson, D. I., Nicoloff, H. & Hjort, K. Mechanisms and clinical relevance of bacterial heteroresistance. Nat. Rev. Microbiol. 17, 479–496 (2019).31235888 10.1038/s41579-019-0218-1
21. Levin BR Theoretical considerations and empirical predictions of the pharmaco- and population dynamics of heteroresistance Proc. Natl. Acad. Sci. USA 2024 121 e2318600121 10.1073/pnas.2318600121 38588431
Levin, B. R. et al. Theoretical considerations and empirical predictions of the pharmaco- and population dynamics of heteroresistance. Proc. Natl. Acad. Sci. USA 121, e2318600121 (2024).38588431 10.1073/pnas.2318600121
22. Baquero F Levin BR Proximate and ultimate causes of the bactericidal action of antibiotics Nat. Rev. Microbiol. 2021 19 123 132 10.1038/s41579-020-00443-1 33024310
Baquero, F. & Levin, B. R. Proximate and ultimate causes of the bactericidal action of antibiotics. Nat. Rev. Microbiol. 19, 123–132 (2021).33024310 10.1038/s41579-020-00443-1
23. Satta G Target for bacteriostatic and bactericidal activities of beta-lactam antibiotics against Escherichia coli resides in different penicillin-binding proteins Antimicrob. Agents Chemother. 1995 39 812 818 10.1128/AAC.39.4.812 7785976
Satta, G. et al. Target for bacteriostatic and bactericidal activities of beta-lactam antibiotics against Escherichia coli resides in different penicillin-binding proteins. Antimicrob. Agents Chemother. 39, 812–818 (1995).7785976 10.1128/AAC.39.4.812
24. Murphy, M. B., Mercer, S. L. & Deweese, J. E. in Advances in Molecular Toxicology Vol. 11 (eds James C. Fishbein & Jacqueline M. Heilman) 203–240 (Elsevier, 2017).
25. Kaur P Agarwal S Datta S Delineating bacteriostatic and bactericidal targets in mycobacteria using IPTG inducible antisense expression PLoS One 2009 4 e5923 10.1371/journal.pone.0005923 19526063
Kaur, P., Agarwal, S. & Datta, S. Delineating bacteriostatic and bactericidal targets in mycobacteria using IPTG inducible antisense expression. PLoS One 4, e5923 (2009).19526063 10.1371/journal.pone.0005923
26. Baudens JG Chabbert YA Rifampicin: a bacteriostatic and bactericidal agent] Pathol. Biol. (Paris) 1969 17 392 397 4307548
Baudens, J. G. & Chabbert, Y. A. Rifampicin: a bacteriostatic and bactericidal agent]. Pathol. Biol. (Paris) 17, 392–397 (1969).4307548
27. Levin, B. R. et al. A Numbers Game: Ribosome Densities, Bacterial Growth, and Antibiotic-Mediated Stasis and Death. mBio 8 (2017). 10.1128/mBio.02253-16.
28. Yang J Barra JT Fung DK Wang JD Bacillus subtilis produces (p)ppGpp in response to the bacteriostatic antibiotic chloramphenicol to prevent its potential bactericidal effect mLife 2022 1 101 113 10.1002/mlf2.12031 38817674
Yang, J., Barra, J. T., Fung, D. K. & Wang, J. D. Bacillus subtilis produces (p)ppGpp in response to the bacteriostatic antibiotic chloramphenicol to prevent its potential bactericidal effect. mLife 1, 101–113 (2022).38817674 10.1002/mlf2.12031
29. Baquero, F., Rodríguez-Beltrán, J. & Levin, B. R. Bacteriostatic cells instead of bacteriostatic antibiotics? Ecoevorxiv (2023). 10.32942/X2R01T.
30. Gomes C Ruiz-Roldán L Mateu J Ochoa TJ Ruiz J Azithromycin resistance levels and mechanisms in Escherichia coli Sci. Rep. 2019 9 6089 10.1038/s41598-019-42423-3 30988366
Gomes, C., Ruiz-Roldán, L., Mateu, J., Ochoa, T. J. & Ruiz, J. Azithromycin resistance levels and mechanisms in Escherichia coli. Sci. Rep. 9, 6089 (2019).30988366 10.1038/s41598-019-42423-3
31. Roggenkamp A Chronic prosthetic hip infection caused by a small-colony variant of Escherichia coli J. Clin. Microbiol 1998 36 2530 2534 10.1128/JCM.36.9.2530-2534.1998 9705387
Roggenkamp, A. et al. Chronic prosthetic hip infection caused by a small-colony variant of Escherichia coli. J. Clin. Microbiol 36, 2530–2534 (1998).9705387 10.1128/JCM.36.9.2530-2534.1998
32. Negishi T Characterization of clinically isolated thymidine-dependent small-colony variants of Escherichia coli producing extended-spectrum β-lactamase J. Med Microbiol. 2018 67 33 39 10.1099/jmm.0.000634 29143727
Negishi, T. et al. Characterization of clinically isolated thymidine-dependent small-colony variants of Escherichia coli producing extended-spectrum β-lactamase. J. Med Microbiol. 67, 33–39 (2018).29143727 10.1099/jmm.0.000634
33. Xia H A yigP mutant strain is a small colony variant of E. coli and shows pleiotropic antibiotic resistance Can. J. Microbiol. 2017 63 961 969 10.1139/cjm-2017-0347 28915357
Xia, H. et al. A yigP mutant strain is a small colony variant of E. coli and shows pleiotropic antibiotic resistance. Can. J. Microbiol. 63, 961–969 (2017).28915357 10.1139/cjm-2017-0347
34. Roggenkamp A Hoffmann H Hornef MW Growth control of small-colony variants by genetic regulation of the hemin uptake system Infect. Immun. 2004 72 2254 2262 10.1128/IAI.72.4.2254-2262.2004 15039350
Roggenkamp, A., Hoffmann, H. & Hornef, M. W. Growth control of small-colony variants by genetic regulation of the hemin uptake system. Infect. Immun. 72, 2254–2262 (2004).15039350 10.1128/IAI.72.4.2254-2262.2004
35. Santos V Hirshfield I The Physiological and Molecular Characterization of a Small Colony Variant of Escherichia coli and Its Phenotypic Rescue PLOS ONE 2016 11 e0157578 10.1371/journal.pone.0157578 27310825
Santos, V. & Hirshfield, I. The Physiological and Molecular Characterization of a Small Colony Variant of Escherichia coli and Its Phenotypic Rescue. PLOS ONE 11, e0157578 (2016).27310825 10.1371/journal.pone.0157578
36. Tashiro Y Eida H Ishii S Futamata H Okabe S Generation of Small Colony Variants in Biofilms by Escherichia coli Harboring a Conjugative F Plasmid Microbes Environ. 2017 32 40 46 10.1264/jsme2.ME16121 28302951
Tashiro, Y., Eida, H., Ishii, S., Futamata, H. & Okabe, S. Generation of Small Colony Variants in Biofilms by Escherichia coli Harboring a Conjugative F Plasmid. Microbes Environ. 32, 40–46 (2017).28302951 10.1264/jsme2.ME16121
37. Hirsch HM Small colony variants of Escherichia coli. Mode of action of copper in variant recovery and population dynamics of cultures containing variants J. Bacteriol. 1961 81 448 458 10.1128/jb.81.3.448-458.1961 13714572
Hirsch, H. M. Small colony variants of Escherichia coli. Mode of action of copper in variant recovery and population dynamics of cultures containing variants. J. Bacteriol. 81, 448–458 (1961).13714572 10.1128/jb.81.3.448-458.1961
38. Dworkin J Understanding the Stringent Response: Experimental Context Matters mBio 2023 14 e0340422 10.1128/mbio.03404-22 36625599
Dworkin, J. Understanding the Stringent Response: Experimental Context Matters. mBio 14, e0340422 (2023).36625599 10.1128/mbio.03404-22
39. Trampari E Functionally distinct mutations within AcrB underpin antibiotic resistance in different lifestyles npj Antimicrobials Resistance 2023 1 2 10.1038/s44259-023-00001-8 38686215
Trampari, E. et al. Functionally distinct mutations within AcrB underpin antibiotic resistance in different lifestyles. npj Antimicrobials Resistance 1, 2 (2023).38686215 10.1038/s44259-023-00001-8
40. Moore SD Sauer RT Revisiting the mechanism of macrolide-antibiotic resistance mediated by ribosomal protein L22 Proc. Natl. Acad. Sci. USA 2008 105 18261 18266 10.1073/pnas.0810357105 19015512
Moore, S. D. & Sauer, R. T. Revisiting the mechanism of macrolide-antibiotic resistance mediated by ribosomal protein L22. Proc. Natl. Acad. Sci. USA 105, 18261–18266 (2008).19015512 10.1073/pnas.0810357105
41. Suryanarayana T Identification by affinity chromatography of Escherichia coli ribosomal proteins that bind erythromycin and chloramphenicol Biochem Int. 1983 7 719 725 6385986
Suryanarayana, T. Identification by affinity chromatography of Escherichia coli ribosomal proteins that bind erythromycin and chloramphenicol. Biochem Int. 7, 719–725 (1983).6385986
42. Matange, N. Highly Contingent Phenotypes of Lon Protease Deficiency in Escherichia coli upon Antibiotic Challenge. J. Bacteriol. 202 (2020). 10.1128/jb.00561-19.
43. Hughes CA Gorabi V Escamilla Y Dean FB Bullard JM Two Forms of Tyrosyl-tRNA Synthetase from Pseudomonas aeruginosa: Characterization and Discovery of Inhibitory Compounds SLAS Discov. 2020 25 1072 1086 10.1177/2472555220934793 32583746
Hughes, C. A., Gorabi, V., Escamilla, Y., Dean, F. B. & Bullard, J. M. Two Forms of Tyrosyl-tRNA Synthetase from Pseudomonas aeruginosa: Characterization and Discovery of Inhibitory Compounds. SLAS Discov. 25, 1072–1086 (2020).32583746 10.1177/2472555220934793
44. Baquero, F. et al. Bacterial Subcellular Architecture, Structural Epistasis, and Antibiotic Resistance. Biology (Basel) 12 (2023). 10.3390/biology12050640.
45. MA W Methods for dilution antimicrobial susceptibility tests for bacteria that grow aerobically: approved standard Clsi (Nccls) 2006 26 M7 A7
MA, W. Methods for dilution antimicrobial susceptibility tests for bacteria that grow aerobically: approved standard. Clsi (Nccls) 26, M7–A7 (2006).
46. Brown DF Brown L Evaluation of the E test, a novel method of quantifying antimicrobial activity J. Antimicrobial Chemother. 1991 27 185 190 10.1093/jac/27.2.185
Brown, D. F. & Brown, L. Evaluation of the E test, a novel method of quantifying antimicrobial activity. J. Antimicrobial Chemother. 27, 185–190 (1991).10.1093/jac/27.2.185
47. Wick RR Judd LM Gorrie CL Holt KE Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads PLoS Comput Biol. 2017 13 e1005595 10.1371/journal.pcbi.1005595 28594827
Wick, R. R., Judd, L. M., Gorrie, C. L. & Holt, K. E. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol. 13, e1005595 (2017).28594827 10.1371/journal.pcbi.1005595
48. Gurevich A Saveliev V Vyahhi N Tesler G QUAST: quality assessment tool for genome assemblies Bioinformatics 2013 29 1072 1075 10.1093/bioinformatics/btt086 23422339
Gurevich, A., Saveliev, V., Vyahhi, N. & Tesler, G. QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075 (2013).23422339 10.1093/bioinformatics/btt086
49. Seemann T Prokka: rapid prokaryotic genome annotation Bioinformatics 2014 30 2068 2069 10.1093/bioinformatics/btu153 24642063
Seemann, T. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069 (2014).24642063 10.1093/bioinformatics/btu153
50. Berryhill BA The book of Lambda does not tell us that naturally occurring lysogens of Escherichia coli are likely to be resistant as well as immune Proc. Natl. Acad. Sci. USA 2023 120 e2212121120 10.1073/pnas.2212121120 36881631
Berryhill, B. A. et al. The book of Lambda does not tell us that naturally occurring lysogens of Escherichia coli are likely to be resistant as well as immune. Proc. Natl. Acad. Sci. USA 120, e2212121120 (2023).36881631 10.1073/pnas.2212121120
51. Chaudhry WN Leaky resistance and the conditions for the existence of lytic bacteriophage PLoS Biol. 2018 16 e2005971 10.1371/journal.pbio.2005971 30114198
Chaudhry, W. N. et al. Leaky resistance and the conditions for the existence of lytic bacteriophage. PLoS Biol. 16, e2005971 (2018).30114198 10.1371/journal.pbio.2005971
