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,274
    views
  • 278
    downloads
  • 30
    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. Chromosomes and Gene Expression
    2. Evolutionary Biology
    Timothy Fuqua, Yiqiao Sun, Andreas Wagner
    Research Article

    Gene regulation is essential for life and controlled by regulatory DNA. Mutations can modify the activity of regulatory DNA, and also create new regulatory DNA, a process called regulatory emergence. Non-regulatory and regulatory DNA contain motifs to which transcription factors may bind. In prokaryotes, gene expression requires a stretch of DNA called a promoter, which contains two motifs called –10 and –35 boxes. However, these motifs may occur in both promoters and non-promoter DNA in multiple copies. They have been implicated in some studies to improve promoter activity, and in others to repress it. Here, we ask whether the presence of such motifs in different genetic sequences influences promoter evolution and emergence. To understand whether and how promoter motifs influence promoter emergence and evolution, we start from 50 ‘promoter islands’, DNA sequences enriched with –10 and –35 boxes. We mutagenize these starting ‘parent’ sequences, and measure gene expression driven by 240,000 of the resulting mutants. We find that the probability that mutations create an active promoter varies more than 200-fold, and is not correlated with the number of promoter motifs. For parent sequences without promoter activity, mutations created over 1500 new –10 and –35 boxes at unique positions in the library, but only ~0.3% of these resulted in de-novo promoter activity. Only ~13% of all –10 and –35 boxes contribute to de-novo promoter activity. For parent sequences with promoter activity, mutations created new –10 and –35 boxes in 11 specific positions that partially overlap with preexisting ones to modulate expression. We also find that –10 and –35 boxes do not repress promoter activity. Overall, our work demonstrates how promoter motifs influence promoter emergence and evolution. It has implications for predicting and understanding regulatory evolution, de novo genes, and phenotypic evolution.

    1. Evolutionary Biology
    Matthew Osmond, Graham Coop
    Research Article Updated

    Spatial patterns in genetic diversity are shaped by individuals dispersing from their parents and larger-scale population movements. It has long been appreciated that these patterns of movement shape the underlying genealogies along the genome leading to geographic patterns of isolation-by-distance in contemporary population genetic data. However, extracting the enormous amount of information contained in genealogies along recombining sequences has, until recently, not been computationally feasible. Here, we capitalize on important recent advances in genome-wide gene-genealogy reconstruction and develop methods to use thousands of trees to estimate per-generation dispersal rates and to locate the genetic ancestors of a sample back through time. We take a likelihood approach in continuous space using a simple approximate model (branching Brownian motion) as our prior distribution of spatial genealogies. After testing our method with simulations we apply it to Arabidopsis thaliana. We estimate a dispersal rate of roughly 60 km2/generation, slightly higher across latitude than across longitude, potentially reflecting a northward post-glacial expansion. Locating ancestors allows us to visualize major geographic movements, alternative geographic histories, and admixture. Our method highlights the huge amount of information about past dispersal events and population movements contained in genome-wide genealogies.