Extinction times in diffusive public good population dynamics

The co-evolutionary dynamics of competing populations can be strongly affected by frequency-dependent selection and population structure in space. As co-evolving populations grow into a spatial domain, their initial spatial arrangement, as well as their growth rate differences determine the dynamics. Here, we are interested in the dynamics of producers and free-rider co-evolution in the context of an ecological public good that is produced by a sub-population but evokes growth benefits to all individuals. We consider the spatial growth dynamics in one, two and three dimensions by modeling producer cell, free-rider cell and public good densities in space, driven by birth, death and diffusion. Typically, one population goes extinct. We find that uncorrelated initial spatial structures do not influence the time to extinction in comparison to the well-mixed system. We derive a slow manifold solution in order to estimate the time to extinction of either free-riders or producers. For invading populations, i.e. for populations that are initially highly segregated, we observe a traveling wave, whose speed can be calculated to improve the extinction time estimate by a simple superposition of the two times. Our results show that local effects of spatial dynamics evolve independently of the dynamics of the mean populations. Our considerations provide quantitative predictions for the transient dynamics of cooperative traits under pressure of extinction, and a potential experiment to derive elusive details of the fitness function of an ecological public goods game through extinction time observations. Author Summary Ecological public goods (PG) relationships emerge in growing cellular populations, for example between bacteria and cancer cells. We study the eco-evolutionary dynamics of a PG in populations that grow in space. In our model, public good-producer cells and free-rider cells can grow according to their own birth and death rates. Co-evolution occurs due to public good-driven surplus in the intrinsic growth rates and a cost to producers. A net growth rate benefit to free-riders leads to the well-known tragedy of the commons in which producers go extinct. What is often omitted from discussions is the time scale on which this extinction can occur, especially in spatial populations. We derive analytical estimates of the time to extinction in different spatial settings, and identify spatial scenarios in which extinction takes long enough such that the tragedy of the commons never occurs within the lifetime of the populations. Using numerical simulations we analyze the deviations from analytical predictions. Our results have direct implications for inferring ecological public good game properties from in vitro and in vivo experimental observations.

[1]  R. Fisher THE WAVE OF ADVANCE OF ADVANTAGEOUS GENES , 1937 .

[2]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947 .

[3]  N. Goel,et al.  Stochastic models in biology , 1975 .

[4]  P. Taylor,et al.  Evolutionarily Stable Strategies and Game Dynamics , 1978 .

[5]  W. Karnes,et al.  Autocrine growth factors and solid tumor malignancy. , 1991, The Western journal of medicine.

[6]  M. Nowak,et al.  Evolutionary games and spatial chaos , 1992, Nature.

[7]  R. Durrett,et al.  The Importance of Being Discrete (and Spatial) , 1994 .

[8]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947, Mathematical Proceedings of the Cambridge Philosophical Society.

[9]  M. Chaplain,et al.  Continuous and discrete mathematical models of tumor-induced angiogenesis , 1998, Bulletin of mathematical biology.

[10]  A. Maroudas,et al.  Synthesis of insulin-like growth factor binding protein 3 in vitro in human articular cartilage cultures. , 2003, Arthritis and rheumatism.

[11]  W. Ewens Mathematical Population Genetics : I. Theoretical Introduction , 2004 .

[12]  M. Nowak,et al.  A symmetry of fixation times in evoultionary dynamics. , 2006, Journal of theoretical biology.

[13]  C. Hauert,et al.  Evolutionary games and population dynamics: maintenance of cooperation in public goods games , 2006, Proceedings of the Royal Society B: Biological Sciences.

[14]  F. Lanni,et al.  Diffusion of insulin-like growth factor-I and ribonuclease through fibrin gels. , 2007, Biophysical journal.

[15]  K. Foster,et al.  Cooperation and conflict in microbial biofilms , 2007, Proceedings of the National Academy of Sciences.

[16]  A. Traulsen,et al.  Fixation times in evolutionary games under weak selection , 2008, 0812.0851.

[17]  C. Plott,et al.  Handbook of Experimental Economics Results , 2008 .

[18]  C. Hauert,et al.  Spatial dynamics of ecological public goods , 2009, Proceedings of the National Academy of Sciences.

[19]  T. Vincent,et al.  Lessons from applied ecology: cancer control using an evolutionary double bind. , 2009, Cancer research.

[20]  B. Frieden,et al.  Adaptive therapy. , 2009, Cancer research.

[21]  R. Gatenby A change of strategy in the war on cancer , 2009, Nature.

[22]  Carey D. Nadell,et al.  Emergence of Spatial Structure in Cell Groups and the Evolution of Cooperation , 2010, PLoS Comput. Biol..

[23]  C. Hauert,et al.  Pattern formation and chaos in spatial ecological public goods games , 2010 .

[24]  Chaitanya S. Gokhale,et al.  Stochastic slowdown in evolutionary processes. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[25]  C. Hauert,et al.  Pattern formation and chaos in spatial ecological public goodsgames. , 2011, Journal of theoretical biology.

[26]  M. Silber,et al.  Transitions between patterned states in vegetation models for semiarid ecosystems. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  M. Archetti Evolutionary dynamics of the Warburg effect: glycolysis as a collective action problem among cancer cells. , 2014, Journal of theoretical biology.

[28]  M. Archetti,et al.  Heterogeneity for IGF-II production maintained by public goods dynamics in neuroendocrine pancreatic cancer , 2015, Proceedings of the National Academy of Sciences.

[29]  K. Polyak,et al.  Tumorigenesis: it takes a village , 2015, Nature Reviews Cancer.

[30]  Complexity and stability in growing cancer cell populations , 2015, Proceedings of the National Academy of Sciences.

[31]  R. Gillies,et al.  Exploiting evolutionary principles to prolong tumor control in preclinical models of breast cancer , 2016, Science Translational Medicine.

[32]  K. Korolev,et al.  Range expansions transition from pulled to pushed waves as growth becomes more cooperative in an experimental microbial population , 2016, Proceedings of the National Academy of Sciences.

[33]  J. Gore,et al.  Self-organized patchiness facilitates survival in a cooperatively growing Bacillus subtilis population , 2016, Nature Microbiology.

[34]  Franziska Wulf,et al.  Mathematical Population Genetics , 2016 .

[35]  A. Bayliss,et al.  Complex predator invasion waves in a Holling–Tanner model with nonlocal prey interaction , 2017 .

[36]  Robert A. Gatenby,et al.  Aggregation Effects and Population-Based Dynamics as a Source of Therapy Resistance in Cancer , 2017, IEEE Transactions on Biomedical Engineering.

[37]  P. Altrock,et al.  Extinction rates in tumour public goods games , 2017, Journal of The Royal Society Interface.

[38]  M. Nowak,et al.  Evolutionary games on cycles with strong selection. , 2017, Physical review. E.

[39]  Joel s. Brown,et al.  Aggregation Effects and Population-Based Dynamics as a Source of Therapy Resistance in Cancer. , 2017, IEEE transactions on bio-medical engineering.

[40]  Joel s. Brown,et al.  Integrating evolutionary dynamics into treatment of metastatic castrate-resistant prostate cancer , 2017, Nature Communications.

[41]  Ying-Cheng Lai,et al.  Transient phenomena in ecology , 2018, Science.

[42]  Gregory J. Kimmel,et al.  Neighborhood size-effects in nonlinear public goods games , 2018, bioRxiv.