Genetic diversity and parasite prevalence in two species of bumblebee

Many bumblebee species have been suffering from significant declines across their ranges in the Northern Hemisphere over the last few decades. The remaining populations of the rare species are now often isolated due to habitat fragmentation and have reduced levels of genetic diversity. The persistence of these populations may be threatened by inbreeding depression, which may result in a higher susceptibility to parasites. Here we investigate the relationship between genetic diversity and prevalence of the parasitic mite Locustacarus buchneri in bumblebees, using the previously-studied system of Bombus muscorum and Bombus jonellus in the Western Isles of Scotland. We recorded L. buchneri prevalence in 17 populations of B. muscorum and 13 populations of B. jonellus and related the results to levels of heterozygosity. For B. muscorum, we found that prevalence of the mite was higher in populations with lower genetic diversity but there was no such relationship in the more genetically diverse B. jonellus. In contrast to population-level measures of genetic diversity, the heterozygosity of individual bees was not correlated with infection status. We suggest population-level genetic homogeneity may facilitate parasite transmission and elevate prevalence, with potential consequences for population persistence.     

Genetic structure in insular and mainland populations of house sparrows (Passer domesticus) and their hemosporidian parasites

Small and isolated populations usually exhibit low levels of genetic variability, and thus, they are expected to have a lower capacity to adapt to changes in environmental conditions, such as exposure to pathogens and parasites. Comparing the genetic variability of selectively neutral versus functional loci allows one to assess the evolutionary history of populations and their future evolutionary potential. The genes of the major histocompatibility complex (MHC) control immune recognition of parasites, and their unusually high diversity is genes which is likely driven by parasite‐mediated balancing selection. Here, we examined diversity and differentiation of neutral microsatellite loci and functional MHC class I genes in house sparrows (Passer domesticus), living in six insular and six mainland populations, and we aimed to determine whether their diversity or differentiation correlates with the diversity and the prevalence of infection of hemosporidian parasites. We found that island bird populations tended to have lower neutral genetic variability, whereas MHC variability gene was similar between island and mainland populations. Similarly, island populations tended to show greater genetic differentiation than mainland populations, especially at microsatellite markers. The maintenance of MHC genetic diversity and its less marked structure in the island populations could be attributed to balancing‐selection. The greater MHC differentiation among populations was negatively correlated with similarity in blood parasites (prevalence and diversity of parasite strains) between populations. Even at low prevalence and small geographical scale, haemosporidian parasites might contribute to structure the variability of immune genes among populations of hosts.     

GENETIC DIVERSITY AND CAPILLARIA HEPATICA (NEMATODA) PREVALENCE IN MICHIGAN DEER MOUSE POPULATIONS

There have been few field tests of the hypothesis that homozygous populations are prone to high levels of disease. I tested for a negative correlation between genetic diversity and parasitism by estimating the allozyme heterozygosity, population density, and proportion of individuals infected by Capillaria hepatica (Nematoda) in nine Michigan populations of deer mice (Peromyscus maniculatus). Parasite prevalence was correlated negatively with heterozygosity when the effects of density were held constant, but was not correlated with population density after controlling for the effects of genetic diversity. These data support the prediction that inbred populations will be more susceptible to parasite infestations.     

Comparative genetic diversity of parasites and their hosts: population structure of an urban cockroach and its haplo-diploid parasite (oxyuroid nematode)

Few studies have investigated the genetic structure of both host and parasite populations at a level of populations and at a level of individuals. We investigated the genetic structure of the urban cockroach Blattella germanica and its oxyuroid parasite Blatticola blattae. Random amplified polymorphic DNA (RAPD) markers were used to quantify genetic diversity between and within four populations (from two cities in France) of the host and its parasite. Diversity based on phenotypic frequencies was calculated for each RAPD marker using Shannon-Wiener's index. We used multivariate analyses to test the significance of genetic differentiation between host and parasite populations. Analysis of molecular variance was also used. Both methods gave similar results. Diversity between pairs of individuals was estimated by Nei & Li's index. Genetic diversity was higher within host or parasite populations (80% and 82%, respectively, of explained diversity) than between host or parasite populations (20% and 18%, respectively, explained diversity). The genetic distances between pairs of parasite populations (or individuals) were not correlated with the genetic distances between the corresponding pairs of host populations (or individuals).     

Population structure of a parasitic plant and its perennial host

Characterization of host and parasite population genetic structure and estimation of gene flow among populations are essential for the understanding of parasite local adaptation and coevolutionary interactions between hosts and parasites. We examined two aspects of population structure in a parasitic plant, the greater dodder (Cuscuta europaea) and its host plant, the stinging nettle (Urtica dioica), using allozyme data from 12 host and eight parasite populations. First, we examined whether hosts exposed to parasitism in the past contain higher levels of genetic variation. Second, we examined whether host and parasite populations differ in terms of population structure and if their population structures are correlated. There was no evidence that host populations differed in terms of gene diversity or heterozygosity according to their history of parasitism. Host populations were genetically more differentiated (F(ST) = 0.032) than parasite populations (F(ST) = 0.009). Based on these F(ST) values, gene flow was high for both host and parasite. Such high levels of gene flow could counteract selection for local adaptation of the parasite. We found no significant correlation between geographic and genetic distance (estimated as pairwise F(ST)), either for the host or for the parasite. Furthermore, host and parasite genetic distance matrices were uncorrelated, suggesting that sites with genetically similar host populations are unlikely to have genetically similar parasite populations.     

Multiple infections by the anther smut pathogen are frequent and involve related strains

Population models of host-parasite interactions predict that when different parasite genotypes compete within a host for limited resources, those that exploit the host faster will be selected, leading to an increase in parasite virulence. When parasites sharing a host are related, however, kin selection should lead to more cooperative host exploitation that may involve slower rates of parasite reproduction. Despite their potential importance, studies that assess the prevalence of multiple genotype infections in natural populations remain rare, and studies quantifying the relatedness of parasites occurring together as natural multiple infections are particularly scarce. We investigated multiple infections in natural populations of the systemic fungal plant parasite Microbotryum violaceum, the anther smut of Caryophyllaceae, on its host, Silene latifolia. We found that multiple infections can be extremely frequent, with different fungal genotypes found in different stems of single plants. Multiple infections involved parasite genotypes more closely related than would be expected based upon their genetic diversity or due to spatial substructuring within the parasite populations. Together with previous sequential inoculation experiments, our results suggest that M. violaceum actively excludes divergent competitors while tolerating closely related genotypes. Such an exclusion mechanism might explain why multiple infections were less frequent in populations with the highest genetic diversity, which is at odds with intuitive expectations. Thus, these results demonstrate that genetic diversity can influence the prevalence of multiple infections in nature, which will have important consequences for their optimal levels of virulence. Measuring the occurrence of multiple infections and the relatedness among parasites within hosts in natural populations may be important for understanding the evolutionary dynamics of disease, the consequences of vaccine use, and forces driving the population genetic structure of parasites.     

THE GEOGRAPHY OF COEVOLUTION: COMPARATIVE POPULATION STRUCTURES FOR A SNAIL AND ITS TREMATODE PARASITE

Gene flow and the genetic structure of host and parasite populations are critical to the coevolutionary process, including the conditions under which antagonistic coevolution favors sexual reproduction. Here we compare the genetic structures of different populations of a freshwater New Zealand snail (Potamopyrgus antipodarum) with its trematode parasite (Microphallus sp.) using allozyme frequency data. Allozyme variation among snail populations was found to be highly structured among lakes; but for the parasite there was little allozyme structure among lake populations, suggesting much higher levels of parasite gene flow. The overall pattern of variation was confirmed with principal component analysis, which also showed that the organization of genetic differentiation for the snail (but not the parasite) was strongly related to the geographic arrangement of lakes. Some snail populations from different sides of the Alps near mountain passes were more similar to each other than to other snail populations on the same side of the Alps. Furthermore, genetic distances among parasite populations were correlated with the genetic distances among host populations, and genetic distances among both host and parasite populations were correlated with “stepping-stone” distances among lakes. Hence, the host snail and its trematode parasite seem to be dispersing to adjacent lakes in a stepping-stone fashion, although parasite dispersal among lakes is clearly greater. High parasite gene flow should help to continuously reintroduce genetic diversity within local populations where strong selection might otherwise isolate “host races.” Parasite gene flow can thereby facilitate the coevolutionary (Red Queen) dynamics that confer an advantage to sexual reproduction by restoring lost genetic variation.     

Immunological change in a parasite-impoverished environment: divergent signals from four island taxa

Dramatic declines of native Hawaiian avifauna due to the human-mediated emergence of avian malaria and pox prompted an examination of whether island taxa share a common altered immunological signature, potentially driven by reduced genetic diversity and reduced exposure to parasites. We tested this hypothesis by characterizing parasite prevalence, genetic diversity and three measures of immune response in two recently-introduced species (Neochmia temporalis and Zosterops lateralis) and two island endemics (Acrocephalus aequinoctialis and A. rimitarae) and then comparing the results to those observed in closely-related mainland counterparts. The prevalence of blood parasites was significantly lower in 3 of 4 island taxa, due in part to the absence of certain parasite lineages represented in mainland populations. Indices of genetic diversity were unchanged in the island population of N. temporalis; however, allelic richness was significantly lower in the island population of Z. lateralis while both allelic richness and heterozygosity were significantly reduced in the two island-endemic species examined. Although parasite prevalence and genetic diversity generally conformed to expectations for an island system, we did not find evidence for a pattern of uniformly altered immune responses in island taxa, even amongst endemic taxa with the longest residence times. The island population of Z. lateralis exhibited a significantly reduced inflammatory cell-mediated response while levels of natural antibodies remained unchanged for this and the other recently introduced island taxon. In contrast, the island endemic A. rimitarae exhibited a significantly increased inflammatory response as well as higher levels of natural antibodies and complement. These measures were unchanged or lower in A. aequinoctialis. We suggest that small differences in the pathogenic landscape and the stochastic history of mutation and genetic drift are likely to be important in shaping the unique immunological profiles of small isolated populations. Consequently, predicting the impact of introduced disease on the many other endemic faunas of the remote Pacific will remain a challenge.     

Molecular Detection of Hematozoa Infections in Tundra Swans Relative to Migration Patterns and Ecological Conditions at Breeding Grounds

Tundra swans (Cygnus columbianus) are broadly distributed in North America, use a wide variety of habitats, and exhibit diverse migration strategies. We investigated patterns of hematozoa infection in three populations of tundra swans that breed in Alaska using satellite tracking to infer host movement and molecular techniques to assess the prevalence and genetic diversity of parasites. We evaluated whether migratory patterns and environmental conditions at breeding areas explain the prevalence of blood parasites in migratory birds by contrasting the fit of competing models formulated in an occupancy modeling framework and calculating the detection probability of the top model using Akaike Information Criterion (AIC). We described genetic diversity of blood parasites in each population of swans by calculating the number of unique parasite haplotypes observed. Blood parasite infection was significantly different between populations of Alaska tundra swans, with the highest estimated prevalence occurring among birds occupying breeding areas with lower mean daily wind speeds and higher daily summer temperatures. Models including covariates of wind speed and temperature during summer months at breeding grounds better predicted hematozoa prevalence than those that included annual migration distance or duration. Genetic diversity of blood parasites in populations of tundra swans appeared to be relative to hematozoa prevalence. Our results suggest ecological conditions at breeding grounds may explain differences of hematozoa infection among populations of tundra swans that breed in Alaska.     

Experimental evolution reveals differences between phenotypic and evolutionary responses to population density

Group living can select for increased immunity, given the heightened risk of parasite transmission. Yet, it also may select for increased male reproductive investment, given the elevated risk of female multiple mating. Trade‐offs between immunity and reproduction are well documented. Phenotypically, population density mediates both reproductive investment and immune function in the Indian meal moth, Plodia interpunctella. However, the evolutionary response of populations to these traits is unknown. We created two replicated populations of P. interpunctella, reared and mated for 14 generations under high or low population densities. These population densities cause plastic responses in immunity and reproduction: at higher numbers, both sexes invest more in one index of immunity [phenoloxidase (PO) activity] and males invest more in sperm. Interestingly, our data revealed divergence in PO and reproduction in a different direction to previously reported phenotypic responses. Males evolving at low population densities transferred more sperm, and both males and females displayed higher PO than individuals at high population densities. These positively correlated responses to selection suggest no apparent evolutionary trade‐off between immunity and reproduction. We speculate that the reduced PO activity and sperm investment when evolving under high population density may be due to the reduced population fitness predicted under increased sexual conflict and/or to trade‐offs between pre‐ and post‐copulatory traits.     

Geographic variation in sterilizing parasite species and the Red Queen

The Red Queen hypothesis predicts that sexual reproduction should be favoured in locations where the risk of infection by virulent parasites is consistently high. When hosts are exposed to multiple parasites over their geographic range, the coevolving parasite species may vary among host populations. We surveyed 26 streams on the South Island of New Zealand to determine whether the frequency of snails ( Potamopyrgus antipodarum ) infected by various sterilizing trematode parasite species was correlated with the frequency of sexual individuals. We compared the results with a survey conducted over 20 years ago to determine whether the associations were consistent. We also evaluated different measures of parasite-mediated selection among populations, including prevalence of the most common local parasite (MCLP) species and parasite diversity to assess the best predictor of sexual reproduction among stream populations. The results showed that the relationship between male frequency and parasite infection is more geographically widespread than previously recorded. Additionally, we found that the prevalence of the MCLP was the best predictor of sex in habitats where hosts populations are infected with multiple parasites (approximately 15 trematode species). This study provides evidence that sexual snails occur more often in environments with high infection levels, and that the pattern of parasite-imposed selection is geographically variable. Support for the Red Queen may be strengthened by focussing on the MCLP, which may vary among host populations.     

Elucidating relationships between P.falciparum prevalence and measures of genetic diversity with a combined genetic-epidemiological model of malaria

There is an abundance of malaria genetic data being collected from the field, yet using these data to understand the drivers of regional epidemiology remains a challenge. A key issue is the lack of models that relate parasite genetic diversity to epidemiological parameters. Classical models in population genetics characterize changes in genetic diversity in relation to demographic parameters, but fail to account for the unique features of the malaria life cycle. In contrast, epidemiological models, such as the Ross-Macdonald model, capture malaria transmission dynamics but do not consider genetics. Here, we have developed an integrated model encompassing both parasite evolution and regional epidemiology. We achieve this by combining the Ross-Macdonald model with an intra-host continuous-time Moran model, thus explicitly representing the evolution of individual parasite genomes in a traditional epidemiological framework. Implemented as a stochastic simulation, we use the model to explore relationships between measures of parasite genetic diversity and parasite prevalence, a widely-used metric of transmission intensity. First, we explore how varying parasite prevalence influences genetic diversity at equilibrium. We find that multiple genetic diversity statistics are correlated with prevalence, but the strength of the relationships depends on whether variation in prevalence is driven by host- or vector-related factors. Next, we assess the responsiveness of a variety of statistics to malaria control interventions, finding that those related to mixed infections respond quickly (∼months) whereas other statistics, such as nucleotide diversity, may take decades to respond. These findings provide insights into the opportunities and challenges associated with using genetic data to monitor malaria epidemiology.     

Population Genetics of the Filarial Worm Wuchereria bancrofti in a Post-treatment Region of Papua New Guinea: Insights into Diversity and Life History

Background

Wuchereria bancrofti (Wb) is the primary causative agent of lymphatic filariasis (LF). Our studies of LF in Papua New Guinea (PNG) have shown that it is possible to reduce the prevalence of Wb in humans and mosquitoes through mass drug administration (MDA; diethylcarbamazine with/without ivermectin). While MDAs in the Dreikikir region through 1998 significantly reduced prevalence of Wb infection, parasites continue to be transmitted in the area.

Methods

We sequenced the Wb mitochondrial Cytochrome Oxidase 1 (CO1) gene from 16 people infected with Wb. Patients were selected from 7 villages encompassing both high and moderate annual transmission potentials (ATP). We collected genetic data with the objectives to (i) document contemporary levels of genetic diversity and (ii) distinguish between populations of parasites and hosts across the study area.

Principle Findings

We discovered 109 unique haplotypes currently segregating in the Wb parasite population, with one common haplotype present in 15 out of 16 infections. We found that parasite diversity was similar among people residing within the same village and clustered within transmission zones. For example, in the high transmission area, diversity tended to be more similar between neighboring villages, while in the moderate transmission area, diversity tended to be less similar.

Conclusions

In the Dreikikir region of PNG there are currently high levels of genetic diversity in populations of Wb. High levels of genetic diversity may complicate future MDAs in this region and the presence of dominant haplotypes will require adjustments to current elimination strategies.     

Clonal diversity of a lizard malaria parasite, Plasmodium mexicanum, in its vertebrate host, the western fence lizard: role of variation in transmission intensity over time and space

Within the vertebrate host, infections of a malaria parasite (Plasmodium) could include a single genotype of cells (single-clone infections) or two to several genotypes (multiclone infections). Clonal diversity of infection plays an important role in the biology of the parasite, including its life history, virulence, and transmission. We determined the clonal diversity of Plasmodium mexicanum, a lizard malaria parasite at a study region in northern California, using variable microsatellite markers, the first such study for any malaria parasite of lizards or birds (the most common hosts for Plasmodium species). Multiclonal infections are common (50-88% of infections among samples), and measures of genetic diversity for the metapopulation (expected heterozygosity, number of alleles per locus, allele length variation, and effective population size) all indicated a substantial overall genetic diversity. Comparing years with high prevalence (1996-1998 = 25-32% lizards infected), and years with low prevalence (2001-2005 = 6-12%) found fewer alleles in samples taken from the low-prevalence years, but no reduction in overall diversity (H = 0.64-0.90 among loci). In most cases, rare alleles appeared to be lost as prevalence declined. For sites chronically experiencing low transmission intensity (prevalence approximately 1%), overall diversity was also high (H = 0.79-0.91), but there were fewer multiclonal infections. Theory predicts an apparent excess in expected heterozygosity follows a genetic bottleneck. Evidence for such a distortion in genetic diversity was observed after the drop in parasite prevalence under the infinite alleles mutation model but not for the stepwise mutation model. The results are similar to those reported for the human malaria parasite, Plasmodium falciparum, worldwide, and support the conclusion that malaria parasites maintain high genetic diversity in host populations despite the potential for loss in alleles during the transmission cycle or during periods/locations when transmission intensity is low.     

Local parasite pressures and host genotype modulate epigenetic diversity in a mixed‐mating fish

Parasite‐mediated selection is one of the main drivers of genetic variation in natural populations. The persistence of long‐term self‐fertilization, however, challenges the notion that low genetic variation and inbreeding compromise the host's ability to respond to pathogens. DNA methylation represents a potential mechanism for generating additional adaptive variation under low genetic diversity. We compared genetic diversity (microsatellites and AFLPs), variation in DNA methylation (MS‐AFLPs), and parasite loads in three populations of Kryptolebias hermaphroditus, a predomintanly self‐fertilizing fish, to analyze the potential adaptive value of DNA methylation in relation to genetic diversity and parasite loads. We found strong genetic population structuring, as well as differences in parasite loads and methylation levels among sampling sites and selfing lineages. Globally, the interaction between parasites and inbreeding with selfing lineages influenced DNA methylation, but parasites seemed more important in determining methylation levels at the local scale.     

Spatio‐temporal variation in parasite communities maintains diversity at the major histocompatibility complex class IIβ in the endangered Rio Grande silvery minnow

Climate change will strongly impact aquatic ecosystems particularly in arid and semi‐arid regions. Fish–parasite interactions will also be affected by predicted altered flow and temperature regimes, and other environmental stressors. Hence, identifying environmental and genetic factors associated with maintaining diversity at immune genes is critical for understanding species’ adaptive capacity. Here, we combine genetic (MHC class IIβ and microsatellites), parasitological and ecological data to explore the relationship between these factors in the remnant wild Rio Grande silvery minnow (Hybognathus amarus) population, an endangered species found in the southwestern United States. Infections with multiple parasites on the gills were observed and there was spatio‐temporal variation in parasite communities and patterns of infection among individuals. Despite its highly endangered status and chronically low genetic effective size, Rio Grande silvery minnow had high allelic diversity at MHC class IIβ with more alleles recognized at the presumptive DAB1 locus compared to the DAB3 locus. We identified significant associations between specific parasites and MHC alleles against a backdrop of generalist parasite prevalence. We also found that individuals with higher individual neutral heterozygosity and higher amino acid divergence between MHC alleles had lower parasite abundance and diversity. Taken together, these results suggest a role for fluctuating selection imposed by spatio‐temporal variation in pathogen communities and divergent allele advantage in maintenance of high MHC polymorphism. Understanding the complex interaction of habitat, pathogens and immunity in protected species will require integrated experimental, genetic and field studies.     

Inbred and outbred honey bees (Apis mellifera) have similar innate immune responses

Social insect colonies provide ideal conditions for the spread of pathogens. It has been proposed that the extreme polyandry and genetic diversity seen in the colonies of some eusocial insect species is central to a colony’s defence against disease. Indeed, empirically, colonies headed by polyandrous queens have lower incidence of pathogens than genetically uniform monoandrous colonies. The mechanisms of improved resistance in genetically diverse colonies could arise from the genetic diversity among worker genotypes or from increased innate immunity arising from heterozygosity at immune gene loci within individual workers. Here, we investigate the effects of heterozygosity on two components of the honey bee (Apis mellifera) innate immune system: encapsulation and phenoloxidase (PO) activity. No significant effect of heterozygosity on immune system activity was evident for either encapsulation or PO activity. Thus, we conclude that while encapsulation and PO activity are important components of the immune response, it seems that they do not underlie the positive effects of genetic diversity on parasite and pathogen resistance in honey bees.     

Can infection by eugregarine parasites mediate species coexistence in Calopteryx damselflies?

1. Parasitism may be an important factor determining the coexistence of closely related species. Although host–parasite interactions can affect the ecology and distribution of the host species, virtually nothing is known about how other interspecific interactions affecting the host, such as competition or predation, relate to the parasite burden of the host. 2. We studied parasite‐mediated competition between two closely related Calopteryx damselflies, C. virgo L. and C. splendens Harris. We investigated a total of 31 populations, including 18 allopatric and 13 sympatric populations. We measured the occurrence of gut parasites, eugregarines. 3. We found that the prevalence of gregarines was higher in C. virgo than in C. splendens. On average, more than half of the C. virgo individuals were infected by eugregarines both in allopatric and sympatric populations. However, hardly any allopatric C. splendens populations had gregarines, but most of sympatric populations had infected individuals. 4. According to our results, co‐existence of the host species affects the likelihood of the subordinate species showing higher levels of parasitism. Interspecific aggression, lower species genetic heterozygosity, and the difference in host species immunity are proposed as possible explanations for greater parasite burdens in the inferior species at sympatric sites.     

The evolutionary significance of parasitism: do parasite‐driven genetic dynamics occur ex silico?

It has long been recognized that reciprocal antagonism might lock host and parasite populations into a process of constant change, adapting and reacting in open‐ended coevolution. A significant body of theory supports this intuition: dynamic genetic polymorphisms are a common outcome of computer simulations of host–parasite coevolution. These in silico experiments have also shown that dynamical interactions could be responsible for high levels of genetic diversity in host populations, and even be the principle determinant of rates of genetic recombination and sexuality. The evolutionary significance of parasitism depends on the strength and prevalence of parasite‐mediated selection in nature. Here I appraise whether parasitism is a pervasive agent of evolutionary change by detailing empirical evidence for selection. Although there is considerable evidence of genetic variation for resistance, and hence the potential for selection, direct observation of parasite‐driven genetic change is lacking.     

Spatially variable coevolution between a haemosporidian parasite and the MHC of a widely distributed passerine

The environment shapes host–parasite interactions, but how environmental variation affects the diversity and composition of parasite‐defense genes of hosts is unresolved. In vertebrates, the highly variable major histocompatibility complex (MHC) gene family plays an essential role in the adaptive immune system by recognizing pathogen infection and initiating the cellular immune response. Investigating MHC‐parasite associations across heterogeneous landscapes may elucidate the role of spatially fluctuating selection in the maintenance of high levels of genetic variation at the MHC. We studied patterns of association between an avian haemosporidian blood parasite and the MHC of rufous‐collared sparrows (Zonotrichia capensis) that inhabit environments with widely varying haemosporidian infection prevalence in the Peruvian Andes. MHC diversity peaked in populations with high infection prevalence, although intra‐individual MHC diversity was not associated with infection status. MHC nucleotide and protein sequences associated with infection absence tended to be rare, consistent with negative frequency‐dependent selection. We found an MHC variant associated with a ~26% decrease in infection probability at middle elevations (1501–3100 m) where prevalence was highest. Several other variants were associated with a significant increase in infection probability in low haemosporidian prevalence environments, which can be interpreted as susceptibility or quantitative resistance. Our study highlights important challenges in understanding MHC evolution in natural systems, but may point to a role of negative frequency‐dependent selection and fluctuating spatial selection in the evolution of Z. capensis MHC.