BACKGROUND
when cooperation is costly to the individual but its benefits are equally shared in a group, one would expect progressively more selfish behaviour to be selected.
this indeed happens when interactions between cooperative and selfish individuals are random.
instead, cooperation can be selected in a population when the interactions among cooperative and selfish individuals are structured, be it genetically, spatially or socially .
population structure favours cooperation when it allows for cooperators to be in contact with each other more frequently than with selfish individuals  <cit> .
the result is a  stable equilibrium configuration in which cooperators persist indefinitely, and selfish individuals may co-exist.
for instance, this can be due to spatial clustering of cooperators  <cit> , or to an inherent structure of the interaction network  <cit> .
then, the conditions for selfish individuals invading and overtaking a group of cooperators represent the limit to the stability of the solutions found.
in general, these conditions state that higher costs of cooperation increase the selective advantage of selfish individuals.
while spatial structure alone can favour cooperators due to population viscosity, a growing body of experimental and theoretical work indicates that self-organised spatial patterns may have profound and complex effects on cooperative interactions, due to emergent heterogeneities in the local distributions and densities of cooperators and selfish individuals .
for instance, as a population of cooperators invade empty space, its expansion front can be enriched in altruistic individuals, while selfish individuals lag behind  <cit> .
alternatively, mutualistic interactions that are favoured in a resident population can be automatically broken on such expansion front  <cit> .
furthermore, spatial self-organisation can sort cooperative strains from selfish ones , thus limiting the spread of the latter.
the customary approach to study the stability of cooperation  consists of fixing the strategies of the interacting cooperative and selfish individuals, and analyse the population dynamics of the system .
because such pre-determined strategies do not mutate over time, their evolutionary stability remains unexplored.
although exceptions to this approach exist , little is known about what strategies evolve  and how they feedback on the spatial self-organisation of a population, even though it is clear that spatial self-organisation affects the population dynamics of cooperative traits .
here, we seek to study this feed-back between evolution and self-organisation with a minimal model where individuals can evolve the degree of cooperation in a spatially extended system.
we model the cooperative trait in terms of public good production , and we let the amount of public good produced mutate in a continuous fashion.
thus, we can study the long term evolutionary dynamics of cooperation without preconceiving the extent to which individuals cooperate or defect.
as we will show, selfish individuals that produce zero public good evolve readily at higher costs, and quickly invade cooperators.
rather than leading to global extinction, this enables cooperators to thrive and selects for a higher degree of public good production over evolutionary time scales.
RESULTS
the model is a straightforward implementation of a population in which individuals replicate depending on the amount of public good produced in their close neighbourhood.
individuals are embedded on a lattice.
they may reproduce, die or move .
competition for reproduction into neighbouring empty nodes is based on fitness, calculated as the difference between benefits and costs .
an individual benefits from the public good produced in its neighbourhood, but pays a cost for producing it.
thus, public good production is a cooperative trait.
we assume that reproductive success is solely based on cooperation, so that individuals do not reproduce if public good in their neighbourhood is insufficient.
mutations slightly change the offspring’s production rate .
fig.
<dig> the model.
a the world is a square lattice with connectivity k= <dig>  and wrapped boundary conditions.
individuals produce public good at rate p, shared in equal parts  among all neighbouring nodes and self.
b individuals compete for reproduction into an adjacent empty spot.
probability of reproduction depends on fitness f
i, which is the difference between benefits and costs.
the sum of the public good an individual collects from itself and the neighbours, if any, is p
sum=/ <dig>  and confers a benefit b·p
sum.
individuals pay a fitness cost proportional to the public good they produce: c·p.
successful reproduction yields a copy of the selected individual.
mutations occur with probability μ and change the public good production rate by a small random number chosen uniformly in the interval .
individuals have a small probability k
move to move to a random adjacent node, and can die with probability k
death, leaving the node empty
high cost leads to the evolution of larger public good production
we set the benefits per unit of public good b= <dig>  and we let the spatial self-organisation and the evolutionary dynamics unfold under different costs c.
when costs are much smaller than benefits , the public good production steadily increases because an individuals’ own production increases its fitness, rather than decreasing it.
moreover, because replication is a local process, mutants with higher than average production rates benefit from each other due to limited dispersal , outcompeting more selfish lineages.
thus, cooperation is maximised in the long run.
for increasing costs, public good production suddenly drops .
in this regime, the clustering advantage of cooperators is insufficient and more selfish individuals replicate the most because, by producing less, they pay a lower cost.
eventually, public good production stabilises at the minimum value needed for survival.
fig.
<dig> an increase in costs results in an evolutionary increase in cooperation.
a evolutionary steady state distribution  and mean  of public good production are plotted for different values of costs .
parameters: benefit  b= <dig>  k
death= <dig> , k
move= <dig> , μ= <dig> , δ= <dig> .
the maximum public good production is set to p
max= <dig>  b snapshots of the lattice at evolutionary steady state.
colour coding depends on public good production rate.
white is background.
lattice dimensions used for the simulations from left to right:  <dig>   <dig>   <dig>   <dig>
strikingly, further increasing costs leads to an increase in cooperation .
the distribution of public good production is bimodal at evolutionary steady state, with most of the population having higher rates of public good production and a minority producing almost no public good at all.
the long term evolutionary dynamics of cooperation at high cost
following the evolutionary and the spatial dynamics of a single case elucidates why cooperation increases and persists for higher costs .
when we initialise the system with highly producing individuals, public good production decreases rapidly due to strong selection for selfishness .
where public good production drops below the minimum for survival, large patches of individuals go extinct .
the surviving individuals can expand into the empty space .
as the expansion progresses, a selfish and cooperative lineage separate from each other .
the selfish strain evolves to zero public good production, becoming incapable of autonomous persistence and relying on the public good produced by cooperators for survival.
fig.
<dig> evolutionary dynamics of public good production.
a at each time point the distribution of public good production in the lattice is plotted as a heat map.
b-g snapshots of the lattice at subsequent time steps .
costs  c= <dig> , other parameters and colour coding in the snapshot as in fig.
<dig>  lattice size = <dig> .
time units are monte carlo steps.
see also additional file 1: section  <dig> for the full snapshots, and movie at  <cit>
while the two strains differentiate from each other, they organise spatially to form travelling waves .
cooperators constitute the front of a wave, and expand into empty space by replicating into it ; selfish individuals invade those cooperators, and constitute the back of the wave.
selfish individuals leave empty space behind a wave after they die, causing the semblance of movement .
the progression of a wave, however, happens on a time scale that is much longer than the life time of an individual, which in turn experiences a fairly constant environment throughout its life.
because waves persist longer than individuals, they can integrate information over several generations.
fig.
<dig> the movement of a travelling wave.
magnification  of the same portion of the lattice at  <dig> time steps distance.
a the wave-front, composed of cooperative individuals , b the back, composed of selfish individuals .
the snapshots are from the same simulation run as in fig.
<dig>  at time steps  <dig>   and  <dig>  .
the dashed line in the right pane marks approximately the position of the wave in the left pane.
colizzi and hogeweg  <cit>  shows one such travelling wave in a video
new waves are “born” from the collisions of older waves.
as cooperators on the front evolve to larger public good production, waves become larger .
the formation of spatial patterns allows the populations in the system to persist indefinitely despite selfish individuals continuously invading cooperators , provided that the lattice is much larger than the spatial patterns .
spatial population dynamics of cooperators and selfish individuals
spatial self-organisation drives the evolution of cooperation in the system.
when spatial patterns are destroyed, e.g.
by mixing, only selfish individuals are selected and public good production decreases, leading to global extinction , in accordance with the result that random interactions favour selfish behaviour.
to unravel the interplay between the two lineages and their spatial organisation, we analysed the spatial population dynamics for cooperative and selfish individuals separately.
to this end, we shaped the lattice into a long, narrow strip, and set the mutation rate to zero .
cooperators expand faster into empty space when they produce more public good, and slower when costs are higher .
when two clustered populations compete at the expansion front , the one with the largest public good production wins because, by replicating faster, it occupies space before the competing one and eventually overtakes the entire wave front.
fig.
<dig> invasion rates and spatial dynamics of competition for cooperative and selfish individuals.
a cooperators’ invasion rate increases with larger public good production and lower costs; selfish invasion rate increases with increasing costs and it is insensitive to the production rates of the cooperators that support them.
the invasion rates for cooperators invading empty space  and selfish individuals invading a population of cooperators  was measured for each combination of cost and public good production rate  in  <dig>  b a population with lower public good production rate out-competes one with larger production in the back of a wave.
parameters: cooperators  p
yellow= <dig> , selfish individuals  p
blue= <dig>  and p
magenta= <dig> .
c the population with a larger production rate out-competes the one with lower production at the front of a wave.
parameters: p
purple= <dig> , p
green= <dig> .
for both a and b, μ= <dig>  b= <dig>  c= <dig> , background in white.
other parameters as in fig.
2
the replication rate of selfish individuals invading a population of cooperators is higher when the cost of public good production is higher, and it is insensitive to how much public good is produced .
clearly, when two strains compete in the back of a wave, the winner is the more selfish one .
the picture emerging from these experiments is that different selection pressures operate depending on the spatial context: a population expanding into empty space  is selected for higher degrees of cooperation , competition in the back  selects for more selfishness.
importantly, even though costs, benefits and fitness function are the same, spatial pattern formation automatically segregates these two opposing evolutionary pressures to spatially different contexts so that they do not balance each other: hence the evolution of a cooperative and a selfish lineage.
in the full system, cooperators and selfish individuals are ecologically and evolutionarily interdependent .
the two lineages establish an evolutionary feedback mediated by their spatial organisation.
empty space is generated by selfish individuals after invading a population of cooperators.
therefore, the condition for increasing public good production, i.e.
the availability of continuous empty space, is mediated by the invasion dynamics of the selfish lineage.
with higher costs, selfish individuals propagate faster, and space is left empty at a higher rate.
the larger the empty space, the more cooperators can increase public good production.
this evolutionary feedback reaches an evolutionary steady state because highly producing cooperators reduce the empty space faster upon faster expansions .
we checked the long term stability of the steady state .
altogether, cooperation evolves to a higher degree for higher costs due to an emergent feedback between self-organised interaction structures  and the evolution of the individuals composing them.
robustness to parameter change
our results are robust when death rates, movement rates and benefits are changed, provided that benefits-to-costs ratio is maintained .
at lower costs, individuals directly benefiting themselves with their own public good are sometimes dubbed weak altruists  <cit>  , with k+1= <dig> the connectivity of the lattice including self), whereas they are considered strong altruists when their public good is only shared among others .
in our spatial model, we observed no qualitative difference in the evolutionary dynamics when individuals did or did not benefit from their own public good .
this could be expected because individuals’ own payoffs at high costs are negative in both models.
indeed, in both cases the evolving populations underwent speciation of a selfish and a cooperative lineage.
weak or strong altruism do make a difference in the corresponding well-mixed systems, where strong altruists go extinct at lower costs, while weak altruists maximise public good production .
discussion
it is known that during population range expansion, cooperation can be promoted on the front of the expansion range  <cit> .
one could argue that in these models cooperation could evolve only as long as empty space is available, and should eventually be out-competed globally by selfish strategies when the invasion dynamics reach an end.
here we have shown that selfish individuals provide the empty space to allow continuous expansion within a limited area.
generally, stable solutions to the problem of cooperation are based on the condition that at equilibrium selfish individuals do not locally invade a population of cooperators.
for instance, it is well known that cooperation can be stabilised in spatially extended systems, as cooperators cluster and segregate selfish individuals to the edges of those clusters  <cit> .
a side effect of these solutions is that higher costs undermines the stability of such clusters.
in the parameter region where travelling waves do not form  our work is in agreement with those results in that there exists an inverse relationship between costs and cooperation , and in particular with the heuristics  <cit>  that cooperation evolves when the benefit-to-cost ratio is larger than the connectivity of the lattice .
however, we have shown that a novel class of solutions exists at high costs, where large degrees of cooperation are maintained in a locally out-of-equilibrium fashion, with selfish individuals always successfully invading cooperators and setting the stage for the evolutionary increase and the global stability of cooperation.
we conclude that spatial self-organisation can reverse the relationship between costs and cooperation, thus extending the evolutionary viability of cooperation to higher costs.
our results rest on two assumptions: population size can vary and some degree of cooperation is necessary for reproduction.
variable population size is obviously realistic, even though it is not often included in evolutionary models of cooperation.
although the assumption of necessary cooperation is not always met, it is reasonable in several cases.
examples in microbiology include, cooperative protection or cooperative virulence in bacterial infections  <cit> ; invertase production in yeasts while growing on sucrose  <cit> ; siderophore production in iron-limited environments , cooperative secretion of digestive enzyme in microbial hunting  <cit> .
outside the microbial world, situations where our model may apply are e.g.
dangerous behaviours in cooperative nest defence  <cit> , and replication in trans in prebiotic evolution  <cit> .
two recent studies have come to conclusions that at first sight are similar to ours  <cit> .
by making the assumption, as we do, that the lack of cooperation leads to death, they observe  static spatial patterns in which cooperation is maintained because despite relatively high costs, clusters of cooperators cannot be invaded by selfish individuals.
in contrast, cheaters can always invade in our system, and do so faster when costs are higher.
this shows that costs are qualitatively higher in our model.
cooperation is maintained despite, and due to the evolution of true cheaters .
furthermore, we show that the amount of public good produced by an individual increases in evolutionary times, whereas the evolutionary stability of the solutions in  <cit>  is left unexplored, and only the long population dynamical transient is analysed.
more in general, the importance of spatial self-organisation for understanding the population dynamics of cooperators and defectors has been highlighted both from a theoretical  <cit>  and from an experimental point of view  <cit> .
here we make a similar point, but with an evolutionary twist: in our case selfish individuals are not merely a burden to cooperators; instead, their emergence as a separate lineage is necessary for the evolution of high degrees of public good production because they generate the spatial conditions in which cooperators thrive and evolve.
allowing mutations to change public good production in a continuous range resulted in the evolution of two separate strains, a selfish and a cooperative one.
the evolution of stable heterogeneity in a population has been observed before in models of cooperation  <cit> .
here, besides stressing that the evolution of two lineages from a single ancestral one might be a rather general feature of models with variable investments , we make the case that true cheating behaviour  is actually functional and beneficial to the long term evolution of cooperation.
CONCLUSIONS
in conclusion, besides extending the theoretical limits of cooperation, our results broaden the search image of cooperative behaviour in nature by suggesting that there need not be a strict trade-off between costs and benefits; rather, a wider view of the self-organised eco-evolutionary processes must be taken into account to understand the occurrence of costly cooperation.

 when cooperation is costly to the individual but its benefits are equally shared in a group, one would expect progressively more selfish behaviour to be selected. this indeed happens when interactions between cooperative and selfish individuals are random. instead, cooperation can be selected in a population when the interactions among cooperative and selfish individuals are structured, be it genetically, spatially or socially. population structure favours cooperation when it allows for cooperators to be in contact with each other more frequently than with selfish individuals. the result is a stable equilibrium configuration in which cooperators persist indefinitely, and selfish individuals may co-exist. for instance, this can be due to spatial clustering of cooperators, or to an inherent structure of the interaction network. then, the conditions for selfish individuals invading and overtaking a group of cooperators represent the limit to the stability of the solutions found. in general, these conditions state that higher costs of cooperation increase the selective advantage of selfish individuals. while spatial structure alone can favour cooperators due to population viscosity, a growing body of experimental and theoretical work indicates that self-organised spatial patterns may have profound and complex effects on cooperative interactions, due to emergent heterogeneities in the local distributions and densities of cooperators and selfish individuals. for instance, as a population of cooperators invade empty space, its expansion front can be enriched in altruistic individuals, while selfish individuals lag behind. alternatively, mutualistic interactions that are favoured in a resident population can be automatically broken on such expansion front. furthermore, spatial self-organisation can sort cooperative strains from selfish ones, thus limiting the spread of the latter. the customary approach to study the stability of cooperation consists of fixing the strategies of the interacting cooperative and selfish individuals, and analyse the population dynamics of the system. because such pre-determined strategies do not mutate over time, their evolutionary stability remains unexplored. although exceptions to this approach exist, little is known about what strategies evolve and how they feedback on the spatial self-organisation of a population, even though it is clear that spatial self-organisation affects the population dynamics of cooperative traits. here, we seek to study this feed-back between evolution and self-organisation with a minimal model where individuals can evolve the degree of cooperation in a spatially extended system. we model the cooperative trait in terms of public good production, and we let the amount of public good produced mutate in a continuous fashion. thus, we can study the long term evolutionary dynamics of cooperation without preconceiving the extent to which individuals cooperate or defect. as we will show, selfish individuals that produce zero public good evolve readily at higher costs, and quickly invade cooperators. rather than leading to global extinction, this enables cooperators to thrive and selects for a higher degree of public good production over evolutionary time scales.  the model is a straightforward implementation of a population in which individuals replicate depending on the amount of public good produced in their close neighbourhood. individuals are embedded on a lattice. they may reproduce, die or move. competition for reproduction into neighbouring empty nodes is based on fitness, calculated as the difference between benefits and costs. an individual benefits from the public good produced in its neighbourhood, but pays a cost for producing it. thus, public good production is a cooperative trait. we assume that reproductive success is solely based on cooperation, so that individuals do not reproduce if public good in their neighbourhood is insufficient. mutations slightly change the offspring’s production rate. fig. the model. a the world is a square lattice with connectivity k= and wrapped boundary conditions. individuals produce public good at rate p, shared in equal parts among all neighbouring nodes and self. b individuals compete for reproduction into an adjacent empty spot. probability of reproduction depends on fitness f i, which is the difference between benefits and costs. the sum of the public good an individual collects from itself and the neighbours, if any, is p sum=/ and confers a benefit b·p sum. individuals pay a fitness cost proportional to the public good they produce: c·p. successful reproduction yields a copy of the selected individual. mutations occur with probability μ and change the public good production rate by a small random number chosen uniformly in the interval. individuals have a small probability k move to move to a random adjacent node, and can die with probability k death, leaving the node empty high cost leads to the evolution of larger public good production we set the benefits per unit of public good b= and we let the spatial self-organisation and the evolutionary dynamics unfold under different costs c. when costs are much smaller than benefits, the public good production steadily increases because an individuals’ own production increases its fitness, rather than decreasing it. moreover, because replication is a local process, mutants with higher than average production rates benefit from each other due to limited dispersal, outcompeting more selfish lineages. thus, cooperation is maximised in the long run. for increasing costs, public good production suddenly drops. in this regime, the clustering advantage of cooperators is insufficient and more selfish individuals replicate the most because, by producing less, they pay a lower cost. eventually, public good production stabilises at the minimum value needed for survival. fig. an increase in costs results in an evolutionary increase in cooperation. a evolutionary steady state distribution and mean of public good production are plotted for different values of costs. parameters: benefit b= k death=, k move=, μ=, δ=. the maximum public good production is set to p max= b snapshots of the lattice at evolutionary steady state. colour coding depends on public good production rate. white is background. lattice dimensions used for the simulations from left to right: strikingly, further increasing costs leads to an increase in cooperation. the distribution of public good production is bimodal at evolutionary steady state, with most of the population having higher rates of public good production and a minority producing almost no public good at all. the long term evolutionary dynamics of cooperation at high cost following the evolutionary and the spatial dynamics of a single case elucidates why cooperation increases and persists for higher costs. when we initialise the system with highly producing individuals, public good production decreases rapidly due to strong selection for selfishness. where public good production drops below the minimum for survival, large patches of individuals go extinct. the surviving individuals can expand into the empty space. as the expansion progresses, a selfish and cooperative lineage separate from each other. the selfish strain evolves to zero public good production, becoming incapable of autonomous persistence and relying on the public good produced by cooperators for survival. fig. evolutionary dynamics of public good production. a at each time point the distribution of public good production in the lattice is plotted as a heat map. b-g snapshots of the lattice at subsequent time steps. costs c=, other parameters and colour coding in the snapshot as in fig. lattice size =. time units are monte carlo steps. see also additional file 1: section for the full snapshots, and movie at while the two strains differentiate from each other, they organise spatially to form travelling waves. cooperators constitute the front of a wave, and expand into empty space by replicating into it; selfish individuals invade those cooperators, and constitute the back of the wave. selfish individuals leave empty space behind a wave after they die, causing the semblance of movement. the progression of a wave, however, happens on a time scale that is much longer than the life time of an individual, which in turn experiences a fairly constant environment throughout its life. because waves persist longer than individuals, they can integrate information over several generations. fig. the movement of a travelling wave. magnification of the same portion of the lattice at time steps distance. a the wave-front, composed of cooperative individuals, b the back, composed of selfish individuals. the snapshots are from the same simulation run as in fig. at time steps and. the dashed line in the right pane marks approximately the position of the wave in the left pane. colizzi and hogeweg shows one such travelling wave in a video new waves are “born” from the collisions of older waves. as cooperators on the front evolve to larger public good production, waves become larger. the formation of spatial patterns allows the populations in the system to persist indefinitely despite selfish individuals continuously invading cooperators, provided that the lattice is much larger than the spatial patterns. spatial population dynamics of cooperators and selfish individuals spatial self-organisation drives the evolution of cooperation in the system. when spatial patterns are destroyed, e.g. by mixing, only selfish individuals are selected and public good production decreases, leading to global extinction, in accordance with the result that random interactions favour selfish behaviour. to unravel the interplay between the two lineages and their spatial organisation, we analysed the spatial population dynamics for cooperative and selfish individuals separately. to this end, we shaped the lattice into a long, narrow strip, and set the mutation rate to zero. cooperators expand faster into empty space when they produce more public good, and slower when costs are higher. when two clustered populations compete at the expansion front, the one with the largest public good production wins because, by replicating faster, it occupies space before the competing one and eventually overtakes the entire wave front. fig. invasion rates and spatial dynamics of competition for cooperative and selfish individuals. a cooperators’ invasion rate increases with larger public good production and lower costs; selfish invasion rate increases with increasing costs and it is insensitive to the production rates of the cooperators that support them. the invasion rates for cooperators invading empty space and selfish individuals invading a population of cooperators was measured for each combination of cost and public good production rate in b a population with lower public good production rate out-competes one with larger production in the back of a wave. parameters: cooperators p yellow=, selfish individuals p blue= and p magenta=. c the population with a larger production rate out-competes the one with lower production at the front of a wave. parameters: p purple=, p green=. for both a and b, μ= b= c=, background in white. other parameters as in fig. 2 the replication rate of selfish individuals invading a population of cooperators is higher when the cost of public good production is higher, and it is insensitive to how much public good is produced. clearly, when two strains compete in the back of a wave, the winner is the more selfish one. the picture emerging from these experiments is that different selection pressures operate depending on the spatial context: a population expanding into empty space is selected for higher degrees of cooperation, competition in the back selects for more selfishness. importantly, even though costs, benefits and fitness function are the same, spatial pattern formation automatically segregates these two opposing evolutionary pressures to spatially different contexts so that they do not balance each other: hence the evolution of a cooperative and a selfish lineage. in the full system, cooperators and selfish individuals are ecologically and evolutionarily interdependent. the two lineages establish an evolutionary feedback mediated by their spatial organisation. empty space is generated by selfish individuals after invading a population of cooperators. therefore, the condition for increasing public good production, i.e. the availability of continuous empty space, is mediated by the invasion dynamics of the selfish lineage. with higher costs, selfish individuals propagate faster, and space is left empty at a higher rate. the larger the empty space, the more cooperators can increase public good production. this evolutionary feedback reaches an evolutionary steady state because highly producing cooperators reduce the empty space faster upon faster expansions. we checked the long term stability of the steady state. altogether, cooperation evolves to a higher degree for higher costs due to an emergent feedback between self-organised interaction structures and the evolution of the individuals composing them. robustness to parameter change our results are robust when death rates, movement rates and benefits are changed, provided that benefits-to-costs ratio is maintained. at lower costs, individuals directly benefiting themselves with their own public good are sometimes dubbed weak altruists, with k+1= the connectivity of the lattice including self), whereas they are considered strong altruists when their public good is only shared among others. in our spatial model, we observed no qualitative difference in the evolutionary dynamics when individuals did or did not benefit from their own public good. this could be expected because individuals’ own payoffs at high costs are negative in both models. indeed, in both cases the evolving populations underwent speciation of a selfish and a cooperative lineage. weak or strong altruism do make a difference in the corresponding well-mixed systems, where strong altruists go extinct at lower costs, while weak altruists maximise public good production. discussion it is known that during population range expansion, cooperation can be promoted on the front of the expansion range. one could argue that in these models cooperation could evolve only as long as empty space is available, and should eventually be out-competed globally by selfish strategies when the invasion dynamics reach an end. here we have shown that selfish individuals provide the empty space to allow continuous expansion within a limited area. generally, stable solutions to the problem of cooperation are based on the condition that at equilibrium selfish individuals do not locally invade a population of cooperators. for instance, it is well known that cooperation can be stabilised in spatially extended systems, as cooperators cluster and segregate selfish individuals to the edges of those clusters. a side effect of these solutions is that higher costs undermines the stability of such clusters. in the parameter region where travelling waves do not form our work is in agreement with those results in that there exists an inverse relationship between costs and cooperation, and in particular with the heuristics that cooperation evolves when the benefit-to-cost ratio is larger than the connectivity of the lattice. however, we have shown that a novel class of solutions exists at high costs, where large degrees of cooperation are maintained in a locally out-of-equilibrium fashion, with selfish individuals always successfully invading cooperators and setting the stage for the evolutionary increase and the global stability of cooperation. we conclude that spatial self-organisation can reverse the relationship between costs and cooperation, thus extending the evolutionary viability of cooperation to higher costs. our results rest on two assumptions: population size can vary and some degree of cooperation is necessary for reproduction. variable population size is obviously realistic, even though it is not often included in evolutionary models of cooperation. although the assumption of necessary cooperation is not always met, it is reasonable in several cases. examples in microbiology include, cooperative protection or cooperative virulence in bacterial infections; invertase production in yeasts while growing on sucrose; siderophore production in iron-limited environments, cooperative secretion of digestive enzyme in microbial hunting. outside the microbial world, situations where our model may apply are e.g. dangerous behaviours in cooperative nest defence, and replication in trans in prebiotic evolution. two recent studies have come to conclusions that at first sight are similar to ours. by making the assumption, as we do, that the lack of cooperation leads to death, they observe static spatial patterns in which cooperation is maintained because despite relatively high costs, clusters of cooperators cannot be invaded by selfish individuals. in contrast, cheaters can always invade in our system, and do so faster when costs are higher. this shows that costs are qualitatively higher in our model. cooperation is maintained despite, and due to the evolution of true cheaters. furthermore, we show that the amount of public good produced by an individual increases in evolutionary times, whereas the evolutionary stability of the solutions in is left unexplored, and only the long population dynamical transient is analysed. more in general, the importance of spatial self-organisation for understanding the population dynamics of cooperators and defectors has been highlighted both from a theoretical and from an experimental point of view. here we make a similar point, but with an evolutionary twist: in our case selfish individuals are not merely a burden to cooperators; instead, their emergence as a separate lineage is necessary for the evolution of high degrees of public good production because they generate the spatial conditions in which cooperators thrive and evolve. allowing mutations to change public good production in a continuous range resulted in the evolution of two separate strains, a selfish and a cooperative one. the evolution of stable heterogeneity in a population has been observed before in models of cooperation. here, besides stressing that the evolution of two lineages from a single ancestral one might be a rather general feature of models with variable investments, we make the case that true cheating behaviour is actually functional and beneficial to the long term evolution of cooperation.  in conclusion, besides extending the theoretical limits of cooperation, our results broaden the search image of cooperative behaviour in nature by suggesting that there need not be a strict trade-off between costs and benefits; rather, a wider view of the self-organised eco-evolutionary processes must be taken into account to understand the occurrence of costly cooperation. 
an emergent feedback between evolution and self-organisation leads to high degrees of cooperation at high costs, under simple and general conditions. a general and intuitive result from theoretical models of cooperative behaviour is that cooperation decreases when its costs are higher, because selfish individuals gain selective advantage. we propose this as a general explanation for the evolution of cooperative behaviours under seemingly prohibitive conditions. we find that a feedback establishes between the evolutionary dynamics of public good production and the spatial self-organisation of the population. this enables selfish individuals to successfully invade cooperators at the expenses of their autonomous replication, and cooperators to increase public good production while expanding in the empty space left behind by cheaters. the ensuing spatial self-organisation automatically diversifies the selection pressure on the two lineages. the evolutionary dynamics lead to the speciation of a cooperative and a selfish lineage. we show that this emergent feedback leads to higher degrees of cooperation when costs are higher. we analyse a minimal model where individuals live on a lattice and evolve the degree of cooperation.
 cooperation is ubiquitous in biological systems, yet its evolution is a long lasting evolutionary problem. a general and intuitive result from theoretical models of cooperative behaviour is that cooperation decreases when its costs are higher, because selfish individuals gain selective advantage.  contrary to this intuition, we show that cooperation can increase with higher costs. we analyse a minimal model where individuals live on a lattice and evolve the degree of cooperation. we find that a feedback establishes between the evolutionary dynamics of public good production and the spatial self-organisation of the population. the evolutionary dynamics lead to the speciation of a cooperative and a selfish lineage. the ensuing spatial self-organisation automatically diversifies the selection pressure on the two lineages. this enables selfish individuals to successfully invade cooperators at the expenses of their autonomous replication, and cooperators to increase public good production while expanding in the empty space left behind by cheaters. we show that this emergent feedback leads to higher degrees of cooperation when costs are higher.  an emergent feedback between evolution and self-organisation leads to high degrees of cooperation at high costs, under simple and general conditions. we propose this as a general explanation for the evolution of cooperative behaviours under seemingly prohibitive conditions. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords evolution of cooperationspatial modelviscous populationstravelling wavesmulti-level evolutionnwo, complexity-net program.004colizzi enrico sandro issue-copyright-statement© the author 2016 