Antibiotic Resistance: Finding the right sequence of drugs

Rapidly switching between similar antibiotics may help to slow down the evolution of resistance.
  1. Anh Huynh
  2. Kevin B Wood  Is a corresponding author
  1. Department of Biophysics, University of Michigan, United States
  2. Department of Physics, University of Michigan, United States

In order to survive, many living organisms need to be able to adapt to their ever-changing environment. These past experiences shape the behavior of creatures big and small – from mutations in single genes to neurological changes that underlie memory formation in primates.

There is substantial evidence that environmental history affects how bacteria respond to antibiotics (Barbosa et al., 2019; Card et al., 2019; Nichol et al., 2019; Santos-Lopez et al., 2019; Yen and Papin, 2017). This has led researchers to suggest that switching between drugs over time could help slow down antibiotic resistance, as this will force the bacteria into a scenario where their existing solutions (resistance to the current drug) cannot protect them from tomorrow’s problem (a new drug). However, this method has led to mixed results (Abel zur Wiesch et al., 2014; Imamovic et al., 2018), and optimizing this approach is challenging, in part because it is unclear which features of the antibiotic sequence are the most important and guarantee the best results.

Now, in eLife, Hinrich Schulenburg (University of Kiel and Max Planck Institute for Evolutionary Biology) and colleagues – including Aditi Batra and Roderich Roemhild as joint first authors – report the results of experiments on the multi-drug resistant bacteria Pseudomonas aeruginosa (Batra et al., 2021). The team (who are based in Austria and Germany) exposed the bacteria to various sequences of three antibiotics that belong to commonly used classes of drugs: one class targets the ribosome, one targets DNA gyrase, and one targets the cell wall. In some cases, three drugs from different classes were used (that is, a heterogeneous sequence), and in some cases all three drugs belonged to the same class (a homogenous sequence). Batra et al. also varied the temporal properties of each sequence by switching between the drugs rapidly, slowly, or in a random order. The growth rate, phenotypic resistance levels and population genetics of the evolved populations were then analyzed to determine which sequences of drugs were the most effective at eliminating the bacteria (Figure 1A).

Experimental and biological features of effective drug sequences.

(A) Batra et al. applied different sequences of antibiotics to 756 populations of P. aeruginosa (top panel). The bacteria were treated with either a single drug (monotherapy; row 1), or three antibiotics which were switched rapidly (row 2), slowly (row 3) or in a random order (row 4): the three drugs were either from the same class (homogeneous) or from different classes (heterogeneous). This experiment revealed that fast (blue line) and random (green line) switching between three homogeneous beta-lactam drugs reduced bacteria growth and resulted in higher levels of extinction (bottom graph). (B) The effects of the different sequences are also impacted by biological features. (Top panel) When sensitive bacteria (shown in purple) are treated with the first drug, some cells will evolve genetic changes that make them resistant to the antibiotic treatment (shown in green). These evolutionary changes can lead to collateral effects that make the bacteria less (top arrow), equally (middle arrow) or more (bottom arrow) resistant to the second drug. (Bottom panel) Treatment with the first drug may also lead to negative hysteresis, when short-term physiological changes enhance the bacteria’s response to the second drug (right), leading to more cell death in the population compared to bacteria not pre-treated with the first drug (left).

Image credit: Anh Huynh.

In addition to these experimental parameters, the impact of different antibiotic sequences could also depend on how the population biologically responds to consecutive drug exposures (Figure 1B). For example, the genetic changes bacteria evolve in response to one antibiotic can lead to collateral effects that increase the population’s resistance or sensitivity to another drug. Because collateral resistance occurs more frequently between drugs of the same class, heterogeneous sequences of antibiotics are thought to be more effective at eliminating bacteria (Imamovic and Sommer, 2013; Lázár et al., 2013; Maltas and Wood, 2019; Pál et al., 2015; Lázár et al., 2014). Treatment with unrelated drugs has also been shown to favor negative hysteresis, which is when short-term physiological changes induced by one antibiotic enhance susceptibility towards another. Indeed, a recent study found that rapid switching between antibiotics from different classes promoted extinction of bacterial populations, even when the drugs were used at sub-inhibitory levels (Roemhild et al., 2018).

However, Batra et al. found that homogenous sequences of beta-lactams (a class of antibiotics that target the cell wall) were surprisingly more effective at clearing bacteria. The experiments also revealed that extinction tended to occur early in the treatment and was less effective when drugs were switched more slowly. A particular heterogeneous set of drugs also tended to not eliminate bacteria, indicating that heterogeneity, alone, does not guarantee success.

So, what are the important characteristics of a ‘good’ antibiotic sequence? To answer this question, Batra et al. used a common class of statistical models to probe for signatures of successful sequences. They found that extinction was strongly favored by two biological properties — low rates of spontaneous resistance and low levels of collateral resistance — and was also enhanced when switching between drugs was fast or random. This suggests that although thefindings of Batra et al. contradict the proposed benefits of using unrelated drugs, they still validate a portion of the underlying logic: using antibiotics with strong collateral effects and hysteresis enhances the impact of sequence therapy. It just turns out, however, that drugs with these characteristics are not always from different classes.

This study offers insight into how past antibiotic exposure shapes the response of bacterial populations to new challenges. In doing so, it provides a roadmap for future studies investigating how even the simplest organisms harbor signatures of past challenges and potential evolutionary solutions.


Article and author information

Author details

  1. Anh Huynh

    Anh Huynh is in the Department of Biophysics, University of Michigan, Ann Arbor, United States

    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8891-4234
  2. Kevin B Wood

    Kevin B Wood is in the Department of Biophysics and Department of Physics, University of Michigan, Ann Arbor, United States

    For correspondence
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0985-7401

Publication history

  1. Version of Record published: September 9, 2021 (version 1)


© 2021, Huynh and Wood

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.


  • 996
    Page views
  • 73
  • 1

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Anh Huynh
  2. Kevin B Wood
Antibiotic Resistance: Finding the right sequence of drugs
eLife 10:e72562.
  1. Further reading

Further reading

    1. Developmental Biology
    2. Evolutionary Biology
    Kwi Shan Seah, Vinodkumar Saranathan
    Research Article

    The study of color patterns in the animal integument is a fundamental question in biology, with many lepidopteran species being exemplary models in this endeavor due to their relative simplicity and elegance. While significant advances have been made in unraveling the cellular and molecular basis of lepidopteran pigmentary coloration, the morphogenesis of wing scale nanostructures involved in structural color production is not well understood. Contemporary research on this topic largely focuses on a few nymphalid model taxa (e.g., Bicyclus, Heliconius), despite an overwhelming diversity in the hierarchical nanostructural organization of lepidopteran wing scales. Here, we present a time-resolved, comparative developmental study of hierarchical scale nanostructures in Parides eurimedes and five other papilionid species. Our results uphold the putative conserved role of F-actin bundles in acting as spacers between developing ridges, as previously documented in several nymphalid species. Interestingly, while ridges are developing in P. eurimedes, plasma membrane manifests irregular mesh-like crossribs characteristic of Papilionidae, which delineate the accretion of cuticle into rows of planar disks in between ridges. Once the ridges have grown, disintegrating F-actin bundles appear to reorganize into a network that supports the invagination of plasma membrane underlying the disks, subsequently forming an extruded honeycomb lattice. Our results uncover a previously undocumented role for F-actin in the morphogenesis of complex wing scale nanostructures, likely specific to Papilionidae.

    1. Evolutionary Biology
    Hironori Funabiki, Isabel E Wassing ... Thomas Carroll
    Research Article

    5-Methylcytosine (5mC) and DNA methyltransferases (DNMTs) are broadly conserved in eukaryotes but are also frequently lost during evolution. The mammalian SNF2 family ATPase HELLS and its plant ortholog DDM1 are critical for maintaining 5mC. Mutations in HELLS, its activator CDCA7, and the de novo DNA methyltransferase DNMT3B, cause immunodeficiency-centromeric instability-facial anomalies (ICF) syndrome, a genetic disorder associated with the loss of DNA methylation. We here examine the coevolution of CDCA7, HELLS and DNMTs. While DNMT3, the maintenance DNA methyltransferase DNMT1, HELLS, and CDCA7 are all highly conserved in vertebrates and green plants, they are frequently co-lost in other evolutionary clades. The presence-absence patterns of these genes are not random; almost all CDCA7 harboring eukaryote species also have HELLS and DNMT1 (or another maintenance methyltransferase, DNMT5). Coevolution of presence-absence patterns (CoPAP) analysis in Ecdysozoa further indicates coevolutionary linkages among CDCA7, HELLS, DNMT1 and its activator UHRF1. We hypothesize that CDCA7 becomes dispensable in species that lost HELLS or DNA methylation, and/or the loss of CDCA7 triggers the replacement of DNA methylation by other chromatin regulation mechanisms. Our study suggests that a unique specialized role of CDCA7 in HELLS-dependent DNA methylation maintenance is broadly inherited from the last eukaryotic common ancestor.