Heterozygote advantage at MHC DRB may influence response to infectious disease epizootics

The effect of MHC polymorphism on individual fitness variation in the wild remains equivocal; however, much evidence suggests that heterozygote advantage is a major determinant. To understand the contribution of MHC polymorphism to individual disease resistance or susceptibility in natural populations, we investigated two MHC class II B loci, DQB and DRB, in the New Zealand sea lion (NZSL, Phocarctos hookeri). The NZSL is a threatened species which is unusually susceptible to death by bacterial infection at an early age; it has suffered three bacterial induced epizootics resulting in high mortality levels of young pups since 1997. The MHC DQB and DRB haplotypes of dead NZSL pups with known cause of death (bacteria, enteritis or trauma) were sequenced and reconstructed, compared to pups that survived beyond 2 months of age, and distinct MHC DRB allele frequency and genotype differences were identified. Two findings were striking: (i) one DRB allele was present only in dead pups, and (ii) one heterozygous DRB genotype, common in live pups, was absent from dead pups. These results are consistent with some functional relationship with these variants and suggest heterozygote advantage is operating at DRB. We found no association between heterozygosity and fitness at 17 microsatellite loci, indicating that general heterozygosity is not responsible for the effect on fitness detected here. This result may be a consequence of recurrent selection by multiple pathogen assault over recent years and highlights the importance of heterozygote advantage at MHC as a potential mechanism for fitness differences in wild populations.     

Introgression from Domestic Goat Generated Variation at the Major Histocompatibility Complex of Alpine Ibex

The major histocompatibility complex (MHC) is a crucial component of the vertebrate immune system and shows extremely high levels of genetic polymorphism. The extraordinary genetic variation is thought to be ancient polymorphisms maintained by balancing selection. However, introgression from related species was recently proposed as an additional mechanism. Here we provide evidence for introgression at the MHC in Alpine ibex (Capra ibex ibex). At a usually very polymorphic MHC exon involved in pathogen recognition (DRB exon 2), Alpine ibex carried only two alleles. We found that one of these DRB alleles is identical to a DRB allele of domestic goats (Capra aegagrus hircus). We sequenced 2489 bp of the coding and non-coding regions of the DRB gene and found that Alpine ibex homozygous for the goat-type DRB exon 2 allele showed nearly identical sequences (99.8%) to a breed of domestic goats. Using Sanger and RAD sequencing, microsatellite and SNP chip data, we show that the chromosomal region containing the goat-type DRB allele has a signature of recent introgression in Alpine ibex. A region of approximately 750 kb including the DRB locus showed high rates of heterozygosity in individuals carrying one copy of the goat-type DRB allele. These individuals shared SNP alleles both with domestic goats and other Alpine ibex. In a survey of four Alpine ibex populations, we found that the region surrounding the DRB allele shows strong linkage disequilibria, strong sequence clustering and low diversity among haplotypes carrying the goat-type allele. Introgression at the MHC is likely adaptive and introgression critically increased MHC DRB diversity in the genetically impoverished Alpine ibex. Our finding contradicts the long-standing view that genetic variability at the MHC is solely a consequence of ancient trans-species polymorphism. Introgression is likely an underappreciated source of genetic diversity at the MHC and other loci under balancing selection.     

Gene duplications and sequence polymorphism of bovine class II DQB genes

The genetic diversity of bovine class II DQB genes was investigated by polymerase chain reaction amplification and DNA sequencing. The first domain exon was amplified from genomic DNA samples representing 14 class II haplotypes, defined by restriction fragment length polymorphism (RFLP) analysis. The presence of a polymorphism in the copy number of DQB genes was confirmed since two DQB sequences were isolated from certain haplotypes. Four subtypes of bovine DQB genes were found. DQB1 is the major type and was found in almost all haplotypes. DQB2 is very similar DQB1 but was found only in the duplicated haplotypes DQ9 to 12. DQB3 and DQB4 are two quite divergent genes only present in certain duplicated haplotypes. The bovine DQB complexity thus resembles that in the human DRB region. Bovine DQB genes were found to be highly polymorphic as ten DQB1 alleles and four DQB2 alleles were identified. The observed sequence polymorphism correlated well with previously defined DQB RFLPs. Bovine and human DQB alleles show striking similarities at the amino acid level. In contrast, the frequency of silent substitutions is much higher in comparisons of DQB alleles between species than within species ruling out the possibility that any of the contemporary DQB alleles have been maintained since the divergence of humans and cattle. The frequency of silent substitutions between DQB alleles was markedly lower in cattle than in humans, in agreement with a previous comparison of human and bovine DRB alleles.     

Diversity and duplication of DQB and DRB-like genes of the MHC in baleen whales (suborder: Mysticeti)

The molecular diversity and phylogenetic relationships of two class II genes of the baleen whale major histocompatibility complex were investigated and compared to toothed whales and out-groups. Amplification of the DQB exon 2 provided sequences showing high within-species and between-species nucleotide diversity and uninterrupted reading frames consistent with functional class II loci found in related mammals (e.g., ruminants). Cloning of amplified products indicated gene duplication in the humpback whale and triplication in the southern right whale, with average nucleotide diversity of 5.9 and 6.3%, respectively, for alleles of each species. Significantly higher nonsynonymous divergence at sites coding for peptide binding (32% for humpback and 40% for southern right) suggested that these loci were subject to positive (overdominant) selection. A population survey of humpback whales detected 23 alleles, differing by up to 21% of their inferred amino acid sequences. Amplification of the DRB exon 2 resulted in two groups of sequences. One was most similar to the DRB3 of the cow and present in all whales screened to date, including toothed whales. The second was most similar to the DRB2 of the cow and was found only in the bowhead and right whales. Both loci showed low diversity among species and apparent loss of function or altered function including interruption of reading frames. Finally, comparison of inferred protein sequence of the DRB3-like locus suggested convergence with the DQB, perhaps resulting from intergenic conversion or recombination.Electronic Supplementary Material Supplementary material is available for this article at     

Limited polymorphism at major histocompatibility complex (MHC) loci in the Swedish moose A. alces

The Swedish moose was analysed for genetic variability at major histocompatibility complex (MHC) class I and class II DQA, DQB and DRB loci using restriction fragment length polymorphism (RFLP) and single strand conformation polymorphism (SSCP) techniques. Both methods revealed limited amounts of polymorphism. Since the SSCP analysis concerned an expressed DRB gene it can be concluded that the level of functional MHC class II polymorphism, at least at the DRB locus, is low in Swedish moose. DNA fingerprinting was used to determine if the unusual pattern of low MHC variability could be explained by a low degree of genome-wide genetic diversity. Hybridizations with two minisatellite probes gave similarity indices somewhat higher than the average for other natural population, but the data suggest that the low MHC variability cannot be explained by a recent population bottleneck. However, since minisatellite sequences evolve more rapidly than MHC sequences, the low levels of MHC diversity may be attributed to a bottleneck of more ancient origin. The selection pressure for MHC variability in moose may also be reduced and we discuss the possibility that its solitary life style may reduce lateral transmission of pathogens in the population.     

Sequence Polymorphism and Evolution of Three Cetacean MHC Genes

Sequence variability at three major histocompatibility complex (MHC) genes (DQB, DRA, and MHC-I) of cetaceans was investigated in order to get an overall understanding of cetacean MHC evolution. Little sequence variation was detected at the DRA locus, while extensive and considerable variability were found at the MHC-I and DQB loci. Phylogenetic reconstruction and sequence comparison revealed extensive sharing of identical MHC alleles among different species at the three MHC loci examined. Comparisons of phylogenetic trees for these MHC loci with the trees reconstructed only based on non-PBR sites revealed that allelic similarity/identity possibly reflected common ancestry and were not due to adaptive convergence. At the same time, trans-species evolution was also evidenced that the allelic diversity of the three MHC loci clearly pre-dated species divergence events according to the relaxed molecular clock. It may be the forces of balancing selection acting to maintain the high sequence variability and identical alleles in trans-specific manner at the MHC-I and DQB loci.     

Characterization of Mhc-DRB allelic diversity in white-tailed deer (Odocoileus virginianus) provides insight into Mhc-DRB allelic evolution within Cervidae

 Although white-tailed deer (Odocoileus virginianus) are one of North America's best studied mammals, no information is available concerning allelic diversity at any locus of the major histocompatibility complex in this taxon. Using the polymerase chain reaction, single-stranded conformation polymorphism analysis, and DNA sequencing techniques, 15 DRB exon 2 alleles were identified among 150 white-tailed deer from a single population in southeastern Oklahoma. These alleles represent a single locus and exhibit a high degree of nucleotide and amino acid polymorphism, with most amino acid variation occurring at positions forming the peptide binding sites. Furthermore, twenty-seven amino acid residues unique to white-tailed deer DRB alleles were detected, with 19 of these occurring at residues forming contact points of the peptide binding region. Significantly higher rates of nonsynonymous than synonymous substitutions were detected among these DRB alleles. In contrast to other studies of Artiodactyla DRB sequences, interallelic recombination does not appear to be playing a significant role in the generation of allelic diversity at this locus in white-tailed deer. To examine evolution of white-tailed deer (Odvi-DRB) alleles within Cervidae, we performed a phylogenetic analysis of all published red deer (Ceel-DRB), roe deer (Caca-DRB), and moose (Alal-DRB) DRB alleles. The phylogenetic tree clearly shows a trans-species persistence of DRB lineages among these taxa. Moreover, this phylogenetic tree provides insight into evolution of DRB allelic lineages within Cervidae and may aid in assignment of red deer DRB alleles to specific loci. Received: 25 June 1998 / Revised: 2 September 1998     

Survey of major histocompatibility complex class II diversity in pig-tailed macaques

Pig-tailed macaques (Macaca nemestrina) serve as important models for human infectious disease research. Major histocompatibility complex (MHC) class II molecules are important to this research since they present peptides to CD4+ T cells. Despite the importance of characterizing the MHC-II alleles expressed in model species like pig-tailed macaques, to date, less than 150 MHC-II alleles have been named for the six most common classical class II loci (DRA, DRB, DQA, DQB, DPA, and DPB) in this population. Additionally, only a small percentage of these alleles are full-length, making it impossible to use the known sequence for reagent development. To address this, we developed a fast, high-throughput method to discover full-length MHC-II alleles and used it to characterize alleles in 32 pig-tailed macaques. By this method, we identified 128 total alleles across all six loci. We also performed an exon 2-based genotyping assay to validate the full-length sequencing results; this genotyping assay could be optimized for use in determining MHC-II allele frequencies in large cohorts of pig-tailed macaques.     

Evolution of MHC-<Emphasis Type="Italic">DRB</Emphasis> class II polymorphism in the genus <Emphasis Type="Italic">Apodemus</Emphasis> and a comparison of <Emphasis Type="Italic">DRB</Emphasis> sequences within the family Muridae (Mammalia: Rodentia)

Allelic diversity at major histocompatibility complex (MHC) genes is thought to be maintained by balancing selection over long periods of time, even across multiple speciation events. Trans-species sharing of MHC alleles among genera has been supported by many studies on mammals and fish, but in rodents, the results are ambiguous. We investigated natural levels of MHC-DRB variability and evolutionary processes in the wood mouse (Apodemus sylvaticus) and the yellow-necked mouse (Apodemus flavicollis), which are common, sympatric murid rodents in European forests. Using single-strand conformation polymorphism analysis and DNA sequencing, 38 DRB exon 2 alleles were detected among 162 A. sylvaticus from nine different locations in Germany and Switzerland, and 15 DRB exon 2 alleles were detected among 60 A. flavicollis from three different locations in northern Germany. There was evidence for balancing selection in both species. Phylogenetic analysis, including additional murid taxa, showed that the DRB exon 2 sequences did not separate according to species, consistent with trans-species evolution of the MHC in these taxa.     

Evolutionary genetics of MHC class II beta genes in the brown hare, <Emphasis Type="Italic">Lepus europaeus</Emphasis>

The genes of the major histocompatibility complex (MHC) are attractive candidates for investigating the link between adaptive variation and individual fitness. High levels of diversity at the MHC are thought to be the result of parasite-mediated selection and there is growing evidence to support this theory. Most studies, however, target just a single gene within the MHC and infer any evidence of selection to be representative of the entire gene region. Here we present data from three MHC class II beta genes (DPB, DQB, and DRB) for brown hares in two geographic regions and compare them against previous results from a class II alpha-chain gene (DQA). We report moderate levels of diversity and high levels of population differentiation in the DQB and DRB genes (Na = 11, D _est = 0.071 and Na = 15, D _est = 0.409, respectively), but not for the DPB gene (Na = 4, D _est = 0.00). We also detected evidence of positive selection within the peptide binding region of the DQB and DRB genes (95% CI, ω > 1.0) but found no signature of selection for DPB. Mutation and recombination were both found to be important processes shaping the evolution of the class II genes. Our findings suggest that while diversifying selection is a significant contributor to the generally high levels of MHC diversity, it does not act in a uniform manner across the entire MHC class II region. The beta-chain genes that we have characterized provide a valuable set of MHC class II markers for future studies of the evolution of adaptive variation in Leporids.     

Genetic diversity of MHC class II <Emphasis Type="Italic">DRB</Emphasis> alleles in the continental and Japanese populations of the sable <Emphasis Type="Italic">Martes zibellina</Emphasis> (Mustelidae,Carnivora, Mammalia)

The sable (Martes zibellina) is a medium-sized mustelid inhabiting forest environments in Siberia, northern China, the Korean Peninsula, and Hokkaido Island, Japan. To further understand the molecular evolution of the major histocompatibility complex (MHC), we sequenced part of exon 2 in MHC class II DRB genes, including codons encoding the antigen binding site, from 33 individuals from continental Eurasia and Japan. We identified 16 MHC class II DRB alleles (Mazi-DRBs), some of which were geographically restricted and others broadly distributed, and eight putative pseudogenes. A single-breakpoint recombination analysis detected a recombination site in the middle of exon 2. A mixed effects model of evolution analysis identified five amino acid sites presumably under positive selection. These sites were all located in the region 3′ to the recombination site, suggesting that positive selection and recombination could be committed to the diversity of the M. zibellina DRB gene. In a Bayesian phylogenetic tree, all Mazi-DRBs and the presumed pseudogenes grouped within a Mustelidae clade. The Mazi-DRBs showed trans-species polymorphism, with some alleles most closely related to alleles from other mustelid species. This result suggests that the sable DRBs have evolved under long-lasting balancing selection.     

Recombination and the origin of sequence diversity in the DRB MHC class II locus in chamois (Rupicapra spp.)

We examined the evolutionary processes contributing to genetic diversity at the major histocompatibility complex (MHC) class II DRB locus in chamois (Rupicapra spp., subfamily Caprinae). We characterised the pattern of intragenic recombination (or homologous gene conversion) and quantified the amount of recombination in the genealogical history of the two chamois species, Pyrenean chamois (Rupicapra pyrenaica) and Alpine chamois (Rupicapra rupicapra). We found evidence for intragenic recombination, and the estimated amount of population recombination suggests that recombination has been a significant process in generating DRB allelic diversity in the genealogical history of the genus Rupicapra. Moreover, positive selection appears to act on the same peptide-binding residues in both analysed chamois species, but not in identical intensity. Recombination coupled with positive selection drives the rapid evolution at the peptide-binding sites in the MHC class II DRB gene. Many chamois MHC class II DRB alleles are thus much younger than previously assumed.     

High MHC DQB Variation and Asymmetric Allelic Distribution in the Endangered Yangtze Finless Porpoise, Neophocaena phocaenoides asiaeorientalis

The endangered Yangtze finless porpoise is found in the middle and lower reaches of the Yangtze River and its adjoining big lakes. To explore the major histocompatibility complex (MHC) genetic diversity and allelic distribution patterns across its range, we investigated variation at DQB exon 2. From 76 porpoises, we identified 18 DQB sequences. The freshwater Yangtze populations had much higher allelic diversity than marine populations. Among these freshwater populations, the middle-reach population had higher allelic diversity than the lower-reach population. The high DQB diversity level, relative to that of a neutral mtDNA locus, suggests that balancing selection is acting at the DQB gene and that rapid evolution and local positive selection play critical roles in generating and retaining high MHC diversity in the freshwater population. As the balancing selection might be driven by environmental pathogens, we suggest that maintaining MHC variation should be a high priority in the conservation and management of this endangered population, especially as an ex situ conservation strategy.     

Low MHC DRB class II diversity in the mountain goat: past bottlenecks and possible role of pathogens and parasites

Major histocompatibility complex (MHC) genes are the most polymorphic in vertebrates and code for molecules playing a central role in pathogen resistance. We studied levels of MHC DRB class II diversity in a long-term study population of mountain goats (Oreamnos americanus) at Caw Ridge, Alberta, and two other populations from British Columbia, Canada. Only two alleles were found among the three populations sampled. The Caw Ridge population was fixed for one of the two MHC DRB alleles, but this lack of variation did not appear to have affected it negatively because the population doubled over two decades and had no history of any apparent infectious diseases. Past population bottlenecks during Pleistocene glaciations are thought to have been the main factor contributing to the low levels of MHC diversity in mountain goats, a hypothesis supported by our previous work reporting low polymorphism at neutral loci. Additionally, the limited MHC variability in mountain goats may be related to its northern distribution as we found that allelic diversity at MHC DRB class II in wild ungulates decreases with increasing latitude, possibly as a result of low parasite diversity at high latitudes. The low MHC variation in mountain goats and other northern ungulates such as muskoxen (Ovibos moschatus) may expose these species to population outbreaks that could be generated by introduced pathogens or northward shifts in the distribution of pathogens with global climate warming.     

Genetic Variation at Exon 2 of the MHC Class II DQB Locus in Blue Whale (Balaenoptera musculus) from the Gulf of California

The genes of the Major Histocompatibility Complex (MHC) play an important role in the vertebrate immune response and are among the most polymorphic genes known in vertebrates. In some marine mammals, MHC genes have been shown to be characterized by low levels of polymorphism compared to terrestrial taxa; this reduction in variation is often explained as a result of lower pathogen pressures in marine habitats. To determine if this same reduction in variation applies to the migratory population of blue whales (Balaenoptera musculus) that occurs in the Gulf of California, we genotyped a 172 bp fragment of exon 2 of the MHC Class II DQB locus for 80 members of this population. Twenty-two putatively functional DQB allotypes were identified, all of which were homologous with DQB sequences from other cetacean species. Up to 5 putative alleles per individual were identified, suggesting that gene duplication has occurred at this locus. Rates of non-synonymous to synonymous substitutions (ω) and maximum likelihood analyses of models of nucleotide variation provided potential evidence of ongoing positive selection at this exon. Phylogenetic analyses of DQB alleles from B. musculus and 16 other species of cetaceans revealed trans-specific conservation of MHC variants, suggesting that selection has acted on this locus over prolonged periods of time. Collectively our findings reveal that immunogenic variation in blue whales is comparable to that in terrestrial mammals, thereby providing no evidence that marine taxa are subject to reduced pathogen-induced selective pressures.