1
|
Moawad A, Abbara A, Bitbol AF. Evolution of cooperation in deme-structured populations on graphs. Phys Rev E 2024; 109:024307. [PMID: 38491653 DOI: 10.1103/physreve.109.024307] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2023] [Accepted: 12/19/2023] [Indexed: 03/18/2024]
Abstract
Understanding how cooperation can evolve in populations despite its cost to individual cooperators is an important challenge. Models of spatially structured populations with one individual per node of a graph have shown that cooperation, modeled via the prisoner's dilemma, can be favored by natural selection. These results depend on microscopic update rules, which determine how birth, death, and migration on the graph are coupled. Recently, we developed coarse-grained models of spatially structured populations on graphs, where each node comprises a well-mixed deme, and where migration is independent from division and death, thus bypassing the need for update rules. Here, we study the evolution of cooperation in these models in the rare-migration regime, within the prisoner's dilemma. We find that cooperation is not favored by natural selection in these coarse-grained models on graphs where overall deme fitness does not directly impact migration from a deme. This is due to a separation of scales, whereby cooperation occurs at a local level within demes, while spatial structure matters between demes.
Collapse
Affiliation(s)
- Alix Moawad
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland and SIB Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland
| | - Alia Abbara
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland and SIB Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland
| | - Anne-Florence Bitbol
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland and SIB Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland
| |
Collapse
|
2
|
Sohel Mondal S, Ray A, Chakraborty S. Hypochaos prevents tragedy of the commons in discrete-time eco-evolutionary game dynamics. CHAOS (WOODBURY, N.Y.) 2024; 34:023122. [PMID: 38377296 DOI: 10.1063/5.0190800] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/08/2023] [Accepted: 01/22/2024] [Indexed: 02/22/2024]
Abstract
While quite a few recent papers have explored game-resource feedback using the framework of evolutionary game theory, almost all the studies are confined to using time-continuous dynamical equations. Moreover, in such literature, the effect of ubiquitous chaos in the resulting eco-evolutionary dynamics is rather missing. Here, we present a deterministic eco-evolutionary discrete-time dynamics in generation-wise non-overlapping population of two types of harvesters-one harvesting at a faster rate than the other-consuming a self-renewing resource capable of showing chaotic dynamics. In the light of our finding that sometimes chaos is confined exclusively to either the dynamics of the resource or that of the consumer fractions, an interesting scenario is realized: The resource state can keep oscillating chaotically, and hence, it does not vanish to result in the tragedy of the commons-extinction of the resource due to selfish indiscriminate exploitation-and yet the consumer population, whose dynamics depends directly on the state of the resource, may end up being composed exclusively of defectors, i.e., high harvesters. This appears non-intuitive because it is well known that prevention of tragedy of the commons usually requires substantial cooperation to be present.
Collapse
Affiliation(s)
- Samrat Sohel Mondal
- Department of Physics, Indian Institute of Technology Kanpur, Uttar Pradesh 208016, India
| | - Avishuman Ray
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, USA
| | - Sagar Chakraborty
- Department of Physics, Indian Institute of Technology Kanpur, Uttar Pradesh 208016, India
| |
Collapse
|
3
|
Park HJ, Hilbe C, Nowak MA, Kim BJ, Jeong HC. Vacancies in growing habitats promote the evolution of cooperation. J Theor Biol 2023; 575:111629. [PMID: 37802182 DOI: 10.1016/j.jtbi.2023.111629] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2023] [Revised: 09/23/2023] [Accepted: 09/28/2023] [Indexed: 10/08/2023]
Abstract
We study evolutionary game dynamics in a growing habitat with vacancies. Fitness is determined by the global effect of the environment and a local prisoner's dilemma among neighbors. We study population growth on a one-dimensional lattice and analyze how the environment affects evolutionary competition. As the environment becomes harsh, an absorbing phase transition from growing populations to extinction occurs. The transition point depends on which strategies are present in the population. In particular, we find a 'cooperative window' in parameter space, where only cooperators can survive. A mutant defector in a cooperative community might briefly proliferate, but over time naturally occurring vacancies separate cooperators from defectors, thereby driving defectors to extinction. Our model reveals that vacancies provide a strong boost for cooperation by spatial selection.
Collapse
Affiliation(s)
- Hye Jin Park
- Department of Physics, Inha University, Incheon, 22212, Republic of Korea.
| | - Christian Hilbe
- Max Planck Research Group 'Dynamics of Social Behavior', Max Planck Institute for Evolutionary Biology, Plön, 24306, Germany
| | - Martin A Nowak
- Department of Mathematics, Harvard University, Cambridge, MA, 02138, United States; Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, 02138, United States
| | - Beom Jun Kim
- Department of Physics, Sungkyunkwan University, Suwon, 16419, Republic of Korea
| | - Hyeong-Chai Jeong
- Department of Physics and Astronomy, Sejong University, Seoul, 05006, Republic of Korea.
| |
Collapse
|
4
|
Arenzon JJ, Peliti L. Emergent cooperative behavior in transient compartments. Phys Rev E 2023; 108:034409. [PMID: 37849208 DOI: 10.1103/physreve.108.034409] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2023] [Accepted: 09/07/2023] [Indexed: 10/19/2023]
Abstract
We introduce a minimal model of multilevel selection on structured populations, considering the interplay between game theory and population dynamics. Through a bottleneck process, finite groups are formed with cooperators and defectors sampled from an infinite pool. After the fragmentation, these transient compartments grow until the maximal number of individuals per compartment is attained. Eventually, all compartments are merged and well mixed, and the whole process is repeated. We show that cooperators, even if interacting only through mean-field intragroup interactions that favor defectors, may perform well because of the intergroup competition and the size diversity among the compartments. These cycles of isolation and coalescence may therefore be important in maintaining diversity among different species or strategies and may help to understand the underlying mechanisms of the scaffolding processes in the transition to multicellularity.
Collapse
Affiliation(s)
- Jeferson J Arenzon
- Instituto de Física, Universidade Federal do Rio Grande do Sul, CP 15051, 91501-970 Porto Alegre RS, Brazil
- Instituto Nacional de Ciência e Tecnologia-Sistemas Complexos, 22290-180 Rio de Janeiro RJ, Brazil
| | - Luca Peliti
- Santa Marinella Research Institute, 00058 Santa Marinella (RM), Italy
| |
Collapse
|
5
|
Park J. Correlation between the formation of new competing group and spatial scale for biodiversity in the evolutionary dynamics of cyclic competition. CHAOS (WOODBURY, N.Y.) 2022; 32:081101. [PMID: 36049957 DOI: 10.1063/5.0102416] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Accepted: 07/15/2022] [Indexed: 06/15/2023]
Abstract
Securing space for species breeding is important in the evolution and maintenance of life in ecological sciences, and an increase in the number of competing species may cause frequent competition and conflict among the population in securing such spaces in a given area. In particular, for cyclically competing species, which can be described by the metaphor of rock-paper-scissors game, most of the previous works in microscopic frameworks have been studied with the initially given three species without any formation of additional competing species, and the phase transition of biodiversity via mobility from coexistence to extinction has never been changed by a change of spatial scale. In this regard, we investigate the relationship between spatial scales and species coexistence in the spatial cyclic game by considering the emergence of a new competing group by mutation. For different spatial scales, our computations reveal that coexistence can be more sensitive to spatial scales and may require larger spaces for frequencies of interactions. By exploiting the calculation of the coexistence probability from Monte-Carlo simulations, we obtain that certain interaction ranges for coexistence can be affected by both spatial scales and mobility, and spatial patterns for coexistence can appear in different ways. Since the issue of spatial scale is important for species survival as competing populations increase, we expect our results to have broad applications in the fields of social and ecological sciences.
Collapse
Affiliation(s)
- Junpyo Park
- Department of Applied Mathematics, Kyung Hee University, Yongin 17104, Republic of Korea
| |
Collapse
|
6
|
Feng TJ, Li C, Zheng XD, Lessard S, Tao Y. Stochastic replicator dynamics and evolutionary stability. Phys Rev E 2022; 105:044403. [PMID: 35590672 DOI: 10.1103/physreve.105.044403] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2021] [Accepted: 03/24/2022] [Indexed: 06/15/2023]
Abstract
To develop the concept of evolutionary stability in a stochastic environment, we investigate the continuous-time dynamics of a two-phenotype linear evolutionary game with generally correlated random payoffs in pairwise interactions. By using the Gram-Schmidt orthogonalization procedure and Itô's formula, we deduce a stochastic differential equation for the phenotype frequencies that extends the replicator equation, called the stochastic replicator equation. We give conditions for stochastic stability of a fixation state or a constant interior equilibrium point with respect to the stochastic dynamics of the two phenotypes. We show that, if a fixation state is stochastically stable, then the pure strategy corresponding to this fixation state must be stochastically evolutionarily stable with respect to mixed strategies. However, this is not the case for a mixed strategy that corresponds to a stochastically stable constant interior equilibrium point with respect to the two phenotypes.
Collapse
Affiliation(s)
- Tian-Jiao Feng
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Cong Li
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an 710072, China
| | - Xiu-Deng Zheng
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Sabin Lessard
- Department of Mathematics and Statistics, University of Montreal, Montreal H3C 3J7, Canada
| | - Yi Tao
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an 710072, China
- Institute of Biomedical Research, Yunnan University, Kunming 650091, China
| |
Collapse
|
7
|
Feng TJ, Mei J, Li C, Zheng XD, Lessard S, Tao Y. Stochastic evolutionary stability in matrix games with random payoffs. Phys Rev E 2022; 105:034303. [PMID: 35428156 DOI: 10.1103/physreve.105.034303] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2021] [Accepted: 02/23/2022] [Indexed: 06/14/2023]
Abstract
Evolutionary game theory and the concept of an evolutionarily stable strategy have been not only extensively developed and successfully applied to explain the evolution of animal behavior, but also widely used in economics and social sciences. Recently, in order to reveal the stochastic dynamical properties of evolutionary games in randomly fluctuating environments, the concept of stochastic evolutionary stability based on conditions for stochastic local stability for a fixation state was developed in the context of a symmetric matrix game with two phenotypes and random payoffs in pairwise interactions [Zheng et al., Phys. Rev. E 96, 032414 (2017)2470-004510.1103/PhysRevE.96.032414]. In this paper, we extend this study to more general situations, namely, multiphenotype symmetric as well as asymmetric matrix games with random payoffs. Conditions for stochastic local stability and stochastic evolutionary stability are established. Conditions for a fixation state to be stochastically unstable and almost everywhere stochastically unstable are distinguished in a multiphenotype setting according to the initial population state. Our results provide some alternative perspective and a more general theoretical framework for a better understanding of the evolution of animal behavior in a stochastic environment.
Collapse
Affiliation(s)
- Tian-Jiao Feng
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Jie Mei
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Cong Li
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an 710072, China
| | - Xiu-Deng Zheng
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Sabin Lessard
- Department of Mathematics and Statistics, University of Montreal, Montreal, Quebec, H3C 3J7, Canada
| | - Yi Tao
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an 710072, China
- Institute of Biomedical Research, Yunnan University, Kunming 650091, China
| |
Collapse
|
8
|
Balasekaran M, Johanis M, Rychtář J, Taylor D, Zhu J. Quasi-neutral evolution in populations under small demographic fluctuations. J Theor Biol 2022; 538:111040. [DOI: 10.1016/j.jtbi.2022.111040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2021] [Revised: 01/11/2022] [Accepted: 01/19/2022] [Indexed: 10/19/2022]
|
9
|
Das Bairagya J, Mondal SS, Chowdhury D, Chakraborty S. Game-environment feedback dynamics in growing population: Effect of finite carrying capacity. Phys Rev E 2021; 104:044407. [PMID: 34781515 DOI: 10.1103/physreve.104.044407] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2021] [Accepted: 09/28/2021] [Indexed: 11/07/2022]
Abstract
The tragedy of the commons (TOC) is an unfortunate situation where a shared resource is exhausted due to uncontrolled exploitation by the selfish individuals of a population. Recently, the paradigmatic replicator equation has been used in conjunction with a phenomenological equation for the state of the shared resource to gain insight into the influence of the games on the TOC. The replicator equation, by construction, models a fixed infinite population undergoing microevolution. Thus, it is unable to capture any effect of the population growth and the carrying capacity of the population although the TOC is expected to be dependent on the size of the population. Therefore, in this paper, we present a mathematical framework that incorporates the density dependent payoffs and the logistic growth of the population in the eco-evolutionary dynamics modeling the game-resource feedback. We discover a bistability in the dynamics: a finite carrying capacity can either avert or cause the TOC depending on the initial states of the resource and the initial fraction of cooperators. In fact, depending on the type of strategic game-theoretic interaction, a finite carrying capacity can either avert or cause the TOC when it is exactly the opposite for the corresponding case with infinite carrying capacity.
Collapse
Affiliation(s)
- Joy Das Bairagya
- Department of Physics, Indian Institute of Technology, Kanpur 208016, India
| | | | | | - Sagar Chakraborty
- Department of Physics, Indian Institute of Technology, Kanpur 208016, India
| |
Collapse
|
10
|
Solving the chemical master equation for monomolecular reaction systems and beyond: a Doi-Peliti path integral view. J Math Biol 2021; 83:48. [PMID: 34635944 DOI: 10.1007/s00285-021-01670-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2019] [Revised: 09/01/2021] [Accepted: 09/09/2021] [Indexed: 10/20/2022]
Abstract
The chemical master equation (CME) is a fundamental description of interacting molecules commonly used to model chemical kinetics and noisy gene regulatory networks. Exact time-dependent solutions of the CME-which typically consists of infinitely many coupled differential equations-are rare, and are valuable for numerical benchmarking and getting intuition for the behavior of more complicated systems. Jahnke and Huisinga's landmark calculation of the exact time-dependent solution of the CME for monomolecular reaction systems is one of the most general analytic results known; however, it is hard to generalize, because it relies crucially on special properties of monomolecular reactions. In this paper, we rederive Jahnke and Huisinga's result on the time-dependent probability distribution and moments of monomolecular reaction systems using the Doi-Peliti path integral approach, which reduces solving the CME to evaluating many integrals. While the Doi-Peliti approach is less intuitive, it is also more mechanical, and hence easier to generalize. To illustrate how the Doi-Peliti approach can go beyond the method of Jahnke and Huisinga, we also find an explicit and exact time-dependent solution to a problem involving an autocatalytic reaction that Jahnke and Huisinga identified as not solvable using their method. Most interestingly, we are able to find a formal exact time-dependent solution for any CME whose list of reactions involves only zero and first order reactions, which may be the most general result currently known. This formal solution also yields a useful algorithm for efficiently computing numerical solutions to CMEs of this type.
Collapse
|
11
|
Pattni K, Overton CE, Sharkey KJ. Evolutionary graph theory derived from eco-evolutionary dynamics. J Theor Biol 2021; 519:110648. [PMID: 33636202 DOI: 10.1016/j.jtbi.2021.110648] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2020] [Revised: 02/15/2021] [Accepted: 02/18/2021] [Indexed: 11/28/2022]
Abstract
A biologically motivated individual-based framework for evolution in network-structured populations is developed that can accommodate eco-evolutionary dynamics. This framework is used to construct a network birth and death model. The evolutionary graph theory model, which considers evolutionary dynamics only, is derived as a special case, highlighting additional assumptions that diverge from real biological processes. This is achieved by introducing a negative ecological feedback loop that suppresses ecological dynamics by forcing births and deaths to be coupled. We also investigate how fitness, a measure of reproductive success used in evolutionary graph theory, is related to the life-history of individuals in terms of their birth and death rates. In simple networks, these ecologically motivated dynamics are used to provide new insight into the spread of adaptive mutations, both with and without clonal interference. For example, the star network, which is known to be an amplifier of selection in evolutionary graph theory, can inhibit the spread of adaptive mutations when individuals can die naturally.
Collapse
Affiliation(s)
- Karan Pattni
- Department of Mathematical Sciences, University of Liverpool, United Kingdom.
| | | | - Kieran J Sharkey
- Department of Mathematical Sciences, University of Liverpool, United Kingdom.
| |
Collapse
|
12
|
Marrec L, Bitbol AF. Adapt or Perish: Evolutionary Rescue in a Gradually Deteriorating Environment. Genetics 2020; 216:573-583. [PMID: 32855198 PMCID: PMC7536851 DOI: 10.1534/genetics.120.303624] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2020] [Accepted: 08/24/2020] [Indexed: 12/31/2022] Open
Abstract
We investigate the evolutionary rescue of a microbial population in a gradually deteriorating environment, through a combination of analytical calculations and stochastic simulations. We consider a population destined for extinction in the absence of mutants, which can survive only if mutants sufficiently adapted to the new environment arise and fix. We show that mutants that appear later during the environment deterioration have a higher probability to fix. The rescue probability of the population increases with a sigmoidal shape when the product of the carrying capacity and of the mutation probability increases. Furthermore, we find that rescue becomes more likely for smaller population sizes and/or mutation probabilities if the environment degradation is slower, which illustrates the key impact of the rapidity of environment degradation on the fate of a population. We also show that our main conclusions are robust across various types of adaptive mutants, including specialist and generalist ones, as well as mutants modeling antimicrobial resistance evolution. We further express the average time of appearance of the mutants that do rescue the population and the average extinction time of those that do not. Our methods can be applied to other situations with continuously variable fitnesses and population sizes, and our analytical predictions are valid in the weak-to-moderate mutation regime.
Collapse
Affiliation(s)
- Loïc Marrec
- Sorbonne Université, CNRS, Institut de Biologie Paris-Seine, Laboratoire Jean Perrin (UMR 8237), 75005 Paris, France
| | - Anne-Florence Bitbol
- Sorbonne Université, CNRS, Institut de Biologie Paris-Seine, Laboratoire Jean Perrin (UMR 8237), 75005 Paris, France
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| |
Collapse
|
13
|
Marcaccioli R, Livan G. Correspondence between temporal correlations in time series, inverse problems, and the spherical model. Phys Rev E 2020; 102:012112. [PMID: 32795068 DOI: 10.1103/physreve.102.012112] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2020] [Accepted: 06/16/2020] [Indexed: 11/07/2022]
Abstract
In this paper we employ methods from statistical mechanics to model temporal correlations in time series. We put forward a methodology based on the maximum entropy principle to generate ensembles of time series constrained to preserve part of the temporal structure of an empirical time series of interest. We show that a constraint on the lag-one autocorrelation can be fully handled analytically and corresponds to the well-known spherical model of a ferromagnet. We then extend such a model to include constraints on more complex temporal correlations by means of perturbation theory, showing that this leads to substantial improvements in capturing the lag-one autocorrelation in the variance. We apply our approach on synthetic data and illustrate how it can be used to formulate expectations on the future values of a data-generating process.
Collapse
Affiliation(s)
- Riccardo Marcaccioli
- Department of Computer Science, University College London, 66-72 Gower Street, London WC1E 6EA, United Kingdom
| | - Giacomo Livan
- Department of Computer Science, University College London, 66-72 Gower Street, London WC1E 6EA, United Kingdom.,Systemic Risk Centre, London School of Economics and Political Sciences, Houghton Street, London WC2A 2AE, United Kingdom
| |
Collapse
|
14
|
Taitelbaum A, West R, Assaf M, Mobilia M. Population Dynamics in a Changing Environment: Random versus Periodic Switching. PHYSICAL REVIEW LETTERS 2020; 125:048105. [PMID: 32794803 DOI: 10.1103/physrevlett.125.048105] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/13/2020] [Revised: 05/13/2020] [Accepted: 06/23/2020] [Indexed: 06/11/2023]
Abstract
Environmental changes greatly influence the evolution of populations. Here, we study the dynamics of a population of two strains, one growing slightly faster than the other, competing for resources in a time-varying binary environment modeled by a carrying capacity switching either randomly or periodically between states of abundance and scarcity. The population dynamics is characterized by demographic noise (birth and death events) coupled to a varying environment. We elucidate the similarities and differences of the evolution subject to a stochastically and periodically varying environment. Importantly, the population size distribution is generally found to be broader under intermediate and fast random switching than under periodic variations, which results in markedly different asymptotic behaviors between the fixation probability of random and periodic switching. We also determine the detailed conditions under which the fixation probability of the slow strain is maximal.
Collapse
Affiliation(s)
- Ami Taitelbaum
- Racah Institute of Physics, Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Robert West
- Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds LS2 9JT, United Kingdom
| | - Michael Assaf
- Racah Institute of Physics, Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Mauro Mobilia
- Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds LS2 9JT, United Kingdom
| |
Collapse
|
15
|
Mittal S, Mukhopadhyay A, Chakraborty S. Evolutionary dynamics of the delayed replicator-mutator equation: Limit cycle and cooperation. Phys Rev E 2020; 101:042410. [PMID: 32422824 DOI: 10.1103/physreve.101.042410] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2020] [Accepted: 04/06/2020] [Indexed: 11/07/2022]
Abstract
Game theory deals with strategic interactions among players and evolutionary game dynamics tracks the fate of the players' populations under selection. In this paper, we consider the replicator equation for two-player-two-strategy games involving cooperators and defectors. We modify the equation to include the effect of mutation and also delay that corresponds either to the delayed information about the population state or in realizing the effect of interaction among players. By focusing on the four exhaustive classes of symmetrical games-the Stag Hunt game, the Snowdrift game, the Prisoners' Dilemma game, and the Harmony game-we analytically and numerically analyze the delayed replicator-mutator equation to find the explicit condition for the Hopf bifurcation bringing forth stable limit cycle. The existence of the asymptotically stable limit cycle imply the coexistence of the cooperators and the defectors; and in the games, where defection is a stable Nash strategy, a stable limit cycle does provide a mechanism for evolution of cooperation. We find that while mutation alone can never lead to oscillatory cooperation state in two-player-two-strategy games, the delay can change the scenario. On the other hand, there are situations when delay alone cannot lead to the Hopf bifurcation in the absence of mutation in the selection dynamics.
Collapse
Affiliation(s)
- Sourabh Mittal
- Department of Physics, Indian Institute of Technology Kanpur, Uttar Pradesh 208016, India
| | - Archan Mukhopadhyay
- Department of Physics, Indian Institute of Technology Kanpur, Uttar Pradesh 208016, India
| | - Sagar Chakraborty
- Department of Physics, Indian Institute of Technology Kanpur, Uttar Pradesh 208016, India
| |
Collapse
|
16
|
Gluzman M, Scott JG, Vladimirsky A. Optimizing adaptive cancer therapy: dynamic programming and evolutionary game theory. Proc Biol Sci 2020; 287:20192454. [PMID: 32315588 DOI: 10.1098/rspb.2019.2454] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Recent clinical trials have shown that adaptive drug therapies can be more efficient than a standard cancer treatment based on a continuous use of maximum tolerated doses (MTD). The adaptive therapy paradigm is not based on a preset schedule; instead, the doses are administered based on the current state of tumour. But the adaptive treatment policies examined so far have been largely ad hoc. We propose a method for systematically optimizing adaptive policies based on an evolutionary game theory model of cancer dynamics. Given a set of treatment objectives, we use the framework of dynamic programming to find the optimal treatment strategies. In particular, we optimize the total drug usage and time to recovery by solving a Hamilton-Jacobi-Bellman equation. We compare MTD-based treatment strategy with optimal adaptive treatment policies and show that the latter can significantly decrease the total amount of drugs prescribed while also increasing the fraction of initial tumour states from which the recovery is possible. We conclude that the use of optimal control theory to improve adaptive policies is a promising concept in cancer treatment and should be integrated into clinical trial design.
Collapse
Affiliation(s)
- Mark Gluzman
- Center for Applied Mathematics, Cornell University, Ithaca, NY, USA
| | - Jacob G Scott
- Department of Translational Hematology and Oncology Research, Cleveland Clinic, Cleveland, OH, USA
| | - Alexander Vladimirsky
- Department of Mathematics and Center for Applied Mathematics, Cornell University, 561 Malott Hall, Ithaca, NY 14853-4201, USA
| |
Collapse
|
17
|
Marrec L, Bitbol AF. Resist or perish: Fate of a microbial population subjected to a periodic presence of antimicrobial. PLoS Comput Biol 2020; 16:e1007798. [PMID: 32275712 PMCID: PMC7176291 DOI: 10.1371/journal.pcbi.1007798] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2019] [Revised: 04/22/2020] [Accepted: 03/19/2020] [Indexed: 12/22/2022] Open
Abstract
The evolution of antimicrobial resistance can be strongly affected by variations of antimicrobial concentration. Here, we study the impact of periodic alternations of absence and presence of antimicrobial on resistance evolution in a microbial population, using a stochastic model that includes variations of both population composition and size, and fully incorporates stochastic population extinctions. We show that fast alternations of presence and absence of antimicrobial are inefficient to eradicate the microbial population and strongly favor the establishment of resistance, unless the antimicrobial increases enough the death rate. We further demonstrate that if the period of alternations is longer than a threshold value, the microbial population goes extinct upon the first addition of antimicrobial, if it is not rescued by resistance. We express the probability that the population is eradicated upon the first addition of antimicrobial, assuming rare mutations. Rescue by resistance can happen either if resistant mutants preexist, or if they appear after antimicrobial is added to the environment. Importantly, the latter case is fully prevented by perfect biostatic antimicrobials that completely stop division of sensitive microorganisms. By contrast, we show that the parameter regime where treatment is efficient is larger for biocidal drugs than for biostatic drugs. This sheds light on the respective merits of different antimicrobial modes of action.
Collapse
Affiliation(s)
- Loïc Marrec
- Sorbonne Université, CNRS, Institut de Biologie Paris-Seine, Laboratoire Jean Perrin (UMR 8237), F-75005 Paris, France
| | - Anne-Florence Bitbol
- Sorbonne Université, CNRS, Institut de Biologie Paris-Seine, Laboratoire Jean Perrin (UMR 8237), F-75005 Paris, France
- Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| |
Collapse
|
18
|
Geyrhofer L, Brenner N. Coexistence and cooperation in structured habitats. BMC Ecol 2020; 20:14. [PMID: 32122337 PMCID: PMC7053132 DOI: 10.1186/s12898-020-00281-y] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Accepted: 02/18/2020] [Indexed: 12/19/2022] Open
Abstract
Background Natural habitats are typically structured, imposing constraints on inhabiting populations and their interactions. Which conditions are important for coexistence of diverse communities, and how cooperative interaction stabilizes in such populations, have been important ecological and evolutionary questions. Results We investigate a minimal ecological framework of microbial population dynamics that exhibits crucial features to show coexistence: Populations repeatedly undergo cycles of separation into compartmentalized habitats and mixing with new resources. The characteristic time-scale is longer than that typical of individual growth. Using analytic approximations, averaging techniques and phase-plane methods of dynamical systems, we provide a framework for analyzing various types of microbial interactions. Population composition and population size are both dynamic variables of the model; they are found to be decoupled both in terms of time-scale and parameter dependence. We present specific results for two examples of cooperative interaction by public goods: collective antibiotics resistance, and enhanced iron-availability by pyoverdine. We find stable coexistence to be a likely outcome. Conclusions The two simple features of a long mixing time-scale and spatial compartmentalization are enough to enable coexisting strains. In particular, costly social traits are often stabilized in such an environment—and thus cooperation established.
Collapse
Affiliation(s)
- Lukas Geyrhofer
- Network Biology Research Laboratories, and Department of Chemical Engineering, Technion-Israel Institute of Technology, Haifa, Israel.
| | - Naama Brenner
- Network Biology Research Laboratories, and Department of Chemical Engineering, Technion-Israel Institute of Technology, Haifa, Israel
| |
Collapse
|
19
|
Li X, Thirumalai D. Share, but unequally: a plausible mechanism for emergence and maintenance of intratumour heterogeneity. J R Soc Interface 2020; 16:20180820. [PMID: 30958159 DOI: 10.1098/rsif.2018.0820] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Intratumour heterogeneity (ITH), referring to the coexistence of different cell subpopulations in a single tumour, has been a major puzzle in cancer research for almost half a century. The lack of understanding of the underlying mechanism of ITH hinders progress in developing effective therapies for cancers. Based on the findings in a recent quantitative experiment on pancreatic cancer, we developed a general evolutionary model for one type of cancer, accounting for interactions between different cell populations through paracrine or juxtacrine factors. We show that the emergence of a stable heterogeneous state in a tumour requires an unequal allocation of paracrine growth factors (public goods) between cells that produce them and those that merely consume them. Our model provides a quantitative explanation of recent in vitro experimental studies in pancreatic cancer in which insulin-like growth factor II (IGF-II) plays the role of public goods. The calculated phase diagrams as a function of exogenous resources and fraction of growth factor producing cells show ITH persists only in a narrow range of concentration of exogenous IGF-II. Remarkably, maintenance of ITH requires cooperation among tumour cell subpopulations in harsh conditions, specified by lack of exogenous IGF-II, whereas surplus exogenous IGF-II elicits competition. Our theory also quantitatively accounts for measured in vivo tumour growth in glioblastoma multiforme (GBM). The predictions for GBM tumour growth as a function of the fraction of tumour cells are amenable to experimental tests. The mechanism for ITH also provides hints for devising efficacious therapies.
Collapse
Affiliation(s)
- Xin Li
- Department of Chemistry, University of Texas at Austin , Austin, TX 78712 , USA
| | - D Thirumalai
- Department of Chemistry, University of Texas at Austin , Austin, TX 78712 , USA
| |
Collapse
|
20
|
Frey F, Ziebert F, Schwarz US. Dynamics of particle uptake at cell membranes. Phys Rev E 2019; 100:052403. [PMID: 31869989 DOI: 10.1103/physreve.100.052403] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2019] [Indexed: 11/07/2022]
Abstract
Receptor-mediated endocytosis requires that the energy of adhesion overcomes the deformation energy of the plasma membrane. The resulting driving force is balanced by dissipative forces, leading to deterministic dynamical equations. While the shape of the free membrane does not play an important role for tensed and loose membranes, in the intermediate regime it leads to an important energy barrier. Here we show that this barrier is similar to but different from an effective line tension and suggest a simple analytical approximation for it. We then explore the rich dynamics of uptake for particles of different shapes and present the corresponding dynamical state diagrams. We also extend our model to include stochastic fluctuations, which facilitate uptake and lead to corresponding changes in the phase diagrams.
Collapse
Affiliation(s)
- Felix Frey
- Institute for Theoretical Physics, Heidelberg University, Philosophenweg 19, 69120 Heidelberg, Germany and BioQuant, Heidelberg University, Im Neuenheimer Feld 267, 69120 Heidelberg, Germany
| | - Falko Ziebert
- Institute for Theoretical Physics, Heidelberg University, Philosophenweg 19, 69120 Heidelberg, Germany and BioQuant, Heidelberg University, Im Neuenheimer Feld 267, 69120 Heidelberg, Germany
| | - Ulrich S Schwarz
- Institute for Theoretical Physics, Heidelberg University, Philosophenweg 19, 69120 Heidelberg, Germany and BioQuant, Heidelberg University, Im Neuenheimer Feld 267, 69120 Heidelberg, Germany
| |
Collapse
|
21
|
Wu Y, Zhang Z, Yan M, Zhang S. Environmental feedback promotes the evolution of cooperation in the structured populations. CHAOS (WOODBURY, N.Y.) 2019; 29:113101. [PMID: 31779368 DOI: 10.1063/1.5120049] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/16/2019] [Accepted: 10/15/2019] [Indexed: 06/10/2023]
Abstract
Environment plays a vital role in individual decision-making. In the game process, employing the strategy of the opponent who behaves better is nontrivial for the evolution and maintenance of cooperation, in that such a behavior may assist the player in achieving higher personal interests and more competitive superiorities. Enlightened by this thought, a coevolutionary model where the mechanisms of dynamic environment and preference selection are introduced in the networked prisoner's dilemma game is considered. Individual preference selection is introduced in the strategy update process to probe how the preferences of the latent strategy sources interfere with individual behaviors. The environment defined in the model is not only heterogeneous, but also evolves with the evolution of individual strategies. Through conducting large-scale Monte Carlo simulations, we draw a conclusion that the introduction of evolutionary environment-related preference selection is an effective promoter of cooperation even under a severe temptation. Our exploration indicates that the coevolutionary model may have a practical significance and provide a new insight into the analysis of the origin of cooperation in structured populations for further research.
Collapse
Affiliation(s)
- Yu'e Wu
- Coordinated Innovation Center for Computable Modeling in Management Science, Tianjin University of Finance and Economics, Tianjin 300222, China
| | - Zhipeng Zhang
- School of Economics and Management, Tianjin University of Technology and Education, Tianjin 300222, China
| | - Ming Yan
- Coordinated Innovation Center for Computable Modeling in Management Science, Tianjin University of Finance and Economics, Tianjin 300222, China
| | - Shuhua Zhang
- Coordinated Innovation Center for Computable Modeling in Management Science, Tianjin University of Finance and Economics, Tianjin 300222, China
| |
Collapse
|
22
|
Cremer J, Melbinger A, Wienand K, Henriquez T, Jung H, Frey E. Cooperation in Microbial Populations: Theory and Experimental Model Systems. J Mol Biol 2019; 431:4599-4644. [PMID: 31634468 DOI: 10.1016/j.jmb.2019.09.023] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2019] [Revised: 09/25/2019] [Accepted: 09/26/2019] [Indexed: 01/07/2023]
Abstract
Cooperative behavior, the costly provision of benefits to others, is common across all domains of life. This review article discusses cooperative behavior in the microbial world, mediated by the exchange of extracellular products called public goods. We focus on model species for which the production of a public good and the related growth disadvantage for the producing cells are well described. To unveil the biological and ecological factors promoting the emergence and stability of cooperative traits we take an interdisciplinary perspective and review insights gained from both mathematical models and well-controlled experimental model systems. Ecologically, we include crucial aspects of the microbial life cycle into our analysis and particularly consider population structures where ensembles of local communities (subpopulations) continuously emerge, grow, and disappear again. Biologically, we explicitly consider the synthesis and regulation of public good production. The discussion of the theoretical approaches includes general evolutionary concepts, population dynamics, and evolutionary game theory. As a specific but generic biological example, we consider populations of Pseudomonas putida and its regulation and use of pyoverdines, iron scavenging molecules, as public goods. The review closes with an overview on cooperation in spatially extended systems and also provides a critical assessment of the insights gained from the experimental and theoretical studies discussed. Current challenges and important new research opportunities are discussed, including the biochemical regulation of public goods, more realistic ecological scenarios resembling native environments, cell-to-cell signaling, and multispecies communities.
Collapse
Affiliation(s)
- J Cremer
- Department of Molecular Immunology and Microbiology, Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, 9747 AG Groningen, the Netherlands
| | - A Melbinger
- Arnold-Sommerfeld-Center for Theoretical Physics and Center for Nanoscience, Ludwig-Maximilians-Universität München, Theresienstrasse 37, D-80333 Munich, Germany
| | - K Wienand
- Arnold-Sommerfeld-Center for Theoretical Physics and Center for Nanoscience, Ludwig-Maximilians-Universität München, Theresienstrasse 37, D-80333 Munich, Germany
| | - T Henriquez
- Microbiology, Department of Biology I, Ludwig-Maximilians-Universität München, Grosshaderner Strasse 2-4, Martinsried, Germany
| | - H Jung
- Microbiology, Department of Biology I, Ludwig-Maximilians-Universität München, Grosshaderner Strasse 2-4, Martinsried, Germany.
| | - E Frey
- Arnold-Sommerfeld-Center for Theoretical Physics and Center for Nanoscience, Ludwig-Maximilians-Universität München, Theresienstrasse 37, D-80333 Munich, Germany.
| |
Collapse
|
23
|
Eliazar I, Metzler R, Reuveni S. Gumbel central limit theorem for max-min and min-max. Phys Rev E 2019; 100:020104. [PMID: 31574728 DOI: 10.1103/physreve.100.020104] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2018] [Indexed: 11/07/2022]
Abstract
The max-min and min-max of matrices arise prevalently in science and engineering. However, in many real-world situations the computation of the max-min and min-max is challenging as matrices are large and full information about their entries is lacking. Here we take a statistical-physics approach and establish limit laws-akin to the central limit theorem-for the max-min and min-max of large random matrices. The limit laws intertwine random-matrix theory and extreme-value theory, couple the matrix dimensions geometrically, and assert that Gumbel statistics emerge irrespective of the matrix entries' distribution. Due to their generality and universality, as well as their practicality, these results are expected to have a host of applications in the physical sciences and beyond.
Collapse
Affiliation(s)
- Iddo Eliazar
- School of Chemistry, The Center for Physics and Chemistry of Living Systems, The Raymond and Beverly Sackler Center for Computational Molecular and Materials Science, and The Mark Ratner Institute for Single Molecule Chemistry, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Ralf Metzler
- University of Potsdam, Institute of Physics and Astronomy, 14476 Potsdam, Germany
| | - Shlomi Reuveni
- School of Chemistry, The Center for Physics and Chemistry of Living Systems, The Raymond and Beverly Sackler Center for Computational Molecular and Materials Science, and The Mark Ratner Institute for Single Molecule Chemistry, Tel Aviv University, Tel Aviv 6997801, Israel
| |
Collapse
|
24
|
Kimmel GJ, Gerlee P, Brown JS, Altrock PM. Neighborhood size-effects shape growing population dynamics in evolutionary public goods games. Commun Biol 2019; 2:53. [PMID: 30729189 PMCID: PMC6363775 DOI: 10.1038/s42003-019-0299-4] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2018] [Accepted: 01/08/2019] [Indexed: 01/19/2023] Open
Abstract
An evolutionary game emerges when a subset of individuals incur costs to provide benefits to all individuals. Public goods games (PGG) cover the essence of such dilemmas in which cooperators are prone to exploitation by defectors. We model the population dynamics of a non-linear PGG and consider density-dependence on the global level, while the game occurs within local neighborhoods. At low cooperation, increases in the public good provide increasing returns. At high cooperation, increases provide diminishing returns. This mechanism leads to diverse evolutionarily stable strategies, including monomorphic and polymorphic populations, and neighborhood-size-driven state changes, resulting in hysteresis between equilibria. Stochastic or strategy-dependent variations in neighborhood sizes favor coexistence by destabilizing monomorphic states. We integrate our model with experiments of cancer cell growth and confirm that our framework describes PGG dynamics observed in cellular populations. Our findings advance the understanding of how neighborhood-size effects in PGG shape the dynamics of growing populations.
Collapse
Affiliation(s)
- Gregory J. Kimmel
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33629 USA
| | - Philip Gerlee
- Department of Mathematical Sciences, Chalmers University of Technology, Gothenburg, SE-412 96 Sweden
- Department of Mathematical Sciences, University of Gothenburg, Gothenburg, SE-412 61 Sweden
| | - Joel S. Brown
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33629 USA
| | - Philipp M. Altrock
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33629 USA
| |
Collapse
|
25
|
Czuppon P, Traulsen A. Fixation probabilities in populations under demographic fluctuations. J Math Biol 2018; 77:1233-1277. [PMID: 29882011 PMCID: PMC6153673 DOI: 10.1007/s00285-018-1251-9] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2017] [Revised: 05/08/2018] [Indexed: 01/09/2023]
Abstract
We study the fixation probability of a mutant type when introduced into a resident population. We implement a stochastic competitive Lotka-Volterra model with two types and intra- and interspecific competition. The model further allows for stochastically varying population sizes. The competition coefficients are interpreted in terms of inverse payoffs emerging from an evolutionary game. Since our study focuses on the impact of the competition values, we assume the same net growth rate for both types. In this general framework, we derive a formula for the fixation probability [Formula: see text] of the mutant type under weak selection. We find that the most important parameter deciding over the invasion success of the mutant is its death rate due to competition with the resident. Furthermore, we compare our approximation to results obtained by implementing population size changes deterministically in order to explore the parameter regime of validity of our method. Finally, we put our formula in the context of classical evolutionary game theory and observe similarities and differences to the results obtained in that constant population size setting.
Collapse
Affiliation(s)
- Peter Czuppon
- Department of Evolutionary Theory, Max-Planck Institute for Evolutionary Biology, Plön, Germany
| | - Arne Traulsen
- Department of Evolutionary Theory, Max-Planck Institute for Evolutionary Biology, Plön, Germany
| |
Collapse
|
26
|
Aguadé‐Gorgorió G, Solé R. Adaptive dynamics of unstable cancer populations: The canonical equation. Evol Appl 2018; 11:1283-1292. [PMID: 30151040 PMCID: PMC6099832 DOI: 10.1111/eva.12625] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2017] [Accepted: 02/15/2018] [Indexed: 12/24/2022] Open
Abstract
In most instances of tumour development, genetic instability plays a role in allowing cancer cell populations to respond to selection barriers, such as physical constraints or immune responses, and rapidly adapt to an always changing environment. Modelling instability is a nontrivial task, since by definition evolving instability leads to changes in the underlying landscape. In this article, we explore mathematically a simple version of unstable tumour progression using the formalism of adaptive dynamics (AD) where selection and mutation are explicitly coupled. Using a set of basic fitness landscapes, the so-called canonical equation for the evolution of genetic instability on a minimal scenario associated with a population of unstable cells is derived. We obtain explicit expressions for the evolution of mutation probabilities, and the implications of the model on further experimental studies and potential mutagenic therapies are discussed.
Collapse
Affiliation(s)
- Guim Aguadé‐Gorgorió
- ICREA‐Complex Systems LabUniversitat Pompeu FabraBarcelonaSpain
- Institut de Biologia Evolutiva (CSIC‐UPF)BarcelonaSpain
| | - Ricard Solé
- ICREA‐Complex Systems LabUniversitat Pompeu FabraBarcelonaSpain
- Institut de Biologia Evolutiva (CSIC‐UPF)BarcelonaSpain
- Santa Fe InstituteSanta FeNMUSA
| |
Collapse
|
27
|
Quantifying the impact of a periodic presence of antimicrobial on resistance evolution in a homogeneous microbial population of fixed size. J Theor Biol 2018; 457:190-198. [PMID: 30172688 DOI: 10.1016/j.jtbi.2018.08.040] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2018] [Revised: 07/21/2018] [Accepted: 08/30/2018] [Indexed: 01/12/2023]
Abstract
The evolution of antimicrobial resistance generally occurs in an environment where antimicrobial concentration is variable, which has dramatic consequences on the microorganisms' fitness landscape, and thus on the evolution of resistance. We investigate the effect of these time-varying patterns of selection within a stochastic model. We consider a homogeneous microbial population of fixed size subjected to periodic alternations of phases of absence and presence of an antimicrobial that stops growth. Combining analytical approaches and stochastic simulations, we quantify how the time necessary for fit resistant bacteria to take over the microbial population depends on the alternation period. We demonstrate that fast alternations strongly accelerate the evolution of resistance, reaching a plateau for sufficiently small periods. Furthermore, this acceleration is stronger in larger populations. For asymmetric alternations, featuring a different duration of the phases with and without antimicrobial, we shed light on the existence of a minimum for the time taken by the population to fully evolve resistance. The corresponding dramatic acceleration of the evolution of antimicrobial resistance likely occurs in realistic situations, and may have an important impact both in clinical and experimental situations.
Collapse
|
28
|
Wienand K, Frey E, Mobilia M. Eco-evolutionary dynamics of a population with randomly switching carrying capacity. J R Soc Interface 2018; 15:20180343. [PMID: 30135263 PMCID: PMC6127162 DOI: 10.1098/rsif.2018.0343] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2018] [Accepted: 07/26/2018] [Indexed: 01/19/2023] Open
Abstract
Environmental variability greatly influences the eco-evolutionary dynamics of a population, i.e. it affects how its size and composition evolve. Here, we study a well-mixed population of finite and fluctuating size whose growth is limited by a randomly switching carrying capacity. This models the environmental fluctuations between states of resources abundance and scarcity. The population consists of two strains, one growing slightly faster than the other, competing under two scenarios: one in which competition is solely for resources, and one in which the slow (cooperating) strain produces a public good (PG) that benefits also the fast (free-riding) strain. We investigate how the coupling of demographic and environmental (external) noise affects the population's eco-evolutionary dynamics. By analytical and computational means, we study the correlations between the population size and its composition, and discuss the social-dilemma-like 'eco-evolutionary game' characterizing the PG production. We determine in what conditions it is best to produce a PG; when cooperating is beneficial but outcompeted by free riding, and when the PG production is detrimental for cooperators. Within a linear noise approximation to populations of varying size, we also accurately analyse the coupled effects of demographic and environmental noise on the size distribution.
Collapse
Affiliation(s)
- Karl Wienand
- Arnold Sommerfeld Center for Theoretical Physics, Department of Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany
| | - Erwin Frey
- Arnold Sommerfeld Center for Theoretical Physics, Department of Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany
| | - Mauro Mobilia
- Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds LS2 9JT, UK
| |
Collapse
|
29
|
Heins AL, Weuster-Botz D. Population heterogeneity in microbial bioprocesses: origin, analysis, mechanisms, and future perspectives. Bioprocess Biosyst Eng 2018. [PMID: 29541890 DOI: 10.1007/s00449-018-1922-3] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Population heterogeneity is omnipresent in all bioprocesses even in homogenous environments. Its origin, however, is only so well understood that potential strategies like bet-hedging, noise in gene expression and division of labour that lead to population heterogeneity can be derived from experimental studies simulating the dynamics in industrial scale bioprocesses. This review aims at summarizing the current state of the different parts of single cell studies in bioprocesses. This includes setups to visualize different phenotypes of single cells, computational approaches connecting single cell physiology with environmental influence and special cultivation setups like scale-down reactors that have been proven to be useful to simulate large-scale conditions. A step in between investigation of populations and single cells is studying subpopulations with distinct properties that differ from the rest of the population with sub-omics methods which are also presented here. Moreover, the current knowledge about population heterogeneity in bioprocesses is summarized for relevant industrial production hosts and mixed cultures, as they provide the unique opportunity to distribute metabolic burden and optimize production processes in a way that is impossible in traditional monocultures. In the end, approaches to explain the underlying mechanism of population heterogeneity and the evidences found to support each hypothesis are presented. For instance, population heterogeneity serving as a bet-hedging strategy that is used as coordinated action against bioprocess-related stresses while at the same time spreading the risk between individual cells as it ensures the survival of least a part of the population in any environment the cells encounter.
Collapse
Affiliation(s)
- Anna-Lena Heins
- Institute of Biochemical Engineering, Technical University of Munich, Boltzmannstr. 15, 85748, Garching, Germany.
| | - Dirk Weuster-Botz
- Institute of Biochemical Engineering, Technical University of Munich, Boltzmannstr. 15, 85748, Garching, Germany
| |
Collapse
|
30
|
Becker F, Wienand K, Lechner M, Frey E, Jung H. Interactions mediated by a public good transiently increase cooperativity in growing Pseudomonas putida metapopulations. Sci Rep 2018; 8:4093. [PMID: 29511247 PMCID: PMC5840296 DOI: 10.1038/s41598-018-22306-9] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2017] [Accepted: 02/21/2018] [Indexed: 01/13/2023] Open
Abstract
Bacterial communities have rich social lives. A well-established interaction involves the exchange of a public good in Pseudomonas populations, where the iron-scavenging compound pyoverdine, synthesized by some cells, is shared with the rest. Pyoverdine thus mediates interactions between producers and non-producers and can constitute a public good. This interaction is often used to test game theoretical predictions on the "social dilemma" of producers. Such an approach, however, underestimates the impact of specific properties of the public good, for example consequences of its accumulation in the environment. Here, we experimentally quantify costs and benefits of pyoverdine production in a specific environment, and build a model of population dynamics that explicitly accounts for the changing significance of accumulating pyoverdine as chemical mediator of social interactions. The model predicts that, in an ensemble of growing populations (metapopulation) with different initial producer fractions (and consequently pyoverdine contents), the global producer fraction initially increases. Because the benefit of pyoverdine declines at saturating concentrations, the increase need only be transient. Confirmed by experiments on metapopulations, our results show how a changing benefit of a public good can shape social interactions in a bacterial population.
Collapse
Affiliation(s)
- Felix Becker
- Microbiology, Department Biology 1, Ludwig-Maximilians-Universität Munich, Grosshaderner Strasse 2-4, D-82152 Martinsried, Germany
| | - Karl Wienand
- Arnold-Sommerfeld-Center for Theoretical Physics and Center for Nanoscience, Ludwig-Maximilians-Universität, Theresienstrasse 37, D-80333, Munich, Germany
| | - Matthias Lechner
- Arnold-Sommerfeld-Center for Theoretical Physics and Center for Nanoscience, Ludwig-Maximilians-Universität, Theresienstrasse 37, D-80333, Munich, Germany
| | - Erwin Frey
- Arnold-Sommerfeld-Center for Theoretical Physics and Center for Nanoscience, Ludwig-Maximilians-Universität, Theresienstrasse 37, D-80333, Munich, Germany.
| | - Heinrich Jung
- Microbiology, Department Biology 1, Ludwig-Maximilians-Universität Munich, Grosshaderner Strasse 2-4, D-82152 Martinsried, Germany.
| |
Collapse
|
31
|
McAvoy A, Fraiman N, Hauert C, Wakeley J, Nowak MA. Public goods games in populations with fluctuating size. Theor Popul Biol 2018; 121:72-84. [PMID: 29408219 DOI: 10.1016/j.tpb.2018.01.004] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2017] [Revised: 01/16/2018] [Accepted: 01/24/2018] [Indexed: 11/18/2022]
Abstract
Many mathematical frameworks of evolutionary game dynamics assume that the total population size is constant and that selection affects only the relative frequency of strategies. Here, we consider evolutionary game dynamics in an extended Wright-Fisher process with variable population size. In such a scenario, it is possible that the entire population becomes extinct. Survival of the population may depend on which strategy prevails in the game dynamics. Studying cooperative dilemmas, it is a natural feature of such a model that cooperators enable survival, while defectors drive extinction. Although defectors are favored for any mixed population, random drift could lead to their elimination and the resulting pure-cooperator population could survive. On the other hand, if the defectors remain, then the population will quickly go extinct because the frequency of cooperators steadily declines and defectors alone cannot survive. In a mutation-selection model, we find that (i) a steady supply of cooperators can enable long-term population survival, provided selection is sufficiently strong, and (ii) selection can increase the abundance of cooperators but reduce their relative frequency. Thus, evolutionary game dynamics in populations with variable size generate a multifaceted notion of what constitutes a trait's long-term success.
Collapse
Affiliation(s)
- Alex McAvoy
- Program for Evolutionary Dynamics, Harvard University, Cambridge, MA 02138, United States.
| | - Nicolas Fraiman
- Department of Statistics and Operations Research, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, United States
| | - Christoph Hauert
- Department of Mathematics, University of British Columbia, 1984 Mathematics Road, Vancouver, BC, Canada V6T 1Z2
| | - John Wakeley
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, United States
| | - Martin A Nowak
- Program for Evolutionary Dynamics, Harvard University, Cambridge, MA 02138, United States; Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, United States; Department of Mathematics, Harvard University, Cambridge, MA 02138, United States
| |
Collapse
|
32
|
Wienand K, Frey E, Mobilia M. Evolution of a Fluctuating Population in a Randomly Switching Environment. PHYSICAL REVIEW LETTERS 2017; 119:158301. [PMID: 29077432 DOI: 10.1103/physrevlett.119.158301] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/16/2017] [Indexed: 06/07/2023]
Abstract
Environment plays a fundamental role in the competition for resources, and hence in the evolution of populations. Here, we study a well-mixed, finite population consisting of two strains competing for the limited resources provided by an environment that randomly switches between states of abundance and scarcity. Assuming that one strain grows slightly faster than the other, we consider two scenarios-one of pure resource competition, and one in which one strain provides a public good-and investigate how environmental randomness (external noise) coupled to demographic (internal) noise determines the population's fixation properties and size distribution. By analytical means and simulations, we show that these coupled sources of noise can significantly enhance the fixation probability of the slower-growing species. We also show that the population size distribution can be unimodal, bimodal, or multimodal and undergoes noise-induced transitions between these regimes when the rate of switching matches the population's growth rate.
Collapse
Affiliation(s)
- Karl Wienand
- Arnold Sommerfeld Center for Theoretical Physics, Department of Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany
| | - Erwin Frey
- Arnold Sommerfeld Center for Theoretical Physics, Department of Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany
| | - Mauro Mobilia
- Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds LS2 9JT, United Kingdom
| |
Collapse
|
33
|
Bauer M, Knebel J, Lechner M, Pickl P, Frey E. Ecological feedback in quorum-sensing microbial populations can induce heterogeneous production of autoinducers. eLife 2017; 6. [PMID: 28741470 PMCID: PMC5526673 DOI: 10.7554/elife.25773] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2017] [Accepted: 06/15/2017] [Indexed: 12/11/2022] Open
Abstract
Autoinducers are small signaling molecules that mediate intercellular communication in microbial populations and trigger coordinated gene expression via ‘quorum sensing’. Elucidating the mechanisms that control autoinducer production is, thus, pertinent to understanding collective microbial behavior, such as virulence and bioluminescence. Recent experiments have shown a heterogeneous promoter activity of autoinducer synthase genes, suggesting that some of the isogenic cells in a population might produce autoinducers, whereas others might not. However, the mechanism underlying this phenotypic heterogeneity in quorum-sensing microbial populations has remained elusive. In our theoretical model, cells synthesize and secrete autoinducers into the environment, up-regulate their production in this self-shaped environment, and non-producers replicate faster than producers. We show that the coupling between ecological and population dynamics through quorum sensing can induce phenotypic heterogeneity in microbial populations, suggesting an alternative mechanism to stochastic gene expression in bistable gene regulatory circuits. DOI:http://dx.doi.org/10.7554/eLife.25773.001 Bacteria and other microbes can communicate with each other using chemical languages. They release small signaling molecules called autoinducers into their surroundings and sense the levels of the autoinducers in the environment. The response to these autoinducers – known as quorum sensing – can regulate how whole communities of microbes grow and behave; for example, autoinducers can alter the ability of microbes to infect humans or enable the microbes to collectively switch on light production. Recent experiments suggest that, in a population of genetically identical microbes, some individuals may produce autoinducers while others do not. The coexistence of these different “phenotypes” in one population may enable different individuals to perform different roles, or act as a “bet-hedging” strategy that helps the population to survive if it is later exposed to a stressful situation. It is not clear how microbes regulate autoinducer production so that only some individuals produce these molecules. Bauer, Knebel et al. developed a theoretical model to address this question. In the model, the microbes shape their environment by producing autoinducers and can respond to this self-shaped environment by changing their level of autoinducer production. Bauer, Knebel et al. found that this establishes a feedback loop that can result in autoinducers being produced by some individuals but not others. The next step following on from this work is to carry out experiments to test the assumptions and predictions made by the theoretical model. These findings may help to understand how the coexistence of different phenotypes affects collective behaviors, and vice versa, in populations of microbes that use quorum-sensing. DOI:http://dx.doi.org/10.7554/eLife.25773.002
Collapse
Affiliation(s)
- Matthias Bauer
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, München, Germany
| | - Johannes Knebel
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, München, Germany
| | - Matthias Lechner
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, München, Germany
| | - Peter Pickl
- Department of Mathematics, Ludwig-Maximilians-Universität München, München, Germany
| | - Erwin Frey
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, München, Germany
| |
Collapse
|
34
|
Ashcroft P, Smith CE, Garrod M, Galla T. Effects of population growth on the success of invading mutants. J Theor Biol 2017; 420:232-240. [DOI: 10.1016/j.jtbi.2017.03.014] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2016] [Revised: 02/28/2017] [Accepted: 03/15/2017] [Indexed: 10/19/2022]
|
35
|
Weber MF, Frey E. Master equations and the theory of stochastic path integrals. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2017; 80:046601. [PMID: 28306551 DOI: 10.1088/1361-6633/aa5ae2] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
This review provides a pedagogic and self-contained introduction to master equations and to their representation by path integrals. Since the 1930s, master equations have served as a fundamental tool to understand the role of fluctuations in complex biological, chemical, and physical systems. Despite their simple appearance, analyses of master equations most often rely on low-noise approximations such as the Kramers-Moyal or the system size expansion, or require ad-hoc closure schemes for the derivation of low-order moment equations. We focus on numerical and analytical methods going beyond the low-noise limit and provide a unified framework for the study of master equations. After deriving the forward and backward master equations from the Chapman-Kolmogorov equation, we show how the two master equations can be cast into either of four linear partial differential equations (PDEs). Three of these PDEs are discussed in detail. The first PDE governs the time evolution of a generalized probability generating function whose basis depends on the stochastic process under consideration. Spectral methods, WKB approximations, and a variational approach have been proposed for the analysis of the PDE. The second PDE is novel and is obeyed by a distribution that is marginalized over an initial state. It proves useful for the computation of mean extinction times. The third PDE describes the time evolution of a 'generating functional', which generalizes the so-called Poisson representation. Subsequently, the solutions of the PDEs are expressed in terms of two path integrals: a 'forward' and a 'backward' path integral. Combined with inverse transformations, one obtains two distinct path integral representations of the conditional probability distribution solving the master equations. We exemplify both path integrals in analysing elementary chemical reactions. Moreover, we show how a well-known path integral representation of averaged observables can be recovered from them. Upon expanding the forward and the backward path integrals around stationary paths, we then discuss and extend a recent method for the computation of rare event probabilities. Besides, we also derive path integral representations for processes with continuous state spaces whose forward and backward master equations admit Kramers-Moyal expansions. A truncation of the backward expansion at the level of a diffusion approximation recovers a classic path integral representation of the (backward) Fokker-Planck equation. One can rewrite this path integral in terms of an Onsager-Machlup function and, for purely diffusive Brownian motion, it simplifies to the path integral of Wiener. To make this review accessible to a broad community, we have used the language of probability theory rather than quantum (field) theory and do not assume any knowledge of the latter. The probabilistic structures underpinning various technical concepts, such as coherent states, the Doi-shift, and normal-ordered observables, are thereby made explicit.
Collapse
Affiliation(s)
- Markus F Weber
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, Theresienstraße 37, 80333 München, Germany
| | - Erwin Frey
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, Theresienstraße 37, 80333 München, Germany
| |
Collapse
|
36
|
Daly AJ, Baetens JM, De Baets B. The impact of resource dependence of the mechanisms of life on the spatial population dynamics of an in silico microbial community. CHAOS (WOODBURY, N.Y.) 2016; 26:123121. [PMID: 28039986 DOI: 10.1063/1.4972788] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Biodiversity has a critical impact on ecosystem functionality and stability, and thus the current biodiversity crisis has motivated many studies of the mechanisms that sustain biodiversity, a notable example being non-transitive or cyclic competition. We therefore extend existing microscopic models of communities with cyclic competition by incorporating resource dependence in demographic processes, characteristics of natural systems often oversimplified or overlooked by modellers. The spatially explicit nature of our individual-based model of three interacting species results in the formation of stable spatial structures, which have significant effects on community functioning, in agreement with experimental observations of pattern formation in microbial communities.
Collapse
Affiliation(s)
- Aisling J Daly
- KERMIT, Department of Mathematical Modelling, Statistics and Bioinformatics, Ghent University, Coupure links 653, Ghent B-9000, Belgium
| | - Jan M Baetens
- KERMIT, Department of Mathematical Modelling, Statistics and Bioinformatics, Ghent University, Coupure links 653, Ghent B-9000, Belgium
| | - Bernard De Baets
- KERMIT, Department of Mathematical Modelling, Statistics and Bioinformatics, Ghent University, Coupure links 653, Ghent B-9000, Belgium
| |
Collapse
|
37
|
|
38
|
Daly AJ, Baetens JM, De Baets B. In silico substrate dependence increases community productivity but threatens biodiversity. Phys Rev E 2016; 93:042414. [PMID: 27176336 DOI: 10.1103/physreve.93.042414] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2016] [Indexed: 11/07/2022]
Abstract
The critical role that biodiversity plays in ecosystem functioning has motivated many studies of the mechanisms that sustain biodiversity, a notable example being cyclic competition. We extend existing models of communities with cyclic competition by incorporating variable community evenness and resource dependence in demographic processes, two features that have generally been neglected. In this way, we align previous approaches more closely with real-world microbial ecosystems. We demonstrate the existence of a trade-off between increasing biomass production and maintaining biodiversity. This supports experimental observations of a net negative biodiversity effect on biomass productivity, due to competition effects suffered by highly productive species in diverse communities. Our results also support the important role assigned by microbial ecologists to evenness in maintaining ecosystem stability, thus far largely overlooked in in silico approaches.
Collapse
Affiliation(s)
- Aisling J Daly
- KERMIT, Department of Mathematical Modelling, Statistics and Bioinformatics, Ghent University, Coupure Links 653, B-9000 Ghent, Belgium
| | - Jan M Baetens
- KERMIT, Department of Mathematical Modelling, Statistics and Bioinformatics, Ghent University, Coupure Links 653, B-9000 Ghent, Belgium
| | - Bernard De Baets
- KERMIT, Department of Mathematical Modelling, Statistics and Bioinformatics, Ghent University, Coupure Links 653, B-9000 Ghent, Belgium
| |
Collapse
|
39
|
Li XY, Pietschke C, Fraune S, Altrock PM, Bosch TCG, Traulsen A. Which games are growing bacterial populations playing? J R Soc Interface 2016; 12:20150121. [PMID: 26236827 PMCID: PMC4528578 DOI: 10.1098/rsif.2015.0121] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Microbial communities display complex population dynamics, both in frequency and absolute density. Evolutionary game theory provides a natural approach to analyse and model this complexity by studying the detailed interactions among players, including competition and conflict, cooperation and coexistence. Classic evolutionary game theory models typically assume constant population size, which often does not hold for microbial populations. Here, we explicitly take into account population growth with frequency-dependent growth parameters, as observed in our experimental system. We study the in vitro population dynamics of the two commensal bacteria (Curvibacter sp. (AEP1.3) and Duganella sp. (C1.2)) that synergistically protect the metazoan host Hydra vulgaris (AEP) from fungal infection. The frequency-dependent, nonlinear growth rates observed in our experiments indicate that the interactions among bacteria in co-culture are beyond the simple case of direct competition or, equivalently, pairwise games. This is in agreement with the synergistic effect of anti-fungal activity observed in vivo. Our analysis provides new insight into the minimal degree of complexity needed to appropriately understand and predict coexistence or extinction events in this kind of microbial community dynamics. Our approach extends the understanding of microbial communities and points to novel experiments.
Collapse
Affiliation(s)
- Xiang-Yi Li
- Department of Evolutionary Theory, Max Planck Institute for Evolutionary Biology, August-Thienemannstraße 2, 24306 Plön, Germany
| | - Cleo Pietschke
- Department of Evolutionary Theory, Max Planck Institute for Evolutionary Biology, August-Thienemannstraße 2, 24306 Plön, Germany
- Zoological Institute, Christian-Albrechts-University Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| | - Sebastian Fraune
- Zoological Institute, Christian-Albrechts-University Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| | - Philipp M. Altrock
- Dana-Farber Cancer Institute, Boston, MA 02215, USA
- Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
- Program for Evolutionary Dynamics, Harvard University, Cambridge, MA 02138, USA
| | - Thomas C. G. Bosch
- Zoological Institute, Christian-Albrechts-University Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| | - Arne Traulsen
- Department of Evolutionary Theory, Max Planck Institute for Evolutionary Biology, August-Thienemannstraße 2, 24306 Plön, Germany
- e-mail:
| |
Collapse
|
40
|
Melbinger A, Vergassola M. The Impact of Environmental Fluctuations on Evolutionary Fitness Functions. Sci Rep 2015; 5:15211. [PMID: 26477392 PMCID: PMC4609966 DOI: 10.1038/srep15211] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2015] [Accepted: 09/21/2015] [Indexed: 01/30/2023] Open
Abstract
The concept of fitness as a measure for a species' success in natural selection is central to the theory of evolution. We here investigate how reproduction rates which are not constant but vary in response to environmental fluctuations, influence a species' prosperity and thereby its fitness. Interestingly, we find that not only larger growth rates but also reduced sensitivities to environmental changes substantially increase the fitness. Thereby, depending on the noise level of the environment, it might be an evolutionary successful strategy to minimize this sensitivity rather than to optimize the reproduction speed. Also for neutral evolution, where species with exactly the same properties compete, variability in the growth rates plays a crucial role. The time for one species to fixate is strongly reduced in the presence of environmental noise. Hence, environmental fluctuations constitute a possible explanation for effective population sizes inferred from genetic data that often are much smaller than the census population size.
Collapse
Affiliation(s)
- Anna Melbinger
- University of California, San Diego, Department of Physics, 9500 Gilman Drive, La Jolla, CA 92093, USA
| | - Massimo Vergassola
- University of California, San Diego, Department of Physics, 9500 Gilman Drive, La Jolla, CA 92093, USA
| |
Collapse
|
41
|
Wienand K, Lechner M, Becker F, Jung H, Frey E. Non-Selective Evolution of Growing Populations. PLoS One 2015; 10:e0134300. [PMID: 26274606 PMCID: PMC4537121 DOI: 10.1371/journal.pone.0134300] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2015] [Accepted: 07/07/2015] [Indexed: 11/24/2022] Open
Abstract
Non-selective effects, like genetic drift, are an important factor in modern conceptions of evolution, and have been extensively studied for constant population sizes (Kimura, 1955; Otto and Whitlock, 1997). Here, we consider non-selective evolution in the case of growing populations that are of small size and have varying trait compositions (e.g. after a population bottleneck). We find that, in these conditions, populations never fixate to a trait, but tend to a random limit composition, and that the distribution of compositions “freezes” to a steady state. This final state is crucially influenced by the initial conditions. We obtain these findings from a combined theoretical and experimental approach, using multiple mixed subpopulations of two Pseudomonas putida strains in non-selective growth conditions (Matthijs et al, 2009) as model system. The experimental results for the population dynamics match the theoretical predictions based on the Pólya urn model (Eggenberger and Pólya, 1923) for all analyzed parameter regimes. In summary, we show that exponential growth stops genetic drift. This result contrasts with previous theoretical analyses of non-selective evolution (e.g. genetic drift), which investigated how traits spread and eventually take over populations (fixate) (Kimura, 1955; Otto and Whitlock, 1997). Moreover, our work highlights how deeply growth influences non-selective evolution, and how it plays a key role in maintaining genetic variability. Consequently, it is of particular importance in life-cycles models (Melbinger et al, 2010; Cremer et al, 2011; Cremer et al, 2012) of periodically shrinking and expanding populations.
Collapse
Affiliation(s)
- Karl Wienand
- Arnold-Sommerfeld-Center for Theoretical Physics and Center for NanoScience, Physics Department, Ludwig-Maximilians-Universität, Munich, Germany
| | - Matthias Lechner
- Arnold-Sommerfeld-Center for Theoretical Physics and Center for NanoScience, Physics Department, Ludwig-Maximilians-Universität, Munich, Germany
| | - Felix Becker
- Department of Biology 1, Microbiology, Ludwig-Maximilians-Universität, Martinsried, Germany
| | - Heinrich Jung
- Department of Biology 1, Microbiology, Ludwig-Maximilians-Universität, Martinsried, Germany
| | - Erwin Frey
- Arnold-Sommerfeld-Center for Theoretical Physics and Center for NanoScience, Physics Department, Ludwig-Maximilians-Universität, Munich, Germany
- * E-mail:
| |
Collapse
|
42
|
Chotibut T, Nelson DR. Evolutionary dynamics with fluctuating population sizes and strong mutualism. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:022718. [PMID: 26382443 DOI: 10.1103/physreve.92.022718] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/18/2014] [Indexed: 06/05/2023]
Abstract
Game theory ideas provide a useful framework for studying evolutionary dynamics in a well-mixed environment. This approach, however, typically enforces a strictly fixed overall population size, deemphasizing natural growth processes. We study a competitive Lotka-Volterra model, with number fluctuations, that accounts for natural population growth and encompasses interaction scenarios typical of evolutionary games. We show that, in an appropriate limit, the model describes standard evolutionary games with both genetic drift and overall population size fluctuations. However, there are also regimes where a varying population size can strongly influence the evolutionary dynamics. We focus on the strong mutualism scenario and demonstrate that standard evolutionary game theory fails to describe our simulation results. We then analytically and numerically determine fixation probabilities as well as mean fixation times using matched asymptotic expansions, taking into account the population size degree of freedom. These results elucidate the interplay between population dynamics and evolutionary dynamics in well-mixed systems.
Collapse
Affiliation(s)
- Thiparat Chotibut
- Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA
| | - David R Nelson
- Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA
| |
Collapse
|
43
|
Melbinger A, Cremer J, Frey E. The emergence of cooperation from a single mutant during microbial life cycles. J R Soc Interface 2015; 12:20150171. [PMID: 26063816 PMCID: PMC4528582 DOI: 10.1098/rsif.2015.0171] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2015] [Accepted: 05/14/2015] [Indexed: 12/23/2022] Open
Abstract
Cooperative behaviour is widespread in nature, even though cooperating individuals always run the risk of being exploited by free-riders. Population structure effectively promotes cooperation given that a threshold in the level of cooperation was already reached. However, the question how cooperation can emerge from a single mutant, which cannot rely on a benefit provided by other cooperators, is still puzzling. Here, we investigate this question for a well-defined but generic situation based on typical life cycles of microbial populations where individuals regularly form new colonies followed by growth phases. We analyse two evolutionary mechanisms favouring cooperative behaviour and study their strength depending on the inoculation size and the length of a life cycle. In particular, we find that population bottlenecks followed by exponential growth phases strongly increase the survival and fixation probabilities of a single cooperator in a free-riding population.
Collapse
Affiliation(s)
- Anna Melbinger
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, München, Germany Department of Physics, UCSD, 9500 Gilman Drive, La Jolla, CA 92093, USA
| | - Jonas Cremer
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, München, Germany Department of Physics, UCSD, 9500 Gilman Drive, La Jolla, CA 92093, USA
| | - Erwin Frey
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, München, Germany
| |
Collapse
|
44
|
Knebel J, Weber MF, Krüger T, Frey E. Evolutionary games of condensates in coupled birth-death processes. Nat Commun 2015; 6:6977. [PMID: 25908384 PMCID: PMC4421816 DOI: 10.1038/ncomms7977] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2014] [Accepted: 03/20/2015] [Indexed: 11/09/2022] Open
Abstract
Condensation phenomena arise through a collective behaviour of particles. They are observed in both classical and quantum systems, ranging from the formation of traffic jams in mass transport models to the macroscopic occupation of the energetic ground state in ultra-cold bosonic gases (Bose-Einstein condensation). Recently, it has been shown that a driven and dissipative system of bosons may form multiple condensates. Which states become the condensates has, however, remained elusive thus far. The dynamics of this condensation are described by coupled birth-death processes, which also occur in evolutionary game theory. Here we apply concepts from evolutionary game theory to explain the formation of multiple condensates in such driven-dissipative bosonic systems. We show that the vanishing of relative entropy production determines their selection. The condensation proceeds exponentially fast, but the system never comes to rest. Instead, the occupation numbers of condensates may oscillate, as we demonstrate for a rock-paper-scissors game of condensates.
Collapse
Affiliation(s)
- Johannes Knebel
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, Theresienstraße 37, 80333 München, Germany
| | - Markus F. Weber
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, Theresienstraße 37, 80333 München, Germany
| | - Torben Krüger
- IST Austria, Am Campus 1, 3400 Klosterneuburg, Austria
| | - Erwin Frey
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, Theresienstraße 37, 80333 München, Germany
| |
Collapse
|
45
|
Zhang L, Zou Y, Guan S, Liu Z. Analytical description for the critical fixations of evolutionary coordination games on finite complex structured populations. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:042807. [PMID: 25974547 DOI: 10.1103/physreve.91.042807] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/25/2014] [Indexed: 06/04/2023]
Abstract
Evolutionary game theory is crucial to capturing the characteristic interaction patterns among selfish individuals. In a population of coordination games of two strategies, one of the central problems is to determine the fixation probability that the system reaches a state of networkwide of only one strategy, and the corresponding expectation times. The deterministic replicator equations predict the critical value of initial density of one strategy, which separates the two absorbing states of the system. However, numerical estimations of this separatrix show large deviations from the theory in finite populations. Here we provide a stochastic treatment of this dynamic process on complex networks of finite sizes as Markov processes, showing the evolutionary time explicitly. We describe analytically the effects of network structures on the intermediate fixations as observed in numerical simulations. Our theoretical predictions are validated by various simulations on both random and scale free networks. Therefore, our stochastic framework can be helpful in dealing with other networked game dynamics.
Collapse
Affiliation(s)
- Liye Zhang
- Department of Physics, East China Normal University, Shanghai 200062, China
| | - Yong Zou
- Department of Physics, East China Normal University, Shanghai 200062, China
- State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China
| | - Shuguang Guan
- Department of Physics, East China Normal University, Shanghai 200062, China
- State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China
| | - Zonghua Liu
- Department of Physics, East China Normal University, Shanghai 200062, China
- State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|
46
|
Ashcroft P, Galla T. Pattern formation in individual-based systems with time-varying parameters. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:062104. [PMID: 24483383 DOI: 10.1103/physreve.88.062104] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/28/2013] [Indexed: 06/03/2023]
Abstract
We study the patterns generated in finite-time sweeps across symmetry-breaking bifurcations in individual-based models. Similar to the well-known Kibble-Zurek scenario of defect formation, large-scale patterns are generated when model parameters are varied slowly, whereas fast sweeps produce a large number of small domains. The symmetry breaking is triggered by intrinsic noise, originating from the discrete dynamics at the microlevel. Based on a linear-noise approximation, we calculate the characteristic length scale of these patterns. We demonstrate the applicability of this approach in a simple model of opinion dynamics, a model in evolutionary game theory with a time-dependent fitness structure, and a model of cell differentiation. Our theoretical estimates are confirmed in simulations. In further numerical work, we observe a similar phenomenon when the symmetry-breaking bifurcation is triggered by population growth.
Collapse
Affiliation(s)
- Peter Ashcroft
- Theoretical Physics, School of Physics and Astronomy, The University of Manchester, Manchester M13 9PL, United Kingdom
| | - Tobias Galla
- Theoretical Physics, School of Physics and Astronomy, The University of Manchester, Manchester M13 9PL, United Kingdom
| |
Collapse
|
47
|
Schröder M, Rahbari SHE, Nagler J. Crackling noise in fractional percolation. Nat Commun 2013; 4:2222. [DOI: 10.1038/ncomms3222] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2012] [Accepted: 07/01/2013] [Indexed: 11/09/2022] Open
|
48
|
Knebel J, Krüger T, Weber MF, Frey E. Coexistence and survival in conservative Lotka-Volterra networks. PHYSICAL REVIEW LETTERS 2013; 110:168106. [PMID: 23679644 DOI: 10.1103/physrevlett.110.168106] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/14/2012] [Indexed: 06/02/2023]
Abstract
Analyzing coexistence and survival scenarios of Lotka-Volterra (LV) networks in which the total biomass is conserved is of vital importance for the characterization of long-term dynamics of ecological communities. Here, we introduce a classification scheme for coexistence scenarios in these conservative LV models and quantify the extinction process by employing the Pfaffian of the network's interaction matrix. We illustrate our findings on global stability properties for general systems of four and five species and find a generalized scaling law for the extinction time.
Collapse
Affiliation(s)
- Johannes Knebel
- Department of Physics, Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany
| | | | | | | |
Collapse
|
49
|
Gelimson A, Cremer J, Frey E. Mobility, fitness collection, and the breakdown of cooperation. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 87:042711. [PMID: 23679453 DOI: 10.1103/physreve.87.042711] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/22/2012] [Revised: 02/02/2013] [Indexed: 06/02/2023]
Abstract
The spatial arrangement of individuals is thought to overcome the dilemma of cooperation: When cooperators engage in clusters, they might share the benefit of cooperation while being more protected against noncooperating individuals, who benefit from cooperation but save the cost of cooperation. This is paradigmatically shown by the spatial prisoner's dilemma model. Here, we study this model in one and two spatial dimensions, but explicitly take into account that in biological setups, fitness collection and selection are separated processes occurring mostly on vastly different time scales. This separation is particularly important to understand the impact of mobility on the evolution of cooperation. We find that even small diffusive mobility strongly restricts cooperation since it enables noncooperative individuals to invade cooperative clusters. Thus, in most biological scenarios, where the mobility of competing individuals is an irrefutable fact, the spatial prisoner's dilemma alone cannot explain stable cooperation, but additional mechanisms are necessary for spatial structure to promote the evolution of cooperation. The breakdown of cooperation is analyzed in detail. We confirm the existence of a phase transition, here controlled by mobility and costs, which distinguishes between purely cooperative and noncooperative absorbing states. While in one dimension the model is in the class of the voter model, it belongs to the directed percolation universality class in two dimensions.
Collapse
Affiliation(s)
- Anatolij Gelimson
- Arnold Sommerfeld Center for Theoretical Physics and Center for NanoScience, Department of Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, D-80333 München, Germany
| | | | | |
Collapse
|
50
|
Requejo RJ, Camacho J. Scarcity may promote cooperation in populations of simple agents. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 87:022819. [PMID: 23496580 DOI: 10.1103/physreve.87.022819] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/28/2012] [Revised: 05/22/2012] [Indexed: 06/01/2023]
Abstract
In the study of the evolution of cooperation, resource limitations are usually assumed just to provide a finite population size. Recently, however, it has been pointed out that resource limitation may also generate dynamical payoffs able to modify the original structure of the games. Here we study analytically a phase transition from a homogeneous population of defectors when resources are abundant to the survival of unconditional cooperators when resources reduce below a threshold. To this end, we introduce a model of simple agents, with no memory or ability of recognition, interacting in well-mixed populations. The result might shed light on the role played by resource constraints on the origin of multicellularity.
Collapse
Affiliation(s)
- R J Requejo
- Departament de Física, Universitat Autònoma de Barcelona, Campus UAB, E-08193 Bellaterra, Spain
| | | |
Collapse
|