Contrasting patterns of variation in MHC loci in the Alpine newt

Major histocompatibility complex (MHC) genes are essential in pathogen recognition and triggering an adaptive immune response. Although they are the most polymorphic genes in vertebrates, very little information on MHC variation and patterns of evolution are available for amphibians, a group known to be declining rapidly worldwide. As infectious diseases are invoked in the declines, information on MHC variation should contribute to devising appropriate conservation strategies. In this study, we examined MHC variation in 149 Alpine newts ( Mesotriton alpestris ) from three allopatric population groups in Poland at the northeastern margin of the distribution of this species. The genetic distinctiveness of the population groups has previously been shown by studies of skin graft rejection, allozymes and microsatellites. Two putative expressed MHC II loci with contrasting levels of variation and clear evidence of gene conversion/recombination between them were detected. The Meal-DAB locus is highly polymorphic (37 alleles), and shows evidence of historical positive selection for amino acid replacements and substantial geographical differentiation in allelic richness. On the contrary, the Meal-DBB locus exhibits low polymorphism (three alleles differing by up to two synonymous substitutions) and a uniform distribution of three alleles among geographical regions. The uniform frequencies of the presumptively neutral Meal-DBB alleles may be explained by linkage to Meal-DAB . We found differences in allelic richness in Meal-DAB between regions, consistent with the hypothesis that genetic drift prevails with increasing distance from glacial refugia. Pseudogene loci appear to have evolved neutrally. The level of DAB variation correlated with variation in microsatellite loci, implying that selection and drift interplayed to produce the pattern of MHC variation observed in marginal populations of the Alpine newt.     

Genetic diversity and differentiation at MHC genes in island populations of tuatara (Sphenodon spp.)

Neutral genetic markers are commonly used to understand the effects of fragmentation and population bottlenecks on genetic variation in threatened species. Although neutral markers are useful for inferring population history, the analysis of functional genes is required to determine the significance of any observed geographical differences in variation. The genes of the major histocompatibility complex (MHC) are well‐known examples of genes of adaptive significance and are particularly relevant to conservation because of their role in pathogen resistance. In this study, we survey diversity at MHC class I loci across a range of tuatara populations. We compare the levels of MHC variation with that observed at neutral microsatellite markers to determine the relative roles of balancing selection, diversifying selection and genetic drift in shaping patterns of MHC variation in isolated populations. In general, levels of MHC variation within tuatara populations are concordant with microsatellite variation. Tuatara populations are highly differentiated at MHC genes, particularly between the northern and Cook Strait regions, and a trend towards diversifying selection across populations was observed. However, overall our results indicate that population bottlenecks and isolation have a larger influence on patterns of MHC variation in tuatara populations than selection.     

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.     

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.     

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.     

Interspecific hybridization increases MHC class II diversity in two sister species of newts

Our understanding of the evolutionary mechanisms generating variation within the highly polymorphic major histocompatibility complex (MHC) genes remains incomplete. Assessing MHC variation across multiple populations, of recent and ancient divergence, may facilitate understanding of geographical and temporal aspects of variation. Here, we applied 454 sequencing to perform a large-scale, comprehensive analysis of MHC class II in the closely related, hybridizing newts, Lissotriton vulgaris (Lv) and Lissotriton montandoni (Lm). Our study revealed an extensive (299 alleles) geographically structured polymorphism. Populations at the southern margin of the Lv distribution, inhabited by old and distinct lineages (southern Lv), exhibited moderate MHC variation and strong population structure, indicating little gene flow or extensive local adaptation. Lissotriton vulgaris in central Europe and the northern Balkans (northern Lv) and almost all Lm populations had a high MHC variation. A much higher proportion of MHC alleles was shared between Lm and northern Lv than between Lm and southern Lv. Strikingly, the average pairwise F(ST) between northern Lv and Lm was significantly lower than between northern and southern Lv for MHC, but not for microsatellites. Thus, high MHC variation in Lm and northern Lv may result from gene flow between species. We hypothesize that the interspecific exchange of MHC genes may be facilitated by frequency-dependent selection. A marginally significant correlation between the MHC and microsatellite allelic richness indicates that demographic factors may have contributed to the present-day pattern of MHC variation, but unequivocal signatures of adaptive evolution in MHC class II sequences emphasize the role of selection on a longer timescale.     

Balancing selection,genetic drift,and human‐mediated introgression interplay to shape MHC (functional) diversity in Mediterranean brown trout

The extraordinary polymorphism of major histocompatibility complex (MHC) genes is considered a paradigm of pathogen‐mediated balancing selection, although empirical evidence is still scarce. Furthermore, the relative contribution of balancing selection to shape MHC population structure and diversity, compared to that of neutral forces, as well as its interaction with other evolutionary processes such as hybridization, remains largely unclear. To investigate these issues, we analyzed adaptive (MHC‐DAB gene) and neutral (11 microsatellite loci) variation in 156 brown trout (Salmo trutta complex) from six wild populations in central Italy exposed to introgression from domestic hatchery lineages (assessed with the LDH gene). MHC diversity and structuring correlated with those at microsatellites, indicating the substantial role of neutral forces. However, individuals carrying locally rare MHC alleles/supertypes were in better body condition (a proxy of individual fitness/parasite load) regardless of the zygosity status and degree of sequence dissimilarity of MHC, hence supporting balancing selection under rare allele advantage, but not heterozygote advantage or divergent allele advantage. The association between specific MHC supertypes and body condition confirmed in part this finding. Across populations, MHC allelic richness increased with increasing admixture between native and domestic lineages, indicating introgression as a source of MHC variation. Furthermore, introgression across populations appeared more pronounced for MHC than microsatellites, possibly because initially rare MHC variants are expected to introgress more readily under rare allele advantage. Providing evidence for the complex interplay among neutral evolutionary forces, balancing selection, and human‐mediated introgression in shaping the pattern of MHC (functional) variation, our findings contribute to a deeper understanding of the evolution of MHC genes in wild populations exposed to anthropogenic disturbance.     

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.     

Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes

Microsatellites and mitochondrial DNA (mtDNA) have traditionally been used in population genetics because of their variability and presumed neutrality, whereas genes of the major histocompatibility complex (MHC) are increasingly of interest because strong selective pressures shape their standing variation. Despite the potential for MHC genes, microsatellites, and mtDNA sequences to complement one another in deciphering population history and demography, the three are rarely used in tandem. Here we report on MHC, microsatellite, and mtDNA variability in a single large population of the eastern tiger salamander (Ambystoma tigrinum tigrinum). We use the mtDNA mismatch distribution and, on microsatellite data, the imbalance index and bottleneck tests to infer aspects of population history and demography. Haplotype and allelic variation was high at all loci surveyed, and heterozygosity was high at the nuclear loci. We find concordance among neutral molecular markers that suggests our study population originated from post-Pleistocene expansions of multiple, fragmented sources that shared few migrants. Differences in N(e) estimates derived from haploid and diploid genetic markers are potentially attributable to secondary contact among source populations that experienced rapid mtDNA divergence and comparatively low levels of nuclear DNA divergence. We find strong evidence of natural selection acting on MHC genes and estimate long-term effective population sizes (N(e)) that are very large, making small selection intensities significant evolutionary forces in this population.     

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.     

Variation in MHC genotypes in two populations of house sparrow (Passer domesticus) with different population histories

Small populations are likely to have a low genetic ability for disease resistance due to loss of genetic variation through inbreeding and genetic drift. In vertebrates, the highest genetic diversity of the immune system is located at genes within the major histocompatibility complex (MHC). Interestingly, parasite‐mediated selection is thought to potentially maintain variation at MHC loci even in populations that are monomorphic at other loci. Therefore, general loss of genetic variation in the genome may not necessarily be associated with low variation at MHC loci. We evaluated inter‐ and intrapopulation variation in MHC genotypes between an inbred (Aldra) and a relatively outbred population (Hestmannøy) of house sparrows (Passer domesticus) in a metapopulation at Helgeland, Norway. Genomic (gDNA) and transcribed (cDNA) alleles of functional MHC class I and IIB loci, along with neutral noncoding microsatellite markers, were analyzed to obtain relevant estimates of genetic variation. We found lower allelic richness in microsatellites in the inbred population, but high genetic variation in MHC class I and IIB loci in both populations. This suggests that also the inbred population could be under balancing selection to maintain genetic variation for pathogen resistance.     

Lack of Spatial Immunogenetic Structure among Wolverine (Gulo gulo) Populations Suggestive of Broad Scale Balancing Selection

Elucidating the adaptive genetic potential of wildlife populations to environmental selective pressures is fundamental for species conservation. Genes of the major histocompatibility complex (MHC) are highly polymorphic, and play a key role in the adaptive immune response against pathogens. MHC polymorphism has been linked to balancing selection or heterogeneous selection promoting local adaptation. However, spatial patterns of MHC polymorphism are also influenced by gene flow and drift. Wolverines are highly vagile, inhabiting varied ecoregions that include boreal forest, taiga, tundra, and high alpine ecosystems. Here, we investigated the immunogenetic variation of wolverines in Canada as a surrogate for identifying local adaptation by contrasting the genetic structure at MHC relative to the structure at 11 neutral microsatellites to account for gene flow and drift. Evidence of historical positive selection was detected at MHC using maximum likelihood codon-based methods. Bayesian and multivariate cluster analyses revealed weaker population genetic differentiation at MHC relative to the increasing microsatellite genetic structure towards the eastern wolverine distribution. Mantel correlations of MHC against geographical distances showed no pattern of isolation by distance (IBD: r = -0.03, p = 0.9), whereas for microsatellites we found a relatively strong and significant IBD (r = 0.54, p = 0.01). Moreover, we found a significant correlation between microsatellite allelic richness and the mean number of MHC alleles, but we did not observe low MHC diversity in small populations. Overall these results suggest that MHC polymorphism has been influenced primarily by balancing selection and to a lesser extent by neutral processes such as genetic drift, with no clear evidence for local adaptation. This study contributes to our understanding of how vulnerable populations of wolverines may respond to selective pressures across their range.     

Inference of selection based on temporal genetic differentiation in the study of highly polymorphic multigene families

The co-evolutionary arms race between host immune genes and parasite virulence genes is known as Red Queen dynamics. Temporal fluctuations in allele frequencies, or the 'turnover' of alleles at immune genes, are concordant with predictions of the Red Queen hypothesis. Such observations are often taken as evidence of host-parasite co-evolution. Here, we use computer simulations of the Major Histocompatibility Complex (MHC) of guppies (Poecilia reticulata) to study the turnover rate of alleles (temporal genetic differentiation, G'(ST)). Temporal fluctuations in MHC allele frequencies can be ≥≤order of magnitude larger than changes observed at neutral loci. Although such large fluctuations in the MHC are consistent with Red Queen dynamics, simulations show that other demographic and population genetic processes can account for this observation, these include: (1) overdominant selection, (2) fluctuating population size within a metapopulation, and (3) the number of novel MHC alleles introduced by immigrants when there are multiple duplicated genes. Synergy between these forces combined with migration rate and the effective population size can drive the rapid turnover in MHC alleles. We posit that rapid allelic turnover is an inherent property of highly polymorphic multigene families and that it cannot be taken as evidence of Red Queen dynamics. Furthermore, combining temporal samples in spatial F(ST) outlier analysis may obscure the signal of selection.     

The mode of evolution of molecular markers in populations of house mice under artificial selection for locomotor behavior

A complete understanding of the mode of evolution of molecular markers is important for making inferences about different population genetic parameters, especially because a number of studies have reported patterns of allelic variation at molecular markers that are not in agreement with neutral evolutionary expectations. In the present study, house mice (Mus domesticus) from the fourteenth generation of a selection experiment for increased voluntary wheel-running activity were used to test how selection on a complex behavior affects the distribution of allelic variation by examining patterns of variation at six microsatellite and four allozyme loci. This population had a hierarchical structure that allowed for simultaneous testing of the effects of selection and genetic drift on the distribution of allelic variation by comparing observed patterns of allele frequencies and estimates of genetic divergence at multiple hierarchical levels to expectations under models of neutral evolution. The levels of genetic divergence among replicate lines and between selection groups, estimated from microsatellite data or pooled microsatellite and allozyme data, were not significantly different from expectations under neutral evolution. Furthermore, the pattern of change of allele frequencies between the base population and generation 14 was largely in agreement with expectations under neutral evolution (although the PGM locus exhibited a pattern of change within populations that was difficult to explain under neutral evolution). Overall the results generally provide support for the neutral evolution of molecular markers.     

Low diversity in the major histocompatibility complex class II DRB1 gene of the Spanish ibex, Capra pyrenaica

During the last two centuries, the Spanish ibex (Capra pyrenaica) has shown a significant demographic decline as a result of the progressive destruction of its natural habitat, disease epidemics, and uncontrolled hunting. Partial sequencing of the class II MHC DRB1 gene revealed that the Spanish ibex has remarkably low levels of genetic variation at this locus, with only six different DRB1 alleles and an observed heterozygosity of 0.429-0.579. The rates of nonsynonymous vs synonymous substitutions were significantly different in the peptide-binding region (dN/dS=5.347, P=0.002), a feature that indicates that the DRB1 gene is under positive selection. A phylogenetic analysis of the Spanish ibex and a set of domestic goat DRB1 alleles revealed that the reported sequences represent four major allelic lineages. The limited allelic repertoire of the DRB1 gene in the Spanish ibex is likely the direct result of the recent history of population bottlenecks and marked demographic decline of this species. A genetic survey of 13 microsatellite loci was consistent with this idea. The Spanish ibex subspecies C. p. hispanica and C. p. victoriae consistently showed considerably lower levels of microsatellite heterozygosity (Ho=0.184-0.231) and allelic diversity (mean number of alleles per locus=2-2.4) than those reported in other wild ruminants. This study demonstrates the significance of both natural selection and the demographic history of populations in determining patterns of genetic variation at MHC loci. In addition, our results emphasize the importance of locally adapted populations for the preservation of genetic diversity.     

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.     

Molecular and quantitative genetic differentiation across Europe in yellow dung flies

Relating geographic variation in quantitative traits to underlying population structure is crucial for understanding processes driving population differentiation, isolation and ultimately speciation. Our study represents a comprehensive population genetic survey of the yellow dung fly Scathophaga stercoraria, an important model organism for evolutionary and ecological studies, over a broad geographic scale across Europe (10 populations from the Swiss Alps to Iceland). We simultaneously assessed differentiation in five quantitative traits (body size, development time, growth rate, proportion of diapausing individuals and duration of diapause), to compare differentiation in neutral marker loci (F(ST)) to that of quantitative traits (Q(ST)). Despite long distances and uninhabitable areas between sampled populations, population structuring was very low but significant (F(ST) = 0.007, 13 microsatellite markers; F(ST) = 0.012, three allozyme markers; F(ST) = 0.007, markers combined). However, only two populations (Iceland and Sweden) showed significant allelic differentiation to all other populations. We estimated high levels of gene flow [effective number of migrants (Nm) = 6.2], there was no isolation by distance, and no indication of past genetic bottlenecks (i.e. founder events) and associated loss of genetic diversity in any northern or island population. In contrast to the low population structure, quantitative traits were strongly genetically differentiated among populations, following latitudinal clines, suggesting that selection is responsible for life history differentiation in yellow dung flies across Europe.     

Wild cyclic voles maintain high neutral and MHC diversity without strong evidence for parasite-mediated selection

The major histocompatibility complex (MHC) is an important component of vertebrate immune defense involved with self/nonself recognition and disease susceptibility. The high variability of genes of the MHC is thought to arise from both parasite-mediated and sexual selection. An outstanding question involves the degree to which balancing selection can oppose genetic drift to maintain high MHC diversity in the face of population bottlenecks. To address this question we examined genetic diversity and population structure at neutral (microsatellite) and MHC genes in montane voles [Microtus montanus (Peale, 1848)] subject to high amplitude population fluctuations, and compared these to measures of infection by common gastrointestinal parasites. We found high neutral and MHC allelic variability, indicating low impacts of genetic drift despite large fluctuations in population size. Greater MHC diversity did not predict lower parasite richness or infection by the two most common endoparasites (cestodes and coccidian protozoa), as might be expected if genotypic composition confers resistance to infection. One specific MHC allele predicted lower cestode intensity, but we found no other associations between MHC and infection measures. Neutral heterozygosity was positively associated with total parasite richness, possibly owing to greater parasite tolerance among heterozygous relative to more inbred hosts. Overall, these results suggest that factors beyond the parasites examined here, such as high inter-patch migration, mate choice, gene conversion or other infectious agents, are likely maintaining the high levels of MHC diversity observed in wild montane voles.     

Contrasting epidemic histories reveal pathogen-mediated balancing selection on class II MHC diversity in a wild songbird

The extent to which pathogens maintain the extraordinary polymorphism at vertebrate Major Histocompatibility Complex (MHC) genes via balancing selection has intrigued evolutionary biologists for over half a century, but direct tests remain challenging. Here we examine whether a well-characterized epidemic of Mycoplasmal conjunctivitis resulted in balancing selection on class II MHC in a wild songbird host, the house finch (Carpodacus mexicanus). First, we confirmed the potential for pathogen-mediated balancing selection by experimentally demonstrating that house finches with intermediate to high multi-locus MHC diversity are more resistant to challenge with Mycoplasma gallisepticum. Second, we documented sequence and diversity-based signatures of pathogen-mediated balancing selection at class II MHC in exposed host populations that were absent in unexposed, control populations across an equivalent time period. Multi-locus MHC diversity significantly increased in exposed host populations following the epidemic despite initial compromised diversity levels from a recent introduction bottleneck in the exposed host range. We did not observe equivalent changes in allelic diversity or heterozygosity across eight neutral microsatellite loci, suggesting that the observations reflect selection rather than neutral demographic processes. Our results indicate that a virulent pathogen can exert sufficient balancing selection on class II MHC to rescue compromised levels of genetic variation for host resistance in a recently bottlenecked population. These results provide evidence for Haldane's long-standing hypothesis that pathogens directly contribute to the maintenance of the tremendous levels of genetic variation detected in natural populations of vertebrates.