Price equation captures the role of drug interactions and collateral effects in the evolution of multidrug resistance

  1. Erida Gjini  Is a corresponding author
  2. Kevin B Wood  Is a corresponding author
  1. University of Lisbon, Portugal, Portugal
  2. University of Michigan, United States

Abstract

Bacterial adaptation to antibiotic combinations depends on the joint inhibitory effects of the two drugs (drug interaction, DI) and how resistance to one drug impacts resistance to the other (collateral effects, CE). Here we model these evolutionary dynamics on two-dimensional phenotype spaces that leverage scaling relations between the drug-response surfaces of drug sensitive (ancestral) and drug resistant (mutant) populations. We show that evolved resistance to the component drugs-and in turn, the adaptation of growth rate-is governed by a Price equation whose covariance terms encode geometric features of both the two-drug response surface (DI) in ancestral cells and the correlations between resistance levels to those drugs (CE). Within this framework, mean evolutionary trajectories reduce to a type of weighted gradient dynamics, with the drug interaction dictating the shape of the underlying landscape and the collateral effects constraining the motion on those landscapes. We also demonstrate how constraints on available mutational pathways can be incorporated into the framework, adding a third key driver of evolution. Our results clarify the complex relationship between drug interactions and collateral effects in multi-drug environments and illustrate how specific dosage combinations can shift the weighting of these two effects, leading to different and temporally-explicit selective outcomes.

Data availability

Data used in this paper was taken from a public repository:Dean, Ziah; Maltas, Jeff; Wood, Kevin (2020), Antibiotic interactions shape short-term evolution of resistance in Enterococcus faecalis, Dryad, Dataset, https://doi.org/10.5061/dryad.j3tx95x92There are no restrictions on any new results.

The following previously published data sets were used

Article and author information

Author details

  1. Erida Gjini

    Center for Computational and Stochastic Mathematics, Instituto Superior Tecnico,, University of Lisbon, Portugal, Lisbon, Portugal
    For correspondence
    erida.gjini@tecnico.ulisboa.pt
    Competing interests
    The authors declare that no competing interests exist.
  2. Kevin B Wood

    Department of Biophysics, University of Michigan, Ann Arbor, United States
    For correspondence
    kbwood@umich.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0985-7401

Funding

National Institutes of Health (1R35GM124875)

  • Kevin B Wood

National Science Foundation (1553028)

  • Kevin B Wood

Fundação Luso-Americana para o Desenvolvimento (274/2016)

  • Erida Gjini

Instituto Gulbenkian de Ciencia

  • Erida Gjini

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Copyright

© 2021, Gjini & Wood

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

Metrics

  • 2,316
    views
  • 286
    downloads
  • 31
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

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. Erida Gjini
  2. Kevin B Wood
(2021)
Price equation captures the role of drug interactions and collateral effects in the evolution of multidrug resistance
eLife 10:e64851.
https://doi.org/10.7554/eLife.64851

Share this article

https://doi.org/10.7554/eLife.64851

Further reading

    1. Ecology
    2. Evolutionary Biology
    Justine Boutry, Océane Rieu ... Fréderic Thomas
    Research Article

    While host phenotypic manipulation by parasites is a widespread phenomenon, whether tumors, which can be likened to parasite entities, can also manipulate their hosts is not known. Theory predicts that this should nevertheless be the case, especially when tumors (neoplasms) are transmissible. We explored this hypothesis in a cnidarian Hydra model system, in which spontaneous tumors can occur in the lab, and lineages in which such neoplastic cells are vertically transmitted (through host budding) have been maintained for over 15 years. Remarkably, the hydras with long-term transmissible tumors show an unexpected increase in the number of their tentacles, allowing for the possibility that these neoplastic cells can manipulate the host. By experimentally transplanting healthy as well as neoplastic tissues derived from both recent and long-term transmissible tumors, we found that only the long-term transmissible tumors were able to trigger the growth of additional tentacles. Also, supernumerary tentacles, by permitting higher foraging efficiency for the host, were associated with an increased budding rate, thereby favoring the vertical transmission of tumors. To our knowledge, this is the first evidence that, like true parasites, transmissible tumors can evolve strategies to manipulate the phenotype of their host.

    1. Evolutionary Biology
    2. Microbiology and Infectious Disease
    Zach Hensel
    Short Report

    Accurate estimation of the effects of mutations on SARS-CoV-2 viral fitness can inform public-health responses such as vaccine development and predicting the impact of a new variant; it can also illuminate biological mechanisms including those underlying the emergence of variants of concern. Recently, Lan et al. reported a model of SARS-CoV-2 secondary structure and its underlying dimethyl sulfate reactivity data (Lan et al., 2022). I investigated whether base reactivities and secondary structure models derived from them can explain some variability in the frequency of observing different nucleotide substitutions across millions of patient sequences in the SARS-CoV-2 phylogenetic tree. Nucleotide basepairing was compared to the estimated ‘mutational fitness’ of substitutions, a measurement of the difference between a substitution’s observed and expected frequency that is correlated with other estimates of viral fitness (Bloom and Neher, 2023). This comparison revealed that secondary structure is often predictive of substitution frequency, with significant decreases in substitution frequencies at basepaired positions. Focusing on the mutational fitness of C→U, the most common type of substitution, I describe C→U substitutions at basepaired positions that characterize major SARS-CoV-2 variants; such mutations may have a greater impact on fitness than appreciated when considering substitution frequency alone.