Molecular characterization of major histocompatibility complex class II alleles in the common frog, Rana temporaria

Major histocompatibility complex (MHC) class II genes, which play a major role in the immune system response, are some of the most polymorphic genes in vertebrates. We developed polymerase chain reaction primers for part of the second exon of an expressed MHC class II gene in the common frog, Rana temporaria. We genotyped this locus in five frog populations in southeast England and detected eight alleles in 215 individuals. Five or six alleles were detected in each population with a maximum of two alleles per individual, indicating that only a single locus was amplified. We also inferred the possible existence of a null allele. There were 23 variable nucleotide sites (out of 136) and 13 variable amino acid sites (out of 44), many of which corresponded to amino acids involved in antigen recognition. We detected a significant excess of nonsynonymous substitutions at antigen binding sites, indicating that this gene is under positive selection. The level of variation we found was similar to that in other amphibian MHC class II loci, such as those in Bombina bombina, Xenopus laevis and Ambystoma tigrinum.     

Genetic Variation of the Major Histocompatibility Complex (MHC Class II B Gene) in the Threatened Hume’s Pheasant,Syrmaticus humiae

Major histocompatibility complex (MHC) genes are the most polymorphic genes in vertebrates and encode molecules that play a crucial role in pathogen resistance. As a result of their diversity, they have received much attention in the fields of evolutionary and conservation biology. Here, we described the genetic variation of MHC class II B (MHCIIB) exon 2 in a wild population of Hume’s pheasant (Syrmaticus humiae), which has suffered a dramatic decline in population over the last three decades across its ranges in the face of heavy exploitation and habitat loss. Twenty-four distinct alleles were found in 73 S. humiae specimens. We found seven shared alleles among four geographical groups as well as six rare MHCIIB alleles. Most individuals displayed between one to five alleles, suggesting that there are at least three MHCIIB loci of the Hume’s pheasant. The d _N ⁄ d _S ratio at putative antigen-binding sites (ABS) was significantly greater than one, indicating balancing selection is acting on MHCIIB exon 2. Additionally, recombination and gene conversion contributed to generating MHCIIB diversity in the Hume’s pheasant. One to three recombination events and seventy-five significant gene conversion events were observed within the Hume’s pheasant MHCIIB loci. The phylogenetic tree and network analysis revealed that the Hume’s pheasant alleles do not cluster together, but are scattered through the tree or network indicating a trans-species evolutionary mode. These findings revealed the evolution of the Hume’s pheasant MHC after suffering extreme habitat fragmentation.     

Genetic variability of MHC class II DQB exon 2 alleles in yak (Bos grunniens)

The major histocompatibility class (MHC) DQ molecules are dimeric glycoproteins revealing antigen presentation to CD⁴⁺ T cells. In the present study, the exon 2 of the MHC class II DQB gene from 32 yaks (Bos grunniens) was cloned, sequenced and compared with previously reported patterns for other bovidae. It was revealed by sequence analyses that there are 25 DQB exon 2 alleles among 32 yaks, all alleles are found to belong to DQB1 loci. These alleles exhibited a high degree of nucleotide and amino acid polymorphisms with most amino acid variations occurring at positions forming the peptide-binding sites. The DQB loci were analyzed for patterns of synonymous (d _S) and non-synonymous (d _N) substitution. The yak was observed to be under strong positive selection in the DQB exon 2 peptide-binding sites (d _N = 0.15, P < 0.001). It appears that this variability among yaks confers the ability to mount immune responses to a wide variety of peptides or pathogens.     

MHC influences infection with parasites and winter survival in the root vole Microtus oeconomus

Selective pressure from parasites is thought to maintain the polymorphism of major histocompatibility complex (MHC) genes. Although a number of studies have shown a relationship between the MHC and parasitic infections, the fitness consequences of such associations are less well documented. In the present paper, we characterised the variation in exon 2 of MHC class II DRB gene in the root vole and examined the effects of that gene on parasite prevalence and winter survival. We identified 18 unique exon 2 sequences, which translated into 10 unique amino acid sequences. Phylogenetic analysis revealed the presence of three distinct clusters, and allele distributions among these individuals suggested that the clusters correspond to three different loci. Although the rate of synonymous substitutions (d_S) exceeded the rate of nonsynonymous substitutions (d_N) across sequences, implying purifying selection, d_N was significantly elevated at antigen-binding sites, suggesting that these sites could be under positive selection. Screening for parasites revealed a moderate prevalence of infection with gastrointestinal parasites (24 % infected), but a high infection rate for blood parasites (56 % infected). Infection with the blood parasite Babesia ssp. decreased survival almost twofold (25.7 vs. 13.9 %). Animals possessing the amino acid sequence AA*08 survived better than others (44.9 vs. 22 %), and they were infected with Babesia ssp. less often (13.9 vs 25.7 %). In contrast, individuals carrying allele AA*05 were infected more often (31.7 vs. 15.3 %). Heterozygosity at one of the putative loci was associated with a lower probability of infection with Babesia ssp., but at the other locus, the association was reversed. The unexpected latter result could be at least partly explained by the increased frequency of the susceptible allele AA*05 among heterozygotes. Overall, we demonstrate that infection with Babesia ssp. is a strong predictor of winter survival and that MHC genes are important predictors of infection status as well as survival in the root vole.     

The transcription of MGAT4A glycosyl transferase is increased in white cells of peripheral blood of Type 2 Diabetes patients

Background

The giant panda (Ailuropoda melanoleuca) is one of the most endangered animals due to habitat fragmentation and loss. Although the captive breeding program for this species is now nearly two decades old, researches on the genetic background of such captive populations, especially on adaptive molecular polymorphism of major histocompatibility complex (MHC), are still limited. In this study, we characterized adaptive variation of the giant panda's MHC DQA gene by PCR amplification of its antigen-recognizing region (i.e. the exon 2) and subsequent single-strand conformational polymorphism (SSCP) and sequence analyses.

Results

The results revealed a low level of DQA exon 2 diversity in this rare animal, presenting 6 alleles from 61 giant panda individuals. The observed polymorphism was restricted to 9 amino acid substitutions, all of which occurred at and adjacent to positions forming the functionally important antigen-binding sites. All the samples were in Hardy-Weinberg proportions. A significantly higher rate of non-synonymous than synonymous substitutions at the antigen-binding sites indicated positive selection for diversity in the locus.

Conclusion

The DQA allelic diversity of giant pandas was low relative to other vertebrates. Nonetheless, the pandas exhibited more alleles in DQA than those in DRB, suggesting the alpha chain genes would play a leading role when coping with certain pathogens and thus should be included in conservation genetic investigation. The microsatellite and MHC loci might predict long-term persistence potential and short-term survival ability, respectively. Consequently, it is recommended to utilize multiple suites of microsatellite markers and multiple MHC loci to detect overall genetic variation in order to design unbiased conservation strategies.     

Evidence of gene orthology and trans‐species polymorphism,but not of parallel evolution,despite high levels of concerted evolution in the major histocompatibility complex of flamingo species

The major histocompatibility complex (MHC) is a cornerstone in the study of adaptive genetic diversity. Intriguingly, highly polymorphic MHC sequences are often not more similar within species than between closely related species. Divergent selection of gene duplicates, balancing selection maintaining trans‐species polymorphism (TSP) that predate speciation and parallel evolution of species sharing similar selection pressures can all lead to higher sequence similarity between species. In contrast, high rates of concerted evolution increase sequence similarity of duplicated loci within species. Assessing these evolutionary models remains difficult as relatedness and ecological similarities are often confounded. As sympatric species of flamingos are more distantly related than allopatric species, flamingos represent an ideal model to disentangle these evolutionary models. We characterized MHC Class I exon 3, Class IIB exon 2 and exon 3 of the six extant flamingo species. We found up to six MHC Class I loci and two MHC Class IIB loci. As all six species shared the same number of MHC Class IIB loci, duplication appears to predate flamingo speciation. However, the high rate of concerted evolution has prevented the divergence of duplicated loci. We found high sequence similarity between all species regardless of codon position. The latter is consistent with balancing selection maintaining TSP, as under this mechanism amino acid sites under pathogen‐mediated selection should be characterized by fewer synonymous codons (due to their common ancestry) than under parallel evolution. Overall, balancing selection maintaining TSP appears to result in high MHC similarity between species regardless of species relatedness and geographical distribution.     

Genetic variation and selection of MHC class I loci differ in two congeneric frogs

Major histocompatibility complex (MHC) genes encode proteins in the acquired immune response pathway that often show distinctive selection-driven patterns in wild vertebrate populations. We examined genetic variation and signatures of selection in the MHC class I alpha 1 (A1)- and alpha 2 (A2)-domain encoding exons of two frog congeners [Agalychnis callidryas (n?=?20) and A. lemur (n?=?20)] from a single locality in Panama. We also investigated how historical demographic processes may have impacted MHC genetic diversity by analyzing a neutral mitochondrial marker. We found that both MHC domains were highly variable in both species, with both species likely expressing three loci. Our analyses revealed different signatures of selection between the two species, most notably that the A. callidryas A2 domain had experienced positive selection while the A2 domain of A. lemur had not. Diversifying selection acted on the same number of A1 and A2 allelic lineages, but on a higher percentage of A1 sites compared to A2 sites. Neutrality tests of mitochondrial haplotypes predominately indicated that the two species were at genetic equilibrium when the samples were collected. In addition, two historical tests of demography indicated both species have had relatively stable population sizes over the past 100,000 years; thus large population size changes are unlikely to have greatly influenced MHC diversity in either species during this time period. In conclusion, our results suggest that the impact of selection on MHC diversity varied between these two closely related species, likely due to a combination of distinct ecological conditions and past pathogenic pressures.     

Population genetic diversity and geographical differentiation of MHC class II DAB genes in the vulnerable Chinese egret (<Emphasis Type="Italic">Egretta eulophotes</Emphasis>)

Major histocompatibility complex (MHC) genes are excellent markers for the study of adaptive genetic variation occurring over different geographical scales. The Chinese egret (Egretta eulophotes) is a vulnerable ardeid species with an estimated global population of 2600–3400 individuals. In this study, we sampled 172 individuals of this egret (approximately 6 % of the global population) from five natural populations that span the entire distribution range of this species in China. We examined their population genetic diversity and geographical differentiation at three MHC class II DAB genes by identifying eight exon 2 alleles at Egeu-DAB1, eight at Egeu-DAB2 and four at Egeu-DAB3. Allelic distributions at each of these three Egeu-DAB loci varied substantially within the five populations, while levels of genetic diversity varied slightly among the populations. Analysis of molecular variance showed low but significant genetic differentiation among five populations at all three Egeu-DAB loci (haplotype-based ?_ST: 0.029, 0.020 and 0.042; and distance-based ?_ST: 0.036, 0.027 and 0.043, respectively; all P < 0.01). The Mantel test suggested that this significant population genetic differentiation was likely due to an isolation-by-distance pattern of MHC evolution. However, the phylogenetic analyses and the Bayesian clustering analysis based on the three Egeu-DAB loci indicated that there was little geographical structuring of the genetic differentiation among five populations. These results provide fundamental population information for the conservation genetics of the vulnerable Chinese egret.     

MHC class I typing in a songbird with numerous loci and high polymorphism using motif-specific PCR and DGGE

The major histocompatibility complex (MHC) has a central role in the specific immune defence of vertebrates. Exon 3 of MHC class I genes encodes the domain that binds and presents peptides from pathogens that trigger immune reactions. Here we develop a fast population screening method for detecting genetic variation in the MHC class I genes of birds. We found evidence of at least 15 exon 3 sequences in the investigated great reed warbler individual. The organisation of the great reed warbler MHC class I genes suggested that a locus-specific screening protocol is impractical due to the high similarity between alleles across loci, including the introns flanking exon 3. Therefore, we used motif-specific PCR to amplify two subsets of alleles (exon 3 sequences) that were separated with by DGGE. The motif-specific primers amplify a substantial proportion of the transcribed class I alleles (2-12 alleles per individual) from as many as six class I loci. Although not exhaustive, this gives a reliable estimate of the class I variation. The method is highly repeatable and more sensitive in detecting genetic variation than the RFLP method. The motif-specific primers also allow us to avoid screening pseudogenes. In our study population of great reed warblers, we found a high level of genetic variation in MHC class I, and no less than 234 DGGE genotypes were detected among 248 screened individuals.     

凹耳蛙MHCII类B基因第二外显子多态性分析

两栖类正经历全球范围内的种群衰退,很多两栖动物集群灭绝事件与环境病原体(如壶菌(Batrachochytrium dendrobatidis)的侵扰有关。MHC基因的表达产物在有颌脊椎动物免疫应答过程中起关键作用,其多态性通常与动物对疾病的抗性或易感性密切相关,因而被认为是研究动物适应性进化的最佳候选基因之一。本文对中国特有的无尾两栖动物凹耳蛙(Odorrana tormota)MHC II类B基因多态性进行初步研究。首先,利用1对通用引物扩增出凹耳蛙MHC II类B基因exon2长约180bp的DNA片段。在此基础上,利用ligation-mediated PCR进一步获取侧翼未知序列,序列拼接后长2,030bp,包含exon2以及intron1和intron2的部分序列。基于上述序列设计出凹耳蛙B基因exon2特异性引物(IIQ1BU/IIQ1BD),对该物种黄山种群32个样品进行PCR扩增和克隆测序,共获得34个不同的等位基因,等位基因序列核苷酸和氨基酸变异位点的比例分别为16.17%(33/204)和26.87%(18/67),大多数氨基酸变异位点位于推测的抗原结合位点(antigen binding sites,ABS)。每个样品包含2-5个等位基因,结合等位基因序列特征以及cDNA表达分析结果,推测凹耳蛙至少拥有3个可表达的B基因座位。与文献报道的蛙科其他物种比较后发现,尽管凹耳蛙目前的分布区非常狭窄,但其MHC II类B基因多态性明显高于蛙科其他动物。等位基因碱基替换模式提示凹耳蛙MHC II类B基因曾经历过强烈的正选择作用,ABS区的dN值显著大于dS(P<0.05),PAML软件包CODEML程序中不同模型的似然比检测(likelihood rate test)结果同样支持上述推论,贝叶斯经验贝叶斯路径(Bayesian Em-pirical Bayes)共检测出5个显著受正选择作用的氨基酸位点。贝叶斯系统树的拓扑结构显示,无尾两栖类不同科的等位基因分别形成单系群,但蛙科不同属的等位基因未能形成单系群,蛙属绿池蛙(Rana clamitans)的1个等位基因与臭蛙属凹耳蛙的部分等位基因享有共同的谱系关系,提示蛙科不同属间的B基因存在跨种多态性。     

The evolution of the major histocompatibility complex in upstream versus downstream river populations of the longnose dace

Populations in upstream versus downstream river locations can be exposed to vastly different environmental and ecological conditions and can thus harbor different genetic resources due to selection and neutral processes. An interesting question is how upstream–downstream directionality in rivers affects the evolution of immune response genes. We used next‐generation amplicon sequencing to identify eight alleles of the major histocompatibility complex (MHC) class II β exon 2 in the cyprinid longnose dace (Rhinichthys cataractae) from three rivers in Alberta, upstream and downstream of municipal and agricultural areas along contaminant gradients. We used these data to test for directional and balancing selection on the MHC. We also genotyped microsatellite loci to examine neutral population processes in this system. We found evidence for balancing selection on the MHC in the form of increased nonsynonymous variation relative to neutral expectations, and selection occurred at more amino acid residues upstream than downstream in two rivers. We found this pattern despite no population structure or isolation by distance, based on microsatellite data, at these sites. Overall, our results suggest that MHC evolution is driven by upstream–downstream directionality in fish inhabiting this system.     

High adaptive variability and virus-driven selection on major histocompatibility complex (MHC) genes in invasive wild rabbits in Australia

The rabbit haemorrhagic disease virus (RHDV) was imported into Australia in 1995 as a biocontrol agent to manage one of the most successful and devastating invasive species, the European rabbit (Oryctolagus cuniculus cuniculus). During the first disease outbreaks, RHDV caused mortality rates of up to 97% and reduced Australian rabbit numbers to very low levels. However, recently increased genetic resistance to RHDV and strong population growth has been reported. Major histocompatibility complex (MHC) class I immune genes are important for immune responses against viruses, and a high MHC variability is thought to be crucial in adaptive processes under pathogen-driven selection. We asked whether strong population bottlenecks and presumed genetic drift would have led to low MHC variability in wild Australian rabbits, and if the retained MHC variability was enough to explain the increased resistance against RHD. Despite the past bottlenecks we found a relatively high number of MHC class I sequences distributed over 2–4 loci. We identified positive selection on putative antigen-binding sites of the MHC. We detected evidence for RHDV-driven selection as one MHC supertype was negatively associated with RHD survival, fitting expectations of frequency-dependent selection. Gene duplication and pathogen-driven selection are possible (and likely) mechanisms that maintained the adaptive potential of MHC genes in Australian rabbits. Our findings not only contribute to a better understanding of the evolution of invasive species, they are also important in the light of planned future rabbit biocontrol in Australia.     

尼罗鳄MHCⅡ类B基因第二外元的序列变异分析

利用一对简并引物扩增了尼罗鳄MHCⅡ类分子B基因第二外元的部分片段,并对PCR产物进行了克隆和测序,结果得到8种不同的序列,序列长度为166 bp;经分析,序列中有56个变异位点,核苷酸的非同义替换多于同义替换,造成30个位点氨基酸的改变,氨基酸的替换趋于集中在假定的抗原结合位点附近.核苷酸和氨基酸序列与已报道的扬子鳄和密河鳄的MHCⅡ类B基因第二外元序列有较高的同源性,利用PAUP4.0软件构建的NJ树显示,鳄类的MHCⅡ类B基因存在跨种多态性现象.     

The salmonid MHC class I: more ancient loci uncovered

An unprecedented level of sequence diversity has been maintained in the salmonid major histocompatibility complex (MHC) class I UBA gene, with between lineage AA sequence identities as low as 34%. The derivation of deep allelic lineages may have occurred through interlocus exon shuffling or convergence of ancient loci with the UBA locus, but until recently, no such ancient loci were uncovered. Herein, we document the existence of eight additional MHC class I loci in salmon (UCA, UDA, UEA, UFA, UGA, UHA, ULA, and ZE), six of which share exon 2 and 3 lineages with UBA, and three of which have not been described elsewhere. Half of the UBA exon 2 lineages and all UBA exon 3 lineages are shared with other loci. Two loci, UGA and UEA, share only a single exon lineage with UBA, likely generated through exon shuffling. Based on sequence homologies, we hypothesize that most exchanges and duplications occurred before or during tetraploidization (50 to 100 Ma). Novel loci that share no relationship with other salmonid loci are also identified (UHA and ZE). Each locus is evaluated for its potential to function as a class Ia gene based on gene expression, conserved residues and polymorphism. UBA is the only locus that can indisputably be classified as a class Ia gene, although three of the eight loci (ZE, UCA, and ULA) conform in three out of four measures. We hypothesize that these additional loci are in varying states of degradation to class Ib genes.     

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.     

MHC class I variation in a natural blue tit population (Cyanistes caeruleus)

The major histocompatibility complex (MHC) is central to the vertebrate immune system and its highly polymorphic genes are considered to influence several life-history traits of individuals. To characterize the MHC in a natural population of blue tits (Cyanistes caeruleus) we investigated the class I exon 3 diversity of more than 900 individuals. We designed two pairs of motif-specific primers that reliably amplify independent subsets of MHC alleles. Applying denaturing gradient gel electrophoresis (DGGE) we obtained 48 independently inherited units of unique band patterns (DGGE-haplogroups), which were validated in a segregation analysis within 105 families. In a second approach, we extensively sequenced 6 unrelated individuals to confirm that DGGE-haplogroup composition reflects individual allelic variation. The highest number of different DGGE-haplogroups in a single individual corresponded in 19 MHC exon 3 sequences, suggesting a minimum of 10 amplified MHC class I loci in the blue tit. In total, we identified 50 unique functional and 3 non-functional sequences. Functional sequences showed high levels of recombination and strong positive selection in the antigen binding region, whereas nucleotide diversity was comparatively low in the range of all passerine species. Finally, in a phylogenetic comparison of passerine MHC class I exon 3 sequences we discuss conflicting evolutionary signals possibly due to recent gene duplication, recombination events and concerted evolution. Our results indicate that the described method is suitable to effectively explore the MHC diversity and its ecological impacts in blue tits in future studies.     

Structure,diversity and evolutionary patterns of expressed MHC class II<Emphasis Type="Italic">B</Emphasis> genes in chub (<Emphasis Type="Italic">Squalius cephalus</Emphasis>), a cyprinid fish species from Europe

The polymorphism of exon 2 of the DAB genes (major histocompatibility complex [MHC] class IIB) was investigated for the first time in the freshwater cyprinid fish species, Squalius cephalus, in the wide range of its distribution in Europe. We identified 111 different MHC class IIB variants in 15 chub populations distributed from Finland to Spain. The sequence analysis showed that many structurally important amino acid sites that were conserved among tetrapods were also conserved in chub. The analysis of recombination indicated that it does not play an important role in producing and maintaining the variation of DAB genes analyzed in the present study. The exon 2 was shown to be subjected to intense positive selection. Phylogenetic analysis and sequence identities suggest the presence of two class IIB loci (DAB1-like and DAB3-like) in chub. Nevertheless, the presence of three DAB3-like sequence variants in several individuals indicates the duplication of the DAB3 gene. A contrasting selection pattern was found in DAB1-like and DAB3-like genes, which suggests the potential functional differences between these genes. Some DAB sequence variants were shared among the populations of different mtDNA lineages. The phylogenetic analyses did not confirm any biogeographical pattern of the genetic structure of MHC IIB in chub, which is in line with balancing selection and trans-species polymorphism in MHC genes. Nevertheless, cluster analysis based on the presence/absence of DAB sequence variants in the populations showed the phylogeophraphical pattern corresponding to the mtDNA lineages, which indicates that neutral selection can partially explain the MHC IIB evolution in chub.     

MHC Class IIB Exon 2 Polymorphism in the Grey Partridge (Perdix perdix) Is Shaped by Selection,Recombination and Gene Conversion

Among bird species, the most studied major histocompatibility complex (MHC) is the chicken MHC. Although the number of studies on MHC in free-ranging species is increasing, the knowledge on MHC variation in species closely related to chicken is required to understand the peculiarities of bird MHC evolution. Here we describe the variation of MHC class IIB (MHCIIB) exon 2 in a population of the Grey partridge (Perdix perdix), a species of high conservation concern throughout Europe and an emerging galliform model in studies of sexual selection. We found 12 alleles in 108 individuals, but in comparison to other birds surprisingly many sites show signatures of historical positive selection. Individuals displayed between two to four alleles both on genomic and complementary DNA, suggesting the presence of two functional MHCIIB loci. Recombination and gene conversion appear to be involved in generating MHCIIB diversity in the Grey partridge; two recombination breakpoints and several gene conversion events were detected. In phylogenetic analysis of galliform MHCIIB, the Grey partridge alleles do not cluster together, but are scattered through the tree instead. Thus, our results indicate that the Grey partridge MHCIIB is comparable to most other galliforms in terms of copy number and population polymorphism.     

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.