Disentangling the roles of natural selection and genetic drift in shaping variation at MHC immunity genes

The major histocompatibility complex (MHC) forms an integral component of the vertebrate immune response and, due to strong selection pressures, is one of the most polymorphic regions of the entire genome. Despite over 15 years of research, empirical studies offer highly contradictory explanations of the relative roles of different evolutionary forces, selection and genetic drift, acting on MHC genes during population bottlenecks. Here, we take a meta-analytical approach to quantify the results of studies into the effects of bottlenecks on MHC polymorphism. We show that the consequences of selection acting on MHC loci prior to a bottleneck event, combined with drift during the bottleneck, will result in overall loss of MHC polymorphism that is ～15% greater than loss of neutral genetic diversity. These results are counter to general expectations that selection should maintain MHC polymorphism, but do agree with the results of recent simulation models and at least two empirical studies. Notably, our results suggest that negative frequency-dependent selection could be more important than overdominance for maintaining high MHC polymorphism in pre-bottlenecked populations.     

Major histocompatibility complex variation in insular populations of the Egyptian vulture: inferences about the roles of genetic drift and selection

Insular populations have attracted the attention of evolutionary biologists because of their morphological and ecological peculiarities with respect to their mainland counterparts. Founder effects and genetic drift are known to distribute neutral genetic variability in these demes. However, elucidating whether these evolutionary forces have also shaped adaptive variation is crucial to evaluate the real impact of reduced genetic variation in small populations. Genes of the major histocompatibility complex (MHC) are classical examples of evolutionarily relevant loci because of their well-known role in pathogen confrontation and clearance. In this study, we aim to disentangle the partial roles of genetic drift and natural selection in the spatial distribution of MHC variation in insular populations. To this end, we integrate the study of neutral (22 microsatellites and one mtDNA locus) and MHC class II variation in one mainland (Iberia) and two insular populations (Fuerteventura and Menorca) of the endangered Egyptian vulture (Neophron percnopterus). Overall, the distribution of the frequencies of individual MHC alleles (n=17 alleles from two class II B loci) does not significantly depart from neutral expectations, which indicates a prominent role for genetic drift over selection. However, our results point towards an interesting co-evolution of gene duplicates that maintains different pairs of divergent alleles in strong linkage disequilibrium on islands. We hypothesize that the co-evolution of genes may counteract the loss of genetic diversity in insular demes, maximize antigen recognition capabilities when gene diversity is reduced, and promote the co-segregation of the most efficient allele combinations to cope with local pathogen communities.     

Genetic drift outweighs balancing selection in shaping post-bottleneck major histocompatibility complex variation in New Zealand robins (Petroicidae)

The Chatham Island black robin, Petroica traversi, is a highly inbred, endangered passerine with extremely low levels of variation at hypervariable neutral DNA markers. In this study we investigated variation in major histocompatibility complex (MHC) class II genes in both the black robin and its nonendangered relative, the South Island robin Petroica australis australis. Previous studies have shown that Petroica have at least four expressed class II B MHC genes. In this study, the sequences of introns flanking exon 2 of these loci were characterized to design primers for peptide-binding region (PBR) sequence analysis. Intron sequences were comprised of varying numbers of repeated units, with highly conserved regions immediately flanking exon 2. Polymerase chain reaction primers designed to this region amplified three or four sequences per black robin individual, and eight to 14 sequences per South Island robin individual. MHC genes are fitness-related genes thought to be under balancing selection, so they may be more likely to retain variation in bottlenecked populations. To test this, we compared MHC variation in the black robin with artificially bottlenecked populations of South Island robin, and with their respective source populations, using restriction fragment length polymorphism analyses and DNA sequencing of the PBR. Our results indicate that the black robin is monomorphic at class II B MHC loci, while both source and bottlenecked populations of South Island robin have retained moderate levels of variation. Comparison of MHC variation with minisatellite DNA variation indicates that genetic drift outweighs balancing selection in determining MHC diversity in the bottlenecked populations. However, balancing selection appears to influence MHC diversity over evolutionary timescales, and the effects of gene conversion are evident.     

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.     

The relative role of drift and selection in shaping the human skull

Human populations across the world vary greatly in cranial morphology. It is highly debated to what extent this variability has accumulated through neutral processes (genetic drift) or through natural selection driven by climate. By taking advantage of recent work showing that geographic distance along landmasses is an excellent proxy for neutral genetic differentiation, we quantify the relative role of drift versus selection in an exceptionally large dataset of human skulls. We show that neutral processes have been much more important than climate in shaping the human cranium. We further demonstrate that a large proportion of the signal for natural selection comes from populations from extremely cold regions. More generally, we show that, if drift is not explicitly accounted for, the effect of natural selection can be greatly overestimated. Am J Phys Anthropol, 2010. © 2009 Wiley‐Liss, Inc.     

Diversifying selection on MHC class I in the house sparrow (Passer domesticus)

Genes of the major histocompatibility complex (MHC) are the most polymorphic loci known in vertebrates. Two main hypotheses have been put forward to explain the maintenance of MHC diversity: pathogen-mediated selection and MHC-based mate choice. Host–parasite interactions can maintain MHC diversity via frequency-dependent selection, heterozygote advantage, and diversifying selection (spatially and/or temporally heterogeneous selection). In this study, we wished to investigate the nature of selection acting on the MHC class I across spatially structured populations of house sparrows ( Passer domesticus ) in France. To infer the nature of the selection, we compared patterns of population differentiation based on two types of molecular markers: MHC class I and microsatellites. This allowed us to test whether the observed differentiation at MHC genes merely reflects demographic and/or stochastic processes. At the global scale, diversifying selection seems to be the main factor maintaining MHC diversity in the house sparrow. We found that (i) overall population differentiation at MHC was stronger than for microsatellites, (ii) MHC marker showed significant isolation by distance. In addition, the slope of the regression of F _ST on geographical distance was significantly steeper for MHC than for microsatellites due to a stronger pairwise differentiation between populations located at large geographical distances. These results are in agreement with the hypothesis that spatially heterogeneous selective pressures maintain different MHC alleles at local scales, possibly resulting in local adaptation.     

A global analysis of selection at the avian MHC

Recent advancements in sequencing technology have resulted in rapid progress in the study of the major histocompatibility complex (MHC) in non‐model avian species. Here, we analyze a global dataset of avian MHC class I and class II sequences (ca. 11,000 sequences from over 250 species) to gain insight into the processes that govern macroevolution of MHC genes in birds. Analysis of substitution rates revealed striking differences in the patterns of diversifying selection between passerine and non‐passerine birds. Non‐passerines showed stronger selection at MHC class II, which is primarily involved in recognition of extracellular pathogens, while passerines showed stronger selection at MHC class I, which is involved in recognition of intracellular pathogens. Positions of positively selected amino‐acid residues showed marked discrepancies with peptide‐binding residues (PBRs) of human MHC molecules, suggesting that using a human classification of PBRs to assess selection patterns at the avian MHC may be unjustified. Finally, our analysis provided evidence that indel mutations can make a substantial contribution to adaptive variation at the avian MHC.     

Demographic processes shaping genetic variation of the solitarious phase of the desert locust

Between plagues, the solitarious desert locust (Schistocerca gregaria) is generally thought to exist as small populations, which are particularly prone to extinction events in arid regions of Africa and Asia. Given the high genetic structuring observed in one geographical area (the Eritrean coast) by former authors, a metapopulation dynamics model involving repeated extinction and colonization events was favoured. In this study, we assessed the validity of a demographic scenario involving temporary populations of the solitarious phase of the desert locust by analysing large‐scale population genetic data. We scored 24 microsatellites in 23 solitarious population samples collected over most of the species range during remission. We found very little genetic structuring and little evidence of declining genetic diversity. A Bayesian clustering method distinguished four genetically differentiated units. Three groups were largely consistent with three population samples which had undergone recent bottleneck events. Nevertheless, the last genetically homogeneous unit included all individuals from the remaining 18 population samples and did not show evidence of demographic disequilibrium. An approximate Bayesian computation treatment indicated a large population size for this main genetic group, moderately reduced between plague and remission but still containing tens of thousands of individuals. Our results diverge from the hypothesis of a classical metapopulation dynamics model. They instead support the scenario in which large populations persist in the solitarious phase of the desert locust.     

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.     

Retracted: MHC diversity and differential exposure to pathogens in kestrels (Aves: Falconidae)

Pathogen diversity is thought to drive major histocompatibility complex (MHC) polymorphism given that host’s immune repertories are dependent on antigen recognition capabilities. Here, we surveyed an extensive community of pathogens (n = 35 taxa) and MHC diversity in mainland versus island subspecies of the Eurasian kestrel Falco tinnunculus and in a sympatric mainland population of the phylogenetically related lesser kestrel Falco naumanni. Insular subspecies are commonly exposed to impoverished pathogen communities whilst different species’ ecologies and contrasting life‐history traits may lead to different levels of pathogen exposure. Although specific host traits may explain differential particular infections, overall pathogen diversity, richness and prevalence were higher in the truly cosmopolitan, euriphagous and long‐distance disperser Eurasian kestrel than in the estenophagous, steppe‐specialist, philopatric but long‐distance migratory lesser kestrel. Accordingly, the continental population of Eurasian kestrels displayed a higher number (64 vs. 49) as well as more divergent alleles at both MHC class I and class II loci. Detailed analyses of amino acid diversity revealed that significant differences between both species were exclusive to those functionally important codons comprising the antigen binding sites. The lowest pathogen burdens and the smallest but still quite divergent set of MHC alleles (n = 16) were found in island Eurasian kestrels, where the rates of allele fixation at MHC loci seem to have occurred faster than at neutral markers. The results presented in this study would therefore support the role of pathogen diversity and abundance in shaping patterns of genetic variation at evolutionary relevant MHC genes.     

Experimental demonstration of a causal relationship between heterogeneity of selection and genetic differentiation in quantitative traits

Comparisons of estimates of genetic differentiation at molecular markers (F(ST)) and at quantitative traits (Q(ST)) are a means of inferring the level and heterogeneity of selection in natural populations. However, such comparisons are questionable because they require that the influence of drift and selection on Q(ST) be detectable over possible background influences of environmental or nonadditive genetic effects on Q(ST)-values. Here we test this using an experimental evolution approach in metapopulations of Arabidopsis thaliana experiencing different levels of drift and selection heterogeneity. We estimated the intensity and heterogeneity of selection on morphological and phenological traits via selection differentials. We demonstrate that Q(ST)-values increased with increasing selection heterogeneity when genetic drift was limited. The effect of selection on Q(ST) was thus detectable despite significant genotype-by-environment interactions that most probably biased the estimates of genetic differentiation. Although they cannot be used as a direct validation of the conclusions of prior studies, our results strongly support both the relevance of Q(ST) as an estimator of genetic differentiation and the role of local selection in shaping the genetic differentiation of natural populations.     

Interplay between extreme drift and selection intensities favors the fixation of beneficial mutations in selfing maize populations

Population and quantitative genetic models provide useful approximations to predict long-term selection responses sustaining phenotypic shifts, and underlying multilocus adaptive dynamics. Valid across a broad range of parameters, their use for understanding the adaptive dynamics of small selfing populations undergoing strong selection intensity (thereafter High Drift-High selection regime, HDHS) remains to be explored. Saclay Divergent Selection Experiments (DSEs) on maize flowering time provide an interesting example of populations evolving under HDHS, with significant selection responses over 20 generations in two directions. We combined experimental data from Saclay DSEs, forward individual-based simulations, and theoretical predictions to dissect the evolutionary mechanisms at play in the observed selection responses. We asked two main questions: How do mutations arise, spread, and reach fixation in populations evolving under HDHS? How does the interplay between drift and selection influence observed phenotypic shifts? We showed that the long-lasting response to selection in small populations is due to the rapid fixation of mutations occurring during the generations of selection. Among fixed mutations, we also found a clear signal of enrichment for beneficial mutations revealing a limited cost of selection. Both environmental stochasticity and variation in selection coefficients likely contributed to exacerbate mutational effects, thereby facilitating selection grasp and fixation of small-effect mutations. Together our results highlight that despite a small number of polymorphic loci expected under HDHS, adaptive variation is continuously fueled by a vast mutational target. We discuss our results in the context of breeding and long-term survival of small selfing populations.     

Rapid loss of MHC class II variation in a bottlenecked population is explained by drift and loss of copy number variation

Population bottlenecks may reduce genetic variation and potentially increase the risk of extinction. Here, we present the first study to use historic samples to analyse loss of variation at the major histocompatibility complex (MHC), which plays a central role in vertebrate disease resistance. Balancing selection acts on the MHC and could moderate the loss of variation expected from drift; however, in a Wisconsin population of greater prairie-chickens (Tympanuchus cupido), the number of MHC class II B alleles per individual declined by 44% following a population bottleneck, compared to a loss of only 8% at microsatellites. Simulations indicate that drift likely reduced MHC variation at the population level, as well as within individuals by reducing the number of gene copies per individual or by fixing the same alleles across multiple loci. These multiple effects of genetic drift on MHC variation could have important implications for immunity and fitness.     

Environmental variation, fluctuating selection and genetic drift in subdivided populations

Although there have many studies of the population genetical consequences of environmental variation, little is known about the combined effects of genetic drift and fluctuating selection in structured populations. Here we use diffusion theory to investigate the effects of temporally and spatially varying selection on a population of haploid individuals subdivided into a large number of demes. Using a perturbation method for processes with multiple time scales, we show that as the number of demes tends to infinity, the overall frequency converges to a diffusion process that is also the diffusion approximation for a finite, panmictic population subject to temporally fluctuating selection. We find that the coefficients of this process have a complicated dependence on deme size and migration rate, and that changes in these demographic parameters can determine both the balance between the dispersive and stabilizing effects of environmental variation and whether selection favors alleles with lower or higher fitness variance.     

Analysis of microsatellite variation in Pinus radiata reveals effects of genetic drift but no recent bottlenecks

Most conifer species occur in large continuous populations, but radiata pine, Pinus radiata, occurs only in five disjunctive natural populations in California and Mexico. The Mexican island populations were presumably colonized from the mainland millions of years ago. According to Axelrod (1981), the mainland populations are relicts of an earlier much wider distribution, reduced some 8,000 years ago, whereas according to Millar (1997, 2000), the patchy metapopulation-like structure is typical of the long-term population demography of the species. We used 19 highly polymorphic microsatellite loci to describe population structure and to search for signs of the dynamics of population demography over space and time. Frequencies of null alleles at microsatellite loci were estimated using an approach based on the probability of identity by descent. Microsatellite genetic diversities were high in all populations [expected heterozygosity (H(e)) = 0.68-0.77], but the island populations had significantly lower estimates. Variation between loci in genetic differentiation (F(ST)) was high, but no locus deviated statistically significantly from the rest at an experiment wide level of 0.05. Thus, all loci were included in subsequent analysis. The average differentiation was measured as F(ST) = 0.14 (SD 0.012), comparable with earlier allozyme results. The island populations were more diverged from the other populations and from an inferred common ancestral gene pool than the mainland ones. All populations showed a deficiency of expected heterozygosity given the number of alleles, the mainland populations more so than the island ones. The results thus do not support a recent important contraction in the mainland range of radiata pine.     

Clonality, genetic diversity and support for the diversifying selection hypothesis in natural populations of a flower-living yeast

Vast amounts of effort have been devoted to investigate patterns of genetic diversity and structuring in plants and animals, but similar information is scarce for organisms of other kingdoms. The study of the genetic structure of natural populations of wild yeasts can provide insights into the ecological and genetic correlates of clonality, and into the generality of recent hypotheses postulating that microbial populations lack the potential for genetic divergence and allopatric speciation. Ninety‐one isolates of the flower‐living yeast Metschnikowia gruessii from southeastern Spain were DNA fingerprinted using amplified fragment length polymorphism (AFLP) markers. Genetic diversity and structuring was investigated with band‐based methods and model‐ and nonmodel‐based clustering. Linkage disequilibrium tests were used to assess reproduction mode. Microsite‐dependent, diversifying selection was tested by comparing genetic characteristics of isolates from bumble bee vectors and different floral microsites. AFLP polymorphism (91%) and genotypic diversity were very high. Genetic diversity was spatially structured, as shown by amova (Φ_st = 0.155) and clustering. The null hypothesis of random mating was rejected, clonality seeming the prevailing reproductive mode in the populations studied. Genetic diversity of isolates declined from bumble bee mouthparts to floral microsites, and frequency of five AFLP markers varied significantly across floral microsites, thus supporting the hypothesis of diversifying selection on clonal lineages. Wild populations of clonal fungal microbes can exhibit levels of genetic diversity and spatial structuring that are not singularly different from those shown by sexually reproducing plants or animals. Microsite‐dependent, divergent selection can maintain high local and regional genetic diversity in microbial populations despite extensive clonality.     

Multiple parasites mediate balancing selection at two MHC class II genes in the fossorial water vole: insights from multivariate analyses and population genetics

We investigated the factors mediating selection acting on two MHC class II genes (DQA and DRB) in water vole (Arvicola scherman) natural populations in the French Jura Mountains. Population genetics showed significant homogeneity in allelic frequencies at the DQA1 locus as opposed to neutral markers (nine microsatellites), indicating balancing selection acting on this gene. Moreover, almost exhaustive screening for parasites, including gastrointestinal helminths, brain coccidia and antibodies against viruses responsible for zoonoses, was carried out. We applied a co-inertia approach to the genetic and parasitological data sets to avoid statistical problems related to multiple testing. Two alleles, Arte-DRB-11 and Arte-DRB-15, displayed antagonistic associations with the nematode Trichuris arvicolae, revealing the potential parasite-mediated selection acting on DRB locus. Selection mechanisms acting on the two MHC class II genes thus appeared different. Moreover, overdominance as balancing selection mechanism was showed highly unlikely in this system.     

Genetic diversity of MHC class I loci in six non-model frogs is shaped by positive selection and gene duplication

Comparative studies of major histocompatibility complex (MHC) genes across vertebrate species can reveal the evolutionary processes that shape the structure and function of immune regulatory proteins. In this study, we characterized MHC class I sequences from six frog species representing three anuran families (Hylidae, Centrolenidae and Ranidae). Using cDNA from our focal species, we amplified a total of 79 unique sequences spanning exons 2-4 that encode the extracellular domains of the functional alpha chain protein. We compared intra- and interspecific nucleotide and amino-acid divergence, tested for recombination, and identified codon sites under selection by estimating the rate of non-synonymous to synonymous substitutions with multiple codon-based maximum likelihood methods. We determined that positive (diversifying) selection was acting on specific amino-acid sites located within the domains that bind pathogen-derived peptides. We also found significant signals of recombination across the physical distance of the genes. Finally, we determined that all the six species expressed two or three putative classical class I loci, in contrast to the single locus condition of Xenopus laevis. Our results suggest that MHC evolution in anurans is a dynamic process and that variation in numbers of loci and genetic diversity can exist among taxa. Thus, the accumulation of genetic data for more species will be useful in further characterizing the relative importance of processes such as selection, recombination and gene duplication in shaping MHC loci among amphibian lineages.     

Interpreting population differentiation in terms of drift and selection

Summary Many empirical studies demonstrate some degree of genetic differentiation among populations of the same species. Understanding the relative importance of the processes causing this genetic differentiation has proven to be a difficult task. In particular, population differentiation can be influenced primarily by selection, genetic drift, and migration. We review the effect of drift and migration on patterns of genetic variation, with special reference to the conditions necessary for population differentiation. Conceptually, selection may be implicated in cases of population differentiation if the effect of drift and migration can be shown to be insufficient to cause the observed patterns. We examine some of the pitfalls of this approach when used with allozyme data, and revise a previous conclusion concerning the relative importance of selection in poulations of scale insects.