Constructing transient amplifiers for death-Birth updating: A case study of cubic and quartic regular graphs

A central question of evolutionary dynamics on graphs is whether or not a mutation introduced in a population of residents survives and eventually even spreads to the whole population, or gets extinct. The outcome naturally depends on the fitness of the mutant and the rules by which mutants and residents may propagate on the network, but arguably the most determining factor is the network structure. Some structured networks are transient amplifiers. They increase for a certain fitness range the fixation probability of beneficial mutations as compared to a well-mixed population. We study a perturbation methods for identifying transient amplifiers for death-Birth updating. The method includes calculating the coalescence times of random walks on graphs and finding the vertex with the largest remeeting time. If the graph is perturbed by removing an edge from this vertex, there is a certain likelihood that the resulting perturbed graph is a transient amplifier. We test all pairwise nonisomorphic cubic and quartic regular graphs up to a certain size and thus cover the whole structural range expressible by these graphs. We carry out a spectral analysis and show that the graphs from which the transient amplifiers can be constructed share certain structural properties. The graphs are path-like, have low conductance and are rather easy to divide into subgraphs by removing edges and/or vertices. This is connected with the subgraphs being identical (or almost identical) building blocks and the frequent occurrence of cut and/or hinge vertices. Identifying spectral and structural properties may promote finding and designing such networks.

[1]  Anirban Banerjee,et al.  Structural distance and evolutionary relationship of networks , 2008, Biosyst..

[2]  P. Stadler,et al.  Spectral classes of regular, random, and empirical graphs , 2014, 1406.6454.

[3]  A. Hartmann,et al.  Bias in generation of random graphs. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[4]  Hendrik Richter,et al.  Evolution of Cooperation for Multiple Mutant Configurations on All Regular Graphs with N≤14 players , 2019, Games.

[5]  F. Alcalde Cuesta,et al.  Evolutionary regime transitions in structured populations , 2018, PloS one.

[6]  Arne Traulsen,et al.  Exact numerical calculation of fixation probability and time on graphs , 2015, Biosyst..

[7]  Eric Lewitus,et al.  Characterizing and comparing phylogenies from their Laplacian spectrum , 2015, bioRxiv.

[8]  Anirban Banerjee,et al.  Graph spectra as a systematic tool in computational biology , 2007, Discret. Appl. Math..

[9]  A. Barabasi,et al.  Spectra of "real-world" graphs: beyond the semicircle law. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[10]  François G. Meyer,et al.  Metrics for graph comparison: A practitioner’s guide , 2019, bioRxiv.

[11]  Anirban Banerjee,et al.  Spectral plots and the representation and interpretation of biological data , 2007, Theory in Biosciences.

[12]  Krishnendu Chatterjee,et al.  Limits on amplifiers of natural selection under death-Birth updating , 2019, PLoS Comput. Biol..

[13]  Michael Doebeli,et al.  Spatial structure often inhibits the evolution of cooperation in the snowdrift game , 2004, Nature.

[14]  Mark Broom,et al.  Evolutionary dynamics on small-order graphs , 2009 .

[15]  Hendrik Richter,et al.  Dynamic landscape models of coevolutionary games , 2016, Biosyst..

[16]  Arne Traulsen,et al.  Exploring and mapping the universe of evolutionary graphs identifies structural properties affecting fixation probability and time , 2018, Communications Biology.

[17]  M. Nowak,et al.  Evolutionary games on isothermal graphs , 2019, Nature Communications.

[18]  M. Nowak,et al.  The spatial ultimatum game , 2000, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[19]  Yue-Li Wang,et al.  A Linear Time Algorithm for Finding all Hinge Vertices of a Permutation Graph , 1996, Inf. Process. Lett..

[20]  Martin A Nowak,et al.  Spatial invasion of cooperation. , 2008, Journal of theoretical biology.

[21]  Arne Traulsen,et al.  Most Undirected Random Graphs Are Amplifiers of Selection for Birth-Death Dynamics, but Suppressors of Selection for Death-Birth Dynamics , 2015, PLoS Comput. Biol..

[22]  B. McKay The expected eigenvalue distribution of a large regular graph , 1981 .

[23]  Krishnendu Chatterjee,et al.  Construction of arbitrarily strong amplifiers of natural selection using evolutionary graph theory , 2018, Communications Biology.

[24]  N. Wormald,et al.  Models of the , 2010 .

[25]  Anirban Banerjee,et al.  Effect on normalized graph Laplacian spectrum by motif attachment and duplication , 2012, Appl. Math. Comput..

[26]  Christoph Hauert,et al.  Fixation probabilities on superstars, revisited and revised. , 2013, Journal of theoretical biology.

[27]  Jacobo Torán,et al.  Isomorphism Testing: Perspective and Open Problems , 2005, Bull. EATCS.

[28]  Martin A. Nowak,et al.  Evolution of cooperation on large networks with community structure , 2019, Journal of the Royal Society Interface.

[29]  Markus Meringer,et al.  Fast generation of regular graphs and construction of cages , 1999, J. Graph Theory.

[30]  Martin A Nowak,et al.  Evolutionary graph theory: breaking the symmetry between interaction and replacement. , 2007, Journal of theoretical biology.

[31]  Karan Pattni,et al.  Evolutionary graph theory revisited: when is an evolutionary process equivalent to the Moran process? , 2015, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[32]  Chiun-Chieh Hsu,et al.  Finding the Set of All Hinge Vertices for Strongly Chordal Graphs in Linear Time , 1997, Inf. Sci..

[33]  Fernando Alcalde Cuesta,et al.  An Accurate Database of the Fixation Probabilities for All Undirected Graphs of Order 10 or Less , 2017, IWBBIO.

[34]  C. Hauert Fundamental clusters in spatial 2×2 games , 2001, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[35]  Benjamin Allen,et al.  Transient amplifiers of selection and reducers of fixation for death-Birth updating on graphs , 2019, PLoS Comput. Biol..

[36]  Hendrik Richter Fixation properties of multiple cooperator configurations on regular graphs , 2019, Theory in Biosciences.

[37]  M. Nowak,et al.  Supplementary Note 1: Population structure determines the tradeoff between fixation probability and fixation time , 2019 .

[38]  Krishnendu Chatterjee,et al.  Amplification on Undirected Population Structures: Comets Beat Stars , 2017, Scientific Reports.

[39]  László Babai,et al.  GROUP, GRAPHS, ALGORITHMS: THE GRAPH ISOMORPHISM PROBLEM , 2019, Proceedings of the International Congress of Mathematicians (ICM 2018).

[40]  Amin Saberi,et al.  A Sequential Algorithm for Generating Random Graphs , 2007, Algorithmica.

[41]  Travis Monk,et al.  Martingales and the fixation probability of high-dimensional evolutionary graphs. , 2018, Journal of theoretical biology.

[42]  Bin Wu,et al.  Computation and Simulation of Evolutionary Game Dynamics in Finite Populations , 2019, Scientific Reports.

[43]  Ping Zhu,et al.  A study of graph spectra for comparing graphs and trees , 2008, Pattern Recognit..

[44]  Hendrik Richter,et al.  Properties of network structures, structure coefficients, and benefit-to-cost ratios , 2019, Biosyst..

[45]  S. Shen-Orr,et al.  Networks Network Motifs : Simple Building Blocks of Complex , 2002 .

[46]  Matthew Kahle,et al.  Spectral Gaps of Random Graphs and Applications , 2012, International Mathematics Research Notices.

[47]  N. Wormald Models of random regular graphs , 2010 .

[48]  Uzy Smilansky,et al.  Trace formulae and spectral statistics for discrete Laplacians on regular graphs (I) , 2009, 1003.1445.