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.

Reviewing Editor

  1. C. Brandon Ogbunugafor, Yale University, United States

Version history

  1. Received: November 13, 2020
  2. Accepted: July 8, 2021
  3. Accepted Manuscript published: July 22, 2021 (version 1)
  4. Accepted Manuscript updated: July 26, 2021 (version 2)
  5. Version of Record published: August 3, 2021 (version 3)

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

  • 1,976
    Page views
  • 242
    Downloads
  • 7
    Citations

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. 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

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.