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.     

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.     

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.     

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.     

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.     

Parasite-mediated evolution of the functional part of the MHC in primates

The major histocompatibility complex (MHC) is a key model of genetic polymorphism, but the mechanisms underlying its extreme variability are debated. Most hypotheses for MHC diversity focus on pathogen-driven selection and predict that MHC polymorphism evolves under the pressure of a diverse parasite fauna. Several studies reported that certain alleles offer protection against certain parasites, yet it remains unclear whether variation in parasite pressure more generally covaries with allelic diversity and rates of molecular evolution of MHC across species. We tested this prediction in a comparative study of 41 primate species. We characterized polymorphism of the exon 2 of DRB region of the MHC class II. Our phylogenetic analyses controlled for the potential effects of neutral mutation rate, population size, geographic origin and body mass and revealed that nematode species richness associates positively with nonsynonymous nucleotide substitution rate at the functional part of the molecule. We failed to find evidence for allelic diversity being strongly related to parasite species richness. Continental distribution was a strong predictor of both allelic diversity and substitution rate, with higher values in Malagasy and Neotropical primates. These results indicate that parasite pressure can influence the different estimates of MHC polymorphism, whereas geography plays an independent role in the natural history of MHC.     

Stabilizing selection maintains exuberant colour polymorphism in the spider Theridion californicum (Araneae, Theridiidae)

Genetically controlled colour polymorphisms provide a physical manifestation of the operation of selection and how this can vary according to the spatial or temporal arrangement of phenotypes, or their frequency in a population. Here, we examine the role of selection in shaping the exuberant colour polymorphism exhibited by the spider Theridion californicum. This species is part of a system in which several distantly related spiders in the same lineage, but living in very different geographical areas, exhibit remarkably convergent polymorphisms. These polymorphisms are characterized by allelic inheritance and the presence of a single common cryptic morph and, in the case of T. californicum and its congener the Hawaiian happy-face spider Theridion grallator, numerous rare patterned morphs. We compare population differentiation estimated from colour phenotypic data to differentiation at neutral amplified fragment length polymorphisms (AFLP) loci and demonstrate that the colour polymorphism appears to be maintained by balancing selection. We also examine the patterns of selection in the genome-wide sample of AFLP loci and compare approaches to detecting signatures of selection in this context. Our results have important implications regarding balancing selection, suggesting that selective agents can act in a similar manner across disparate taxa in globally disjunct locales resulting in parallel evolution of exuberant polymorphism.     

Low genetic diversity and local adaptive divergence of Dracaena cambodiana (Liliaceae) populations associated with historical population bottlenecks and natural selection: an endangered long‐lived tree endemic to Hainan Island,China

Historical population bottlenecks and natural selection have important effects on the current genetic diversity and structure of long‐lived trees. Dracaena cambodiana is an endangered, long‐lived tree endemic to Hainan Island, China. Our field investigations showed that only 10 populations remain on Hainan Island and that almost all have been seriously isolated and grow in distinct habitats. A considerable amount of genetic variation at the species level, but little variation at the population level, and a high level of genetic differentiation among the populations with limited gene flow in D. cambodiana were detected using inter‐simple sequence repeat (ISSR) and random amplified polymorphic DNA (RAPD) analyses. No significant correlation was found between genetic diversity and actual population size, as the genetic diversities were similar regardless of population size. The Mantel test revealed that there was no correlation between genetic and geographic distances among the 10 populations. The UPGMA, PCoA and Bayesian analyses showed that local adaptive divergence has occurred among the D. cambodiana populations, which was further supported by habitat‐private fragments. We suggest that the current genetic diversity and population differentiation of D. cambodiana resulted from historical population bottlenecks and natural selection followed by historical isolation. However, the lack of natural regeneration of D. cambodiana indicates that former local adaptations with low genetic diversity may have been genetically weak and are unable to adapt to the current ecological environments.     

Microsatellite and morphometric variation in Drosophila gouveai: the relative importance of historical and current factors in shaping the genetic population structure

Drosophila gouveai is a cactus-breeding species with a naturally fragmented distribution in central-western and southeastern Brazil. In this study, a neutral genetic marker (microsatellite DNA) and a quantitative trait (wing morphology) were used to investigate the population structure of eight populations of D. gouveai . Quantitative analysis was done using a morphometric approach and multivariate analysis of 17 wing parameters. Drosophila gouveai showed a large degree of genetic population differentiation ( F _ST = 0.245). A molecular analysis of variance showed that the population structure was attributable mainly to genetic differences among the population groups, which were correlated with regional groupings. A Bayesian cluster analysis also identified the same regional groupings as contributing to the population structure. Assignment tests revealed that the current gene flow was very restricted. The divergence in wing morphology among populations was also high, and revealed a geographical pattern that conformed to a latitudinal cline. In contrast to current factors such as migration-drift equilibrium, these results indicated that the genetic population structure of D. gouveai was shaped predominantly by historical factors.     

Loss of MHC and neutral variation in Peary caribou: genetic drift is not mitigated by balancing selection or exacerbated by MHC allele distributions

Theory and empirical results suggest that the rate of loss of variation at Mhc and neutral microsatellite loci may differ because selection influences Mhc genes, and because a high proportion of rare alleles at Mhc loci may result in high rates of loss via drift. Most published studies compare Mhc and microsatellite variation in various contemporary populations to infer the effects of population size on genetic variation, even though different populations are likely to have different demographic histories that may also affect contemporary genetic variation. We directly compared loss of variation at Mhc and microsatellite loci in Peary caribou by comparing historical and contemporary samples. We observed that similar proportions of genetic variation were lost over time at each type of marker despite strong evidence for selection at Mhc genes. These results suggest that microsatellites can be used to estimate genome-wide levels of variation, but also that adaptive potential is likely to be lost following population bottlenecks. However, gene conversion and recombination at Mhc loci may act to increase variation following bottlenecks.     

Extreme environmental variation sharpens selection that drives the evolution of a mutualism

The importance of infrequent events for both adaptive evolution and the evolution of species interactions is largely unknown. We investigated how the infrequent production of large seed crops (masting) of a bird-dispersed tree (whitebark pine, Pinus albicaulis) influenced phenotypic selection exerted by its primary avian seed predator-disperser, the Clark's nutcracker (Nucifraga columbiana). Selection was not evident during common years of low seed abundance, whereas it was replicated among areas and favoured traits facilitating seed dispersal during infrequent years of high seed abundance. Since nutcrackers act mostly as seed predators during small seed crops but as seed dispersers during the largest seed crops, trees experienced strong selection from nutcrackers only during infrequent years when the interaction was most strongly mutualistic. Infrequent events can thus be essential to both adaptive evolution and the evolutionary dynamics of species interactions.     

Colonizing the world in spite of reduced MHC variation

The major histocompatibility complex (MHC), which harbours the most polymorphic vertebrate genes, plays a critical role in the host-pathogen coevolutionary arms race. However, the extent to which MHC diversity determines disease susceptibility and long-term persistence of populations is currently under debate, as recent studies have demonstrated that low MHC variability does not necessarily hamper population viability. However, these studies typically assayed small and decimated populations in species with restricted distribution, thereby making inferences about the evolutionary potential of these populations difficult. Here, we show that MHC impoverishment has not constrained the ecological radiation and flourishing of falcons (Aves: Falconidae) worldwide. We found two remarkably different patterns of MHC variation within the genus Falco. Whereas MHC variation in kestrels (the basal group within the genus) is very high, falcons exhibit ancestrally low intra- and interspecific MHC variability. This pattern is not due to the inadvertent survey of paralogous genes or pseudogenes. Further, patterns of variation in mitochondrial or other nuclear genes do not indicate a generalized low level of genome-wide variability among falcons. Although a relative contribution of genetic drift cannot be completely ruled out, we propose the falcons went through an evolutionary transition, driven and maintained by natural selection, from primarily highly variable towards low polymorphic and slow-evolving MHC genes with a very specific immune function. This study highlights that the importance of MHC diversity cannot be generalized among vertebrates, and hints at the evolution of compensatory immune mechanisms in falcons to cope with emerging and continuously evolving pathogens.