Heterozygote advantage fails to explain the high degree of polymorphism of the MHC

Major histocompatibility (MHC) molecules are encoded by extremely polymorphic genes and play a crucial role in vertebrate immunity. Natural selection favors MHC heterozygous hosts because individuals heterozygous at the MHC can present a larger diversity of peptides from infectious pathogens than homozygous individuals. Whether or not heterozygote advantage is sufficient to account for a high degree of polymorphism is controversial, however. Using mathematical models we studied the degree of MHC polymorphism arising when heterozygote advantage is the only selection pressure. We argue that existing models are misleading in that the fitness of heterozygotes is not related to the MHC alleles they harbor. To correct for this, we have developed novel models in which the genotypic fitness of a host directly reflects the fitness contributions of its MHC alleles. The mathematical analysis suggests that a high degree of polymorphism can only be accounted for if the different MHC alleles confer unrealistically similar fitnesses. This conclusion was confirmed by stochastic simulations, including mutation, genetic drift, and a finite population size. Heterozygote advantage on its own is insufficient to explain the high population diversity of the MHC.Electronic Supplementary Material Supplementary material is available in the online version of this article at     

MHC diversity in bottlenecked populations: a simulation model

The depletion of variation at MHC loci, which play a crucial role in pathogen recognition, has been postulated to be one of important extinction risk factors for endangered populations. Thus, it is important to understand how selection affects the level of polymorphism in these genes when populations undergo a reduction in size. We followed MHC diversity in computer simulations of population bottlenecks. The fates of MHC alleles in the simulations were determined either by drift, or by balancing selection resulting from host–parasite coevolution. We found that the impact of selection on MHC polymorphism in bottlenecked populations was dependent upon the timescales involved. Initially, selection maintained lower number of alleles than drift, but after ~40 generations of hosts selection maintained higher MHC diversity, as compared to drift. The adverse effects of decreased MHC polymorphism on population viability may be, to some extent, compensated for if selection helps to retain MHC alleles which show high functional diversity, which should allow protection against a broader range of pathogens. Our simulation shows, however, that the mean divergence of alleles retained under selection in bottlenecked populations is not, on average, significantly higher than the divergence due to drift.     

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.     

Isolation over 35 years in a heated biotest basin causes selection on MHC class IIß genes in the European perch (Perca fluviatilis L.)

Genes that play key roles in host immunity such as the major histocompatibility complex (MHC) in vertebrates are expected to be major targets of selection. It is well known that environmental conditions can have an effect on host–parasite interactions and may thus influence the selection on MHC. We analyzed MHC class IIß variability over 35 years in a population of perch (Perca fluviatilis) from the Baltic Sea that was split into two populations separated from each other. One population was subjected to heating from cooling water of a nuclear power plant and was isolated from the surrounding environment in an artificial lake, while the other population was not subjected to any change in water temperature (control). The isolated population experienced a change of the allelic composition and a decrease in allelic richness of MHC genes compared to the control population. The two most common MHC alleles showed cyclic patterns indicating ongoing parasite–host coevolution in both populations, but the alleles that showed a cyclic behavior differed between the two populations. No such patterns were observed at alleles from nine microsatellite loci, and no genetic differentiation was found between populations. We found no indications for a genetic bottleneck in the isolated population during the 35 years. Additionally, differences in parasitism of the current perch populations suggest that a change of the parasite communities has occurred over the isolation period, although the evidence in form of in‐depth knowledge of the change of the parasite community over time is lacking. Our results are consistent with the hypothesis of a selective sweep imposed by a change in the parasite community.     

Evidence for evolutionary convergence at MHC in two broadly distributed mesocarnivores

Variation within major histocompatibility complex (MHC) genes is important in recognizing pathogens and initiating an immune response. These genes are relevant in enhancing our understanding of how species cope with rapid environmental changes and concomitant fluctuations in selective pressures such as invasive, infectious diseases. Disease-based models suggest that diversity at MHC is maintained through balancing selection arising from the coevolution of hosts and pathogens. Despite intensive balancing selection, sequence motifs or even identical MHC alleles can be shared across multiple species; three potential mechanisms have been put forth to explain this phenomenon: common ancestry, convergent evolution, and random chance. To understand the processes that maintain MHC similarity across divergent species, we examined the variation at two orthologous MHC-DRB genes in widespread North American Musteloid species, striped skunks (Mephitis mephitis), and raccoons (Procyon lotor). These species are often sympatric and exposed to a similar suite of diseases (e.g., rabies, canine distemper, and parvovirus). Given their exposure to similar selective pressures from pathogens, we postulated that similar DRB alleles may be present in both species. Our results indicated that similar motifs are present within both species, at functionally relevant polymorphic sites. However, based on phylogenetic analyses that included previously published MHC sequences of several closely related carnivores, the respective MHC-DRB alleles do not appear to have been maintained through common ancestry and unlikely through random chance. Instead, the similarities observed between the two mesocarnivore species may rather be due to evolutionary convergence.     

Parasite load and MHC diversity in undisturbed and agriculturally modified habitats of the ornate dragon lizard

Major histocompatibility complex (MHC) gene polymorphism is thought to be driven by host–parasite co‐evolution, but the evidence for an association between the selective pressure from parasites and the number of MHC alleles segregating in a population is scarce and inconsistent. Here, we characterized MHC class I polymorphism in a lizard whose habitat preferences (rock outcrops) lead to the formation of well‐defined and stable populations. We investigated the association between the load of ticks, which were used as a proxy for the load of pathogens they transmit, and MHC class I polymorphism across populations in two types of habitat: undisturbed reserves and agricultural land. We hypothesized that the association would be positive across undisturbed reserve populations, but across fragmented agricultural land populations, the relationship would be distorted by the loss of MHC variation due to drift. After controlling for habitat, MHC diversity was not associated with tick number, and the habitats did not differ in this respect. Neither did we detect a difference between habitats in the relationship between MHC and neutral diversity, which was positive across all populations. However, there was extensive variation in the number of MHC alleles per individual, and we found that tick number was positively associated with the average number of alleles carried by lizards across reserve populations, but not across populations from disturbed agricultural land. Our results thus indicate that local differences in selection from parasites may contribute to MHC copy number variation within species, but habitat degradation can distort this relationship.     

Neutral and selective processes shape MHC gene diversity and expression in stocked brook charr populations (Salvelinus fontinalis)

The capacity of an individual to battle infection is an important fitness determinant in wild vertebrate populations. The major histocompatibility complex (MHC) genes are crucial for a host's adaptive immune system to detect pathogens. However, anthropogenic activities may disrupt natural cycles of co‐evolution between hosts and pathogens. In this study, we investigated the dynamic sequence and expression variation of host parasite interactions in brook charr (Salvelinus fontinalis) in a context of past human disturbance via population supplementation from domestic individuals. To do so, we developed a new method to examine selection shaping MHC diversity within and between populations and found a complex interplay between neutral and selective processes that varied between lakes that were investigated. We provided evidence for a lower introgression rate of domestic alleles and found that parasite infection increased with domestic genomic background of individuals. We also documented an association between individual MHC alleles and parasite taxa. Finally, longer cis‐regulatory minisatellites were positively correlated with MHC II down‐regulation and domestic admixture, suggesting that inadvertent selection during domestication resulted in a lower immune response capacity, through a trade‐off between growth and immunity, which explained the negative selection of domestic alleles at least under certain circumstances.     

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.     

Two-step infection processes can lead to coevolution between functionally independent infection and resistance pathways

There is growing evidence that successful infection of hosts by pathogens requires a series of independent steps. However, how multistep infection processes affect host-pathogen coevolution is unclear. We present a coevolutionary model, inspired by empirical observations from a range of host-pathogen systems, where the infection process consists of the following two steps: the first is for the pathogen to recognize and locate a suitable host, and the second is to exploit the host while evading immunity. Importantly, these two steps conform to different models of infection genetics: inverse-gene-for-gene (IGFG) and gene-for-gene (GFG), respectively. We show that coevolution under this scenario can lead to coupled gene frequency changes across these two systems. In particular, selection often favors pathogens that are infective at the first, IGFG, step and hosts that are resistant at the second, GFG, step. Hence, there may be signals of positive selection between functionally independent systems whenever there are multistep processes determining resistance and infectivity. Such multistep infection processes are a fundamental, but overlooked feature of many host-pathogen interactions, and have important consequences for our understanding of host-pathogen coevolution.     

Bacteria-phage coevolution and the emergence of generalist pathogens

Understanding the genetic constraints on pathogen evolution will help to predict the emergence of generalist pathogens that can infect a range of different host genotypes. Here we show that generalist viral pathogens are more likely to emerge during coevolution between the bacterium Pseudomonas fluorescens and the lytic phage SBW25Φ2 than when the same pathogen is challenged to adapt to a nonevolving population of novel hosts. When phages were able to adapt to nonevolving novel hosts, the resulting phenotypes had relatively narrow host ranges compared with coevolved phages. Evolved (rather than coevolved) phages also had lower virulence, although they attained virulence similar to that of coevolved phages after continued adaptation to a nonevolving population of the same host. We explain these results by using sequence data showing that the evolution of broad host range is associated with several different amino acid substitutions and therefore occurs only through repeated rounds of selection for novel infectivity alleles. These findings suggest that generalist bacteriophages are more likely to emerge through long-term coevolution with their hosts than through spontaneous adaptation to a single novel host. These results are likely to be relevant to host-parasite systems where parasite generalism can evolve through the acquisition of multiple mutations or alleles, as appears to be the case for many plant-bacteria and bacteria-virus interactions.     

Long-Term and Short-Term Evolutionary Impacts of Transposable Elements on Drosophila

Transposable elements (TEs) are considered to be genomic parasites and their interactions with their hosts have been likened to the coevolution between host and other nongenomic, horizontally transferred pathogens. TE families, however, are vertically inherited as integral segments of the nuclear genome. This transmission strategy has been suggested to weaken the selective benefits of host alleles repressing the transposition of specific TE variants. On the other hand, the elevated rates of TE transposition and high incidences of deleterious mutations observed during the rare cases of horizontal transfers of TE families between species could create at least a transient process analogous to the influence of horizontally transmitted pathogens. Here, we formally address this analogy, using empirical and theoretical analysis to specify the mechanism of how host–TE interactions may drive the evolution of host genes. We found that host TE-interacting genes actually have more pervasive evidence of adaptive evolution than immunity genes that interact with nongenomic pathogens in Drosophila. Yet, both our theoretical modeling and empirical observations comparing Drosophila melanogaster populations before and after the horizontal transfer of P elements, which invaded D. melanogaster early last century, demonstrated that horizontally transferred TEs have only a limited influence on host TE-interacting genes. We propose that the more prevalent and constant interaction with multiple vertically transmitted TE families may instead be the main force driving the fast evolution of TE-interacting genes, which is fundamentally different from the gene-for-gene interaction of host–pathogen coevolution.     

Trans-species polymorphism and evidence of selection on class II MHC loci in tuco-tucos (Rodentia: Ctenomyidae)

Balancing selection acting over the evolutionary history of a lineage can result in the retention of alleles among species for longer than expected under neutral evolution. The associated pattern of trans-species polymorphism, in which similar or even identical alleles are shared among species, is often used to infer that balancing selection has occurred. The genes of the major histocompatibility complex (MHC) are thought to be subject to balancing selection that maintains alleles associated with response to specific pathogens. To explore the role of balancing selection in shaping MHC diversity in ctenomyid rodents, we examined allelic variability at the class II DRB and DQA loci in 18 species in the genus Ctenomys. Previous studies of four of these species had revealed significant within-population evidence of positive selection on MHC loci. The current study expands upon these analyses to (1) evaluate among-species evidence of positive selection and (2) explore the potential for balancing selection on MHC genes. Interspecific nucleotide sequence variation revealed significant evidence of positive selection on the DRB and DQA loci. At the same time, comparisons of phylogenetic trees for these MHC loci with a putative species tree based on mitochondrial sequence data revealed multiple examples of trans-specific polymorphism, including sharing of identical DRB and DQA alleles among distantly related species of Ctenomys. These findings suggest that MHC genes in these animals have historically been subject to balancing selection and yield new insights into the complex suite of forces shaping MHC diversity in free-living vertebrates.     

Major histocompatibility complex variation and age‐specific endoparasite load in subadult European rabbits

Genes of the major histocompatibility complex (MHC) play a fundamental role in the vertebrate immune response and are amongst the most polymorphic genes in vertebrate genomes. It is generally agreed that the highly polymorphic nature of the MHC is maintained through host–parasite co‐evolution. Two nonexclusive mechanisms of selection are supposed to act on MHC genes: superiority of MHC heterozygous individuals (overdominance) and an advantage for rare MHC alleles. However, the precise mechanisms and their relative importance are still unknown. Here, we examined MHC dependent parasite load in European rabbits (Oryctolagus cuniculus) from a distinct population with low MHC diversity (three alleles, six genotypes). Using a multivariate approach, we tested for associations of individual MHC class II DRB constitution and the rabbits’ intestinal burden with nematodes and coccidia. Rabbits having a particular allele showed lower infestations with hepatic coccidia (E. stiedai). However, a comparison of all six genotypes in the population revealed that carriers of this allele only benefit when they are heterozygous, and furthermore, MHC heterozygosity in general did not affect individual parasite load. In conclusion, this study suggests an immunogenetic basis of European rabbit resistance to hepatic coccidiosis, which can strongly limit survival to maturity in this species. Our study gives a complex picture of MHC–parasite correlations, unveiling the limits of the classical hypotheses of how MHC polymorphism is maintained in natural systems.     

Sympatric and Allopatric Divergence of MHC Genes in Threespine Stickleback

Parasites can strongly affect the evolution of their hosts, but their effects on host diversification are less clear. In theory, contrasting parasite communities in different foraging habitats could generate divergent selection on hosts and promote ecological speciation. Immune systems are costly to maintain, adaptable, and an important component of individual fitness. As a result, immune system genes, such as those of the Major Histocompatability Complex (MHC), can change rapidly in response to parasite-mediated selection. In threespine stickleback (Gasterosteus aculeatus), as well as in other vertebrates, MHC genes have been linked with female mating preference, suggesting that divergent selection acting on MHC genes might influence speciation. Here, we examined genetic variation at MHC Class II loci of sticklebacks from two lakes with a limnetic and benthic species pair, and two lakes with a single species. In both lakes with species pairs, limnetics and benthics differed in their composition of MHC alleles, and limnetics had fewer MHC alleles per individual than benthics. Similar to the limnetics, the allopatric population with a pelagic phenotype had few MHC alleles per individual, suggesting a correlation between MHC genotype and foraging habitat. Using a simulation model we show that the diversity and composition of MHC alleles in a sympatric species pair depends on the amount of assortative mating and on the strength of parasite-mediated selection in adjacent foraging habitats. Our results indicate parallel divergence in the number of MHC alleles between sympatric stickleback species, possibly resulting from the contrasting parasite communities in littoral and pelagic habitats of lakes.     

Major histocompatibility complex polymorphism: dynamics and consequences of parasite-mediated local adaptation in fishes

Parasitism is a common form of life and represents a strong selective pressure for host organisms. In response to this evolutionary pressure, vertebrates have developed genetically coded defences such as the major histocompatibility complex (MHC). Mechanisms of parasite-mediated selection not only maintain outstanding polymorphism in these genes but have also been proposed to further promote host population divergence and ultimately speciation because it can drive evolution of local adaptation in which MHC genes play a crucial role. This review first highlights the dynamics and complexity of parasite-mediated selection in natural systems, which not only depends on dominating parasite strategies and on the taxonomic diversity of the parasite community but also includes the differences in parasite communities between habitats and niches, creating divergent selection on locally adapted populations. Then the different ways in which MHC genes potentially allow vertebrates to respond to these dynamics and to adapt locally are outlined. Finally, it is proposed that varying selection strength in time and space may lead to variation in the strength of precopulatory reproductive isolation which has evolved to maintain local adaptation.     

Coevolutionary genetics of plants and pathogens

Summary The genetic polymorphism maintained by host-pathogen coevolution is analysed in a multilocus model. The model assumes gene-for-gene interactions of the type commonly observed between host plants and their fungal pathogens. Unstable (epidemic) systems maintain more resistance genes, fewer virulence genes, and less overall genetic diversity than stable (endemic) diseases. The stability of the system depends primarily on demographic parameters, such as the pathogen's intrinsic rate of increase, rather than genetic parameters, such as the costs of resistance and virulence. At equilibrium the model predicts that the number of resistance alleles in each host plant follows a binomial distribution that depends on the cost to the pathogen for carrying virulence alleles. Similarly, the number of virulence alleles in each pathogen spore follows a binomial distribution that depends on one minus the cost to the host for carrying resistance alleles. Data from wild populations match the predicted binomial distributions.     

Pathogen burden,co‐infection and major histocompatibility complex variability in the European badger (Meles meles)

Pathogen‐mediated selection is thought to maintain the extreme diversity in the major histocompatibility complex (MHC) genes, operating through the heterozygote advantage, rare‐allele advantage and fluctuating selection mechanisms. Heterozygote advantage (i.e. recognizing and binding a wider range of antigens than homozygotes) is expected to be more detectable when multiple pathogens are considered simultaneously. Here, we test whether MHC diversity in a wild population of European badgers (Meles meles) is driven by pathogen‐mediated selection. We examined individual prevalence (infected or not), infection intensity and co‐infection of 13 pathogens from a range of taxa and examined their relationships with MHC class I and class II variability. This population has a variable, but relatively low, number of MHC alleles and is infected by a variety of naturally occurring pathogens, making it very suitable for the investigation of MHC–pathogen relationships. We found associations between pathogen infections and specific MHC haplotypes and alleles. Co‐infection status was not correlated with MHC heterozygosity, but there was evidence of heterozygote advantage against individual pathogen infections. This suggests that rare‐allele advantages and/or fluctuating selection, and heterozygote advantage are probably the selective forces shaping MHC diversity in this species. We show stronger evidence for MHC associations with infection intensity than for prevalence and conclude that examining both pathogen prevalence and infection intensity is important. Moreover, examination of a large number and diversity of pathogens, and both MHC class I and II genes (which have different functions), provide an improved understanding of the mechanisms driving MHC diversity.     

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.     

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.     

Unique Pathogen Peptidomes Facilitate Pathogen-Specific Selection and Specialization of MHC Alleles

A key component of pathogen-specific adaptive immunity in vertebrates is the presentation of pathogen-derived antigenic peptides by major histocompatibility complex (MHC) molecules. The excessive polymorphism observed at MHC genes is widely presumed to result from the need to recognize diverse pathogens, a process called pathogen-driven balancing selection. This process assumes that pathogens differ in their peptidomes—the pool of short peptides derived from the pathogen’s proteome—so that different pathogens select for different MHC variants with distinct peptide-binding properties. Here, we tested this assumption in a comprehensive data set of 51.9 Mio peptides, derived from the peptidomes of 36 representative human pathogens. Strikingly, we found that 39.7% of the 630 pairwise comparisons among pathogens yielded not a single shared peptide and only 1.8% of pathogen pairs shared more than 1% of their peptides. Indeed, 98.8% of all peptides were unique to a single pathogen species. Using computational binding prediction to characterize the binding specificities of 321 common human MHC class-I variants, we investigated quantitative differences among MHC variants with regard to binding peptides from distinct pathogens. Our analysis showed signatures of specialization toward specific pathogens especially by MHC variants with narrow peptide-binding repertoires. This supports the hypothesis that such fastidious MHC variants might be maintained in the population because they provide an advantage against particular pathogens. Overall, our results establish a key selection factor for the excessive allelic diversity at MHC genes observed in natural populations and illuminate the evolution of variable peptide-binding repertoires among MHC variants.