Polymorphism in multilocus host parasite coevolutionary interactions

Numerous loci in host organisms are involved in parasite recognition, such as major histocompatibility complex (MHC) genes in vertebrates or genes involved in gene-for-gene (GFG) relationships in plants. Diversity is commonly observed at such loci and at corresponding loci encoding antigenic molecules in parasites. Multilocus theoretical models of host-parasite coevolution predict that polymorphism is more likely than in single-locus interactions because recurrent coevolutionary cycles are sustained by indirect frequency-dependent selection as rare genotypes have a selective advantage. These cycles are stabilized by direct frequency-dependent selection, resulting from repeated reinfection of the same host by a parasite, a feature of most diseases. Here, it is shown that for realistically small costs of resistance and virulence, polycyclic disease and high autoinfection rates, stable polymorphism of all possible genotypes is obtained in parasite populations. Two types of epistatic interactions between loci tend to increase the parameter space in which stable polymorphism can occur with all possible host and parasite genotypes. In the parasite, the marginal cost of each additional virulence allele should increase, while in the host, the marginal cost of each additional resistance allele should decrease. It is therefore predicted that GFG polymorphism will be stable (and hence detectable) when there is partial complementation of avirulence genes in the parasite and of resistance genes in the host.     

The long-term evolution of multilocus traits under frequency-dependent disruptive selection

Frequency-dependent disruptive selection is widely recognized as an important source of genetic variation. Its evolutionary consequences have been extensively studied using phenotypic evolutionary models, based on quantitative genetics, game theory, or adaptive dynamics. However, the genetic assumptions underlying these approaches are highly idealized and, even worse, predict different consequences of frequency-dependent disruptive selection. Population genetic models, by contrast, enable genotypic evolutionary models, but traditionally assume constant fitness values. Only a minority of these models thus addresses frequency-dependent selection, and only a few of these do so in a multilocus context. An inherent limitation of these remaining studies is that they only investigate the short-term maintenance of genetic variation. Consequently, the long-term evolution of multilocus characters under frequency-dependent disruptive selection remains poorly understood. We aim to bridge this gap between phenotypic and genotypic models by studying a multilocus version of Levene's soft-selection model. Individual-based simulations and deterministic approximations based on adaptive dynamics theory provide insights into the underlying evolutionary dynamics. Our analysis uncovers a general pattern of polymorphism formation and collapse, likely to apply to a wide variety of genetic systems: after convergence to a fitness minimum and the subsequent establishment of genetic polymorphism at multiple loci, genetic variation becomes increasingly concentrated on a few loci, until eventually only a single polymorphic locus remains. This evolutionary process combines features observed in quantitative genetics and adaptive dynamics models, and it can be explained as a consequence of changes in the selection regime that are inherent to frequency-dependent disruptive selection. Our findings demonstrate that the potential of frequency-dependent disruptive selection to maintain polygenic variation is considerably smaller than previously expected.     

ADAPTIVE DYNAMICS IN ALLELE SPACE: EVOLUTION OF GENETIC POLYMORPHISM BY SMALL MUTATIONS IN A HETEROGENEOUS ENVIRONMENT

We demonstrate how a genetic polymorphism of distinctly different alleles can develop during long-term frequency-dependent evolution in an initially monomorphic diploid population, if mutations have only small phenotypic effect. As a specific example, we use a version of Levene's (1953) soft selection model, where stabilizing selection acts on a continuous trait within each of two habitats. If the optimal phenotypes within the habitats are sufficiently different, then two distinctly different alleles evolve gradually from a single ancestral allele. In a wide range of parameter values, the two locally optimal phenotypes will be realized by one of the homozygotes and the heterozygote, rather than by the two homozygotes. Unlike in the haploid analogue of the model, there can be multiple polymorphic evolutionary attractors with different probabilities of convergence. Our results differ from the population genetic models of short-term evolution in two aspects: (1) a polymorphism that is population genetically stable may be invaded by a new mutant allele and, as a consequence, the population may fall back to monomorphism, (2) long-term evolution by allele substitutions may lead from a population where polymorphism is not possible into one where polymorphism is possible.     

Stability of genetic polymorphism in host-parasite interactions

Allelic diversity is common at host loci involved in parasite recognition, such as the major histocompatibility complex in vertebrates or gene-for-gene relationships in plants, and in corresponding loci encoding antigenic molecules in parasites. Diverse factors have been proposed in models to account for genetic polymorphism in host-parasite recognition. Here, a simple but general theory of host-parasite coevolution is developed. Coevolution implies the existence of indirect frequency-dependent selection (FDS), because natural selection on the host depends on the frequency of a parasite gene, and vice versa. It is shown that polymorphism can be maintained in both organisms only if there is negative, direct FDS, such that the strength of natural selection for the host resistance allele, the parasite virulence allele or both declines with increasing frequency of that allele itself. This condition may be fulfilled if the parasite has more than one generation in the same host individual, a feature which is common to most diseases. It is argued that the general theory encompasses almost all factors previously proposed to account for polymorphism at corresponding host and parasite loci, including those controlling gene-for-gene interactions.     

Ploidally antagonistic selection maintains stable genetic polymorphism

Understanding the maintenance of genetic variation in the face of selection remains a key issue in evolutionary biology. One potential mechanism for the maintenance of genetic variation is opposing selection during the diploid and haploid stages of biphasic life cycles universal among eukaryotic sexual organisms. If haploid and diploid gene expression both occur, selection can act in each phase, potentially in opposing directions. In addition, sex-specific selection during haploid phases is likely simply because male and female gametophytes/gametes tend to have contrasting life histories. We explored the potential for the maintenance of a stable polymorphism under ploidally antagonistic as well as sex-specific selection. Furthermore, we examined the role of the chromosomal location of alleles (autosomal or sex-linked). Our analyses show that the most permissible conditions for the maintenance of polymorphism occur under negative ploidy-by-sex interactions, where stronger selection for an allele in female than male diploids is coupled with weaker selection against the allele in female than male haploids. Such ploidy-by-sex interactions also promote allele frequency differences between the sexes. With constant fitness, ploidally antagonistic selection can maintain stable polymorphisms for autosomal and X-linked genes but not for Y-linked genes. We discuss the implications of our results and outline a number of biological settings where the scenarios modeled may apply.     

The evolution of genetic diversity.

The existence within natural populations of large amounts of genetic variation in molecules and morphology presents an evolutionary problem. The 'neutralist' solution to this problem, that the variation is usually unimportant to the organism displaying it, has now lost much of its strength. Interpretations that assume widespread heterozygous advantage also face serious difficulties. A resolution is possible in terms of frequency-dependent selection by predators, parasites and competitors. The evidence for pervasive frequency-dependent selection is now very strong. It appears to follow naturally from the behaviour of predators, from the evolutionary lability of parasites, from the ecology of competition and, at the molecular level, from the phenomena of enzyme kinetics. Such selection can explain the maintenance not only of conventional polymorphism but also of continuous variation in both molecular and morphological characters. It can account for the occurrence of diversity within groups of haploid and self-fertilizing organisms, and for the evolution of differences between individuals in their systems of genetic control.     

Discrete polymorphisms due to disruptive selection on a continuous trait--I: the one-locus case

We have investigated, numerically and analytically, long-term evolution under frequency-dependent disruptive selection of a continuous trait varying in a finite range and controlled by one diploid mendelian locus. We found that evolution converges towards a unique long-term equilibrium where only two extreme phenotypes are present with frequencies identical to those of the mixed strategy that would be the unique ESS of the game defined by the basic fitness function of the model. As long as this precise phenotypic composition is preserved, any genetic configuration of the polymorphism is equally acceptable (selectively neutral) at the equilibrium. Thus the number of alleles and their dominance pattern may vary considerably among different equilibrium populations. If genetic expression of the trait is variable but the amount of variability is genetically modifiable, disruptive selection, acting on such modifiers, produces a steady increase of expression variability before the equilibrium is attained. In this case a population at the long-term equilibrium might even be genetically monomorphic, with the phenotypic dimorphism resulting from purely random individual variation.     

Pathogen resistance and genetic variation at MHC loci

Abstract.— Balancing selection in the form of heterozygote advantage, frequency-dependent selection, or selection that varies in time and/or space, has been proposed to explain the high variation at major histocompatibility complex (MHC) genes. Here the effect of variation of the presence and absence of pathogens over time on genetic variation at multiallelic loci is examined. In the basic model, resistance to each pathogen is conferred by a given allele, and this allele is assumed to be dominant. Given that s is the selective disadvantage for homozygotes (and heterozygotes) without the resistance allele and the proportion of generations, which a pathogen is present, is e , fitnesses for homozygotes become (1 — s )^(n-1)e and the fitnesses for heterozygotes become (1 — s )^(n-2)e, where n is the number of alleles. In this situation, the conditions for a stable, multiallelic polymorphism are met even though there is no intrinsic heterozygote advantage. The distribution of allele frequencies and consequently heterozygosity are a function of the autocorrelation of the presence of the pathogen in subsequent generations. When there is a positive autocorrelation over generations, the observed heterozygosity is reduced. In addition, the effects of lower levels of selection and dominance and the influence of genetic drift were examined. These effects were compared to the observed heterozygosity for two MHC genes in several South American Indian samples. Overall, resistance conferred by specific alleles to temporally variable pathogens may contribute to the observed polymorphism at MHC genes and other similar host defense loci.     

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.     

COMPETITION AND GENOTYPE-BY-ENVIRONMENT INTERACTION IN NATURAL BREEDING SUBSTRATES OF DROSOPHILA

Although empirical studies frequently suggest that genotype-by-environment (G X E) interaction can maintain genetic variation, very few data are available to test for the specific conditions necessary for the existence of a protected polymorphism (i.e., the property of persistence of an allele even when initially rare). Drosophila species live in patchy environments and their local population structure may be characterized to some extent by Levene's migration pattern, namely by a single pool of individuals that presumably mate at random and breed on discrete and ephemeral resources. We present here a field experiment that links Drosophila ecology and population genetics, which used the alcohol dehydrogenase (Adh) and α-glycerophosphate dehydrogenase (αGpdh) polymorphic loci in D. melanogaster flies raised from Opuntia ficus-indica fruits (prickly pears). The results show that there is density-dependent mortality in those fruits with a relatively high number of larvae (i.e., selection is “soft”) and suggest that there is differential viability for αGpdh genotypes. Additionally, a pattern of G X E interaction for fitness values, which is fully compatible with the theoretical conditions required for the existence of a protected polymorphism, was found after weighting the fitness estimates by the relative contribution that each fruit makes to the total adult population. The strong association between Adh^S and αGpdh^F alleles suggests that the occurrence of the common cosmopolitan inversion In(2L)t in the population might be responsible for the negative frequency-dependent selection predicted by Levene's model when genetic variation persists in heterogeneous environments.     

MHC allele frequency distributions under parasite-driven selection: A simulation model

Background  

The extreme polymorphism that is observed in major histocompatibility complex (MHC) genes, which code for proteins involved in recognition of non-self oligopeptides, is thought to result from a pressure exerted by parasites because parasite antigens are more likely to be recognized by MHC heterozygotes (heterozygote advantage) and/or by rare MHC alleles (negative frequency-dependent selection). The Ewens-Watterson test (EW) is often used to detect selection acting on MHC genes over the recent history of a population. EW is based on the expectation that allele frequencies under balancing selection should be more even than under neutrality. We used computer simulations to investigate whether this expectation holds for selection exerted by parasites on host MHC genes under conditions of heterozygote advantage and negative frequency-dependent selection acting either simultaneously or separately.     

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.     

The evolution of genetic architecture under frequency-dependent disruptive selection

We propose a model to analyze a quantitative trait under frequency-dependent disruptive selection. Selection on the trait is a combination of stabilizing selection and intraspecific competition, where competition is maximal between individuals with equal phenotypes. In addition, there is a density-dependent component induced by population regulation. The trait is determined additively by a number of biallelic loci, which can have different effects on the trait value. In contrast to most previous models, we assume that the allelic effects at the loci can evolve due to epistatic interactions with the genetic background. Using a modifier approach, we derive analytical results under the assumption of weak selection and constant population size, and we investigate the full model by numerical simulations. We find that frequency-dependent disruptive selection favors the evolution of a highly asymmetric genetic architecture, where most of the genetic variation is concentrated on a small number of loci. We show that the evolution of genetic architecture can be understood in terms of the ecological niches created by competition. The phenotypic distribution of a population with an adapted genetic architecture closely matches this niche structure. Thus, evolution of the genetic architecture seems to be a plausible way for populations to adapt to regimes of frequency-dependent disruptive selection. As such, it should be seen as a potential evolutionary pathway to discrete polymorphisms and as a potential alternative to other evolutionary responses, such as the evolution of sexual dimorphism or assortative mating.     

Maintenance of handedness polymorphism in humans: a frequency-dependent selection model

Frequency-dependent selection is an important process in the maintenance of genetic variation in fitness. In humans, it has been proposed that the polymorphism of handedness is maintained by negative frequency-dependent selection, through a strategic advantage of left-handers in fighting interactions. Using simple mathematical models, we explore: (1) whether it is possible to predict the range of left-handedness frequencies observed in human populations by the frequency and the violence of fighting interactions; (2) the consequences of the sex differences in the probability of transmission of hand preference to offspring. We show that a wide range of values of the frequency of left-handers can be obtained with realistic changes of the parameters values. Our models reinforce the idea that negative frequency-dependence may have played a role in maintaining left-handedness in human populations, and provide further support for the importance of fighting interactions in the evolution of hand preference. Moreover, they suggest an explanation for the occurrence of left-handedness among women in this context, namely an indirect selective advantage through their male offspring.     

Allelic Genealogy under Overdominant and Frequency-Dependent Selection and Polymorphism of Major Histocompatibility Complex Loci

To explain the long-term persistence of polymorphic alleles (trans-specific polymorphism) at the major histocompatibility complex (MHC) loci in rodents and primates, a computer simulation study was conducted about the coalescence time of different alleles sampled under various forms of selection. At the same time, average heterozygosity, the number of alleles in a sample, and the rate of codon substitution were examined to explain the mechanism of maintenance of polymorphism at the MHC loci. The results obtained are as follows. (1) The coalescence time for neutral alleles is too short to explain the trans-specific polymorphism at the MHC loci. (2) Under overdominant selection, the coalescence time can be tens of millions of years, depending on the parameter values used. The average heterozygosity and the number of alleles observed are also high enough to explain MHC polymorphism. (3) The pathogen adaptation model proposed by Snell is incapable of explaining MHC polymorphism, since the coalescence time for this model is too short and the expected heterozygosity and the expected number of alleles are too small. (4) From the mathematical point of view, the minority advantage model of frequency-dependent selection is capable of explaining a high degree of polymorphism and trans-specific polymorphism. (5) The molecular mimicry hypothesis also gives a sufficiently long coalescence time when the mutation rate is low in the host but very high in the parasite. However, the expected heterozygosity and the expected number of alleles tend to be too small. (6) Consideration of the molecular mechanism of the function of MHC molecules and other biological observations suggest that the most important factor for the maintenance of MHC polymorphism is overdominant selection. However, some experiments are necessary to distinguish between the overdominance and frequency-dependent selection hypotheses.     

Sexually antagonistic cytonuclear fitness interactions in Drosophila melanogaster

Theoretical and empirical studies have shown that selection cannot maintain a joint nuclear-cytoplasmic polymorphism within a population except under restrictive conditions of frequency-dependent or sex-specific selection. These conclusions are based on fitness interactions between a diploid autosomal locus and a haploid cytoplasmic locus. We develop a model of joint transmission of X chromosomes and cytoplasms and through simulation show that nuclear-cytoplasmic polymorphisms can be maintained by selection on X-cytoplasm interactions. We test aspects of the model with a "diallel" experiment analyzing fitness interactions between pairwise combinations of X chromosomes and cytoplasms from wild strains of Drosophila melanogaster. Contrary to earlier autosomal studies, significant fitness interactions between X chromosomes and cytoplasms are detected among strains from within populations. The experiment further demonstrates significant sex-by-genotype interactions for mtDNA haplotype, cytoplasms, and X chromosomes. These interactions are sexually antagonistic--i.e., the "good" cytoplasms in females are "bad" in males--analogous to crossing reaction norms. The presence or absence of Wolbachia did not alter the significance of the fitness effects involving X chromosomes and cytoplasms but tended to reduce the significance of mtDNA fitness effects. The negative fitness correlations between the sexes demonstrated in our empirical study are consistent with the conditions that maintain cytoplasmic polymorphism in simulations. Our results suggest that fitness interactions with the sex chromosomes may account for some proportion of cytoplasmic variation in natural populations. Sexually antagonistic selection or reciprocally matched fitness effects of nuclear-cytoplasmic genotypes may be important components of cytonuclear fitness variation and have implications for mitochondrial disease phenotypes that differ between the sexes.     

A theoretical investigation of sympatric evolution of temporal reproductive isolation as illustrated by marine broadcast spawners

Recent theory suggests that frequency-dependent disruptive selection in combination with assortative mating can lead to the establishment of reproductive isolation in sympatry. Here we explore how temporal variation in reproduction might simultaneously generate both disruptive selection and assortative mating, and result in sympatric speciation. The conceptual framework of the model may be applicable to biological systems with negative frequency-dependent selection, such as marine broadcast spawners or systems with pollinator limitation. We present a model that is motivated by recent findings in marine broadcast spawners and is parameterized with data from the Montastraea annularis species complex. Broadcast spawners reproduce via external fertilization and synchronous spawning is required to increase the probability of successful fertilization, but empirical evidence shows that as density increases, so does the risk of polyspermy. Polyspermy is the fusion of multiple sperm with an egg at fertilization, a process that makes the embryo unviable. Synchrony can therefore also act as a source of negative density-dependent disruptive selection. Model analysis shows that the interaction between polyspermy and spawning synchrony can lead to temporal reproductive isolation in sympatry and that, more generally, increased density promotes maintenance of genetic variation.     

Role of pleiotropy in the evolution of a cryptic developmental variation in Caenorhabditis elegans

Robust biological systems are expected to accumulate cryptic genetic variation that does not affect the system output in standard conditions yet may play an evolutionary role once phenotypically expressed under a strong perturbation. Genetic variation that is cryptic relative to a robust trait may accumulate neutrally as it does not change the phenotype, yet it could also evolve under selection if it affects traits related to fitness in addition to its cryptic effect. Cryptic variation affecting the vulval intercellular signaling network was previously uncovered among wild isolates of Caenorhabditis elegans. Using a quantitative genetic approach, we identify a non-synonymous polymorphism of the previously uncharacterized nath-10 gene that affects the vulval phenotype when the system is sensitized with different mutations, but not in wild-type strains. nath-10 is an essential protein acetyltransferase gene and the homolog of human NAT10. The nath-10 polymorphism also presents non-cryptic effects on life history traits. The nath-10 allele carried by the N2 reference strain leads to a subtle increase in the egg laying rate and in the total number of sperm, a trait affecting the trade-off between fertility and minimal generation time in hermaphrodite individuals. We show that this allele appeared during early laboratory culture of N2, which allowed us to test whether it may have evolved under selection in this novel environment. The derived allele indeed strongly outcompetes the ancestral allele in laboratory conditions. In conclusion, we identified the molecular nature of a cryptic genetic variation and characterized its evolutionary history. These results show that cryptic genetic variation does not necessarily accumulate neutrally at the whole-organism level, but may evolve through selection for pleiotropic effects that alter fitness. In addition, cultivation in the laboratory has led to adaptive evolution of the reference strain N2 to the laboratory environment, which may modify other phenotypes of interest.     

Phenotypic assortative mating among the Peruvian Cashinahua

The persistence of behaviorally deleterious genes in the human population poses an interesting question for population genetics: If certain alleles at these loci are deleterious, why have they survived in the population? We consider evidence for phenotypic capacitance and/or frequency-dependent selection for an allele that has been putatively shown to have negative associations with human behaviors (the “short” 5-HTT promoter region allele) yet has persisted in human and nonhuman primate populations. Using data from the National Longitudinal Study of Adolescent Health, we compare sibling and twin variation in depression by 5-HTT genotype (specified in several ways) and investigate sibship-level cross-person gene-gene interactions. In support of the “orchid/dandelion” hypothesis, we find evidence that the short allele increases variation in phenotypes in response to environmental (or genetic) differences (i.e., acts as a perturbation of a phenotypic capacitor). Further, we also find some evidence that the effects of allelic variation at this locus are moderated by the genetic environment of the sibship unit (i.e., effects may be susceptible to frequency-dependent selection). We discuss implications of these findings for genetic models in general, specifically with respect to stable unit treatment value assumption violations (i.e., nonindependence of units of analysis).