# The Silent Cooperator: An Epigenetic Model for Emergence of Altruistic Traits in Biological Systems

Authors: I. Hashem; D. Telen; P. Nimmegeers; J. Van Impe
Journal: Complexity
 (2018)

---

## Abstract

Spatial evolutionary game theory explains how cooperative traits can survive the intense competition in biological systems. If the spatial distribution allows cooperators to interact with each other frequently, the benefits of cooperation will outweigh the losses due to exploitation by selfish organisms. However, for a cooperative behavior to get established in a system, it needs to be found initially in a sufficiently large cluster to allow a high frequency of intracooperator interactions. Since mutations are rare events, this poses the question of how cooperation can arise in a biological system in the first place. We present a simple model which captures two concepts from genetics that can explain how evolution overcomes the emergence problem. The first concept is, often in nature, a gene may not express its phenotype except under specific environmental conditions, rendering it to be a “silent” gene. The second key idea is that a neutral gene, one that does not harm or improve an organism’s survival chances, can still spread through a population if it is physically near to another gene that is positively selected. Through these two ideas, our model offers a possible solution to the fundamental problem of emergence of cooperation in biological systems.


---

## Body

## 1. Introduction

Cooperation is an unexpected product of evolution. In a struggle and race for existence between living organisms, it is intriguing that cooperative phenomena are widely spread in nature: bacteria share enzymes among each other to boast their overall growth; altruistic insects defend their colonies and, in the process, waste their shot at reproduction; and humans engage in a wide spectrum of altruistic activities, from paying taxes to donating to complete strangers. Cooperation is indeed widespread across all levels of biological societies [1, 2]. However, despite this fact, classic evolutionary game theory could not explain how cooperators in nature can survive the risk of being exploited by “cheaters,” organisms that aim to maximize their self-interest. This problem can be illustrated mathematically using the prisoner’s dilemma as a game of evolution. The prisoner’s dilemma is a standard example in game theory that can model a wide range of cooperative phenomena in a variety of biological systems, including viruses [3], cancer cells [4], bacterial life [5], and mammals [6].The payoffs for a prisoner’s dilemma (PD) game are summarized in Table1. A cooperator provides a gain g and endures a cooperation cost c. Defectors receive nothing when encountering each other. In contrast, if both agents are cooperative, each of them gets g−c. A cooperator interacting with a defector gets a payoff of −c as it gets exploited. On the other hand, the defector reaps the benefit of exploitation g without bearing any cost. Hence, it is seen that while cooperation is better for the group, defection is always better for the individual. A possible solution to this problem has been offered in the framework of spatial evolutionary game theory [2]. Including the effects of spatial arrangements can make the higher frequency of interactions among cooperators compensate the losses due to encounters with cheaters. The results of including spatial dimensions in the study of evolution of cooperation have been demonstrated by numerous studies [2, 7–9]. In [2], cells playing a PD game with each other are distributed on a square lattice. In a biological model, cells do not choose their strategies. Instead, they are genetically programmed to be either defectors or cooperators as they model living organisms. At each round, every cell plays a PD game with each of its neighbors. After that, a cell’s fitness, the sum of the games’ payoffs, is compared with its neighbors. If a neighbor has higher fitness than the focal cell, the focal cell dies and its location gets repopulated by a cell with the same strategy of the most successful neighbor. Similar studies have been performed using different update rules [7, 9] and on other social games, a public good game, for instance [10], showing that cooperative traits can survive on a lattice under certain parameter ranges.Table 1
Payoffs for agenta against agent b in a prisoner’s dilemma game; g>0 and c>0.
b
C
D
a
C
g−c
−c
D
g
0A feature such models share is that for cooperators to do a successful invasion in a world of defectors, an initial cluster of cooperators is needed [2, 7, 8, 11]. Spatial effects can only come to play if a cooperator interacts frequently with other cooperators. Hence, there is always a need of a critical population size of altruists so that they can survive and, in some situations, invade a lattice of selfish cells. The problem is even more severe in a public good game in which cooperators do not reap the benefits of cooperation except when their population reaches a critical threshold [12]. This raises the problem of how a simple cooperative trait can first emerge in a competitive biological system. We define a simple trait as one that does not get adjusted depending on the specific interaction it is involved in. Thus, it does not rely on complex behaviors such as punishment, memory, or communication. One possible answer to how simple cooperative traits appear is that with high enough population, and very long time, a highly unlikely coincidence could happen resulting in the formation of a cluster of cooperators exceeding the critical size. This can occur due to a number of mutations happening very close to each other or by introducing uncertainty in the game rules such that there is a tiny probability that a cooperator with low payoffs will invade a defector with higher payoffs. A common view currently in explaining the emergence of a simple cooperative strategy is that in a biological system, given enough time, even a highly unlikely event can happen [13].In this paper, we aim to provide an alternative solution to the problem of emergence of cooperation. By means of known genetic phenomena, we show how a cooperative mutation can overcome the hurdle of the initial cluster size to establish itself in a biological system. The paper is structured as follows: Section2 describes the concept and formulation of the spatial cooperation model examined in this work. The results and a discussion of the key parameters of the model are provided in Section 3. Finally, Section 4 summarizes the paper’s conclusions.


## 2. Materials and Methods

The model presented in this section is inspired by how social bacteria overcome the problem of establishing cooperation when growing in a biofilm mode, a highly competitive biological system. Hence, before introducing our model, we first take a closer look into the social life of one of the most successful communities in earth: biofilms.


### 2.1. Cooperation in Biofilms

A biofilm is a community of bacterial cells living together on a surface enclosed by a self-produced polymeric matrix. The biofilm lifecycle starts with the colonization of a surface by motile microorganisms which get attached to the surface. Subsequently, they grow by cell division till a mature biofilm is formed. Finally, the dispersal of cells from its top leads to the colonization of new surfaces [14]. During the growth phase of a biofilm, bacteria often engage in cooperative behaviors, such as the production of extracellular enzymes to break down the nutrients for the benefit of the whole group and growing at a slower mode with a higher yield to increase the nutrient availability to the population as a whole. The distinctive feature of biofilms which allows cooperative traits to persist is the extracellular matrix, since it causes diffusion limitations in the biofilm and hence intensifies the spatial effects of segregation. In a study by [15], it is shown that cooperative cells can survive in a lowly mixed biofilm owing to the diffusion constraints of the polymeric matrix which reduce the effects of exploitation of cooperators by surrounding selfish cells. This result parallels the conclusion of [2, 11] at the more abstract setting of a spatial PD game.On the other hand, in a highly mixed biofilm, such constitutive cooperative cells cannot survive as they will be outgrown by selfish cells. Hence, the need for more complex cooperative strategies arises. One technique by which cooperation can emerge in a highly mixed biofilm isquorum sensing. Quorum sensing is the process of producing chemical signals by bacterial cells to accumulate in the biofilm. When the concentration of the signals sensed by the cell reaches a certain threshold, quorum, it induces a change in the behavior/phenotype of the cell. In a study by [16], when colonizing a new surface, it is shown that cooperative cells using quorum sensing can dominate highly mixed biofilms by delaying the expression of their cooperative gene. In the beginning, the cooperative gene is unexpressed. However, as the biofilm grows, the density of cooperative quorum-sensing cells increases as well as the concentration of their chemical signals. The quorum signal concentration eventually reaches the threshold required for activating the cooperative gene. By that time, quorum-sensing cells are spatially segregated from selfish cells. Hence, they are less likely to be exploited, and they dominate the biofilm. Quorum sensing effectively acts as a timer to delay the expression of an altruistic trait till the cells having the “silent” cooperative genes become spatially separated from selfish cells [16]. Thus, it reduces the risks of cooperative cells being exploited and increases their benefits from cooperation. This represents a highly sophisticated behavior where bacterial cells deliberately delay the expression of the cooperative gene till segregation occurs. By the means of quorum sensing, the silent gene gets activated by a change in environment that the bacteria themselves induce. This way, cooperative cells overcome the early highly mixed phase of the colonization process.It should be noted that cooperative biofilm species have developed such sophisticated behavior as a result of immense selection pressure and the crucial need to establish cooperation in newborn communities repeatedly as they colonize new surfaces [16]. Using quorum sensing, they efficiently fine-tune the timing of the expression of their cooperative genes. Nevertheless, the concept of emergence of cooperation by a silent gene that gets activated at a future point of time due to possible changes in environmental conditions can still be a route for emergence of cooperation in general biological systems, as we will explain in the next section.


### 2.2. Model Concept

Once it is established in a population, a cooperative gene could be evolutionary stable if the population structure allows high frequency of intracooperator interactions [2, 11]. However, a mutation resulting in a cooperative trait, such as the production of public good or sharing vital resources, faces steep odds in order to spread in a biological system in the first place. To explain a possible route for the emergence of cooperation, our model relies on two genetic phenomena. Firstly, it is common in nature that gene expression could rely on environmental factors [17, 18]. A gene could give rise to a certain phenotype only beyond certain temperature [19], concentration of a chemical such as oxygen [20], or the pH level of the environment [21]. This alteration of gene expression due to nongenetic influences is referred to as epigenetic activation [22]. In the absence of the appropriate level of the environmental factor to trigger it, a gene could stay unexpressed (silent) with neither positive nor negative impact on the organism’s fitness [18]. Hence, as illustrated in Figure 1, for a cooperative trait to first appear in a population via a mutation, it can be either through a gene that has been instantly expressed, as a consequence of being at the required environmental conditions/needing no environmental trigger to be active, or as a silent gene that will be activated at under the appropriate environmental conditions at a point in the future. In the latter case, the silent cooperative trait has a time interval at which it exerts no influence on its carriers’ fitness. During this interval, if it happened to spread to sufficient extent in the population, the cooperative trait could be evolutionary stable in the population by the time it gets expressed.Figure 1
Silent genes: gene expression is often dependent on environmental factors. At their absence, the gene will stay unexpressed, although it can still be inherited and, in certain cases, spread in a population. The initial and terminal states here could represent low/high temperature, pH, or any combination of environmental triggers required for the activation of the gene in hand.
(a)
Appearance of a cooperative gene (C) in the population that is always expressed, at any environmental condition
(b)
Appearance of a gene (S) that requires a future state of the environment to be expressed. Till then, the gene carriers act as defective (D)This brings us to the question of how a gene that does not alter the fitness of its carrier, such as silent genes, could possibly spread in a population. A possible mechanism is genetic hitchhiking. First proposed by [23], it is the increase in the frequency of a gene in the population not because it is itself under positive selection but due to being physically close, on the same DNA chain, of another gene that is beneficial to the organism and, hence, under positive selection [24], as illustrated in Figure 2.Figure 2
Genetic hitchhiking: the frequency of a gene could increase in the population due to lying at the same chromosome of another advantageous gene. In these “domino organisms,” the top gene, the number of dots, represents a trait that is advantageous to its carrier, such as resistance to toxins or diseases. Hence, as the domino organisms with the highest dot number get positively selected, their bottom genes, which have no influence on their fitness, also spread in the population.As a consequence of combining the two ideas, a possible route for the emergence of an altruistic trait could be suggested, as illustrated in Figure3. If a mutation produced a cooperative trait in the absence of the required abundance of the environmental factors necessary for its expression, the silent gene will not get exploited. Additionally, if it happened to hitchhike a selective sweep of another beneficial gene till it exceeds a critical spatial distribution, it could get established in the population once it gets triggered as the environment evolves to a more favorable state for its expression.Figure 3
Silent genes and genetic hitchhiking: a cooperative gene in a world of defectors will be heavily exploited with no other cooperators to interact with. Hence, it has low chances of spreading in the population. A silent gene, on the other hand, does not affect its carriers’ fitness till activation. And, it might hitchhike a selective sweep if it happened to be in an organism with relatively higher fitness, due to influence of the rest of the genes.


### 2.3. Model Description

In this paper, we adopt the spatial model from [7]; however, the main conclusions can be generalized to other models in which a critical size of cooperator clusters is needed for a successful invasion. The model from [7] consists of a L×L two-dimensional square lattice with periodic boundary conditions, where each position is occupied by a cell that can be either cooperative or defective. Each cell engages in a pairwise PD game with its neighbors. We consider the trait of taking part in a PD interaction as gene Y, which can be either cooperative C or defective D, Y∈C,D. The payoffs of a single interaction are normalized to depend on a single parameter u, as shown in Table 2 [25].Table 2
Normalized payoffs for agenta against agent b in a prisoner’s dilemma game; 0<u<1.
b
C
D
a
C
1
0
D
1+u
uThe population is updated as follows: a focal cellj is selected at random order. The cell interacts with all neighbor cells bi in its extended Moore’s neighborhood of radius r=3, that is, all cells reached within three chess king’s moves, with i∈1,2,…,Nm. Nm is the number of cells in Moore’s neighborhood of the focal cell, with Nm=2r+12−1. The fitness of cell j in any generation t, Πjt, is equal to the fitness resulting from gene Y, πj,Yt, which is defined as the sum of the payoffs from the pairwise PD interactions with neighbor cells weighted according to the distance between the neighbor cell bi from the focal cell a. This can be expressed by the following equation, which is adapted from [7]:
(1)Πjt=πj,Yt=∑i=1Nm1dbiEbit,with Ebit as the payoff resulting from an encounter with neighbor cell bi, which takes one of the values in Table 2 depending on the strategies of both j and bi. dbi is the distance between bi and the focal cell j, dbi∈1,2,…,r. The payoffs are weighted depending on the distance as illustrated in Figure 4(a), such that the further the cell, the less the weight of its effect on a focal cell. After the fitness of all cells in the population is calculated, each cell’s fitness is compared with the fitness of a random direct neighbor. As shown in Figure 4(b), for a focal cell j, a neighbor cell bi within its simple Moore’s neighborhood is chosen at random. The fitness of the two cells Πjt and Πbit is compared. If cell j has higher fitness, it keeps its position, repopulating it with a cell that inherits its own gene Y and has the same strategy of the parent. On the other hand, if bi has higher fitness, it could invade the position of j, repopulating it with a cell having the strategy of bi. The probability of such invasion depends on the difference between the fitness of the two cells, expressed as follows:
(2)fΠbit−Πjt=Πbit−Πjtk,ifΠbit>Πjt,0,otherwise,with fΠbit−Πjt as the invasion probability and k as a normalization constant such that fΠbit−Πjt∈0,1, here k=∑i=1Nm1/dbi1+u=241+u. Using these population update rules, four extreme scenarios are shown in Figure 5. In Figure 5(b), a cooperator in a world of defectors gets a zero payoff and, in consequence, has a zero chance of invading any other cell. As explained in [7], a cluster of cooperators needs to be larger than a critical size to invade such a system. Hence, a major question is posed by such models: How could cooperation possibly emerge by rare mutations in a system dominated initially by defectors?Figure 4
An illustration of the update rules. A focal cell plays the prisoner’s dilemma game with cells in its extended Moore neighborhood; scores are weighted depending on the distance between the neighbor and the focal cell. After all cells in the grid are assigned fitness values, each cell competes with one random neighbor.
(a)
The weights of the interactions between a focal cell and neighbor cells in an extended Moore’s neighborhood of radius 3
(b)
A cell can have its position invaded by a randomly chosen cell from its simple Moore’s neighborhood, one chess king’s move from the focal cell, depending on the difference in fitness between the two cellsFigure 5
An illustration of a focal cell interacting with neighbors in its extended Moore neighborhood with range = 3, under 4 extreme scenarios (withu<1).
(a)
Cooperator among cooperators, focal cell payoffs = 24
(b)
Cooperator among defectors, focal cell payoffs = 0
(c)
Defector among cooperators, focal cell payoffs = 24(1 +u)
(d)
Defector among defectors, focal cell payoffs = 24uTo address this question, we introduce two changes in the model to abstract two concepts from genetics that could offer a possible route to the emergence of cooperation in spatial systems. The first phenomenon is that genes in nature often require specific environmental factors to be active. In the absence of these environmental triggers, the gene stays unexpressed (silent), although it can be still inherited. We model this by extending the definition of gene Y such that it can be either cooperativeC, defective D, or silent S, Y∈C,D,S. S aims to model the behavior of a cooperative mutation that only gets triggered after the environment reaches a favorable state for its expression (see Figure 1). This could be the result of the accumulation of a certain chemical or a certain degree of a physical parameter such as the temperature or the pH of the environment. Thus, S can be formulated as follows:
(3)S≡D,t<t′,C,t>t′,with t′ as a random variable that models the time till activation of the gene, and the equivalence symbol ≡ here means “expresses the trait of”. Equation (3) means that a silent gene acts as a defector, expressing the defector phenotype, till the moment when a random environmental event activates the cooperative phenotype. Hence, while the environment is not explicitly modeled here, each independent mutation could be considered to be activated by a specific state of the environment with a waiting time t′. For example, if the environmental trigger that activates the silent gene is high temperatures, each independent mutation giving rise to this gene is assumed to be activated by a specific high temperature which will be realized by the environment after time t′.Hence, when a silent gene appears in a population of defectors, its carrier keeps acting as a defector. And in a world of defectors, till the gene gets expressed, it has neither negative nor positive effect on its carrier’s fitness, although it can still be inherited and, in some cases, spread in the population. The second genetic phenomena to be modeled, known as genetic hitchhiking, are related to how the frequency of an unexpressed gene could increase in the population if it happened to be in the same DNA chain of another gene that is positively selected. To model the hitchhiking effect, the other modification to be introduced to the model of [7] is that the fitness of any cell j in the grid, Πjt, will not only be a result of the PD (gene Y) interactions with neighbors, but will also include another component. Hence, the fitness equation used in our model will extend (1) by including an additional fitness term as follows:
(4)Πjt=πj,Xt+πj,Yt,with πj,Yt as the gain in fitness resulting from the action of gene Y, πj,Yt∈0,241+u, depending on the strategy of the cell and its neighbors, calculated identically to (1). Besides, we model the gain in fitness πj,Xt resulting from the action of all other genes of the cell, aggregated as gene X, which contribute to the cell fitness independently of the PD game. The initial generation of cells is assigned a base fitness due to gene X that is equal to zero for all cells in the population, πj,Xt=0=0, for any cell j in the grid. However, each cell has a probability PXμ of a mutation in gene X leading to a slight increase in its fitness due to gene X, μ, such that for any cell j in generation t, the fitness due to gene X, πj,Xt, is as follows:
(5)πj,Xt=πj,Xt−1,ifnobeneficialmutationoccurred,πj,Xt−1+μ,ifabeneficialmutationoccurred.In this way, gene X aims to model the increase in fitness of the cell resulting from any other gene rather than the one taking part in the PD game. If all cells have the same gene Y strategy, all of them are defectors or all are cooperators, a cellj with higher πj,Xt can still invade the population (see Figures 2 and 3). A flowchart of the model is presented in Figure 6 and an illustration of the fitness function is presented in Figure 7. A summary of notations used is presented in Table 3. All model simulations are executed using the Repast Simphony tool kit, an agent-based simulation platform created by [26].Figure 6
A flowchart of the spatial prisoner’s dilemma model.Figure 7
Each cell in the grid has two sources of fitness.πj,Yt is the familiar fitness function resulting from the PD interactions with neighbor cells. Additionally, there is πj,Xt which represents the contribution in fitness from all the other genes of the cell, aggregated as gene X. All cells in the grid are assigned the same background fitness at the beginning of the simulation, with πj,Xt=0=0 for any cell j in the grid.Table 3
Notation system.
Notation
Meaning
Πjt
Fitness of cellj in generation t
πj,Yt
Fitness component of cellj resulting from PD interactions
πj,Xt
Fitness component of cellj resulting from background genes, aggregated as gene X
Gene Y
The trait of taking part in PD interactions; it can be either C, D, or S
C
A variant of gene Y which always expresses cooperative behavior
D
A variant of gene Y which always expresses defector’s behavior
S
A variant of gene Y which expresses cooperative behavior only when triggered by a future state of the environment
t′
Time, in generations, between appearance ofS in the population and its activation
Gene X
An abstraction of all genes not taking part in the PD game that contribute to the cell’s fitness
μ
The increase in the background fitness,πj,Xt, of cell j due to a beneficial mutation in the background genes, unrelated to the PD interactions


## 2.1. Cooperation in Biofilms

A biofilm is a community of bacterial cells living together on a surface enclosed by a self-produced polymeric matrix. The biofilm lifecycle starts with the colonization of a surface by motile microorganisms which get attached to the surface. Subsequently, they grow by cell division till a mature biofilm is formed. Finally, the dispersal of cells from its top leads to the colonization of new surfaces [14]. During the growth phase of a biofilm, bacteria often engage in cooperative behaviors, such as the production of extracellular enzymes to break down the nutrients for the benefit of the whole group and growing at a slower mode with a higher yield to increase the nutrient availability to the population as a whole. The distinctive feature of biofilms which allows cooperative traits to persist is the extracellular matrix, since it causes diffusion limitations in the biofilm and hence intensifies the spatial effects of segregation. In a study by [15], it is shown that cooperative cells can survive in a lowly mixed biofilm owing to the diffusion constraints of the polymeric matrix which reduce the effects of exploitation of cooperators by surrounding selfish cells. This result parallels the conclusion of [2, 11] at the more abstract setting of a spatial PD game.On the other hand, in a highly mixed biofilm, such constitutive cooperative cells cannot survive as they will be outgrown by selfish cells. Hence, the need for more complex cooperative strategies arises. One technique by which cooperation can emerge in a highly mixed biofilm isquorum sensing. Quorum sensing is the process of producing chemical signals by bacterial cells to accumulate in the biofilm. When the concentration of the signals sensed by the cell reaches a certain threshold, quorum, it induces a change in the behavior/phenotype of the cell. In a study by [16], when colonizing a new surface, it is shown that cooperative cells using quorum sensing can dominate highly mixed biofilms by delaying the expression of their cooperative gene. In the beginning, the cooperative gene is unexpressed. However, as the biofilm grows, the density of cooperative quorum-sensing cells increases as well as the concentration of their chemical signals. The quorum signal concentration eventually reaches the threshold required for activating the cooperative gene. By that time, quorum-sensing cells are spatially segregated from selfish cells. Hence, they are less likely to be exploited, and they dominate the biofilm. Quorum sensing effectively acts as a timer to delay the expression of an altruistic trait till the cells having the “silent” cooperative genes become spatially separated from selfish cells [16]. Thus, it reduces the risks of cooperative cells being exploited and increases their benefits from cooperation. This represents a highly sophisticated behavior where bacterial cells deliberately delay the expression of the cooperative gene till segregation occurs. By the means of quorum sensing, the silent gene gets activated by a change in environment that the bacteria themselves induce. This way, cooperative cells overcome the early highly mixed phase of the colonization process.It should be noted that cooperative biofilm species have developed such sophisticated behavior as a result of immense selection pressure and the crucial need to establish cooperation in newborn communities repeatedly as they colonize new surfaces [16]. Using quorum sensing, they efficiently fine-tune the timing of the expression of their cooperative genes. Nevertheless, the concept of emergence of cooperation by a silent gene that gets activated at a future point of time due to possible changes in environmental conditions can still be a route for emergence of cooperation in general biological systems, as we will explain in the next section.


## 2.2. Model Concept

Once it is established in a population, a cooperative gene could be evolutionary stable if the population structure allows high frequency of intracooperator interactions [2, 11]. However, a mutation resulting in a cooperative trait, such as the production of public good or sharing vital resources, faces steep odds in order to spread in a biological system in the first place. To explain a possible route for the emergence of cooperation, our model relies on two genetic phenomena. Firstly, it is common in nature that gene expression could rely on environmental factors [17, 18]. A gene could give rise to a certain phenotype only beyond certain temperature [19], concentration of a chemical such as oxygen [20], or the pH level of the environment [21]. This alteration of gene expression due to nongenetic influences is referred to as epigenetic activation [22]. In the absence of the appropriate level of the environmental factor to trigger it, a gene could stay unexpressed (silent) with neither positive nor negative impact on the organism’s fitness [18]. Hence, as illustrated in Figure 1, for a cooperative trait to first appear in a population via a mutation, it can be either through a gene that has been instantly expressed, as a consequence of being at the required environmental conditions/needing no environmental trigger to be active, or as a silent gene that will be activated at under the appropriate environmental conditions at a point in the future. In the latter case, the silent cooperative trait has a time interval at which it exerts no influence on its carriers’ fitness. During this interval, if it happened to spread to sufficient extent in the population, the cooperative trait could be evolutionary stable in the population by the time it gets expressed.Figure 1
Silent genes: gene expression is often dependent on environmental factors. At their absence, the gene will stay unexpressed, although it can still be inherited and, in certain cases, spread in a population. The initial and terminal states here could represent low/high temperature, pH, or any combination of environmental triggers required for the activation of the gene in hand.
(a)
Appearance of a cooperative gene (C) in the population that is always expressed, at any environmental condition
(b)
Appearance of a gene (S) that requires a future state of the environment to be expressed. Till then, the gene carriers act as defective (D)This brings us to the question of how a gene that does not alter the fitness of its carrier, such as silent genes, could possibly spread in a population. A possible mechanism is genetic hitchhiking. First proposed by [23], it is the increase in the frequency of a gene in the population not because it is itself under positive selection but due to being physically close, on the same DNA chain, of another gene that is beneficial to the organism and, hence, under positive selection [24], as illustrated in Figure 2.Figure 2
Genetic hitchhiking: the frequency of a gene could increase in the population due to lying at the same chromosome of another advantageous gene. In these “domino organisms,” the top gene, the number of dots, represents a trait that is advantageous to its carrier, such as resistance to toxins or diseases. Hence, as the domino organisms with the highest dot number get positively selected, their bottom genes, which have no influence on their fitness, also spread in the population.As a consequence of combining the two ideas, a possible route for the emergence of an altruistic trait could be suggested, as illustrated in Figure3. If a mutation produced a cooperative trait in the absence of the required abundance of the environmental factors necessary for its expression, the silent gene will not get exploited. Additionally, if it happened to hitchhike a selective sweep of another beneficial gene till it exceeds a critical spatial distribution, it could get established in the population once it gets triggered as the environment evolves to a more favorable state for its expression.Figure 3
Silent genes and genetic hitchhiking: a cooperative gene in a world of defectors will be heavily exploited with no other cooperators to interact with. Hence, it has low chances of spreading in the population. A silent gene, on the other hand, does not affect its carriers’ fitness till activation. And, it might hitchhike a selective sweep if it happened to be in an organism with relatively higher fitness, due to influence of the rest of the genes.


## 2.3. Model Description

In this paper, we adopt the spatial model from [7]; however, the main conclusions can be generalized to other models in which a critical size of cooperator clusters is needed for a successful invasion. The model from [7] consists of a L×L two-dimensional square lattice with periodic boundary conditions, where each position is occupied by a cell that can be either cooperative or defective. Each cell engages in a pairwise PD game with its neighbors. We consider the trait of taking part in a PD interaction as gene Y, which can be either cooperative C or defective D, Y∈C,D. The payoffs of a single interaction are normalized to depend on a single parameter u, as shown in Table 2 [25].Table 2
Normalized payoffs for agenta against agent b in a prisoner’s dilemma game; 0<u<1.
b
C
D
a
C
1
0
D
1+u
uThe population is updated as follows: a focal cellj is selected at random order. The cell interacts with all neighbor cells bi in its extended Moore’s neighborhood of radius r=3, that is, all cells reached within three chess king’s moves, with i∈1,2,…,Nm. Nm is the number of cells in Moore’s neighborhood of the focal cell, with Nm=2r+12−1. The fitness of cell j in any generation t, Πjt, is equal to the fitness resulting from gene Y, πj,Yt, which is defined as the sum of the payoffs from the pairwise PD interactions with neighbor cells weighted according to the distance between the neighbor cell bi from the focal cell a. This can be expressed by the following equation, which is adapted from [7]:
(1)Πjt=πj,Yt=∑i=1Nm1dbiEbit,with Ebit as the payoff resulting from an encounter with neighbor cell bi, which takes one of the values in Table 2 depending on the strategies of both j and bi. dbi is the distance between bi and the focal cell j, dbi∈1,2,…,r. The payoffs are weighted depending on the distance as illustrated in Figure 4(a), such that the further the cell, the less the weight of its effect on a focal cell. After the fitness of all cells in the population is calculated, each cell’s fitness is compared with the fitness of a random direct neighbor. As shown in Figure 4(b), for a focal cell j, a neighbor cell bi within its simple Moore’s neighborhood is chosen at random. The fitness of the two cells Πjt and Πbit is compared. If cell j has higher fitness, it keeps its position, repopulating it with a cell that inherits its own gene Y and has the same strategy of the parent. On the other hand, if bi has higher fitness, it could invade the position of j, repopulating it with a cell having the strategy of bi. The probability of such invasion depends on the difference between the fitness of the two cells, expressed as follows:
(2)fΠbit−Πjt=Πbit−Πjtk,ifΠbit>Πjt,0,otherwise,with fΠbit−Πjt as the invasion probability and k as a normalization constant such that fΠbit−Πjt∈0,1, here k=∑i=1Nm1/dbi1+u=241+u. Using these population update rules, four extreme scenarios are shown in Figure 5. In Figure 5(b), a cooperator in a world of defectors gets a zero payoff and, in consequence, has a zero chance of invading any other cell. As explained in [7], a cluster of cooperators needs to be larger than a critical size to invade such a system. Hence, a major question is posed by such models: How could cooperation possibly emerge by rare mutations in a system dominated initially by defectors?Figure 4
An illustration of the update rules. A focal cell plays the prisoner’s dilemma game with cells in its extended Moore neighborhood; scores are weighted depending on the distance between the neighbor and the focal cell. After all cells in the grid are assigned fitness values, each cell competes with one random neighbor.
(a)
The weights of the interactions between a focal cell and neighbor cells in an extended Moore’s neighborhood of radius 3
(b)
A cell can have its position invaded by a randomly chosen cell from its simple Moore’s neighborhood, one chess king’s move from the focal cell, depending on the difference in fitness between the two cellsFigure 5
An illustration of a focal cell interacting with neighbors in its extended Moore neighborhood with range = 3, under 4 extreme scenarios (withu<1).
(a)
Cooperator among cooperators, focal cell payoffs = 24
(b)
Cooperator among defectors, focal cell payoffs = 0
(c)
Defector among cooperators, focal cell payoffs = 24(1 +u)
(d)
Defector among defectors, focal cell payoffs = 24uTo address this question, we introduce two changes in the model to abstract two concepts from genetics that could offer a possible route to the emergence of cooperation in spatial systems. The first phenomenon is that genes in nature often require specific environmental factors to be active. In the absence of these environmental triggers, the gene stays unexpressed (silent), although it can be still inherited. We model this by extending the definition of gene Y such that it can be either cooperativeC, defective D, or silent S, Y∈C,D,S. S aims to model the behavior of a cooperative mutation that only gets triggered after the environment reaches a favorable state for its expression (see Figure 1). This could be the result of the accumulation of a certain chemical or a certain degree of a physical parameter such as the temperature or the pH of the environment. Thus, S can be formulated as follows:
(3)S≡D,t<t′,C,t>t′,with t′ as a random variable that models the time till activation of the gene, and the equivalence symbol ≡ here means “expresses the trait of”. Equation (3) means that a silent gene acts as a defector, expressing the defector phenotype, till the moment when a random environmental event activates the cooperative phenotype. Hence, while the environment is not explicitly modeled here, each independent mutation could be considered to be activated by a specific state of the environment with a waiting time t′. For example, if the environmental trigger that activates the silent gene is high temperatures, each independent mutation giving rise to this gene is assumed to be activated by a specific high temperature which will be realized by the environment after time t′.Hence, when a silent gene appears in a population of defectors, its carrier keeps acting as a defector. And in a world of defectors, till the gene gets expressed, it has neither negative nor positive effect on its carrier’s fitness, although it can still be inherited and, in some cases, spread in the population. The second genetic phenomena to be modeled, known as genetic hitchhiking, are related to how the frequency of an unexpressed gene could increase in the population if it happened to be in the same DNA chain of another gene that is positively selected. To model the hitchhiking effect, the other modification to be introduced to the model of [7] is that the fitness of any cell j in the grid, Πjt, will not only be a result of the PD (gene Y) interactions with neighbors, but will also include another component. Hence, the fitness equation used in our model will extend (1) by including an additional fitness term as follows:
(4)Πjt=πj,Xt+πj,Yt,with πj,Yt as the gain in fitness resulting from the action of gene Y, πj,Yt∈0,241+u, depending on the strategy of the cell and its neighbors, calculated identically to (1). Besides, we model the gain in fitness πj,Xt resulting from the action of all other genes of the cell, aggregated as gene X, which contribute to the cell fitness independently of the PD game. The initial generation of cells is assigned a base fitness due to gene X that is equal to zero for all cells in the population, πj,Xt=0=0, for any cell j in the grid. However, each cell has a probability PXμ of a mutation in gene X leading to a slight increase in its fitness due to gene X, μ, such that for any cell j in generation t, the fitness due to gene X, πj,Xt, is as follows:
(5)πj,Xt=πj,Xt−1,ifnobeneficialmutationoccurred,πj,Xt−1+μ,ifabeneficialmutationoccurred.In this way, gene X aims to model the increase in fitness of the cell resulting from any other gene rather than the one taking part in the PD game. If all cells have the same gene Y strategy, all of them are defectors or all are cooperators, a cellj with higher πj,Xt can still invade the population (see Figures 2 and 3). A flowchart of the model is presented in Figure 6 and an illustration of the fitness function is presented in Figure 7. A summary of notations used is presented in Table 3. All model simulations are executed using the Repast Simphony tool kit, an agent-based simulation platform created by [26].Figure 6
A flowchart of the spatial prisoner’s dilemma model.Figure 7
Each cell in the grid has two sources of fitness.πj,Yt is the familiar fitness function resulting from the PD interactions with neighbor cells. Additionally, there is πj,Xt which represents the contribution in fitness from all the other genes of the cell, aggregated as gene X. All cells in the grid are assigned the same background fitness at the beginning of the simulation, with πj,Xt=0=0 for any cell j in the grid.Table 3
Notation system.
Notation
Meaning
Πjt
Fitness of cellj in generation t
πj,Yt
Fitness component of cellj resulting from PD interactions
πj,Xt
Fitness component of cellj resulting from background genes, aggregated as gene X
Gene Y
The trait of taking part in PD interactions; it can be either C, D, or S
C
A variant of gene Y which always expresses cooperative behavior
D
A variant of gene Y which always expresses defector’s behavior
S
A variant of gene Y which expresses cooperative behavior only when triggered by a future state of the environment
t′
Time, in generations, between appearance ofS in the population and its activation
Gene X
An abstraction of all genes not taking part in the PD game that contribute to the cell’s fitness
μ
The increase in the background fitness,πj,Xt, of cell j due to a beneficial mutation in the background genes, unrelated to the PD interactions


## 3. Results and Discussion

In this section, we start by examining how the need for a sufficiently large cluster of cooperators hampers the possibilities of inception of cooperation in a spatial system. Next, the effects of emergence of a cooperative trait as a silent gene are explored. The rest of the section explains the biological intuition behind the results and analyses the impact of key parameters.


### 3.1. Emergence of Cooperation in a Biological System

First, to illustrate the problem of emergence of cooperation, we reproduce earlier results in literature by examining the effect of the initial cooperators’ cluster size on their fate, with no mutations in the simulation. In Figures8 and 9, we seed a 50×50 lattice with clusters of cooperators of sizes 9 and 25, respectively, which are cells with gene Y being cooperative (i.e., Y=C). The rest of the positions in the lattice are occupied by defectors, hence gene Y being defective (i.e., Y=D), while the value of the normalized PD variable u is set to u=0.09. In Figures 8 and 9, the interactions among cooperators are not sufficient to compensate for the losses suffered by the exploitation from selfish agents and the cooperative gene goes extinct rapidly.Figure 8
Evolution of cooperator frequency (black) when starting with an initial3×3 cooperator cluster within a 50×50 grid of defectors (yellow).
(a)
Start of simulation
(b)
End of simulation
(c)
Cooperator frequencyFigure 9
Evolution of cooperator frequency (black) when starting with an initial5×5 cooperator cluster within a 50×50 grid of defectors (yellow).
(a)
Start of simulation
(b)
End of simulation
(c)
Cooperator frequencyOn the other hand, when starting with large cluster sizes of6×6 and 8×8 cooperators, as shown in Figures 10 and 11, the initial cooperator population undergoes aggressive expansion at the beginning. The number of cells with Y=C increases quadratically till it reaches a relatively stable value, around 65% of the total population. Simulations are stopped after 105 generations. A dynamic equilibrium state is reached where compact clusters of cooperators are surrounded by stripes of defectors, in accordance with previous research [2, 7, 11, 27]. These results confirm that in a spatial system, cooperators in a PD setting can dominate a world of defectors as their interactions among each other compensate losses due to encounters with cheaters. However, this will happen only if they started from a large enough initial cluster.Figure 10
Evolution of cooperator frequency (black) when starting with an initial6×6 cooperator cluster within a 50×50 grid of defectors (yellow).
(a)
Start of simulation
(b)
End of simulation
(c)
Cooperator frequencyFigure 11
Evolution of cooperator frequency (black) when starting with an initial8×8 cooperator cluster within a 50×50 grid of defectors (yellow).
(a)
Start of simulation
(b)
End of simulation
(c)
Cooperator frequencyWe proceed to study the problem of emergence of cooperation in a world of defectors. The50×50 lattice is initialized by filling it completely with defectors, agents with Y=D, as in Figure 12(a). However, mutations are allowed to occur in the model. The probability that a cell with gene Y=D will mutate to Y=C is PYD→C=1×10−4. Additionally, the probability of a beneficial mutation in gene X, PXμ, is set to be PXμ=1×10−4, with μ=1.0. As observed in Figures 12 and 13, the instantly expressed gene Y=C fails to make a successful invasion of the system. This has been expected as a cooperator in a sea of defectors is very prone to exploitation as not enough interactions occur between cooperators to out-weight exploitation losses. Also, even if a cooperative mutation occurred in a cell with slightly higher X than its neighbors, the losses suffered from exploitation by defectors will be much larger than the gain in fitness due to gene X. Finally, the probability that cooperative mutations will coappear in a narrow spatial range such that they form a cluster of altruists is extremely small.Figure 12
Illustration of a system where mutations of instantly expressed cooperator genes (Y=C) occur in a 50×50 system dominated by defector genes (Y=D).
(a)
Start of simulation
(b)
After 5 × 105 generations
(c)
End of simulationFigure 13
Evolution of the frequency of cooperators in a50×50 system of defectors, where cooperation arises by instantly expressed mutations of gene Y=C.Do mutations of a silent cooperator gene have better chances for invading the system? We test that by allowing, with a probability ofPYD→S=1×10−4, that a cell with gene Y=D will mutate into Y=S, which in turn will get activated after passing of t′ generations after the mutation. For this simulation, t′ is assumed to be exponentially distributed, t′∼Expλ, with a mean value of 200 generations. The simulation is again initialized by filling the lattice with defectors. In the early beginning, no cooperators at all appear, as mutations occurring in the system are still inactivated. After that, the system evolves in a way similar to the previous case of instantly expressed cooperator genes, with cooperators appearing isolated from each other and doomed to extinction. Then, after around 500 generations in this simulation, a sudden outbreak of cooperation appears with a cluster size that is large enough to eventually perform a successful invasion of the system. Thus, cooperation prevailed here. But why exactly did cooperator genes with delayed expression succeed while instantly expressed genes failed? We take a closer look at the two key features of the model which allowed these results to occur in the next section.


### 3.2. Silent Genes and Genetic Hitchhiking

Epigenetics refer to changes in gene expression that are not resulting from changes in the DNA code [28]. A gene could produce a phenotype only when activated by certain environmental conditions. Also, it has been found that a gene could be tuned by an inheritance mechanism to get expressed only after a number of generations, a phenomenon known as transgenerational gene silencing [29, 30].In our model, when a silent gene appears by a mutation, it remains unexpressed for a random period associated with the physical nature of each mutation. Then, after passing oft′ generations, it gets activated at all cells in the system carrying this mutation. This silent interval is however essential for cooperation to break out as the gene can still spread through the population by a mechanism known as genetic hitchhiking. Genetic hitchhiking is the process by which a gene that offers no evolutionary advantage to its carrier can spread through a population due to being associated with another gene that is under a strong selection [31]. Due to this phenomenon, the frequency of a neutral gene can increase within the population if it is on the same DNA chain of another beneficial gene. In the model, mutations of gene X play the role of the beneficial gene. A cell with higher fitness value due to the contribution of gene X can slowly invade the population. As illustrated in Figure 14, by combining the concepts of silent genes and genetic hitchhiking, the results of our model can be easily explained:
(i)
When a mutation produces a silent cooperative trait, it will still act as a defector for a number of generations. Hence, in a world of defectors, it does not get instantly exploited.
(ii)
If this mutation happened in a cellj that has higher fitness due to gene X, πj,Xt, than its neighbors, it will spread through the population by the hitchhiking effect. If it stayed silent till the number of cells carrying the cooperative mutation exceeds a critical cluster size, it will get fixed in the population once it is expressed.
(iii)
If it is expressed too early, while the number of cells carrying the cooperative mutation is still less than the critical cluster size, the damage from exploitation by defectors is most probably much higher than the relative advantage of fitness due to gene X and the mutation dies out.
(iv)
If the mutation happened in a cell that has low or same level of gene X fitness compared to its neighbors, it will not spread. And when it is expressed, it will die out.Figure 14
An illustration of the different scenarios for emergence of a cooperative behavior in a world of defectors in the model. The dot density here represents the background fitness (gene X fitness) of a cell.
(a)
Case 1: a mutation resulting in an instantly expressed cooperative trait,Y=C, gets heavily exploited in a world of defectors and most probably dies out, even if it occurred in a cell with higher gene X fitness than its neighbors
(b)
Case 2: a mutation resulting in a silent trait,Y=S, and the cell carrying the mutation still acts as a defector, S≡D, till time t′ after the mutation. Then, it expresses the cooperative behavior, S≡C. Here, the silent gene did not spread, and it will die out once it is expressed, similar to Case 1
(c)
Case 3: a mutation resulting in a silent trait,Y=S, that happened to be on a cell with higher background fitness (gene X fitness) than its neighbors. The mutation spreads in the system via hitchhiking, and if it gets activated after exceeding a critical cluster size, it stays in the systemIn Figure15, we notice that for the first 500 generations, no silent gene has been “lucky” enough to spread through the population, and so once they are expressed, they get eliminated. Then, a vertical increase in the frequency of cooperators is observed, associated with Figure 16(b). This marks the event that a silent cooperator has hitchhiked a selective sweep of gene X till the number of cells carrying the silent gene S has exceeded the critical cluster size. Cooperators expressed do not die out, and they eventually carry out a successful invasion of the lattice.Figure 15
Evolution of the frequency of cooperators in a50×50 system of defectors, where cooperation arises by silent mutations of gene Y=S.Figure 16
Illustration of a system where mutations of silent cooperator genes (Y=S) occur in a 50×50 system dominated by defector genes (Y=D). Silent genes get activated later at a random point of time t′ to express the cooperative trait.
(a)
Start of simulation
(b)
After 500 generations
(c)
After 4500 generations
(d)
End of simulationTo summarize the factors influencing the emergence of altruistic traits, in a spatial biological system, the probability of a successful invasionPIS for a silent cooperative gene S can be formulated as follows:
(6)PIS=PS⋅PH⋅Pt,with PS as the probability of a mutation producing a silent cooperative gene; PH as the probability that the organism getting the silent mutation has another gene that is going on a selective sweep, thus hitchhiking occurs; and Pt is the probability that the silent gene gets activated only after the elapse of enough time so that its spatial distribution exceeds the critical cluster size by the action of the hitchhiking effect.On the other hand, for comparison, the probability of a successful invasionPIC for an instantly expressed cooperator C is
(7)PIC=PCnwith PC as the probability of a mutation producing an instantly expressed cooperative gene and n as the number of cells required to exceed a critical cluster size of cooperators. If, by chance, mutations happened at the same time interval in cells very close to each other, successful invasion will occur. Hence, the larger the required cluster size for a successful invasion n, the more likely that PIC<<PIS.Finally, it should be noted that the value used for initializing the parameterπj,Xt=0, in this work, it is πj,Xt=0=0 for any cell j in the grid, plays absolutely no part in the results, as long as the same value is used for all the cells in the grid. That is because the invasion probability fΠbit−Πjt is a function of the difference in fitness between two neighboring cells Πbit−Πjt. Hence, the initial value of πj,Xt=0 is not a part of the probability equation. Nonetheless, a mutation in gene X happening to a cell j will lead to a slight increase of its fitness due to gene X, πj,Xt, compared to its neighbors. However, such mutations will spread in the system, and all the cells will always have very close values of fitness due to gene X, making the PD interactions the decisive factor in the competition. In this model, the major influence of beneficial mutations in gene X is simulating the hitchhiking effect. Additionally, such background mutations will also lead to an interesting effect once the cooperators and defectors in the system become spatially segregated, as we will examine in the next section.


### 3.3. Spatial Segregation of Social Traits Is Inherently Unstable

A common feature of spatial cooperation models is that a successful invasion of a cooperative trait results in a state of dynamic equilibrium where cooperator clusters are surrounded by thin patches of defectors [2, 7, 27]. The spatial segregation happens between the two traits not due to any geographical constraints but because of the dynamics of cooperative and selfish interactions. This result makes an intuitive sense as cooperators enjoy the benefits of being surrounded by other cooperators while defectors thrive by exploiting cooperators on the borders of the clusters. So why does our model produce a different fate where the cooperative trait fully invades the system?This is due to beneficial mutations being more likely to arise and get fixed as the population grows [32, 33]. After the initial expansion, a large population of cooperators will have higher frequency of beneficial gene X mutations than the small defector patches. This leads to an accumulated increase in the fitness of cooperators due to gene X compared to defectors, which in turn results in improving the chances for cooperators to expand in the system and, consequently, increase their population, thereby completing a positive-feedback loop. To examine this effect, we calculate the mean of the fitness due to gene X for the defectors, π¯D,Xt, and the cooperators, π¯C,Xt, in this simulation, defined as follows:
(8)π¯D,Xt=∑j=1nDtπj,xtnDt,π¯C,Xt=∑j=1nCtπj,xtnCt,with nDt and nCt as the number of cells in any generation t that are expressing defector and cooperative behavior in the PD interactions, respectively. In Figure 17, as cooperative cells become dominant in the population, the rate of increase of their fitness due to gene X exceeds the defectors’. Till the end of the simulation, cooperators enjoy significantly higher gene X fitness, allowing them to complete a full invasion of the grid. The last defector in the system, d, will have πd,Yt=241+u (see Figure 5(c)). For u=0.09, to be invaded by a neighbor cooperator bi with πbi,Yt=23, the difference in π¯D,Xt and π¯C,Xt will always approach 3 by the end of the simulation.Figure 17
Evolution of the mean of fitness due to gene X for cooperators and defectors during the simulation.It has to be stressed though that this result of a total invasion by the cooperative trait of the system cannot be used to draw general conclusions about the stability of cooperative traits in spatial systems. Since the model is concerned mainly with the problem of emergence of cooperation through rare mutations in a system dominated by defectors, only mutations giving rise to cooperative genes have been considered. In reality, it is likely that mutations can go in both directions, which will generally affect the stability of the cooperative trait in the system depending on the rates of mutations in both directions. The interested reader could refer to the works of [34, 35] for an analysis of this subject.Next, we explore the effects of the activation time of the silent gene,t′, and the gain in fitness due to gene X mutations, μ, on the invasion time of the altruistic trait. To get established in a system, a silent cooperative mutation has to spread via hitchhiking a beneficial mutation, to the extent that when it gets activated, a cooperator cluster larger than the critical size appears. Therefore, the faster a silent mutation gets activated, the less likely that its spatial distribution allows it to invade the system. Also, if the propagation of beneficial mutations is slow within the system, a silent mutation is more likely to get activated before its concentration at a region in the grid reaches the required threshold. The next two sections explore the impact of t′ and μ on the invasion dynamics in the simulation.


### 3.4. Effect of Activation Time

Activation timet′ refers here to how many generations the silent cooperative gene will stay unexpressed. In reality, this could depend on the rarity, or the severity, of the environmental events that could lead to the activation of a particular cooperative mutation. As observed in Figure 18, at low activation times, the invasion time of the cooperative trait is high. When using t′=80 generations, 11% of the simulations did not even reach equilibrium within 105 generations. This number increased to 72% and 100% at t′=40 and t′=0, respectively. On the other hand, as the activation time increases, the invasion time decreases. This is a consequence of the fact that the survival of a cooperative mutation rests on it remaining unexpressed, thus increasing its chances of spreading through the population via hitchhiking till it exceeds the critical cluster size. For a certain critical cluster size, this effect becomes less strong as the activation time becomes excessively long.Figure 18
Relationship between invasion time of a silent cooperative gene and the length of its activation time. For every value of the activation time, the simulation has been run 100 times and the invasion time mean has been calculated. The vertical bars represent the standard deviation of the results.


### 3.5. Effect of Gene X Gain

The second factor that influences the outbreak of a silent cooperative gene via hitchhiking is the speed by which beneficial mutations spread within the system. This will depend on the magnitude of the increase in fitness induced by a mutation and the frequency of interactions between the cells. We study how the change in the fitness gain resulting from gene X mutations,μ, affects the invasion time of the silent cooperative gene. In Figure 19, it is shown that at low values of μ, resulting in slow selective sweeps, the invasion time is high. When using μ=0.2 and μ=0, 79% and 100%, respectively, of the simulations did not reach a successful invasion of the cooperative trait within 105 generations. If beneficial mutations spread slowly within the system, chances are higher that the hitchhiking silent genes will get expressed before reaching a critical distribution which allows it to perform a successful invasion. The higher the μ, the faster genes spread through the population and the lower the invasion time. This effect becomes weaker as μ increases, for fixed activation time and critical cluster size.Figure 19
Relationship between invasion time of a silent cooperative gene and the magnitude of the gain in fitness due to gene X mutations. For every value of gene X gain, the simulation has been run 100 times and the invasion time mean has been calculated. The vertical bars represent the standard deviation of the results.A parameter map visualizing the effects of varying the two parameters on the rate of successful invasions is presented in Figure20.Figure 20
Rate of successful invasions within105 generations in 20 simulations, for different combinations of activation time and gene X gain.


### 3.6. Effect of Moore’s Neighborhood Radius

Biological systems are characterized by their spatial diffusivity, as the effects of interactions between an organism and its neighbors depend on the distance between them. For example, when a living cell shares public good, the surrounding cells benefit from it in different degrees, depending on their distance to the focal cell and the diffusivity of the system. A system with low spatial diffusivity, where cells interact exclusively with their direct neighbors, is expected to offer higher chances for emergence of cooperation than systems with high interaction radius, as a smaller initial cluster of cooperators is needed for cooperation to emerge. Here, this effect is investigated by varying the Moore’s neighborhood radiusr in the simulations. As shown in Figure 21 for fixed parameter values of t′=240 and μ=1.4, the invasion time mean increases as the neighborhood radius increases. A higher neighborhood radius means that the silent gene is more likely to get activated before spreading enough in the system for cooperation to be stable. For comparison, simulations were also carried out for instantly expressed cooperative mutations. For r=1, 29% of the simulations resulted in successful invasions of the cooperative trait within 105 generations. This ratio has dropped to only 5% for r=2 and 0% for larger neighborhood radii. Emergence of cooperation in both cases is significantly easier for smaller values of Moore’s neighborhood radius. However, as r increases, it becomes more difficult for cooperation to emerge by a relatively high number of instantly expressed mutations occurring in close proximity in time and space.Figure 21
Relationship between invasion time of a silent cooperative gene and Moore’s neighborhood radius. Simulations are run for 100 times for each point with the vertical bars representing the standard deviation.We end this section with a short discussion of how silent cooperative genes can influence our understanding of natural selection.


### 3.7. Selection for Quietness

In [11], it is shown that cooperation can be stable in the setting of spatial PD interactions. Reference [7] has constructed a similar model in order to investigate the invasion dynamics of a cooperative trait in spatial PD and snowdrift games. There, it has been shown that in the case of the PD game, the model reaches a state of spatial segregation between cooperators and defectors with cooperators breaking up into isolated clusters in order to reduce the exploitation effects from defectors. In Section 3.3, we have shown that including the occurrence of background beneficial mutations into the model leads to the instability of the spatial segregation between traits at the long term. More importantly, a common feature in the works of [2, 7] is that for a cooperative trait to survive and expand in the system, a critically sized cluster of cells carrying this trait has to exist at the beginning. Both works have thoroughly investigated the invasion dynamics, the stability of a cooperative trait, and the resulting spatial patterns in spatial PD systems. This paper focuses on the problem of how an initial cluster of cooperators could appear by random mutations in the first place.Our model showed that the delayed expression of a cooperative gene could improve its chances of a successful invasion. Also, the longer the time till switching on the gene, the better. Similar observations have been made by [16] when studying how quorum sensing can enable a cooperative bacterial strain to dominate a biofilm, by delaying the expression of the cooperative gene till enough segregation occurs. Furthermore, it has been shown that in case of strong competition (highly mixed biofilms), only bacterial strains with a high quorum-sensing threshold can dominate the biofilm as segregation occurs later in time. Natural selection in this case selects for the strains that delay expression of their cooperative trait.In a general setting, the larger the critical size of a cooperator cluster required for a successful invasion, the stronger the need for a delayed expression of the cooperator gene. In other words, natural selection here acts to select the cooperative genes which require more aggressive change in the environment to be expressed.A typical view of natural selection is that it acts on the currently expressed phenotypes. We fix a point in time, and we look at the genes expressing phenotypes that have higher fitness compared to their alleles. Based on that, one can predict that the frequency of those more adapted alleles will increase in the environment as the biological system evolves. On the other hand, natural selection, in our work and in [16], selects genes not on their expressed phenotypes but on their ability to not express them. That is, it favors genes that have devised the most devious ways to endure a possibly long hitchhiking trip.


## 3.1. Emergence of Cooperation in a Biological System

First, to illustrate the problem of emergence of cooperation, we reproduce earlier results in literature by examining the effect of the initial cooperators’ cluster size on their fate, with no mutations in the simulation. In Figures8 and 9, we seed a 50×50 lattice with clusters of cooperators of sizes 9 and 25, respectively, which are cells with gene Y being cooperative (i.e., Y=C). The rest of the positions in the lattice are occupied by defectors, hence gene Y being defective (i.e., Y=D), while the value of the normalized PD variable u is set to u=0.09. In Figures 8 and 9, the interactions among cooperators are not sufficient to compensate for the losses suffered by the exploitation from selfish agents and the cooperative gene goes extinct rapidly.Figure 8
Evolution of cooperator frequency (black) when starting with an initial3×3 cooperator cluster within a 50×50 grid of defectors (yellow).
(a)
Start of simulation
(b)
End of simulation
(c)
Cooperator frequencyFigure 9
Evolution of cooperator frequency (black) when starting with an initial5×5 cooperator cluster within a 50×50 grid of defectors (yellow).
(a)
Start of simulation
(b)
End of simulation
(c)
Cooperator frequencyOn the other hand, when starting with large cluster sizes of6×6 and 8×8 cooperators, as shown in Figures 10 and 11, the initial cooperator population undergoes aggressive expansion at the beginning. The number of cells with Y=C increases quadratically till it reaches a relatively stable value, around 65% of the total population. Simulations are stopped after 105 generations. A dynamic equilibrium state is reached where compact clusters of cooperators are surrounded by stripes of defectors, in accordance with previous research [2, 7, 11, 27]. These results confirm that in a spatial system, cooperators in a PD setting can dominate a world of defectors as their interactions among each other compensate losses due to encounters with cheaters. However, this will happen only if they started from a large enough initial cluster.Figure 10
Evolution of cooperator frequency (black) when starting with an initial6×6 cooperator cluster within a 50×50 grid of defectors (yellow).
(a)
Start of simulation
(b)
End of simulation
(c)
Cooperator frequencyFigure 11
Evolution of cooperator frequency (black) when starting with an initial8×8 cooperator cluster within a 50×50 grid of defectors (yellow).
(a)
Start of simulation
(b)
End of simulation
(c)
Cooperator frequencyWe proceed to study the problem of emergence of cooperation in a world of defectors. The50×50 lattice is initialized by filling it completely with defectors, agents with Y=D, as in Figure 12(a). However, mutations are allowed to occur in the model. The probability that a cell with gene Y=D will mutate to Y=C is PYD→C=1×10−4. Additionally, the probability of a beneficial mutation in gene X, PXμ, is set to be PXμ=1×10−4, with μ=1.0. As observed in Figures 12 and 13, the instantly expressed gene Y=C fails to make a successful invasion of the system. This has been expected as a cooperator in a sea of defectors is very prone to exploitation as not enough interactions occur between cooperators to out-weight exploitation losses. Also, even if a cooperative mutation occurred in a cell with slightly higher X than its neighbors, the losses suffered from exploitation by defectors will be much larger than the gain in fitness due to gene X. Finally, the probability that cooperative mutations will coappear in a narrow spatial range such that they form a cluster of altruists is extremely small.Figure 12
Illustration of a system where mutations of instantly expressed cooperator genes (Y=C) occur in a 50×50 system dominated by defector genes (Y=D).
(a)
Start of simulation
(b)
After 5 × 105 generations
(c)
End of simulationFigure 13
Evolution of the frequency of cooperators in a50×50 system of defectors, where cooperation arises by instantly expressed mutations of gene Y=C.Do mutations of a silent cooperator gene have better chances for invading the system? We test that by allowing, with a probability ofPYD→S=1×10−4, that a cell with gene Y=D will mutate into Y=S, which in turn will get activated after passing of t′ generations after the mutation. For this simulation, t′ is assumed to be exponentially distributed, t′∼Expλ, with a mean value of 200 generations. The simulation is again initialized by filling the lattice with defectors. In the early beginning, no cooperators at all appear, as mutations occurring in the system are still inactivated. After that, the system evolves in a way similar to the previous case of instantly expressed cooperator genes, with cooperators appearing isolated from each other and doomed to extinction. Then, after around 500 generations in this simulation, a sudden outbreak of cooperation appears with a cluster size that is large enough to eventually perform a successful invasion of the system. Thus, cooperation prevailed here. But why exactly did cooperator genes with delayed expression succeed while instantly expressed genes failed? We take a closer look at the two key features of the model which allowed these results to occur in the next section.


## 3.2. Silent Genes and Genetic Hitchhiking

Epigenetics refer to changes in gene expression that are not resulting from changes in the DNA code [28]. A gene could produce a phenotype only when activated by certain environmental conditions. Also, it has been found that a gene could be tuned by an inheritance mechanism to get expressed only after a number of generations, a phenomenon known as transgenerational gene silencing [29, 30].In our model, when a silent gene appears by a mutation, it remains unexpressed for a random period associated with the physical nature of each mutation. Then, after passing oft′ generations, it gets activated at all cells in the system carrying this mutation. This silent interval is however essential for cooperation to break out as the gene can still spread through the population by a mechanism known as genetic hitchhiking. Genetic hitchhiking is the process by which a gene that offers no evolutionary advantage to its carrier can spread through a population due to being associated with another gene that is under a strong selection [31]. Due to this phenomenon, the frequency of a neutral gene can increase within the population if it is on the same DNA chain of another beneficial gene. In the model, mutations of gene X play the role of the beneficial gene. A cell with higher fitness value due to the contribution of gene X can slowly invade the population. As illustrated in Figure 14, by combining the concepts of silent genes and genetic hitchhiking, the results of our model can be easily explained:
(i)
When a mutation produces a silent cooperative trait, it will still act as a defector for a number of generations. Hence, in a world of defectors, it does not get instantly exploited.
(ii)
If this mutation happened in a cellj that has higher fitness due to gene X, πj,Xt, than its neighbors, it will spread through the population by the hitchhiking effect. If it stayed silent till the number of cells carrying the cooperative mutation exceeds a critical cluster size, it will get fixed in the population once it is expressed.
(iii)
If it is expressed too early, while the number of cells carrying the cooperative mutation is still less than the critical cluster size, the damage from exploitation by defectors is most probably much higher than the relative advantage of fitness due to gene X and the mutation dies out.
(iv)
If the mutation happened in a cell that has low or same level of gene X fitness compared to its neighbors, it will not spread. And when it is expressed, it will die out.Figure 14
An illustration of the different scenarios for emergence of a cooperative behavior in a world of defectors in the model. The dot density here represents the background fitness (gene X fitness) of a cell.
(a)
Case 1: a mutation resulting in an instantly expressed cooperative trait,Y=C, gets heavily exploited in a world of defectors and most probably dies out, even if it occurred in a cell with higher gene X fitness than its neighbors
(b)
Case 2: a mutation resulting in a silent trait,Y=S, and the cell carrying the mutation still acts as a defector, S≡D, till time t′ after the mutation. Then, it expresses the cooperative behavior, S≡C. Here, the silent gene did not spread, and it will die out once it is expressed, similar to Case 1
(c)
Case 3: a mutation resulting in a silent trait,Y=S, that happened to be on a cell with higher background fitness (gene X fitness) than its neighbors. The mutation spreads in the system via hitchhiking, and if it gets activated after exceeding a critical cluster size, it stays in the systemIn Figure15, we notice that for the first 500 generations, no silent gene has been “lucky” enough to spread through the population, and so once they are expressed, they get eliminated. Then, a vertical increase in the frequency of cooperators is observed, associated with Figure 16(b). This marks the event that a silent cooperator has hitchhiked a selective sweep of gene X till the number of cells carrying the silent gene S has exceeded the critical cluster size. Cooperators expressed do not die out, and they eventually carry out a successful invasion of the lattice.Figure 15
Evolution of the frequency of cooperators in a50×50 system of defectors, where cooperation arises by silent mutations of gene Y=S.Figure 16
Illustration of a system where mutations of silent cooperator genes (Y=S) occur in a 50×50 system dominated by defector genes (Y=D). Silent genes get activated later at a random point of time t′ to express the cooperative trait.
(a)
Start of simulation
(b)
After 500 generations
(c)
After 4500 generations
(d)
End of simulationTo summarize the factors influencing the emergence of altruistic traits, in a spatial biological system, the probability of a successful invasionPIS for a silent cooperative gene S can be formulated as follows:
(6)PIS=PS⋅PH⋅Pt,with PS as the probability of a mutation producing a silent cooperative gene; PH as the probability that the organism getting the silent mutation has another gene that is going on a selective sweep, thus hitchhiking occurs; and Pt is the probability that the silent gene gets activated only after the elapse of enough time so that its spatial distribution exceeds the critical cluster size by the action of the hitchhiking effect.On the other hand, for comparison, the probability of a successful invasionPIC for an instantly expressed cooperator C is
(7)PIC=PCnwith PC as the probability of a mutation producing an instantly expressed cooperative gene and n as the number of cells required to exceed a critical cluster size of cooperators. If, by chance, mutations happened at the same time interval in cells very close to each other, successful invasion will occur. Hence, the larger the required cluster size for a successful invasion n, the more likely that PIC<<PIS.Finally, it should be noted that the value used for initializing the parameterπj,Xt=0, in this work, it is πj,Xt=0=0 for any cell j in the grid, plays absolutely no part in the results, as long as the same value is used for all the cells in the grid. That is because the invasion probability fΠbit−Πjt is a function of the difference in fitness between two neighboring cells Πbit−Πjt. Hence, the initial value of πj,Xt=0 is not a part of the probability equation. Nonetheless, a mutation in gene X happening to a cell j will lead to a slight increase of its fitness due to gene X, πj,Xt, compared to its neighbors. However, such mutations will spread in the system, and all the cells will always have very close values of fitness due to gene X, making the PD interactions the decisive factor in the competition. In this model, the major influence of beneficial mutations in gene X is simulating the hitchhiking effect. Additionally, such background mutations will also lead to an interesting effect once the cooperators and defectors in the system become spatially segregated, as we will examine in the next section.


## 3.3. Spatial Segregation of Social Traits Is Inherently Unstable

A common feature of spatial cooperation models is that a successful invasion of a cooperative trait results in a state of dynamic equilibrium where cooperator clusters are surrounded by thin patches of defectors [2, 7, 27]. The spatial segregation happens between the two traits not due to any geographical constraints but because of the dynamics of cooperative and selfish interactions. This result makes an intuitive sense as cooperators enjoy the benefits of being surrounded by other cooperators while defectors thrive by exploiting cooperators on the borders of the clusters. So why does our model produce a different fate where the cooperative trait fully invades the system?This is due to beneficial mutations being more likely to arise and get fixed as the population grows [32, 33]. After the initial expansion, a large population of cooperators will have higher frequency of beneficial gene X mutations than the small defector patches. This leads to an accumulated increase in the fitness of cooperators due to gene X compared to defectors, which in turn results in improving the chances for cooperators to expand in the system and, consequently, increase their population, thereby completing a positive-feedback loop. To examine this effect, we calculate the mean of the fitness due to gene X for the defectors, π¯D,Xt, and the cooperators, π¯C,Xt, in this simulation, defined as follows:
(8)π¯D,Xt=∑j=1nDtπj,xtnDt,π¯C,Xt=∑j=1nCtπj,xtnCt,with nDt and nCt as the number of cells in any generation t that are expressing defector and cooperative behavior in the PD interactions, respectively. In Figure 17, as cooperative cells become dominant in the population, the rate of increase of their fitness due to gene X exceeds the defectors’. Till the end of the simulation, cooperators enjoy significantly higher gene X fitness, allowing them to complete a full invasion of the grid. The last defector in the system, d, will have πd,Yt=241+u (see Figure 5(c)). For u=0.09, to be invaded by a neighbor cooperator bi with πbi,Yt=23, the difference in π¯D,Xt and π¯C,Xt will always approach 3 by the end of the simulation.Figure 17
Evolution of the mean of fitness due to gene X for cooperators and defectors during the simulation.It has to be stressed though that this result of a total invasion by the cooperative trait of the system cannot be used to draw general conclusions about the stability of cooperative traits in spatial systems. Since the model is concerned mainly with the problem of emergence of cooperation through rare mutations in a system dominated by defectors, only mutations giving rise to cooperative genes have been considered. In reality, it is likely that mutations can go in both directions, which will generally affect the stability of the cooperative trait in the system depending on the rates of mutations in both directions. The interested reader could refer to the works of [34, 35] for an analysis of this subject.Next, we explore the effects of the activation time of the silent gene,t′, and the gain in fitness due to gene X mutations, μ, on the invasion time of the altruistic trait. To get established in a system, a silent cooperative mutation has to spread via hitchhiking a beneficial mutation, to the extent that when it gets activated, a cooperator cluster larger than the critical size appears. Therefore, the faster a silent mutation gets activated, the less likely that its spatial distribution allows it to invade the system. Also, if the propagation of beneficial mutations is slow within the system, a silent mutation is more likely to get activated before its concentration at a region in the grid reaches the required threshold. The next two sections explore the impact of t′ and μ on the invasion dynamics in the simulation.


## 3.4. Effect of Activation Time

Activation timet′ refers here to how many generations the silent cooperative gene will stay unexpressed. In reality, this could depend on the rarity, or the severity, of the environmental events that could lead to the activation of a particular cooperative mutation. As observed in Figure 18, at low activation times, the invasion time of the cooperative trait is high. When using t′=80 generations, 11% of the simulations did not even reach equilibrium within 105 generations. This number increased to 72% and 100% at t′=40 and t′=0, respectively. On the other hand, as the activation time increases, the invasion time decreases. This is a consequence of the fact that the survival of a cooperative mutation rests on it remaining unexpressed, thus increasing its chances of spreading through the population via hitchhiking till it exceeds the critical cluster size. For a certain critical cluster size, this effect becomes less strong as the activation time becomes excessively long.Figure 18
Relationship between invasion time of a silent cooperative gene and the length of its activation time. For every value of the activation time, the simulation has been run 100 times and the invasion time mean has been calculated. The vertical bars represent the standard deviation of the results.


## 3.5. Effect of Gene X Gain

The second factor that influences the outbreak of a silent cooperative gene via hitchhiking is the speed by which beneficial mutations spread within the system. This will depend on the magnitude of the increase in fitness induced by a mutation and the frequency of interactions between the cells. We study how the change in the fitness gain resulting from gene X mutations,μ, affects the invasion time of the silent cooperative gene. In Figure 19, it is shown that at low values of μ, resulting in slow selective sweeps, the invasion time is high. When using μ=0.2 and μ=0, 79% and 100%, respectively, of the simulations did not reach a successful invasion of the cooperative trait within 105 generations. If beneficial mutations spread slowly within the system, chances are higher that the hitchhiking silent genes will get expressed before reaching a critical distribution which allows it to perform a successful invasion. The higher the μ, the faster genes spread through the population and the lower the invasion time. This effect becomes weaker as μ increases, for fixed activation time and critical cluster size.Figure 19
Relationship between invasion time of a silent cooperative gene and the magnitude of the gain in fitness due to gene X mutations. For every value of gene X gain, the simulation has been run 100 times and the invasion time mean has been calculated. The vertical bars represent the standard deviation of the results.A parameter map visualizing the effects of varying the two parameters on the rate of successful invasions is presented in Figure20.Figure 20
Rate of successful invasions within105 generations in 20 simulations, for different combinations of activation time and gene X gain.


## 3.6. Effect of Moore’s Neighborhood Radius

Biological systems are characterized by their spatial diffusivity, as the effects of interactions between an organism and its neighbors depend on the distance between them. For example, when a living cell shares public good, the surrounding cells benefit from it in different degrees, depending on their distance to the focal cell and the diffusivity of the system. A system with low spatial diffusivity, where cells interact exclusively with their direct neighbors, is expected to offer higher chances for emergence of cooperation than systems with high interaction radius, as a smaller initial cluster of cooperators is needed for cooperation to emerge. Here, this effect is investigated by varying the Moore’s neighborhood radiusr in the simulations. As shown in Figure 21 for fixed parameter values of t′=240 and μ=1.4, the invasion time mean increases as the neighborhood radius increases. A higher neighborhood radius means that the silent gene is more likely to get activated before spreading enough in the system for cooperation to be stable. For comparison, simulations were also carried out for instantly expressed cooperative mutations. For r=1, 29% of the simulations resulted in successful invasions of the cooperative trait within 105 generations. This ratio has dropped to only 5% for r=2 and 0% for larger neighborhood radii. Emergence of cooperation in both cases is significantly easier for smaller values of Moore’s neighborhood radius. However, as r increases, it becomes more difficult for cooperation to emerge by a relatively high number of instantly expressed mutations occurring in close proximity in time and space.Figure 21
Relationship between invasion time of a silent cooperative gene and Moore’s neighborhood radius. Simulations are run for 100 times for each point with the vertical bars representing the standard deviation.We end this section with a short discussion of how silent cooperative genes can influence our understanding of natural selection.


## 3.7. Selection for Quietness

In [11], it is shown that cooperation can be stable in the setting of spatial PD interactions. Reference [7] has constructed a similar model in order to investigate the invasion dynamics of a cooperative trait in spatial PD and snowdrift games. There, it has been shown that in the case of the PD game, the model reaches a state of spatial segregation between cooperators and defectors with cooperators breaking up into isolated clusters in order to reduce the exploitation effects from defectors. In Section 3.3, we have shown that including the occurrence of background beneficial mutations into the model leads to the instability of the spatial segregation between traits at the long term. More importantly, a common feature in the works of [2, 7] is that for a cooperative trait to survive and expand in the system, a critically sized cluster of cells carrying this trait has to exist at the beginning. Both works have thoroughly investigated the invasion dynamics, the stability of a cooperative trait, and the resulting spatial patterns in spatial PD systems. This paper focuses on the problem of how an initial cluster of cooperators could appear by random mutations in the first place.Our model showed that the delayed expression of a cooperative gene could improve its chances of a successful invasion. Also, the longer the time till switching on the gene, the better. Similar observations have been made by [16] when studying how quorum sensing can enable a cooperative bacterial strain to dominate a biofilm, by delaying the expression of the cooperative gene till enough segregation occurs. Furthermore, it has been shown that in case of strong competition (highly mixed biofilms), only bacterial strains with a high quorum-sensing threshold can dominate the biofilm as segregation occurs later in time. Natural selection in this case selects for the strains that delay expression of their cooperative trait.In a general setting, the larger the critical size of a cooperator cluster required for a successful invasion, the stronger the need for a delayed expression of the cooperator gene. In other words, natural selection here acts to select the cooperative genes which require more aggressive change in the environment to be expressed.A typical view of natural selection is that it acts on the currently expressed phenotypes. We fix a point in time, and we look at the genes expressing phenotypes that have higher fitness compared to their alleles. Based on that, one can predict that the frequency of those more adapted alleles will increase in the environment as the biological system evolves. On the other hand, natural selection, in our work and in [16], selects genes not on their expressed phenotypes but on their ability to not express them. That is, it favors genes that have devised the most devious ways to endure a possibly long hitchhiking trip.


## 4. Conclusion

Cooperation plays a central role in countless biological phenomena, and explaining its emergence and preservation in living systems remains a major challenge. Spatial competition dynamics have provided an intuitive answer to the problem of how cooperation can persist in a biological system. Since individuals are more likely to interact with their neighbors, cooperative cells in a spatial PD game can dominate the system by existing in clusters. Consequently, the high frequency of mutually beneficial interactions among altruists can balance the exploitation by defectors on the boundaries of the cooperative segments. This explanation however leads us to consider the fundamental problem of the lonely cooperator. Rare mutations producing individuals with simple cooperative traits, not involving memory or communication, in a system dominated by defectors are doomed to extinction, as they are prone to heavy exploitation by their neighbors. In hostile environments, emergence of such traits is highly unlikely.In this paper, we propose the idea that silent genes could provide a route for emergence of cooperation in biological systems. A silent mutation is one that does not get immediately expressed. Instead, it only gets activated at a random point of the future, possibly due to specific environmental conditions. A rare mutant carrying a silent cooperative trait will still act as a defector, making it safe from exploitation, till its latent altruistic trait gets activated. Subsequently, a silent gene has a chance to spread through the biological system, which occurs in our model by hitchhiking to another beneficial gene under positive selection. Afterwards, when a change in environment happens that can trigger this specific mutation, cells carrying this gene could happen to exist in the required spatial distribution to survive, and even invade, the system. In short, two established genetic concepts, silent genes and genetic hitchhiking, could explain how a simple cooperative trait emerges in a biological system.It is also interesting to note that our model suggests that spatial segregation between social traits, an equilibrium state frequently produced by spatial cooperation models, could be evolutionary unstable. Two types of the same organism which are spatially segregated and genetically identical except in a single social trait can be assumed to have the same mutation frequency. Since the rate of beneficial mutations in a population of any type is proportional to its size, and the population size itself expands as the fitness of its individuals increases, the coexistence of two competing simple social traits in a biological system is bound to be an unstable equilibrium point. That said, it has to be emphasized that the fourth-mentioned result is qualified by the noninclusion of mutations giving rise to the defector trait in the model, which will significantly affect the general problem of the stability of cooperation in biological systems. Additionally, it should be also highlighted that this model did not include background deleterious mutations;μ is always positive. However, this will not affect the results as such mutations will not spread in the system. Only beneficial mutations play a significant role as they carry out the hitchhiking process. Also, it is assumed in this work that when the environment reaches a state that could trigger a silent gene, it gets synchronously activated in all the organisms carrying it in the population. This is clearly a simplifying assumption. In reality, the environment never changes at exactly the same degree for all the population. Nevertheless, it is natural to expect that organisms which are spatially close to each other will experience similar environmental conditions. This could be enough for cooperation to emerge as the trait needs to be expressed only by a cluster of organisms at the same time interval to be stable.Finally, we would like to highlight that in the model presented in this paper, for simplicity, genes are assumed to only transfer vertically, from a single parent to a single offspring. In real life, genes will also transfer horizontally. In bacteria, this happens by the means of horizontal gene transfer, and in sexual organisms, this happens by the process of sexual reproduction. Hence, for a silent mutation to hitchhike a gene undergoing a selective sweep, it needs to be physically close to the beneficial gene on the DNA. It has been shown that hitchhiking events can be detected via statistical analysis of genetic variation along chromosomes [36]. Therefore, when considering real-life situations, our model predicts that cooperative genes are more likely to be close in the DNA chain to “competitive” genes, ones that may have gone through a selective sweep in the past, than being close to other cooperative/public good genes or silent segments of the DNA.



---

*Source: 2082037-2018-09-18.xml*