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     

The effects of historical fragmentation on major histocompatibility complex class II β and microsatellite variation in the Aegean island reptile,Podarcis erhardii

The major histocompatibility complex (MHC) plays a key role in disease resistance and is the most polymorphic gene region in vertebrates. Although habitat fragmentation is predicted to lead to a loss in MHC variation through drift, the impact of other evolutionary forces may counter this effect. Here we assess the impact of selection, drift, migration, and recombination on MHC class II and microsatellite variability in 14 island populations of the Aegean wall lizard Podarcis erhardii. Lizards were sampled from islands within the Cyclades (Greece) formed by rising sea levels as the last glacial maximum approximately 20,000 before present. Bathymetric data were used to determine the area and age of each island, allowing us to infer the corresponding magnitude and timing of genetic bottlenecks associated with island formation. Both MHC and microsatellite variation were positively associated with island area, supporting the hypothesis that drift governs neutral and adaptive variation in this system. However, MHC but not microsatellite variability declined significantly with island age. This discrepancy is likely due to the fact that microsatellites attain mutation‐drift equilibrium more rapidly than MHC. Although we detected signals of balancing selection, recombination and migration, the effects of these evolutionary processes appeared negligible relative to drift. This study demonstrates how land bridge islands can provide novel insights into the impact of historical fragmentation on genetic diversity as well as help disentangle the effects of different evolutionary forces on neutral and adaptive diversity.     

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.     

Geographic heterogeneity in natural selection on an MHC locus in sockeye salmon

Balancing selection maintains high levels of polymorphism and heterozygosity in genes of the MHC (major histocompatibility complex) of vertebrate organisms, and promotes long evolutionary persistence of individual alleles and strongly differentiated allelic lineages. In this study, genetic variation at the MHC class II DAB-beta1 locus was examined in 31 populations of sockeye salmon (Oncorhynchus nerka) inhabiting the Fraser River drainage of British Columbia, Canada. Twenty-five percent of variation at the locus was partitioned among sockeye populations, as compared with 5% at neutral genetic markers. Geographic heterogeneity of balancing selection was detected among four regions in the Fraser River drainage and among lake systems within regions. High levels of beta1 allelic diversity and heterozygosity, as well as distributions of alleles and allelic lineages that were more even than expected for a neutral locus, indicated the presence of balancing selection in populations throughout much of the interior Fraser drainage. However, proximate populations in the upper Fraser region, and four of six populations from the lower Fraser drainage, exhibited much lower levels of genetic diversity and had beta1 allele frequency distributions in conformance with those expected for a neutral locus, or a locus under directional selection. Pair-wise FST values for beta1 averaged 0.19 and tended to exceed the corresponding values estimated for neutral loci at all levels of population structure, although they were lower among populations experiencing balancing selection than among other populations. The apparent heterogeneity in selection resulted in strong genetic differentiation between geographically proximate populations with and without detectable levels of balancing selection, in stark contrast to observations at neutral loci. The strong partitioning and complex structure of beta1 diversity within and among sockeye populations on a small geographic scale illustrates the value of incorporating adaptive variation into conservation planning for the species.     

Adaptive tolerance to a pathogenic fungus drives major histocompatibility complex evolution in natural amphibian populations

Amphibians have been affected globally by the disease chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), and we are just now beginning to understand how immunogenetic variability contributes to disease susceptibility. Lineages of an expressed major histocompatibility complex (MHC) class II locus involved in acquired immunity are associated with chytridiomycosis susceptibility in controlled laboratory challenge assays. Here, we extend these findings to natural populations that vary both in exposure and response to Bd. We find that MHC alleles and supertypes associated with Bd survival in the field show a molecular signal of positive selection, while those associated with susceptibility do not, supporting the hypothesis that heritable Bd tolerance is rapidly evolving. We compare MHC supertypes to neutral loci to demonstrate where selection versus demography is shaping MHC variability. One population with Bd tolerance in nature shows a significant signal of directional selection for the same allele (allele Q) that was significantly associated with survival in an earlier laboratory study. Our findings indicate that selective pressure for Bd survival drives rapid immunogenetic adaptation in some natural populations, despite differences in environment and demography. Our field-based analysis of immunogenetic variation confirms that natural amphibian populations have the evolutionary potential to adapt to chytridiomycosis.     

On the relative roles of selection and genetic drift in shaping MHC variation

Genes of the major histocompatibility complex (MHC) have provided some of the clearest examples of how natural selection generates discordances between adaptive and neutral variation in natural populations. The type and intensity of selection as well as the strength of genetic drift are believed to be important in shaping the resulting pattern of MHC diversity. However, evaluating the relative contribution of multiple microevolutionary forces is challenging, and empirical studies have reported contrasting results. For instance, balancing selection has been invoked to explain high levels of MHC diversity and low population differentiation in comparison with other nuclear markers. Other studies have shown that genetic drift can sometimes overcome selection and then patterns of genetic variation at adaptive loci cannot be discerned from those occurring at neutral markers. Both empirical and simulated data also indicate that loss of genetic diversity at adaptive loci can occur faster than at neutral loci when selection and population bottlenecks act simultaneously. Diversifying selection, on the other hand, explains accelerated MHC divergence as the result of spatial variation in pathogen‐mediated selective regimes. Because of all these possible scenarios and outcomes, collecting information from as many study systems as possible, is crucial to enhance our understanding about the evolutionary forces driving MHC polymorphism. In this issue, Miller and co‐workers present an illuminating contribution by combining neutral markers (microsatellites) and adaptive MHC class I loci during the investigation of genetic differentiation across island populations of tuatara Sphenodon punctatus. Their study of geographical variation reveals a major role of genetic drift in shaping MHC variation, yet they also discuss some support for diversifying selection.     

Are there Ubiquitous Parasite-driven Major Histocompatibility Complex Selection Mechanisms in Gray Mouse Lemurs?

A major goal of evolutionary biology is to understand how selection drives local adaptation. For example, the major histocompatibility complex (MHC) plays an important role in the immune system, and high levels of MHC variation are thought to be a form of adaptation in natural populations. Individual MHC composition may influence parasite resistance via advantages associated with 1) heterozygosity, because heterozygotes recognize a broader range of different antigens than homozygotes (heterozygote advantage); 2) highly variable amino acid sequences in MHC alleles, allowing individuals to bind a broader spectrum of parasite-derived peptides (divergent-alleles advantage, a mechanistic variant of the heterozygote advantage model); or 3) specific MHC alleles (rare allele advantage or frequency dependent selection). We investigated relationships between gastrointestinal nematode burden and both adaptive immune gene variability (MHC class II DRB) and neutral microsatellites in free-living gray mouse lemurs (Microcebus murinus) native to a dry deciduous forest population in western Madagascar to test these hypotheses. The individual MHC composition was related to parasite infestation. Specific MHC alleles were involved in parasite resistance and the presence of common alleles negatively influenced infestation intensity. We found no support for the heterozygote advantage hypothesis, but we did find support for the divergent-MHC allele advantage hypothesis: Individuals with very divergent MHC alleles carried fewer and less intense nematode infestations than individuals with more similar alleles in the more variable dry deciduous forest population. These results indicate that intestinal parasites are important selection pressures under natural conditions and suggest that different selection mechanisms are not mutually exclusive. In contrast, we detected no association between neutral overall individual genetic diversity (measured via 17 microsatellites) and parasite load. Finally, we investigated the ubiquity of parasite-driven selection mechanisms by comparing our results with a previous study of a mouse lemur population from the climatically different littoral forest in southeastern Madagascar, ca. 500 km away. This revealed that different specific MHC alleles were involved in parasite resistance in the 2 habitats, showing that gene-parasite associations are not consistent between populations.     

Pathogen richness and abundance predict patterns of adaptive major histocompatibility complex variation in insular amphibians

The identification of the factors responsible for genetic variation and differentiation at adaptive loci can provide important insights into the evolutionary process and is crucial for the effective management of threatened species. We studied the impact of environmental viral richness and abundance on functional diversity and differentiation of the MHC class Ia locus in populations of the black‐spotted pond frog (Pelophylax nigromaculatus), an IUCN‐listed species, on 24 land‐bridge islands of the Zhoushan Archipelago and three nearby mainland sites. We found a high proportion of private MHC alleles in mainland and insular populations, corresponding to 32 distinct functional supertypes, and strong positive selection on MHC antigen‐binding sites in all populations. Viral pathogen diversity and abundance were reduced at island sites relative to the mainland, and islands housed distinctive viral communities. Standardized MHC diversity at island sites exceeded that found at neutral microsatellites, and the representation of key functional supertypes was positively correlated with the abundance of specific viruses in the environment (Frog virus 3 and Ambystoma tigrinum virus). These results indicate that pathogen‐driven diversifying selection can play an important role in maintaining functionally important MHC variation following island isolation, highlighting the importance of considering functionally important genetic variation and host–pathogen associations in conservation planning and management.     

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake,the Eastern Massasauga (Sistrurus catenatus)

Genetic diversity is fundamental to maintaining the long‐term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.     

Patterns of selection and allele diversity of class I and class II major histocompatibility loci across the species range of sockeye salmon (Oncorhynchus nerka)

The major histocompatibility complex (MHC), an important component of the vertebrate immune system, provides an important suite of genes to examine the role of genetic diversity at non‐neutral loci for population persistence. We contrasted patterns of diversity at the two classical MHC loci in sockeye salmon (Oncorhynchus nerka), MHC class I (UBA) and MHC class II (DAB), and neutral microsatellite loci across 70 populations spanning the species range from Washington State to Japan. There was no correlation in allelic richness or heterozygosity between MHC loci or between MHC loci and microsatellites. The two unlinked MHC loci may be responding to different selective pressures; the distribution of F_ST values for the two loci was uncorrelated, and evidence for both balancing and directional selection on alleles and lineages of DAB and UBA was observed in populations throughout the species range but rarely on both loci within a population. These results suggest that fluctuating selection has resulted in the divergence of MHC loci in contemporary populations.     

Detecting the signature of selection on immune genes in highly structured populations of wild sheep (Ovis dalli)

The confounding effects of population structure complicate efforts to identify regions of the genome under the influence of selection in natural populations. Here we test for evidence of selection in three genes involved in vertebrate immune function - the major histocompatibility complex (MHC), interferon gamma (IFNG) and natural resistance associated macrophage polymorphism (NRAMP) - in highly structured populations of wild thinhorn sheep (Ovis dalli). We examined patterns of variation at microsatellite loci linked to these gene regions and at the DNA sequence level. Simple Watterson's tests indicated balancing selection at all three gene regions. However, evidence for selection was confounded by population structure, as the Watterson's test statistics from linked markers were not outside of the range of values from unlinked and presumably neutral microsatellites. The translated coding sequences of thinhorn IFNG and NRAMP are fixed and identical to those of domestic sheep (Ovis aries). In contrast, the thinhorn MHC DRB locus shows significant evidence of overdominance through both an excess of nonsynonymous substitution and trans-species polymorphism. The failure to detect balancing selection at microsatellite loci linked to the MHC is likely the result of recombination between the markers and expressed gene regions.     

“Balancing” balancing selection? Assortative mating at the major histocompatibility complex despite molecular signatures of balancing selection

In vertebrate animals, genes of the major histocompatibility complex (MHC) determine the set of pathogens to which an individual's adaptive immune system can respond. MHC genes are extraordinarily polymorphic, often showing elevated nonsynonymous relative to synonymous sequence variation and sharing presumably ancient polymorphisms between lineages. These patterns likely reflect pathogen‐mediated balancing selection, for example, rare‐allele or heterozygote advantage. Such selection is often reinforced by disassortative mating at MHC. We characterized exon 2 of MHC class II, corresponding to the hypervariable peptide‐binding region, in song sparrows (Melospiza melodia). We compared nonsynonymous to synonymous sequence variation in order to identify positively selected sites; assessed evidence for trans‐species polymorphisms indicating ancient balancing selection; and compared MHC similarity of socially mated pairs to expectations under random mating. Six codons showed elevated ratios of nonsynonymous to synonymous variation, consistent with balancing selection, and we characterized several alleles similar to those occurring in at least four other avian families. Despite this evidence for historical balancing selection, mated pairs were significantly more similar at MHC than were randomly generated pairings. Nonrandom mating at MHC thus appears to partially counteract, not reinforce, pathogen‐mediated balancing selection in this system. We suggest that in systems where individual fitness does not increase monotonically with MHC diversity, assortative mating may help to avoid excessive offspring heterozygosity that could otherwise arise from long‐standing balancing selection.     

Balancing selection on MHC class I in wild brown trout Salmo trutta

Evidence is reported for balancing selection acting on variation at major histocompatibility complex (MHC) in wild populations of brown trout Salmo trutta. First, variation at an MHC class I (satr-uba)-linked microsatellite locus (mhc1) is retained in small S. trutta populations isolated above waterfalls although variation is lost at neutral microsatellite markers. Second, populations across several catchments are less differentiated at mhc1 than at neutral markers, as predicted by theory. The population structure of these fish was also elucidated.     

Natural selection on MHC IIβ in parapatric lake and stream stickleback: Balancing,divergent, both or neither?

Major histocompatibility complex (MHC) genes encode proteins that play a central role in vertebrates' adaptive immunity to parasites. MHC loci are among the most polymorphic in vertebrates' genomes, inspiring many studies to identify evolutionary processes driving MHC polymorphism within populations and divergence between populations. Leading hypotheses include balancing selection favouring rare alleles within populations, and spatially divergent selection. These hypotheses do not always produce diagnosably distinct predictions, causing many studies of MHC to yield inconsistent or ambiguous results. We suggest a novel strategy to distinguish balancing vs. divergent selection on MHC, taking advantage of natural admixture between parapatric populations. With divergent selection, individuals with immigrant alleles will be more infected and less fit because they are susceptible to novel parasites in their new habitat. With balancing selection, individuals with locally rare immigrant alleles will be more fit (less infected). We tested these contrasting predictions using three‐spine stickleback from three replicate pairs of parapatric lake and stream habitats. We found numerous positive and negative associations between particular MHC IIβ alleles and particular parasite taxa. A few allele–parasite comparisons supported balancing selection, and others supported divergent selection between habitats. But, there was no overall tendency for fish with immigrant MHC alleles to be more or less heavily infected. Instead, locally rare MHC alleles (not necessarily immigrants) were associated with heavier infections. Our results illustrate the complex relationship between MHC IIβ allelic variation and spatially varying multispecies parasite communities: different hypotheses may be concurrently true for different allele–parasite combinations.     

Comparative analysis of population structure across environments and geographical scales at major histocompatibility complex and microsatellite loci in Atlantic salmon (Salmo salar)

Evidence of selection acting on major histocompatibility complex (MHC) genes has been illustrated with the analysis of their nucleotide sequences and allele frequency distribution. Comparing the patterns of population differentiation at neutral markers and MHC genes in the wild may provide further insights about the relative role of selection and neutrality in shaping their diversity. In this study, we combine both methods to assess the role of selection on a MHC gene in Atlantic salmon. We compare variation at a MHC class II B locus and microsatellites among 14 samples from seven different rivers and seven subpopulations within a single river system covering a variety of habitats and different geographical scales. We show that diversifying selection is acting on the sites involved in antigen presentation and that balancing selection maintains a high level of polymorphism within populations. Despite important differences in habitat type, the comparison of the population structure at MHC and microsatellites on large geographical scales reveals a correlation between patterns of differentiation, indicating that drift and migration have been more important than selection in shaping population differentiation at the MHC locus. In contrast, strong discrepancies between patterns of population differentiation at the two types of markers provides support for the role of selection in shaping population structure within rivers. Together, these results confirm that natural selection is influencing MHC gene diversity in wild Atlantic salmon although neutral forces may also be important in their evolution.     

Selection from parasites favours immunogenetic diversity but not divergence among locally adapted host populations

The unprecedented polymorphism in the major histocompatibility complex (MHC) genes is thought to be maintained by balancing selection from parasites. However, do parasites also drive divergence at MHC loci between host populations, or do the effects of balancing selection maintain similarities among populations? We examined MHC variation in populations of the livebearing fish Poecilia mexicana and characterized their parasite communities. Poecilia mexicana populations in the Cueva del Azufre system are locally adapted to darkness and the presence of toxic hydrogen sulphide, representing highly divergent ecotypes or incipient species. Parasite communities differed significantly across populations, and populations with higher parasite loads had higher levels of diversity at class II MHC genes. However, despite different parasite communities, marked divergence in adaptive traits and in neutral genetic markers, we found MHC alleles to be remarkably similar among host populations. Our findings indicate that balancing selection from parasites maintains immunogenetic diversity of hosts, but this process does not promote MHC divergence in this system. On the contrary, we suggest that balancing selection on immunogenetic loci may outweigh divergent selection causing divergence, thereby hindering host divergence and speciation. Our findings support the hypothesis that balancing selection maintains MHC similarities among lineages during and after speciation (trans‐species evolution).     

Population fragmentation and major histocompatibility complex variation in the spotted suslik,Spermophilus suslicus

The fragmentation of populations typically enhances depletion of genetic variation, but highly polymorphic major histocompatibility complex (MHC) genes are thought to be under balancing selection and therefore retain polymorphism despite population bottlenecks. In this study, we investigate MHC DRB (class II) exon 2 variation in 14 spotted suslik populations from two regions differing in their degree of habitat fragmentation and gene flow. We found 16 alleles that segregated in a sample of 248 individuals. The alleles were highly divergent and revealed the hallmark signs of positive selection acting on them in the past, showing a significant excess of nonsynonymous substitutions. This excess was concentrated in putative antigen‐binding sites, which suggests that past selection was driven by pathogens. MHC diversity was significantly lower in fragmented western populations than in the eastern populations, characterized by significant gene flow. In contrast to neutral variation, amova did not reveal genetic differentiation between the two regions. This may indicate similar selective pressures shaping MHC variation in both regions until the recent past. However, MHC allelic richness within a population was correlated with that for microsatellites. F_ST outlier analyses have shown that population differentiation at DRB was neither higher nor lower than expected under neutrality. The results suggest that selection on MHC is not strong enough to counteract drift that results from recent fragmentation of spotted suslik populations.     

TEMPORAL VARIATION AT THE MHC CLASS IIB IN WILD POPULATIONS OF THE GUPPY (POECILIA RETICULATA)

Understanding genetic diversity in natural populations is a fundamental objective of evolutionary biology. The immune genes of the major histocompatibility complex (MHC) are excellent candidates to study such diversity because they are highly polymorphic in populations. Although balancing selection may be responsible for maintaining diversity at these functionally important loci, temporal variation in selection pressure has rarely been examined. We examine temporal variation in MHC class IIB diversity in nine guppy (Poecilia reticulata) populations over two years. We found that five of the populations changed significantly more at the MHC than at neutral (microsatellite) loci as measured by F_ST, which suggests that the change at the MHC was due to selection and not neutral processes. Additionally, pairwise population differentiation measures at the MHC were higher in 2007 than in 2006, with the signature of selection changing from homogenizing to diversifying selection or neutral evolution. Interestingly, within the populations the magnitude of the change at the MHC between years was related to the change in the proportion of individuals infected by a common parasite, indicating a link between genetic structure and the parasite. Our data thereby implicate temporal variation in selective pressure as an important mechanism maintaining diversity at the MHC in wild populations.     

Signals of major histocompatibility complex overdominance in a wild salmonid population

The major histocompatibility complex (MHC) contains the most variable genes in vertebrates, but despite extensive research, the mechanisms maintaining this polymorphism are still unresolved. One hypothesis is that MHC polymorphism is a result of balancing selection operating by overdominance, but convincing evidence for overdominant selection in natural populations has been lacking. We present strong evidence consistent with MHC-specific overdominance in a free-living population of Arctic charr (Salvelinus alpinus) in northernmost Europe. In this population, where just two MHC alleles were observed, MHC heterozygous fish had a lower parasite load, were in better condition (as estimated by a fatness indicator) and had higher survival under stress than either of the homozygotes. Conversely, there was no consistent association between these fitness measures and assumedly neutral microsatellite variability, indicating an MHC-specific effect. Our results provide convincing empirical evidence consistent with the notion that overdominance can be an important evolutionary mechanism contributing to MHC polymorphism in wild animal populations. They also support a recent simulation study indicating that the number of alleles expected to be maintained at an MHC loci can be low, even under strong heterozygote advantage.     

Major histocompatibility complex variation in insular populations of the Egyptian vulture: inferences about the roles of genetic drift and selection

Insular populations have attracted the attention of evolutionary biologists because of their morphological and ecological peculiarities with respect to their mainland counterparts. Founder effects and genetic drift are known to distribute neutral genetic variability in these demes. However, elucidating whether these evolutionary forces have also shaped adaptive variation is crucial to evaluate the real impact of reduced genetic variation in small populations. Genes of the major histocompatibility complex (MHC) are classical examples of evolutionarily relevant loci because of their well-known role in pathogen confrontation and clearance. In this study, we aim to disentangle the partial roles of genetic drift and natural selection in the spatial distribution of MHC variation in insular populations. To this end, we integrate the study of neutral (22 microsatellites and one mtDNA locus) and MHC class II variation in one mainland (Iberia) and two insular populations (Fuerteventura and Menorca) of the endangered Egyptian vulture (Neophron percnopterus). Overall, the distribution of the frequencies of individual MHC alleles (n=17 alleles from two class II B loci) does not significantly depart from neutral expectations, which indicates a prominent role for genetic drift over selection. However, our results point towards an interesting co-evolution of gene duplicates that maintains different pairs of divergent alleles in strong linkage disequilibrium on islands. We hypothesize that the co-evolution of genes may counteract the loss of genetic diversity in insular demes, maximize antigen recognition capabilities when gene diversity is reduced, and promote the co-segregation of the most efficient allele combinations to cope with local pathogen communities.