Inclusive fitness analysis on mathematical groups

The evolution of environmentally mediated social interactions and posthumous spite under isolation by distance

Many social interactions happen indirectly via modifications of the environment, e.g. through the secretion of functional compounds or the depletion of renewable resources. Here, we derive the selection gradient on a quantitative trait affecting dynamical environmental variables that feed back on reproduction and survival in a finite patch-structured population subject to isolation by distance. Our analysis shows that the selection gradient depends on how a focal individual influences the fitness of all future individuals in the population through modifications of the environmental variables they experience, weighted by the neutral relatedness between recipients and the focal. The evolutionarily relevant trait-driven environmental modifications are formalized as the extended phenotypic effects of an individual, quantifying how a trait change in an individual in the present affects the environmental variables in all patches at all future times. When the trait affects reproduction and survival through a payoff function, the selection gradient can be expressed in terms of extended phenotypic effects weighted by scaled relatedness. We show how to compute extended phenotypic effects, relatedness, and scaled relatedness using Fourier analysis, which allow us to investigate a broad class of environmentally mediated social interactions in a tractable way. We use our approach to study the evolution of a trait controlling the costly production of some lasting commons (e.g. a common-pool resource or a toxic compound) that can diffuse in space and persist in time. We show that indiscriminate posthumous spite readily evolves in this scenario. More generally, whether selection favours environmentally mediated altruism or spite is determined by the spatial correlation between an individual's lineage and the commons originating from its patch. The sign of this correlation depends on interactions between dispersal patterns and the commons' renewal dynamics. More broadly, we suggest that selection can favour a wide range of social behaviours when these have carry-over effects in space and time.

Symmetry in models of natural selection

Symmetry arguments are frequently used-often implicitly-in mathematical modelling of natural selection. Symmetry simplifies the analysis of models and reduces the number of distinct population states to be considered. Here, I introduce a formal definition of symmetry in mathematical models of natural selection. This definition applies to a broad class of models that satisfy a minimal set of assumptions, using a framework developed in previous works. In this framework, population structure is represented by a set of sites at which alleles can live, and transitions occur via replacement of some alleles by copies of others. A symmetry is defined as a permutation of sites that preserves probabilities of replacement and mutation. The symmetries of a given selection process form a group, which acts on population states in a way that preserves the Markov chain representing selection. Applying classical results on group actions, I formally characterize the use of symmetry to reduce the states of this Markov chain, and obtain bounds on the number of states in the reduced chain.

Multiple social encounters can eliminate Crozier's paradox and stabilise genetic kin recognition

Crozier's paradox suggests that genetic kin recognition will not be evolutionarily stable. The problem is that more common tags (markers) are more likely to be recognised and helped. This causes common tags to increase in frequency, and hence eliminates the genetic variability that is required for genetic kin recognition. It has therefore been assumed that genetic kin recognition can only be stable if there is some other factor maintaining tag diversity, such as the advantage of rare alleles in host-parasite interactions. We show that allowing for multiple social encounters before each social interaction can eliminate Crozier's paradox, because it allows individuals with rare tags to find others with the same tag. We also show that rare tags are better indicators of relatedness, and hence better at helping individuals avoid interactions with non-cooperative cheats. Consequently, genetic kin recognition provides an advantage to rare tags that maintains tag diversity, and stabilises itself.

Evolutionary games on isothermal graphs

Population structure affects the outcome of natural selection. These effects can be modeled using evolutionary games on graphs. Recently, conditions were derived for a trait to be favored under weak selection, on any weighted graph, in terms of coalescence times of random walks. Here we consider isothermal graphs, which have the same total edge weight at each node. The conditions for success on isothermal graphs take a simple form, in which the effects of graph structure are captured in the ‘effective degree’—a measure of the effective number of neighbors per individual. For two update rules (death-Birth and birth-Death), cooperative behavior is favored on a large isothermal graph if the benefit-to-cost ratio exceeds the effective degree. For two other update rules (Birth-death and Death-birth), cooperation is never favored. We relate the effective degree of a graph to its spectral gap, thereby linking evolutionary dynamics to the theory of expander graphs. Surprisingly, we find graphs of infinite average degree that nonetheless provide strong support for cooperation.

The spatial structure of a population is often critical for the evolution of cooperation. Here, Allen and colleagues show that when spatial structure is represented by an isothermal graph, the effective number of neighbors per individual determines whether or not cooperation can evolve.

Fidelity of parent-offspring transmission and the evolution of social behavior in structured populations

The theoretical investigation of how spatial structure affects the evolution of social behavior has mostly been done under the assumption that parent-offspring strategy transmission is perfect, i.e., for genetically transmitted traits, that mutation is very weak or absent. Here, we investigate the evolution of social behavior in structured populations under arbitrary mutation probabilities. We consider populations of fixed size N, structured such that in the absence of selection, all individuals have the same probability of reproducing or dying (neutral reproductive values are the all same). Two types of individuals, A and B, corresponding to two types of social behavior, are competing; the fidelity of strategy transmission from parent to offspring is tuned by a parameter μ. Social interactions have a direct effect on individual fecundities. Under the assumption of small phenotypic differences (implying weak selection), we provide a formula for the expected frequency of type A individuals in the population, and deduce conditions for the long-term success of one strategy against another. We then illustrate our results with three common life-cycles (Wright-Fisher, Moran Birth-Death and Moran Death-Birth), and specific population structures (graph-structured populations). Qualitatively, we find that some life-cycles (Moran Birth-Death, Wright-Fisher) prevent the evolution of altruistic behavior, confirming previous results obtained with perfect strategy transmission. We also show that computing the expected frequency of altruists on a regular graph may require knowing more than just the graph's size and degree.

How Life History Can Sway the Fixation Probability of Mutants

In this work, we study the effects of demographic structure on evolutionary dynamics when selection acts on reproduction, survival, or both. In contrast to the previously discovered pattern that the fixation probability of a neutral mutant decreases while the population becomes younger, we show that a mutant with a constant selective advantage may have a maximum or a minimum of the fixation probability in populations with an intermediate fraction of young individuals. This highlights the importance of life history and demographic structure in studying evolutionary dynamics. We also illustrate the fundamental differences between selection on reproduction and selection on survival when age structure is present. In addition, we evaluate the relative importance of size and structure of the population in determining the fixation probability of the mutant. Our work lays the foundation for also studying density- and frequency-dependent effects in populations when demographic structures cannot be neglected.

The causal meaning of Hamilton's rule

Hamilton's original derivation of his rule for the spread of an altruistic gene (rb>c) assumed additivity of costs and benefits. Recently, it has been argued that an exact version of the rule holds under non-additive pay-offs, so long as the cost and benefit terms are suitably defined, as partial regression coefficients. However, critics have questioned both the biological significance and the causal meaning of the resulting rule. This paper examines the causal meaning of the generalized Hamilton's rule in a simple model, by computing the effect of a hypothetical experiment to assess the cost of a social action and comparing it to the partial regression definition. The two do not agree. A possible way of salvaging the causal meaning of Hamilton's rule is explored, by appeal to R. A. Fisher's 'average effect of a gene substitution'.

Social evolution and genetic interactions in the short and long term

The evolution of social traits remains one of the most fascinating and feisty topics in evolutionary biology even after half a century of theoretical research. W.D. Hamilton shaped much of the field initially with his 1964 papers that laid out the foundation for understanding the effect of genetic relatedness on the evolution of social behavior. Early theoretical investigations revealed two critical assumptions required for Hamilton's rule to hold in dynamical models: weak selection and additive genetic interactions. However, only recently have analytical approaches from population genetics and evolutionary game theory developed sufficiently so that social evolution can be studied under the joint action of selection, mutation, and genetic drift. We review how these approaches suggest two timescales for evolution under weak mutation: (i) a short-term timescale where evolution occurs between a finite set of alleles, and (ii) a long-term timescale where a continuum of alleles are possible and populations evolve continuously from one monomorphic trait to another. We show how Hamilton's rule emerges from the short-term analysis under additivity and how non-additive genetic interactions can be accounted for more generally. This short-term approach reproduces, synthesizes, and generalizes many previous results including the one-third law from evolutionary game theory and risk dominance from economic game theory. Using the long-term approach, we illustrate how trait evolution can be described with a diffusion equation that is a stochastic analogue of the canonical equation of adaptive dynamics. Peaks in the stationary distribution of the diffusion capture classic notions of convergence stability from evolutionary game theory and generally depend on the additive genetic interactions inherent in Hamilton's rule. Surprisingly, the peaks of the long-term stationary distribution can predict the effects of simple kinds of non-additive interactions. Additionally, the peaks capture both weak and strong effects of social payoffs in a manner difficult to replicate with the short-term approach. Together, the results from the short and long-term approaches suggest both how Hamilton's insight may be robust in unexpected ways and how current analytical approaches can expand our understanding of social evolution far beyond Hamilton's original work.

The molecular clock of neutral evolution can be accelerated or slowed by asymmetric spatial structure

Over time, a population acquires neutral genetic substitutions as a consequence of random drift. A famous result in population genetics asserts that the rate, K, at which these substitutions accumulate in the population coincides with the mutation rate, u, at which they arise in individuals: K = u. This identity enables genetic sequence data to be used as a "molecular clock" to estimate the timing of evolutionary events. While the molecular clock is known to be perturbed by selection, it is thought that K = u holds very generally for neutral evolution. Here we show that asymmetric spatial population structure can alter the molecular clock rate for neutral mutations, leading to either Ku. Our results apply to a general class of haploid, asexually reproducing, spatially structured populations. Deviations from K = u occur because mutations arise unequally at different sites and have different probabilities of fixation depending on where they arise. If birth rates are uniform across sites, then K ≤ u. In general, K can take any value between 0 and Nu. Our model can be applied to a variety of population structures. In one example, we investigate the accumulation of genetic mutations in the small intestine. In another application, we analyze over 900 Twitter networks to study the effect of network topology on the fixation of neutral innovations in social evolution.

Collapse

Environmental evolutionary graph theory

Understanding the influence of an environment on the evolution of its resident population is a major challenge in evolutionary biology. Great progress has been made in homogeneous population structures while heterogeneous structures have received relatively less attention. Here we present a structured population model where different individuals are best suited to different regions of their environment. The underlying structure is a graph: individuals occupy vertices, which are connected by edges. If an individual is suited for their vertex, they receive an increase in fecundity. This framework allows attention to be restricted to the spatial arrangement of suitable habitat. We prove some basic properties of this model and find some counter-intuitive results. Notably, (1) the arrangement of suitable sites is as important as their proportion, and (2) decreasing the proportion of suitable sites may result in a decrease in the fixation time of an allele.

A mathematical description of the inclusive fitness theory

Recent developments in the inclusive fitness theory have revealed that the direction of evolution can be analytically predicted in a wider class of models than previously thought, such as those models dealing with network structure. This paper aims to provide a mathematical description of the inclusive fitness theory. Specifically, we provide a general framework based on a Markov chain that can implement basic models of inclusive fitness. Our framework is based on the probability distribution of "offspring-to-parent map", from which the key concepts of the theory, such as fitness function, relatedness and inclusive fitness, are derived in a straightforward manner. We prove theorems showing that inclusive fitness always provides a correct prediction on which of two competing genes more frequently appears in the long run in the Markov chain. As an application of the theorems, we prove a general formula of the optimal dispersal rate in the Wright's island model with recurrent mutations. We also show the existence of the critical mutation rate, which does not depend on the number of islands and below which a positive dispersal rate evolves. Our framework can also be applied to lattice or network structured populations.

Measures of success in a class of evolutionary models with fixed population size and structure

We investigate a class of evolutionary models, encompassing many established models of well-mixed and spatially structured populations. Models in this class have fixed population size and structure. Evolution proceeds as a Markov chain, with birth and death probabilities dependent on the current population state. Starting from basic assumptions, we show how the asymptotic (long-term) behavior of the evolutionary process can be characterized by probability distributions over the set of possible states. We then define and compare three quantities characterizing evolutionary success: fixation probability, expected frequency, and expected change due to selection. We show that these quantities yield the same conditions for success in the limit of low mutation rate, but may disagree when mutation is present. As part of our analysis, we derive versions of the Price equation and the replicator equation that describe the asymptotic behavior of the entire evolutionary process, rather than the change from a single state. We illustrate our results using the frequency-dependent Moran process and the birth-death process on graphs as examples. Our broader aim is to spearhead a new approach to evolutionary theory, in which general principles of evolution are proven as mathematical theorems from axioms.

An inclusive fitness analysis of synergistic interactions in structured populations

We study the evolution of a pair of competing behavioural alleles in a structured population when there are non-additive or 'synergistic' fitness effects. Under a form of weak selection and with a simple symmetry condition between a pair of competing alleles, Tarnita et al. provide a surprisingly simple condition for one allele to dominate the other. Their condition can be obtained from an analysis of a corresponding simpler model in which fitness effects are additive. Their result uses an average measure of selective advantage where the average is taken over the long-term--that is, over all possible allele frequencies--and this precludes consideration of any frequency dependence the allelic fitness might exhibit. However, in a considerable body of work with non-additive fitness effects--for example, hawk-dove and prisoner's dilemma games--frequency dependence plays an essential role in the establishment of conditions for a stable allele-frequency equilibrium. Here, we present a frequency-dependent generalization of their result that provides an expression for allelic fitness at any given allele frequency p. We use an inclusive fitness approach and provide two examples for an infinite structured population. We illustrate our results with an analysis of the hawk-dove game.

Does synergy rescue the evolution of cooperation? An analysis for homogeneous populations with non-overlapping generations

Recent developments of social evolution theory have revealed conditions under which cooperation is favored by natural selection. Effects of population structure on the evolution of cooperation have been one of the central questions in this issue, and inclusive fitness analyses have unveiled two different selective forces that favor cooperation; the direct fitness effect to the helper and the indirect fitness benefit to the helper via its kin. Although these theoretical frameworks have made a significant contribution to our understanding of cooperative traits, there is still one factor to be taken into account, synergy. Synergy means a nonlinear effect that arises when two individuals help each other. In other words, it represents deviation from additivity, to which inclusive fitness theory has paid relatively little attention. Here I provide a theoretical result on the possibility that synergy favors the evolution of cooperation. For homogeneously structured populations with non-overlapping generations, I show that incorporating synergistic effects does not rescue the evolution of cooperation. Potential factors that could enable synergy to rescue the evolution of cooperation are also discussed.

Resistance and relatedness on an evolutionary graph

When investigating evolution in structured populations, it is often convenient to consider the population as an evolutionary graph-individuals as nodes, and whom they may act with as edges. There has, in recent years, been a surge of interest in evolutionary graphs, especially in the study of the evolution of social behaviours. An inclusive fitness framework is best suited for this type of study. A central requirement for an inclusive fitness analysis is an expression for the genetic similarity between individuals residing on the graph. This has been a major hindrance for work in this area as highly technical mathematics are often required. Here, I derive a result that links genetic relatedness between haploid individuals on an evolutionary graph to the resistance between vertices on a corresponding electrical network. An example that demonstrates the potential computational advantage of this result over contemporary approaches is provided. This result offers more, however, to the study of population genetics than strictly computationally efficient methods. By establishing a link between gene transfer and electric circuit theory, conceptualizations of the latter can enhance understanding of the former.