Red Queen Processes Drive Positive Selection on Major Histocompatibility Complex (MHC) Genes

Major Histocompatibility Complex (MHC) genes code for proteins involved in the incitation of the adaptive immune response in vertebrates, which is achieved through binding oligopeptides (antigens) of pathogenic origin. Across vertebrate species, substitutions of amino acids at sites responsible for the specificity of antigen binding (ABS) are positively selected. This is attributed to pathogen-driven balancing selection, which is also thought to maintain the high polymorphism of MHC genes, and to cause the sharing of allelic lineages between species. However, the nature of this selection remains controversial. We used individual-based computer simulations to investigate the roles of two phenomena capable of maintaining MHC polymorphism: heterozygote advantage and host-pathogen arms race (Red Queen process). Our simulations revealed that levels of MHC polymorphism were high and driven mostly by the Red Queen process at a high pathogen mutation rate, but were low and driven mostly by heterozygote advantage when the pathogen mutation rate was low. We found that novel mutations at ABSs are strongly favored by the Red Queen process, but not by heterozygote advantage, regardless of the pathogen mutation rate. However, while the strong advantage of novel alleles increased the allele turnover rate, under a high pathogen mutation rate, allelic lineages persisted for a comparable length of time under Red Queen and under heterozygote advantage. Thus, when pathogens evolve quickly, the Red Queen is capable of explaining both positive selection and long coalescence times, but the tension between the novel allele advantage and persistence of alleles deserves further investigation.     

Sexual selection and the evolutionary dynamics of the major histocompatibility complex

The genes of the major histocompatibility complex (MHC) are a key component of the adaptive immune system and among the most variable loci in the vertebrate genome. Pathogen-mediated natural selection and MHC-based disassortative mating are both thought to structure MHC polymorphism, but their effects have proven difficult to discriminate in natural systems. Using the first model of MHC dynamics incorporating both survival and reproduction, we demonstrate that natural and sexual selection produce distinctive signatures of MHC allelic diversity with critical implications for understanding host–pathogen dynamics. While natural selection produces the Red Queen dynamics characteristic of host–parasite interactions, disassortative mating stabilizes allele frequencies, damping major fluctuations in dominant alleles and protecting functional variants against drift. This subtle difference generates a complex interaction between MHC allelic diversity and population size. In small populations, the stabilizing effects of sexual selection moderate the effects of drift, whereas pathogen-mediated selection accelerates the loss of functionally important genetic diversity. Natural selection enhances MHC allelic variation in larger populations, with the highest levels of diversity generated by the combined action of pathogen-mediated selection and disassortative mating. MHC-based sexual selection may help to explain how functionally important genetic variation can be maintained in populations of conservation concern.     

Balancing selection versus allele and supertype turnover in MHC class II genes in guppies

Selection pressure from parasites is thought to be a major force shaping the extreme polymorphism of the major histocompatibility complex (MHC) genes, but the modes and consequences of selection remain unclear. Here, we analyse MHC class II and microsatellite diversity in 16 guppy (Poecilia reticulata) populations from two islands (Trinidad and Tobago) that have been separated for at least 10 ky. Within-population MHC diversity was high, but allele sharing was limited within islands and even lower between islands, suggesting relatively fast turnover of alleles. Allelic lineages strongly supported in phylogenetic analyses tended to be island-specific, suggesting rapid lineage sorting, and an expansion of an allelic lineage private to Tobago was observed. New alleles appear to be generated locally at a detectably high frequency. We did not detect a consistent signature of local adaptation, but F_ST outlier analysis suggested that balancing selection may be the more general process behind spatial variation in MHC allele frequencies in this system, particularly within Trinidad. We found no evidence for divergent allele advantage within populations, or for decreased genetic structuring of MHC supertypes compared to MHC alleles. The dynamic and complex nature of MHC evolution we observed in guppies, coupled with some evidence for balancing selection shaping MHC allele frequencies, are consistent with Red Queen processes of host-parasite coevolution.Subject terms: Population genetics, Evolutionary genetics     

Spatially and temporally fluctuating selection at non-MHC immune genes: evidence from TAP polymorphism in populations of brown trout (Salmo trutta, L.)

Temporal samples of Danish brown trout (Salmo trutta) from populations representing varying geographical scales were analysed using eight putatively neutral microsatellite loci and two microsatellite loci embedded in TAP genes (Transporter associated with Antigen Processing). These genes encode molecules that are central to the major histocompatibility complex (MHC) class I restricted antigen presentation and thus integral components in the adaptive immune system. As such, they could be influenced by selection, driven by pathogens and parasites in a manner similar to MHC genes. Analysis of allele frequencies at presumably neutral microsatellite loci revealed a temporally unstable population structure within regions, while the population structure was stable over time among regions. Analyses of the two TAP markers indicated an effect of selection at both a regional and micro-geographical spatial scale. Moreover, signals of divergent selection among temporal samples within localities suggest that selection also might fluctuate at a temporal scale. These results suggest that immune genes other than the classical MHC class I and II might be subject to selection and warrant further studies of functional polymorphism of such genes in natural populations.     

Between‐year variation of MHC allele frequencies in great reed warblers: selection or drift?

The major histocompatibility complex (MHC) genes are extremely polymorphic and this variation is assumed to be maintained by balancing selection. Cyclic interactions between pathogens and their hosts could generate such selection, and specific MHC alleles or heterozygosity at certain MHC loci have been shown to confer resistance against particular pathogens. Here we compare the temporal variation in allele frequencies of 23 MHC class I alleles with that of 23 neutral microsatellite markers in adult great reed warblers (a passerine bird) in nine successive cohorts. Overall, the MHC alleles showed a significantly higher variation in allele frequencies between cohorts than the microsatellite alleles, using a multi-variate genetic analysis (amova). The frequency of two specific MHC alleles, A3e (P = 0.046) and B4b (P = 0.0018), varied more between cohorts than expected from random, whereas none of the microsatellite alleles showed fluctuations exceeding the expectation from stochastic variation. These results imply that the variation in MHC allele frequencies between cohorts is not a result of demographic events, but rather an effect of selection favouring different MHC alleles in different years.     

Red Queen Dynamics with Non-Standard Fitness Interactions

Antagonistic coevolution between hosts and parasites can involve rapid fluctuations of genotype frequencies that are known as Red Queen dynamics. Under such dynamics, recombination in the hosts may be advantageous because genetic shuffling can quickly produce disproportionately fit offspring (the Red Queen hypothesis). Previous models investigating these dynamics have assumed rather simple models of genetic interactions between hosts and parasites. Here, we assess the robustness of earlier theoretical predictions about the Red Queen with respect to the underlying host-parasite interactions. To this end, we created large numbers of random interaction matrices, analysed the resulting dynamics through simulation, and ascertained whether recombination was favoured or disfavoured. We observed Red Queen dynamics in many of our simulations provided the interaction matrices exhibited sufficient ‘antagonicity’. In agreement with previous studies, strong selection on either hosts or parasites favours selection for increased recombination. However, fast changes in the sign of linkage disequilibrium or epistasis were only infrequently observed and do not appear to be a necessary condition for the Red Queen hypothesis to work. Indeed, recombination was often favoured even though the linkage disequilibrium remained of constant sign throughout the simulations. We conclude that Red Queen-type dynamics involving persistent fluctuations in host and parasite genotype frequencies appear to not be an artefact of specific assumptions about host-parasite fitness interactions, but emerge readily with the general interactions studied here. Our results also indicate that although recombination is often favoured, some of the factors previously thought to be important in this process such as linkage disequilibrium fluctuations need to be reassessed when fitness interactions between hosts and parasites are complex.     

Major histocompatibility complex variability in the clonal Amazon molly,Poecilia formosa: is copy number less important than genotype?

The evolution of sex is still a major unsolved puzzle in biology. One of the most promising theoretical models to answer this question is the Red Queen hypothesis. The Red Queen hypothesis proposes a fast adaptation of pathogens to common genotypes and therefore a negative frequency-dependent selection against common genotypes. Clonal organisms should be especially endangered when co-occurring with closely related sexual species. In this context, major histocompatibility (MHC) genes have been discussed to be auspicious candidates that could provide the genetic basis on which selection for immune competence could act. In this study, we investigated MHC variability in a clonal teleost fish: the Amazon molly, Poecilia formosa . The Amazon molly is an ideal candidate to test the Red Queen hypothesis as it is a clonal species but co-occurs with a closely related sexual species and should therefore be especially susceptible to pathogen infection. We found that allele numbers did in general not differ between sexual and clonal 'species' but that genotypic variability is reduced in the clonally reproducing fish, especially in the polyploids. We conclude that in clonal organisms, genotype frequency might be more important for immune competence than MHC allele number. Amazon mollies and their co-occurring parental species clearly fulfil a prerequisite of the Red Queen hypothesis and should therefore provide an ideal system to experimentally test this basic principle probably underlying the evolution of sex.     

Selective pressures on MHC class II genes in the guppy (Poecilia reticulata) as inferred by hierarchical analysis of population structure

Genes of the major histocompatibility complex, which are the most polymorphic of all vertebrate genes, are a pre‐eminent system for the study of selective pressures that arise from host–pathogen interactions. Balancing selection capable of maintaining high polymorphism should lead to the homogenization of MHC allele frequencies among populations, but there is some evidence to suggest that diversifying selection also operates on the MHC. However, the pattern of population structure observed at MHC loci is likely to depend on the spatial and/or temporal scale examined. Here, we investigated selection acting on MHC genes at different geographic scales using Venezuelan guppy populations inhabiting four regions. We found a significant correlation between MHC and microsatellite allelic richness across populations, which suggests the role of genetic drift in shaping MHC diversity. However, compared to microsatellites, more MHC variation was explained by differences between populations within larger geographic regions and less by the differences between the regions. Furthermore, among proximate populations, variation in MHC allele frequencies was significantly higher compared to microsatellites, indicating that selection acting on MHC may increase population structure at small spatial scales. However, in populations that have significantly diverged at neutral markers, the population‐genetic signature of diversifying selection may be eradicated in the long term by that of balancing selection, which acts to preserve rare alleles and thus maintain a common pool of MHC alleles.     

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.     

MHC polymorphism under host-pathogen coevolution

The genes encoding major histocompatibility (MHC) molecules are among the most polymorphic genes known for vertebrates. Since MHC molecules play an important role in the induction of immune responses, the evolution of MHC polymorphism is often explained in terms of increased protection of hosts against pathogens. Two selective pressures that are thought to be involved are (1) selection favoring MHC heterozygous hosts, and (2) selection for rare MHC alleles by host-pathogen coevolution. We have developed a computer simulation of coevolving hosts and pathogens to study the relative impact of these two mechanisms on the evolution of MHC polymorphism. We found that heterozygote advantage per se is insufficient to explain the high degree of polymorphism at the MHC, even in very large host populations. Host-pathogen coevolution, on the other hand, can easily account for realistic polymorphisms of more than 50 alleles per MHC locus. Since evolving pathogens mainly evade presentation by the most common MHC alleles in the host population, they provide a selective pressure for a large variety of rare MHC alleles. Provided that the host population is sufficiently large, a large set of MHC alleles can persist over many host generations under host-pathogen coevolution, despite the fact that allele frequencies continuously change.Electronic Supplementary Material Supplementary material is available in the online version of this article at     

Density-related changes in selection pattern for major histocompatibility complex genes in fluctuating populations of voles

Host-pathogen interactions are of particular interest in studies of the interplay between population dynamics and natural selection. The major histocompatibility complex (MHC) genes of demographically fluctuating species are highly suitable markers for such studies, because they are involved in initiating the immune response against pathogens and display a high level of adaptive genetic variation. We investigated whether two MHC class II genes (DQA1, DRB) were subjected to contemporary selection during increases in the density of fossorial water vole (Arvicola terrestris) populations, by comparing the neutral genetic structure of seven populations with that estimated from MHC genes. Tests for heterozygosity excess indicated that DQA1 was subject to intense balancing selection. No such selection operated on neutral markers. This pattern of selection became more marked with increasing abundance. In the low-abundance phase, when populations were geographically isolated, both overall differentiation and isolation-by-distance were more marked for MHC genes than for neutral markers. Model-based simulations identified DQA1 as an outlier (i.e. under selection) in a single population, suggesting the action of local selection in fragmented populations. The differences between MHC and neutral markers gradually disappeared with increasing effective migration between sites. In the high-abundance year, DQA1 displayed significantly lower levels of overall differentiation than the neutral markers. This gene therefore displayed stronger homogenization than observed under drift and migration alone. The observed signs of selection were much weaker for DRB. Spatial and temporal fluctuations in parasite pressure and locus-specific selection are probably the most plausible mechanisms underlying the observed changes in selection pattern during the demographic cycle.     

Consequences of population structure on genes under balancing selection

This paper describes a new approach to modeling population structure for genes under strong balancing selection of the type seen in plant self-incompatibility systems and the major histocompatibility complex (MHC) system of vertebrates. Simple analytic solutions for the number of alleles maintained at equilibrium and the expected proportion of alleles shared between demes at various levels are derived and checked against simulation results. The theory accurately captures the dynamics of allele number in a subdivided population and identifies important values of m (migration rate) at which allele number and distribution change qualitatively. Starting from a panmictic population, as migration among demes decreases a qualitative change in dynamics is seen at approximately m(crit) approximately equal to the square root of(s/4piNT) where NT is the total population size and s is a measure of the strength of selection. At this point, demes can no longer maintain their panmictic allele number, due to increasing isolation from the total population. Another qualitative change occurs at a migration rate on the same order of magnitude as the mutation rate, mu. At this point, the demes are highly differentiated for allele complement, and the total number of alleles in the population is increased. Because in general u < m<(crit) at intermediate migration rates slightly fewer alleles may be maintained in the total population than are maintained at panmixia. Within this range, total allele number may not be the best indicator of whether a population is effectively panmictic, and some caution should be used when interpreting samples from such populations. The theory presented here can help to analyze data from genes under balancing selection in subdivided populations.     

The consequences of fluctuating selection for isozyme polymorphisms in Daphnia

Temporal sequences of allele frequencies in natural populations of Daphnia are analyzed to obtain the mean and variance of the selection coefficient for both asexual and sexual phases. In general, the alleles at enzyme loci appear to be quasi-neutral. Although significant variation exists for the estimated selection coefficients, the means are in all cases close to zero. Estimates of the variance of selection intensity are applied to existing models to demonstrate the implications of fluctuating selection for the spatial and temporal distribution of gene frequencies in Daphnia. The empirical and analytical results are shown to provide a possible solution to some previously puzzling aspects of Daphnia population genetic surveys. Neither genetic drift nor diversifying selection are necessary conditions for the local diversification of gene frequencies.     

A comparison of variability and population structure for major histocompatibility complex and microsatellite loci in California coastal steelhead (Oncorhynchus mykiss Walbaum)

The major histocompatibility complex (MHC) contains genes integral to immune response in vertebrates. MHC genes have been shown to be under selection in a number of vertebrate taxa, making them intriguing for population genetic studies. We have conducted a survey of genetic variation in an MHC class II gene for steelhead trout from 24 sites in coastal California and compared this variation to that observed at 16 presumably neutral microsatellite loci. A high amount of allelic variation was observed at the MHC when compared to previously published studies on other Pacific salmonids. Elevated nonsynonymous substitutions, relative to synonymous substitutions, were detected at the MHC gene, indicating the signature of historical balancing selection. The MHC data were tested for correlations to and deviations from the patterns found with the microsatellite data. Estimates of allelic richness for the MHC gene and for the microsatellites were positively correlated, as were estimates of population differentiation (F(ST)). An analysis for F(ST) outliers indicates that the MHC locus has an elevated F(ST) relative to the neutral expectation, although a significant result was found for only one particular geographical subgroup. Relatively uniform allele frequency distributions were detected in four populations, although this finding may be partially due to recent population bottlenecks. These results indicate that, at the scale studied here, drift and migration play a major role in the observed geographical variability of MHC genes in steelhead, and that contemporary selection is relatively weak and difficult to detect.     

Microsatellite null alleles and estimation of population differentiation

Microsatellite null alleles are commonly encountered in population genetics studies, yet little is known about their impact on the estimation of population differentiation. Computer simulations based on the coalescent were used to investigate the evolutionary dynamics of null alleles, their impact on F(ST) and genetic distances, and the efficiency of estimators of null allele frequency. Further, we explored how the existing method for correcting genotype data for null alleles performed in estimating F(ST) and genetic distances, and we compared this method with a new method proposed here (for F(ST) only). Null alleles were likely to be encountered in populations with a large effective size, with an unusually high mutation rate in the flanking regions, and that have diverged from the population from which the cloned allele state was drawn and the primers designed. When populations were significantly differentiated, F(ST) and genetic distances were overestimated in the presence of null alleles. Frequency of null alleles was estimated precisely with the algorithm presented in Dempster et al. (1977). The conventional method for correcting genotype data for null alleles did not provide an accurate estimate of F(ST) and genetic distances. However, the use of the genetic distance of Cavalli-Sforza and Edwards (1967) corrected by the conventional method gave better estimates than those obtained without correction. F(ST) estimation from corrected genotype frequencies performed well when restricted to visible allele sizes. Both the proposed method and the traditional correction method have been implemented in a program that is available free of charge at http://www.montpellier.inra.fr/URLB/. We used 2 published microsatellite data sets based on original and redesigned pairs of primers to empirically confirm our simulation results.     

Drift and selection influence geographic variation at immune loci of prairie-chickens

Previous studies of immunity in wild populations have focused primarily on genes of the major histocompatibility complex (MHC); however, studies of model species have identified additional immune-related genes that also affect fitness. In this study, we sequenced five non-MHC immune genes in six greater prairie-chicken (Tympanuchus cupido) populations that have experienced varying degrees of genetic drift as a consequence of population bottlenecks and fragmentation. We compared patterns of geographic variation at the immune genes with six neutral microsatellite markers to investigate the relative effects of selection and genetic drift. Global F(ST) outlier tests identified positive selection on just one of five immune genes (IAP-1) in one population. In contrast, at other immune genes, standardized G'(ST) values were lower than those at microsatellites for a majority of pairwise population comparisons, consistent with balancing selection or with species-wide positive or purifying selection resulting in similar haplotype frequencies across populations. The effects of genetic drift were also evident as summary statistics (e.g., Tajima's D) did not differ from neutrality for the majority of cases, and immune gene diversity (number of haplotypes per gene) was correlated positively with population size. In summary, we found that both genetic drift and selection shaped variation at the five immune genes, and the strength and type of selection varied among genes. Our results caution that neutral forces, such as drift, can make it difficult to detect current selection on genes.     

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.     

Sequence diversity of the MHC DRB gene in the Eurasian beaver (Castor fiber)

Major histocompatibility complex (MHC) genes, coding molecules which play an important role in immune response, are the most polymorphic genes known in vertebrates. However, MHC polymorphism in some species is limited. MHC monomorphism at several MHC class I and II loci was previously reported for two neighbouring northern European populations of the Eurasian beaver (Castor fiber) and reduced selection for polymorphism has been hypothesized. Here, we analysed a partial sequence of the second exon of the MHC II DRB locus from seven relict European and Asian beaver populations. We detected 10 unique alleles among 76 beavers analysed. Only a western Siberian population was polymorphic, with four alleles detected in 10 individuals. Each of the remaining populations was fixed for a different allele. Sequences showed considerable divergence, suggesting the long persistence of allelic lineages. A significant excess of nonsynonymous substitutions was detected at the antigen binding sites, indicating that sequence evolution of beaver DRB was driven by positive selection. Current MHC monomorphism in the majority of populations may be the result of the superimposition of the recent bottleneck on pre-existing genetic structure resulting from population subdivision and differential pathogen pressure.     

Parasite transmission among relatives halts Red Queen dynamics

The theory that coevolving hosts and parasites create a fluctuating selective environment for one another (i.e., produce Red Queen dynamics) has deep roots in evolutionary biology; yet empirical evidence for Red Queen dynamics remains scarce. Fluctuating coevolutionary dynamics underpin the Red Queen hypothesis for the evolution of sex, as well as hypotheses explaining the persistence of genetic variation under sexual selection, local parasite adaptation, the evolution of mutation rate, and the evolution of nonrandom mating. Coevolutionary models that exhibit Red Queen dynamics typically assume that hosts and parasites encounter one another randomly. However, if related individuals aggregate into family groups or are clustered spatially, related hosts will be more likely to encounter parasites transmitted by genetically similar individuals. Using a model that incorporates familial parasite transmission, we show that a slight degree of familial parasite transmission is sufficient to halt coevolutionary fluctuations. Our results predict that evidence for Red Queen dynamics, and its evolutionary consequences, are most likely to be found in biological systems in which hosts and parasites mix mainly at random, and are less likely to be found in systems with familial aggregation. This presents a challenge to the Red Queen hypothesis and other hypotheses that depend on coevolutionary cycling.