Inclusive Fitness from Multitype Branching Processes

The evolution of manipulative cheating

A social cheat is typically assumed to be an individual that does not perform a cooperative behaviour, or performs less of it, but can still exploit the cooperative behaviour of others. However, empirical data suggests that cheating can be more subtle, involving evolutionary arms races over the ability to both exploit and resist exploitation. These complications have not been captured by evolutionary theory, which lags behind empirical studies in this area. We bridge this gap with a mixture of game-theoretical models and individual-based simulations, examining what conditions favour more elaborate patterns of cheating. We found that as well as adjusting their own behaviour, individuals can be selected to manipulate the behaviour of others, which we term 'manipulative cheating'. Further, we found that manipulative cheating can lead to dynamic oscillations (arms races), between selfishness, manipulation, and suppression of manipulation. Our results can help explain both variation in the level of cheating, and genetic variation in the extent to which individuals can be exploited by cheats.

Metacommunities, fitness and gradual evolution

We analyze the evolution of a multidimensional quantitative trait in a class-structured focal species interacting with other species in a wider metacommunity. The evolutionary dynamics in the focal species as well as the ecological dynamics of the whole metacommunity is described as a continuous-time process with birth, physiological development, dispersal, and death given as rates that can depend on the state of the whole metacommunity. This can accommodate complex local community and global metacommunity environmental feedbacks owing to inter- and intra-specific interactions, as well as local environmental stochastic fluctuations. For the focal species, we derive a fitness measure for a mutant allele affecting class-specific trait expression. Using classical results from geometric singular perturbation theory, we provide a detailed proof that if the effect of the mutation on phenotypic expression is small ("weak selection"), the large system of dynamical equations needed to describe selection on the mutant allele in the metacommunity can be reduced to a single ordinary differential equation on the arithmetic mean mutant allele frequency that is of constant sign. This invariance on allele frequency entails the mutant either dies out or will out-compete the ancestral resident (or wild) type. Moreover, the directional selection coefficient driving arithmetic mean allele frequency can be expressed as an inclusive fitness effect calculated from the resident metacommunity alone, and depends, as expected, on individual fitness differentials, relatedness, and reproductive values. This formalizes the Darwinian process of gradual evolution driven by random mutation and natural selection in spatially and physiologically class-structured metacommunities.

The components of directional and disruptive selection in heterogeneous group-structured populations

We derive how directional and disruptive selection operate on scalar traits in a heterogeneous group-structured population for a general class of models. In particular, we assume that each group in the population can be in one of a finite number of states, where states can affect group size and/or other environmental variables, at a given time. Using up to second-order perturbation expansions of the invasion fitness of a mutant allele, we derive expressions for the directional and disruptive selection coefficients, which are sufficient to classify the singular strategies of adaptive dynamics. These expressions include first- and second-order perturbations of individual fitness (expected number of settled offspring produced by an individual, possibly including self through survival); the first-order perturbation of the stationary distribution of mutants (derived here explicitly for the first time); the first-order perturbation of pairwise relatedness; and reproductive values, pairwise and three-way relatedness, and stationary distribution of mutants, each evaluated under neutrality. We introduce the concept of individual k-fitness (defined as the expected number of settled offspring of an individual for which k-1 randomly chosen neighbors are lineage members) and show its usefulness for calculating relatedness and its perturbation. We then demonstrate that the directional and disruptive selection coefficients can be expressed in terms individual k-fitnesses with k=1,2,3 only. This representation has two important benefits. First, it allows for a significant reduction in the dimensions of the system of equations describing the mutant dynamics that needs to be solved to evaluate explicitly the two selection coefficients. Second, it leads to a biologically meaningful interpretation of their components. As an application of our methodology, we analyze directional and disruptive selection in a lottery model with either hard or soft selection and show that many previous results about selection in group-structured populations can be reproduced as special cases of our model.

The evolution of eusociality: no risk-return tradeoff but the ecology matters

The origin of eusociality in the Hymenoptera is a question of major interest. Theory has tended to focus on genetic relatedness, but ecology can be just as important a determinant of whether eusociality evolves. Using the model of Fu et al. (2015), we show how ecological assumptions critically affect the conclusions drawn. Fu et al. inferred that eusociality rarely evolves because it faces a fundamental 'risk-return tradeoff'. The intuitive logic was that worker production represents an opportunity cost because it delays realising a reproductive payoff. However, making empirically justified assumptions that (1) workers take over egg-laying following queen death and (2) productivity increases gradually with each additional worker, we find that the risk-return tradeoff disappears. We then survey Hymenoptera with more specialised morphological castes, and show how the interaction between two common features of eusociality - saturating birth rates and group size-dependent helping decisions - can determine whether eusociality outperforms other strategies.

An evolutionary quantitative genetics model for phenotypic (co)variances under limited dispersal, with an application to socially synergistic traits

Darwinian evolution consists of the gradual transformation of heritable traits due to natural selection and the input of random variation by mutation. Here, we use a quantitative genetics approach to investigate the coevolution of multiple quantitative traits under selection, mutation, and limited dispersal. We track the dynamics of trait means and of variance-covariances between traits that experience frequency-dependent selection. Assuming a multivariate-normal trait distribution, we recover classical dynamics of quantitative genetics, as well as stability and evolutionary branching conditions of invasion analyses, except that due to limited dispersal, selection depends on indirect fitness effects and relatedness. In particular, correlational selection that associates different traits within-individuals depends on the fitness effects of such associations between-individuals. We find that these kin selection effects can be as relevant as pleiotropy for the evolution of correlation between traits. We illustrate this with an example of the coevolution of two social traits whose association within-individuals is costly but synergistically beneficial between-individuals. As dispersal becomes limited and relatedness increases, associations between-traits between-individuals become increasingly targeted by correlational selection. Consequently, the trait distribution goes from being bimodal with a negative correlation under panmixia to unimodal with a positive correlation under limited dispersal.

Environmental fitness heterogeneity in the Moran process

Many mathematical models of evolution assume that all individuals experience the same environment. Here, we study the Moran process in heterogeneous environments. The population is of finite size with two competing types, which are exposed to a fixed number of environmental conditions. Reproductive rate is determined by both the type and the environment. We first calculate the condition for selection to favour the mutant relative to the resident wild-type. In large populations, the mutant is favoured if and only if the mutant's spatial average reproductive rate exceeds that of the resident. But environmental heterogeneity elucidates an interesting asymmetry between the mutant and the resident. Specifically, mutant heterogeneity suppresses its fixation probability; if this heterogeneity is strong enough, it can even completely offset the effects of selection (including in large populations). By contrast, resident heterogeneity has no effect on a mutant's fixation probability in large populations and can amplify it in small populations.

Public goods games in populations with fluctuating size

Many mathematical frameworks of evolutionary game dynamics assume that the total population size is constant and that selection affects only the relative frequency of strategies. Here, we consider evolutionary game dynamics in an extended Wright-Fisher process with variable population size. In such a scenario, it is possible that the entire population becomes extinct. Survival of the population may depend on which strategy prevails in the game dynamics. Studying cooperative dilemmas, it is a natural feature of such a model that cooperators enable survival, while defectors drive extinction. Although defectors are favored for any mixed population, random drift could lead to their elimination and the resulting pure-cooperator population could survive. On the other hand, if the defectors remain, then the population will quickly go extinct because the frequency of cooperators steadily declines and defectors alone cannot survive. In a mutation-selection model, we find that (i) a steady supply of cooperators can enable long-term population survival, provided selection is sufficiently strong, and (ii) selection can increase the abundance of cooperators but reduce their relative frequency. Thus, evolutionary game dynamics in populations with variable size generate a multifaceted notion of what constitutes a trait's long-term success.

Invasion fitness, inclusive fitness, and reproductive numbers in heterogeneous populations

How should fitness be measured to determine which phenotype or "strategy" is uninvadable when evolution occurs in a group-structured population subject to local demographic and environmental heterogeneity? Several fitness measures, such as basic reproductive number, lifetime dispersal success of a local lineage, or inclusive fitness have been proposed to address this question, but the relationships between them and their generality remains unclear. Here, we ascertain uninvadability (all mutant strategies always go extinct) in terms of the asymptotic per capita number of mutant copies produced by a mutant lineage arising as a single copy in a resident population ("invasion fitness"). We show that from invasion fitness uninvadability is equivalently characterized by at least three conceptually distinct fitness measures: (i) lineage fitness, giving the average individual fitness of a randomly sampled mutant lineage member; (ii) inclusive fitness, giving a reproductive value weighted average of the direct fitness costs and relatedness weighted indirect fitness benefits accruing to a randomly sampled mutant lineage member; and (iii) basic reproductive number (and variations thereof) giving lifetime success of a lineage in a single group, and which is an invasion fitness proxy. Our analysis connects approaches that have been deemed different, generalizes the exact version of inclusive fitness to class-structured populations, and provides a biological interpretation of natural selection on a mutant allele under arbitrary strength of selection.

Evolutionary Stability of Jointly Evolving Traits in Subdivided Populations

The evolutionary stability of quantitative traits depends on whether a population can resist invasion by any mutant. While uninvadability is well understood in well-mixed populations, it is much less so in subdivided populations when multiple traits evolve jointly. Here, we investigate whether a spatially subdivided population at a monomorphic equilibrium for multiple traits can withstand invasion by any mutant or is subject to diversifying selection. Our model also explores the correlations among traits arising from diversifying selection and how they depend on relatedness due to limited dispersal. We find that selection tends to favor a positive (negative) correlation between two traits when the selective effects of one trait on relatedness is positively (negatively) correlated to the indirect fitness effects of the other trait. We study the evolution of traits for which this matters: dispersal that decreases relatedness and helping that has positive indirect fitness effects. We find that when dispersal cost is low and the benefits of helping accelerate faster than its costs, selection leads to the coexistence of mobile defectors and sessile helpers. Otherwise, the population evolves to a monomorphic state with intermediate helping and dispersal. Overall, our results highlight the effects of population subdivision for evolutionary stability and correlations among traits.

Does evolution lead to maximizing behavior?

A long-standing question in biology and economics is whether individual organisms evolve to behave as if they were striving to maximize some goal function. We here formalize this "as if" question in a patch-structured population in which individuals obtain material payoffs from (perhaps very complex multimove) social interactions. These material payoffs determine personal fitness and, ultimately, invasion fitness. We ask whether individuals in uninvadable population states will appear to be maximizing conventional goal functions (with population-structure coefficients exogenous to the individual's behavior), when what is really being maximized is invasion fitness at the genetic level. We reach two broad conclusions. First, no simple and general individual-centered goal function emerges from the analysis. This stems from the fact that invasion fitness is a gene-centered multigenerational measure of evolutionary success. Second, when selection is weak, all multigenerational effects of selection can be summarized in a neutral type-distribution quantifying identity-by-descent between individuals within patches. Individuals then behave as if they were striving to maximize a weighted sum of material payoffs (own and others). At an uninvadable state it is as if individuals would freely choose their actions and play a Nash equilibrium of a game with a goal function that combines self-interest (own material payoff), group interest (group material payoff if everyone does the same), and local rivalry (material payoff differences).

Inclusive-fitness logic of cooperative breeding with benefits of natal philopatry

In cooperatively breeding species, individuals help to raise offspring that are not their own. We use two inclusive-fitness models to study the advantage of this kind of helpful behaviour in social groups with high reproductive skew. Our first model does not allow for competition among relatives to occur but our second model does. Specifically, our second model assumes a competitive hierarchy among nest-mates, with non-breeding helpers ranked higher than their newborn siblings. For each model, we obtain an expression for the change in inclusive fitness experienced by a helpful individual in a selfish population. The prediction suggested by each expression is confirmed with computer simulation. When model predictions are compared to one another, we find that helping emerges under a broader range of conditions in the second model. Although competition among kin occurs in our second model, we conclude that the life-history features associated with this competition also act to promote the evolutionary transition from solitary to cooperative breeding.