1
|
Cirne D, Campos PRA. Rate of environmental variation impacts the predictability in evolution. Phys Rev E 2022; 106:064408. [PMID: 36671169 DOI: 10.1103/physreve.106.064408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2022] [Accepted: 12/09/2022] [Indexed: 12/24/2022]
Abstract
In the two last decades, we have improved our understanding of the adaptive evolution of natural populations under constant and stable environments. For instance, experimental methods from evolutionary biology have allowed us to explore the structure of fitness landscapes and survey how the landscape properties can constrain the adaptation process. However, understanding how environmental changes can affect adaptation remains challenging. Very little progress has been made with respect to time-varying fitness landscapes. Using the adaptive-walk approximation, we survey the evolutionary process of populations under a scenario of environmental variation. In particular, we investigate how the rate of environmental variation influences the predictability in evolution. We observe that the rate of environmental variation not only changes the duration of adaptive walks towards fitness peaks of the fitness landscape, but also affects the degree of repeatability of both outcomes and evolutionary paths. In general, slower environmental variation increases the predictability in evolution. The accessibility of endpoints is greatly influenced by the ecological dynamics. The dependence of these quantities on the genome size and number of traits is also addressed. To our knowledge, this contribution is the first to use the predictive approach to quantify and understand the impact of the speed of environmental variation on the degree of parallelism of the evolutionary process.
Collapse
Affiliation(s)
- Diego Cirne
- Departamento de Física, Universidade Federal de Pernambuco, 50740-560 Recife-PE, Brazil
| | - Paulo R A Campos
- Departamento de Física, Universidade Federal de Pernambuco, 50740-560 Recife-PE, Brazil
| |
Collapse
|
2
|
Sidhom L, Galla T. Ecological communities from random generalized Lotka-Volterra dynamics with nonlinear feedback. Phys Rev E 2020; 101:032101. [PMID: 32289927 DOI: 10.1103/physreve.101.032101] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2019] [Accepted: 01/28/2020] [Indexed: 12/11/2022]
Abstract
We investigate the outcome of generalized Lotka-Volterra dynamics of ecological communities with random interaction coefficients and nonlinear feedback. We show in simulations that the saturation of nonlinear feedback stabilizes the dynamics. This is confirmed in an analytical generating-functional approach to generalized Lotka-Volterra equations with piecewise linear saturating response. For such systems we are able to derive self-consistent relations governing the stable fixed-point phase and to carry out a linear stability analysis to predict the onset of unstable behavior. We investigate in detail the combined effects of the mean, variance, and covariance of the random interaction coefficients, and the saturation value of the nonlinear response. We find that stability and diversity increases with the introduction of nonlinear feedback, where decreasing the saturation value has a similar effect to decreasing the covariance. We also find cooperation to no longer have a detrimental effect on stability with nonlinear feedback, and the order parameters mean abundance and diversity to be less dependent on the symmetry of interactions with stronger saturation.
Collapse
Affiliation(s)
- Laura Sidhom
- Theoretical Physics, Department of Physics and Astronomy, School of Natural Sciences, The University of Manchester, Manchester M13 9PL, United Kingdom
| | - Tobias Galla
- Theoretical Physics, Department of Physics and Astronomy, School of Natural Sciences, The University of Manchester, Manchester M13 9PL, United Kingdom and Instituto de Física Interdisciplinar y Sistemas Complejos, IFISC (CSIC-UIB), Campus Universitat Illes Balears, E-07122 Palma de Mallorca, Spain
| |
Collapse
|
3
|
Galla T. Dynamically evolved community size and stability of random Lotka-Volterra ecosystems
(a). ACTA ACUST UNITED AC 2018. [DOI: 10.1209/0295-5075/123/48004] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
4
|
The prevalence of chaotic dynamics in games with many players. Sci Rep 2018; 8:4902. [PMID: 29559641 PMCID: PMC5861132 DOI: 10.1038/s41598-018-22013-5] [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: 04/07/2017] [Accepted: 01/23/2018] [Indexed: 11/23/2022] Open
Abstract
We study adaptive learning in a typical p-player game. The payoffs of the games are randomly generated and then held fixed. The strategies of the players evolve through time as the players learn. The trajectories in the strategy space display a range of qualitatively different behaviours, with attractors that include unique fixed points, multiple fixed points, limit cycles and chaos. In the limit where the game is complicated, in the sense that the players can take many possible actions, we use a generating-functional approach to establish the parameter range in which learning dynamics converge to a stable fixed point. The size of this region goes to zero as the number of players goes to infinity, suggesting that complex non-equilibrium behaviour, exemplified by chaos, is the norm for complicated games with many players.
Collapse
|
5
|
Bairey E, Kelsic ED, Kishony R. High-order species interactions shape ecosystem diversity. Nat Commun 2016; 7:12285. [PMID: 27481625 PMCID: PMC4974637 DOI: 10.1038/ncomms12285] [Citation(s) in RCA: 174] [Impact Index Per Article: 21.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2016] [Accepted: 06/20/2016] [Indexed: 02/02/2023] Open
Abstract
Classical theory shows that large communities are destabilized by random interactions among species pairs, creating an upper bound on ecosystem diversity. However, species interactions often occur in high-order combinations, whereby the interaction between two species is modulated by one or more other species. Here, by simulating the dynamics of communities with random interactions, we find that the classical relationship between diversity and stability is inverted for high-order interactions. More specifically, while a community becomes more sensitive to pairwise interactions as its number of species increases, its sensitivity to three-way interactions remains unchanged, and its sensitivity to four-way interactions actually decreases. Therefore, while pairwise interactions lead to sensitivity to the addition of species, four-way interactions lead to sensitivity to species removal, and their combination creates both a lower and an upper bound on the number of species. These findings highlight the importance of high-order species interactions in determining the diversity of natural ecosystems.
Collapse
Affiliation(s)
- Eyal Bairey
- Department of Physics, Technion—Israel Institute of Technology, Haifa 3200003, Israel
| | - Eric D. Kelsic
- Department of Systems Biology, Harvard Medical School, Boston, Massachusetts 02115, USA
| | - Roy Kishony
- Department of Systems Biology, Harvard Medical School, Boston, Massachusetts 02115, USA
- Department of Biology and Department of Computer Science, Technion—Israel Institute of Technology, Haifa 3200003, Israel
| |
Collapse
|
6
|
Obuchi T, Kabashima Y, Tokita K. Multiple peaks of species abundance distributions induced by sparse interactions. Phys Rev E 2016; 94:022312. [PMID: 27627322 DOI: 10.1103/physreve.94.022312] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2015] [Indexed: 06/06/2023]
Abstract
We investigate the replicator dynamics with "sparse" symmetric interactions which represent specialist-specialist interactions in ecological communities. By considering a large self-interaction u, we conduct a perturbative expansion which manifests that the nature of the interactions has a direct impact on the species abundance distribution. The central results are all species coexistence in a realistic range of the model parameters and that a certain discrete nature of the interactions induces multiple peaks in the species abundance distribution, providing the possibility of theoretically explaining multiple peaks observed in various field studies. To get more quantitative information, we also construct a non-perturbative theory which becomes exact on tree-like networks if all the species coexist, providing exact critical values of u below which extinct species emerge. Numerical simulations in various different situations are conducted and they clarify the robustness of the presented mechanism of all species coexistence and multiple peaks in the species abundance distributions.
Collapse
Affiliation(s)
- Tomoyuki Obuchi
- Department of Mathematical and Computing Science, Tokyo Institute of Technology, Yokohama 226-8502, Japan
| | - Yoshiyuki Kabashima
- Department of Mathematical and Computing Science, Tokyo Institute of Technology, Yokohama 226-8502, Japan
| | - Kei Tokita
- Graduate School of Information Science, Nagoya University, Nagoya 464-8601, Japan
| |
Collapse
|
7
|
Wagner N, Ashkenasy G. How Catalytic Order Drives the Complexification of Molecular Replication Networks. Isr J Chem 2015. [DOI: 10.1002/ijch.201400198] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
|
8
|
Analytical theory of species abundance distributions of a random community model. POPUL ECOL 2015. [DOI: 10.1007/s10144-014-0476-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
9
|
Helbing D, Johansson A. Cooperation, norms, and revolutions: a unified game-theoretical approach. PLoS One 2010; 5:e12530. [PMID: 20967256 PMCID: PMC2953489 DOI: 10.1371/journal.pone.0012530] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2010] [Accepted: 08/02/2010] [Indexed: 12/04/2022] Open
Abstract
Background Cooperation is of utmost importance to society as a whole, but is often challenged by individual self-interests. While game theory has studied this problem extensively, there is little work on interactions within and across groups with different preferences or beliefs. Yet, people from different social or cultural backgrounds often meet and interact. This can yield conflict, since behavior that is considered cooperative by one population might be perceived as non-cooperative from the viewpoint of another. Methodology and Principal Findings To understand the dynamics and outcome of the competitive interactions within and between groups, we study game-dynamical replicator equations for multiple populations with incompatible interests and different power (be this due to different population sizes, material resources, social capital, or other factors). These equations allow us to address various important questions: For example, can cooperation in the prisoner's dilemma be promoted, when two interacting groups have different preferences? Under what conditions can costly punishment, or other mechanisms, foster the evolution of norms? When does cooperation fail, leading to antagonistic behavior, conflict, or even revolutions? And what incentives are needed to reach peaceful agreements between groups with conflicting interests? Conclusions and Significance Our detailed quantitative analysis reveals a large variety of interesting results, which are relevant for society, law and economics, and have implications for the evolution of language and culture as well.
Collapse
|
10
|
Wagner N, Tannenbaum E, Ashkenasy G. Second order catalytic quasispecies yields discontinuous mean fitness at error threshold. PHYSICAL REVIEW LETTERS 2010; 104:188101. [PMID: 20482213 DOI: 10.1103/physrevlett.104.188101] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/24/2009] [Indexed: 05/29/2023]
Abstract
The quasispecies model describes processes related to the origin of life and viral evolutionary dynamics. We discuss how the error catastrophe that reflects the transition from localized to delocalized quasispecies population is affected by catalytic replication of different reaction orders. Specifically, we find that second order mechanisms lead to a discontinuity in the mean fitness of the population at the error threshold. This is in contrast to the behavior of the first order, autocatalytic replication mechanism considered in the standard quasispecies model. This suggests that quasispecies models with higher order replication mechanisms produce discontinuities in the mean fitness, and hence the viable population fraction as well, at the error threshold, while lower order replication mechanisms yield a continuous mean fitness function. We discuss potential implications for understanding replication in the RNA world and in virology.
Collapse
Affiliation(s)
- Nathaniel Wagner
- Department of Chemistry, Ben-Gurion University of the Negev, Be'er-Sheva, Israel
| | | | | |
Collapse
|
11
|
Helbing D, Johansson A. Evolutionary dynamics of populations with conflicting interactions: classification and analytical treatment considering asymmetry and power. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 81:016112. [PMID: 20365437 DOI: 10.1103/physreve.81.016112] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/29/2009] [Indexed: 05/29/2023]
Abstract
Evolutionary game theory has been successfully used to investigate the dynamics of systems, in which many entities have competitive interactions. From a physics point of view, it is interesting to study conditions under which a coordination or cooperation of interacting entities will occur, be it spins, particles, bacteria, animals, or humans. Here, we analyze the case, where the entities are heterogeneous, particularly the case of two populations with conflicting interactions and two possible states. For such systems, explicit mathematical formulas will be determined for the stationary solutions and the associated eigenvalues, which determine their stability. In this way, four different types of system dynamics can be classified and the various kinds of phase transitions between them will be discussed. While these results are interesting from a physics point of view, they are also relevant for social, economic, and biological systems, as they allow one to understand conditions for (1) the breakdown of cooperation, (2) the coexistence of different behaviors ("subcultures"), (3) the evolution of commonly shared behaviors ("norms"), and (4) the occurrence of polarization or conflict. We point out that norms have a similar function in social systems that forces have in physics.
Collapse
|
12
|
|
13
|
Yoshino Y, Galla T, Tokita K. Rank abundance relations in evolutionary dynamics of random replicators. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 78:031924. [PMID: 18851082 DOI: 10.1103/physreve.78.031924] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/21/2008] [Revised: 07/09/2008] [Indexed: 05/26/2023]
Abstract
We present a nonequilibrium statistical mechanics description of rank abundance relations (RAR) in random community models of ecology. Specifically, we study a multispecies replicator system with quenched random interaction matrices. We here consider symmetric interactions as well as asymmetric and antisymmetric cases. RARs are obtained analytically via a generating functional analysis, describing fixed-point states of the system in terms of a small set of order parameters, and in dependence on the symmetry or otherwise of interactions and on the productivity of the community. Our work is an extension of Tokita [Phys. Rev. Lett. 93, 178102 (2004)], where the case of symmetric interactions was considered within an equilibrium setup. The species abundance distribution in our model come out as truncated normal distributions or transformations thereof and, in some case, are similar to left-skewed distributions observed in ecology. We also discuss the interaction structure of the resulting food-web of stable species at stationarity, cases of heterogeneous cooperation pressures as well as effects of finite system size and of higher-order interactions.
Collapse
Affiliation(s)
- Yoshimi Yoshino
- Graduate School of Science and Cybermedia Center, Osaka University, Toyonaka, Osaka 560-0043, Japan.
| | | | | |
Collapse
|
14
|
|
15
|
|
16
|
Michaelian K. Thermodynamic stability of ecosystems. J Theor Biol 2005; 237:323-35. [PMID: 15978624 DOI: 10.1016/j.jtbi.2005.04.019] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2004] [Revised: 03/29/2005] [Accepted: 04/25/2005] [Indexed: 11/23/2022]
Abstract
The stability of ecosystems during periods of stasis in their macro-evolutionary trajectory is studied from a non-equilibrium thermodynamic perspective. Individuals of the species are considered as units of entropy production and entropy exchange in an open thermodynamic system. Within the framework of the classical theory of irreversible thermodynamics, and under the condition of constant external constraints, such a system will naturally evolve toward a globally stable thermodynamic stationary state. It is thus suggested that the ecological steady state, or stasis, is a particular case of the thermodynamic stationary state, and that the evolution of community stability through natural selection is a manifestation of non-equilibrium thermodynamic directives. Furthermore, it is argued that punctuation of stasis leading to ecosystem succession, may be a manifestation of non-equilibrium "phase transitions" brought on by a change of external constraints through a thermodynamic critical point.
Collapse
Affiliation(s)
- K Michaelian
- Instituto de Física, Universidad Nacional Autónoma de México, A.P. 20-364, 01000 México D.F., Mexico.
| |
Collapse
|
17
|
Santos DOC, Fontanari JF. Model ecosystems with random nonlinear interspecies interactions. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2004; 70:061914. [PMID: 15697409 DOI: 10.1103/physreve.70.061914] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/29/2004] [Indexed: 05/24/2023]
Abstract
The principle of competitive exclusion in ecology establishes that two species living together cannot occupy the same ecological niche. Here we present a model ecosystem in which the species are described by a series of phenotypic characters and the strength of the competition between two species is given by a nondecreasing (modulating) function of the number of common characters. Using analytical tools of statistical mechanics we find that the ecosystem diversity, defined as the fraction of species that coexist at equilibrium, decreases as the complexity (i.e., number of characters) of the species increases, regardless of the modulating function. By considering both selective and random elimination of the links in the community web, we show that ecosystems composed of simple species are more robust than those composed of complex species. In addition, we show that the puzzling result that there exists either rich or poor ecosystems for a linear modulating function is not typical of communities in which the interspecies interactions are determined by a complementarity rule.
Collapse
Affiliation(s)
- Danielle O C Santos
- Instituto de Física de São Carlos, Universidade de São Paulo, Caixa Postal 369, 13560-970 São Carlos, São Paulo, Brazil
| | | |
Collapse
|
18
|
Tokita K. Species abundance patterns in complex evolutionary dynamics. PHYSICAL REVIEW LETTERS 2004; 93:178102. [PMID: 15525129 DOI: 10.1103/physrevlett.93.178102] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/03/2004] [Indexed: 05/24/2023]
Abstract
An analytic theory of species abundance patterns (SAPs) in biological networks is presented. The theory is based on multispecies replicator dynamics equivalent to the Lotka-Volterra equation, with diverse interspecies interactions. Various SAPs observed in nature are derived from a single parameter. The abundance distribution is formed like a widely observed left-skewed lognormal distribution. As the model has a general form, the result can be applied to similar patterns in other complex biological networks, e.g., gene expression.
Collapse
Affiliation(s)
- Kei Tokita
- Program for Evolutionary Dynamics, Harvard University, One Brattle Square, Cambridge, Massachusetts 02138, USA.
| |
Collapse
|
19
|
de Oliveira VM, Fontanari JF. Complementarity and diversity in a soluble model ecosystem. PHYSICAL REVIEW LETTERS 2002; 89:148101. [PMID: 12366077 DOI: 10.1103/physrevlett.89.148101] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/05/2002] [Indexed: 05/23/2023]
Abstract
Complementarity among species with different traits is one of the basic processes affecting biodiversity, defined as the number of species in the ecosystem. We present here a soluble model of microbial-based ecosystems in which the species are characterized by binary traits and their pairwise interactions follow a complementarity principle. Manipulation of the species composition, and so the study of its effects on the species diversity, is achieved through the introduction of a bias parameter favoring one of the traits. Using statistical mechanics tools we find explicit expressions for the allowed values of the equilibrium species concentrations in terms of the control parameters of the model.
Collapse
Affiliation(s)
- Viviane M de Oliveira
- Instituto de Física de São Carlos, Universidade de São Paulo, Caixa Postal 369, 13560-970 São Carlos SP, Brazil
| | | |
Collapse
|
20
|
Sayama H, de Aguiar MAM, Bar-Yam Y, Baranger M. Spontaneous pattern formation and genetic invasion in locally mating and competing populations. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2002; 65:051919. [PMID: 12059605 DOI: 10.1103/physreve.65.051919] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2001] [Indexed: 05/23/2023]
Abstract
We present a theoretical model of evolution of spatially distributed populations in which organisms mate with and compete against each other only locally. We show using both analysis and numerical simulation that the typical dynamics of population density variation is a spontaneous formation of isolated groups due to competition for resources. The resulting spatial separation between groups strongly affects the process of genetic invasion by local reproductive mixing, and spatially inhomogeneous genetic distributions are possible in the final states. We then consider a specific version of this model in the presence of disruptive selection, favoring two fittest types against their genetic intermediates. This case can be simplified to a system that involves just two nonconserved order parameters: population density and type difference. Since the coexistence of two fittest types is unstable in this case, symmetry breaking and coarsening occur in type difference, implying eventual dominance by one type over another for finite populations. However, such coarsening patterns may be pinned by the spontaneously generated spatial separation between isolated groups. The long-term evolution of genetic composition is found to be sensitive to the ratio of the mating and competition ranges, and other parameters. Our model may provide a theoretical basis for consideration of various properties of spatially extended evolutionary processes, including spontaneous formation of subpopulations and lateral invasion of different types.
Collapse
Affiliation(s)
- Hiroki Sayama
- New England Complex Systems Institute, Cambridge, Massachusetts 02138, USA
| | | | | | | |
Collapse
|
21
|
de Oliveira VM, Fontanari JF. Extinctions in the random replicator model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2001; 64:051911. [PMID: 11735972 DOI: 10.1103/physreve.64.051911] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/23/2001] [Indexed: 05/23/2023]
Abstract
The statistical properties of an ecosystem composed of species interacting via pairwise, random interactions and deterministic, concentration limiting self-interactions are studied analytically with tools of equilibrium statistical mechanics of disordered systems. Emphasis is given to the effects of externally induced extinction of a fixed fraction of species at the outset of the evolutionary process. The manner the ecosystem copes with the initial extinction event depends on the degree of competition among the species as well as on the strength of that event. For instance, in the regime of high competition the ecosystem diversity, given by the fraction of surviving species, is practically insensitive to the strength of the initial extinction provided it is not too large, while in the less competitive regime the diversity decreases linearly with the size of the event. In the case of large extinction events we find that no further biotic extinctions take place and, furthermore, that rare species become very unlikely to be found in the ecosystem at equilibrium. In addition, we show that the reciprocal of the Edwards-Anderson order parameter yields a good measure of the diversity of the model ecosystem.
Collapse
Affiliation(s)
- V M de Oliveira
- Instituto de Física de São Carlos, Universidade de São Paulo, Caixa Postal 369, 13560-970 São Carlos São Paolo, Brazil
| | | |
Collapse
|