Effect of balancing selection on spatial genetic structure within populations: theoretical investigations on the self-incompatibility locus and empirical studies in Arabidopsis halleri

The effect of selection on patterns of genetic structure within and between populations may be studied by contrasting observed patterns at the genes targeted by selection with those of unlinked neutral marker loci. Local directional selection on target genes will produce stronger population genetic structure than at neutral loci, whereas the reverse is expected for balancing selection. However, theoretical predictions on the intensity of this signal under precise models of balancing selection are still lacking. Using negative frequency-dependent selection acting on self-incompatibility systems in plants as a model of balancing selection, we investigated the effect of such selection on patterns of spatial genetic structure within a continuous population. Using numerical simulations, we tested the effect of the type of self-incompatibility system, the number of alleles at the self-incompatibility locus and the dominance interactions among them, the extent of gene dispersal, and the immigration rate on spatial genetic structure at the selected locus and at unlinked neutral loci. We confirm that frequency-dependent selection is expected to reduce the extent of spatial genetic structure as compared to neutral loci, particularly in situations with low number of alleles at the self-incompatibility locus, high frequency of codominant interactions among alleles, restricted gene dispersal and restricted immigration from outside populations. Hence the signature of selection on spatial genetic structure is expected to vary across species and populations, and we show that empirical data from the literature as well as data reported here on three natural populations of the herb Arabidopsis halleri confirm these theoretical results.     

Evolution of foraging behaviour in response to chronic malnutrition in Drosophila melanogaster

Chronic exposure to food of low quality may exert conflicting selection pressures on foraging behaviour. On the one hand, more active search behaviour may allow the animal to find patches with slightly better, or more, food; on the other hand, such active foraging is energetically costly, and thus may be opposed by selection for energetic efficiency. Here, we test these alternative hypotheses in Drosophila larvae. We show that populations which experimentally evolved improved tolerance to larval chronic malnutrition have shorter foraging path length than unselected control populations. A behavioural polymorphism in foraging path length (the rover-sitter polymorphism) exists in nature and is attributed to the foraging locus (for). We show that a sitter strain (for(s2)) survives better on the poor food than the rover strain (for(R)), confirming that the sitter foraging strategy is advantageous under malnutrition. Larvae of the selected and control populations did not differ in global for expression. However, a quantitative complementation test suggests that the for locus may have contributed to the adaptation to poor food in one of the selected populations, either through a change in for allele frequencies, or by interacting epistatically with alleles at other loci. Irrespective of its genetic basis, our results provide two independent lines of evidence that sitter-like foraging behaviour is favoured under chronic larval malnutrition.     

Parasite-mediated selection drives an immunogenetic trade-off in plains zebras (Equus quagga)

Pathogen evasion of the host immune system is a key force driving extreme polymorphism in genes of the major histocompatibility complex (MHC). Although this gene family is well characterized in structure and function, there is still much debate surrounding the mechanisms by which MHC diversity is selectively maintained. Many studies have investigated relationships between MHC variation and specific pathogens, and have found mixed support for and against the hypotheses of heterozygote advantage, frequency-dependent or fluctuating selection. Few, however, have focused on the selective effects of multiple parasite types on host immunogenetic patterns. Here, we examined relationships between variation in the equine MHC gene, ELA-DRA, and both gastrointestinal (GI) and ectoparasitism in plains zebras (Equus quagga). Specific alleles present at opposing population frequencies had antagonistic effects, with rare alleles associated with increased GI parasitism and common alleles with increased tick burdens. These results support a frequency-dependent mechanism, but are also consistent with fluctuating selection. Maladaptive GI parasite ‘susceptibility alleles’ were reduced in frequency, suggesting that these parasites may play a greater selective role at this locus. Heterozygote advantage, in terms of allele mutational divergence, also predicted decreased GI parasite burden in genotypes with a common allele. We conclude that an immunogenetic trade-off affects resistance/susceptibility to parasites in this system. Because GI and ectoparasites do not directly interact within hosts, our results uniquely show that antagonistic parasite interactions can be indirectly modulated through the host immune system. This study highlights the importance of investigating the role of multiple parasites in shaping patterns of host immunogenetic variation.     

Impact of density and disease on frequency-dependent selection and genetic polymorphism: experiments with stripe rust and wheat

Frequency-dependent disease impacts may contribute to the maintenance of genetic diversity and sexual reproduction in plant populations. In earlier work with experimental wheat (Triticum aestivum) populations at a single density, we found that stripe rust (caused by Puccinia striiformis) created frequency-dependent selection on its host but competitive interactions between host genotypes reduced the potential for disease to maintain genetic polymorphisms in this highly self-pollinated species; the weaker competitor actually exhibited positive disease-mediated frequency-dependent selection. Based on these results we predicted that at low density, where the overall level of competition is lower, disease would have a stronger impact relative to competition and thus be more likely to maintain genetic polymorphisms; at low densities the greatest effect of disease for negative frequency-dependent selection should be seen in the weak competitor. Here we report on results with wheat stripe rust in which we altered both the frequency and density of host genotypes in factorial combinations of two-way mixtures where each host genotype was attacked by its own specialized race of rust. Within each density disease levels increased with genotype frequencies, creating frequency-dependent disease attack at all densities. Similarly, disease created negative frequency-dependent selection on its host at all densities, as a genotype’s fitness was often greater at low than high frequency when disease was present. Disease levels increased with plant density in 1997 but decreased in 1998. While increasing plant density reduced absolute fitness, presumably as a result of increased competition, a genetic polymorphism was not more likely to be maintained at low than high density as we had predicted. Within each density, the impact of disease was insufficient to reverse the slope of the relationship between absolute fitness and planted frequency from positive to negative for the less competitive host genotype, thus preventing the maintenance of a genetic polymorphism.     

Models of Frequency-Dependent Selection with Mutation from Parental Alleles

Frequency-dependent selection (FDS) remains a common heuristic explanation for the maintenance of genetic variation in natural populations. The pairwise-interaction model (PIM) is a well-studied general model of frequency-dependent selection, which assumes that a genotype’s fitness is a function of within-population intergenotypic interactions. Previous theoretical work indicated that this type of model is able to sustain large numbers of alleles at a single locus when it incorporates recurrent mutation. These studies, however, have ignored the impact of the distribution of fitness effects of new mutations on the dynamics and end results of polymorphism construction. We suggest that a natural way to model mutation would be to assume mutant fitness is related to the fitness of the parental allele, i.e., the existing allele from which the mutant arose. Here we examine the numbers and distributions of fitnesses and alleles produced by construction under the PIM with mutation from parental alleles and the impacts on such measures due to different methods of generating mutant fitnesses. We find that, in comparison with previous results, generating mutants from existing alleles lowers the average number of alleles likely to be observed in a system subject to FDS, but produces polymorphisms that are highly stable and have realistic allele-frequency distributions.     

Molecular population genetics of the OBP83 genomic region in Drosophila subobscura and D. guanche: contrasting the effects of natural selection and gene arrangement expansion in the patterns of nucleotide variation

Chromosomal inversion polymorphism play a major role in the evolutionary dynamics of populations and species because of their effects on the patterns of genetic variability in the genomic regions within inversions. Though there is compelling evidence for the adaptive character of chromosomal polymorphisms, the mechanisms responsible for their maintenance in natural populations is not fully understood. For this type of analysis, Drosophila subobscura is a good model species as it has a rich and extensively studied chromosomal inversion polymorphism system. Here, we examine the patterns of DNA variation in two natural populations segregating for chromosomal arrangements that differentially affect the surveyed genomic region; in particular, we analyse both nucleotide substitutions and insertion/deletion variations in the genomic region encompassing the odorant-binding protein genes Obp83a and Obp83b (Obp83 region). We show that the two main gene arrangements are genetically differentiated, but are consistent with a monophyletic origin of inversions. Nevertheless, these arrangements interchange some genetic information, likely by gene conversion. We also find that the frequency spectrum-based tests indicate that the pattern of nucleotide variation is not at equilibrium; this feature probably reflects the rapid increase in the frequency of the new gene arrangement promoted by positive selection (that is an adaptive change). Furthermore, a comparative analysis of polymorphism and divergence patterns reveals a relaxation of the functional constraints at the Obp83b gene, which might be associated with particular ecological or demographic features of the Canary island endemic species D. guanche     

Negative frequency-dependent selection is intensified at higher population densities in protist populations

Natural populations of free-living protists often exhibit high-levels of intraspecific diversity, yet this is puzzling as classic evolutionary theory predicts dominance by genotypes with high fitness, particularly in large populations where selection is efficient. Here, we test whether negative frequency-dependent selection (NFDS) plays a role in the maintenance of diversity in the marine flagellate Oxyrrhis marina using competition experiments between multiple pairs of strains. We observed strain-specific responses to frequency and density, but an overall signature of NFDS that was intensified at higher population densities. Because our strains were not selected a priori on the basis of particular traits expected to exhibit NFDS, these data represent a relatively unbiased estimate of the role for NFDS in maintaining diversity in protist populations. These findings could help to explain how bloom-forming plankton, which periodically achieve exceptionally high population densities, maintain substantial intraspecific diversity.     

Evolutionary Genetics of Nodule Bacteria: Molecular and Population Aspects

The molecular analysis of the genetic systems controlling the main stages of nodule bacteria (rhizobia) interaction with a legume host (signaling at early stages and symbiotic nitrogen fixation) has shown that the widespread recombination of genetic material in free-living ancestors of rhizobia was an important factor in the evolution of these systems. These recombinations could be conditioned by a high content of repeated DNA sequences and the IS elements in the rhizobial genome. A high recombination activity of rhizobia is manifested in the panmictic structure of their populations, which is associated with frequency-dependent selection favoring rare recombinants. This selection is realized through the competition of virulent strains for the nodule formation and can be controlled by the genes whose expression depends on population density (via the quorum sensing mechanism). A high degree of panmixia in rhizobial populations is associated with their ecotypic polymorphism, manifested as the coexistence of symbiotic and nonsymbiotic strains. This type of polymorphism is caused by individual selection during the periodic changes of ecological niches (soil–plant host) in the rhizobia life cycle. The rhizobia–plant interaction stimulates selection in bacterial populations, which results in the increased levels of their heterogeneity and panmixia. The combination of individual and frequency-dependent selection types resulted in the high rates of symbiosis evolution and polyphyletic origin of diverse rhizobial species.     

Are reward polymorphisms subject to frequency‐ and density‐dependent selection? Evidence from a monoecious species pollinated by deceit

Abstract Most deceit‐pollinated species involve floral dimorphisms characterized by the presence of rewarding male flowers and nonrewarding female flowers. It has been proposed that this polymorphism establishes the conditions for the action of frequency dependent selection (FDS). The tendency of foraging animals to aggregate in areas of high resource density suggests that pollination efficiency and fruit production may be positively influenced by flower density (density dependent selection, DDS). In this paper we offer a graphical model describing the effects of FDS and DDS on a monoecious species pollinated by deceit. We test the FDS and DDS assumptions and the predictions of the model using field observations and experimental populations of Begonia gracilis in which population sex ratio and flower density were controlled. We found a marked effect of both FDS and DDS on pollinator visitation, fruit‐set, and on the probability of female flowers to setting fruits. We conclude that these two types of selection have had a strong influence on the evolution of deceit‐pollinated species.     

Impact of negative frequency-dependent selection on mating pattern and genetic structure: a comparative analysis of the S-locus and nuclear SSR loci in Prunus lannesiana var. speciosa

Mating processes of local demes and spatial genetic structure of island populations at the self-incompatibility (S-) locus under negative frequency-dependent selection (NFDS) were evaluated in Prunus lannesiana var. speciosa in comparison with nuclear simple sequence repeat (SSR) loci that seemed to be evolutionarily neutral. Our observations of local mating patterns indicated that male-female pair fecundity was influenced by not only self-incompatibility, but also various factors, such as kinship, pollen production and flowering synchrony. In spite of the mating bias caused by these factors, the NFDS effect on changes in allele frequencies from potential mates to mating pollen was detected at the S-locus but not at the SSR loci, although the changes from adult to juvenile cohorts were not apparent at any loci. Genetic differentiation and isolation-by-distance over various spatial scales were smaller at the S-locus than at the SSR loci, as expected under the NFDS. Allele-sharing distributions among the populations also had a unimodal pattern at the S-locus, indicating the NFDS effect except for alleles unique to individual populations probably due to isolation among islands, although this pattern was not exhibited by the SSR loci. Our results suggest that the NFDS at the S-locus has an impact on both the mating patterns and the genetic structure in the P. lannesiana populations studied.     

Extensive intragenic recombination and patterns of linkage disequilibrium at the CSN3 locus in European rabbit

Kappa-casein (CSN3) plays an important role in stabilising the Ca-sensitive caseins in the micelle. The European rabbit (Oryctolagus cuniculus) CSN3 has previously been shown to possess two alleles (A and B), which differ deeply in their intronic regions (indels of 100 and 1550 nucleotides in introns 1 and 4, respectively). Furthermore, a correlation between several reproductive performance traits and the different alleles was described. However, all these data were exclusively collected in rabbit domestic breeds, preventing a deeper understanding of the extensive polymorphism observed in the CSN3 gene. Additionally, the techniques available for the typing of both indel polymorphisms were until now not suitable for large-scale studies. In this report, we describe a simple, PCR-based typing method to distinguish rabbit CSN3 alleles. We analyse both ancient wild rabbit populations from the Iberian Peninsula and France, and the more recently derived English wild rabbits and domestic stocks. A new allele (C) showing another major indel (250 bp) in intron 1 was found, but exclusively detected in Iberian wild rabbits. In addition, our survey revealed the occurrence of new haplotypes in wild populations, suggesting that intragenic recombination is important in creating genetic diversity at this locus. This easy and low cost single-step PCR-based method results in an improvement over previous described techniques, can be easily set up in a routine molecular laboratory and would probably be a valuable tool in the management of rabbit domestic breeds.     

The influence of stochastic and selective forces in the population divergence of female colour polymorphism in damselflies of the genus Ischnura

Disentangling the relative importance and potential interactions of selection and genetic drift in driving phenotypic divergence of species is a classical research topic in population genetics and evolutionary biology. Here, we evaluate the role of stochastic and selective forces on population divergence of a colour polymorphism in seven damselfly species of the genus Ischnura, with a particular focus on I. elegans and I. graellsii. Colour-morph frequencies in Spanish I. elegans populations varied greatly, even at a local scale, whereas more similar frequencies were found among populations in eastern Europe. In contrast, I. graellsii and the other five Ischnura species showed little variation in colour-morph frequencies between populations. F(ST)-outlier analyses revealed that the colour locus deviated strongly from neutral expectations in Spanish populations of I. elegans, contrasting the pattern found in eastern European populations, and in I. graellsii, where no such discrepancy between morph divergence and neutral divergence could be detected. This suggests that divergent selection has been operating on the colour locus in Spanish populations of I. elegans, whereas processes such as genetic drift, possibly in combination with other forms of selection (such as negative frequency-dependent selection), appear to have been present in other regions, such as eastern Europe. Overall, the results indicate that both selective and stochastic processes operate on these colour polymorphisms, and suggest that the relative importance of factors varies between geographical regions.     

Global allele polymorphism indicates a high rate of allele genesis at a locus under balancing selection

When selection favours rare alleles over common ones (balancing selection in the form of negative frequency-dependent selection), a locus may maintain a large number of alleles, each at similar frequency. To better understand how allelic richness is generated and maintained at such loci, we assessed 201 sequences of the complementary sex determiner (csd) of the Asian honeybee (Apis cerana), sampled from across its range. Honeybees are haplodiploid; hemizygotes at csd develop as males and heterozygotes as females, while homozygosity is lethal. Thus, csd is under strong negative frequency-dependent selection because rare alleles are less likely to end up in the lethal homozygous form. We find that in A. cerana, as in other Apis, just a few amino acid differences between csd alleles in the hypervariable region are sufficient to trigger female development. We then show that while allelic lineages are spread across geographical regions, allelic differentiation is high between populations, with most csd alleles (86.3%) detected in only one sample location. Furthermore, nucleotide diversity in the hypervariable region indicates an excess of recently arisen alleles, possibly associated with population expansion across Asia since the last glacial maximum. Only the newly invasive populations of the Austral-Pacific share most of their csd alleles. In all, the geographic patterns of csd diversity in A. cerana indicate that high mutation rates and balancing selection act together to produce high rates of allele genesis and turnover at the honeybee sex locus, which in turn leads to its exceptionally high local and global polymorphism.Subject terms: Evolutionary genetics, Rare variants, Ecological genetics     

Plasmodium relictum infection and MHC diversity in the house sparrow (Passer domesticus)

Antagonistic coevolution between hosts and parasites has been proposed as a mechanism maintaining genetic diversity in both host and parasite populations. In particular, the high level of genetic diversity usually observed at the major histocompatibility complex (MHC) is generally thought to be maintained by parasite-driven selection. Among the possible ways through which parasites can maintain MHC diversity, diversifying selection has received relatively less attention. This hypothesis is based on the idea that parasites exert spatially variable selection pressures because of heterogeneity in parasite genetic structure, abundance or virulence. Variable selection pressures should select for different host allelic lineages resulting in population-specific associations between MHC alleles and risk of infection. In this study, we took advantage of a large survey of avian malaria in 13 populations of the house sparrow (Passer domesticus) to test this hypothesis. We found that (i) several MHC alleles were either associated with increased or decreased risk to be infected with Plasmodium relictum, (ii) the effects were population specific, and (iii) some alleles had antagonistic effects across populations. Overall, these results support the hypothesis that diversifying selection in space can maintain MHC variation and suggest a pattern of local adaptation where MHC alleles are selected at the local host population level.     

Genetic variation in the major histocompatibility complex of the European brown hare (Lepus europaeus) across distinct phylogeographic areas

The major histocompatibility complex is one of the best studied systems in vertebrates providing evidence for the long-term action of selection. Here, we examined the intra- and inter-population genetic diversity of the MHC class II DRB locus in European brown hare (Lepus europaeus) and correlated the results with genetic variability already estimated from the MHC DQA locus and from maternally (mitochondrial DNA (mtDNA)) and biparentally (allozymes, microsatellites) inherited loci. L. europaeus showed remarkable genetic polymorphism in both DQA and DRB1 loci. The Anatolian populations exhibited the highest genetic polymorphism for both loci. Balancing selection has established increased variability in the European populations despite the founder effects after the last glaciation. Different evolutionary rates were traced for DRB1 and DQA loci, as evidenced by the higher number of common DRB1 than DQA alleles and the greater differences between DRB1 alleles with common origin in comparison with DQA alleles. The high number of rare alleles with low frequencies detected implies that frequency-dependent selection drives MHC evolution in the brown hare through the advantage of rare alleles. Both loci were under the influence of positive selection within the peptide-binding region. The functional polymorphism, recorded as amino acid substitutions within the binding pockets, fell also within distinct geographic patterns, yet it was much narrower than the genetic polymorphism. We hypothesize that certain structural and functional characteristics of the binding pockets set limitations to the actual shape of genetic polymorphism in MHC.     

A trade-off between oxidative stress resistance and DNA repair plays a role in the evolution of elevated mutation rates in bacteria

The dominant paradigm for the evolution of mutator alleles in bacterial populations is that they spread by indirect selection for linked beneficial mutations when bacteria are poorly adapted. In this paper, we challenge the ubiquity of this paradigm by demonstrating that a clinically important stressor, hydrogen peroxide, generates direct selection for an elevated mutation rate in the pathogenic bacterium Pseudomonas aeruginosa as a consequence of a trade-off between the fidelity of DNA repair and hydrogen peroxide resistance. We demonstrate that the biochemical mechanism underlying this trade-off in the case of mutS is the elevated secretion of catalase by the mutator strain. Our results provide, to our knowledge, the first experimental evidence that direct selection can favour mutator alleles in bacterial populations, and pave the way for future studies to understand how mutation and DNA repair are linked to stress responses and how this affects the evolution of bacterial mutation rates.     

Evolutionary genetics of rhizobia: molecular and population aspects

The molecular analysis of the genetic systems controlling the main stages of nodule bacteria (rhizobia) interaction with a legume host (signaling at early stages and symbiotic nitrogen fixation) has shown that the widespread recombination of genetic material in free-living ancestors of rhizobia was an important factor in the evolution of these systems. These recombinations could be conditioned by a high content of repeated DNA sequences and the IS elements in the rhizobial genome. A high recombination activity of rhizobia is manifested in the panmictic structure of their populations, which is associated with frequency-dependent selection favoring rare recombinants. This selection is realized through the competition of virulent strains for the nodule formation and can be controlled by the genes whose expression depends on population density (via the quorum sensing mechanism). A high degree of panmixia in rhizobial populations is associated with their ecotypic polymorphism, manifested as the coexistence of symbiotic and nonsymbiotic strains. This type of polymorphism is caused by individual selection during the periodic changes of ecological niches (soil-plant host) in the rhizobia life cycle. The rhizobia-plant interaction stimulates selection in bacterial populations, which results in the increased levels of their heterogeneity and panmixia. The combination of individual and frequency-dependent selection types resulted in the high rates of symbiosis evolution and polyphyletic origin of diverse rhizobial species.     

Disease, frequency-dependent selection, and genetic polymorphisms: experiments with stripe rust and wheat

Abstract.— Pathogens have the potential to maintain genetic polymorphisms by creating frequency-dependent selection on their host. This can occur when a rare host genotype is less likely to be attacked by a pathogen (frequency-dependent disease attack) and has higher fitness at low frequency (negative frequency-dependent selection). In this study, we used wheat genotypes that were susceptible to different races of the pathogen Puccinia striiformis to test whether disease created frequency-selection on its host and whether such selection could maintain polymorphisms for resistance genes in the wheat populations. Four different two-way mixtures of wheat genotypes were planted at different frequencies in both the presence and absence of disease. Disease created frequency-dependent selection on its host in some populations. Unknown factors other than disease also created frequency-dependent selection in this system because, in some instances, rare genotype advantage was observed in the absence of disease. Although the pathogen created frequency-dependent selection on its host, this selection was not sufficient to maintain genetic polymorphism in the host populations. In all cases where frequency-dependent selection occurred only in the diseased plots, one of the two genotypes was predicted to dominate in the population and the same genotype was predicted to dominate in both the presence and absence of disease. Only in cases where frequency-dependent selection was not caused by disease was there evidence that genetic polymorphisms would be maintained in the population. The frequency-dependent selection described in this study is a consequence of epidemiological effects of disease and differs from the time-lagged frequency-dependent selection resulting from coevolution between hosts and parasites. The impact of this direct frequency-dependent selection on the maintenance of genetic polymorphisms in the host population is discussed.     

Interplay of Darwinian and frequency-dependent selection in the host-associated microbial populations

In order to analyze the microevolutionary processes in host-associated microorganisms, we simulated the dynamics of rhizobia populations composed of a parental strain and its mutants possessing the altered fitness within "plant-soil" system. The population dynamics was presented as a series of cycles (each one involves "soil-->rhizosphere-->nodules-->soil" succession) described using recurrent equations. For representing the selection and mutation pressures, we used a universal approach based on calculating the shifts in the genetic ratios of competing bacterial genotypes within the particular habitats and across several habitats. Analysis of the model demonstrated that a balanced polymorphism may be established in rhizobia population: mutants with an improved fitness do not supplant completely the parental strain while mutants with a decreased fitness may be maintained stably. This polymorphism is caused by a rescue of low-fitted genotypes via negative frequency-dependent selection (FDS) that is implemented during inoculation of nodules and balances the Darwinian selection that occurs during multiplication or extinction of bacteria at different habitats. The most diverse populations are formed if the rhizobia are equally successful in soil and nodules, while a marked preference for any of these habitats results in the decrease of diversity. Our simulation suggests that FDS can maintain the mutualistic rhizobia-legume interactions under the stress conditions deleterious for surviving the bacterial strains capable for intensive N2 fixation. Genetic consequences of releasing the modified rhizobia strains may be addressed using the presented model.     

FREQUENCY-DEPENDENT SUCCESS OF CHEATERS DURING FORAGING BOUTS MIGHT LIMIT THEIR SPREAD WITHIN COLONIES OF A SOCIALLY POLYMORPHIC SPIDER

Although several studies have demonstrated that frequency-dependent effects can promote the maintenance of cooperative behavior in microbes, experimental evidence of frequency-dependent effects in cooperative animal societies is rare. We staged mixed phenotype feeding bouts in the spider Anelosimus studiosus , which shows a within-population social polymorphism, to determine how phenotype frequency affects the foraging success of the social (cooperative) and asocial (cheater) phenotypes. Foraging performance was inferred from average change in percent mass for the respective phenotypes after staged group foraging events. We then performed a field census of multifemale colonies of A. studiosus to determine the phenotypic composition of naturally occurring colonies. Our data indicate that asocial (i.e., cheater) individuals experience negative frequency-dependent foraging success in staged foraging contests. Asocial individuals outperform social individuals when their representation is low, but lose this competitive advantage as their relative numbers increase. Naturally occurring colonies, on average, contained 58.33% social and 41.67% asocial individuals.