Genetic Variation of Microsatellite Loci in the Major Histocompatibility Complex (MHC) Region in the Southeast Asian House Mouse (<Emphasis Type="Italic">Mus musculus castaneus</Emphasis>)

Major histocompatibility complex (MHC) genes are the most polymorphic loci known for vertebrates. Here we employed five microsatellite loci closely linked to the MHC region in an attempt to study the amount of genetic variation in 19 populations of the southeast Asian house mouse (Mus musculus castaneus) in Taiwan. The overall polymorphism at the five loci was high (H_e = 0.713), and the level of polymorphism varied from locus to locus. Furthermore, in order to investigate if selection is operating on MHC genes in natural mouse populations, we compared the extent and pattern of genetic variation for the MHC-linked microsatellite loci (the MHC loci) with those for the microsatellite loci located outside the MHC region (the non-MHC loci). The number of alleles and the logarithm of variance in repeat number were significantly higher for the MHC loci than for the non-MHC loci, presumably reflecting linkage to a locus under balancing selection. Although three statistical tests used do not provide support for selection, their lack of support may be due to low statistical power of the tests, to weakness of selection, or to a profound effect of genetic drift reducing the signature of balancing selection. Our results also suggested that the populations in the central and the southwestern regions of Taiwan might be one part of a metapopulation structure.     

Selection on MHC-linked microsatellite loci in sheep populations

Microsatellites within the major histocompatibility complex (MHC) region have received increasing attention as proxy measures of the level of polymorphism at the Mhc genes themselves. We assessed the diversity of microsatellite loci within or in close proximity of the Mhc genes in several breeds of domestic sheep (Ovis aries) and the wild Mouflon (Ovis orientalis musimon). This was compared to variation at other microsatellite loci scattered throughout the sheep genome. Significantly higher number of alleles were observed at the MHC microsatellites. The sheep breeds studied fell into high- and low-diversity group. This grouping is not related to the agricultural use of the breeds, whether for milk, meat or wool. It is, however, correlated with the geographic origins of the breeds. Southern breeds are genetically more diverse than northern breeds. The observed heterozygosity was in most cases lower than Hardy-Weinberg expectations. The potential impact of selective breeding by man on this is discussed. Neutrality tests indicated that for most of the breeds, the distribution of alleles at the MHC-linked microsatellites are more even than would be expected if the genes were neutral and sampled from populations under drift-mutation equilibrium. Hitchhiking due to tight linkage with alleles at the MHC loci that are under balancing selection is proposed as a possible explanation for this pattern.     

Globally dispersed Y chromosomal haplotypes in wild and domestic sheep

To date, investigations of genetic diversity and the origins of domestication in sheep have utilised autosomal microsatellites and variation in the mitochondrial genome. We present the first analysis of both domestic and wild sheep using genetic markers residing on the ovine Y chromosome. Analysis of a single nucleotide polymorphism (oY1) in the SRY promoter region revealed that allele A-oY1 was present in all wild bighorn sheep (Ovis canadensis), two subspecies of thinhorn sheep (Ovis dalli), European Mouflon (Ovis musimon) and the Barbary (Ammontragis lervia). A-oY1 also had the highest frequency (71.4%) within 458 domestic sheep drawn from 65 breeds sampled from Africa, Asia, Australia, the Caribbean, Europe, the Middle East and Central Asia. Sequence analysis of a second locus, microsatellite SRYM18, revealed a compound repeat array displaying fixed differences, which identified bighorn and thinhorn sheep as distinct from the European Mouflon and domestic animals. Combined genotypic data identified 11 male-specific haplotypes that represented at least two separate lineages. Investigation of the geographical distribution of each haplotype revealed that one (H6) was both very common and widespread in the global sample of domestic breeds. The remaining haplotypes each displayed more restricted and informative distributions. For example, H5 was likely founded following the domestication of European breeds and was used to trace the recent transportation of animals to both the Caribbean and Australia. A high rate of Y chromosomal dispersal appears to have taken place during the development of domestic sheep as only 12.9% of the total observed variation was partitioned between major geographical regions.     

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

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

Re-sequencing regions of the ovine Y chromosome in domestic and wild sheep reveals novel paternal haplotypes

The male-specific region of the ovine Y chromosome (MSY) remains poorly characterized, yet sequence variants from this region have the potential to reveal the wild progenitor of domestic sheep or examples of domestic and wild paternal introgression. The 5' promoter region of the sex-determining gene SRY was re-sequenced using a subset of wild sheep including bighorn ( Ovis canadensis ), thinhorn ( Ovis dalli spp.), urial ( Ovis vignei ), argali ( Ovis ammon ), mouflon ( Ovis musimon ) and domestic sheep ( Ovis aries ). Seven novel SNPs ( oY 2– oY 8) were revealed; these were polymorphic between but not within species. Re-sequencing and fragment analysis was applied to the MSY microsatellite SRYM18 . It contains a complex compound repeat structure and sequencing of three novel size fragments revealed that a pentanucleotide element remained fixed, whilst a dinucleotide element displayed variability within species. Comparison of the sequence between species revealed that urial and argali sheep grouped more closely to the mouflon and domestic breeds than the pachyceriforms (bighorn and thinhorn). SNP and microsatellite data were combined to define six previously undetected haplotypes. Analysis revealed the mouflon as the only species to share a haplotype with domestic sheep, consistent with its status as a feral domesticate that has undergone male-mediated exchange with domestic animals. A comparison of the remaining wild species and domestic sheep revealed that O. aries is free from signatures of wild sheep introgression.     

Population genetic structure of North American thinhorn sheep (Ovis dalli)

The thinhorn sheep (Ovis dalli ssp.) provides a rare example of a North American large mammal that occupies most of its native range and maintains close to ancestral population size. There are currently two recognized subspecies, Dall's sheep (O. d. dalli) and Stone's sheep (O. d. stonei), the validity of which remains uncertain. We investigated the spatial genetic structure of thinhorn sheep populations representing both subspecies by genotyping individuals (n = 919) from across the species range at 12 variable microsatellite loci. We found high levels of genetic diversity within (HE = 0.722) and significant genetic structure among the 24 sampled areas (FST = 0.160). Genetic distance measures and Bayesian clustering analyses revealed the presence of at least eight subpopulations that are delineated by mountain range topology. A strong overall pattern of isolation-by-distance is evident across the sampling range (r = 0.75, P < 0.001) suggesting limited dispersal and extensive philopatry. Partial Mantel tests of this relationship showed mountain range distinctions represent significant barriers to gene flow (P = 0.0001), supporting the Bayesian analyses. Genetic structure was more strongly pronounced in southern Yukon and Alaska than elsewhere. We also show evidence for genetic differences between the two currently recognized thinhorn subspecies.     

Identifying footprints of directional and balancing selection in marine and freshwater three-spined stickleback (Gasterosteus aculeatus) populations

Natural selection is expected to leave an imprint on the neutral polymorphisms at the adjacent genomic regions of a selected gene. While directional selection tends to reduce within-population genetic diversity and increase among-population differentiation, the reverse is expected under balancing selection. To identify targets of natural selection in the three-spined stickleback ( Gasterosteus aculeatus ) genome, 103 microsatellite and two indel markers including expressed sequence tags (EST) and quantitative trait loci (QTL)-associated loci, were genotyped in four freshwater and three marine populations. The results indicated that a high proportion of loci (14.7%) might be affected by balancing selection and a lower proportion (2.8%) by directional selection. The strongest signatures of directional selection were detected in a microsatellite locus and two indel markers located in the intronic regions of the Eda-gene coding for the number of lateral plates. Yet, other microsatellite loci previously found to be informative in QTL-mapping studies revealed no signatures of selection. Two novel microsatellite loci ( Stn12 and Stn90 ) located in chromosomes I and VIII, respectively, showed signals of directional selection and might be linked to genomic regions containing gene(s) important for adaptive divergence. Although the coverage of the total genomic content was relatively low, the predominance of balancing selection signals is in agreement with the contention that balancing, rather than directional selection is the predominant mode of selection in the wild.     

Selective pressures on MHC class II genes in the guppy (Poecilia reticulata) as inferred by hierarchical analysis of population structure

Genes of the major histocompatibility complex, which are the most polymorphic of all vertebrate genes, are a pre‐eminent system for the study of selective pressures that arise from host–pathogen interactions. Balancing selection capable of maintaining high polymorphism should lead to the homogenization of MHC allele frequencies among populations, but there is some evidence to suggest that diversifying selection also operates on the MHC. However, the pattern of population structure observed at MHC loci is likely to depend on the spatial and/or temporal scale examined. Here, we investigated selection acting on MHC genes at different geographic scales using Venezuelan guppy populations inhabiting four regions. We found a significant correlation between MHC and microsatellite allelic richness across populations, which suggests the role of genetic drift in shaping MHC diversity. However, compared to microsatellites, more MHC variation was explained by differences between populations within larger geographic regions and less by the differences between the regions. Furthermore, among proximate populations, variation in MHC allele frequencies was significantly higher compared to microsatellites, indicating that selection acting on MHC may increase population structure at small spatial scales. However, in populations that have significantly diverged at neutral markers, the population‐genetic signature of diversifying selection may be eradicated in the long term by that of balancing selection, which acts to preserve rare alleles and thus maintain a common pool of MHC alleles.     

Genetic Variation of Major Histocompatibility Complex and Microsatellite Loci: A Comparison in Bighorn Sheep

Examining and comparing genetic variation for major histocompatibility complex (MHC) and microsatellite (MS) loci in the same individuals provides an opportunity to understand the forces influencing genetic variation. We examined five MHC and three MS loci in 235 bighorn sheep (Ovis canadensis) from 14 populations and found that both types of loci were highly variable and were in Hardy-Weinberg proportions. Mean F(ST) values for both markers were very similar and MHC and MS genetic variability was predominantly distributed within rather than among populations. However, analyses of genetic distances and tree topologies revealed different spatial patterns of variation for the two types of loci. Collectively, these results indicated that neutral forces substantially influenced MS and MHC variation, and they provided limited evidence for selection acting on the MHC.     

Signals of major histocompatibility complex overdominance in a wild salmonid population

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

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

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

No evidence for major histocompatibility complex-dependent mating patterns in a free-living ruminant population.

Conventionally, the extraordinary diversity of the vertebrate major histocompatibility complex (MHC is thought to have evolved in response to parasites and pathogens affecting fitness. More recently, reproductive mechanisms such as disassortative mating have been suggested as alternative mechanisms maintaining MHC diversity. A large unmanaged population of Soay sheep (Ovis aries L.) was used to investigate reproductive mechanisms in the maintenance of MHC diversity. Animals were sampled as new-born lambs and between 887 and 1209 individuals were typed at each of five microsatellite markers located either within or flanking the ovine MHC. All loci were in Hardy-Weinberg proportions. A novel likelihood-based approach was developed to analyse mating patterns using paternity data. No evidence for non-random mating with respect to MHC markers was found using this technique. We conclude that MHC diversity in the St Kildan Soay sheep population is unlikely to be maintained by mating preferences and that, in contrast with evidence from experimental mice populations, MHC variation plays no role in the mating structure of this population.     

Spatially and temporally fluctuating selection at non-MHC immune genes: evidence from TAP polymorphism in populations of brown trout (Salmo trutta, L.)

Temporal samples of Danish brown trout (Salmo trutta) from populations representing varying geographical scales were analysed using eight putatively neutral microsatellite loci and two microsatellite loci embedded in TAP genes (Transporter associated with Antigen Processing). These genes encode molecules that are central to the major histocompatibility complex (MHC) class I restricted antigen presentation and thus integral components in the adaptive immune system. As such, they could be influenced by selection, driven by pathogens and parasites in a manner similar to MHC genes. Analysis of allele frequencies at presumably neutral microsatellite loci revealed a temporally unstable population structure within regions, while the population structure was stable over time among regions. Analyses of the two TAP markers indicated an effect of selection at both a regional and micro-geographical spatial scale. Moreover, signals of divergent selection among temporal samples within localities suggest that selection also might fluctuate at a temporal scale. These results suggest that immune genes other than the classical MHC class I and II might be subject to selection and warrant further studies of functional polymorphism of such genes in natural populations.     

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

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

Evolution of HLA class II molecules: Allelic and amino acid site variability across populations

Analysis of the highly polymorphic beta1 domains of the HLA class II molecules encoded by the DRB1, DQB1, and DPB1 loci reveals contrasting levels of diversity at the allele and amino acid site levels. Statistics of allele frequency distributions, based on Watterson's homozygosity statistic F, reveal distinct evolutionary patterns for these loci in ethnically diverse samples (26 populations for DQB1 and DRB1 and 14 for DPB1). When examined over all populations, the DQB1 locus allelic variation exhibits striking balanced polymorphism (P < 10(-4)), DRB1 shows some evidence of balancing selection (P < 0.06), and while there is overall very little evidence for selection of DPB1 allele frequencies, there is a trend in the direction of balancing selection (P < 0.08). In contrast, at the amino acid level all three loci show strong evidence of balancing selection at some sites. Averaged over polymorphic amino acid sites, DQB1 and DPB1 show similar deviation from neutrality expectations, and both exhibit more balanced polymorphic amino acid sites than DRB1. Across ethnic groups, polymorphisms at many codons show evidence for balancing selection, yet data consistent with directional selection were observed at other codons. Both antigen-binding pocket- and non-pocket-forming amino acid sites show overall deviation from neutrality for all three loci. Only in the case of DRB1 was there a significant difference between pocket- and non-pocket-forming amino acid sites. Our findings indicate that balancing selection at the MHC occurs at the level of polymorphic amino acid residues, and that in many cases this selection is consistent across populations.     

Extensive polymorphism and geographical variation at a positively selected MHC class II B gene of the lesser kestrel (Falco naumanni)

Understanding the selective forces that shape genetic variation in natural populations remains a high priority in evolutionary biology. Genes at the major histocompatibility complex (MHC) have become excellent models for the investigation of adaptive variation and natural selection because of their crucial role in fighting off pathogens. Here we present one of the first data sets examining patterns of MHC variation in wild populations of a bird of prey, the lesser kestrel, Falco naumanni . We report extensive polymorphism at the second exon of a putatively functional MHC class II gene, Fana- DAB*1. Overall, 103 alleles were isolated from 121 individuals sampled from Spain to Kazakhstan. Bayesian inference of diversifying selection suggests that several amino acid sites may have experienced strong positive selection (ω = 4.02 per codon). The analysis also suggests a prominent role of recombination in generating and maintaining MHC diversity (ρ = 4 Nc  = 0.389 per codon, θ = 0.017 per codon). Both the Fana -DAB*1 locus and a set of eight polymorphic microsatellite markers revealed an isolation-by-distance pattern across the Western Palaearctic ( r  = 0.67; P  = 0.01 and r  = 0.50; P  = 0.04, respectively). Nonetheless, geographical variation at the MHC contrasts with relatively uniform distributions in the frequencies of microsatellite alleles. In addition, we found lower fixation rates in the MHC than those predicted by genetic drift after controlling for neutral mitochondrial sequences. Our results therefore underscore the role of balancing selection as well as spatial variations in parasite-mediated selection regimes in shaping MHC diversity when gene flow is limited.     

Genome‐wide set of SNPs reveals evidence for two glacial refugia and admixture from postglacial recolonization in an alpine ungulate

Past glaciation events have played a major role in shaping the genetic diversity and distribution of wild sheep in North America. The advancement of glaciers can isolate populations in ice‐free refugia, where they can survive until the recession of ice sheets. The major Beringian refugium is thought to have held thinhorn sheep (Ovis dalli) populations during times of glacial advance. While isolation in the major refugium can account for much of the genetic and morphological diversity seen in extant thinhorn sheep populations, mounting evidence suggests the persistence of populations in smaller minor refugia. We investigated the refugial origins of thinhorn sheep using ~10 000 SNPs obtained via a cross‐species application of the domestic sheep ovine HD BeadChip to genotype 52 thinhorn sheep and five bighorn sheep (O. canadensis) samples. Phylogenetic inference revealed a distinct lineage of thinhorn sheep inhabiting British Columbia, which is consistent with the survival of a group of thinhorn sheep in a minor refugium separate from the Beringian refugium. Isolation in separate glacial refugia probably mediated the evolution of the two thinhorn sheep subspecies, the white Dall's sheep (O. d. dalli), which persisted in Beringia, and the dark Stone's sheep (O. d. stonei), which utilized the minor refugium. We also found the first genetic evidence for admixture between sheep from different glacial refugia in south‐central Yukon as a consequence of post glacial expansion and recolonization. These results show that glaciation events can have a major role in the evolution of species inhabiting previously glaciated habitats and the need to look beyond established refugia when examining the evolutionary history of such species.     

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

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

Role of selection versus neutral processes determining genetic variation in a small mammal along a climatic gradient in southern Africa

Empirical evidence is accumulating that pathogens drive selection and explain common patterns of high immune gene (major histocompatibility complex, MHC) polymorphism. While most previous studies have identified that selection has acted over large time scales on the MHC, there still is a paucity of information in mammal species that demonstrates how processes operate on MHC genes in extant generations. Here we investigated 439 striped mouse individuals (Rhabdomys pumilio), trapped across seven different locations along a climatic gradient in southern Africa. Data from a previous study, conducted in the same study system, revealed that gastro-intestinal nematode infections were higher in individuals from study sites located within wetter climates compared to those from drier ones. In order to improve our understanding about the role of parasite-driven selection on the MHC in contemporary generations we tested for population divergences based on seven neutral microsatellite markers and the MHC DRB exon II locus. If divergences exist, we wanted to know if they are influenced by the spatial variation in parasite pressure mediated by different climatic conditions along the study site transect. Our analysis revealed an extensive polymorphism of 249 different MHC alleles and isolation-by-distance showed significant correlations at the microsatellite loci but not at the MHC. Nematode pressure was lowest at the driest site (Fish River Canyon, Namibia) and specifically this population revealed the highest divergence between MHC and microsatellite loci. We conclude that spatial variation in parasite pressure can facilitate local immune gene adaptations and thus mediate interactions of directional and balancing selection shaping MHC polymorphism in contemporary generations.