The Efficacy of Group Selection is Increased by Coexistence Dynamics within Groups
Selection on the level of loosely associated groups has been suggested as a route towards the evolution of cooperation between individuals and the subsequent formation of higher-level biological entities. Such group selection explanations remain problematic, however, due to the narrow range of parameters under which they can overturn within-group selection that favours selfish behaviour. In principle, individual selection could act on such parameters so as to strengthen the force of between-group selection and hence increase cooperation and individual fitness, as illustrated in our previous work. However, such a process cannot operate in parameter regions where group selection effects are totally absent, since there would be no selective gradient to follow. One key parameter, which when increased often rapidly causes group selection effects to tend to zero, is initial group size, for when groups are formed randomly then even moderately sized groups lack significant variance in their composition. However, the consequent restriction of any group selection effect to small sized groups is derived from models that assume selfish types will competitively exclude their more cooperative counterparts at within-group equilibrium. In such cases, diversity in the migrant pool can tend to zero and accordingly variance in group composition cannot be generated. In contrast, we show that if within-group dynamics lead to a stable coexistence of selfish and cooperative types, then the range of group sizes showing some effect of group selection is much larger.
The Efficacy of Group Selection is Increased by Coexistence Dynamics within Groups
Simon T. Powers, Alexandra S. Penn and Richard A. Watson School of Electronics & Computer Science, University of Southampton, Highfield, Southampton SO17 1BJ, UK Simon.Powers@Unil.ch
The evolution of cooperation between biological individuals, both generally and as a vital part of the formation of new higher-level composite individuals, is an important and much discussed open question in both evolutionary biology (Maynard Smith and Szathmary, 1995; Keller, 1999; Michod, 1999; Hammerstein, 2003; Okasha, 2006) and artificial life (Bedau et al., 2000). The fundamental problem is that any group of cooperative types, whose members donate some component of individual fitness in order to benefit their group, is vulnerable to invasion by selfish cheats that reap the group benefits of cooperation without paying the individual cost. Various mechanisms by which cooperation can nevertheless evolve have been suggested, including so-called ‘group selection’ models in which population structure exists such that individuals spend part of their time in groups (rather than freely mixed in the whole population) before mixing in a migrant pool from which new groups are formed (Wilson, 1980). Although selection within any given group will always favour selfish individuals, groups with a higher proportion of cooperators are more productive and hence contribute more individuals to the migrant pool and the next generation of groups.
Such models have been found to allow cooperation to evolve or be preserved in a population, but only within certain narrow parameter ranges (discussed below). This clearly presents a problem when appealing to such mechanisms as a route towards the evolution of higher-level individuals. As we have argued in previous work (Powers et al., 2007), the conditions that allow group selection to be effective and control its strength need not be externally imposed. Instead, they may be products of individual characters and hence subject to individual selection. Specifically, if key parameters such as group size are subject to individual adaptation (for example, via production of extracellular matrix in a bacterial biofilm), then a process akin to niche construction (Odling-Smee et al., 2003) supporting the evolution of cooperation may occur whereby cooperative traits and those affecting the strength of group selection evolve concurrently.
However, this process would require the existence of an adaptive gradient, such that a small parameter change could increase the strength of group selection and consequently cooperation and individual fitness. In this paper we consider one important parameter, group size. Where groups are formed randomly, existing models have shown that increasing group size rapidly causes the measurable effect of group selection to reach zero (Wade, 1978). This then means that for a large region of parameter space, a small decrease in group size, e.g. by individual mutation, would have no effect on fitness, i.e. there would be no selective gradient towards smaller groups and increased individual fitness. However, the model developed and presented in this paper suggests that the rapid tendency of group selection effects to zero is a consequence of an assumption of directional within-group selection. Specifically, classical group selection models attempt to explain the global promotion of a cooperative allele that is driven extinct at within-group equilibrium through directional selection for a rival selfish allele. While this assumption of directional selection and hence competitive exclusion of types is commonplace in population genetics models, a competitive coexistence of two types is instead often permitted in ecological models. For example, the classical Lotka-Volterra competition equations allow for coexistence as well as exclusion (May, 1976). Such dynamics are of relevance to potential multi-species group selection scenarios; for example, during egalitarian major evolutionary transitions in which different unrelated individuals eventually form a new level of selection, and within extant multi-species consortia such as bacterial biofilms (Burmolle et al., 2006) in which group selection effects may be pertinent. 111Coexistence dynamics may also apply to single-species scenarios under certain regimes such as balancing selection (see below) and are hence relevant to group selection discussions more generally. Results presented in this paper show that where a within-group stable coexistence of types exists, measurable group selection effects are sustained over a much larger range of group sizes, potentially providing an individual adaptive gradient towards smaller groups that enhance group selection over a much larger range of parameter space.
The Limits of Group Selection
Group selection theory is typically understood as the idea that the differential productivity of groups can lead to changes in allele frequency in a gene pool, in a manner analogous to individual selection acting through the differential productivity of individuals (Wade, 1978; Wilson, 1980). In this paper we are concerned with type 1 group selection, which defines group fitness as the average fitness of the group members 222Type 2 defines group fitness as number of offspring groups. (Okasha, 2006, chap. 2). Models of this process therefore seek to investigate the effect of group structure on the evolution of individual social traits such as cooperation. In particular, a trait that is individually disadvantageous may nevertheless evolve if it has a positive effect on the group as a whole, an idea that dates back to Darwin (1871). Using the language of multilevel selection theory, such cooperative traits are selected against within-group, since they confer a fitness disadvantage relative to other group members, but are favoured under between-group selection, since they increase group productivity, i.e. average absolute group member fitness. Whether or not the trait spreads in the global gene pool then depends on the balance of these two selective forces.
Three factors are pertinent in determining the outcome of such a scenario. The first is the individual cost to group benefit ratio of the cooperative act; a large individual cost strengthens within-group selection against cooperation, while a large group benefit strengthens between-group selection. Secondly, a group mixing mechanism must exist so that the increased productivity of more cooperative groups may affect the global allele frequencies. We consider here a multi-generational variant of D.S. Wilson’s (1980) trait-group model, where a global mixing stage occurs every generations in which the progeny of all groups disperse, join a global migrant pool, and then reform new groups of random composition. Since cooperative groups will contain more individuals prior to dispersal, they will constitute a larger fraction of the global migrant pool, thereby biasing the global allele frequencies. How often global mixing occurs is therefore a crucial parameter, since groups must be mixed before the within-group equilibrium has been reached, otherwise there will be no difference in group productivity for selection to act on (assuming selection within all groups leads to the same group equilibrium (Wilson, 1992)).
The third factor, and the one which is most often commented on in the group selection literature, is the requirement for there to be variation in the groups’ allelic composition. The larger this variance, the greater the difference in group productivity (since within-group dynamics are assumed to be deterministic) and likewise the effect of group selection. Many theoretical models have shown that if group composition is random then very small initial group sizes are needed in order to produce the between-group variance necessary for a measurable effect (see (Wade, 1978) for a classical review). It is therefore usually concluded that some kind of non-random group formation is required in order for group selection to have any significant effect. The most common way that such assortative grouping is believed to occur in nature is through kin grouping, where the group members are related by descent from a common ancestor and are hence more similar to each other than to members of other groups. Consequently, kin selection (Hamilton, 1964a, b) is commonly seen as the only pertinent force in social evolution. However, even with only simple random sampling, the range of group sizes which produce a group selection effect may be strongly affected if the assumption of within-group competitive exclusion is relaxed. This is because possible between-group variance from sampling error is controlled by the frequency of the least frequent type in the migrant pool. In the competitive exclusion case the possible variance rapidly tends to zero as the selfish type approaches fixation. If, however, coexistence of the two types is possible at equilibrium then we would expect the possible variance to be maintained above zero even at within-group equilibrium frequencies.
Group Competition: Aggregation and Dispersal Through a Migrant Pool
In this paper, we consider a model of a population structure where individuals reproduce in groups for a number of generations. After this period of reproduction within groups, all groups disperse and their progeny mix freely together in a migrant pool. Thereafter, new groups are created from the individuals in the migrant pool, and the process repeats. This process is a multi-generational variant of the trait-group model of Wilson (1980, 1987), and corresponds closely to the “Haystack” model of Maynard Smith (1964). Examples of natural populations that fit this model particularly well include soil-dwelling populations of micro-organisms that are occasionally mixed together during rainstorms (Wilson, 1980, p. 22), and a type of desert leaf cutting ant which lives in colonies that periodically disband and which are founded by unrelated females from a mating swarm (Rissing et al., 1989). However, the issues explored by the general form of the model are broadly applicable to a wide range of population structures where individuals have the majority of their interactions with a subset of the population. An algorithmic description of such an aggregation and dispersal process is as follows:
Initialisation: Initialise the migrant pool with individuals;
Group formation (aggregation): Assign individuals in the migrant pool randomly to groups of size ;
Reproduction: Perform reproduction within groups for time-steps, as described in the next section;
Migrant pool formation (dispersal): Return the progeny of each group to the migrant pool;
Iteration: Repeat from step 2 onwards for a number of generations, .
Group competition occurs in a population structured in this fashion during the dispersal stage, since groups that have grown to a larger size will contribute more individuals to the migrant pool. Crucially, this means that cooperative traits that are individually disadvantageous but that benefit the group can potentially increase in frequency (see Figure 1), depending on the balance of within- and between-group selection.
A key factor in determining the balance between the levels of selection in any model of group selection is the variance in group composition (Wilson, 1980), for if there is no variance then there is nothing for selection to act on (Darwin, 1859). In particular, there must be variance in initial group composition which causes a variance in group size prior to dispersal. In the above model, this variance is generated through random sampling of individuals from the migrant pool, where the sample size corresponds to the initial group size. Since an increased sample size causes a decrease in between-sample variance, increasing the initial group size decreases between-group variance and hence the efficacy of group selection (Wade, 1978). It is therefore often assumed that the upper limits of group size which produce a non-zero group selection effect are very small.
However, in this paper we present results which suggest that this follows from an assumption of within-group dynamics that lead to the competitive exclusion of a cooperative type by its selfish counterpart at within-group equilibrium. Specifically, we are able to show that where within-group dynamics instead lead to a stable coexistence of types, then the range of initial group sizes over which an effect of group selection can be seen is much larger. This is due to the fact that since the cooperative type cannot be driven extinct, variance in group composition when sampling from the migrant pool is always possible. In the next section, we describe how both competitive exclusion and coexistence dynamics can be modelled within groups.
Within-group Dynamics: Competitive Exclusion Versus Coexistence
Classical models of group selection consider a scenario where a selfish type ultimately drives its cooperative counterpart extinct at within-group equilibrium. In particular, fitness functions of the following form, first proposed by Wright (1945) but subsequently used in a plethora of other models (Williams and Williams, 1957; Maynard Smith, 1964; Charnov and Krebs, 1975; Wilson, 1980, 1987), are typically used to model within-group selection:
In the above equations, and denote the per capita fitness of selfish and cooperative individuals within a group, respectively. Cooperators, whose proportion within the group is denoted by , confer a fitness benefit on every group member. Crucially, both types receive this benefit, while only cooperators pay a cost, represented by the selection coefficient against cooperation, . It is then clear that if these equations are iterated until equilibrium is reached then the selfish type will be driven to fixation within the group.
Both competitive exclusion and stable coexistence within-group dynamics can instead be modelled using the standard two-species symmetric Lotka-Volterra competition equations (e.g. (May, 1976)). For implementation purposes, we use the following difference equation as a discrete approximation:
In the above equation, is the biomass of species at time , and the intrinsic per capita growth rate. Each species has an intrinsic carrying capacity , which is then modified through interspecific interactions. Specifically, the per capita effect of species on species is given by , the coefficient of interaction. All such interactions are competitive in the above equation, as ensured by the negative sign and the stipulation that all . Similarly, denotes the negative density-dependant effect of species on itself that prevents unbounded exponential growth. This coefficient can be seen as representing crowding and can vary for different species.
We define selfish () and cooperative () strategies in the above equation through settings of the within- and between-type interaction coefficients. Specifically, a selfish type is defined as having a large negative per capita effect on both itself () and the other type (). A cooperative type is then defined as having correspondingly smaller per capita negative effects ( and ). A pure group of cooperators will therefore grow to a larger size than a pure group of selfish individuals, creating a group productivity differential on which selection can potentially act. However, within mixed-groups the selfish type will reach the larger frequency (provided that is not too large), since . In other words, cooperators are favoured by between-group selection, while selfish individuals are favoured under within-group selection, exactly as in a classic group selection scenario. It should be noted that our definition of cooperative behaviour corresponds to weak, rather than strong, altruism (Wilson, 1980). This follows because although cooperation confers a relative fitness disadvantage compared to a selfish individual within the same group, it nevertheless increases the absolute fitness of all group members, including the cooperator.
It is well known that a stable coexistence of both types occurs in such a model when competition for resources (space, food, etc.) between individuals of the same type is stronger than competition between individuals of different types. Such a case corresponds at the ecological level to species occupying different niches, i.e. only partially overlapping in their resource requirements (May, 1976). Conversely, if between-type competition is stronger than within-type competition then competitive exclusion of one type will occur. Between- and within-type competition are both modelled in the Lotka-Volterra equations through the settings of the interaction coefficients. Throughout this paper, we assume the following:
and , i.e. that cooperators have lower negative density-dependant effects on themselves and others;
and , the intrinsic per capita growth rates and carrying capacities respectively, are the same for both types.
Given these assumptions, competitive exclusion of the cooperative type occurs when , producing qualitatively similar dynamics to those of the traditional within-group selection equations (1) and (2). However, when then the cooperative type is maintained at within-group equilibrium at an above-zero frequency, i.e. a stable coexistence of types occurs. When the interaction coefficients of the cooperative strategy are fixed, the equilibrium frequency at which it is maintained then depends upon the settings of and , i.e. the magnitude of the negative density-dependant effects of selfish individuals on themselves and cooperators, respectively. In addition, the within-group equilibrium is reached more quickly the greater the effects of the selfish type. Although a Lotka-Volterra model such as this is typically interpreted at the ecological level as representing species interactions, it could also be interpreted as a model of allelic competition dynamics within a single species group. In particular, coexistence Lotka-Volterra dynamics are analogous to balancing selection for a stable allelic polymorphism within a group. Conversely, competitive exclusion of one species by another is analogous to directional selection driving one allele to fixation. The motivation for using the language of allelic competition is to facilitate comparison with classical group selection models, which consider competition between selfish and cooperative alleles.
Our use of the Lotka-Volterra equations in this paper should be contrasted from their use in community or ecosystem selection models (Wilson, 1992; Penn, 2003). Such models do not consider explicit cooperative and selfish types in the fashion of traditional group selection models. Instead, they examine the complex within-group dynamics that arise when a larger number of types are present. These complex dynamics can give rise to multiple within-group attractors, which can then provide a source of variation in their own right upon which selection can act (Penn, 2003). By contrast, in this paper we consider simple two-type within-group dynamics, where only a single group attractor exists (either coexistence or competitive exclusion, as discussed above). As far as we are aware, our use of the Lotka-Volterra equations to define explicit selfish and cooperative strategies is novel.
|Parameter||Value (competitive exclusion)||Value (coexistance)|
The parameter settings used for the Lotka-Volterra equations throughout this paper are shown in Table 1. Changing between competitive exclusion and coexistence within-group dynamics is achieved by simply switching over the values of and , since that determines whether and hence whether competitive exclusion occurs. The values of the interaction coefficients in Table 1 produce representative within- and between-group dynamics; other settings produce the same qualitative trends. In this section, we first present results using classical competitive exclusion dynamics, and then contrast these to results from the coexistence case.
Group Selection Dynamics in the Competitive Exclusion Case
The within-group dynamics for a group initialised with unit biomass of each type are shown in Figures 2(a) and 2(b), for the competitive exclusion case. Initially, both types are in their growth phase; their biomass is below the intrinsic type carrying capacity of 100. However, the selfish type grows at a faster rate, despite the fact that their intrinsic growth rates, , are the same. This is because of the greater negative density-dependant effect of the selfish type on cooperators, i.e. . Finally, since , the cooperative type is driven to extinction. Furthermore, as Figure 2(b) shows, the proportion of selfish individuals increases monotonically. Such behaviour is qualitatively identical to that of directional within-group selection for a selfish allele in classical group selection models (e.g. (Wright, 1945; Wilson, 1980)).
Now let us consider global dynamics under group selection in this competitive exclusion case. In order for group selection to operate through an aggregation and dispersal process, a difference in group size at the dispersal stage must exist. Figure 3 illustrates how final group size varies as a function of the time spent in the group prior to dispersal, for various starting frequencies of cooperators in groups of initial size 10. It can be seen from this graph that, using the parameters described in Table 1, groups with a greater proportion of cooperators do indeed grow to a larger size. In addition, the results for the coexistence case, where and are swapped, are quantitatively similar. These results therefore confirm that group selection can in principal operate, since there is a variation in group productivity on which selection can act.
To determine the magnitude of the effect of group selection, the aggregation and dispersal process was executed for 5000 iterations, which preliminary experimentation had shown to be a sufficient length of time for a global equilibrium to be reached, using the within-group parameter settings described in Table 1. Equation 3 was iterated 30 times in the reproduction stage, while the global population size was maintained at 5000. Initial group size was then varied from 1 to 100 inclusive, while the migrant pool was initialised with 50% of each type. The result of this process after 5000 aggregation and dispersal cycles is shown in Figure 4, where ‘effect of group selection’ on the -axis is defined as the difference between the frequency of the selfish type at within-group equilibrium and the global frequency of the selfish type after 5000 aggregation and dispersal cycles. Since the within-group equilibrium in the competitive exclusion case is the selfish type at 100%, the -axis equivilantely shows the global frequency of cooperators in this case. Furthermore, it should also be stressed that the within-group equilibrium is the equilibrium that would be reached in an unstructured population where there were no groups. The -axis therefore shows the effect that group structure is having on the outcome of evolution compared to that in an unstructured population.
There are two points to note from Figure 4. Firstly, increasing the initial group size decreases the effect of group selection, and consequently the global proportion of cooperators. In particular, for small group sizes, the cooperative type reaches global fixation (and remains there because we do not reintroduce types by mutation). However, for group sizes above 10, it is driven extinct. This follows from the fact that the between-group variance necessary for group selection to act is generated by random sampling from the migrant pool, and therefore rests on the existence of a small initial group size, as previously discussed.
The second, and a key point for this paper, is that the effect of group selection rapidly tends to zero as initial group size increases. Specifically, above a group size of 10, there is no measurable effect at all. Such a result may therefore make the idea of group selection acting on randomly formed groups seem rather implausible as a significant evolutionary pathway. However, the above results only consider the competitive exclusion case; in the coexistence case, the results are somewhat different, as shown in the following section.
The Efficacy of Group Selection under Coexistence Dynamics
Let us now consider the coexistence dynamics that arise from redefining the selfish strategy as and , i.e. by swapping the interaction coefficients over from the competitive exclusion case. Figures 5(a) and 5(b) show how the cooperative type is no longer driven extinct at the within-group equilibrium. In particular, the change in the frequency of the selfish type from an initialisation of 50% shows clear balancing selective dynamics resulting in the maintance of cooperation at an above-zero frequency. In other words, the result is a stable coexistence of cooperative and selfish types within a group.
Group selection dynamics under the aggregation and dispersal process are now as shown in the black curve in Figure 6. Crucially, in contrast to the competitive exclusion case (shown again in the dotted line), an effect of group selection can be seen over the entire range of group sizes examined. For example, in groups of initial size 50, group selection can be seen to still increase the global frequency of cooperation above the within-group equilibrium. The significance of this observation is that since the within-group equilibrium is the same equilibrium that would be reached in an unstructured population, these results show that group structure is having an effect on population dynamics across a wide range of group sizes.
Finally, to verify that this result is not an artefact of the particular values of and used, the same curves were plotted for a variety of other parameters. Figure 7 provides an example of this, where and in the competitive exclusion case, vice versa for the coexistence case. These parameters were chosen since they represent stronger within-group selection towards selfish behaviour in the coexistence case than in the previous example. Specifically, the within-group equilibrium frequency of the selfish type in the coexistence case is 98.04%, compared to 83.3% previously. The results in Figure 7 show that while an effect of group selection is still seen over a larger range of group sizes in the coexistence case, the magnitude of the effect is reduced compared to Figure 6. The reason for this is that variance in group composition is proportional to the frequency of cooperators in the migrant pool in this case, and hence to the corresponding within-group equilibrium frequency, as discussed in detail in the following section.
The results in the previous section demonstrate that where a stable coexistence of types occurs at within-group equilibrium, an effect of group selection on global frequencies can be seen over a much larger range of initial group sizes than in the competitive exclusion case. In particular, as group size increases in the competitive exclusion case, any measurable effect of group selection on the global frequency of cooperation rapidly tends to zero. By contrast, in the coexistence case, some effect on global frequencies is seen over the entire range of group sizes examined. It must be stressed that we do not make a particular claim about the magnitude of the effect for large group sizes. Rather, our model implies that there is some measurable effect on frequencies over a large range of group sizes; how large this effect may be will depend on the properties of the natural system under consideration. However, the fact that any effect of group selection still exists over a large range of parameters is significant, since it suggests that where within-group dynamics in nature are of the coexistence type, some effect of a group population structure may always be acting.
Coexistence dynamics allow group selection effects to be sustained over a larger range of group sizes because of the effect of migrant pool frequencies on between-group variance. In particular, because group formation constitutes random sampling from the migrant pool, initial between-group variance can be approximated by the binomial distribution, and is then given by , where is the proportion of the cooperative type in the migrant pool, the proportion of the selfish type, and S the initial group size (Wilson, 1980, p. 27). Since , it follows that between-group variance is proportional to the frequency of the least frequent type, i.e. variance is maximal when both types are of equal frequency, and zero when one type is at fixation. Therefore, where one type reaches global fixation then there can be no variance and hence no group selection. However, in the coexistence case, where one type cannot reach fixation, it follows that there must always be some variance and hence some possible effect of group selection. The fact that the variance is proportional to the frequency of the least frequent type is illustrated by the difference between Figures 6 and 7, where the lower within-group equilibrium frequency of cooperators in Figure 7 results in a reduced effect of group selection for large group sizes.
A further observation from Figure 6 is that a gradient towards an increased effect of group selection also exists over a larger range of group sizes in the coexistence case. Specifically, decreasing group size by a small amount yields an increase in the effect of group selection for groups of size 20 in the coexistence case. However, there is no gradient at this size in the competitive exclusion case. The significance of this is that increasing the effect of group selection increases average absolute individual fitness in the population, due to an increased global level of cooperation. If group size can be partly determined by individual traits (Powers et al., 2007) then this may provide an adaptive gradient towards smaller groups, increased levels of cooperation, and greater fitness. In the competitive exclusion case, however, such a gradient only exists over a much smaller range of group sizes. While numerical experimentation investigating whether either gradient can be followed by a series of small mutations will be the subject of a future study, the results presented here do suggest that the concurrent evolution of group size and cooperation is more plausible in cases where a stable coexistence of types within groups exists.
Any group selection process requires there to be a variation in group composition. In aggregation and dispersal style models, this variation arises through the random assignment of individuals from the migrant pool into groups. Consequently, it is often suggested that an effect of group selection on the global frequency of types will only be seen for very small initial group sizes. However, the models on which this claim is based typically only consider within-group dynamics that lead to the competitive exclusion of a cooperative type by its selfish counterpart.
In this paper, we suggest that within-group competitive exclusion dynamics may be an unnecessary assumption in a number of situations, including the modelling of multi-species consortia in bacterial biofilms and during egalitarian major transitions. A model has been presented which shows that where such coexistence dynamics are present within groups, the range of initial group sizes over which an effect of group selection can be seen is much larger. Consequently, the potential for an adaptive gradient towards smaller groups and increased cooperation exists over a much larger range of group sizes under coexistence dynamics. This is in sharp contrast to the competitive exclusion case, where the effects of group selection rapidly reach zero as initial group size increases, excluding the possibility of such a gradient for a large range of parameters. Our results suggest that, where group size can be influenced by individual traits, the evolution of smaller groups and increased cooperation is more plausible under coexistence dynamics. Such increased group cooperation is a vital component of many major transitions in evolution (Maynard Smith and Szathmary, 1995; Michod, 1999).
We have shown in this paper that the conventional conclusion that group selection effects can only be seen for very small groups rests on the assumption that within-group dynamics lead to competitive exclusion. If a within-group coexistence of competing types is instead permitted, then the range of group sizes over which an effect can be seen is much larger. This result follows from the fact that the variance in group composition upon which group selection acts is dependant not only on group size but also on the frequencies of types in the migrant pool. In particular, since neither type can be driven extinct under coexistence dynamics, there will always be some variance in group composition when sampling from the migrant pool, which can then be acted on by group selection. Thus, since it is not necessary to assume that within-group dynamics lead to competitive exclusion, this result shows that group selection can operate in a wider range of conditions than previously realised.
- Bedau et al. (2000) Bedau, M., McCaskill, J., Packard, N., Rasmussen, S., Adami, C., Green, D., Ikegami, T., Kaneko, K., and Ray, T. (2000). Open questions in artificial life. Artificial Life, 6:352–357.
- Burmolle et al. (2006) Burmolle, M., Webb, J. S., Rao, D., Hansen, L. H., Sorensen, S. J., and Kjelleberg, S. (2006). Enhanced biofilm formation and increased resistance towards antimicrobial agents and bacterial invasion caused by synergistic interactions in multi-species biofilms. Appl Environ Microbiol, 72:3916–23.
- Charnov and Krebs (1975) Charnov, E. L. and Krebs, J. R. (1975). The evolution of alarm calls: Altruism or manipulation? The American Naturalist, 109(965):107–112.
- Darwin (1859) Darwin, C. (1859). On the origin of species by means of natural selection, or the preservation of favoured races in the struggle for life. John Murray, London.
- Darwin (1871) Darwin, C. (1871). The Descent of Man, and Selection in Relation to Sex. John Murray, London.
- Hamilton (1964a) Hamilton, W. D. (1964a). The genetical evolution of social behaviour. i. Journal of Theoretical Biology, 7(1):1–16.
- Hamilton (1964b) Hamilton, W. D. (1964b). The genetical evolution of social behaviour. ii. Journal of Theoretical Biology, 7(1):17–52.
- Hammerstein (2003) Hammerstein, P. (2003). Genetic and Cultural Evolution of Cooperation. MIT Press.
- Keller (1999) Keller, L., editor (1999). Levels of Selection in Evolution. Monographs in behavior and ecology. Princeton University Press.
- May (1976) May, R. M., editor (1976). Theoretical Ecology: Principles and Applications. Blackwell Scientific Publications, Oxford.
- Maynard Smith (1964) Maynard Smith, J. (1964). Group selection and kin selection. Nature, 201:1145–1147.
- Maynard Smith and Szathmary (1995) Maynard Smith, J. and Szathmary, E. (1995). Major Transitions in Evolution. Spektrum.
- Michod (1999) Michod, R. E. (1999). Darwinian Dynamics: Evolutionary Transitions in Fitness and Individuality. Princeton University Press.
- Odling-Smee et al. (2003) Odling-Smee, F. J., Laland, K. N., and Feldman, M. W. (2003). Niche construction: the neglected process in evolution. Monographs in population biology; no. 37. Princeton University Press.
- Okasha (2006) Okasha, S. (2006). Evolution and the Levels of Selection. Clarendon Press.
- Penn (2003) Penn, A. S. (2003). Modelling artificial ecosystem selection: A preliminary investigation. In Proceedings of the 7th European Conference on Artificial Life (ECAL 2003).
- Powers et al. (2007) Powers, S. T., Penn, A. S., and Watson, R. A. (2007). Individual selection for cooperative group formation. In Proceedings of the 9th European Conference on Artificial Life (ECAL 2007).
- Rissing et al. (1989) Rissing, S. W., Pollock, G. B., Higgins, M. R., Hagen, R. H., and Smith, D. R. (1989). Foraging specialization without relatedness or dominance among co-founding ant queens. Nature, 338:420–422.
- Wade (1978) Wade, M. J. (1978). A critical review of the models of group selection. The Quarterly Review of Biology, 53:101–114.
- Williams and Williams (1957) Williams, G. C. and Williams, D. C. (1957). Natural selection of individually harmful social adaptations among sibs with special reference to social insects. Evolution, 11:32–39.
- Wilson (1980) Wilson, D. S. (1980). The Natural Selection of Populations and Communities. Benjamin/Cummings, California.
- Wilson (1987) Wilson, D. S. (1987). Altruism in mendelian populations derived from sibling groups: The haystack model revisited. Evolution, 41(5):1059–1070.
- Wilson (1992) Wilson, D. S. (1992). Complex interactions in metacommunities, with implications for biodiversity and higher levels of selection. Ecology, 73(6):1984–2000.
- Wright (1945) Wright, S. (1945). Tempo and mode in evolution: A critical review. Ecology, 26:415–419.