Variation on islands: major histocompatibility complex (Mhc) polymorphism in populations of the Australian bush rat

The role of demographic history and selection in shaping genetic diversity of the Galápagos penguin (Spheniscus mendiculus)

Although many studies have documented the effects of demographic bottlenecks on the genetic diversity of natural populations, there is conflicting evidence of the roles that genetic drift and selection may play in driving changes in genetic variation at adaptive loci. We analyzed genetic variation at microsatellite and mitochondrial loci in conjunction with an adaptive MHC class II locus in the Galápagos penguin (Spheniscus mendiculus), a species that has undergone serial demographic bottlenecks associated with El Niño events through its evolutionary history. We compared levels of variation in the Galápagos penguin to those of its congener, the Magellanic penguin (Spheniscus magellanicus), which has consistently maintained a large population size and thus was used as a non-bottlenecked control. The comparison of neutral and adaptive markers in these two demographically distinct species allowed assessment of the potential role of balancing selection in maintaining levels of MHC variation during bottleneck events. Our analysis suggests that the lack of genetic diversity at both neutral and adaptive loci in the Galápagos penguin likely resulted from its restricted range, relatively low abundance, and history of demographic bottlenecks. The Galápagos penguin revealed two MHC alleles, one mitochondrial haplotype, and six alleles across five microsatellite loci, which represents only a small fraction of the diversity detected in Magellanic penguins. Despite the decreased genetic diversity in the Galápagos penguin, results revealed signals of balancing selection at the MHC, which suggest that selection can mitigate some of the effects of genetic drift during bottleneck events. Although Galápagos penguin populations have persisted for a long time, increased frequency of El Niño events due to global climate change, as well as the low diversity exhibited at immunological loci, may put this species at further risk of extinction.

Variation in Major Histocompatibility Complex diversity in invasive cane toad populations

Context The cane toad (Rhinella marina), a native species of central and southern America, was introduced to Australia in 1935 as a biocontrol agent after a complex history of prior introductions. The population rapidly expanded and has since spread through much of the Australian landmass, with severe impacts on the endemic wildlife, primarily via toxicity to predators. The invasion process has taken its toll on the cane toad, with changes in the immunological capacity across the Australian invasive population. Aims To investigate the immunogenetic underpinnings of these changes, we studied the diversity of the Major Histocompatiblity Complex (MHC) genes in introduced cane toad populations. Methods We studied the diversity of two MHC genes (the classical class I UA locus and a class II DAB locus) and compared these with neutral microsatellite markers in toads from the Australian site of introduction and the Australian invasion front. We also included toads from Hawai’i, the original source of the Australian toads, to infer founder effect. Key results Diversity across all markers was low across Australian and Hawai’ian samples, consistent with a reduction in genetic diversity through multiple founder effects during the course of the successive translocations. In Australia, allelic diversity at the microsatellite markers and the UA locus was reduced at the invasion front, whereas all three alleles at the DAB locus were maintained in the invasion-front toads. Conclusions Loss of allelic diversity observed at the microsatellite markers and the UA locus could be the result of drift and bottlenecking along the invasion process, however, the persistence of DAB diversity warrants further investigation to disentangle the evolutionary forces influencing this locus. Implications Through the use of different molecular markers, we provide a preliminary description of the adaptive genetic processes occurring in this invasive population. The extremely limited MHC diversity may represent low immunogenetic competence across the Australian population, which could be exploited for invasive species management.

Genetic diversity comparison of the DQA gene in European rabbit (Oryctolagus cuniculus) populations

The European rabbit (Oryctolagus cuniculus) natural populations within the species native region, the Iberian Peninsula, are considered a reservoir of genetic diversity. Indeed, the Iberia was a Pleistocene refuge to the species and currently two subspecies are found in the peninsula (Oryctolagus cuniculus cuniculus and Oryctolagus cuniculus algirus). The genes of the major histocompatibility complex (MHC) have been substantially studied in wild populations due to their exceptional variability, believed to be pathogen driven. They play an important function as part of the adaptive immune system affecting the individual fitness and population viability. In this study, the MHC variability was assessed by analysing the exon 2 of the DQA gene in several European rabbit populations from Portugal, Spain and France and in domestic breeds. Twenty-eight DQA alleles were detected, among which 18 are described for the first time. The Iberian rabbit populations are well differentiated from the French population and domestic breeds. The Iberian populations retained the higher allelic diversity with the domestic breeds harbouring the lowest; in contrast, the DQA nucleotide diversity was higher in the French population. Signatures of positive selection were detected in four codons which are putative peptide-binding sites and have been previously detected in other mammals. The evolutionary relationships showed instances of trans-species polymorphism. Overall, our results suggest that the DQA in European rabbits is evolving under selection and genetic drift.

MHC class II diversity of koala (Phascolarctos cinereus) populations across their range

Major histocompatibility complex class II (MHCII) genes code for proteins that bind and present antigenic peptides and trigger the adaptive immune response. We present a broad geographical study of MHCII DA β1 (DAB) and DB β1 (DBB) variants of the koala (Phascolarctos cinereus; n=191) from 12 populations across eastern Australia, with a total of 13 DAB and 7 DBB variants found. We identified greater MHCII variation and, possibly, additional gene copies in koala populations in the north (Queensland and New South Wales) relative to the south (Victoria), confirmed by STRUCTURE analyses and genetic differentiation using analysis of molecular variance. The higher MHCII diversity in the north relative to south could potentially be attributed to (i) significant founder effect in Victorian populations linked to historical translocation of bottlenecked koala populations and (ii) increased pathogen-driven balancing selection and/or local genetic drift in the north. Low MHCII genetic diversity in koalas from the south could reduce their potential response to disease, although the three DAB variants found in the south had substantial sequence divergence between variants. This study assessing MHCII diversity in the koala with historical translocations in some populations contributes to understanding the effects of population translocations on functional genetic diversity.

Major histocompatibility complex alleles associated with parasite susceptibility in wild giant pandas

Major histocompatibility complex (MHC) polymorphism is thought to be driven by antagonistic coevolution between pathogens and hosts, mediated through either overdominance or frequency-dependent selection. However, investigations under natural conditions are still rare for endangered mammals which often exhibit depleted variation, and the mechanism of selection underlying the maintenance of characteristics remains a considerable debate. In this study, 87 wild giant pandas were used to investigate MHC variation associated with parasite load. With the knowledge of the MHC profile provided by the genomic data of the giant panda, seven DRB1, seven DQA1 and eight DQA2 alleles were identified at each single locus. Positive selection evidenced by a significantly higher number of non-synonymous substitutions per non-synonymous codon site relative to synonymous substitutions per synonymous codon site could only be detected at the DRB1 locus, which leads to the speculation that DRB1 may have a more important role in dealing with parasite infection for pandas. Coprological analyses revealed that 55.17% of individuals exhibited infection with 1-2 helminthes and 95.3% of infected pandas carried Baylisascaris shroederi. Using a generalized linear model, we found that Aime-DRB1*10 was significantly associated with parasite infection, but no resistant alleles could be detected. MHC heterozygosity of the pandas was found to be uncorrelated with the infection status or the infection intensity. These results suggested that the possible selection mechanisms in extant wild pandas may be frequency dependent rather than being determined by overdominance selection. Our findings could guide the candidate selection for the ongoing reintroduction or translocation of pandas.

How closely does genetic diversity in finite populations conform to predictions of neutral theory? Large deficits in regions of low recombination

Levels of genetic diversity in finite populations are crucial in conservation and evolutionary biology. Genetic diversity is required for populations to evolve and its loss is related to inbreeding in random mating populations, and thus to reduced population fitness and increased extinction risk. Neutral theory is widely used to predict levels of genetic diversity. I review levels of genetic diversity in finite populations in relation to predictions of neutral theory. Positive associations between genetic diversity and population size, as predicted by neutral theory, are observed for microsatellites, allozymes, quantitative genetic variation and usually for mitochondrial DNA (mtDNA). However, there are frequently significant deviations from neutral theory owing to indirect selection at linked loci caused by balancing selection, selective sweeps and background selection. Substantially lower genetic diversity than predicted under neutrality was found for chromosomes with low recombination rates and high linkage disequilibrium (compared with 'normally' recombining chromosomes within species and adjusted for different copy numbers and mutation rates), including W (median 100% lower) and Y (89% lower) chromosomes, dot fourth chromosomes in Drosophila (94% lower) and mtDNA (67% lower). Further, microsatellite genetic and allelic diversity were lost at 12 and 33% faster rates than expected in populations adapting to captivity, owing to widespread selective sweeps. Overall, neither neutral theory nor most versions of the genetic draft hypothesis are compatible with all empirical results.

Selection maintains MHC diversity through a natural population bottleneck

A perceived consequence of a population bottleneck is the erosion of genetic diversity and concomitant reduction in individual fitness and evolutionary potential. Although reduced genetic variation associated with demographic perturbation has been amply demonstrated for neutral molecular markers, the effective management of genetic resources in natural populations is hindered by a lack of understanding of how adaptive genetic variation will respond to population fluctuations, given these are affected by selection as well as drift. Here, we demonstrate that selection counters drift to maintain polymorphism at a major histocompatibility complex (MHC) locus through a population bottleneck in an inbred island population of water voles. Before and after the bottleneck, MHC allele frequencies were close to balancing selection equilibrium but became skewed by drift when the population size was critically low. MHC heterozygosity generally conformed to Hardy-Weinberg expectations except in one generation during the population recovery where there was a significant excess of heterozygous genotypes, which simulations ascribed to strong differential MHC-dependent survival. Low allelic diversity and highly skewed frequency distributions at microsatellite loci indicated potent genetic drift due to a strong founder affect and/or previous population bottlenecks. This study is a real-time examination of the predictions of fundamental evolutionary theory in low genetic diversity situations. The findings highlight that conservation efforts to maintain the genetic health and evolutionary potential of natural populations should consider the genetic basis for fitness-related traits, and how such adaptive genetic diversity will vary in response to both the demographic fluctuations and the effects of selection.

Impact of historical founder effects and a recent bottleneck on MHC variability in Commander Arctic foxes (Vulpes lagopus)

Populations of Arctic foxes (Vulpes lagopus) have been isolated on two of the Commander Islands (Bering and Mednyi) from the circumpolar distributed mainland population since the Pleistocene. In 1970-1980, an epizootic outbreak of mange caused a severe population decline on Mednyi Island. Genes of the major histocompatibility complex (MHC) play a primary role in infectious disease resistance. The main objectives of our study were to compare contemporary variation of MHC class II in mainland and island Arctic foxes, and to document the effects of the isolation and the recent bottleneck on MHC polymorphism by analyzing samples from historical and contemporary Arctic foxes. In 184 individuals, we found 25 unique MHC class II DRB and DQB alleles, and identified evidence of balancing selection maintaining allelic lineages over time at both loci. Twenty different MHC alleles were observed in mainland foxes and eight in Bering Island foxes. The historical Mednyi population contained five alleles and all contemporary individuals were monomorphic at both DRB and DQB. Our data indicate that despite positive and diversifying selection leading to elevated rates of amino acid replacement in functionally important antigen-binding sites, below a certain population size, balancing selection may not be strong enough to maintain genetic diversity in functionally important genes. This may have important fitness consequences and might explain the high pathogen susceptibility in some island populations. This is the first study that compares MHC diversity before and after a bottleneck in a wild canid population using DNA from museum samples.

Major Histocompatibility Complex (MHC) markers in conservation biology

Human impacts through habitat destruction, introduction of invasive species and climate change are increasing the number of species threatened with extinction. Decreases in population size simultaneously lead to reductions in genetic diversity, ultimately reducing the ability of populations to adapt to a changing environment. In this way, loss of genetic polymorphism is linked with extinction risk. Recent advances in sequencing technologies mean that obtaining measures of genetic diversity at functionally important genes is within reach for conservation programs. A key region of the genome that should be targeted for population genetic studies is the Major Histocompatibility Complex (MHC). MHC genes, found in all jawed vertebrates, are the most polymorphic genes in vertebrate genomes. They play key roles in immune function via immune-recognition and -surveillance and host-parasite interaction. Therefore, measuring levels of polymorphism at these genes can provide indirect measures of the immunological fitness of populations. The MHC has also been linked with mate-choice and pregnancy outcomes and has application for improving mating success in captive breeding programs. The recent discovery that genetic diversity at MHC genes may protect against the spread of contagious cancers provides an added impetus for managing and protecting MHC diversity in wild populations. Here we review the field and focus on the successful applications of MHC-typing for conservation management. We emphasize the importance of using MHC markers when planning and executing wildlife rescue and conservation programs but stress that this should not be done to the detriment of genome-wide diversity.

Genetic impacts of Anacapa deer mice reintroductions following rat eradication

The Anacapa deer mouse is an endemic subspecies that inhabits Anacapa Island, part of Channel Islands National Park, California. We used mitochondrial DNA cytochrome c oxidase subunit II gene (COII) and 10 microsatellite loci to evaluate the levels of genetic differentiation and variation in ~1400 Anacapa deer mice sampled before and for 4 years after a black rat (Rattus rattus) eradication campaign that included trapping, captive holding and reintroduction of deer mice. Both mitochondrial and microsatellite analyses indicated significant differentiation between Anacapa deer mice and mainland mice, and genetic variability of mainland mice was significantly higher than Anacapa mice even prior to reintroduction. Bayesian cluster analysis and Principal Coordinates Analysis indicated that East, Middle and West Anacapa mice were genetically differentiated from each other, but translocation of mice among islands resulted in the East population becoming less distinct as a result of management. Levels of heterozygosity were similar before and after management. However, numerous private alleles in the founder populations were not observed after reintroduction and shifts in allele frequencies occurred, indicating that the reintroduced populations experienced substantial genetic drift. Surprisingly, two mitochondrial haplotypes observed in an earlier study of Anacapa deer mice were lost in the 20 years prior to the rat eradication program, leaving only a single haplotype in Anacapa deer mice. This study demonstrates how genetic monitoring can help to understand the re-establishment of endemic species after the eradication of invasive species and to evaluate the effectiveness of the management strategies employed.

Reduced MHC and neutral variation in the Galápagos hawk, an island endemic

BACKGROUND

Genes at the major histocompatibility complex (MHC) are known for high levels of polymorphism maintained by balancing selection. In small or bottlenecked populations, however, genetic drift may be strong enough to overwhelm the effect of balancing selection, resulting in reduced MHC variability. In this study we investigated MHC evolution in two recently diverged bird species: the endemic Galápagos hawk (Buteo galapagoensis), which occurs in small, isolated island populations, and its widespread mainland relative, the Swainson's hawk (B. swainsoni).

RESULTS

We amplified at least two MHC class II B gene copies in each species. We recovered only three different sequences from 32 Galápagos hawks, while we amplified 20 unique sequences in 20 Swainson's hawks. Most of the sequences clustered into two groups in a phylogenetic network, with one group likely representing pseudogenes or nonclassical loci. Neutral genetic diversity at 17 microsatellite loci was also reduced in the Galápagos hawk compared to the Swainson's hawk.

CONCLUSIONS

The corresponding loss in neutral diversity suggests that the reduced variability present at Galápagos hawk MHC class II B genes compared to the Swainson's hawk is primarily due to a founder event followed by ongoing genetic drift in small populations. However, purifying selection could also explain the low number of MHC alleles present. This lack of variation at genes involved in the adaptive immune response could be cause for concern should novel diseases reach the archipelago.

Contrasting responses to selection in class I and class IIα major histocompatibility-linked markers in salmon

Comparison of levels and patterns of genetic variation in natural populations either across loci or against neutral expectation can yield insight into locus-specific differences in the strength and direction of evolutionary forces. We used both approaches to test the hypotheses on patterns of selection on major histocompatibility (MH)-linked markers. We performed temporal analyses of class I and class IIα MH-linked markers and eight microsatellite loci in two Atlantic salmon populations in Ireland on two temporal scales: over six decades and 9 years in the rivers Burrishoole and Delphi, respectively. We also compared contemporary Burrishoole and Delphi samples with nearby populations for the same loci. On comparing patterns of temporal and spatial differentiation among classes of loci, the class IIα MH-linked marker was consistently identified as an outlier compared with patterns at the other microsatellite loci or neutral expectation. We found higher levels of temporal and spatial heterogeneity in heterozygosity (but not in allelic richness) for the class IIα MH-linked marker compared with microsatellites. Tests on both within- and among-population differentiation are consistent with directional selection acting on the class IIα-linked marker in both temporal and spatial comparisons, but only in temporal comparisons for the class I-linked marker. Our results indicate a complex pattern of selection on MH-linked markers in natural populations of Atlantic salmon. These findings highlight the importance of considering selection on MH-linked markers when using these markers for management and conservation purposes.

The impact of isolation and bottlenecks on genetic diversity in the Pearson Island population of the black-footed rock-wallaby (Petrogale lateralis pearsoni; Marsupialia:Macropodidae)

A morphologically distinct subspecies of black-footed rock-wallaby (Petrogale lateralis pearsoni), is naturally found only on North Pearson Island (~160 ha) in the Investigator Group, SA, where it was isolated by rising sea levels ~10 500 years ago. Subsequent translocations have seen additional populations established on Middle-South Pearson Island (~53 ha) in 1960 and Wedge Island (96 ha) in 1975. We have used 10 hypervariable microsatellite loci to examine the levels of genetic diversity in the endemic (n = 38) and translocated (n = 45–77) P. l. pearsoni populations compared with mainland P. lateralis populations (n = 19–52). Results show that all sampled P. l. pearsoni populations have very low levels of genetic diversity (A = 1.5–1.9; HE = 0.02–0.13) compared with mainland populations (A = 3.5–12.7; HE = 0.54–0.87). Intriguingly, more diversity was detected in the translocated Middle-South Pearson population than in its source population from North Pearson Island. In contrast, the Wedge Island population was almost monomorphic. Overall, the severe loss of genetic diversity (up to 98%) in P. l. pearsoni populations appears to result from random genetic drift on a small isolated population, exacerbated by some subsequent one-off translocation events. Although additional supplementary translocations are recommended to enhance genetic diversity, populations of P. l. pearsoni are likely to remain inherently vulnerable to extinction and therefore of special conservation concern.

Allelic diversity of the MHC class II DRB genes in brown bears (Ursus arctos) and a comparison of DRB sequences within the family Ursidae

The allelic diversity of the DRB locus in major histocompatibility complex (MHC) genes was analyzed in the brown bear (Ursus arctos) from the Hokkaido Island of Japan, Siberia, and Kodiak of Alaska. Nineteen alleles of the DRB exon 2 were identified from a total of 38 individuals of U. arctos and were highly polymorphic. Comparisons of non-synonymous and synonymous substitutions in the antigen-binding sites of deduced amino acid sequences indicated evidence for balancing selection on the bear DRB locus. The phylogenetic analysis of the DRB alleles among three genera (Ursus, Tremarctos, and Ailuropoda) in the family Ursidae revealed that DRB allelic lineages were not separated according to species. This strongly shows trans-species persistence of DRB alleles within the Ursidae.

Major histocompatibility complex heterozygote advantage and widespread bacterial infections in populations of Chinook salmon (Oncorhynchus tshawytscha)

Despite growing evidence for parasite-mediated selection on the vertebrate major histocompatibility complex (MHC), little is known about variation in the bacterial parasite community within and among host populations or its influence on MHC evolution. In this study, we characterize variation in the parasitic bacterial community associated with Chinook salmon (Oncorhynchus tshawytscha) fry in five populations in British Columbia (BC), Canada across 2 years, and examine whether bacterial infections are a potential source of selection on the MHC. We found an unprecedented diversity of bacteria infecting fry with a total of 55 unique bacteria identified. Bacterial infection rates varied from 9% to 29% among populations and there was a significant isolation by distance relationship in bacterial community phylogenetic similarity across the populations. Spatial variation in the frequency of infections and in the phylogenetic similarity of bacterial communities may result in differential parasite-mediated selection at the MHC across populations. Across all populations, we found evidence of a heterozygote advantage at the MHC class II, which may be a source of balancing selection on this locus. Interestingly, a co-inertia analysis indicated only susceptibility associations between a few of the MHC class I and II alleles and specific bacterial parasites; there was no evidence that any of the alleles provided resistance to the bacteria. Our results reveal a complex bacterial community infecting populations of a fish and underscore the importance of considering the role of multiple pathogens in the evolution of host adaptations.

MHC genetic structure and divergence across populations of Chinook salmon (Oncorhynchus tshawytscha)

The major histocompatibility complex (MHC) is thought to be under strong selection pressure because of its integral role in pathogen recognition. Consequently, patterns of MHC genetic variation should reflect selection pressures across the landscape. We examined genetic variation and population genetic structure at the MHC class I-A1 and class II-B1 exons in five Chinook salmon (Oncorhynchus tshawytscha) populations from two geographic regions in British Columbia, Canada. We then compared estimates of population structure at the MHC genes with neutral estimates based on microsatellites to examine the potential for local adaptation at the MHC. Chinook salmon are in decline throughout much of their native range and understanding the degree of local adaptation exhibited by the MHC may be important in conservation planning. Comparisons among populations yielded higher G'(ST) estimates for the MHC class I than expected under neutrality based on the microsatellites. In contrast, the MHC class II tended to exhibit lower G'(ST) values than did the microsatellites. These results suggest that across populations unique selection pressures are driving allele frequency differences at the MHC class I but that the MHC class II may be the subject of homogenizing selection. Rates of nonsynonymous versus synonymous substitutions found in codons associated within the MHC class I and II peptide-binding regions provided strong evidence of positive selection. Together, these results support the hypothesis that selection is influencing genetic variation at the MHC, but suggest that selection pressures may vary at the two classes of loci both at the sequence and population levels.

Tail shedding in island lizards [Lacertidae, Reptilia]: decline of antipredator defenses in relaxed predation environments

The ability of an animal to shed its tail is a widespread antipredator strategy among lizards. The degree of expression of this defense is expected to be shaped by prevailing environmental conditions including local predation pressure. We test these hypotheses by comparing several aspects of caudal autotomy in 15 Mediterranean lizard taxa existing across a swath of mainland and island localities that differ in the number and identity of predator species present. Autotomic ease varied substantially among the study populations, in a pattern that is best explained by the presence of vipers. Neither insularity nor the presence of other types of predators explain the observed autotomy rates. Final concentration of accumulated tail muscle lactate and duration of movement of a shed tail, two traits that were previously thought to relate to predation pressure, are in general not shaped by either predator diversity or insularity. Under conditions of relaxed predation selection, an uncoupling of different aspects of caudal autotomy exists, with some elements (ease of autotomy) declining faster than others (duration of movement, lactate concentration). We compared rates of shed tails in the field against rates of laboratory autotomies conducted under standardized conditions and found very high correlation values (r > 0.96). This suggests that field autotomy rates, rather than being a metric of predatory attacks, merely reflect the innate predisposition of a taxon to shed its tail.

Balancing and directional selection at exon-2 of the MHC DQB1 locus among populations of odontocete cetaceans

The diversity of exon-2 (peptide-binding region) of the DQB1 locus (Class II, major histocompatibility complex, MHC) was investigated on an extended sample of populations of three focal cetacean species (two sibling delphinid species and another in the same family). We tested the hypothesis that dolphin populations with a worldwide distribution across different habitats and geographic regions will be under differential selective pressure by comparing DQB1 variation with variation at neutral markers and by investigating putative functional residues within the exon-2 sequence at the population level. Variation at the DQB1 locus was not correlated to neutral differentiation (assessed by comparison with microsatellite DNA markers), and overall F(ST) values were significantly lower for the MHC locus, consistent with expectations for balancing selection. Measures of heterozygosity and d(n)/d(s) ratios were also consistent with balancing selection. However, outliers in the F(ST) comparisons and the analysis of putative functional residues suggested incidences of directional selection in local populations.

Evolutionary history of an MHC gene in two leporid species: characterisation of Mhc-DQA in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning-sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism.

Diversity and evolutionary history of the MHC DQA gene in leporids

The European rabbit (Oryctolagus cuniculus) is used as a model for many human diseases, yet comparatively little is known of its genetics, particularly at important loci such as the major histocompatibility complex (MHC). This study investigated genetic diversity and evolutionary history of the DQA gene in a range of leporid species by analysing coding sequence diversity of exon 2 and intron 2 in 53 individuals of 16 different species. Fifty leporid DQA alleles were detected, including 13 novel European rabbit alleles. In the rabbit, the highest levels of diversity were observed in wild rabbits from Portugal, with wild rabbits from England and domestic rabbits showing less diversity. Within the sample, several recombination events were detected and trans-specific evolution of alleles was evidenced, both being general characteristics of mammalian MHC genes. Positive selection is implicated as operating on six codons within exon 2, which are also subject to positive selection in other mammals. Some of these positions are putative antigen recognition sites and underline the importance of pathogen-driven selection on these MHC genes.

Genetic drift vs. natural selection in a long-term small isolated population: major histocompatibility complex class II variation in the Gulf of California endemic porpoise (Phocoena sinus)

Although many studies confirm long-term small isolated populations (e.g. island endemics) commonly sustain low neutral genetic variation as a result of genetic drift, it is less clear how selection on adaptive or detrimental genes interplay with random forces. We investigated sequence variation at two major histocompatibility complex (Mhc) class II loci on a porpoise endemic to the upper Gulf of California, México (Phocoena sinus, or vaquita). Its unique declining population is estimated around 500 individuals. Single-strand conformation polymorphism analysis revealed one putative functional allele fixed at the locus DQB (n = 25). At the DRB locus, we found two presumed functional alleles (n = 29), differing by a single nonsynonymous nucleotide substitution that could increase the stability at the dimer interface of alphabeta-heterodimers on heterozygous individuals. Identical trans-specific DQB1 and DRB1 alleles were identified between P. sinus and its closest relative, the Burmeister's porpoise (Phocoena spinipinnis). Comparison with studies on four island endemic mammals suggests fixation of one allele, due to genetic drift, commonly occurs at the DQA or DQB loci (effectively neutral). Similarly, deleterious alleles of small effect are also effectively neutral and can become fixed; a high frequency of anatomical malformations on vaquita gave empirical support to this prediction. In contrast, retention of low but functional polymorphism at the DRB locus was consistent with higher selection intensity. These observations indicated natural selection could maintain (and likely also purge) some crucial alleles even in the face of strong and prolonged genetic drift and inbreeding, suggesting long-term small populations should display low inbreeding depression. Low levels of Mhc variation warn about a high susceptibility to novel pathogens and diseases in vaquita.

Low MHC variation in the endangered Galápagos penguin (Spheniscus mendiculus)

The major histocompatibility complex (MHC) is one of the most polymorphic regions of the genome, likely due to balancing selection acting to maintain alleles over time. Lack of MHC variability has been attributed to factors such as genetic drift in small populations and relaxed selection pressure. The Galápagos penguin (Spheniscus mendiculus), endemic to the Galápagos Islands, is the only penguin that occurs on the equator. It relies upon cold, nutrient-rich upwellings and experiences severe population declines when ocean temperatures rise during El Niño events. These bottlenecks, occurring in an already small population, have likely resulted in reduced genetic diversity in this species. In this study, we used MHC class II exon 2 sequence data from a DRB1-like gene to characterize the amount of genetic variation at the MHC in 30 Galápagos penguins, as well as one Magellanic penguin (S. magellanicus) and two king penguins (Aptenodytes patagonicus), and compared it to that in five other penguin species for which published data exist. We found that the Galápagos penguin had the lowest MHC diversity (as measured by number of polymorphic sites and average divergence among alleles) of the eight penguin species studied. A phylogenetic analysis showed that Galápagos penguin MHC sequences are most closely related to Humboldt penguin (Spheniscus humboldti) sequences, its putative sister species based on other loci. An excess of non-synonymous mutations and a pattern of trans-specific evolution in the neighbor-joining tree suggest that selection is acting on the penguin MHC.

Duplication, balancing selection and trans-species evolution explain the high levels of polymorphism of the DQA MHC class II gene in voles (Arvicolinae)

Major histocompatibility complex (MHC) genes play important role in host-parasite interactions and parasites are crucial factors influencing the population dynamics of hosts. We described the structure and diversity of exon 2 of the MHC class II DQA gene in three species of voles (Arvicolinae) exhibiting regular multi-annual fluctuations of population density and analysed the processes leading to the observed MHC polymorphism. By using cloning-sequencing methodology and capillary electrophoresis-single strand conformation polymorphism, we described seven sequences in the water, eight in the common, and seven in the bank voles coming from an area of 70 km(2) around the Nozeroy canton in the Jura Mountains (Franche Comté, France). All exon 2 sequences translate to give unique amino acid sequences and positive selection was found to act very intensively on antigen binding sites. We documented the presence of recombination at vole DQA region but its importance in generating allelic polymorphism seems to be relatively limited. For the first time within rodents, we documented the duplication of the DQA gene in all three species with both copies being transcriptionally active. Phylogenetic analysis of allelic sequences revealed extensive trans-species polymorphism within the subfamily although no alleles were shared between species in our data set. We discuss possible role of parasites in forming the recent polymorphism pattern of the DQA locus in voles.

Genetic Variation of the MHCDQBLocus 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.

MHC gene configuration variation in geographically disparate populations of California sea lions (Zalophus californianus)

Major histocompatibility complex (MHC) class II DRB genotypes were examined in two geographically isolated populations of California sea lions (Zalophus californianus) (Gulf of California and California coastal Pacific Ocean). Genomic DNA from 227 California sea lions was examined using eight sequence-specific primer (SSP) pairs flanking the putative peptide-binding site. A total of 40 different Zaca-DRB genotype configurations were identified among the 227 individuals. Using SSP-PCR, significant differences were found between coastal California and Gulf of California Zalophus populations in numbers of DRB sequences per individual and configuration of sequences within individuals. Additionally, unique local patterns of MHC diversity were identified among the Midriff Island animals. These population differences are consistent with either ecologically distinct patterns of selection pressures and/or geographical isolation. The consequences of these partitioned MHC configurations at the population level are as yet unknown; however, the worldwide increase in emerging marine diseases lends urgency to their examination.

Sequence diversity of the MHC DRB gene in the Eurasian beaver (Castor fiber)

Major histocompatibility complex (MHC) genes, coding molecules which play an important role in immune response, are the most polymorphic genes known in vertebrates. However, MHC polymorphism in some species is limited. MHC monomorphism at several MHC class I and II loci was previously reported for two neighbouring northern European populations of the Eurasian beaver (Castor fiber) and reduced selection for polymorphism has been hypothesized. Here, we analysed a partial sequence of the second exon of the MHC II DRB locus from seven relict European and Asian beaver populations. We detected 10 unique alleles among 76 beavers analysed. Only a western Siberian population was polymorphic, with four alleles detected in 10 individuals. Each of the remaining populations was fixed for a different allele. Sequences showed considerable divergence, suggesting the long persistence of allelic lineages. A significant excess of nonsynonymous substitutions was detected at the antigen binding sites, indicating that sequence evolution of beaver DRB was driven by positive selection. Current MHC monomorphism in the majority of populations may be the result of the superimposition of the recent bottleneck on pre-existing genetic structure resulting from population subdivision and differential pathogen pressure.

The importance of immune gene variability (MHC) in evolutionary ecology and conservation

Genetic studies have typically inferred the effects of human impact by documenting patterns of genetic differentiation and levels of genetic diversity among potentially isolated populations using selective neutral markers such as mitochondrial control region sequences, microsatellites or single nucleotide polymorphism (SNPs). However, evolutionary relevant and adaptive processes within and between populations can only be reflected by coding genes. In vertebrates, growing evidence suggests that genetic diversity is particularly important at the level of the major histocompatibility complex (MHC). MHC variants influence many important biological traits, including immune recognition, susceptibility to infectious and autoimmune diseases, individual odours, mating preferences, kin recognition, cooperation and pregnancy outcome. These diverse functions and characteristics place genes of the MHC among the best candidates for studies of mechanisms and significance of molecular adaptation in vertebrates. MHC variability is believed to be maintained by pathogen-driven selection, mediated either through heterozygote advantage or frequency-dependent selection. Up to now, most of our knowledge has derived from studies in humans or from model organisms under experimental, laboratory conditions. Empirical support for selective mechanisms in free-ranging animal populations in their natural environment is rare. In this review, I first introduce general information about the structure and function of MHC genes, as well as current hypotheses and concepts concerning the role of selection in the maintenance of MHC polymorphism. The evolutionary forces acting on the genetic diversity in coding and non-coding markers are compared. Then, I summarise empirical support for the functional importance of MHC variability in parasite resistance with emphasis on the evidence derived from free-ranging animal populations investigated in their natural habitat. Finally, I discuss the importance of adaptive genetic variability with respect to human impact and conservation, and implications for future studies.

Major histocompatibility complex and microsatellite variation in two populations of wild gorillas

In comparison to their close relatives the chimpanzees and humans, very little is known concerning the amount and structure of genetic variation in gorillas. Two species of gorillas are recognized and while the western gorillas number in the tens of thousands, only several hundred representatives of the mountain gorilla subspecies of eastern gorillas survive. To analyse the possible effects of these different population sizes, this study compares the variation observed at microsatellite and major histocompatibility complex (MHC) loci in samples of wild western and mountain gorillas, collected using a sampling scheme that targeted multiple social groups within defined geographical areas. Noninvasive samples proved a viable source of DNA for sequence analysis of the second exon of the DRB loci of the MHC. Observed levels of variation at the MHC locus were similar between the two gorilla species and were comparable to those in other primates. Comparison of results from analysis of variation at multiple microsatellite loci found only a slight reduction in heterozygosity for the mountain gorillas despite the relatively smaller population size.

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.