Disentangling the effects of recombination,selection, and demography on the genetic variation at a major histocompatibility complex class II gene in the alpine chamois

The major histocompatibility complex (MHC) harbours some of the most polymorphic loci in vertebrate genomes. MHC genes are thought to be subject to some form of balancing selection, most likely pathogen‐mediated selection. Hence, MHC genes are excellent candidates for exploring adaptive processes. In this study, we investigated the genetic variation at exon 2 of the DRB class II MHC locus in 191 alpine chamois (Rupicapra rupicapra) from 10 populations in the eastern Alps of Italy. In particular, we were interested in distinguishing and estimating the relative impact of selective and demographic factors, while taking into account the confounding effect of recombination. The extremely high d_n/d_s ratio and the presence of trans‐species polymorphisms suggest that a strong long‐term balancing selection effect has been operating at this locus throughout the evolutionary history of this species. We analysed patterns of genetic variation within and between populations, and the mitochondrial D‐loop polymorphism patterns were analysed to provide a baseline indicator of the effects of demographic processes. These analyses showed that (i) the chamois experienced a demographic decline in the last 5000–30 000 years, most likely related to the postglacial elevation in temperature; (ii) this demographic process can explain the results of neutrality tests applied to MHC variation within populations, but cannot justify the much weaker divergence between populations implied by MHC as opposed to mitochondrial DNA; (iii) similar sets of divergent alleles are probably maintained with similar frequencies by balancing selection in different populations, and this mechanism is also operating in small isolated populations, which are strongly affected by drift.     

A new test suggests hundreds of amino acid polymorphisms in humans are subject to balancing selection

The role that balancing selection plays in the maintenance of genetic diversity remains unresolved. Here, we introduce a new test, based on the McDonald–Kreitman test, in which the number of polymorphisms that are shared between populations is contrasted to those that are private at selected and neutral sites. We show that this simple test is robust to a variety of demographic changes, and that it can also give a direct estimate of the number of shared polymorphisms that are directly maintained by balancing selection. We apply our method to population genomic data from humans and provide some evidence that hundreds of nonsynonymous polymorphisms are subject to balancing selection.

What maintains genetic variation remains an unresolved mystery. This study describes the development of a new test and its application to human population genomic data, suggesting that natural selection may have a much more important role than previously thought, with hundreds of non-synonymous polymorphisms subject to balancing selection.     

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

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

Evidence of amino acid diversity-enhancing selection within humans and among primates at the candidate sperm-receptor gene PKDREJ

Sperm-egg interaction is a crucial step in fertilization, yet the identity of most interacting sperm-egg proteins that mediate this process remains elusive. Rapid evolution of some fertilization proteins has been observed in a number of species, including evidence of positive selection in the evolution of components of the mammalian egg coat. The rapid evolution of the egg-coat proteins could strongly select for changes on the sperm receptor, to maintain the interaction. Here, we present evidence that positive selection has driven the evolution of PKDREJ, a candidate sperm receptor of mammalian egg-coat proteins. We sequenced PKDREJ from a panel of 14 primates, including humans, and conducted a comparative maximum-likelihood analysis of nucleotide changes and found evidence of positive selection. An additional panel of 48 humans was surveyed for nucleotide polymorphisms at the PKDREJ locus. The regions predicted to have been subject to adaptive evolution among primates show several amino acid polymorphisms within humans. The distribution of polymorphisms suggests that balancing selection may maintain diverse PKDREJ alleles in some populations. It remains unknown whether there are functional differences associated with these diverse alleles, but their existence could have consequences for human fertility.     

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

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

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

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

Can balancing selection on MHC loci counteract genetic drift in small fragmented populations of black grouse?

The ability of natural populations to adapt to new environmental conditions is crucial for their survival and partly determined by the standing genetic variation in each population. Populations with higher genetic diversity are more likely to contain individuals that are better adapted to new circumstances than populations with lower genetic diversity. Here, we use both neutral and major histocompatibility complex (MHC) markers to test whether small and highly fragmented populations hold lower genetic diversity than large ones. We use black grouse as it is distributed across Europe and found in populations with varying degrees of isolation and size. We sampled 11 different populations; five continuous, three isolated, and three small and isolated. We tested patterns of genetic variation in these populations using three different types of genetic markers: nine microsatellites and 21 single nucleotide polymorphisms (SNPs) which both were found to be neutral, and two functional MHC genes that are presumably under selection. The small isolated populations displayed significantly lower neutral genetic diversity compared to continuous populations. A similar trend, but not as pronounced, was found for genotypes at MHC class II loci. Populations were less divergent at MHC genes compared to neutral markers. Measures of genetic diversity and population genetic structure were positively correlated among microsatellites and SNPs, but none of them were correlated to MHC when comparing all populations. Our results suggest that balancing selection at MHC loci does not counteract the power of genetic drift when populations get small and fragmented.     

A major genetic component of BSE susceptibility

Background  

Coding variants of the prion protein gene (PRNP) have been shown to be major determinants for the susceptibility to transmitted prion diseases in humans, mice and sheep. However, to date, the effects of polymorphisms in the coding and regulatory regions of bovine PRNP on bovine spongiform encephalopathy (BSE) susceptibility have been considered marginal or non-existent. Here we analysed two insertion/deletion (indel) polymorphisms in the regulatory region of bovine PRNP in BSE affected animals and controls of four independent cattle populations from UK and Germany.     

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.     

Balancing selection and heterozygote advantage in major histocompatibility complex loci of the bottlenecked Finnish wolf population

Maintaining effective immune response is an essential factor in the survival of small populations. One of the most important immune gene regions is the highly polymorphic major histocompatibility complex (MHC). We investigated how a population bottleneck and recovery have influenced the diversity and selection in three MHC class II loci, DLA‐DRB1, DLA‐DQA1 and DLA‐DQB1, in the Finnish wolf population. We studied the larger Russian Karelian wolf population for comparison and used 17 microsatellite markers as reference loci. The Finnish and Karelian wolf populations did not differ substantially in their MHC diversities ( = 0.047, P = 0.377), but differed in neutral microsatellite diversities ( = 0.148, P = 0.008). MHC allele frequency distributions in the Finnish population were more even than expected under neutrality, implying balancing selection. In addition, an excess of nonsynonymous compared to synonymous polymorphisms indicated historical balancing selection. We also studied association between helminth (Trichinella spp. and Echinococcus canadensis) prevalence and MHC diversity at allele and SNP level. MHC‐heterozygous wolves were less often infected by Trichinella spp. and carriers of specific MHC alleles, SNP haplotypes and SNP alleles had less helminth infections. The associated SNP haplotypes and alleles were shared by different MHC alleles, which emphasizes the necessity of single‐nucleotide‐level association studies also in MHC. Here, we show that strong balancing selection has had similar effect on MHC diversities in the Finnish and Russian Karelian wolf populations despite significant genetic differentiation at neutral markers and small population size in the Finnish population.     

A new theory of MHC evolution: beyond selection on the immune genes

The major histocompatibility complex (MHC) is a dense region of immune genes with high levels of polymorphism, which are arranged in haplotype blocks. Traditional models of balancing selection (i.e. overdominance and negative frequency dependence) were developed to study the population genetics of single genes. However, the MHC is a multigene family surrounded by linked (non-neutral) polymorphisms, and not all of its features are well explained by these models. For example, (i) the high levels of polymorphism in small populations, (ii) the unexpectedly large genetic differentiation between populations, (iii) the shape of the allelic genealogy associated with trans-species evolution, and (iv) the close associations between particular MHC (human leucocyte antigen, HLA) haplotypes and the approximately 100 pathologies in humans. Here, I propose a new model of MHC evolution named Associative Balancing Complex evolution that can explain these phenomena. The model proposes that recessive deleterious mutations accumulate as a 'sheltered load' nearby MHC genes. These mutations can accumulate because (i) they are rarely expressed as homozygotes given the high MHC gene diversity and (ii) purifying selection is inefficient with low recombination rates (cf. Muller's ratchet). Once fixed, these mutations add to balancing selection and further reinforce linkage through epistatic selection against recombinants.     

Trans-species polymorphism and evidence of selection on class II MHC loci in tuco-tucos (Rodentia: Ctenomyidae)

Balancing selection acting over the evolutionary history of a lineage can result in the retention of alleles among species for longer than expected under neutral evolution. The associated pattern of trans-species polymorphism, in which similar or even identical alleles are shared among species, is often used to infer that balancing selection has occurred. The genes of the major histocompatibility complex (MHC) are thought to be subject to balancing selection that maintains alleles associated with response to specific pathogens. To explore the role of balancing selection in shaping MHC diversity in ctenomyid rodents, we examined allelic variability at the class II DRB and DQA loci in 18 species in the genus Ctenomys. Previous studies of four of these species had revealed significant within-population evidence of positive selection on MHC loci. The current study expands upon these analyses to (1) evaluate among-species evidence of positive selection and (2) explore the potential for balancing selection on MHC genes. Interspecific nucleotide sequence variation revealed significant evidence of positive selection on the DRB and DQA loci. At the same time, comparisons of phylogenetic trees for these MHC loci with a putative species tree based on mitochondrial sequence data revealed multiple examples of trans-specific polymorphism, including sharing of identical DRB and DQA alleles among distantly related species of Ctenomys. These findings suggest that MHC genes in these animals have historically been subject to balancing selection and yield new insights into the complex suite of forces shaping MHC diversity in free-living vertebrates.     

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.     

Genetic variation of the MHC DQB locus in the finless porpoise (Neophocaena phocaenoides)

The Major Histocompatibility Complex (MHC) is a large multigene coding for glycoproteins that play a key role in the initiation of immune responses in vertebrates. The exon 2 region of the MHC DQB locus was analyzed using 160 finless porpoises from 5 populations in Japanese waters. The 5 populations were based on a previous mitochondrial DNA control region analysis, which showed distinct geographical separation. Eight DQB alleles were detected, and the geographical distribution of the alleles indicated that most of them are shared among the populations. Heterozygosity of the DQB alleles in each population ranged from 0.55 to 0.78, and for all 5 populations was 0.78. Low MHC variability is not a common feature in marine mammals, but the finless porpoise populations inhabiting coastal waters had a relatively high MHC heterozygosity. Balancing selection in the MHC DQB alleles of the finless porpoise was indicated by the higher rate of nonsynonymous than synonymous substitutions for PBR; however, an excess of hetrozygotes compared to expectation was not observed. This suggests that the MHC DQB locus in the finless porpoise may have been under balancing selection for a long evolutionary time period, and is influenced by genetic drift beyond the effect of balancing selection for short time periods in small local populations.     

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.     

Density-related changes in selection pattern for major histocompatibility complex genes in fluctuating populations of voles

Host-pathogen interactions are of particular interest in studies of the interplay between population dynamics and natural selection. The major histocompatibility complex (MHC) genes of demographically fluctuating species are highly suitable markers for such studies, because they are involved in initiating the immune response against pathogens and display a high level of adaptive genetic variation. We investigated whether two MHC class II genes (DQA1, DRB) were subjected to contemporary selection during increases in the density of fossorial water vole (Arvicola terrestris) populations, by comparing the neutral genetic structure of seven populations with that estimated from MHC genes. Tests for heterozygosity excess indicated that DQA1 was subject to intense balancing selection. No such selection operated on neutral markers. This pattern of selection became more marked with increasing abundance. In the low-abundance phase, when populations were geographically isolated, both overall differentiation and isolation-by-distance were more marked for MHC genes than for neutral markers. Model-based simulations identified DQA1 as an outlier (i.e. under selection) in a single population, suggesting the action of local selection in fragmented populations. The differences between MHC and neutral markers gradually disappeared with increasing effective migration between sites. In the high-abundance year, DQA1 displayed significantly lower levels of overall differentiation than the neutral markers. This gene therefore displayed stronger homogenization than observed under drift and migration alone. The observed signs of selection were much weaker for DRB. Spatial and temporal fluctuations in parasite pressure and locus-specific selection are probably the most plausible mechanisms underlying the observed changes in selection pattern during the demographic cycle.     

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

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

Computer simulation analysis suggests weak balancing selection operative at the MICA locus

A high degree of polymorphism has been reported at the major histocompatibility class I chain-related gene A (MICA) locus, which is located 46 kb away from HLA-Bin the human major histocompatibility complex (MHC) class I region. Although it is known that the polymorphisms at the conventional MHC class I loci have been maintained by balancing selection, it is unclear whether positive natural selection is also operative in maintaining the polymorphism at the MICA locus. In order to explain the degree of polymorphism at the MICA locus, a computer simulation study was carried out. The high degree of polymorphism at the MICA locus (heterozygosity and number of polymorphic residues) could not be explained solely by balancing selection at the HLA-B locus even if no recombination was assumed between MICA and HLA-B. Although there is no definite evidence indicating that balancing selection is operative at the MICA locus, our results suggest that the MICA gene is subject to weak balancing selection.     

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

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