Parasite-mediated selection of major histocompatibility complex variability in wild brandt’s voles (Lasiopodomys brandtii) from Inner Mongolia, China
© Zhang and He; licensee BioMed Central Ltd. 2013
Received: 1 March 2013
Accepted: 27 June 2013
Published: 12 July 2013
Genes of the major histocompatibility complex (MHC) exhibit high levels of variability, which is believed to have arisen through pathogen-mediated selection. We investigated the relationship between parasite load and genetic diversity at selectively neutral, non-coding markers (microsatellites) and adaptive genetic variation at a functionally important part of the MHC in six independent natural populations of Brandt’s voles (Lasiopodomys brandtii) from two regions of the Xilingol Grassland area of Inner Mongolia.
Two-hundred and fifty-two voles were screened for gastrointestinal parasites, and were assessed for genetic variation. Parasite screening was done through non-invasive fecal egg counts, while allelic diversity was determined via single-stranded conformation polymorphism and DNA sequencing. We detected eight distinct helminth egg morphotypes. A total of 10 microsatellite loci were genotyped and 19 unique MHC class II B alleles were isolated. The rate of nonsynonymous substitutions (dN) exceeded the rate of synonymous substitutions (dS) at putative antigen binding sites of DRB. Neutral and adaptive genetic diversity differed between the six vole populations. To test the main pathogen-driven selection hypotheses for the maintenance of host MHC diversity and parasite species-specific co-evolutionary effects, multivariate approaches (generalized linear mixed models) were used to test for associations between the MHC class II DRB genotype and infections with nematodes. We found no evidence for heterozygote advantage, and overall heterozygosity was lower than expected in the MHC alleles. We identified an association between the parasite load and specific MHC alleles in the voles, and this pattern varied between geographic regions.
The results suggest that MHC variability in Brandt’s voles is maintained by rare allele advantage and fluctuating selection, but the data failed to show any heterozygote advantage effect. Our results add to a growing body of evidence showing that the mode and relative strength of pathogen-driven selection acting on MHC diversity varies within specific wild populations. In addition, our study contributes to the understanding of what maintains MHC diversity, of host-pathogen coevolution and of how genetic diversity is maintained in voles.
KeywordsMajor histocompatibility complex Genetic diversity Parasite-driven selection Fluctuating selection Heterozygote advantage Rare allele advantage
Genetic diversity is widely believed to influence the evolutionary and adaptive potential of populations and species . Analysis of patterns and levels of genetic variation at neutral markers, such as autosomal microsatellites and mitochondrial DNA regions, has been widely used in the last decades to infer historical events (e.g. past demographic expansions or contractions) [2–5] and geographical features (e.g. fragmentation) [6–8] in natural populations. However, studying molecular polymorphism at loci under selection is the only direct way to understand the genetics of adaptive processes . Pathogens represent very powerful agents of selection that have the potential to drive rapid changes in the genetic composition of natural host populations. In the co-evolutionary host-pathogen interplay pathogens are particularly important for maintaining host genetic variation . The role that genetic variation plays in buffering host populations from pathogens has been emphasized in several studies. These studies found associations between low levels of genetic diversity, increased pathogen susceptibility, and high pathogen loads [11–14].
In vertebrates, the genes of the major histocompatibility complex (MHC) are among the most debated candidates in the co-evolutionary process of host-parasite interactions at the molecular level . They have been studied extensively in model species under laboratory conditions, but, because of their functional importance in the immune system and mate choice, they have also become the focus of an increasing number of studies on natural populations [15–18]. MHC genes code for cell surface molecules that present self and nonself antigens to T-cells. This function enables them to play a vital role in the recognition of pathogens invading the body. The region of the molecule responsible for binding antigens, the so-called antigen-binding sites (ABS), show particularly high levels of variation, not only in the number of alleles, but also in the extent of sequence divergence between alleles. In particular, the ABS sites display more non-synonymous than synonymous substitutions that change the amino acid sequence of the peptide and thus allow binding of a diverse array of antigens . This indicates that selection processes maintain polymorphism in the functionally important regions of the MHC. Importantly, genetic diversity in the ABS facilitates binding of a diverse array of antigens to the molecule .
The exceptionally high allelic polymorphism found in the MHC loci is believed to be maintained by pathogen-mediated selection, although the relative importance of a number of nonexclusive hypotheses explaining the potential selection mechanisms that enhance or maintain adaptive genetic variation is debated [14, 21, 22]. The heterozygote advantage hypothesis  proposes that individuals heterozygous at MHC loci are able to respond to a greater range of pathogen peptides than homozygotes and, consequently, benefit from increased resistance to pathogens. Heterozygotes are, therefore, more likely to have higher relative fitness and, as a result, more MHC alleles persist, on average, in the population. This hypothesis has been used to explain the persistence of highly divergent MHC alleles over millions of years . Support for heterozygote advantage hypothesis comes mainly from mate choice studies [25, 26]; however, researchers have rarely investigated the possible advantages of MHC heterozygosity in one individual in the context of infectious diseases, such as intestinal parasite infestations (but see [15, 27–29]).
The second selection mechanism is described by the rare allele advantage hypothesis (also known as the negative frequency dependent selection hypothesis) . In this scenario, the selection pressure exerted by common parasites favors rare resistant host alleles. As those host alleles become more common, the host population exerts a reciprocal selection pressure on the parasite population, favoring other parasite genotypes to which the host has not yet adapted. With time, this could lead to continual cycling of host and parasite genotype frequencies within the population, which could maintain high levels of MHC variability [27, 28, 31].
Finally, the fluctuating selection hypothesis  proposes that spatial and temporal heterogeneity in the type and abundance of pathogens may maintain diversity at the MHC. In short-term field studies, the detection of associations between specific MHC alleles and parasite load is usually presumed to be an indicator of this selection mechanism. However, host-parasite interactions are also shaped by environmental conditions, which play an important regulating role in the distribution, transmission, and developmental success of parasites and pathogens . These conditions can influence parasite species richness, as well as the intensity of infestation in the host species. Therefore, co-evolutionary selection processes should vary not only in time but also in space, and different specific MHC alleles should have an advantage in different environments .
Over the last decades, there have been ample findings of associations between MHC alleles and parasite load, even in studies on free-ranging species under constant challenge by a diverse range of pathogens [14, 21]. However, under homogeneous parasite selection, these mechanisms alone do not explain the observed large allelic diversity at the metapopulation level. Therefore, the idea of parasites exerting divergent selection on locally adapted MHC allele pools in heterogeneous environments has been put forth to suggest how this unparalleled genetic diversity is maintained. Several recent field surveys investigating MHC variation on different geographical scales and in heterogeneous habitats have proposed that contrasting parasite communities may shape MHC composition (e.g., in mammals , in birds [36, 37] and in fish [38–40]).
Brandt's vole (Lasiopodomys brandtii) is the dominant rodent species of the typical steppe habitat extending from the central part of Inner Mongolia through the middle and east of Dornod Aimag in the Republic of Mongolia, to the southern borders of Mongolia in Transbaikalia, Russia . The distribution of this species is discontinuous in Inner Mongolia [42, 43]. However, it is a widespread species living in a diverse environment, which presumably necessitates constant adaptation to environmental change, such as encountering new parasites, but also persistence to challenges from ancient pathogen . There is significant variation in the density of these voles between years . The presumably high parasite pressure in these habitats could lead to pronounced and therefore detectable signs of otherwise more subtle mechanisms of selection .
In this study, we examined the role of parasite-mediated MHC polymorphism in six independent natural populations of Brandt’s voles, from Maodeng Livestock Farm (MD) and East Ujimqin (DWQ) of the Xilingol Grassland area of Inner Mongolia, to understand the selective mechanisms that act on MHC in response to parasitism. Our specific aims were to test for an association between: (1) individual parasite load and MHC heterozygosity (indicating heterozygote advantage) and (2) individual pathogen load and specific MHC alleles (rare allele advantage and fluctuating selection) .
We detected eight distinct helminth egg morphotypes in 252 Brandt’s voles’ fecal samples. Five of these were classified as nematodes and, among them, two nematode morphotypes were identified as Syphacia obvelata and Aspiculuris tetraptera. The remaining three morphotypes belonged to the Trichostrongylidae family. Two different cestode morphotypes were detected, which were identified as Schizorchis ochotonae and Hymenolepis nana. In addition, one trematode morphotype from the Echinostomatidae family was detected. Among the individuals examined, 94.5% had infections with one to four helminths, with most of the infections caused by nematodes (99.2% of infected individuals), whereas only 6.3% and 3.4% of the infections were caused by cestodes and trematodes, respectively. Because of the high frequency of nematode infections found in this study and the minor prevalence of cestode and trematode infections, the latter two helminths were excluded from subsequent analyses.
Mean parasite prevalence, mean species richness, and mean infection intensity for Brandt’s voles
Maodeng livestock farm (MD)
East Ujimqin (DWQ)
Population (Sample size)
M1 (n = 41)
M2 (n = 43)
M3 (n = 41)
D1 (n = 44)
D2 (n = 41)
D3 (n = 42)
Mean species richness
1.78 ± 0.80
1.74 ± 0.83
1.75 ± 0.72
2.11 ± 0.75
2.39 ± 0.97
2.03 ± 0.72
Mean infection intensity
3.95 ± 0.29
3.96 ± 0.31
3.97 ± 0.31
4.21 ± 0.10
4.24 ± 0.12
4.12 ± 0.26
Evidence for positive selection
The rates of synonymous (dS) and nonsynonymous (dN) substitutions were calculated separately for ABS and non-ABS. For the ABS, dN (0.390) was significantly higher than dS (0.072), resulting in a dN/dS ratio of 5.43 (Z > 2.613; P < 0.01). In contrast, the non-ABS ratio between nonsynonymous (dN = 0.057) and synonymous substitutions (dS = 0.038) did not significantly deviate from unity (dN/dS = 1.49, Z < 1.165; n.s.). Moreover, dN was 6.84 times higher in the ABS than in the non-ABS. This indicates positive selection processes that maintain polymorphism in the functionally important regions of the MHC.
Genetic diversity for Brandt’s voles
Estimation of pairwise genetic distance (F ST ) between Brandt’s vole populations
Association between parasite load and MHC heterozygosity
Genetic diversity effects on nematode load in Brandt’s voles as calculated by generalized linear mixed models
a) Nematode species richness
β ± SE
2.576 ± 1.369
0.038 ± 0.008
MHC H obs
0.023 ± 0.015
MHC allelic richness
−0.206 ± 0.091
b) Nematode infection intensity
β ± SE
0.176 ± 0.054
0.031 ± 0.009
MHC H obs
0.446 ± 0.172
MHC allelic richness
−0.332 ± 0.015
Similarly, the two separate GLMMs, which included either prevalence or infection intensity as response variables (from all of the nematode species combined), in addition to MHC genotypes (homozygote, heterozygote) and MLH as predictors, also did not reveal any support for heterozygote advantage. There was no support for the hypothesis that MHC heterozygous individuals are less infected than homozygotes (prevalence: β ± SE =17.683 ± 0.439, t =15.383, P = 0.328). The same applied to MLH (β ± SE = 0.635 ± 0.611, t =1.874, P = 0.267). Additionally, restricting the data to each of the three most common nematodes (S. obvelata, A. tetraptera, Trichostrongylidae) did not reveal any evidence for heterozygote advantage (all P > 0.10).
Association between parasite load and specific MHC alleles
Effects of the most abundant vole Labr-DRB* alleles on nematode prevalence and infection intensity
a) Nematode prevalence
β ± SE
Prev (Syphacia obvelata)
−0.699 ± 0.156
0.981 ± 0.214
1.343 ± 0.262
Prev (Aspiculuris tetraptera)
−0.955 ± 0.247
0.748 ± 0.250
−2.512 ± 0.208
−0.897 ± 0.362
1.194 ± 0.191
3.256 ± 0.814
b) Nematode infection intensity
β ± SE
FEC (Syphacia obvelata)
−0.592 ± 0.178
0.643 ± 0.439
3.565 ± 0.255
FEC (Aspiculuris tetraptera)
−0.678 ± 0.176
1.190 ± 0.182
2.409 ± 0.914
−0.833 ± 0.286
0.639 ± 0.342
3.042 ± 0.570
Our GLMMs also revealed that alleles associated with high or low infection intensity differed between sampling regions. We found a significant region-specific effect (MD: t = −2.35, p = 0.021; DWQ: t = 1.56, p = 0.014) of the Labr-DRB*03 allele on the intensity of infection with A. tetraptera. In MD, animals carrying Labr-DRB*03 had fewer parasites than animals without it, but in DWQ the association was the opposite.
The aim of our study was to investigate whether parasite-mediated selection could explain MHC variability in free-ranging Brandt’s voles. MHC class II-DRB exon2 polymorphism and nematode burden were surveyed and selective mechanisms that may be acting on the MHC in the presence of nematodes were tested.
In 252 individuals of L. brandtii, 19 distinct Labr-DRB alleles were detected. The alleles showed high levels of nucleotide and amino acid sequence divergence. Polymorphism was highest in the functionally important antigen recognition and binding sites of the MHC. In the ABS, significantly more nonsynonymous substitutions than synonymous substitutions were found. This is considered a clear indication of positive selection  and characteristic of proteins with antigen-presenting functions . High levels of polymorphism given by the number of alleles as well as by the sequence divergence, especially at ABS sites, are common findings in MHC genes, and were found in a variety of studies [48–50]. Because the MHC plays a major role in vertebrate immune systems, pathogen-driven selection processes are thought to be involved in the maintenance of diversity at MHC loci . Many studies of wild mammals have demonstrated a correlation between MHC diversity and resistance against pathogens. In a large survey of avian malaria in 13 populations of the house sparrow (Passer domesticus), variable selection pressures were observed to select for different host allelic lineages resulting in population-specific associations between MHC alleles and risk of infection . Similarly, rodent species that face a rich ectoparasite fauna also maintain increased allelic polymorphism at the MHC [51, 52].
It is important to note that while we monitored the community of gastrointestinal parasites in the fecal samples, we concentrated our attention on helminth parasite species because of their prevalence [53–55], their impact on fitness and mortality in a wide range of wild animal species [56–58], and because their infestation intensity can be assessed non-invasively by fecal sampling. Overall, we detected eight distinct helminth egg morphotypes and an infestation rate of 94.5%. To date, the studies of individual MHC and parasite loads in natural populations, particularly for mammals, have mostly focused on helminths [15, 49, 59–61]. However, most organisms are faced with enormous numbers of pathogens, and identifying and measuring the vast community of parasites and pathogens that can infect a natural population will be crucial . Whether results from highly simplified study systems (a single pathogen species) are applicable to more complex systems is questionable. Therefore, extending the scope of studies across a broader range of parasite taxa would enhance our understanding of MHC-parasite dynamics in natural populations. A serious challenge lies in fully characterizing the MHC and pathogen load. This is unlikely to be possible in most study systems, and even if it were, statistical analysis may be intractable. The best study systems will probably be characterized by intermediate levels of pathogen diversity and simple, well-characterized, MHC structures, thus avoiding oversimplification while retaining statistical tractability .
In our study, we found neither support for the heterozygote advantage hypothesis on the population nor on the individual level, and heterozygosity of MHC alleles was lower than expected in all of the study populations. These findings suggest little or no direct selection for MHC heterozygosity in populations at our study sites, at least for the generation of voles we sampled. Consistent low MHC heterozygosity may have arisen from underdominance . The lack of associations between MHC genetic diversity and parasite load at the population level in our study adds to the mixed results of former studies [28, 50, 63]. Studies have indicated that if hosts and pathogens share a long-term co-evolutionary history, selection via diverse pathogens causes high MHC polymorphism in a species or population, whereas low MHC polymorphism indicates the presence of relaxed pathogenic selection pressure [38, 64, 65]. Alternatively, it might have been because MHC diversity was not fully characterized as a result of methodological errors, such as null alleles. However, as the MHC sequences could be amplified from all of the study individuals, and on the basis of at least two independent polymerase chain reaction (PCR) and single-stranded conformation polymorphism (SSCP) assays, as well as forward and reverse sequence analyses, respectively, the obstacle of null alleles is improbable.
We detected an association between parasite load and specific MHC alleles in the individuals. The Labr-DRB*11 and Labr-DRB*13 alleles were significantly associated with the status of S. obvelata infection. As for A. tetraptera, voles that carried the Labr-DRB*04 allele were significantly less infected than animals without it, while Labr-DRB*12 alleles were associated with an increased probability of a higher infection intensity. Positive associations of the subfamily Trichostrongylidae and Labr-DRB*19 alleles were revealed for both parasite abundance and infection intensity. In MD, animals carrying Labr-DRB*03 had fewer parasites than animals without it, but in DWQ the association was the opposite. Differences in vole susceptibility to parasitism were not explained by differences in the geographic distribution of alleles; all above alleles were detected in at least five populations and in both sampling regions. The association between certain alleles and susceptibility, or resistance to certain parasites and spatial variation in resistance, matches the predictions for rare allele advantage, but also for fluctuating selection . Rare allele advantage arises as the evolution of new parasite and pathogen phenotypes reduces the relative fitness of common host genotypes, thereby providing a selective advantage to rare host . Evidence in favor of this hypothesis comes from multiple studies demonstrating associations between specific MHC alleles and resistance to viral (e.g., hepatitis , Puumala virus ) and bacterial (e.g., tuberculosis ) infections as well as blood-borne parasites(e.g., malaria ) and helminths [10, 71]. In reality, however, natural populations are exposed to fluctuating environmental conditions and, subsequently, host-pathogen interactions are similarly expected to vary spatiotemporally . Charbonnel and Pemberton  detected fluctuating selection at an MHC locus during a 13-year survey of a population of feral sheep, perhaps driven by interactions with parasitic nematodes. Empirical genetic evidence for geographic heterogeneity in selective pressures has been shown in fish , birds [36, 75] and mammals [35, 76]. A selection model demonstrated that temporal variation in pathogen resistance may be sufficient to maintain polymorphism in the absence of both heterozygote and rare-allele advantage . While it is generally accepted that any or both of these proposed mechanisms can play a role in shaping the distribution of MHC variation, distinguishing between them in natural populations can be challenging due to similarities in the expected genetic outcomes, requiring long-term studies of multiple populations to examine temporal changes in parasite resistance to alleles and spatiotemporal variation in the forces driving parasite abundance . Under rare-allele advantage, one would expect to see different alleles conferring resistance to the same pathogen in different populations, and for resistance to change with time, so that different alleles become associated with resistance. Under fluctuating selection, one would expect to observe external biotic and/or abiotic forces driving spatio-temporal variation in pathogen abundance, leading to distinct subsets of alleles being selected for in different populations and/or different time periods , and higher population structure at MHC relative to neutral loci.
In wild rodents, parasite prevalence and intensity are known to underlie spatial and temporal variation, as well as intrinsic factors such as the immune status, sex, and age of the host [28, 63]. By using the individual multivariate GLMM approach we were able to focus on ‘pure’ parasite driven selection mechanisms because our models allowed us to include confounding factors that could obscure the detection of MHC effects. In 252 adult Brandt’s voles, male voles were significantly heavier than females (males: 35.4 ± 1.2 g, females: 33.5 ± 1.1 g, Analysis of covariance, F = 1.2, P = 0.028). There were no significant associations between parasite load and host body mass, however, vole gender had a significant effect on nematode prevalence and intensity (Table 4). Furthermore, some aspects of the biology of Brandt’s voles suggest that a stable equilibrium between these voles and their parasites, which is a prerequisite for rare allele advantage, may rarely be achieved. Mixing of voles and parasite populations over an extensive area is possible, and the dynamic nature of this environment may mitigate the evolution of a stable equilibrium between MHC polymorphism in the host and parasite abundance and diversity. In addition, most of the parasites hosted by voles do not have direct lifecycles  and variability among populations may be dependent on the abundance of intermediate or final hosts of these parasites. In the case of Strongyloides, which have been found to be the most abundant and widespread parasites in wildlife, the abundance of intermediate hosts can affect population level differences in parasite abundance . Thus, a significant association between 1) parasite infection and specific MHC alleles, 2) marked spatial variation in parasite infection and an association with specific MHC alleles, and 3) the dynamic nature of the environment in which voles are found, all implicate a role for fluctuating selection in maintaining MHC polymorphism in voles.
A potential criticism of the present study is that we must more rigorously characterise the diversity that is being detected. The ease with which MHC variation can now be resolved could paradoxically impede progress with our understanding of MHC dynamics . A large number of studies are characterizing MHC variation from genomic DNA extracts with little or incomplete knowledge of (1) the number of loci that are actually being simultaneously PCR amplified; (2) whether any or all of these genes are actually expressed and (3) whether the variation that is resolved reflects sequence differences at structurally important regions such as the peptide-binding amino acids . Ultimately, an inability to determine the actual number of expressed loci, identify alleles as being structural MHC variants, assign alleles to loci and determine whether an individual is heterozygous or homozygous for a particular locus compromises rigorous statistical testing of association between MHC and parasite load, either by introducing variation that is not directly under selection, or precluding detection of overdominant selection . This is the first time that field population of Brandt’s voles have been typed at the MHC region, we have made every effort to initially ensure that PCR primers do actually amplify single expressed products to clarify that subsequent analysis based on genomic DNA is appropriate . Detailed molecular groundwork is required to ensure that the full complement of MHC variation is accurately assessed in future experiment. The use of next-generation sequencing for MHC screening  is likely to be a great help in terms of more accurately characterizing MHC diversity. Despite all this, it will be a particular challenge to be able to completely validate studies from mRNA extracts, especially from natural populations.
Our results are consistent with pathogen-mediated selection operating through rare allele advantage and fluctuating selection, but these two mechanisms could not be differentiated. Our results failed to show any effect of heterozygote advantage, and MHC polymorphism in wild Brandt’s voles may be constrained through underdominance. Furthermore, our results add to a growing body of evidence showing that the mode and relative strength of pathogen-driven selection acting on MHC diversity varies within specific wild populations, while the understanding of what maintains MHC diversity will also feed into our general understanding of host-pathogen coevolution and the maintenance of genetic diversity.
Traps were set before 6 a.m. and were collected after 8 a.m. in the morning. The details of each trapped vole were individually recorded (gender, body condition and body mass). Fecal samples were collected from each trap (no feces were gathered when traps contained more than one individual) during each trapping session. We sampled ear tissue with a biopsy ear punch and stored the tissue immediately in 95% ethanol until DNA isolation. Animals were handled and immediately released at the point of capture. The animal handling and sampling protocol followed the guidelines approved by the Institutional Animal Use and Care Committee, the Institute of Zoology, Chinese Academy of Sciences (CAS IAUCC).
Parasite identification and counts
Fecal material collected from each individual was stored at 4°C overnight in Petri dishes on damp blotting paper to standardize the humidity content. Thereafter, each sample was weighed and the gastrointestinal parasite load measured; we measured the fecal egg counts (FEC; number of eggs per gram feces) using a McMaster floatation technique. Following the method of Schad, Ganzhorn and Sommer , we counted two chambers of a McMaster slide and used the mean values of the individual samples. We used potassium iodide in the flotation dilution, which enhances the detection of eggs with a higher specific weight, as proposed by Meyer-Lucht and Sommer . We classified helminth eggs into morphotypes based on size and appearance and photographed them for taxonomic identification at a later stage.
For each vole, DNA extraction from the ear tissue sample was conducted using a TIANamp Genomic DNA Kit (Tiangen Biotech, Co., Ltd., Beijing, China) following the manufacturer’s protocol. All voles were typed for 10 microsatellites to assess neutral genetic diversity. These loci have previously been described for Brandt’s voles  and were amplified following the authors’ protocols .
Screening of MHC variation
We examined variation in a highly polymorphic 200 bp fragment of exon 2 of the MHC class II DRB gene, which includes the functionally important antigen binding and recognition sites. As there was no sequence data available for Brandt's voles, primers were designed after alignment of the published sequences for small animals [17, 48, 85]; these were designated BVF (up): 5’-ATTACAACAACGGGACGCA-3’ and BVR (down): 5’-CTCGTAGTTGTGTCTGCA-3’. To identify suitable primers, an extensive blast search (http://www.ncbi.nlm.nih.gov) was carried out and DRB sequences from a wide range of animal species from different phylogenetic radiations were aligned. Amplifications were conducted in a final reaction volume of 50 μl, which included 15–50 ng of DNA, 0.25 mM of each primer, 200 mM of dNTPs, 5 μl of a 10 × reaction buffer solution and 0.5 U of Taq DNA polymerase (Beijing CoWin Bioscience Co., Ltd.). The PCR instrument used was a TaKaRa Thermal Cycler Dice TP600 (TaKaRa Bio Inc., Japan) and thermal cycling started with 3 min denaturation at 95°C, followed by 30 cycles at 95°C, 54°C and 72°C for 30 s each and ended with an elongation step at 72°C for 7 min.
Alleles were separated using SSCP . SSCP is a sensitive method that can distinguish minimal allele differences; it has been widely used in human genetics and is popular in population genetics and evolutionary ecology [11, 20, 29]. PCR products were denatured at 98°C for 10 min and immediately transferred to ice for snap-cooling to produce single-strands and hinder reannealing. The ssDNA was then mixed with loading dye and loaded onto a non-denaturing 15% polyacrylamide gel according to the manufacturer’s protocol and with the following modifications: 12°C running temperature, pre-run for 5 min at 200 V maintained for 4 h at 100 V.
Gels were silver stained, scanned and processed with Quantity One (Bio-Rad Laboratories Inc., CA, USA) to align the individual band patterns. We counted bands with the same mobility as the same alleles and each of these alleles was sequenced at least once from both directions, and where possible twice, to confirm this assumption. Appropriate bands were cut from the polyacrylamide gel, eluted in Tris/borate/EDTA (TBE) buffer, and amplified using the same protocol described above. PCR products were gel purified and then sequenced.
Estimates of genetic variation
To measure the overall neutral genetic diversity per population we used MLH , and mean microsatellite d2. DRB sequences were revised manually using the BIOEDIT Sequence Alignment Editor  and aligned in GENEDOC version 2.6 . We verified the sequence identity through homology with the published MHC alleles of other species using BLAST from NCBI with a cutoff E-value of 10-6. MEGA 5  was employed to construct a phylogenetic tree of the DRB alleles based on the exon2 sequence, using the minimum evolutionary criteria , and to calculate the relative rates of nonsynonymous (dN) and synonymous (dS) substitutions according to the model of Nei and Gojobory  with the correction of Jukes and Cantor  for multiple hits. The dN/dS rates were tested for significant differences with a Z-test. Calculations were carried out separately for ABS and non-ABS, assuming concordance with antigen binding sites in the human HLA class II molecule, DR1 .
MHC genetic diversity was described by the observed heterozygosity and the allelic richness. As the observed number of alleles in a sample is highly dependent on the number of individuals sampled, we calculated the allelic richness corrected in the different sample sizes using a rarefaction index implemented in FSTAT : thereby, the expected number of alleles in each sub-sample is calculated for the number of individuals present in the smallest sample.
For both markers, differentiation across all populations and between population pairs was tested using GenePop 4.0  for microsatellites and Arlequin 3.0  for MHC. In addition, both global and pairwise estimate of FST were estimated using Arlequin following Weir and Cockerham . Observed and expected heterozygosity for both markers were calculated by Arlequin. Linkage disequilibrium between pairs of loci and deviations from Hardy–Weinberg equilibrium for each locus were also tested in Arlequin.
We used three measurements to describe the parasite burden: 1) parasite prevalence (for individuals, presence/absence of a given species; for populations, the percentage of the animals infected); 2) species richness, which was defined as the number of all of the parasite species present in one host; and 3) infection intensity, which was estimated using nematode fecal egg counts (FEC = log10 EPG; EPG: eggs/g feces).
Dissimilarities in parasite community composition between all host populations were assessed by calculating Hellinger distances using the R library vegan . Hellinger distances are based on square-rooted proportional abundances  and, therefore, reflect relative differences in parasite community composition. A permutational multivariate analysis of variance on Hellinger distances was performed using the Adonis function in vegan . Statistical significance was obtained through 1000 permutations of the data. This analysis was then repeated for each pair of host populations among MD or DWQ separately, applying a significance level corrected for multiple comparisons (3 pairwise comparisons, α = 0.017).
To investigate associations between multiple nematode infestation and the gene constitutions considering both type of markers on the population level, we used GLMMs. Models were fitted for overall mean species richness and overall mean parasite infection intensity. The models for species richness were calculated using a Poisson distribution and log-linear-link function. In the mean infection intensity models we applied a Gaussian error distribution with an identity link function. Owing to the small number of populations (n = 6), and to avoid colinearity, the four predictors of genetic diversity (microsatellite MLH, microsatellite d2, MHC heterozygosity, MHC alleles richness) were included in separate, but otherwise identical, models.
In order to test for possible associations and interactions between MHC gene and parasite loads, as well as finding support for parasite-driven selection mechanisms at the individual level, we also used GLMMs following the methods of Froeschke and Sommer . We took (a) the prevalence (presence/absence) and (b) the infection intensity (FEC) data as response variables. For the prevalence data, logistic regression models were applied with a binomial error distribution and logit link function. For the log-transformed infection intensity data, we used a Gaussian error distribution with an identity link function. To analyze relationships between genetic constitution (heterozygote advantage) and intestinal nematode susceptibility, we included the MHC genotype (homozygote or heterozygote) and microsatellite MLH for each individual as fixed factors in our model. Heterozygote host individuals and animals with a higher allele divergence should be able to recognize a broader spectrum of parasites and thus potential lower prevalence and FEC rates would be interpreted as an advantage. To test the impact of specific MHC alleles (rare allele advantage and fluctuating selection) on nematode burden, we continued with the three most abundant nematode (S. obvelata, A. tetraptera and Trichostrongylidae family). As predictors we included the presence/absence of specific MHC alleles observed in more than five individuals as fixed factors. Simplification carried out by removing variables in the order of non-significance derived the model: Parasite load (prevalence or infection intensity) ~ specific MHC allele (present or absent) + Sex + Body mass. Negative associations may be interpreted as indicating alleles conferring resistance to the parasite species, whereas positive associations indicate susceptibility to the parasite species. To examine spatial variation in individual pathogen load and specific MHC alleles, the above analysis was then repeated for each region (MD or DWQ) separately. Because parasite load is probably influenced by individual sex and body mass, we included them as explanatory variables in all of the above GLMM analyses. We included ‘population’ as a random factor in our models to consider extra sources of variation in variances through the influences of different populations and, accordingly, geographical position of each individual.
Statistical analyses were performed using the R 2.14 statistical package . We applied a modified false discovery rate procedure  to estimate the critical p value for the effect of MHC alleles. This procedure is an alternative to the Bonferroni correction and regarded as the best practical solution to the problem of multiple comparisons .
Number of synonymous substitutions per synonymous site
Number of nonsynonymous substitutions per non synonymous site
Fecal egg counts (log10 EPG)
General linear mixed models
Maodeng livestock farm
Major histocompatibility complex
Single-stranded conformation polymorphism.
Fieldwork depended on the authorization and support of the Maodeng Ecology Experimental Station. We thank all the staff from Maodeng for their assistance collecting samples. Thanks to Joselito Acosta and two anonymous reviewers for helpful comments on the manuscript. We also express our appreciation to Dr. Sandra Cheesman from Edanz Editing and Alexander J. Pritchard for helping with the language expression of the manuscript. This work was partially supported by grants from The National Key Basic Research and Development Program of China (9732007BC109103), the Science & Technology support project of the eleventh-5-year plan of China (2009BAI83B01), the General Program of National Nature Science foundation of China (31072126), and NWRC USDA-IOZ CAS Joint Project (0760621234).
- Alcaide M, Lemus JA, Blanco G, Tella JL, Serrano D, Negro JJ, Rodriguez A, Garcia-Montijano M: MHC diversity and differential exposure to pathogens in kestrels (Aves: Falconidae). Mol Ecol. 2010, 19: 691-705. 10.1111/j.1365-294X.2009.04507.x.PubMedView ArticleGoogle Scholar
- Campagna L, Benites P, Lougheed SC, Lijtmaer DA, Di Giacomo AS, Eaton MD, Tubaro PL: Rapid phenotypic evolution during incipient speciation in a continental avian radiation. P R Soc B. 2012, 279: 1847-1856. 10.1098/rspb.2011.2170.View ArticleGoogle Scholar
- Freedman AH, Thomassen HA, Buermann W, Smith TB: Genomic signals of diversification along ecological gradients in a tropical lizard. Mol Ecol. 2010, 19: 3773-3788. 10.1111/j.1365-294X.2010.04684.x.PubMedView ArticleGoogle Scholar
- Martins FM, Templeton AR, Pavan ACO, Kohlbach BC, Morgante JS: Phylogeography of the common vampire bat (Desmodus rotundus): Marked population structure, Neotropical Pleistocene vicariance and incongruence between nuclear and mtDNA markers. BMC Evol Biol. 2009, 9: 294-10.1186/1471-2148-9-294.PubMed CentralPubMedView ArticleGoogle Scholar
- Mirol PM, Routtu J, Hoikkala A, Butlin RK: Signals of demographic expansion in Drosophila virilis. BMC Evol Biol. 2008, 8: 59-10.1186/1471-2148-8-59.PubMed CentralPubMedView ArticleGoogle Scholar
- Willi Y, Van Buskirk J, Schmid B, Fischer M: Genetic isolation of fragmented populations is exacerbated by drift and selection. J Evolution Biol. 2007, 20: 534-542. 10.1111/j.1420-9101.2006.01263.x.View ArticleGoogle Scholar
- Johansson M, Primmer CR, Merila J: Does habitat fragmentation reduce fitness and adaptability? A case study of the common frog (Rana temporaria). Mol Ecol. 2007, 16: 2693-2700. 10.1111/j.1365-294X.2007.03357.x.PubMedView ArticleGoogle Scholar
- Mhemmed G, Kamel H, Chedly A: How selection fashions morphological variation in Cakile maritima: a comparative analysis of population structure using random amplified polymorphic DNA and quantitative traits. J Syst Evol. 2012, 50: 109-118. 10.1111/j.1759-6831.2011.00176.x.View ArticleGoogle Scholar
- Mona S, Crestanello B, Bankhead-Dronnet S, Pecchioli E, Ingrosso S, D'Amelio S, Rossi L, Meneguz PG, Bertorelle G: Disentangling the effects of recombination, selection, and demography on the genetic variation at a major histocompatibility complex class II gene in the alpine chamois. Mol Ecol. 2008, 17: 4053-4067. 10.1111/j.1365-294X.2008.03892.x.PubMedView ArticleGoogle Scholar
- Meyer-Lucht Y, Otten C, Puttker T, Pardini R, Metzger JP, Sommer S: Variety matters: adaptive genetic diversity and parasite load in two mouse opossums from the Brazilian Atlantic forest. Conserv Genet. 2010, 11: 2001-2013. 10.1007/s10592-010-0093-9.View ArticleGoogle Scholar
- Meyer-Lucht Y, Sommer S: Number of MHC alleles is related to parasite loads in natural populations of yellow necked mice, Apodemus flavicollis. Evol Ecol Res. 2009, 11: 1085-1097.Google Scholar
- Turner AK, Begon M, Jackson JA, Paterson S: Evidence for selection at cytokine loci in a natural population of field voles (Microtus agrestis). Mol Ecol. 2012, 21: 1632-1646. 10.1111/j.1365-294X.2012.05501.x.PubMedView ArticleGoogle Scholar
- Turner AK, Begon M, Jackson JA, Bradley JE, Paterson S: Genetic diversity in cytokines associated with immune variation and resistance to multiple pathogens in a natural rodent population. PLoS Genet. 2011, 7: e1002343-10.1371/journal.pgen.1002343.PubMed CentralPubMedView ArticleGoogle Scholar
- Sommer S: The importance of immune gene variability (MHC) in evolutionary ecology and conservation. Frontiers in Zoology. 2005, 2: 1-18. 10.1186/1742-9994-2-1.View ArticleGoogle Scholar
- Meyer-Lucht Y, Sommer S: MHC diversity and the association to nematode parasitism in the yellow-necked mouse (Apodemus flavicollis). Mol Ecol. 2005, 14: 2233-2243. 10.1111/j.1365-294X.2005.02557.x.PubMedView ArticleGoogle Scholar
- Eizaguirre C, Lenz TL, Kalbe M, Milinski M: Divergent selection on locally adapted major histocompatibility complex immune genes experimentally proven in the field. Ecol Lett. 2012, 15: 723-731. 10.1111/j.1461-0248.2012.01791.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Lenz TL, Wells K, Pfeiffer M, Sommer S: Diverse MHC IIB allele repertoire increases parasite resistance and body condition in the long-tailed giant rat (Leopoldamys sabanus). BMC Evol Biol. 2009, 9: 1-13. 10.1186/1471-2148-9-1.View ArticleGoogle Scholar
- Meyer-Lucht Y, Otten C, Puettker T, Sommer S: Selection, diversity and evolutionary patterns of the MHC class II DAB in free-ranging Neotropical marsupials. BMC Genet. BMC Genetics 2008, 9: 39-PubMed CentralPubMedView ArticleGoogle Scholar
- Brown JH, Jardetzky TS, Gorga JC, Stern LJ: Three-dimensional structure of the human class II histocompatibility antigen HLA-DR1. Nature. 1993, 364: 33-39. 10.1038/364033a0.PubMedView ArticleGoogle Scholar
- Axtner J, Sommer S: Gene duplication, allelic diversity, selection processes and adaptive value of MHC class II DRB genes of the bank vole, Clethrionomys glareolus. Immunogenetics. 2007, 59: 417-426. 10.1007/s00251-007-0205-y.PubMedView ArticleGoogle Scholar
- Spurgin LG, Richardson DS: How pathogens drive genetic diversity: MHC, mechanisms and misunderstandings. P R Soc B. 2010, 277: 979-988. 10.1098/rspb.2009.2084.View ArticleGoogle Scholar
- Bernatchez L, Landry C: MHC studies in nonmodel vertebrates: what have we learned about natural selection in 15 years?. J Evol Biol. 2003, 16: 363-377. 10.1046/j.1420-9101.2003.00531.x.PubMedView ArticleGoogle Scholar
- Doherty PC, Zinkernagel RM: Enhanced immunological surveillance in mice heterozygous at the H-2 gene complex. Nature. 1975, 256: 50-52. 10.1038/256050a0.PubMedView ArticleGoogle Scholar
- She J, Wakeland E, Boehme S: The generation and maintenance of MHC class II gene polymorphism in rodents. Immunol Rev. 1992, 113: 207-226.Google Scholar
- Sommer S: Major histocompatibility complex and mate choice in a monogamous rodent. Behav Ecol Sociobiol. 2005, 58: 181-189. 10.1007/s00265-005-0909-7.View ArticleGoogle Scholar
- Bahr A, Wilson AB: The impact of sex-role reversal on the diversity of the major histocompatibility complex: insights from the seahorse (Hippocampus abdominalis). BMC Evol Biol. 2011, 11: 121-10.1186/1471-2148-11-121.PubMed CentralPubMedView ArticleGoogle Scholar
- Schwensow N, Fietz J, Dausmann KH, Sommer S: Neutral versus adaptive genetic variation in parasite resistance: importance of major histocompatibility complex supertypes in a free-ranging primate. Heredity. 2007, 99: 265-277. 10.1038/sj.hdy.6800993.PubMedView ArticleGoogle Scholar
- Froeschke G, Sommer S: Insights into the complex associations between MHC class II DRB polymorphism and multiple gastrointestinal parasite infestations in the striped mouse. PLoS One. 2012, 7 (2): e31820-10.1371/journal.pone.0031820.PubMed CentralPubMedView ArticleGoogle Scholar
- Schwensow N, Eberle M, Sommer S: Are there ubiquitous parasite-driven major histocompatibility complex selection mechanisms in gray mouse lemurs?. Int J Primatol. 2010, 31: 519-537. 10.1007/s10764-010-9411-9.View ArticleGoogle Scholar
- Slade RW, McCallum HI: Overdominant vs. frequency-dependent selection at MHC loci. Genetics. 1992, 132: 861-862.PubMed CentralPubMedGoogle Scholar
- Westerdahl H, Asghar M, Hasselquist D, Bensch S: Quantitative disease resistance: to better understand parasite-mediated selection on major histocompatibility complex. P R Soc B. 2012, 279: 577-584. 10.1098/rspb.2011.0917.View ArticleGoogle Scholar
- Hill A: HLA associations with malaria in Africa: some implications for MHC evolution. Molecular Evolution of the Major Histocompatibility Complex. Edited by: Klein J, Klein D. 1991, Berlin: Springer-Verlag, 403-434.View ArticleGoogle Scholar
- De Bellocq JG, Morand S, Feliu C: Patterns of parasite species richness of Western Palaeartic micro-mammals: island effects. Ecography. 2002, 25: 173-183. 10.1034/j.1600-0587.2002.250205.x.View ArticleGoogle Scholar
- Hedrick PW: Pathogen resistance and genetic variation at MHC loci. Evolution. 2002, 56: 1902-1908.PubMedView ArticleGoogle Scholar
- Oliver MK, Lambin X, Cornulier T, Piertney SB: Spatio-temporal variation in the strength and mode of selection acting on major histocompatibility complex diversity in water vole (Arvicola terrestris) metapopulations. Mol Ecol. 2009, 18: 80-92.PubMedGoogle Scholar
- Ekblom R, Saether SA, Jacobsson P, Fiske P, Sahlman T, Grahn M, Kalas JA, Hoglund J: Spatial pattern of MHC class II variation in the great snipe (Gallinago media). Mol Ecol. 2007, 16: 1439-1451. 10.1111/j.1365-294X.2007.03281.x.PubMedView ArticleGoogle Scholar
- Loiseau C, Zoorob R, Robert A, Chastel O, Julliard R, Sorci G: Plasmodium relictum infection and MHC diversity in the house sparrow (Passer domesticus). P R Soc B. 2011, 278: 1264-1272. 10.1098/rspb.2010.1968.View ArticleGoogle Scholar
- Wegner K, Reusch T, Kalbe M: Multiple parasites are driving major histocompatibility complex polymorphism in the wild. J Evol Biol. 2003, 16: 224-232. 10.1046/j.1420-9101.2003.00519.x.PubMedView ArticleGoogle Scholar
- Fraser BA, Neff BD: Parasite mediated homogenizing selection at the MHC in guppies. Genetica. 2010, 138: 273-278. 10.1007/s10709-009-9402-y.PubMedView ArticleGoogle Scholar
- Eizaguirre C, Lenz TL, Sommerfeld RD, Harrod C, Kalbe M, Milinski M: Parasite diversity, patterns of MHC II variation and olfactory based mate choice in diverging three-spined stickleback ecotypes. Evol Ecol. 2011, 25: 605-622. 10.1007/s10682-010-9424-z.View ArticleGoogle Scholar
- Shenbrot GI, Krasnov BR: An atlas of the geographic distribution of the Arvicoline rodents of the world (Rodentia, Muridae: Arvicolinae). 2005, Moscow: Pensoft pressGoogle Scholar
- Shi DZ: A preliminary study on Brandt's voles' distribution region in China and its relation to vegetation and water-temperature condition. Acta Theriologica Sinica. 1988, 8: 299-306.Google Scholar
- Zhong WQ, Zhou QQ, Sun CL: The basic characteristics of the rodent pests on the pasture in Inner Mongolia and the ecolocical strategies of controlling. Acta Theriologica Sinica. 1985, 5: 241-249.Google Scholar
- Wan X, Liu W, Wang G, WANG M, ZHONG, W: Seasonal changes of the activity patterns of Brandt. s vole (Lasiopodomys brandtii) in the typical steppe in Inner Mongolia. Acta Theriologica Sinica. 2006, 26: 226-234.Google Scholar
- Wan XR, Zhang XJ, Liu W, Wang GH, Wang MJ, Zhong WQ: Social hierarchy and its seaonal changes of marked Lasiopodomys brandtii population. Chinese Journal of Ecology. 2007, 26: 359-362.Google Scholar
- Klein J, Bontrop RE, Dawkins RL: Nomenclature for the major histocompatibility complexes of different species, a proposal. Immunogenetics. 1990, 31: 217-219.PubMedGoogle Scholar
- Hughes AL, Nei M: Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection. Nature. 1988, 335: 167-170. 10.1038/335167a0.PubMedView ArticleGoogle Scholar
- Schad J, Sommer S, Ganzhorn JU: MHC variability of a small lemur in the littoral forest fragments of southeastern Madagascar. Conserv Genet. 2004, 5: 299-309.View ArticleGoogle Scholar
- Harf R, Sommer S: Association between major histocompatibility complex class II DRB alleles and parasite load in the hairy-footed gerbil, Gerbillurus paeba, in the southern Kalahari. Mol Ecol. 2005, 14: 85-91.PubMedView ArticleGoogle Scholar
- Kloch A, Baran K, Buczek M, Konarzewski M, Radwan J: MHC influences infection with parasites and winter survival in the root vole Microtus oeconomus. Evol Ecol. 2012, 27: 635-653.View ArticleGoogle Scholar
- Oliver MK, Telfer S, Piertney SB: Major histocompatibility complex (MHC) heterozygote superiority to natural multi-parasite infections in the water vole (Arvicola terrestris). Proceedings Biological sciences / The Royal Society. 2009, 276: 1119-1128. 10.1098/rspb.2008.1525.PubMed CentralPubMedView ArticleGoogle Scholar
- Cutrera AP, Zenuto RR, Lacey EA: MHC variation, multiple simultaneous infections and physiological condition in the subterranean rodent Ctenomys talarum. Infect Genet Evol. 2011, 11: 1023-1036. 10.1016/j.meegid.2011.03.016.PubMedView ArticleGoogle Scholar
- Behnke JM, Lewis JW, Zain SN, Gilbert FS: Helminth infections in Apodemus sylvaticus in southern England: interactive effects of host age, sex and year on the prevalence and abundance of infections. J Helminthol. 1999, 73: 31-44.PubMedGoogle Scholar
- Montgomery SS, Montgomery WI: Cyclic and non-cyclic dynamics in populations of the helminth parasites of wood mice, Apodemus sylvaticus. J Helminthol. 1988, 62: 78-90. 10.1017/S0022149X00011263.PubMedView ArticleGoogle Scholar
- Stefancikova A, Gajdos O, Macko JK, Tomasovicova O: Helminth fauna of small mammals in the urban and suburban area of Kosice. Biologia. 1994, 49: 147-152.Google Scholar
- Albon SD, Stien A, Irvine RJ, Langvatn R, Ropstad E, Halvorsen O: The role of parasites in the dynamics of a reindeer population. Proceedings of the Royal Society B - Biological sciences. 2002, 269: 1625-1632. 10.1098/rspb.2002.2064.PubMed CentralView ArticleGoogle Scholar
- Gulland FMD, Fox M: Epidemiology of nematode infections of Soay sheep on St Kilda. Parasitology. 1992, 105: 481-492. 10.1017/S0031182000074667.PubMedView ArticleGoogle Scholar
- Tompkins DM, Begon M: Parasites can regulate wildlife populations. Parasitol Today. 1999, 15: 311-313. 10.1016/S0169-4758(99)01484-2.PubMedView ArticleGoogle Scholar
- Paterson S, Wilson K, Pemberton JM: Major histocompatibility complex variation associated with juvenile survival and parasite resistance in a large unmanaged ungulate population (Ovis aries L.). P Natl Acad Sci USA. 1998, 95: 3714-3719. 10.1073/pnas.95.7.3714.View ArticleGoogle Scholar
- Froeschke G, Sommer S: MHC class II DRB variability and parasite load in the striped mouse (Rhabdomys pumilio) in the southern Kalahari. Mol Biol Evol. 2005, 22: 1254-1259. 10.1093/molbev/msi112.PubMedView ArticleGoogle Scholar
- Schad J, Ganzhorn JU, Sommer S: Parasite burden and constitution of major histocompatibility complex in the malagasy mouse lemur, Microcebus murinus. Evolution. 2005, 59: 439-450.PubMedView ArticleGoogle Scholar
- Pedersen AB, Fenton A: Emphasizing the ecology in parasite community ecology. Trends Ecol Evol. 2007, 22: 133-139. 10.1016/j.tree.2006.11.005.PubMedView ArticleGoogle Scholar
- Oliver MK, Telfer S, Piertney SB: Major histocompatibility complex (MHC) heterozygote superiority to natural multi-parasite infections in the water vole (Arvicola terrestris). P R Soc B. 2010, 276: 1119-1128.View ArticleGoogle Scholar
- Prugnolle F, Manica A, Charpentier M, Guegan JF, Guernier V, Balloux F: Pathogen-driven selection and worldwide HLA class I diversity. Curr Biol. 2005, 15: 1022-1027. 10.1016/j.cub.2005.04.050.PubMedView ArticleGoogle Scholar
- Alcaide M: On the relative roles of selection and genetic drift in shaping MHC variation. Mol Ecol. 2010, 19: 3842-3844. 10.1111/j.1365-294X.2010.04772.x.PubMedView ArticleGoogle Scholar
- Lively CM, Dybdahl MF: Parasite adaptation to locally common host genotypes. Nature. 2000, 405: 679-681. 10.1038/35015069.PubMedView ArticleGoogle Scholar
- Thursz MR, Kwiatkowski D, Allsopp CE, Greenwood BM, Thomas HC, Hill AV: Association between an MHC class II allele and clearance of hepatitis B virus in the Gambia. 1995, 332: 1065-1069.Google Scholar
- Deter J, Bryja J, Chaval Y, Galan M, Henttonen H, Laakkonen J, Voutilainen L, Vapalahti O, Vaheri A, Salvador AR, et al: Association between the DQA MHC class II gene and Puumala virus infection in Myodes glareolus, the bank vole. Infection, genetics and evolution: journal of molecular epidemiology and evolutionary genetics in infectious diseases. 2008, 8: 450-458. 10.1016/j.meegid.2007.07.003.PubMedView ArticleGoogle Scholar
- Sousa AO, Mazzaccaro RJ, Russell RG, Lee FK, Turner OC, Hong S, Van Kaer L, Bloom BR: Relative contributions of distinct MHC class I-dependent cell populations in protection to tuberculosis infection in mice. Proc Natl Acad Sci USA. 2000, 97: 4204-4208. 10.1073/pnas.97.8.4204.PubMed CentralPubMedView ArticleGoogle Scholar
- Bonneaud C, Perez-Tris J, Federici P, Chastel O, Sorci G: Major histocompatibility alleles associated with local resistance to malaria in a passerine. Evolution. 2006, 60: 383-389.PubMedView ArticleGoogle Scholar
- Tollenaere C, Bryja J, Galan M, Cadet P, Deter J, Chaval Y, Berthier K, Salvador AR, Voutilainen L, Laakkonen J, et al: Multiple parasites mediate balancing selection at two MHC class II genes in the fossorial water vole: insights from multivariate analyses and population genetics. J Evolution Biol. 2008, 21: 1307-1320. 10.1111/j.1420-9101.2008.01563.x.View ArticleGoogle Scholar
- Jensen LF, Hansen MM, Mensberg KL, Loeschcke V: Spatially and temporally fluctuating selection at non-MHC immune genes: evidence from TAP polymorphism in populations of brown trout (Salmo trutta, L.). Heredity. 2008, 100: 79-91. 10.1038/sj.hdy.6801067.PubMedView ArticleGoogle Scholar
- Charbonnel N, Pemberton J: A long-term genetic survey of an ungulate population reveals balancing selection acting on MHC through spatial and temporal fluctuations in selection. Heredity (Edinb). 2005, 95: 377-388. 10.1038/sj.hdy.6800735.View ArticleGoogle Scholar
- Miller KM, Kaukinen KH, Beacham TD, Withler RE: Geographic heterogeneity in natural selection on an MHC locus in sockeye salmon. Genetica. 2001, 111: 237-257. 10.1023/A:1013716020351.PubMedView ArticleGoogle Scholar
- Alcaide M, Edwards SV, Negro JJ, Serrano D, Tella JL: Extensive polymorphism and geographical variation at a positively selected MHC class IIB gene of the lesser kestrel (Falco naumanni). Mol Ecol. 2008, 17: 2652-2665. 10.1111/j.1365-294X.2008.03791.x.PubMedView ArticleGoogle Scholar
- Kloch A, Babik W, Bajer A, Sinski E, Radwan J: Effects of an MHC-DRB genotype and allele number on the load of gut parasites in the bank vole Myodes glareolus. Mol Ecol. 2010, 19 (Suppl 1): 255-265.PubMedView ArticleGoogle Scholar
- Friberg IM, Bradley JE, Jackson JA: Macroparasites, innate immunity and immunoregulation: developing natural models. Trends Parasitol. 2010, 26: 540-549. 10.1016/j.pt.2010.06.010.PubMedView ArticleGoogle Scholar
- Raharivololona BM, Ganzhorn JU: Seasonal variations in gastrointestinal parasites excreted by the gray mouse lemur Microcebus murinus in Madagascar. Endangered Species Research. 2010, 11: 113-122.View ArticleGoogle Scholar
- Piertney SB, Oliver MK: The evolutionary ecology of the major histocompatibility complex. Heredity (Edinb). 2006, 96: 7-21.Google Scholar
- Piertney SB: Major histocompatibility complex B-LB gene variation in red grouse Lagopus lagopus scoticus. Wildlife Biol. 2003, 9: 251-259.Google Scholar
- Babik W, Taberlet P, Ejsmond MJ, Radwan J: New generation sequencers as a tool for genotyping of highly polymorphic multilocus MHC system. Mol Ecol Resour. 2009, 9: 713-719. 10.1111/j.1755-0998.2009.02622.x.PubMedView ArticleGoogle Scholar
- Ren XT, Shen G, Wang ZL, LU JQ: Effects of road and grazing on spatiotemporal distribution of Brandt's vole population in Xilin Gol grassland of Inner Mongolia. Chinese Journal of Ecology. 2011, 30: 2245-2249.Google Scholar
- Wang D, Shi DZ: Isolation and characterization of polymorphic microsatellite loci from Brandt's voles (Lasiopodomys brandtii). Molecular Ecology Notes. 2007, 7: 671-673. 10.1111/j.1471-8286.2006.01673.x.View ArticleGoogle Scholar
- Wang D, Guo Y, Shi D: Genetic structure of Brandt's vole (Lasiopodomys brandtii) populations in Inner Mongolia, China, based on microsatellite analysis. Conservation Genetics. 2011, 12: 659-10.1007/s10592-010-0172-y.View ArticleGoogle Scholar
- Sommer S, Schwab D, Ganzhorn JU: MHC diversity of endemic Malagasy rodents in relation to geographic range and social system. Behav Ecol Sociobiol. 2002, 51: 214-221. 10.1007/s00265-001-0432-4.View ArticleGoogle Scholar
- Coltman D, Pilkington J, Smith J, Pemberton J: Parasite-mediated selection against inbred soay sheep in a free-living, island population. Evolution. 1999, 53: 1259-1267. 10.2307/2640828.View ArticleGoogle Scholar
- Coulson T, Pemberton J, Albon S, Beaumont M, Marshall T: Microsatellites reveal heterosis in red deer. Proc R Soc Lond B Biol Sci. 1998, 265: 489-495. 10.1098/rspb.1998.0321.View ArticleGoogle Scholar
- Hall TA: BIOEDIT: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Nicholas KB, HBJ Nicholas: GENEDOC: a tool for editing and annotating multiple sequence alignments. 1997, http://www.nrbsc.org/gfx/genedoc/,Google Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.PubMedView ArticleGoogle Scholar
- Nei M, Gojobory T: Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol. 1986, 3: 418-426.PubMedGoogle Scholar
- Jukes TH, Cantor CR: Evolution of protein molecules. 1969, New York: Academic PressView ArticleGoogle Scholar
- Goudet J: FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). http://www2.unil.ch/popgen/softwares/fstat.htm,
- Rousset F: enepop'007: a complete re-implementation of the genepop software for Windows and Linux. Molecular ecology resources. 2008, 8: 103-106. 10.1111/j.1471-8286.2007.01931.x.PubMedView ArticleGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform. 2005, 1: 47-50.Google Scholar
- Weir BS, Cockerham CC: Estimating F-Statistics for the Analysis of Population-Structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O'Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H: vegan: Community Ecology Package. R package version 2.0-7. http://CRAN.R-project.org/package=vegan,
- Legendre P, Gallagher E: Ecologically meaningful transformations for ordination of species data. Oecologia. 2001, 129: 271-280. 10.1007/s004420100716.View ArticleGoogle Scholar
- Anderson MJ: A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001, 26: 32-46.Google Scholar
- R Development Core Team: R: A language and environment for statistical computing. http://www.R-project.org,
- Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001, 29: 1165-1188. 10.1214/aos/1013699998.View ArticleGoogle Scholar
- Narum SR: Beyond Bonferroni: less conservative analyses for conservation genetics. Conservation Genetics. 2006, 7: 783-787. 10.1007/s10592-005-9056-y.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.